5

R でシェープファイルを開く 2 つの基本的な方法を見つけました - と を使用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よろしいですか?

4

2 に答える 2

4

次の と の比較はmaptools、 National Center for Ecological Analysis rgdaland PBSMappingSynthesis からのもので、興味があるかもしれません。

「R で ESRI シェープファイルを読み書きする」

于 2013-10-14T21:07:51.353 に答える