如何在Python中为伊朗TIF地图添加流域边界Shapefile

要在TIF地图上叠加清晰的流域边界Shapefile，你可以借助geopandas库读取并绘制Shapefile，然后将其叠加到现有的matplotlib轴上。以下是具体实现步骤：

步骤1：安装依赖 #

如果还未安装geopandas，先执行下方安装命令：

pip install geopandas matplotlib

步骤2：修改代码叠加流域边界 #

from pyidw import idw
from matplotlib import pyplot as plt
import geopandas as gpd

# 生成TIF底图
fig, ax, cbar = idw.show_map(
    input_raster="Anomaly (2051-2075)_SSP126.tif",
    colormap="brg",
    image_size=1.5,
    return_figure=True)

# 保留你原有的颜色范围调整逻辑
image = ax.get_images()[0]
image.set_clim(10, -7)
cbar.update_normal(image)

# 读取流域边界Shapefile
basin_shp = gpd.read_file("你的流域边界文件路径.shp")

# 关键：确保坐标系匹配（若TIF与ShapefileCRS不同，必须转换）
# 示例：若TIF用WGS84（EPSG:4326），将Shapefile转换为对应CRS
# basin_shp = basin_shp.to_crs(epsg=4326)

# 绘制清晰的流域边界：通过参数控制样式，确保边界在底图之上
basin_shp.boundary.plot(ax=ax, linewidth=1.8, color="#000000", zorder=2)

# 设置地图标注
ax.set_title("Maximum temperature map")
ax.set_xlabel("Longitude")
ax.set_ylabel("Latitude")

plt.show()

核心细节说明 #

• 坐标系对齐：如果TIF和Shapefile的坐标参考系统（CRS）不一致，必须用to_crs()方法统一坐标系，否则边界会出现错位。可通过basin_shp.crs查看Shapefile的CRS，TIF的CRS可通过GDAL工具或读取TIF元数据获取。
• 边界样式优化：linewidth控制线条粗细，color设置边界颜色（黑色是最通用的醒目选择），zorder=2确保边界图层覆盖在TIF底图之上（数值越大层级越高），这是让边界清晰可见的关键。
• 路径替换：将代码中的"你的流域边界文件路径.shp"替换为实际的Shapefile文件路径。

内容的提问来源于stack exchange，提问作者mahyar