A non-existent mesh cord?

[uribo]

Slack (r-wakalang) より。

library(jpmesh)
library(sf)
#> Linking to GEOS 3.7.1, GDAL 2.4.0, PROJ 5.2.0

# 以下の結果が同じになる (500mメッシュ)
all.equal(
  export_mesh("513583021") %>% 
    st_as_sf(), 
  export_mesh("523503021") %>% 
    st_as_sf())
#> [1] TRUE

# 問題は10kmメッシュ以下で発生している
all.equal(
  export_mesh("513583") %>% 
    st_as_sf(), 
  export_mesh("523503") %>% 
    st_as_sf())
#> [1] TRUE

all.equal(
  export_mesh("5135") %>% 
    st_as_sf(), 
  export_mesh("5235") %>% 
    st_as_sf())
#> [1] "Component \"x\": Attributes: < Component \"bbox\": Mean relative difference: 0.01941743 >"
#> [2] "Component \"x\": Component 1: Component 1: Mean relative difference: 0.01945523"

^{Created on 2020-04-05 by the reprex package (v0.3.0)}

10kmメッシュの下2桁は0~7の範囲なので、 513583 の段階で処理を停止しなければいけない。

[uribo]

同じ問題は neighbor_mesh()でも

library(jpmesh)

neighbor_mesh(513474923, contains = TRUE)
#> [1] "513474912" "513474914" "513474921" "513474922" "513474923" "513474924"
#> [7] "513484012" "513484021" "513484022"

^{Created on 2020-04-05 by the reprex package (v0.3.0)}

xxxx8xxx が出力されるのはおかしい。