Python绘制水汽通量图问题求助

Tao_Xu

1. from netCDF4 import Dataset
   
2. import numpy as np
   
3. import matplotlib.pyplot as plt
   
4. from mpl_toolkits.basemap import Basemap
   
5. import matplotlib.cm as cm
   
6. 
   
7. meteo_file = "/Users/ACER-PC/PycharmProjects/ERA-interim_Vertical integral of divergence of moisture flux/Data/ERA-interim_Vertical integral of divergence of moisture flux201812.nc"
   
8. fh = Dataset(meteo_file, mode='r')
   
9. 
   
10. # 获取每个变量的值
    
11. lons = fh.variables['longitude'][:]
    
12. lats = fh.variables['latitude'][:]
    
13. tlml = fh.variables['p84.162'][:]
    
14. elml = fh.variables['p71.162'][:]
    
15. nlml = fh.variables['p72.162'][:]
    
16. 
    
17. 
    
18. 
    
19. m = Basemap(llcrnrlon=50., llcrnrlat=-10.,
    
20.              urcrnrlon=160., urcrnrlat=60,
    
21.              projection='cyl', resolution='l')
    
22. lon, lat = np.meshgrid(lons, lats)
    
23. xi, yi = m(lon, lat)
    
24. 
    
25. 
    
26. # 数据只绘制第1个时间的（tlml_0）
    
27. tlml_0 = tlml[0:1:, ::, ::]
    
28. cs = m.contourf(xi, yi, np.squeeze(tlml_0), range (-21,1), cmap=cm.bone)
    
29. elml_0 = elml[0:1:, ::, ::]
    
30. nlml_0 = nlml[0:1:, ::, ::]
    
31. wind= m.quiver(lon, lat, np.squeeze(elml_0), np.squeeze(nlml_0), color='red', width=0.0005, headwidth=60, headlength=100)
    
32. 
    
33. 
    
34. # Add Grid Lines
    
35. # 绘制经纬线
    
36. m.drawparallels(np.arange(-90., 91., 20.),labels=[True,False,True,True],color='dimgrey',dashes=[1, 3])
    
37. m.drawmeridians(np.arange(-180., 181., 30.),labels=[True,True,False,True],color='dimgrey',dashes=[1, 3])
    
38. 
    
39. # Add Coastlines
    
40. m.drawcoastlines(color='grey', linewidth=0.3)
    
41. 
    
42. 
    
43. 
    
44. # Add Colorbar
    
45. cbar = m.colorbar(cs, location='bottom', pad="10%")
    
46. cbar.set_label('Divergence of moisture flux (kg m**-2 s**-1)')
    
47. 
    
48. # Add Title
    
49. plt.title('Vertical integral of divergence of moisture flux', fontsize= 12)
    
50. plt.savefig("Vertical integral of divergence of moisture flux.png", dpi=500, bbox_inches='tight')
    
51. plt.show()
    

复制代码

Tao_Xu

好像帖子没显示，问题是：最近在利用ERA-interim中计算好的整层水汽通量散度（Vertical integral of divergence of moisture flux）和水汽通量（Vertical integral of eastward water vapour flux及Vertical integral of northward water vapour flux）数据来画图，图中的箭头总是很小（我想用箭头表示水汽通量值的大小），需要图放很大才能看见，不知道是哪里出了问题。刚开始学习Python，不是很熟练，希望好心人帮忙看一下脚本，谢谢，不胜感激！！

Tao_Xu

Tao_Xu 发表于 2019-9-1 21:14
好像帖子没显示，问题是：最近在利用ERA-interim中计算好的整层水汽通量散度（Vertical integral of diverg ...

from netCDF4 import Dataset
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
import matplotlib.cm as cm

meteo_file = "/Users/ACER-PC/PycharmProjects/ERA-interim_Vertical integral of divergence of moisture flux/Data/ERA-interim_Vertical integral of divergence of moisture flux201812.nc"
fh = Dataset(meteo_file, mode='r')

# 获取每个变量的值
lons = fh.variables['longitude'][:]
lats = fh.variables['latitude'][:]
tlml = fh.variables['p84.162'][:]
elml = fh.variables['p71.162'][:]
nlml = fh.variables['p72.162'][:]



m = Basemap(llcrnrlon=50., llcrnrlat=-10.,
             urcrnrlon=160., urcrnrlat=60,
             projection='cyl', resolution='l')
lon, lat = np.meshgrid(lons, lats)
xi, yi = m(lon, lat)


# 数据只绘制第1个时间的（tlml_0）
tlml_0 = tlml[0:1:, ::, ::]
cs = m.contourf(xi, yi, np.squeeze(tlml_0), range (-21,1), cmap=cm.bone)
elml_0 = elml[0:1:, ::, ::]
nlml_0 = nlml[0:1:, ::, ::]
wind= m.quiver(lon, lat, np.squeeze(elml_0), np.squeeze(nlml_0), color='red', width=0.0005, headwidth=60, headlength=100)


# Add Grid Lines
# 绘制经纬线
m.drawparallels(np.arange(-90., 91., 20.),labels=[True,False,True,True],color='dimgrey',dashes=[1, 3])
m.drawmeridians(np.arange(-180., 181., 30.),labels=[True,True,False,True],color='dimgrey',dashes=[1, 3])

# Add Coastlines
m.drawcoastlines(color='grey', linewidth=0.3)



# Add Colorbar
cbar = m.colorbar(cs, location='bottom', pad="10%")
cbar.set_label('Divergence of moisture flux (kg m**-2 s**-1)')

# Add Title
plt.title('Vertical integral of divergence of moisture flux', fontsize= 12)
plt.savefig("Vertical integral of divergence of moisture flux.png", dpi=500, bbox_inches='tight')
plt.show()

Luu_

完全没有在图里看到箭头耶

Luu_

试一下用颜色来表示大小？

yeti

厉害，学习了

小其其格

虽然不是很懂python，但是看你脚本里绘制风场的语句：wind= m.quiver(lon, lat, np.squeeze(elml_0), np.squeeze(nlml_0), color='red', width=0.0005, headwidth=60, headlength=100)

这里控制箭头大小的三个参数是 width, headwidth, headlength，你把数值都调整大一些试试，都放大10倍看看呗

对方的凤飞飞

涛哥，建议你的颜色棒换成彩虹颜色，这个最小的白色和你底图都重合了

Tao_Xu

对方的凤飞飞 发表于 2019-9-6 14:48
涛哥，建议你的颜色棒换成彩虹颜色，这个最小的白色和你底图都重合了

其实是故意的，为了让大于0的数据不显示出来

Tao_Xu

Tao_Xu 发表于 2019-9-1 21:14
好像帖子没显示，问题是：最近在利用ERA-interim中计算好的整层水汽通量散度（Vertical integral of diverg ...

找到原因了，将数据的分辨率从0.75*0.75变到2.5*2.5后，箭头就变大了，图上数据点也变的没有那么密集了