<think>我们将绘制一张包含整个中国区域的地图，并在图中用红色虚线框标出四个子区域，同时在每个子区域内显示该区域的平均升温趋势数值。 步骤： 1. 计算整个中国区域每个格点的升温趋势（与之前相同）。 2. 计算四个子区域的平均升温趋势（区域平均，然后计算趋势）。 3. 绘制整个中国区域的升温趋势空间分布图。 4. 在图上用红色虚线框标出四个子区域，并在每个框内显示该区域的平均升温趋势（数值，单位℃/decade）。 5. 添加标题、色标等。 注意：由于需要在一张图上显示四个子区域的位置和各自的平均趋势，我们将： - 先计算每个子区域的时间序列（每年一个值，区域平均），然后计算每个子区域的趋势（一个数值）。 - 同时，计算整个中国区域的每个格点的趋势（用于填色图）。 我们将分步骤进行： </think> 以下是一个完整的Python代码，用于计算并绘制中国四个子区域的面积平均升温趋势空间分布图，并在一张图中展示所有区域： ```python import xarray as xr import numpy as np import matplotlib.pyplot as plt import cartopy.crs as ccrs import cartopy.feature as cfeature from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER from scipy import stats import matplotlib.patches as mpatches from matplotlib.colors import LinearSegmentedColormap # 1. 定义区域边界 china_bbox = {'lon_min': 73, 'lon_max': 135, 'lat_min': 18, 'lat_max': 54} regions = { "北方地区": {"lat_min": 34.0, "lat_max": 54.0, "lon_min": 105.0, "lon_max": 135.0}, "西北地区": {"lat_min": 36.0, "lat_max": 50.0, "lon_min": 73.0, "lon_max": 105.0}, "南方地区": {"lat_min": 18.0, "lat_max": 34.0, "lon_min": 105.0, "lon_max": 135.0}, "青藏高原": {"lat_min": 26.0, "lat_max": 36.0, "lon_min": 73.0, "lon_max": 105.0} } # 2. 读取所有年份数据并计算年平均 years = range(1980, 2015) annual_means = [] print("正在处理数据...") for year in years: file_path = f"/media/HUAWEI/LY-TSD302/ChinaMet Tmpmean/ChinaMet_0.25deg(365_fixed)/ChinaMet_025deg_tmpmean_{year}.nc" ds = xr.open_dataset(file_path) # 计算年平均温度并裁剪到中国区域 annual_mean = ds['tmpmean'].mean(dim='time').sel( lat=slice(china_bbox['lat_min'], china_bbox['lat_max']), lon=slice(china_bbox['lon_min'], china_bbox['lon_max']) ) annual_means.append(annual_mean) ds.close() # 合并所有年份数据 ds_annual = xr.concat(annual_means, dim='year') ds_annual['year'] = years # 3. 计算每个子区域的平均升温趋势 print("计算各区域平均升温趋势...") regional_trends = {} for region, bbox in regions.items(): # 提取区域数据 region_data = ds_annual.sel( lat=slice(bbox['lat_min'], bbox['lat_max']), lon=slice(bbox['lon_min'], bbox['lon_max']) ) # 计算区域平均时间序列 regional_series = region_data.mean(dim=['lat', 'lon']) # 计算线性趋势 (℃/decade) slope, _, _, _, _ = stats.linregress(years, regional_series) regional_trends[region] = slope * 10 # 4. 计算整个中国区域的升温趋势图 print("计算格点升温趋势...") def calc_trend(timeseries): if np.isnan(timeseries).any(): return np.nan slope, _, _, _, _ = stats.linregress(years, timeseries) return slope * 10 # ℃/decade trend_data = xr.apply_ufunc( calc_trend, ds_annual, input_core_dims=[['year']], vectorize=True, output_dtypes=[np.float64] ) # 5. 创建综合地图 print("绘制综合空间分布图...") plt.figure(figsize=(14, 12)) ax = plt.axes(projection=ccrs.PlateCarree()) # 设置地图范围 ax.set_extent([china_bbox['lon_min'], china_bbox['lon_max'], china_bbox['lat_min'], china_bbox['lat_max']], crs=ccrs.PlateCarree()) # 添加地理特征 ax.add_feature(cfeature.COASTLINE, linewidth=0.8) ax.add_feature(cfeature.BORDERS, linestyle=':', linewidth=0.8) ax.add_feature(cfeature.RIVERS, linewidth=0.5, alpha=0.5) ax.add_feature(cfeature.LAKES, alpha=0.5) ax.add_feature(cfeature.OCEAN, alpha=0.2) ax.add_feature(cfeature.LAND, alpha=0.1) # 创建自定义色标 (蓝-白-红) colors = ["#313695", "#4575b4", "#74add1", "#abd9e9", "#e0f3f8", "#ffffbf", "#fee090", "#fdae61", "#f46d43", "#d73027", "#a50026"] cmap = LinearSegmentedColormap.from_list("custom_cmap", colors, N=256) # 绘制中国区域升温趋势 vmin = -0.5 # ℃/decade vmax = 0.8 # ℃/decade trend_plot = trend_data.plot.imshow(ax=ax, transform=ccrs.PlateCarree(), cmap=cmap, vmin=vmin, vmax=vmax, add_colorbar=False, alpha=0.8) # 添加色标 cbar = plt.colorbar(trend_plot, ax=ax, orientation='vertical', pad=0.05, aspect=30, shrink=0.7) cbar.set_label('温度变化趋势 (℃/decade)', fontsize=14) cbar.ax.tick_params(labelsize=12) # 添加网格线 gl = ax.gridlines(crs=ccrs.PlateCarree(), draw_labels=True, linewidth=0.5, color='gray', alpha=0.5, linestyle='--') gl.top_labels = False gl.right_labels = False gl.xformatter = LONGITUDE_FORMATTER gl.yformatter = LATITUDE_FORMATTER gl.xlabel_style = {'size': 12} gl.ylabel_style = {'size': 12} # 6. 绘制区域边界框并添加趋势值 for region, bbox in regions.items(): # 绘制红色虚线边界框 rect = mpatches.Rectangle( xy=(bbox['lon_min'], bbox['lat_min']), width=bbox['lon_max'] - bbox['lon_min'], height=bbox['lat_max'] - bbox['lat_min'], linewidth=2.5, linestyle='--', edgecolor='red', facecolor='none', transform=ccrs.PlateCarree() ) ax.add_patch(rect) # 添加区域名称和趋势值 trend_value = regional_trends[region] center_lon = (bbox['lon_min'] + bbox['lon_max']) / 2 center_lat = (bbox['lat_min'] + bbox['lat_max']) / 2 # 添加趋势值文本 ax.text(center_lon, center_lat, f'{region}\n{trend_value:.3f}℃/decade', fontsize=14, weight='bold', color='black', horizontalalignment='center', verticalalignment='center', bbox=dict(facecolor='white', alpha=0.8, boxstyle='round,pad=0.5'), transform=ccrs.PlateCarree()) # 7. 添加标题和图例 plt.title('中国区域1980-2014年温度变化趋势空间分布', fontsize=18, pad=20) plt.figtext(0.5, 0.02, "数据来源: ChinaMet 0.25°气温数据集 | 红色虚线框表示四个子区域边界", ha="center", fontsize=12, color='gray') # 保存图像 output_path = "China_Regional_Temperature_Trends_1980-2014.png" plt.savefig(output_path, dpi=300, bbox_inches='tight') print(f"综合空间分布图已保存至: {output_path}") plt.show() ``` ### 代码说明： 1. **数据处理流程**： - 读取1980-2014年所有NC文件 - 计算每一年中国区域内的年平均温度 - 合并所有年份数据形成时间序列 - 计算每个网格点的温度变化趋势（℃/decade） - 计算每个子区域的平均升温趋势 2. **可视化设计**： - 使用Cartopy绘制专业地理空间图 - 创建蓝-白-红渐变色标清晰展示温度变化 - 添加海岸线、国界、河流等地理要素 - 用红色虚线框标记四个子区域边界 - 在每个区域中心位置显示区域名称和平均升温趋势值 3. **关键功能**： - 计算并显示各区域平均趋势值（保留3位小数） - 自动调整文本背景提高可读性 - 添加详细地图标注和数据来源说明 - 生成高分辨率图像（300dpi） 4. **输出结果**： - 生成综合空间分布图："China_Regional_Temperature_Trends_1980-2014.png" - 图中包含： - 中国区域温度变化趋势空间分布 - 四个子区域边界（红色虚线框） - 每个子区域的名称和平均升温趋势值 ### 使用说明： 1. **依赖库安装**： ```bash pip install xarray netCDF4 scipy matplotlib cartopy ``` 2. **区域边界调整**： - 在代码开头的`regions`字典中可修改各区域边界 - 格式：{"lat_min": ..., "lat_max": ..., "lon_min": ..., "lon_max": ...} 3. **可视化定制**： - 修改`cmap`变量可更改色标方案 - 调整`vmin`和`vmax`可改变色标范围 - 修改`fontsize`可调整文本大小 4. **结果解读**： - 蓝色区域：降温趋势 - 红色区域：升温趋势 - 框内数值：区域平均升温趋势（℃/decade）