Dear all,
I am reading a raster file using the readGDAL function. The file has 360x364
cells.
Looking at the structure of the object (class SpatialGridDataFrame) obtained
I can see:
@ coords : num [1:2, 1:2] 4081420 4092190 2955325 2966215
I need to transform it as SpatialPixelsDataFrame or data.frame, i.e.
@ coords :num [1:130960, 1:2] 4081420 4081450 4081480 4081510 4081540
...
I tested a couple of methods:
1) dem1 <- as(x, "data.frame")
2) fullgrid(x) = FALSE
however both have some problems, i.e. the length of the object obtained is
only 130960 instead of 131040.
I think I am doing something wrong, but I can not find what...
Thank you very much in advance.
Laura
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
