The bounding box value returned by `mesh_to_coords()` is incorrect

[uribo]

r-wakalang slackでの報告から

library(dplyr, warn.conflicts = FALSE)
library(tidyr)
library(purrr)
library(jpmesh)

df <- 
  tibble(longitude = 139.71475, 
       latitude = 35.70078, 
       mesh_size = c("80km", "10km", "1km", "500m", "250m", "125m")
) %>% 
  mutate(meshcode = pmap(list(longitude, latitude, mesh_size), 
                         coords_to_mesh)) %>% 
  unnest()

df %>% 
  mutate(mesh = pmap(list(meshcode), 
                     mesh_to_coords)) %>% 
  unnest() %>% 
  select(meshcode, mesh_size, lng_error, lat_error)
#> # A tibble: 6 x 4
#>   meshcode    mesh_size lng_error lat_error
#>   <chr>       <chr>         <dbl>     <dbl>
#> 1 5339        80km       0.5       0.333   
#> 2 533945      10km       0.0625    0.0417  
#> 3 53394547    1km        0.00625   0.00417 
#> 4 533945471   500m       0.00312   0.00208 
#> 5 5339454711  250m       0.000781  0.000521
#> 6 53394547112 125m       0.000195  0.000130

^{Created on 2019-05-15 by the reprex package (v0.2.1)}

250m, 125mのlng_errorおよびlat_errorの値は 500mのそれの1/2, 1/4でなければならないが出力がおかしい。

[uribo]

fix.

df %>% 
  mutate(mesh = pmap(list(meshcode), 
                     mesh_to_coords)) %>% 
  unnest() %>% 
  select(meshcode, mesh_size, lng_error, lat_error)
#> # A tibble: 6 x 4
#>   meshcode    mesh_size lng_error lat_error
#>   <chr>       <chr>         <dbl>     <dbl>
#> 1 5339        80km       0.5       0.333   
#> 2 533945      10km       0.0625    0.0417  
#> 3 53394547    1km        0.00625   0.00417 
#> 4 533945471   500m       0.00312   0.00208 
#> 5 5339454711  250m       0.00156   0.00104 
#> 6 53394547112 125m       0.000781  0.000521