Les données SIG utilisées dans «Statistiques sur cartes (SIG statistique)» des statistiques du gouvernement japonais vous permettent de télécharger les mêmes données dans différents formats.
shp est un ancien format binaire et GML est un format XML récemment développé. Lors de la lecture avec python geopandas (en interne fiona), utilisez shp. Il semble que le «GDAL» utilisé en interne par «fiona» ne supporte pas entièrement les spécifications XML requises.
Après cela, en raison de la différence de projection, quand j'ai regardé la définition de proj4, c'était comme ça ([Explication] de geopandas (http://geopandas.org/projections.html)). Il a été sorti dans la ville de Kobe et le système de coordonnées orthogonales du plan devrait être un peu différent selon la zone.
| Nom | chaîne proj4 |
|---|---|
| Latitude et longitude du système géodésique du Japon | +ellps=bessel +no_defs +proj=longlat |
| Latitude et longitude du système mondial de levés | +ellps=GRS80 +no_defs +proj=longlat |
| Japon système géodésique plan système de coordonnées orthogonales | +lon_0=134.3333333333333 +y_0=0 +lat_0=36 +ellps=bessel +no_defs +units=m +proj=tmerc +x_0=0 +k=0.9999 |
| Système de coordonnées orthogonales du plan du système géodésique mondial | +lon_0=134.3333333333333 +y_0=0 +lat_0=36 +ellps=GRS80 +no_defs +units=m +proj=tmerc +x_0=0 +k=0.9999 |
Le système de coordonnées orthogonales du plan est un nombre qui ressemble au sol. ʻUnitsest maintenantm`. Cela signifie-t-il que la sphère elliptique est différente entre le système géographique japonais et le système géographique mondial dans le système de coordonnées rectangulaires plan? Ce qu'on appelle le système de mesure du Japon est l'ellipse du navire, c'est-à-dire l'ancien système de mesure du Japon. Puisqu'il y a beaucoup de données WGS84 de nos jours, il semble préférable d'utiliser GRS80 du "World Survey System" qui est 2000 du Japan Survey System.
Si vous souhaitez effectuer divers calculs, vous devez télécharger "Latitude et longitude du système World Survey", le calculer et le convertir en un système de coordonnées rectangulaires plan lorsque vous faites une figure finale. Vous pouvez le faire d'un seul coup avec to_crs ().
data.to_crs("+lon_0=134.3333333333333 +y_0=0 +lat_0=36 +ellps=GRS80 +no_defs +units=m +proj=tmerc +x_0=0 +k=0.9999").plot()
Lors de l'attribution de valeurs à l'aide de l'aire, dois-je calculer dans le système de coordonnées orthogonales plan?