Displaying 1 result from an estimated 1 matches for "soilpt".
Did you mean:
soilpts
2010 Jun 25
1
variograms and kriging
Hello
Trying to develop variograms and kriged surfaces from a point file. Here is
what I've done so far.
library(gstat) # also loads library(sp)
library(lattice)
soilpts$x <- soilpts$UTM_X
soilpts$y <- soilpts$UTM_Y
soil.dat <- subset(soilpts, select=c(x, y, Area, BulkDensity, LOI, TP, TN,
TC, Total_Mg))
dim(soil.dat)
[1] 1292 7
coordinates(soil.dat) <- ~ x+y
gridded(soil.dat) <- TRUE
Warning messages:
1: In points2grid(points, tolerance,...