Hi everybody,
I have a 3D krigged data that I want to slice it into 2D layers and open
them as a grid in Arcview. In order to do that I tried to subset the 3D into
2D :
mckrig<-as.data.frame(mc.krig)
class(mckrig)
mc5<-subset(mckrig, x==262762.1)
mc5<-subset(mc5, select=c(var1.pred,y,z))
than since write.asciigrid needs spatialgriddataframe I did:
grid<-GridTopology(c(3296542.8, 0), c(10,1),c(920, 30))
mc5=SpatialGridDataFrame(grid, data=mc5)
But when I tried it, I get this error:
Error in validityMethod(object) : unequal number of objects in full grid and
data slot
Could you please help me where is my mistake or is there a better way to
export my krigged data?
I appriciate alot if anyone can help me.
--
Figen Esen Ramirez
[[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
