RとRcppで緯度経度から都道府県と行政区域を高速に検索(したかった) - Qiita

前回(Rで緯度経度から都道府県を求める - Qiita)のやり方でも行政区域まで10ms程度で検索できるが、たくさん検索したかったのでRcppで高速化を試みた。結果は 前回のやり方の概要 行政区域のセントロイドと目標座標の間の距離を求めておいて、最も距離の短かった行政区域の情報を参照する。 library(sf) # 日本地図のシェープファイルを読み込む f <- "data_source/N03-160101_GML/N03-16_160101.shp" japan <- st_read(f) jpn_cent <- st_centroid(japan) # セントロイド計算 cent_mat <- t(matrix(unlist(jpn_cent$geometry), nrow=2)) # 地点をmatrixに find_admin_old <- function(point){ #