如何解决R中外部漂移的误差克里金？答案

【问题标题】：How to solve error kriging with external drift in R?如何解决R中外部漂移的误差克里金？
【发布时间】：2020-07-15 01:23:45
【问题描述】：

我正在尝试使用 Hydrostm 和 Hydrokrige 将克里金法与外部漂移结合使用。虽然，试图插值，我不能很好地插值。我正在使用 dem 30 m 数据。代码如下：

library(zoo)
library(hydroTSM)
library(xts)
library(sp)
library(automap)
library(gstat)
library(sp)
library(rgdal)
library(raster)
bena.ts<- read.csv("EbroPPtsDaily.csv")
bena.gis = read.csv("EbroPPgis.csv")
bena.rs <- raster("DEM-UTM.tif")
bena.sa = readOGR("subBASIN-UTM.shp")
# Covert Raster to spatial grid data frame 
y <- as(bena.rs, "SpatialGridDataFrame")
class(y)
# Make Projection according to study area
p4s <- CRS("+proj=utm +zone=43 +datum=WGS84 +units=m +no_defs +ellps=WGS84
+towgs84=0,0,0")
bena.ts1 <- as.numeric(bena.ts[180, 2:ncol(bena.ts)])
names(bena.ts1) = colnames(bena.ts[180,2:ncol(bena.ts)])
# Computing OK, over the spatial grid defined by the DEM
## Not run: 
x.ok <- hydrokrige(x.ts= bena.ts1, x.gis=bena.gis, 
                   X="EAST_ED50", Y="NORTH_ED50", sname="ID", 
                   formula=value~1,
                   subcatchments= bena.sa,
                   p4s= p4s, 
                   predictors=y,
                   ColorRamp= "Precipitation", 
                   main= "OK Precipitation on the Ebro",
                   arrow.plot= TRUE, 
                   arrow.offset= c(900000,4750000), arrow.scale= 20000,
                   scalebar.plot= TRUE, 
                   sb.offset= c(400000,4480000), sb.scale= 100000)

【问题讨论】：

请创建一个最小的自包含可重现示例。您可以使用所用方法的帮助文件中的示例。
我正在使用我自己的数据。
是的，这是个问题。这就是你现在应该在这个论坛上提问的方式。

标签： r interpolation rastervis hydrotsm

【解决方案1】：

bena.ts show daily precipitation values at different stations bena.gis shows coordinates in UTM, elevation of points and subbasins name bena.sa tells about shapefile of study area in UTM bena.rs show DEM 30 m raster file

我的问题是我可以使用高程对其进行插值。我还将我的栅格 DEM 30 m 文件转换为空间网格数据框，但它仍然没有显示任何结果。请指导我。

【讨论】：