Rで解析：地図などのエリアプロットを簡単にメッシュで分割できます「geoaxe」パッケージ

地図などの空間を管理する「SpatialPolygons」classのプロットにメッシュ分割を追加するパッケージの紹介です。

なお、紹介では本パッケージの他に空間描写に便利な「rgeos」と「sp」パッケージを利用しています。使用環境によるかもしれませんが、これらパッケージは「geoaxe」パッケージのインストール時に自動でインストールされます。

パッケージバージョンは0.1.2.91。実行コマンドはwindows 11のR version 4.1.3で確認しています。

パッケージのインストール

下記、コマンドを実行してください。

#パッケージのインストール
install.packages("devtools")
devtools::install_github("ropenscilabs/geoaxe")

実行コマンド

詳細はコメント、パッケージのヘルプを確認してください。

#&#12497;&#12483;&#12465;&#12540;&#12472;&#12398;&#35501;&#12415;&#36796;&#12415;
library("geoaxe")
#install.packages("rgeos")
library("rgeos")
#install.packages("sp")
library("sp")

#&#25551;&#20889;&#12487;&#12540;&#12479;&#12398;&#20316;&#25104;:&#12300;rgeos&#12301;&#12497;&#12483;&#12465;&#12540;&#12472;;readWKT&#12467;&#12510;&#12531;&#12489;
#&#21442;&#32771;:&#25551;&#20889;&#20869;&#23481;&#12398;&#25351;&#23450;:text&#12458;&#12503;&#12471;&#12519;&#12531;;&#20869;&#23481;&#12399;&#12486;&#12461;&#12473;&#12488;&#12391;&#25351;&#23450;&#12375;&#12414;&#12377;
#&#12394;&#12362;,&#25551;&#20889;&#20869;&#23481;&#12399;"POINT","LINESTRING","POLYGON","MULTIPOINT","MULTILINESTRING","MULTIPOLYGON"&#12398;&#35373;&#23450;&#12364;&#21487;&#33021;&#12391;&#12377;
#&#12300;geoaxe&#12301;&#12497;&#12483;&#12465;&#12540;&#12472;&#12398;chop&#12467;&#12510;&#12531;&#12489;&#21033;&#29992;&#12398;&#12383;&#12417;&#12395;"POLYGON"&#12434;&#25351;&#23450;
PolyData <- readWKT(text = "POLYGON((-180 -20, -140 155, 10 0, -140 -60, -180 -20))")
#&#12503;&#12525;&#12483;&#12488;
plot(PolyData)

#&#12509;&#12522;&#12468;&#12531;&#12487;&#12540;&#12479;&#12434;&#20998;&#21106;:chop&#12467;&#12510;&#12531;&#12489;
#&#20998;&#21106;&#12469;&#12452;&#12474;&#12398;&#25351;&#23450;:size&#12458;&#12503;&#12471;&#12519;&#12531;;&#21021;&#26399;&#20516;10
#&#20998;&#21106;&#25968;&#12398;&#25351;&#23450;:n&#12458;&#12503;&#12471;&#12519;&#12531;;&#21021;&#26399;&#20516;20
polys <- chop(x = PolyData, size = 15, n = 15)
#&#12503;&#12525;&#12483;&#12488;
plot(polys, col = "#e1e6ea")

#size,n&#12458;&#12503;&#12471;&#12519;&#12531;&#20516;&#12395;&#21512;&#12431;&#12379;&#12390;&#12509;&#12522;&#12468;&#12531;&#12487;&#12540;&#12479;&#12398;&#12503;&#12525;&#12483;&#12488;&#12399;&#22793;&#21270;&#12375;&#12414;&#12377;
#size&#12458;&#12503;&#12471;&#12519;&#12531;&#12398;&#22793;&#21270;&#20363;
layout(matrix(c(1, 2, 3, 4), 2, 2, byrow = TRUE))
par(mar = c(1, 0, 1, 0))
plot(chop(polys, size = 5, n = 15), col = "#e1e6ea")
plot(chop(polys, size = 10, n = 15), col = "#505457")
plot(chop(polys, size = 15, n = 15), col = "#4b61ba")
plot(chop(polys, size = 50, n = 15), col = "#a87963")

#n&#12458;&#12503;&#12471;&#12519;&#12531;&#12398;&#22793;&#21270;&#20363;
layout(matrix(c(1, 2, 3, 4), 2, 2, byrow = TRUE))
par(mar = c(1, 0, 1, 0))
plot(chop(polys, n = 5, size = 15), col = "#e1e6ea")
plot(chop(polys, n = 10, size = 15), col = "#505457")
plot(chop(polys, n = 15, size = 15), col = "#4b61ba")
plot(chop(polys, n = 50, size = 15), col = "#a87963")

出力例

・chopコマンド

・sizeオプションの変化例

・nオプションの変化例

少しでも、あなたのウェブや実験の解析が楽になりますように！！