I would like to save the information generated from an ordinary kriging (processed with the function autoKrige from the automap R package), into a file extension readable in ArcGis for desktop 9.3.1 (more specifically, an .asc file - aka ESRI ASCII raster format).
I tried to search for functions inside the raster (raster) and rgdal (writeOGR) packages but without success.
require(gstat)require(automap)require(sp)#Reading data to provide a reproducible example.data(meuse)coordinates(meuse) =~ x+ydata(meuse.grid)gridded(meuse.grid) =~ x+y#Perform ordinary kriging and store results on object of type "autoKrige" "list" kriging_result = autoKrige(zinc~1, meuse, meuse.grid)plot(kriging_result) #I clipped one panel from the entire output plot.
#Extract coordinates and variable of interest to object of type "data.frame"kriging_variable = data.frame(kriging_result[[1]])head(kriging_variable)#Example of observations for the variable of interest "var1.pred"x y var1.pred var1.var var1.stdev181180 333740 777.5909 96809.44 311.1421181140 333700 839.7388 78629.08 280.4088181180 333700 773.9795 84108.76 290.0151181220 333700 710.5693 90272.17 300.4533181100 333660 911.2036 61045.30 247.0735181140 333660 836.4262 65746.42 256.4107Does anyone have a clue?
أكثر...
I tried to search for functions inside the raster (raster) and rgdal (writeOGR) packages but without success.
require(gstat)require(automap)require(sp)#Reading data to provide a reproducible example.data(meuse)coordinates(meuse) =~ x+ydata(meuse.grid)gridded(meuse.grid) =~ x+y#Perform ordinary kriging and store results on object of type "autoKrige" "list" kriging_result = autoKrige(zinc~1, meuse, meuse.grid)plot(kriging_result) #I clipped one panel from the entire output plot.
#Extract coordinates and variable of interest to object of type "data.frame"kriging_variable = data.frame(kriging_result[[1]])head(kriging_variable)#Example of observations for the variable of interest "var1.pred"x y var1.pred var1.var var1.stdev181180 333740 777.5909 96809.44 311.1421181140 333700 839.7388 78629.08 280.4088181180 333700 773.9795 84108.76 290.0151181220 333700 710.5693 90272.17 300.4533181100 333660 911.2036 61045.30 247.0735181140 333660 836.4262 65746.42 256.4107Does anyone have a clue?
أكثر...