Displaying 1 result from an estimated 1 matches for "coordinatevalues".
2008 Aug 07
0
3d kriging et al
...dimensions(PAHdata)
hist(logTotalPAH)
lpah.vgm = variogram(logTotalPAH ~1, PAHdata)
lpah.vgm
plot(lpah.vgm )
lpah.fit = fit.variogram(lpah.vgm, model = vgm( 6, "Sph", 80, 1.5))
lpah.fit
plot(lpah.vgm, lpah.fit)
x = GridTopology(c(534500,3531400, 20), c(2,2, 0.5), c(176,126, 31))
class(x)
coordinatevalues(x)
x = SpatialGrid(x)
class(x)
grd = as(x, "SpatialPixels")
lpah.kriged = krige(logTotalPAH~1, PAHdata, x, model = lpah.fit)
class(lpah.kriged)
lpah.kriged = as(lpah.kriged, "SpatialPixelsDataFrame")
Any help would be appreciated.
...using R 2.7.0 on a windows XP desktop compu...