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

Question

我正在尝试自动映射。我想要一张不确定性图，isolines显示哪些不确定性。我想导出这些isolines以及与它们关联的值，然后通过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

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