如何用Numpy reshape()将二维气象数组转为匹配经纬度的三维格式

现有一组按经纬度顺序排列的气象数据数组，形状为(1038240,4)——每一行对应一个经纬度点位的4项气象数据，需要将其重塑为(4,721,1440)的三维数组，以此对应4个气象变量在721×1440的全球网格上的分布。

我试过用下面的代码：

newarr = t_new.reshape(4,721,1440)

虽然能得到目标形状，但没法让前两列的经纬度坐标匹配到预期的网格格式。

举个简化的例子：

• 原始示例数据是形状为(6,4)的数组，前两列是经纬度值，后两列是两个气象变量
• 预期转换后得到形状为(2,3,2)的三维数组，对应2个气象变量在3×2的网格上的分布

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