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]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
