Question

我正在尝试使用Gstat包在R中进行通用cokriging。我有一个脚本，我得到了帮助，但现在我卡住了，不能从原始来源请求帮助。问题是我无法改变cokriged数据的输出分辨率。我想将插值地图导入ArcMap，而点到光栅会让我的分辨率非常低。

我的脚本如下：

  library(raster)
library(gstat)
library(sp)
library(rgdal)
library(FitAR)

加载我的数据集，其中包含坐标和采样值：

kova<-read.table("katvus_point_modif3.txt",sep="    ",header=T)
coordinates(kova)=~POINT_X+POINT_Y

在与前一个坐标相同的坐标处加载深度值，这是我的协变量：

Sygavus<-read.table("sygavus_point_cokrig.txt",sep="    ",header=T)
coordinates(Sygavus)=~POINT_X+POINT_Y

overlay <- over(kova,Sygavus)
kova$Sygavus <- overlay$Sygavus

这应该为我的插值设置边界，该文件是从ArcMap导出的shapefile：

border <- shapefile("area_2014.shp")
projection(kova)=projection(border)

这应该为cokriging创建一个网格，res =应该让我指定我希望输出的分辨率，但是无论我使用什么数字输出都不会改变。

grid <- spsample(border,type="regular",res=25)

我删除重叠点：

zero <- zerodist(kova)
kova <- kova[-zero[,2],]

我加载深度协变量光栅文件。这是从ArcMap到ascii格式的深度栅格导出：

depth <- raster("htp_depth_covar.asc")
projection(depth)=projection(border)

overlay <- extract(depth,kova)
kova$depth <- overlay

我删除了na！来自覆盖深度值的值（这些值应与相应坐标处先前加载的深度协变量表相同，但如果我将该部分保留，则脚本将停止运行）

kova <- kova[!is.na(kova$depth),]


kova.gstat <- gstat(id="Kova",formula=kova~depth,data=kova)
kova.gstat <- gstat(kova.gstat,id="Sygavus",formula=Sygavus~depth,data=kova)

var.kova <- variogram(kova.gstat)
plot(var.kova)

kova.gstat <- gstat(kova.gstat,id=c("Kova","Sygavus"),model=vgm(psill=cov(kova$kova,kova$Sygavus),model="Mat",range=12000,nugget=0))
kova.gstat <- fit.lmc(var.kova,kova.gstat,model=vgm(psill=cov(kova$kova,kova$Sygavus),model="Mat",range=12000,nugget=0))

plot(var.kova,kova.gstat$model)

overlay <- extract(depth,grid)
grid <- as.data.frame(grid)
grid$depth <- overlay
coordinates(grid)=~x1+x2
projection(grid)=projection(border)

krige <- predict.gstat(kova.gstat,grid)

spplot(krige,c("Kova.pred"))

write.table(krige, "kova.raster1.ck.csv", sep=";", dec=",", row.names=F)

任何帮助理解gstat cokriging和整个脚本将不胜感激！

Answer 1

因为您没有提供可重现的示例，我只能猜测，但我认为spsample忽略了res=25参数。请尝试n=1000，然后增加该值以获得更高的分辨率。

R：Gstat通用cokriging分辨率

1 个答案: