如何对xarray数据的每个经纬度点计算小时温度最大值并绘图

解决方案 #

你之前直接调用mx2t.max()未指定聚合维度，xarray会默认对数组的所有维度（time、latitude、longitude）做聚合，因此只会返回单个全局最大值。你只需要在调用max方法时指定仅对时间维度做聚合，就能保留经纬度维度得到每个格点的最大值：

# 仅在time维度上求最大值，输出结果保留latitude、longitude两个维度
point_max = mx2t.max(dim="time")

得到的point_max就是每个经纬度格点对应的所有时刻温度最大值，可直接用于后续地图绘制。
如果后续有其他聚合需求，也可以传入维度列表同时对多个维度聚合，示例如下：

# 同时对时间、高度层两个维度求最大值
mx2t.max(dim=["time", "level"])

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