Открывать шейп-файлы в R с помощью rgdal всегда лучше, чем с помощью maptools?

Я нашел два основных способа открывать шейп-файлы в R — используя rgdal и maptools:

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

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

Структуры данных кажутся одинаковыми в обоих случаях (класс SpatialLinesDataFrame, пакет sp). Однако в то время как rgdal правильно считывает проекцию, maptools нет (вам, возможно, придется назначить CRS вручную):

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

Итак, зачем мне использовать maptools для открытия шейп-файлов? Я могу только ошибиться, назначив CRS вручную!

Могу ли я сделать вывод, что оба способа эквивалентны, но использование rgdal всегда безопаснее для открытия шейп-файла?


r gis shapefile rgdal maptools
person Tomas    schedule 20.07.2013    source источник

Ответы (2)


arrow_upward
4
arrow_downward

Следующее сравнение maptools, rgdal и PBSMapping из Национального центра экологического анализа и синтеза может вас заинтересовать:

"Чтение и запись шейп-файлов ESRI с помощью R"

person maj    schedule 14.10.2013

arrow_upward
1
arrow_downward

maptools::readShapeLines() получает проблему с именем файла, содержащим пробел, например. Немецкий умляут.

rgdal может справиться с таким случаем

person satart    schedule 26.07.2013