読者です 読者をやめる 読者になる 読者になる

RでのGIS演算時には座標参照系の確認が大切

r gis

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

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

ポリゴン内のラインデータの抽出(切取り)

r gis

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