如何将样本数据集从R包"spatstat"转换为shapefile

Question

我在Java中编写了一个内核密度估计器,它以ESRI shapefile的形式输入并输出估计表面的GeoTIFF图像.为了测试这个模块,我需要一个示例shapefile,无论出于什么原因,我被告知要从R中包含的示例数据中检索一个.问题是没有任何示例数据是shapefile ...

所以我试图使用shapefiles包的convert.to.shapefile(4)函数将包含在R中的spatstat包中的bei数据集转换为shapefile.不幸的是,事实证明这比我想象的要难.有没有人有这方面的经验？如果你非常友好地帮助我,我会非常感激.

谢谢,瑞恩

参考文献: spatstat, shapefile

Answer 1

可以使用和包中的Spatial对象的转换器函数.shapefile至少由每个对象的点(或线或多边形)和属性组成.spatstatmaptools

library(spatstat)
library(sp)
library(maptools)
data(bei)

强制bei的Spatial对象,因为没有上"标记"在这里只是点无属性的ppp对象.

spPoints <- as(bei, "SpatialPoints")

shapefile至少需要一列属性数据,因此创建一个虚拟数据.

dummyData <- data.frame(dummy = rep(0, npoints(bei)))

使用SpatialPoints对象和虚拟数据,生成一个SpatialPointsDataFrame.

spDF <- SpatialPointsDataFrame(spPoints, dummyData)

此时,您一定要考虑使用的坐标系是什么,bei以及是否可以用WKT CRS(众所周知的文本坐标参考系统)表示.您可以将其Spatial作为另一个参数分配给对象SpatialPointsDataFrame,或者在创建后使用proj4string(spDF) <- CRS("+proj=etc...")(但这是一个完整的问题,我们可以自己编写页面).

加载rgdal包(这是最常用的选项,因为它支持许多格式并使用GDAL库,但由于系统依赖性而可能无法使用.

library(rgdal)

(如果没有,请writePolyShape在maptools包装中rgdal使用).

语法是对象,然后是"数据源名称"(这里是当前目录,这可以是.shp或文件夹的完整路径),然后是图层(用于shapefile文件名而不带扩展名),以及然后是输出驱动程序的名称.

writeOGR(obj = spDF, dsn = ".", layer = "bei", driver = "ESRI Shapefile")

请注意,如果"bei.shp"已经存在,写入将失败,因此必须先删除unlink("bei.shp").

列出以"bei"开头的所有文件:

list.files(pattern = "^bei")

[1] "bei.dbf" "bei.shp" "bei.shx"

请注意,ppp对象没有通用的"as.Spatial"转换器,因为必须决定这是否是带有标记的点模式等等 - 尝试编写一个可能很有趣,报告虚拟数据是否为要求等等.

有关这些数据表示之间差异的更多信息和详细信息,请参阅以下插图:

库(SP); vignette("sp")库(spatstat); 晕影( "spatstat")