[R] Geostatitics 3D Variogram Map

Murphy, Mark P (AU)
Fri, 13 Apr 2012 23:53:44 -0000

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