[数据分析与可视化] Python绘制数据地图2-GeoPandas地图可视化

本文所有代码见：Python-Study-Notes

pip install geopandas

conda install geopandas

conda install --channel conda-forge geopandas

# jupyter notebook环境去除warning
import warnings
warnings.filterwarnings("ignore"

# 查看geopandas版本
import geopandas as gpd

gpd.__version__

'0.10.2'

1 基础绘图

1.1 绘图接口说明

GeoPandas基于matplotlib库封装高级接口用于制作地图，GeoSeries或GeoDataFrame对象都提供了plot函数以进行对象可视化。与GeoSeries对象相比，GeoDataFrame对象提供的plot函数在参数上更为复杂，也更为常用。
GeoDataFrame对象提供的plot函数的常用输入参数如下：

def GeoDataFrame.plot(
	column: str, np.array, pd.Series (default None, # 用于绘图的列名或数据列
	kind: str, # 绘图类型
	cmap: str, # matplotlib的颜色图Colormaps
	color: str, np.array, pd.Series (default None, # 指定所有绘图对象的统一颜色
	ax: matplotlib.pyplot.Artist (default None, # 指定matplotlib的绘图轴
	cax: matplotlib.pyplot Artist (default None, # 设置图例的坐标轴
	categorical: bool (default False, # 是否按照类别设置对象颜色
	legend: bool (default False, # 是否显示图例，如果column或color参数未赋值，则此参数无效
	scheme: str (default None, # 用于控制分层设色
	k：int (default 5, # scheme的层次数
	vmin：None or float (default None, # 图例cmap的最小值
	vmax：None or float (default None, # 图例cmap的最大值
	markersize：str or float or sequence (default None, # 绘图点的大小
	figsize: tuple of integers (default None, # 用于控制matplotlib.figure.Figure
	legend_kwds: dict (default None, # matplotlib图例参数
	missing_kw: dsdict (default None, # 缺失值区域绘制参数
	aspect:‘auto’, ‘equal’, None or float (default ‘auto’, # 设置绘图比例
	**style_kwds: dict, # 其他参数，如对象边界色edgecolor, 对象填充色facecolor, 边界宽linewidth,透明度alpha
->ax: matplotlib axes instance

GeoSeries对象提供的plot函数的常用输入参数如下：

def GeoSeries.plot(
	s: Series, # GeoSeries对象
	cmap: str (default None,
	color: str, np.array, pd.Series, List (default None,
	ax: matplotlib.pyplot.Artist (default None,
	figsize: pair of floats (default None,
	aspect: ‘auto’, ‘equal’, None or float (default ‘auto’,
	**style_kwds: dict
->ax: matplotlib axes instance

想要更好使用以上参数最好熟悉matplotlib，如有不懂可以查阅matplotlib文档。下面简要介绍GeoPandas中绘图函数的使用。

读取数据集

import geopandas as gpd
# 读取自带世界地图数据集
world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'
# 人口，大洲，国名，国家缩写，ISO国家代码，gdp，地理位置数据
world.head(

world.plot(

<matplotlib.axes._subplots.AxesSubplot at 0x7f17b270fe10>

# 读取自带世界各大城市数据集
cities = gpd.read_file(gpd.datasets.get_path('naturalearth_cities'
# 名字，地理位置数据
cities.head(

分区统计图

# 去除南极地区
world = world[(world.pop_est>0 & (world.name!="Antarctica"]
# 计算人均gpd
world['gdp_per_cap'] = world.gdp_md_est / world.pop_est
# 绘图
world.plot(column='gdp_per_cap'

<matplotlib.axes._subplots.AxesSubplot at 0x7f17b241f390>

图例设置

import matplotlib.pyplot as plt

fig, ax = plt.subplots(1, 1
# 根据人口绘图数绘制
world.plot(column='pop_est', ax=ax, legend=True

<matplotlib.axes._subplots.AxesSubplot at 0x7f17b23585d0>

from mpl_toolkits.axes_grid1 import make_axes_locatable

fig, ax = plt.subplots(1, 1

divider = make_axes_locatable(ax

# 设置图例，right表示位置，size=5%表示图例宽度，pad表示图例离图片间距
cax = divider.append_axes("right", size="5%", pad=0.2

world.plot(column='pop_est', ax=ax, legend=True, cax=cax

<matplotlib.axes._subplots.AxesSubplot at 0x7f17b21de650>

import matplotlib.pyplot as plt

fig, ax = plt.subplots(1, 1

world.plot(column='pop_est',
           ax=ax,
           legend=True,
           legend_kwds={'label': "Population by Country and Area",
                        'orientation': "horizontal"}

<matplotlib.axes._subplots.AxesSubplot at 0x7f17b20ec110>

颜色设置

# 设置颜色图为tab20
world.plot(column='gdp_per_cap', cmap='tab20'

<matplotlib.axes._subplots.AxesSubplot at 0x7f17b205a490>

world.boundary.plot(

<matplotlib.axes._subplots.AxesSubplot at 0x7f15ac1cecd0>

缺失数据绘制

import numpy as np

world.loc[world[world.continent=="Africa"].index, 'gdp_per_cap'] = np.nan
world.plot(column='gdp_per_cap'

<matplotlib.axes._subplots.AxesSubplot at 0x7f15ac0afe50>

# 将缺失值地区颜色设置为lightgrey
world.plot(column='gdp_per_cap', missing_kwds={'color': 'lightgrey'}

<matplotlib.axes._subplots.AxesSubplot at 0x7f15ac094250>

# 更复杂的例子，hatch表示内部图像样式
ax = world.plot(column='gdp_per_cap', missing_kwds={'color': 'lightgrey','edgecolor':'red',"hatch": "///","label": "Missing values"}
# 不显示坐标轴
ax.set_axis_off(;

图层设置

# 以world数据集为准，对齐crs
cities.plot(marker='*', color='green', markersize=5
cities = cities.to_crs(world.crs

# base为绘图轴
base = world.plot(color='white', edgecolor='black'

cities.plot(ax=base, marker='o', color='red', markersize=5;

import matplotlib.pyplot as plt


# 设置绘图轴
fig, ax = plt.subplots(
ax.set_aspect('equal'

world.plot(ax=ax, color='white', edgecolor='black'
cities.plot(ax=ax, marker='o', color='red', markersize=5

plt.show(

图层顺序控制

# cities先绘制，则图层顺序更靠下，导致一些数据点被world图层遮盖。
ax = cities.plot(color='k'

world.plot(ax=ax;

# 调整绘图顺序，完整显示cities数据点
ax = world.plot(

cities.plot(ax=ax,color='k'

<matplotlib.axes._subplots.AxesSubplot at 0x7f159d535f10>

# 通过zorder参数设置图层顺序，完整显示cities数据点
ax = cities.plot(color='k',zorder=2

world.plot(ax=ax,zorder=1

<matplotlib.axes._subplots.AxesSubplot at 0x7f159d4e8a10>

Pandas Plots

gdf = world.head(10

gdf.plot(kind='scatter', x="pop_est", y="gdp_md_est"

<matplotlib.axes._subplots.AxesSubplot at 0x7f159d5d69d0>

1.2 绘图实例之中国地图绘制

基于geopandas的空间数据分析之基础可视化
• geopandas 中国地图绘制

1 读取地图数据

import geopandas as gpd
import matplotlib.pyplot as plt

# 读取中国地图数据，数据来自DataV.GeoAtlas，将其投影到EPSG:4573
data = gpd.read_file('https://geo.datav.aliyun.com/areas_v3/bound/100000_full.json'.to_crs('EPSG:4573'
# 中国有34个省级行政区，最后一条数据为南海九段线
data.shape

(35, 8

# 保存各个省级行政区的面积，单位万平方公里
data['area'] = data.area/1e6/1e4
# 中国疆域总面积，和官方数据会有差距
data['area'].sum(

989.1500769770855

# 查看最后五条数据
data.tail(

# 拆分数据
nine_dotted_line = data.iloc[-1]
data = data[:-1]
nine_dotted_line

adcode                                                     100000_JD
name                                                                
adchar                                                            JD
childrenNum                                                      NaN
level                                                            NaN
parent                                                           NaN
subFeatureIndex                                                  NaN
geometry           MULTIPOLYGON (((20309263.208229765 2708129.154...
area                                                        0.280824
Name: 34, dtype: object

2 地图绘制

# 创建画布matplotlib
fig, ax = plt.subplots(figsize=(12, 9
# 绘制主要区域
ax = data.plot(ax=ax
# 绘制九段线
ax = gpd.GeoSeries(nine_dotted_line.geometry.plot(ax=ax,edgecolor='red',linewidth=3
# 保存结果
fig.savefig('res.png', dpi=300, bbox_inches='tight'

3 绘图自定义

# 创建画布matplotlib
fig, ax = plt.subplots(figsize=(12, 9
# 绘制主要区域
ax = data.plot(ax=ax,facecolor='grey',edgecolor='lightgrey',alpha=0.8,linewidth=1
# 绘制九段线
ax = gpd.GeoSeries(nine_dotted_line.geometry.plot(ax=ax,edgecolor='grey',linewidth=3
# 强调首都
ax = data[data.name=="北京市"].representative_point(.plot(ax=ax, facecolor='red',marker='*', markersize=150 

# 强调港澳台
# 设置字体
fontdict = {'family':'FZSongYi-Z13S', 'size':8, 'color': "blue",'weight': 'bold'}
for index in data[data.adcode.isin(['710000','810000','820000']].index:
    if data.iloc[index]['name'] == "台湾省":
        x = data.iloc[index].geometry.centroid.x
        y = data.iloc[index].geometry.centroid.y
        name = "台湾省"
        ax.text(x, y, name, ha="center", va="center", fontdict=fontdict
    elif data.iloc[index]['name'] == "香港特别行政区":
        x = data.iloc[index].geometry.centroid.x
        y = data.iloc[index].geometry.centroid.y
        name = "香港"
        ax.text(x, y, name, ha="left", va="top", fontdict=fontdict
        gpd.GeoSeries(data.iloc[index].geometry.centroid.plot(ax=ax, facecolor='black', markersize=5
    elif data.iloc[index]['name'] == "澳门特别行政区":
        x = data.iloc[index].geometry.centroid.x
        y = data.iloc[index].geometry.centroid.y
        name = "澳门"
        ax.text(x, y, name, ha="right", va="top", fontdict=fontdict
        gpd.GeoSeries(data.iloc[index].geometry.centroid.plot(ax=ax, facecolor='black', markersize=5

# 移除坐标轴
ax.axis('off'

(15508157.566510478, 20969512.421140186, 92287.20109282521, 6384866.68627745

4 图例设置

# 创建画布matplotlib
fig, ax = plt.subplots(figsize=(12, 9
# 绘制主要区域
ax = data.plot(ax=ax,facecolor='grey',edgecolor='lightgrey',alpha=0.8,linewidth=1
# 绘制九段线
ax = gpd.GeoSeries(nine_dotted_line.geometry.plot(ax=ax,edgecolor='grey',linewidth=3
# 强调首都
ax = data[data.name=="北京市"].representative_point(.plot(ax=ax, facecolor='red',marker='*', markersize=150 

# 强调港澳台
# 设置字体
fontdict = {'family':'FZSongYi-Z13S', 'size':8, 'color': "blue",'weight': 'bold'}
# 这一段代码可能因为不同matplotlib版本出现不同结果
for index in data[data.adcode.isin(['710000','810000','820000']].index:
    if data.iloc[index]['name'] == "台湾省":
        x = data.iloc[index].geometry.centroid.x
        y = data.iloc[index].geometry.centroid.y
        name = "台湾省"
        ax.text(x, y, name, ha="center", va="center", fontdict=fontdict
    elif data.iloc[index]['name'] == "香港特别行政区":
        x = data.iloc[index].geometry.centroid.x
        y = data.iloc[index].geometry.centroid.y
        name = "香港"
        ax.text(x, y, name, ha="left", va="top", fontdict=fontdict
        gpd.GeoSeries(data.iloc[index].geometry.centroid.plot(ax=ax, facecolor='black', markersize=5
    elif data.iloc[index]['name'] == "澳门特别行政区":
        x = data.iloc[index].geometry.centroid.x
        y = data.iloc[index].geometry.centroid.y
        name = "澳门"
        ax.text(x, y, name, ha="right", va="top", fontdict=fontdict
        gpd.GeoSeries(data.iloc[index].geometry.centroid.plot(ax=ax, facecolor='black', markersize=5

# 移除坐标轴
ax.axis('off'

# 单独绘制图例
plt.rcParams["font.family"] = 'FZSongYi-Z13S'
ax.scatter([], [], c='red', s=80,  marker='*', label='首都'
ax.plot([], [], c='grey',linewidth=3, label='南海九段线' 
# 设置图例顺序
handles, labels = ax.get_legend_handles_labels(
ax.legend(handles[::-1], labels[::-1], title="图例",frameon=True, shadow=True, loc="lower left",fontsize=10
# ax.legend(title="图例",frameon=True, shadow=True, loc="lower left",fontsize=10

<matplotlib.legend.Legend at 0x7f159cbe7510>

5 小地图绘制

这种绘图方式非常不专业，也不推荐，建议只是学习使用思路。具体使用看如下代码。

# 创建画布
fig = plt.figure(figsize=(6,3
# 创建一个填充整个画布的子图
ax = fig.add_axes((0,0,1,1
# 从画布宽40%，高40%处绘制子图 
ax_child = fig.add_axes((0.4,0.4,0.2,0.2

from shapely.geometry import Point

# 设中国大陆区域范围和南海范围，估计得到
bound = gpd.GeoDataFrame({
    'x': [80, 140, 106.5, 123],
    'y': [15, 50, 2.8, 24.5]
}
bound.geometry = bound.apply(lambda row: Point([row['x'], row['y']], axis=1
# 初始化CRS
bound.crs = 'EPSG:4326'
bound = bound.to_crs('EPSG:4573'
bound

# 创建画布matplotlib
fig, ax = plt.subplots(figsize=(12, 9
# 绘制主要区域
ax = data.plot(ax=ax,facecolor='grey',edgecolor='lightgrey',alpha=0.8,linewidth=1
# 绘制九段线
ax = gpd.GeoSeries(nine_dotted_line.geometry.plot(ax=ax,edgecolor='grey',linewidth=3
# 强调首都
ax = data[data.name=="北京市"].representative_point(.plot(ax=ax, facecolor='red',marker='*', markersize=150 

# 强调港澳台
# 设置字体
fontdict = {'family':'FZSongYi-Z13S', 'size':8, 'color': "blue",'weight': 'bold'}
# 这一段代码可能因为不同matplotlib版本出现不同结果
for index in data[data.adcode.isin(['710000','810000','820000']].index:
    if data.iloc[index]['name'] == "台湾省":
        x = data.iloc[index].geometry.centroid.x
        y = data.iloc[index].geometry.centroid.y
        name = "台湾省"
        ax.text(x, y, name, ha="center", va="center", fontdict=fontdict
    elif data.iloc[index]['name'] == "香港特别行政区":
        x = data.iloc[index].geometry.centroid.x
        y = data.iloc[index].geometry.centroid.y
        name = "香港"
        ax.text(x, y, name, ha="left", va="top", fontdict=fontdict
        gpd.GeoSeries(data.iloc[index].geometry.centroid.plot(ax=ax, facecolor='black', markersize=5
    elif data.iloc[index]['name'] == "澳门特别行政区":
        x = data.iloc[index].geometry.centroid.x
        y = data.iloc[index].geometry.centroid.y
        name = "澳门"
        ax.text(x, y, name, ha="right", va="top", fontdict=fontdict
        gpd.GeoSeries(data.iloc[index].geometry.centroid.plot(ax=ax, facecolor='black', markersize=5

# 设置大陆区域范围
ax.set_xlim(bound.geometry[0].x, bound.geometry[1].x
ax.set_ylim(bound.geometry[0].y, bound.geometry[1].y

# 移除坐标轴
ax.axis('off'

# 单独绘制图例
plt.rcParams["font.family"] = 'FZSongYi-Z13S'
ax.scatter([], [], c='red', s=80,  marker='*', label='首都'
ax.plot([], [], c='grey',linewidth=3, label='南海九段线' 
# 设置图例顺序
handles, labels = ax.get_legend_handles_labels(
ax.legend(handles[::-1], labels[::-1], title="图例",frameon=True, shadow=True, loc="lower left",fontsize=10
# ax.legend(title="图例",frameon=True, shadow=True, loc="lower left",fontsize=10

# 创建南海插图对应的子图，调整这些参数以调整地图位置
ax_child = fig.add_axes([0.75, 0.15, 0.2, 0.2]
ax_child = data.plot(ax=ax_child,facecolor='grey',edgecolor='lightgrey',alpha=0.8,linewidth=1
# 绘制九段线
ax_child = gpd.GeoSeries(nine_dotted_line.geometry.plot(ax=ax_child,edgecolor='grey',linewidth=3

# 设置子图显示范围
ax_child.set_xlim(bound.geometry[2].x, bound.geometry[3].x
ax_child.set_ylim(bound.geometry[2].y, bound.geometry[3].y
ax_child.text(0.98,0.02,'Produced by luohenyueji',transform = ax.transAxes,
        ha='center', va='center',fontsize = 12,color='black'
# 移除子图坐标轴
ax_child.set_xticks([]
ax_child.set_yticks([]

fig.savefig('res.png', dpi=300, bbox_inches='tight'

2 分层设色

2.1 分层设色基本介绍

# 读取2019江苏省各市GDP数据

import geopandas as gpd
import matplotlib.pyplot as plt
import pandas as pd
plt.rcParams["font.family"] = 'FZSongYi-Z13S'
# 数据来自互联网
gdp = pd.read_csv("2022江苏省各市GDP.csv"
gdp

# 读取江苏地图数据，数据来自DataV.GeoAtlas，将其投影到EPSG:4573
data = gpd.read_file('https://geo.datav.aliyun.com/areas_v3/bound/320000_full.json'.to_crs('EPSG:4573'
# 合并数据
data = data.join(gdp.set_index('地级市'["2022年GDP（亿元）"],on='name'
# 修改列名
data.rename(columns={'2022年GDP（亿元）':'GDP'},inplace=True
data.head(

fig, ax = plt.subplots(figsize=(9, 9

# legend_kwds设置matplotlib的legend参数
data.plot(ax=ax,column='GDP', cmap='coolwarm', legend=True,legend_kwds={'label': "GDP（亿元）", 'shrink':0.5}
ax.axis('off'

# 设置
fontdict = {'family':'FZSongYi-Z13S', 'size':8, 'color': "black",'weight': 'bold'}

# 设置标题
ax.set_title('江苏省地级市2022年GDP数据可视化', fontsize=24 
for index in data.index:
    x = data.iloc[index].geometry.centroid.x
    y = data.iloc[index].geometry.centroid.y
    name = data.iloc[index]["name"]
    if name in ["苏州市","无锡市"]:
        x = x*1.001
    ax.text(x, y, name, ha="center", va="center", fontdict=fontdict
# 保存图片
fig.savefig('res.png', dpi=300, bbox_inches='tight'

在本文通过Python模块mapclassify用于分层设色和数据可视化。使用mapclassify之前需要输入以下命令安装相关模块：

mapclassify官方仓库见：mapclassify。mapclassify提供了多种分组方法，可以帮助我们更好地理解数据的分布情况，mapclassify提供的方法包括：

BoxPlot: 基于箱线图的分类方法。这种分类方法适用于数据分布比较规律的情况。
• EqualInterval: 等距离分类方法。这种分类方法将数据划分为等距离的若干区间。适用于数据分布比较均匀的情况。
• FisherJenks: 基于Fisher-Jenks算法的分类方法。这种分类方法将数据划分为若干区间，使得每个区间内部的差异最小，不同区间之间的差异最大。适用于数据分布比较不规律的情况。
• HeadTailBreaks: 基于Head-Tail算法的分类方法。这种分类方法将给定的数据集分为两部分：头部和尾部。头部通常包含出现频率最高的值，而尾部包含出现频率较低的值。适用于识别数据集中的异常值和离群值。
• JenksCaspall: 基于Jenks-Caspall算法的分类方法。这种分类方法根据数据中发现的自然分组将数据集划分为类。适用于需要将数据分类为几个具有明显含义的区间的情况。
• JenksCaspallForced: 强制基于Jenks-Caspall算法的分类方法。与JenksCaspall算法类似，但是它对区间的数量和大小有更强的控制力。适用于需要精确控制区间数量和大小的情况。
• JenksCaspallSampled: 采样基于Jenks-Caspall算法的分类方法。该方法对数据进行采样，然后使用Jenks-Caspall算法对采样后的数据进行分类，适用于数据量比较大的情况。
• MaxP: 基于最大界限的分类方法。这种分类方法将数据划分为几个区间，使得不同区间之间的差异最大。适用于需要将数据分类为几个具有明显差异的区间的情况。
• MaximumBreaks: 基于最大间隔的分类方法。这种分类方法与MaxP算法类似，但是它更加注重区间的可理解性。适用于需要将数据分类为几个具有明显含义的区间的情况。
• NaturalBreaks: 基于自然间隔的分类方法。这种分类方法将数据划分为几个区间，使得每个区间内部的差异最小，不同区间之间的差异最大。适用于数据分布比较不规律的情况
• Quantiles: 基于分位数的分类方法。
• Percentiles: 基于百分位数的分类方法。
• StdMean: 基于标准差分组的分类方法。
• UserDefined: 基于自定义分组的分类方法。

示例1

import mapclassify
# 导入示例数据
y = mapclassify.load_example(
print(type(y
print(y.mean(,y.min(, y.max(
y.head(

<class 'pandas.core.series.Series'>
125.92810344827588 0.13 4111.45





0    329.92
1      0.42
2      5.90
3     14.03
4      2.78
Name: emp/sq km, dtype: float64

mapclassify.EqualInterval(y

EqualInterval

     Interval        Count
--------------------------
[   0.13,  822.39] |    57
( 822.39, 1644.66] |     0
(1644.66, 2466.92] |     0
(2466.92, 3289.19] |     0
(3289.19, 4111.45] |     1

示例2

y = [1,2,3,4,5,6,7,8,9,0]
# 分为四个区间
mapclassify.JenksCaspall(y, k=4

JenksCaspall

  Interval     Count
--------------------
[0.00, 2.00] |     3
(2.00, 4.00] |     2
(4.00, 6.00] |     2
(6.00, 9.00] |     3

示例3

y = [1,2,3,4,5,6,7,8,9,0]
# 自定义区间
mapclassify.UserDefined(y, bins=[5, 8, 9]

UserDefined

  Interval     Count
--------------------
[0.00, 5.00] |     6
(5.00, 8.00] |     3
(8.00, 9.00] |     1

示例4

import mapclassify 
import pandas
from numpy import linspace as lsp
demo = [lsp(3,8,num=10, lsp(10, 0, num=10, lsp(-5, 15, num=10]
demo = pandas.DataFrame(demo.T
demo.head(

# 使用apply函数应用分层，rolling表示是否进行滑动窗口计算以消除随机波动
demo.apply(mapclassify.Quantiles.make(rolling=True.head(

2.2 绘图实例之用于地图的分层设色

方法1

fig, ax = plt.subplots(figsize=(9, 9

# 使用分层设色后,legend_kwds要进行相应修改
data.plot(ax=ax,column='GDP', cmap='coolwarm', legend=True,scheme='JenksCaspall',k=4,legend_kwds={
                                                     'loc': 'lower left',
                                                     'title': 'GDP数据分级（亿元）',
                                                 }
ax.axis('off'

# 设置
fontdict = {'family':'FZSongYi-Z13S', 'size':8, 'color': "black",'weight': 'bold'}

# 设置标题
ax.set_title('江苏省地级市2022年GDP数据可视化', fontsize=24 
for index in data.index:
    x = data.iloc[index].geometry.centroid.x
    y = data.iloc[index].geometry.centroid.y
    name = data.iloc[index]["name"]
    if name in ["苏州市","无锡市"]:
        x = x*1.001
    ax.text(x, y, name, ha="center", va="center", fontdict=fontdict
# 保存图片
fig.savefig('res.png', dpi=300, bbox_inches='tight'

方法2

# 创建分类器
classifier = mapclassify.HeadTailBreaks(data['GDP']
classifier

HeadTailBreaks

      Interval         Count
----------------------------
[ 4005.00,  9473.76] |     8
( 9473.76, 15329.34] |     3
(15329.34, 20433.10] |     1
(20433.10, 23958.30] |     1

# 赋值数据
data['GDP_class'] = data['GDP'].apply(classifier
data['GDP_class'] = data['GDP_class'].apply(lambda x : int(x
data.head(

import matplotlib.colors as colors
fig, ax = plt.subplots(figsize=(9, 9

# 设置分层颜色条
cmap = plt.cm.get_cmap('coolwarm', len(set(data['GDP_class']
# vmax和vmin设置是为了让等级值居中
data.plot(ax=ax,column='GDP_class', cmap=cmap, legend=False,vmin=-0.5,vmax=3.5
ax.axis('off'
# 设置Colorbar的刻度
cbar = ax.get_figure(.colorbar(ax.collections[0],shrink=0.5
cbar.set_ticks([0,1,2,3]
cbar.set_label('GDP数据分级'
cbar.set_ticklabels(['等级0','等级1','等级2','等级3']
# 隐藏刻度线
ticks = cbar.ax.get_yaxis(.get_major_ticks(
for tick in ticks:
    tick.tick1line.set_visible(False
    tick.tick2line.set_visible(False


# 设置
fontdict = {'family':'FZSongYi-Z13S', 'size':8, 'color': "black",'weight': 'bold'}

# 设置标题
ax.set_title('江苏省地级市2022年GDP数据可视化', fontsize=24 
for index in data.index:
    x = data.iloc[index].geometry.centroid.x
    y = data.iloc[index].geometry.centroid.y
    name = data.iloc[index]["name"]
    if name in ["苏州市","无锡市"]:
        x = x*1.001
    ax.text(x, y, name, ha="center", va="center", fontdict=fontdict
# 保存图片
fig.savefig('res.png', dpi=300, bbox_inches='tight'

方法3

# 创建分类器
classifier = mapclassify.MaximumBreaks(data['GDP'], k=3
classifier

MaximumBreaks

      Interval         Count
----------------------------
[ 4005.00, 13115.20] |    10
(13115.20, 20433.10] |     2
(20433.10, 23958.30] |     1

# 赋值数据
data['GDP_class'] = data['GDP'].apply(classifier
data['GDP_class'] = data['GDP_class'].apply(lambda x : int(x
data.head(

import matplotlib.patches as mpatches

# 设置图案列表
patterns = ["///", "",  "*", "\\\\",".", "o",  "O",]
cmap = plt.cm.get_cmap('coolwarm', len(set(data['GDP_class']
color_list = cmap([0,1,2]
fig, ax = plt.subplots(figsize=(9, 9
# 自定义图示
legend_list = []

# 按层次设置legend
for i in set(data['GDP_class']:
    tmp = data[data['GDP_class']==i]
    tmp.plot(ax=ax,column='GDP_class', legend=False,hatch=patterns[i],edgecolor='black',color=color_list[i], linestyle='-',linewidth=2
    legend_list.append(
        mpatches.Patch(facecolor=color_list[i], edgecolor='black',linestyle='-', linewidth=2,hatch=patterns[i], label='等级{}'.format(i
    
ax.axis('off'
# 设置标题
ax.set_title('江苏省地级市2022年GDP数据可视化', fontsize=24 

# 自定义图示
ax.legend(handles = legend_list, loc='best', fontsize=12, title='GDP数据分级', shadow=True

# 保存图片
fig.savefig('res.png', dpi=300, bbox_inches='tight'

3 参考

GeoPandas
• GeoPandas-doc
• [数据分析与可视化] Python绘制数据地图1-GeoPandas入门指北
• Chart visualization
• 基于geopandas的空间数据分析之基础可视化
• geopandas 中国地图绘制
• mapclassify
• 基于geopandas的空间数据分析——深入浅出分层设色

声明：本文由用户投稿上传，本站不保证内容的真实与正确性，并且不承担相关法律责任，如有侵权请提供版权资料并联系删除！

编程笔记 » [数据分析与可视化] Python绘制数据地图2-GeoPandas地图可视化