run <- FALSE
if (nchar(Sys.getenv("GISRC")) > 0 &&
read.dcf(Sys.getenv("GISRC"))[1, "LOCATION_NAME"] == "nc_basic_spm_grass7") {
run <- TRUE
}
if (run) {
G <- gmeta()
print(G)
}
if (run) {
cat(getLocationProj(), "\n")
cat(getLocationProj(g.proj_WKT = FALSE), "\n")
}
if (run) {
grd <- gmeta2grd()
print(grd)
}
if (run) {
ncells <- prod(slot(grd, "cells.dim"))
df <- data.frame(k = rep(1, ncells))
mask_SG <- sp::SpatialGridDataFrame(grd, data = df)
print(summary(mask_SG))
}
Run the code above in your browser using DataLab