use_sp()
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")
}
if (run) {
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