Closed yokoyamt closed 6 years ago
ちなみにこっちは大丈夫でした
> spdf_jpn_cities(12, admin_name = "大網白里市")
class : SpatialPolygonsDataFrame
features : 1
extent : 140.2706, 140.4305, 35.47162, 35.56063 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +ellps=GRS80 +no_defs
variables : 5
# A tibble: 1 × 5
pref_name city_name_ city_name city_name_full city_code
* <chr> <chr> <chr> <chr> <fctr>
1 千葉県 <NA> 大網白里市 大網白里市 12239
> spdf_jpn_cities(12, admin_name = "大網白里町")
class : SpatialPolygonsDataFrame
features : 0
coord. ref. : +proj=longlat +ellps=GRS80 +no_defs
variables : 5
# A tibble: 0 × 5
# ... with 5 variables: pref_name <chr>, city_name_ <chr>, city_name <chr>, city_name_full <chr>,
# city_code <fctr>
報告ありがとうございます。
元データの方に更新があったようです。
こちらで確認しましたところ、旧市町村コードの 12402
では該当するデータがなく、新市町村コード 12239
に引っかかりました。
library(magrittr)
library(jpndistrict)
#> Loading required package: jpmesh
#> This package provide map data is based on the Digital Map
#> 25000(Map Image) published by Geospatial Information Authorityof
#> Japan (Approval No.603FY2017 information usage
#> <http://www.gsi.go.jp>)
df_pref12 <-
jpn_admins(jis_code = 12)
#> options: ENCODING=cp932
#> Reading layer `P34-14_12' from data source `/private/var/folders/0x/mb63hycs4k30_7httqyxh2rr0000gn/T/Rtmplw47gK/P34-14_12_GML/P34-14_12.shp' using driver `ESRI Shapefile'
#> Simple feature collection with 203 features and 4 fields
#> geometry type: POINT
#> dimension: XY
#> bbox: xmin: 139.8206 ymin: 34.90821 xmax: 140.8268 ymax: 36.0259
#> epsg (SRID): NA
#> proj4string: +proj=longlat +ellps=GRS80 +no_defs
# 旧市町村コード
df_pref12 %>%
dplyr::filter(jis_code == 12402)
#> [1] jis_code type name address geometry
#> <0 rows> (or 0-length row.names)
# 2013-01-01以降の市町村コード
df_pref12 %>%
dplyr::filter(jis_code == 12239)
#> jis_code type name address
#> 1 12239 1 大網白里市役所 大網白里市大網115-2
#> 2 12239 2 白里出張所 大網白里市南今泉1088-1
#> geometry
#> 1 POINT (140.321 35.52171)
#> 2 POINT (140.4003 35.49253)
# 内部データからの参照
jpndistrict:::read_ksj_p34(pref_code = 12) %>%
dplyr::filter(jis_code == 12239)
#> options: ENCODING=cp932
#> Reading layer `P34-14_12' from data source `/private/var/folders/0x/mb63hycs4k30_7httqyxh2rr0000gn/T/Rtmplw47gK/P34-14_12_GML/P34-14_12.shp' using driver `ESRI Shapefile'
#> Simple feature collection with 203 features and 4 fields
#> geometry type: POINT
#> dimension: XY
#> bbox: xmin: 139.8206 ymin: 34.90821 xmax: 140.8268 ymax: 36.0259
#> epsg (SRID): NA
#> proj4string: +proj=longlat +ellps=GRS80 +no_defs
#> jis_code type name address
#> 1 12239 1 大網白里市役所 大網白里市大網115-2
#> 2 12239 2 白里出張所 大網白里市南今泉1088-1
#> geometry
#> 1 POINT (140.321 35.52171)
#> 2 POINT (140.4003 35.49253)
Created on 2018-05-01 by the reprex package (v0.2.0).
パッケージは全く悪くないんですが気づいたので書かせてください。 (パッケージ作者の判断でcloseしちゃって結構です)
http://nlftp.mlit.go.jp/ksj/gml/datalist/KsjTmplt-P34.html
なので、それ以降に市制施行された場合(2016/10/10市制施行した宮城県富谷市とか)はしょうがない気がするのですが、 それ以前に市制施行された場合も、read_ksj_p34で得られるspdfの値が若干古いことがありそうです。
例)千葉県大網白里市 旧市町村コード12402→2013/1/1市制施行により12239 addressは大網白里市になってるのにnameが町のまま