本文将介绍通过Python绘制地形图的方法,所需第三方Python相关模块包括 rasterio、geopandas、cartopy 等,可通过 pip 等方式安装。
本文将介绍通过Python绘制地形图的方法,所需第三方Python相关模块包括 rasterio、geopandas、cartopy 等,可通过 pip 等方式安装。
1 示例代码
1.1 导入相关模块
| import rasterio |
| import geopandas as gpd |
| import numpy as np |
| import cartopy.crs as ccrs |
| import matplotlib.pyplot as plt |
| from matplotlib.colors import LinearSegmentedColormap |
1.2 设置地图字体及字号
| plt.rcParams["font.family"] = "Times New Roman" |
| plt.rcParams["font.size"] = 14 |
1.3 绘制地图
| |
| |
| projection = ccrs.Mercator() |
| |
| |
| fig, ax = plt.subplots(figsize=(20, 10), subplot_kw={'projection': projection}) |
| |
| |
| ax.set_extent([13530000, 14630000, 4960000, 5850000], crs=projection) |
| |
| |
| shp = gpd.read_file("Data/Jilin_Mercator.shp") |
| shp.plot(ax=ax, transform=projection, edgecolor="black", linewidth=1, facecolor="none") |
| |
| |
| colors = ["#369121", "#95C769", "#FFFFBF", "#E6865A", "#D14E30", "#BA1414"] |
| n_bins = 100 |
| cmap_name = "green_brown" |
| cm = LinearSegmentedColormap.from_list(cmap_name, colors, N=n_bins) |
| |
| |
| dataset = rasterio.open("Data/DEM_Jilin_Mercator.tif") |
| data = dataset.read(1) |
| nodata_value = dataset.nodata |
| |
| data = np.ma.masked_where(data == nodata_value, data) |
| |
| |
| extent = [dataset.bounds.left, dataset.bounds.right, |
| dataset.bounds.bottom, dataset.bounds.top] |
| im = ax.imshow(data, origin="upper", extent=extent, |
| transform=projection, cmap=cm) |
| |
| |
| gl = ax.gridlines(draw_labels=True, linestyle="--", color="#4F4F4F") |
| gl.xlocator = plt.FixedLocator(range(120, 135, 3)) |
| gl.ylocator = plt.FixedLocator(range(40, 50, 2)) |
| |
| |
| cbar = plt.colorbar(im, ax=ax, orientation="horizontal", shrink=0.4, pad=0.1) |
| cbar.set_label("Elevation (m)", labelpad=10) |
| cbar.ax.xaxis.set_label_position("top") |
| |
| plt.savefig("Pic.jpg", dpi=600) |
| plt.show() |
2 结果图展示

【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 10亿数据,如何做迁移?
· 推荐几款开源且免费的 .NET MAUI 组件库
· 清华大学推出第四讲使用 DeepSeek + DeepResearch 让科研像聊天一样简单!
· c# 半导体/led行业 晶圆片WaferMap实现 map图实现入门篇
· 易语言 —— 开山篇