台風の軌跡アニメーションを作る

動機

• 関西の実家へのwarningとして、東から西に動く台風の影響を過去から分析したかった。

結果

• 確かに今回の動きは稀であるのがわかった。
• 台風の速度と中心気圧についての関係も分析できる可能性がある。
• また他の気圧配置データと降水雨量のデータをあわせてみる必要がある。

台風の軌跡表示

2001年以降の台風軌跡

2015年以降の台風軌跡

Data

• 日本の地理情報: (TopoJSON 形式) GitHub - dataofjapan/land
• 台風位置情報: 気象庁 (csv形式) 気象庁｜台風位置表

Code

データ取得

• bash - wget

地理情報オブジェクト

• GeoPandas
• shapely

データ整理と表示

• pandas
• matplotlib

方針

• 日本の都道府県区分を polygon to multipolygon に変えて、表示。
• 台風名、時間ごとにソートして、5 data point ごとにプロット。

感想

• geojson をpython で触ったことがなかったので、いい経験になった。
• shapely object に少し戸惑った。神奈川オブジェクトなど離島ある子はMultiPolygonだぞ！

japan_geo = gpd.read_file('japan.geojson')
Allpolygons = []
for _, row in japan_geo.iterrows():
    polygons = row['geometry']
    if isinstance(polygons, shapely.geometry.multipolygon.MultiPolygon):
        polygons = [x for x in polygons]
    elif isinstance(polygons, shapely.geometry.polygon.Polygon):
        polygons = [polygons]
    else:
        raise ValueError
    Allpolygons.append(polygons)
Allpolygons = sum(Allpolygons,[])
m = MultiPolygon(Allpolygons)