如何在 autoKrige 中插入 x 和 y 标签和刻度线？

Question

使用了两个数据集：

1. 3 列中的空间数据（x、y、数据）

2. 2 列 (x, y) 中的网格数据

automap 包 autoKrige 进行克里金的计算，可以在没有 x 和 y 刻度线和标签的情况下绘制：

plot(kriging_result)
automapPlot(kriging_result$krige_output, "var1.pred", sp.layout = list("sp.points", shimadata), main="OK without grids", xlab="x", ylab="y")

当我使用 ggplot2 包时，它会显示错误，但是它会计算克里金法：

mydata<-read.table("D:/.../mydata.txt",header=T,sep=",")
#Renaming desired columns:
x<-mydata[,1]
y<-mydata[,2]
waterelev<-mydata[,3]
library(gstat)
coordinates(mydata)=~x+y
library(ggplot2)
theme_set(theme_bw())
library(scales)
library(automap)
grids<-read.table("D:/.../grids.txt",header=T,sep=",")
gridded(grids)=~x+y
kriging_result = autoKrige(log(waterelev)~1, mydata)
#This line turns the log(data) back to the original data:
kriging_result$krige_output$var1.pred<-exp(kriging_result$krige_output$var1.pred)
library(reshape2)
ggplot_data = as.data.frame(kriging_result$krige_output)
ggplot(ggplot_data, aes(x = x, y = y, fill = var1.pred)) + 
   geom_raster() + coord_fixed() + 
   scale_fill_gradient(low = 'white', high = muted('blue'))

错误：

错误：美学长度必须为 1，或与 dataProblems:x,y 长度相同

Answer 1

automapPlot 函数是spplot 的包装，因此适用于spplot 的任何修复也适用于automapPlot。您可以从 spplot 文档开始。

但是，我通常使用ggplot2 包进行任何绘图工作，包括空间数据。下面是一个重现automapPlot结果的例子：

library(ggplot2)
theme_set(theme_bw())
library(scales)
library(automap)
data(meuse)
coordinates(meuse) =~ x+y
data(meuse.grid)
gridded(meuse.grid) =~ x+y
kriging_result = autoKrige(zinc~1, meuse, meuse.grid)

# Cast the Spatial object to a data.frame
library(reshape2)
ggplot_data = as.data.frame(kriging_result$krige_output)

ggplot(ggplot_data, aes(x = x, y = y, fill = var1.pred)) + 
    geom_raster() + coord_fixed() + 
    scale_fill_gradient(low = 'white', high = muted('blue'))

Answer 2

我使用了您提供的数据集，并将结果像图片一样发布：

我只是运行您的代码，直到 ggplot_data 从 spPointsDataFrame（网格化）转换为 data.frame（使用 reshape2）。 然后我一步一步构建了ggplot对象，但没有发现错误：

g=ggplot(data=ggplot_data,aes(x=x1,y=x2,fill=var1.pred))
g=g+geom_raster()
g=g+coord_fixed()
g=g+scale_fill_gradient(low = 'white', high = muted('blue'))
print(g)

这是你所期望的吗？