data-raw/jpmesh_simple.R

In uribo/jpmesh: Utilities for Japanese Mesh Code

#################################################
# 47都道府県の地図を単純に表現したもの
#################################################
# Load Employed Packages --------------------------------------------------
pkgload::load_all()
library(dplyr)
library(tidyr)
library(testthat)
library(sf)
# mesh_rectangeがおかしい ------------------------------------------------------
# 各都道府県の位置を示す4点を含むデータセット
jpnrect <- 
  tibble::tibble(
  mesh4 = meshcode(
    c(
      # 1 - 6
      6041, # 北海道
      5840, # 青森
      5740, # 岩手
      5640, # 宮城
      5739, # 秋田
      5639, # 山形
      5540, # 福島
      5539, # 栃木
      5438, # 群馬
      5440, # 茨城
      5439, # 埼玉
      5340, # 千葉
      5339, # 東京都
      5239, # 神奈川県
      # 15 - 20
      5538,  # 新潟県
      5437, # 富山県
      5436, # 石川県
      5336, # 福井県
      5338, # 山梨
      5337, # 長野県
      # 21 - 25
      5236, # 岐阜県
      5238, # 静岡
      5237, # 愛知県
      5136, # 三重県
      5235, # 滋賀県
      # 26 - 30
      5335, # 京都
      5234, # 大阪
      5334, # 兵庫
      5135, # 奈良
      5035, # 和歌山
      # 31 - 35
      5333, # 島根
      5332, # 鳥取
      5233, # 岡山
      5232, # 広島
      5231, # 山口
      # # 36 - 39
      4933, # 徳島
      5033, # 香川
      5032, # 愛媛県
      4932, # 高知
      # # 40 - 47
      5130, # 福岡
      5129, # 佐賀
      5029, # 長崎
      5030, # 熊本
      5031, # 大分
      4931, # 宮崎
      4930, # 鹿児島
      4728 # 沖縄県
    )
  )) %>% 
  mutate(jis_code = stringr::str_pad(seq_len(47), width = 2, side = "left", pad = "0"),
         abb_name = c("HKD", "AOM", "IWT", "MYG", "AKT",
                      "YGT", "FKS", "IBR", "TCG", "GNM",
                      "SIT", "CHB", "TKY", "KNG", "NGT",
                      "TYM", "ISK", "FKI", "YMN", "NGN",
                      "GIF", "SZO", "AIC", "MIE", "SIG",
                      "KYT", "OSK", "HYG", "NAR", "WKY",
                      "TTR", "SMN", "OKY", "HRS", "YGC",
                      "TKS", "KGW", "EHM", "KUC", "FKO",
                      "SAG", "NGS", "KMM", "OIT", "MYZ",
                      "KGS", "OKN")) %>%
  meshcode_sf(mesh4) %>% 
  transmute(jis_code,
            abb_name,
            meshcode = mesh4) %>% 
  st_sf(crs = 4326)

# check -------------------------------------------------------------------
expect_s3_class(jpnrect, c("sf", "tbl", "data.frame"))
expect_equal(dim(jpnrect), c(47, 4))
expect_named(jpnrect, c("jis_code", "abb_name", "meshcode", "geometry"))
# plot(jpnrect["abb_name"])

usethis::use_data(jpnrect, overwrite = TRUE)