HF.092 | 请把你画好的全球要素分布图甩给导师，然后等待夸奖！

#加载包import rasterioimport geopandas as gpdimport matplotlib.pyplot as pltfrom rasterio.warp import calculate_default_transform, reproject, Resamplingimport cartopy.crs as ccrsimport cartopy.feature as cfeatureimport numpy as npfrom matplotlib.font_manager import FontProperties
# 文件路径raster_path = "C:/Users/ASUS/Desktop/t2m/t2m_202001.tif"shapefile_path = "C:/Users/ASUS/Desktop/t2m/World_countries.shp"
# 加载 shapefileworld = gpd.read_file(shapefile_path)
# 设置 Equal Earth 投影equal_earth_crs = ccrs.EqualEarth()
# 打开栅格文件with rasterio.open(raster_path) as src:    # 计算变换并重投影    transform, width, height = calculate_default_transform(        src.crs, equal_earth_crs, src.width, src.height, *src.bounds)        kwargs = src.meta.copy()    kwargs.update({        'crs': equal_earth_crs,        'transform': transform,        'width': width,        'height': height    })        with rasterio.open('reprojected.tif', 'w', **kwargs) as dst:        for i in range(1, src.count + 1):            reproject(                source=rasterio.band(src, i),                destination=rasterio.band(dst, i),                src_transform=src.transform,                src_crs=src.crs,                dst_transform=transform,                dst_crs=equal_earth_crs,                resampling=Resampling.nearest)
# 绘制fig, ax = plt.subplots(figsize=(12, 8), subplot_kw={'projection': equal_earth_crs})ax.set_global()
# 打开重投影后的栅格文件并使用 imshow 绘制with rasterio.open('reprojected.tif') as src:    data = src.read(1)    extent = [src.bounds.left, src.bounds.right, src.bounds.bottom, src.bounds.top]    img = ax.imshow(data, cmap='RdPu', extent=extent, transform=equal_earth_crs, vmin=220, vmax=320)
# 添加 shapefileworld.boundary.plot(ax=ax, transform=ccrs.PlateCarree(), linewidth=0.5, edgecolor='black')
# 设置字体以支持中文