使用 idw 的空间插值误差答案

【问题标题】：spatial interpolation error using idw使用 idw 的空间插值误差
【发布时间】：2017-01-28 05:02:12
【问题描述】：

我正在尝试按纬度和经度对海水 pH 值数据集进行空间插值：

sample<-read.csv(file="Station locations 2016.csv", header=TRUE, sep=",", strip.white=T)

head(sample)    
  Station       lat     long       pH
1     B17 -23.49174 152.0718 8.222411
2     B23 -23.49179 152.0718 8.199310
3     B26 -23.49182 152.0717 8.140428
4     B28 -23.49183 152.0717 8.100752
5     B30 -23.49185 152.0717 8.068141
6     B31 -23.49187 152.0717 8.048852

我已经根据 lat/long 数据中的现有范围创建了一个网格，并希望对 pH 值进行插值，以便生成一个颜色编码的 pH 值图。以下代码一直有效，直到使用 idw 进行空间集成步骤，此时我收到以下错误：

library(ggplot2)
library(gstat)
library(sp)

x.range <- range(sample$long)
y.range <- range(sample$lat)


x<-seq(x.range[1], x.range[2], length.out=20)
y<-seq(y.range[1], y.range[2], length.out=20)
grd<-expand.grid(x,y)


coordinates(sample) = ~long+lat
coordinates(grd) <- ~ Var1+Var2
gridded(grd) <- TRUE

plot(grd, cex=1.5)

dat.idw<-idw(formula=pH ~ 1, data=sample, newdata=grd, idp=2.0)

Error in (function (classes, fdef, mtable)  : unable to find an inherited method for function ‘idw’ for signature ‘"formula", "missing"’

我也尝试过 krige 并得到类似的错误。有关如何解决此问题的任何建议？谢谢你。

【问题讨论】：

试试idw(formula=pH ~ 1, locations=sample, newdata=grid, idp=2.0)
我得到一个类似的错误：“（函数（类，fdef，mtable）中的错误：无法为签名“公式”，“data.frame”找到函数“idw”的继承方法'

标签： interpolation spatial sp

【解决方案1】：

library(gstat)
library(sp)


lat <-  c(-23.49174, -23.49179, -23.49182, -23.49183, -23.49185, -23.49187)
long <- c(152.0718, 152.0718, 152.0717, 152.0717, 152.0717, 152.0717)
pH <- c(8.222411, 8.19931, 8.140428, 8.100752, 8.068141, 8.048852)
sample <- data.frame(lat, long, pH)


x.range <- range(sample$long)
y.range <- range(sample$lat)


x<-seq(x.range[1], x.range[2], length.out=20)
y<-seq(y.range[1], y.range[2], length.out=20)
grd<-expand.grid(x,y)


coordinates(sample) = ~long+lat
coordinates(grd) <- ~ Var1+Var2
gridded(grd) <- TRUE

proj4string(sample) <- CRS("+proj=longlat +datum=WGS84")
proj4string(grd) <- CRS("+proj=longlat +datum=WGS84")

plot(grd, cex=1.5)

dat.idw <- idw(formula=pH ~ 1, locations = sample, newdata = grd, idp = 2.0)
#> [inverse distance weighted interpolation]

plot(dat.idw)

【讨论】：