r - 导出轮廓线以通过 WMS 为它们提供服务

Question

我正在尝试进行自动映射。我想要一张不确定性地图，这些不确定性由等值线显示。我想导出这些等值线以及与它们关联的值，然后通过 web 地图服务显示它们。可能吗？通过以下代码，您可以生成轮廓线

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

但我现在无法弄清楚如何提取它。也许这可以帮助你一点。