求问python画温度场图的时候出现错误怎么解决

wxylndd

import netCDF4
from netCDF4 import Dataset
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
nc_obj=Dataset('D://Anaconda//envs//py27//pressure-Part2//pressure//air.2009.nc')
print(nc_obj)
lon=nc_obj.variables['lon'][:]
lon = np.array(lon)
lat=nc_obj.variables['lat'][:]
lat = np.array(lat)
temperature = nc_obj.variables['air'][:]
输入到这一步时显示：MemoryError: Unable to allocate 995. MiB for an array with shape (1460, 17, 73, 144) and data type float32
我以为是空间不够，修改了虚拟内存大小，还是不行
求助各位大神

帅帅f

不应该啊 代码没问题啊  你把你数据给我看看

wxylndd

帅帅f 发表于 2020-7-10 17:12
不应该啊 代码没问题啊  你把你数据给我看看

好像是我的python版本问题，我之前是32位，用了64位的试了下，没有这个提示了，谢谢你

kingjames

完整代码方便分享学习下吗？

wxylndd

kingjames 发表于 2020-7-11 21:32
完整代码方便分享学习下吗？

1. import netCDF4 as nc
   
2. from netCDF4 import Dataset
   
3. import numpy as np
   
4. import matplotlib.pyplot as plt
   
5. from mpl_toolkits.basemap import Basemap
   
6. nc_obj=Dataset('D://Anaconda//pressure-Part2//pressure//air.2009.nc')
   
7. print(nc_obj.variables.keys())#打印变量信息
   
8. lon=nc_obj.variables['lon']
   
9. lon=np.array(lon)
   
10. lat=nc_obj.variables['lat']
    
11. lat=np.array(lat)
    
12. temperature=nc_obj.variables['air']
    
13. temperature=np.array(temperature)
    
14. level=nc_obj.variables['level']
    
15. level = np.array(level)
    
16. level_need = 5
    
17. temperature_need = temperature[429,level_need,5:34,13:66]
    
18. plt.figure(figsize=(16,8))
    
19. map=Basemap(projection='cyl',llcrnrlat=10,urcrnrlat=90,llcrnrlon=25,urcrnrlon=165,resolution='l')
    
20. map.drawcoastlines()
    
21. map.drawcountries(linewidth=1.5)
    
22. map.drawparallels(np.arange(10.,90.,5.),labels=[1,0,0,0],fontsize=15)
    
23. map.drawmeridians(np.arange(25.,165.,15.),labels=[0,0,0,1],fontsize=15)
    
24. map.drawlsmask()
    
25. Lon,Lat=np.meshgrid(lon[13:66],lat[5:34])
    
26. x,y=map(Lon,Lat)
    
27. c=plt.contour(x,y,temperature_need,colors='red')
    
28. plt.title('20090418 00h 500hpa',size=20)
    
29. plt.savefig('air1.png')
    
30. plt.show()
    

复制代码