Python读取nc数据和绘图

亲三儿

平流层的萝卜 发表于 2020-7-17 10:16
u和v，分别平方，相加，再根号

好的，谢谢

wxylndd

平流层的萝卜 发表于 2020-7-1 09:26
[0,35:76,70:141] 0是指第一个时次，35:76是指在南北方向上，从北到南第34-74个格点，70:141是指东西方向 ...

不应该是第35~75个格点吗。。不太懂

平流层的萝卜

wxylndd 发表于 2020-7-24 11:10
不应该是第35~75个格点吗。。不太懂

在python里，索引是从0开始的，0,1,2,3......

wxylndd

平流层的萝卜 发表于 2020-7-24 11:27
在python里，索引是从0开始的，0,1,2,3......

好的，谢谢

wxylndd

楼主你好，我想请问等温线不显示怎么解决。。代码如下

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,42:71,13:66]
    
18. plt.figure(figsize=(16,8))
    
19. map=Basemap(projection='cyl',llcrnrlat=5,urcrnrlat=55,llcrnrlon=70,urcrnrlon=140,resolution='l')
    
20. map.drawcoastlines()
    
21. map.drawcountries(linewidth=1.5)
    
22. map.drawcounties
    
23. map.drawparallels(np.arange(3.,55.,5.),labels=[1,0,0,0],fontsize=15)
    
24. map.drawmeridians(np.arange(70.,140.,10.),labels=[0,0,0,1],fontsize=15)
    
25. map.drawlsmask()
    
26. Lon,Lat=np.meshgrid(lon[13:66],lat[42:71])
    
27. x,y=map(Lon,Lat)
    
28. c=plt.contour(x,y,temperature_need,colors='red')
    
29. plt.title('20090418 00h 500hpa',size=20)
    
30. plt.savefig('air1.png')
    
31. plt.show()
    

复制代码

出来的图没有等温线。。。

平流层的萝卜

temperature_need = temperature[429,level_need,42:71,13:66]，你要保证这几个值在图上是可以画出来的。意思就是值在可以画的区间内。

shi_8517

牛人，好好学习学习

少年易

感谢楼主分享

来自苹果APP客户端

天边的那边

真的是太有用

不戒香烟

前来学习，膜拜大佬