r - 在 GML 中获取克里金结果

Question

是否有可能获得 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 格式，对于每个点我都可以有坐标和值。

Answer 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>