Pyhthon-foliumのコロプレス図の作成にて、色を付けられない

Question

地理区画情報の入ったjsonファイルにcsvファイルの数値を反映させる

以下のリンクを元に、次のようなコロプレス図を作成しようと考えてました。

https://akatak.hatenadiary.jp/entry/2018/08/26/093839

使用したjsonファイルは東京２３区のポリゴンの情報に加えて
行政区画コードでkeyが設定されています。以下の場所から拝借。
https://github.com/niiyz/JapanCityGeoJson/tree/master/geojson/13

該当するコードがこちらです。

import folium
geojson_23 = r"\tokyo23.json" #ファイルの相対パス指定
#Stamen Tonerのコピーライト
copyright_st = '© ' \
            'Map tiles by <a href="http://stamen.com">Stamen Design</a>,' \
            ' under <a href="http://creativecommons.org/licenses/by/3.0">CC BY 3.0</a>.' \
            'Data by <a href="http://openstreetmap.org">OpenStreetMap</a>,' \
            'under <a href="http://www.openstreetmap.org/copyright">ODbL</a>.'

tokyo23_location = [35.658593, 139.745441]   # 地図の中心点を指定。

f = folium.Map(location=tokyo23_location,tiles='cartodbpositron',zoom_start=10,) 
import pandas as pd
tokyo_data=pd.read_csv(\tokyo_data.csv",header=0)#相対パス指定
tokyo_data['行政区分コード'] = tokyo_data['行政区分コード'].astype('str')

f.choropleth(
    geo_data=geojson_23,
    name='choropleth',
    data=tokyo_data,# DataFrameまたはSeriesを指定
    columns=['行政区分コード', '物件数'],  # 行政区分コードと表示データ
    key_on='feature.id',# GeoJSONのキー（行政区分コード）
    fill_color='YlGn',# 色パレットを指定（※）
    threshold_scale=[0,4000, 8000, 12000, 16000, 20000], # 境界値を指定
    fill_opacity=0.7, # 透明度（色塗り）
    line_opacity=0.2,# 透明度（境界） 
    legend_name='物件数'# 凡例表示名
)
f.save('sample_23.html')

色がつかないエラー

以上のことを決行すると、コード自体は実行してくれますが、細かいエラー文と共に
求めていたものと違うものが出てきました。

・・・\folium.py:415: FutureWarning: The choropleth  method has been deprecated.
 Instead use the new Choropleth class, which has the same arguments. 
　See the example notebook 'GeoJSON_and_choropleth' for how to do this.
  FutureWarning

これの意味が少々わかっていません。

試したこと

何に問題があるのかを明確にするために、次のようなことを考えました。
①csvファイルに問題がある
②jsonファイルに問題がある。
③コードに問題がある。

①について
csvファイルがおかしいならば、数値を読み込んでくれずに着色できないと考え、
他のjsonファイルで同じことを試した。
すると、色がつかないでもないことが分かりました。

これを踏まえて、表したいものに対して複数回実行して重ねれば
上記のコードで行う必要はないと考えたのですが試行してみると
ゲージが二つ表示されるなどで、使用できないことが分かりました。

となると、コードはとりあえず動いてるのかなとも考えて
jsonファイルを覗いたのですが、ちんぷんかんぷんでした。

補足情報（FW/ツールのバージョンなど）

OS windows10

使用ツールjupyter notebook version　4.4.0です。

。

Answer 1

https://github.com/niiyz/JapanCityGeoJson/tree/master/geojson/13

からのGEOJSONデータを確認してみたところ、
一部の区において、同じIDで複数個エリア（Polygon）が飛地のように定義されているようです。

こちらの仕様書にて仕様を確認してみたところ、

For Polygons with more than one of these rings, the first MUST be the exterior ring, and any others MUST be interior rings.

複数のポリゴンの場合、最初のリングは外周でなければならず、他は内部のリングでなければなりません。（意訳）

との記述がありましたので、どうもこのデータはGeoJsonの仕様に合致していないような気がしています。
（専門家ではないので明言はできませんが・・）

そこで、仕様書の次の項目にある "MultiPolygon" にこれらのデータを書き換えてみました。
私が変換に使ったコードはこんな感じです。

import json

with open('tokyo23.json') as f:
    d = json.load(f)

features = d['features']
# idリスト
ids = sorted(set([d['id'] for d in features]))
new_features = []
for id in ids:
    data = [d['geometry']['coordinates'] for d in features if d['id']==id]
    if len(data) == 1:
        feature = dict(
            type = "Feature",
            geometry = dict(
                type = "Polygon",
                coordinates = data[0],
            ),
            id = id,
        )
    else:
        feature = dict(
            type = "Feature",
            geometry = dict(
                type = "MultiPolygon",
                coordinates = data,
            ),
            id = id,
        )
    new_features.append(feature)
d['features'] = new_features

with open('new_tokyo23.json', 'w') as f:
    json.dump(d, f)

とりあえず私の環境では、このコードにて吐き出された変換されたデータにて問題なく動作するようです。

データは乱数で生成しております

Comments

[takumammoth]

横からすいません。
私も質問者の方と同じ問題で困っていまして、
>仕様書の次の項目にある "MultiPolygon" にこれらのデータを書き換えてみました。
の部分が分かりません。
どこをどう書き換えているのかもう少し詳細に教えていただけないでしょうか。
お手数お掛け致しますが、よろしくお願い申し上げます。
←解決しました。ありがとうございます。

[magichan]

既に解決済みとのことですが、一応。
元のデータが
{
"type": "FeatureCollection",
"features": [
{
"type": "Feature",
"geometry": {
"type": "Polygon",
"coordinates": [(AREA1のデータ)]
},
"id": ""ID1"
},
{
"type": "Feature",
"geometry": {
"type": "Polygon",
"coordinates": [(AREA2のデータ)]
},
"id": "ID1" ;;<- １つ上と同じID
},
(以下省略)
]
}
のように同じIDの　Polygon データが並んでいる部分を
{
"type": "FeatureCollection",
"features": [
{
"type": "Feature",
"geometry": {
"type": "MultiPolygon",
"coordinates": [[(AREA1のデータ)],[(AREA2のデータ)]]
},
"id": "同じID"
},
(以下省略)
]
}
のように MultiPolygonでまとめたということです

[magichan]

あっコメントはインデントがなくなるんだった・・