在 GML 中获取克里金结果答案

【问题标题】：Get kriging result in GML在 GML 中获取克里金结果
【发布时间】：2012-12-12 01:38:06
【问题描述】：

是否可以得到 GML 格式的坐标、插值和预测方差的克里金结果？

require(gstat)
data(meuse)
coordinates(meuse) = ~x+y
data(meuse.grid)
gridded(meuse.grid) = ~x+y
m <- vgm(.59, "Sph", 874, .04)
# ordinary kriging:
x <- krige(log(zinc)~1, meuse, meuse.grid, model = m)

x 包含坐标和插值锌值和相关的方差。我希望它采用 GML 格式，对于每个点我都可以有坐标和值。

【问题讨论】：

可重现的例子？ library("sos"); findFn("{GML format}") ?
@Ben，我加了一个例子
finzi.psych.upenn.edu/R/library/rgdal/html/writeOGR.html

标签： r

【解决方案1】：

类似：

require(rgdal)
writeOGR(x,"x.gml","x",driver="GML")

为您提供具有此类功能的 GML：

 <gml:featureMember>
    <ogr:x fid="x.0">
      <ogr:geometryProperty><gml:Point><gml:coordinates>181180,333740</gml:coordinates></gml:Point></ogr:geometry
Property>
      <ogr:var1.pred>6.49662446268745</ogr:var1.pred>
      <ogr:var1.var>0.310842119817784</ogr:var1.var>
    </ogr:x>
  </gml:featureMember>

【讨论】：