Dear R Helpers
I'm investigation the geostatistics tools in R and have found the package
'gstat', which looks to be useful four two dimensional data. However, I
usually deal with three dimensional information. I would like to compute a 3D
variogram map using the variogram tool but cannot seem to get it to recognise
the third (z) dimension
The data file reads ok and the spatial coordinates seem to be correctly applied
(see result below)
Dt<-read.table(DtFl,header = TRUE,sep = ",",na.strings =
"-" )>
>
> #set spatial coordinates
> coordinates(Dt)= c("X","Y","Z")
> summary(Dt)
Object of class SpatialPointsDataFrame
Coordinates:
min max
X 5070.7241 6427.053
Y 12324.2617 13226.178
Z 602.1896 782.880
However when I attempt to compute the variogram map I get an error:
> Vrgrm.Mp.StdNi1<- variogram(Dt$StdREC1~1, Dt, cutoff = 1500, width = 50,
map = TRUE)
Error: dimensions do not match: locations 14397 and data 3780
While the variogram tool is touted to be able to deal with three D data perhaps
the map option does not?
Mark Murphy
Technical Director Mining and Geology
AMEC - Perth
________________________________
This email contains confidential information. The contents must
not be disclosed to anyone else except with the authority of the sender.
Unauthorised recipients are requested to maintain this confidentiality and
immediately advise the sender of any error or misdirection in transmission.
[[alternative HTML version deleted]]