使用 Geopandas 的基礎製圖

使用 Python 分析美國 Census 資料

Lee Hachadoorian

Asst. Professor of Instruction, Temple University

地理空間資料－延伸學習

載入地理空間資料

import geopandas as gpd

# Load a geospatial file
geo_state = gpd.read_file("state_computer_use.gpkg")

type(geo_state)

geopandas.geodataframe.GeoDataFrame

Geopandas DataFrame

print(geo_state.columns)

Index(['state', 'postal', 'name', 'geometry', 'total', 'has_computer',
       'desktop_laptop', 'desktop_laptop_only', 'portable_device',
       'portable_device_only', 'no_computer'],
      dtype='object')

Geopandas DataFrame

print(geo_state.head())

  state postal  ... portable_device_only no_computer
0    06     CA  ...              1052406     1263635
1    08     CO  ...               148749      168639
2    11     DC  ...                23554       32916
3    16     ID  ...                46565       67454
4    17     IL  ...               415840      640062

[5 rows x 11 columns]

Geopandas 的 geometry 欄

print(geo_state["geometry"].head())

0    (POLYGON ((-1716661.572795197 -1091585.2669098...
1    (POLYGON ((-787764.8711845676 -679501.90042785...
2    (POLYGON ((1950824.825659711 -402441.066172522...
3    (POLYGON ((-1357254.591159301 77405.1771214068...
4    (POLYGON ((720437.1241238854 -496471.759394428...
Name: geometry, dtype: object

Geopandas 繪圖

geo_state.plot()

美國本土地圖。

等值色階圖（Choropleth）

geo_state.plot(column = "has_computer")

美國本土地圖，各州以擁有電腦的家戶數著色。

等值色階圖（Choropleth）

geo_state.plot(column = "has_computer", cmap = "YlOrRd")

美國本土地圖，各州以擁有電腦的家戶數著色。

等值色階圖（Choropleth）

geo_state["pct_has_computer"] = 100 * geo_state["has_computer"]/geo_state["total"]

geo_state.plot(column = "pct_has_computer", cmap = "YlOrRd")

美國本土地圖，各州以擁有電腦家戶比例著色。

matplotlib 可用的 18 種連續型由淺至深色彩映射。

Matplotlib 連續型色彩映射

https://matplotlib.org/users/colormaps.html

一起來練習吧！

使用 Python 分析美國 Census 資料