はじめに

Pythonで可視化する地理空間データ

Mary van Valkenburg

Data Science Program Manager, Nashville Software School

位置情報

1854年ロンドンのコレラ流行
600人以上の死者

瘴気の雲

スノーのドットマップ

1864年のコレラ死者のドットマップ

https://myweb.rollins.edu/jsiry/London-Cholera-John_Snow-copy.jpg

このコースで学ぶこと

地理空間ポイントを散布図としてプロットする方法
geopandas を使ってジオメトリをプロットする方法
DataFrame から GeoDataFrame を構築する方法
データセットの空間結合の方法
プロットに街路地図を追加する方法
コロプレスマップの作成方法と適切な場面

経度と緯度

父と息子の身長の散布図

経線と緯線が描かれた地球

plt.scatter(schools.Longitude, 
            schools.Latitude, 
            c = 'darkgreen', 
            marker = 'p')
plt.show()

学校の位置の単純な散布図

plt.scatter(schools.Longitude, schools.Latitude, 
            c = 'darkgreen', marker = 'p')
plt.xlabel('Longitude')
plt.ylabel('Latitude')
plt.title('Nashville Public Schools')
plt.grid()
plt.show()

ラベル・タイトル・グリッド線付きの学校の位置の散布図

経度と緯度の抽出

bus_stops.head()

Stop ID    StopName       Location
4431      MCC5_11    (36.16659, -86.781996)
588       CHA7AWN    (36.165, -86.78406)
590       CHA8AWN    (36.164393, -86.785451)
541       CXONGULC   (36.162249, -86.790464)
5231      7AVUNISM   (36.163822, -86.783791)

経度と緯度の抽出

bus_stops['lat'] = [loc[0] for loc in bus_stops.Location]
bus_stops['lng'] = [loc[1] for loc in bus_stops.Location]
bus_stops.head()

Stop ID    StopName  Location                  lat         lng
4431      MCC5_11   (36.16659, -86.781996)    36.16659    -86.781996    
588       CHA7AWN  (36.165, -86.78406)        36.165      -86.78406
590       CHA8AWN  (36.164393, -86.785451)    36.164393   -86.785451
541       CXONGULC  (36.162249, -86.790464)   36.162249   -86.790464
5231      7AVUNISM  (36.163822, -86.783791)   36.163822   -86.783791

正規表現による経度・緯度の抽出

bus_stops2.head()

Stop ID        Location            
4431       MCC - BAY 11\nNashville, TN\n(36.16659, -86.78199)
588        CHARLOTTE AVE\nNashville, TN\n(36.165, -86.78406)
590        CHARLOTTE AV\nNashville, TN\n(36.164393, -86.785451)
541        CHARLOTTE\nNashville, TN\n(36.162249, -86.790464)
5231       Nashville, TN\n(36.163822, -86.783791)

正規表現による経度・緯度の抽出

lat_lng_pattern = re.compile(r'\((.*),\s*(.*)\)', flags=re.MULTILINE)
def extract_lat_lng(address):
    try:
        lat_lng_match = lat_lng_pattern.search(address)
        lat = float(lat_lng_match.group(1))
        lng = float(lat_lng_match.group(2))
        return (lat, lng)
    except:
        return (np.NaN, np.NaN)

lat_lngs = [extract_lat_lng(location)for location in \
           bus_stops2.loc[:, 'Location']]
bus_stops2['lat'] = [lat for lat, lng in lat_lngs]
bus_stops2['lng'] = [lng for lat, lng in lat_lngs]

ナッシュビルのオープンデータ

nashville.data.govのスクリーンショット

巨大な鶏の頭

練習しましょう！

Pythonで可視化する地理空間データ