if (nchar(Sys.getenv("GISRC")) > 0) {
G <- gmeta6()
print(G)
CRS(getLocationProj())
grd <- gmeta2grd()
print(grd)
ncells <- prod(slot(grd, "cells.dim"))
df <- data.frame(k=rep(1, ncells))
mask_SG <- SpatialGridDataFrame(grd, data=df)
print(summary(mask_SG))
}
Run the code above in your browser using DataLab