0

イタリアのロンバルディア地方の地図が必要です。

ただし、実行するmap("italy")と、データセットが地域レベルではなく州レベルで分割されていることがわかります。ロンバルディア州は 12 の州で構成されています。を実行するとミラノ (都市と州) の地図を取得できますが、map("italy","Milan")実行するとエラーが発生しますmap("italy","Lombardy")。

私の質問は、1 つのマップに複数の州 (ロンバルディア州の州) を配置するにはどうすればよいですか? 内部境界線を削除することはできますか?

4

1 に答える 1

2

mapsヘルプから:

このイタリアのデータベースは、国連環境計画 (UNEP) の GRID-Geneva データ セットの NUTS III (欧州共同体の第三行政単位) データベースに由来します。これらは 1989 年頃に作成されたものであり、多少古くなっている可能性があります。

当時、ロンバルディアにはプロヴィンスが 10 しかなかったので、それらを使用します。

library(maps)

lombardy_provs <- map('italy', c('Milan', 'Bergamo', 'Brescia', 'Como', 
                                 'Cremona', 'Lecco', 'Mantova', 'Varese', 
                                 'Pavia', 'Sondrio'), 
                      resolution = 0)

引数は、resolution = 0利用可能な最大解像度を取得するために使用されます。そうしないと、アーティファクトが発生します。

新しいsf( simple features) パッケージを使用して、これらの領域を統合します。

まず、mapオブジェクトをsf

library(sf)
#> Linking to GEOS 3.5.0, GDAL 2.1.1, proj.4 4.9.3

lombardy_sf_provs <- st_as_sf(lombardy_provs)
plot(lombardy_sf_provs)

次に、次を使用してそれらを「結合」できますst_union()

lombardy_sf <- st_union(lombardy_sf_provs)
plot(lombardy_sf, col = 'lightblue')

この形式は、YMMV形式sfよりもあらゆる用途に適しています。maps

于 2017-07-10T12:42:54.773 に答える