r - 使用 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"