python小白在线求助

不会发光i

使用python绘制平均场，距平图和纬偏图代码如下：
不知道错误在哪，求帮助
import numpy as np
import matplotlib.pyplot as plt
f=open('CL01\hgt500.grd','rb')
a=np.fromfile(f,dtype=np.float32).reshape(144,73,12,65)
f.close()
ave = np.mean(a[:,:,0,0:60],axis=2)
anom = a[:,:,0,60]-ave
avelon = np.mean(a[:,:,1,60],axis=0)
ho = a[:,:,1,60] - avelon
lon = np.arange(0,360,2.5)
lat = np.arange(-90,92.5,2.5)
lats,lons= np.meshgrid(lat,lon)

fig, ax = plt.subplots()
CS = ax.contour(lons,lats,ho)

雨落森林

python重组数据是从右往左重组的，也就是说，你重组成一个(a,b,c,d)的形状，它是先满足每d个数据作为一小组，然后每c个小组数据又是一个大组，然后每b个大组又是一个大大组，然后一共有a个大大组。而grd数据是按经度、纬度、高度、变量、时间存放的，那么你在python，应该是(时间、变量、高度、纬度、经度)这个顺序，其中，时间你可以(年份，月份)这样子。所以应该是65*12*73*144

一山止水

我觉得你的维度放错了，应该是先时间然后再空间。。a(65,12,144,73)

Masterpiece

可能维度有些问题，最好再确认一下

好久不见

维度问题吧。感觉应该是65*12*73*144。144*73*12*65应该是MATLAB的习惯吧。。。

不会发光i

问题已经解决，确实是几位指出的错误。十分感谢楼上几位。