所以我试图使用shapefiles包的funciton convert.to.shapefile(4)将包含在R中的spatstat包中的bei数据集转换为shapefile.不幸的是,事实证明这比我想象的要难.有没有人有这方面的经验？如果你非常友好地帮助我,我会非常感激.
谢谢,
瑞安
参考文献：
spatstat,
shapefiles
解决方法
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)
(如果rgdal不可用,请在maptools包中使用writePolyShape).
语法是对象,然后是“数据源名称”(这里是当前目录,这可以是.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);小品( “SP”)库(spatstat);晕影( “spatstat”)
