1

座標と補間値、予測分散であるクリギング結果を 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 形式にして、各ポイントに座標と値を設定できるようにします。

4

1 に答える 1

4

何かのようなもの：

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>
于 2012-12-11T17:31:45.763 に答える