使用rgdal在R中打开shapefile始终比使用maptools更好？

Question

我找到了两种在R中打开shapefile的基本方法 - 使用rgdal和maptools：

# 1
require(maptools)
shape_maptools <- readShapeLines("file.shp")

# 2
require(rgdal)
shape_rgdal <- readOGR(".", "file")

两种情况下的数据结构看起来完全相同（类SpatialLinesDataFrame，包sp）。但是，虽然rgdal reads the projection properly, maptools does not（您可能需要手动分配CRS）：

> proj4string(shape_maptools)
[1] NA
> proj4string(shape_rgdal)
[1] "+proj=utm +zone=31 +ellps=intl +units=m +no_defs"

那么，为什么我会使用maptools打开形状文件？我可能只会犯错误 CRS手动！

我可以断定两种方式都是等价的，但使用rgdal总是更安全的方式来打开shapefile吗？

Answer 1

以下对国家生态分析与综合中心的maptools，rgdal和PBSMapping的比较可能会让您感兴趣：

"Read and write ESRI Shapefiles with R"

Answer 2

maptools :: readShapeLines（）会出现包含空格字符的文件名问题，例如。德国变形金刚。

rgdal可以处理此类案件