MeteoInfoLab示例：三维等值线填色切片

MeteoInfo

contourfslice 函数可以绘制三维数组在 x, y, z 方向的多个等值线填色切片。

1. fn = os.path.join(migl.get_sample_folder(), 'GrADS', 'model.ctl')
   
2. f = addfile(fn)
   
3. data = f['U'][0,:,:,'0:180']
   
4. pres = data.dimvalue(0)
   
5. z = meteolib.pressure_to_height_std(pres)
   
6. y = data.dimvalue(1)
   
7. x = data.dimvalue(2)
   
8. 
   
9. #Plot
   
10. ax = axes3d()
    
11. geoshow('continent', facecolor=[255,231,177], edgecolor='b')
    
12. levs = arange(-16, 57, 8)
    
13. contourfslice(x, y, z, data, levs, edgecolor='lightgray', alpha=0.8,
    
14.     xslice=[120], yslice=[40], zslice=[z[3]])
    
15. colorbar()
    
16. xlim(0, 180)
    
17. title('contourfslice example')

fn = os.path.join(migl.get_sample_folder(), 'GrADS', 'model.ctl')<br /> f = addfile(fn)<br /> data = f['U'][0,:,:,'0:180']<br /> pres = data.dimvalue(0)<br /> z = meteolib.pressure_to_height_std(pres)<br /> y = data.dimvalue(1)<br /> x = data.dimvalue(2)<br /> <br /> #Plot<br /> ax = axes3d()<br /> geoshow('continent', facecolor=[255,231,177], edgecolor='b')<br /> levs = arange(-16, 57, 8)<br /> contourfslice(x, y, z, data, levs, edgecolor='lightgray', alpha=0.8, <br />     xslice=[120], yslice=[40], zslice=[z[3]])<br /> colorbar()<br /> xlim(0, 180)<br /> title('contourfslice example')

随缘

{:eb502:}

VV09

最近在学习剖切  真不错~~感谢分享

dlxieyr

感谢楼主，怎么运行提示NameError: name 'contourfslice' is not defined。如何解决，谢谢！

MeteoInfo

dlxieyr 发表于 2021-8-21 21:31
感谢楼主，怎么运行提示NameError: name 'contourfslice' is not defined。如何解决，谢谢！

使用MeteoInfo最新版本

朱府鸿

求助各位大神提供MeteoInfo调用MICAPS数据的具体实例PPT更好，若提供本人将不胜感激。

freekiller

王老师好！可以斜剖面吗？

MeteoInfo

freekiller 发表于 2021-9-3 15:10
王老师好！可以斜剖面吗？

可以，比较麻烦一些，有空了可以考虑做些简化。

1. fn = os.path.join(migl.get_sample_folder(), 'GrADS', 'model.ctl')
   
2. f = addfile(fn)
   
3. u = f['U'][0,:,'20','0']
   
4. lev1 = meteo.pressure_to_height_std(u.dimvalue(0)) / 1000
   
5. 
   
6. lon1 = 0.
   
7. lon2 = 180.
   
8. lat1 = -50.
   
9. lat2 = 50.
   
10. lon = lon1
    
11. lat_c = []
    
12. tdata = []
    
13. #tdata_terrain = []
    
14. while lon <= lon2:
    
15.     lat = lat1 + (lat2 - lat1) * (lon - lon1) / (lon2 - lon1)
    
16.     lat_c.append(lat)
    
17.     print lon
    
18.     tdata.append(f['U'][0,:,'%f'%lat,'%f'%lon])
    
19.     lon = lon + 2.5
    
20.    
    
21. alldata = concatenate(tdata, axis=0)
    
22. alldata = reshape(alldata, (len(tdata), len(lev1)))   
    
23. alldata = alldata.T
    
24. x = arange(lon1, lon2 + 1, 2.5)
    
25. 
    
26. #Plot
    
27. ax = axes3d(tickfontsize=14)
    
28. geoshow('continent', color='c', edgecolor='b')
    
29. ls = ax.contourf(x, lev1, alldata, 10, edgecolor='k', zdir='xy', alpha=0.8, \
    
30.     sepoint=[lon1,lat1,lon2,lat2])
    
31. colorbar(ls)
    
32. xlabel('Longitude (degrees)', fontsize=16)
    
33. ylabel('Latitude (degrees)', fontsize=16)
    
34. zlabel('Altitude (km)', fontsize=16)
    
35. xlim(0, 180)
    
36. title('vertical slice plot example')

fn = os.path.join(migl.get_sample_folder(), 'GrADS', 'model.ctl')<br /> f = addfile(fn)<br /> u = f['U'][0,:,'20','0']<br /> lev1 = meteo.pressure_to_height_std(u.dimvalue(0)) / 1000<br /> <br /> lon1 = 0.<br /> lon2 = 180.<br /> lat1 = -50.<br /> lat2 = 50.<br /> lon = lon1<br /> lat_c = []<br /> tdata = []<br /> #tdata_terrain = []<br /> while lon <= lon2:<br />     lat = lat1 + (lat2 - lat1) * (lon - lon1) / (lon2 - lon1)<br />     lat_c.append(lat)<br />     print lon<br />     tdata.append(f['U'][0,:,'%f'%lat,'%f'%lon])<br />     lon = lon + 2.5<br />     <br /> alldata = concatenate(tdata, axis=0)<br /> alldata = reshape(alldata, (len(tdata), len(lev1)))    <br /> alldata = alldata.T<br /> x = arange(lon1, lon2 + 1, 2.5)<br /> <br /> #Plot<br /> ax = axes3d(tickfontsize=14)<br /> geoshow('continent', color='c', edgecolor='b')<br /> ls = ax.contourf(x, lev1, alldata, 10, edgecolor='k', zdir='xy', alpha=0.8, \<br />     sepoint=[lon1,lat1,lon2,lat2])<br /> colorbar(ls)<br /> xlabel('Longitude (degrees)', fontsize=16)<br /> ylabel('Latitude (degrees)', fontsize=16)<br /> zlabel('Altitude (km)', fontsize=16)<br /> xlim(0, 180)<br /> title('vertical slice plot example')

freekiller

MeteoInfo 发表于 2021-9-3 16:57
可以，比较麻烦一些，有空了可以考虑做些简化。

谢谢王老师！这下就灵活多了

oydl1985

请问下王老师这个三维剖面可以叠加地形数据吗？