r - 导出等高线以通过 WMS 为他们服务

Question

我正在尝试进行自动映射。我想要一张不确定性 map ，其中的不确定性由等值线显示。我想导出这些等值线和与它们关联的值，然后通过网络 map 服务显示它们。是否可以？通过以下代码，您可以生成轮廓线

require(gstat)
require(lattice)
data(meuse)
data(meuse.grid)
coordinates(meuse) = ~x+y
gridded(meuse.grid) = ~x+y
m <- vgm(.59, "Sph", 874, .04)
# ordinary kriging:
x <- krige(log(zinc)~1, meuse, meuse.grid, model = m)
l=as.data.frame(x)
contourplot(var1.var ~ x+y,data=l,cuts=10)

Answer 1

这不是您问题的直接答案，但也许可以为您指明正确的方向:

你可以从你的轮廓图中生成一个对象

CP<- contourplot(var1.var ~ x+y,data=l,cuts=10)

并试着看看它的结构

fix(CP)

关于等值线的数据存储在这里:

x <- CP$panel.args.common[1]$x
y <- CP$panel.args.common[2]$y
z <- CP$panel.args.common[3]$z

但我现在不知道如何提取它。也许这可以帮助你一点。