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]]


