How to use mesh cord in R

概要

jpmeshパッケージは、日本国内で利用される「地域メッシュ（メッシュコード）」をRから容易に利用可能にするパッケージです。地域メッシュとは、日本国土を緯度・経度により方形の小地域区画に細分することで、同一の基準による統計調査などの結果を得ることが可能となります。80kmから100mまでのメッシュコードに対応し、メッシュコードと緯度経度座標との相互変換を行います。

jpmeshの主な機能として、「緯度経度からの地域メッシュへの変換」、「地域メッシュからの緯度経度の取得」、「都道府県単位やleaflet上へのマッピング」があります。なお地域メッシュと緯度経度の変換の基礎は竹中明夫さんのRコードが基礎となっています。ここではその使い方を紹介します。

パッケージをインストールしたら環境内で利用可能にするため、読み込みを行いましょう。また、このサンプルページで用いるパッケージも同時に読み込むようにしておきます。

library(jpmesh)
library(sf)
library(purrr)

mesh_to_coords(5133)

## # A tibble: 1 × 5
##     meshcode lng_center lat_center lng_error lat_error
##   <meshcode>      <dbl>      <dbl>     <dbl>     <dbl>
## 1       5133       134.       34.3       0.5     0.333

mesh_to_coords(513377)

## # A tibble: 1 × 5
##     meshcode lng_center lat_center lng_error lat_error
##   <meshcode>      <dbl>      <dbl>     <dbl>     <dbl>
## 1     513377       134.       34.6    0.0625    0.0417

mesh_to_coords(51337783)

## # A tibble: 1 × 5
##     meshcode lng_center lat_center lng_error lat_error
##   <meshcode>      <dbl>      <dbl>     <dbl>     <dbl>
## 1   51337783       134.       34.7   0.00625   0.00417

export_mesh(5133778311)

## Geometry set for 1 feature 
## Geometry type: POLYGON
## Dimension:     XY
## Bounding box:  xmin: 133.9125 ymin: 34.65 xmax: 133.9156 ymax: 34.65208
## Geodetic CRS:  WGS 84

## POLYGON ((133.9125 34.65, 133.9156 34.65, 133.9...

coords_to_mesh(133, 34)

## <meshcode[1]>
## [1] 51330000

coords_to_mesh(133, 34, to_mesh_size = 80)

## <meshcode[1]>
## [1] 5133

coords_to_mesh(133, 34, to_mesh_size = 0.125)

## <meshcode[1]>
## [1] 51330000111

# 80kmメッシュに含まれる10kmメッシュを返します
coords_to_mesh(133, 34, to_mesh_size = 80) %>% 
  fine_separate()

## <meshcode[64]>
##  [1] 513300 513301 513302 513303 513304 513305 513306 513307 513310 513311
## [11] 513312 513313 513314 513315 513316 513317 513320 513321 513322 513323
## [21] 513324 513325 513326 513327 513330 513331 513332 513333 513334 513335
## [31] 513336 513337 513340 513341 513342 513343 513344 513345 513346 513347
## [41] 513350 513351 513352 513353 513354 513355 513356 513357 513360 513361
## [51] 513362 513363 513364 513365 513366 513367 513370 513371 513372 513373
## [61] 513374 513375 513376 513377

# 隣接するメッシュコードを同じスケールで返します
coords_to_mesh(133, 34, to_mesh_size = 80) %>% 
  neighbor_mesh()

## <meshcode[9]>
## [1] 5032 5033 5034 5132 5133 5134 5232 5233 5234

coords_to_mesh(133, 34, 0.5) %>% 
  neighbor_mesh()

## <meshcode[9]>
## [1] 503277994 503370903 503370904 513207092 513207094 513300001 513300002
## [8] 513300003 513300004

administration_mesh(code = 33, to_mesh_size = 80) %>% 
  head() %>% 
  knitr::kable(format = "markdown")

Session Information

sessionInfo()

## R version 4.4.2 (2024-10-31)
## Platform: x86_64-pc-linux-gnu
## Running under: Ubuntu 24.04.1 LTS
## 
## Matrix products: default
## BLAS:   /usr/lib/x86_64-linux-gnu/openblas-pthread/libblas.so.3 
## LAPACK: /usr/lib/x86_64-linux-gnu/openblas-pthread/libopenblasp-r0.3.26.so;  LAPACK version 3.12.0
## 
## locale:
##  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
##  [3] LC_TIME=en_US.UTF-8        LC_COLLATE=C              
##  [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
##  [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
##  [9] LC_ADDRESS=C               LC_TELEPHONE=C            
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       
## 
## time zone: Etc/UTC
## tzcode source: system (glibc)
## 
## attached base packages:
## [1] stats     graphics  grDevices utils     datasets  methods   base     
## 
## other attached packages:
## [1] purrr_1.0.2       sf_1.0-19         jpmesh_2.1.0.9000 rmarkdown_2.29   
## 
## loaded via a namespace (and not attached):
##  [1] s2_1.1.7           sass_0.4.9         utf8_1.2.4         class_7.3-22      
##  [5] KernSmooth_2.23-24 digest_0.6.37      magrittr_2.0.3     evaluate_1.0.1    
##  [9] grid_4.4.2         fastmap_1.2.0      jsonlite_1.8.9     e1071_1.7-16      
## [13] DBI_1.2.3          promises_1.3.2     fansi_1.0.6        crosstalk_1.2.1   
## [17] jquerylib_0.1.4    cli_3.6.3          shiny_1.9.1        rlang_1.1.4       
## [21] units_0.8-5        cachem_1.1.0       yaml_2.3.10        tools_4.4.2       
## [25] memoise_2.0.1      httpuv_1.6.15      buildtools_1.0.0   vctrs_0.6.5       
## [29] R6_2.5.1           mime_0.12          proxy_0.4-27       lifecycle_1.0.4   
## [33] classInt_0.4-10    leaflet_2.2.2      htmlwidgets_1.6.4  miniUI_0.1.1.1    
## [37] pkgconfig_2.0.3    bslib_0.8.0        pillar_1.9.0       later_1.4.1       
## [41] glue_1.8.0         Rcpp_1.0.13-1      xfun_0.49          tibble_3.2.1      
## [45] sys_3.4.3          knitr_1.49         xtable_1.8-4       htmltools_0.5.8.1 
## [49] maketools_1.3.1    wk_0.9.4           compiler_4.4.2