Question

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

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

谢谢，莱恩

参考文献： spatstat， shapefiles

Answer 1

Spatial和spatstat包中有maptools个对象的转换器函数可用于此目的。 shapefile至少由每个对象的点（或线或多边形）和属性组成。

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）;小品（ “SP”）库（spatstat）;晕影（ “spatstat”）

Answer 2

一般解决方案是：

将"ppp"或"owin"分类对象转换为sp包中的适当分类对象
使用包writeOGR()中的rgdal函数编写Shapefile

例如，如果我们考虑来自hamster的{{1}}数据集：

spatstat

首先将此对象转换为require(spatstat) require(maptools) require(sp) require(rgdal) data(hamster)对象：

SpatialPointsDataFrame

这为我们提供了一个ham.sp <- as.SpatialPointsDataFrame.ppp(hamster)对象：

sp

此对象在> str(ham.sp, max = 2) Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots ..@ data :'data.frame': 303 obs. of 1 variable: ..@ coords.nrs : num(0) ..@ coords : num [1:303, 1:2] 6 10.8 25.8 26.8 32.5 ... .. ..- attr(*, "dimnames")=List of 2 ..@ bbox : num [1:2, 1:2] 0 0 250 250 .. ..- attr(*, "dimnames")=List of 2 ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots广告位中有一个变量：

@data

所以说我们现在想把这个变量写成ESRI Shapefile，我们使用> head(ham.sp@data) marks 1 dividing 2 dividing 3 dividing 4 dividing 5 dividing 6 dividing

writeOGR()

这将在当前工作目录中创建的目录writeOGR(ham.sp, "hamster", "marks", driver = "ESRI Shapefile")中创建多个marks.xxx文件。那组文件是 ShapeFile 。

我没有使用hamster数据集执行上述操作的原因之一是它不包含任何数据，因此我们无法将其强制转换为bei对象。我们可以使用数据，SpatialPointsDataFrame（同时加载bei.extra），但是这些额外的数据或在常规网格上。所以我们必须

将bei转换为bei.extra个对象（比如SpatialGridDataFrame）
将bei.spg转换为bei个对象（比如SpatialPoints）
bei.sp overlay()指向bei.sp网格，为bei.spg
应该为我们提供一个bei，可以使用上面的SpatialPointsDataFrame写出

正如您所看到的，只是为了给您一个Shapefile而涉及更多。我显示的writeOGR()数据示例是否足够？如果没有，我明天就可以查看Bivand et al {{}}}并执行hamster的步骤。

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

2 个答案: