# 回遊データをラインポリゴン化 kaiyuu_line_polygons <- postings_with_kaiyuu_id %>% sf::st_as_sf( coords = c('longitude', 'latitude') ) %>% group_by( # group_by + summariseをすると個々の投稿座標がMULTIPOINTでまとまる kaiyuu_id ) %>% summari…

postings_with_kaiyuu_id %>% sf::st_as_sf( coords = c('longitude', 'latitude') ) %>% group_by( # group_by + summariseをすると個々の投稿座標がMULTIPOINTでまとまる kaiyuu_id ) %>% summarise( dummy = n(), do_union = F # 元のデータの順番でMULTI…

# 日本地図 japan <- sf::st_read( 'https://oku.edu.mie-u.ac.jp/~okumura/stat/data/japan.geojson', stringsAsFactors = F )

japanmesh でメッシュコードから緯度経度をだしたい その作業を tibble でやりたいとき data %>% mutate( latitude = map_dbl( .x = mesh_code, .f = ~ {jpgrid::grid_to_coords(jpgrid::grid_auto(.x))[[2]]} ), longitude = map_dbl( .x = mesh_code, .f =…

地理空間データの形式として .shp 形式をよく利用してた .shp はどうも重い？ wkt はオブジェクトとそのオブジェクトの頂点座標の集合を1つのデータセットとしてもつ形式 ファイルの形式そのものは .csv になる あとでちゃんと調べよう

gisをしていると、接触するポリゴンを1つのポリゴンとしてまとめたいときがある。sfパッケージでできるのだが、どうも思ったようにはいまのところ動かないときがある。 例えば、geom_sfで可視化したときに、接触していないようなポリゴンが接触したと判定さ…

japanmeshはjpmeshより何かと処理速度が速くて助かる。mesh_codeさえわかれば、それをポリゴンにしてsf形式のtibbleにすることもできる。こんな感じ。 japanmeshがなくなってjpgridになったver (data with meshcode) %>% mutate( geometry = map( .x = mesh_…

描画する緯度経度で範囲指定したい coord_sf( xlim = lonlim, ylim = latlim ) 緯度経度の刻みを調整したい scale_x_continuous(breaks = seq( trunc(lonlim[1] / 10) * 10, floor(lonlim[2] / 10) * 10, 5)) + scale_y_continuous(breaks = seq( trunc(latl…

distances::nearest_neighbor_search()をつかってみた。引数で、distances::distances()で作った距離行列、検索する近傍点の数、検索する範囲を指定できる。 distances::nearest_neighbor_search( distances = distances::distances(座標データフレーム), k …

rstudioでgisを触る機会が増えてきたけど、勉強してはすぐ忘れてしまう。メモ。 raster::rasterFromXYZ(long, lat, value) 緯度経度つきの環境データをラスターに変換するとき使うことあり。 crs関係はここをよく使う JGD座標系とEPSGの一覧表 | OpenなGISの…