R : How to write an XYZ file from a SpatialPointsDataFrame?

Question

我有一个 SpatialPointsDataFrame，它有一个属性(为了方便起见，我们称之为 z)以及纬度/经度坐标。

我想将其写入 XYZ 文件(即具有三列的 ASCII 文件)。

最初我试过

write.table(spdf, filename, row.names=FALSE)

但这首先在每一行上写入 z 值，然后是坐标。所以它是 ZXY 格式而不是 XYZ。也许没什么大不了的，但对于必须使用该文件的其他人来说很烦人。

目前，我正在使用一种感觉非常糟糕的 bodge 来执行此操作(如下所示)，但我的问题是:是否有一种很好且直接的方法可以将 SPDF 写为 XYZ，并且列的顺序正确？似乎应该很容易!

感谢您的任何建议。

博奇:

dfOutput <- data.frame(x = coordinates(spdf)[,1], y = coordinates(spdf)[,2])
dfOutput$z <- data.frame(spdf)[,1]
write.table(dfOutput, filename, row.names=FALSE)

Answer 1

为什么不只是

library(sp)
spdf <- SpatialPointsDataFrame(coords=matrix(rnorm(30), ncol = 2), 
                               data=data.frame(z = rnorm(15)))


write.csv(cbind(coordinates(spdf), spdf@data), file = "example.csv", 
          row.names = FALSE)