Hi All,

 

I have this question. I wish to create a kriging map with R but I haven't a
auxiliary map. I have only one txt file with X, Y and Z records. 

When I use this code: 

 

 

plot(variogram(Z~sqrt(X+Y), subground, cutoff=1800, width=10))

subground.or <- variogram(Z~sqrt(X+Y), subground, cutoff=1800, width=10)

subground.ovgm <- fit.variogram(subground.or, vgm(1, "Gau", 1000, 1))

subground.ok <- krige (Z~sqrt(X+Y), subground, subground.ovgm, nmax=40) *

 

*[the original code is elevations.ok = krige(Z~1, elevations, dem,
elevations.ovgm, nmax=80) ]

 

when I use this one there is this message (with my data setting), because I
have only txt file and I haven't a DEM ascii file: 

 

"Error in proj4string(newdata) :  proj4string only works for class(es
extending) Spatial." 

 

Is it possible use only one file with krige?

 

Thank you

 

Ale

 


	[[alternative HTML version deleted]]

