简介

在 Python 中可视化地理空间数据

Mary van Valkenburg

Data Science Program Manager, Nashville Software School

位置

  • 1854 年伦敦霍乱爆发
  • 超过 600 人死亡

瘴气云示意图

在 Python 中可视化地理空间数据

斯诺的点图

1864 年霍乱死亡分布点图

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

在 Python 中可视化地理空间数据

本课程将学到什么

  • 如何将地理点绘制为散点图
  • 如何用 geopandas 绘制几何对象
  • 如何从 DataFrame 构建 GeoDataFrame
  • 如何进行空间连接
  • 如何为图添加街道底图
  • 何时以及如何创建分级填色图
在 Python 中可视化地理空间数据

经度与纬度

父子身高散点图

带经纬度网格的地球

在 Python 中可视化地理空间数据
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()

带标签、标题、网格的学校散点图

在 Python 中可视化地理空间数据

提取经度和纬度

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)
在 Python 中可视化地理空间数据

提取经度和纬度

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
在 Python 中可视化地理空间数据

用正则表达式提取经纬度

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)
在 Python 中可视化地理空间数据

用正则表达式提取经纬度

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]
在 Python 中可视化地理空间数据

纳什维尔开放数据

nashville.data.gov 截图

巨型鸡头

在 Python 中可视化地理空间数据

让我们来练习！

在 Python 中可视化地理空间数据

Preparing Video For Download...