【可视化】绘制中国标准行政区划地图（2024版）

晨昏三变

于 2024-08-27 16:26:16 发布

阅读量1.2k

点赞数 7

分类专栏：可视化文章标签： python 信息可视化

本文链接：https://blog.csdn.net/m0_38065162/article/details/141606151

版权

可视化专栏收录该内容

8 篇文章 1 订阅

订阅专栏

【可视化】绘制中国标准行政区划地图（2024版）

在科研投稿论文过程中，多数期刊杂志要求我们制作涉及国界的图件，必须使用自然资源部地图技术审查中心发布的标准地图底图。凡涉国界图件（国内部分地区、全国、世界部分地区、全球）必须使用
自然资源部标准地图服务系统或国家测绘局地理信息局标准地图底图，所用底图边界要完全无修改（包括南海诸岛位置）。本文结合Geopandas和Cartopy两个强大的地理数据处理库进行展示。

代码实现

导入相关库

# -*- coding: UTF-8 -*-
# 导入相关库
import numpy as np
import xarray as xr
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
from copy import copy
from cartopy.mpl.gridliner import LATITUDE_FORMATTER, LONGITUDE_FORMATTER
import shapely.geometry as sgeom

加载数据

**国家地理信息公共服务平台（天地图）**新增了矢量地图底图下载服务，用户可以根据需求下载不同审图号的标准地图矢量底图。矢量地图相较于像素地图，具有更高的清晰度和可编辑性，可以无限放大而不失真，更适合进行地图编辑和二次创作。获取方式见-天地图。

读取中国区域的省、市、县的区划信息：

# 读取省级行政区划
province = gpd.read_file('中国_省.geojson')
# 读取市级行政区划
city = gpd.read_file('中国_市.geojson')
# 读取县级行政区划
county = gpd.read_file('中国_县.geojson')

市级行政区划展示

# 设置地图投影
proj = ccrs.LambertConformal(central_longitude=105, central_latitude=90,standard_parallels=(25, 47))

# 设置图片大小
fig = plt.figure(figsize=[10, 8],frameon=True)

# Set projection and plot the main figure
ax = fig.add_axes([0.08, 0.05, 0.8, 0.94], projection=proj)
# 设置地图范围
ax.set_extent([80, 130, 15, 55],crs=ccrs.PlateCarree())
# 添加地理要素
ax.add_geometries(city.geometry, 
                  crs=ccrs.PlateCarree(), 
                  edgecolor='k',
                  facecolor='none',
                  zorder=1)

# 添加海洋、陆地、河流、湖泊
ax.add_feature(cfeature.OCEAN.with_scale('50m'),zorder=0)
ax.add_feature(cfeature.LAND.with_scale('50m'),zorder=0)
ax.add_feature(cfeature.RIVERS.with_scale('50m'),zorder=0)
ax.add_feature(cfeature.LAKES.with_scale('50m'),zorder=0)

# 设置并添加经纬度标签
xticks = [55, 65, 75, 85, 95, 105, 115, 125, 135, 145, 155, 165]
yticks = [0 , 5 , 10, 15, 20, 25 , 30 , 35 , 40 , 45 , 50 , 55 , 60 , 65]

gl = ax.gridlines(xlocs=xticks, 
                  ylocs=yticks,
                  draw_labels=True,
                  x_inline=False, 
                  y_inline=False,
                  linewidth=1,
                  color='dimgrey',
                  alpha=0.5,
                  linestyle='--')

gl.top_labels = False 
gl.right_labels = False 
gl.xformatter = LONGITUDE_FORMATTER  
gl.yformatter = LATITUDE_FORMATTER  

# 添加南海附图
sub_ax = fig.add_axes([0.75, 0.102, 0.15, 0.2],
                      projection=ccrs.LambertConformal(central_latitude=90,
                                                       central_longitude=105))

# 添加海洋、陆地、河流、湖泊
sub_ax.add_feature(cfeature.OCEAN.with_scale('50m'))
sub_ax.add_feature(cfeature.LAND.with_scale('50m'))
sub_ax.add_feature(cfeature.RIVERS.with_scale('50m'))
sub_ax.add_feature(cfeature.LAKES.with_scale('50m'))

# 设置南海附图范围
sub_ax.set_extent([106, 120, 2, 24], crs=ccrs.PlateCarree())
sub_ax.add_geometries(city.geometry, crs=ccrs.PlateCarree(), edgecolor='k', facecolor='none')

plt.show()