Python

2024-07-17 00:38| 来源: 网络整理| 查看: 265

看过很多大佬绘制过中国地图，有R-ggplot的，也有python-geopandas等，还有其他cartopy等的作图包，我就不一一举例了，今天用proplot这个包，这个包的简介可以参考链接proplot，主要的知识点如下：

proplot绘制中国地图 ax.inset添加南海子图导入需要的包 import os import numpy as np import xarray as xr import proplot as plot import cartopy.crs as ccrs from cartopy.io.shapereader import Reader import maskout shp_path = 'D:/China_shp' 导入数据以历史气温数据为例（historical） ##### 数据读取 ds = xr.open_dataset('*historical*.nc') tas = ds['tas'][75] # 随意选择的年份 #截取一段时间内的降水作为绘图的数据源 tas -= 273.15 lon = ds.lon lat = ds.lat lon_range = lon[(lon > 70) & (lon < 140)] lat_range = lat[(lat > 0) & (lat < 60)] tas_region = tas.sel(lon=lon_range, lat=lat_range) levels = plot.arange(-5,30,5) 主图&白化[1] fig, axs = plot.subplots(width=6,proj=ccrs.PlateCarree()) # First axis ax = axs[0] m = ax.contourf(tas_region, cmap='CoolWarm', levels=levels, extend='both') # 搜白化教程 # http://bbs.06climate.com/forum.php?mod=viewthread&tid=42437&highlight=%CD%EA%C3%C0%B0%D7%BB%AF con_mask = maskout.shp2clip(m, ax, r'E:\maskout\country1.shp', ['China']) ax.add_geometries(Reader(os.path.join(shp_path, 'bou2_4l.shp')).geometries(), ccrs.PlateCarree(), facecolor='none', edgecolor='black', linewidth=0.5) ax.format( coast=False, labels=True, borders=False, latlim=(10, 60), lonlim=(70, 140), lonlines=10, latlines=10, )

[1]地图白化

添加南海子图 # Create inset axes representing a "zoom-in" # https://proplot.readthedocs.io/en/latest/insets_panels.html iax = ax.inset([0.83,0,0.175,0.3]) # 不同的主图参数，这个参数设置也不一样，自己摸索 n = iax.contourf(tas_region, cmap='CoolWarm', levels=levels, extend='both') con_mask = maskout.shp2clip(n, iax, r'E:\maskout\country1.shp', ['China']) iax.add_geometries(Reader(os.path.join(shp_path, 'bou2_4l.shp')).geometries(), ccrs.PlateCarree(), facecolor='none', edgecolor='black', linewidth=0.5) iax.format( coast=False, labels=False, borders=False, latlim=(0, 26), lonlim=(105, 125), lonlines=10, latlines=10, ) fig.colorbar(m, loc='r') 出图 China 总结

Proplot这个包太强大了，在这里就不细说了！共同学习，那就这样吧！(时间紧，任务重/奋斗脸)，写这个就当是抛砖引玉吧！希望大家多多指教！

【本文地址】

公司简介

联系我们