gis操作について

rstudioでgisを触る機会が増えてきたけど、勉強してはすぐ忘れてしまう。メモ。

 

raster::rasterFromXYZ(long, lat, value)

緯度経度つきの環境データをラスターに変換するとき使うことあり。

 

crs関係はここをよく使う

JGD座標系とEPSGの一覧表 | OpenなGISのこと

 

メッシュ⇔緯度経度の変換

japanmeshライブラリが良き。jpmeshはかなり遅い…

 

元の.shpファイルにイカれたオブジェクトがあることがあるので、まずはじめにsf::st_make_validでそのオブジェクトを「大丈夫」にしてあげる。たまにそれすらできないオブジェクトがあるのでそのときはsf::st_is_validでNAになったオブジェクトをこの世から消し去ると良い。