RでGIS演算（例えば面積を算出）を行う場合には座標参照系を確認しておく必要がある。 adm_area <- sf::st_read("N03/2016/N03-16_17_160101.shp", crs = 4612, options = "ENCODING=UTF-8") adm_area %>% as(., "Spatial") %>% rgeos::gArea(byid = T) %>% …

GISデータでメッシュ単位でデータを公開している場合があるが、 メッシュデータにはいくつかの決まり事があるらしい。 詳しくは総務省統計局、パスコのHPやArcGIS Tipsがわかりやすい。 代表的なものとして第1次～第3次メッシュがあり、以下のようになってい…

何回も調べたりしているため、メモ。 ポリゴン内のライン長さを出したいときがある。 その際はラインをポリゴンでインターセクトする必要がある。 関数はrgeos::gIntersection()を利用する。 sp::over()関数で行うと意図と異なる結果となるため注意。 （適切…