Python绘图：绘制地形图

jack_ln

1、Python绘图：使用Basemap绘制全球地形图_菩提有树的博客-CSDN博客

                               
登录/注册后可看大图

2、Python绘制中南半岛与中国南海DEM地形图 (qq.com)
参考了文章1中的地形配色，综合以上链接，代码如下。

1. # -*- encoding: utf-8 -*-
   
2. import numpy as np
   
3. import matplotlib.pyplot as plt
   
4. import cartopy.crs as ccrs
   
5. import cartopy.feature as cfeature
   
6. import cartopy.mpl.ticker as cticker
   
7. import xarray as xr
   
8. from matplotlib import rcParams
   
9. from cartopy.io.shapereader import Reader
   
10. 
    
11. config = {"font.family":'Times New Roman',"font.size": 12,"mathtext.fontset":'stix'}
    
12. rcParams.update(config)
    
13. filename=r'ETOPO2v2c_f4.nc'
    
14. proj=ccrs.PlateCarree()
    
15. #########设置投影及经纬度范围#####################################################
    
16. extent = [70, 140, 10, 60]
    
17. lonmin, lonmax, latmin, latmax = extent
    
18. fig = plt.figure(figsize=(9, 8),dpi=300)
    
19. ax = fig.add_axes([0.1, 0.1, 0.8, 0.6], projection=ccrs.PlateCarree())
    
20. ax.set_extent(extent, crs=ccrs.PlateCarree())
    
21. ax.set_xticks(np.arange(lonmin,lonmax+10,10), crs=ccrs.PlateCarree())
    
22. ax.set_yticks(np.arange(latmin,latmax+5,5), crs=ccrs.PlateCarree())
    
23. lon_formatter = cticker.LongitudeFormatter()
    
24. lat_formatter = cticker.LatitudeFormatter()
    
25. ax.xaxis.set_major_formatter(lon_formatter)
    
26. ax.yaxis.set_major_formatter(lat_formatter)
    
27. ##########读取.NC格式的高程数据##################################################
    
28. f=xr.open_dataset(filename)
    
29. lon=f['x'][:]
    
30. lat=f['y'][:]
    
31. height=f['z'][:]
    
32. 
    
33. levels = [-8000, -6000, -4000, -2000, -1000, -200, -50, 0, 50, 200, 500, 1000, 1500, 2000, 3000, 4000,
    
34.         5000, 6000, 7000, 8000]  # 创建分级
    
35. color = ['#084594', '#2171b5', '#4292c6', '#6baed6', '#9ecae1', '#c6dbef', '#deebf7', '#006837',
    
36.         '#31a354', '#78c679', '#addd8e', '#d9f0a3', '#f7fcb9', '#c9bc87', '#a69165', '#856b49',
    
37.         '#664830', '#ad9591', '#d7ccca']  # 设置色带
    
38. 
    
39. cf=ax.contourf(lon,lat,height,levels=levels,colors=color,extend='neither')
    
40. font3={'family':'SimHei','size':12,'color':'k'}
    
41. ax.add_feature(cfeature.OCEAN)
    
42. ##########设置colorbar属性######################################################
    
43. position=fig.add_axes([0.91, 0.12, 0.025,0.46])
    
44. cb=fig.colorbar(cf,cax=position)
    
45. ax_1=cb.ax
    
46. ax_1.set_title('海拔',fontdict=font3)
    
47. ax_1.tick_params(which='major',direction='in',labelsize=10,length=3.5)
    
48. ################加载shp文件######################################################
    
49. ax.add_feature(cfeature.LAKES.with_scale('50m'))
    
50. ax.add_geometries(Reader(r'country1.shp').geometries(),ccrs.PlateCarree(),facecolor='none',edgecolor='k',linewidth=1)
    
51. ax.add_geometries(Reader(r'province.shp').geometries(),ccrs.PlateCarree(),facecolor='none',edgecolor='k',linewidth=1)
    
52. plt.savefig('china.png',dpi=300,bbox_inches='tight',pad_inches=0)

复制代码

                               
登录/注册后可看大图

jack_ln

guoyumin 发表于 2023-3-23 09:09
你好 这个数据可以共享下吗filename=r'ETOPO2v2c_f4.nc'

ETOPO2v2c_f4.nc下载地址：https://www.ngdc.noaa.gov/mgg/gl ... O2v2c_f4_netCDF.zip

river

学习了，好看！

dengfenglai573

好看！赞一个~~~

雨落森林

这个好！！

likailing2008

这配色配的是真的好看

对方的凤飞飞

不错

努力的叮宝

好看{:eb502:}{:eb502:}

guoyumin

你好 这个数据可以共享下吗filename=r'ETOPO2v2c_f4.nc'

翘课小公主

好看！博主太牛啦！感谢！