R 导出轮廓线,通过WMS为其提供服务
我正在尝试做自动映射。我想要一张不确定性地图,不确定性由等值线显示。我想导出这些等值线及其相关值,然后通过web地图服务显示它们。可能吗?通过以下代码,可以生成等高线R 导出轮廓线,通过WMS为其提供服务,r,mapping,wms,R,Mapping,Wms,我正在尝试做自动映射。我想要一张不确定性地图,不确定性由等值线显示。我想导出这些等值线及其相关值,然后通过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
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)
require(gstat)
需要(格子)
数据(meuse)
数据(meuse.grid)
坐标(meuse)=~x+y
网格化(meuse.grid)=~x+y
m这不是对你问题的直接回答,但也许它可以为你指明正确的方向:
可以从轮廓图生成对象
CP<- contourplot(var1.var ~ x+y,data=l,cuts=10)
有关等值线的数据存储在此处:
x <- CP$panel.args.common[1]$x
y <- CP$panel.args.common[2]$y
z <- CP$panel.args.common[3]$z
x亲爱的Andrey,谢谢你的回答,但这些值不是轮廓值。这些是对象x的输出。为什么要将等高线转换为多边形?为什么不通过您的网络地图服务将它们作为线路使用?@Spacedman。我编辑了标题。
x <- CP$panel.args.common[1]$x
y <- CP$panel.args.common[2]$y
z <- CP$panel.args.common[3]$z