Have you tried
Map1 <- Map[1:100,1:100]
?
Note that this works differently for SpatialGridDataFrames and
SpatialPixelsDataFrames; the first knows about [rows,cols,bands], the
second about [pixels, bands].
--
Edzer
Yong Li wrote:
Dear All,
I want to subset a large spatial grid (2000*2000) into many smaller
grids for further processing in my simulation models. Could anyone do it
before and any tips will be highly appreciated?
Library(rgdal)
Map <- readGDAL("yong.img") # one image from erdras imagine
Map$va <- Map$band1
Map$band1 <- NULL
????
Map1 <- a subset of Map$va based on selected cols and rows.
????......
Cheers
Yong Li
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
