I used the rgdal package to open a geoTiff and then manipulate the
data. The result is a 3-dimensional array that I want to convert back
to a SpatialGridDataFrame in order to save it as a geoTiff. I used
the following code:
grd <- GridTopology(gt[,1],gt[,2],gt[,3])
SGDF <-
SpatialGridDataFrame(grd,data=data.frame(c(data)),proj4string=CRS(proj))
Which works fine for a single band image (a 2-Dimensional matrix).
When I used this code on a 3D array, I get the following error:
Error in validityMethod(object) :
unequal number of objects in full grid and data slot
Does anyone have a solution? Thanks.
Jonathan B. Thayn
Ph.D. Candidate Geography and Remote Sensing
University of Kansas
2101 Constant Avenue, Higuchi Hall, Rm 121
Lawrence, Kansas 66047-3759
[EMAIL PROTECTED]
people.ku.edu/~jonthayn
[[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
