MeteoInfoLab脚本示例：自动标注气压高低中心值

MeteoInfo

可以用minimum_filter和maximum_filter函数来定位气压低值和高值中心。

1. f = addfile('D:/Temp/GrADS/model.ctl')
   
2. pres = f['PS'][0,'-60:70','50:280']
   
3. lat = pres.dimvalue(0)
   
4. lon = pres.dimvalue(1)
   
5. pres = imagelib.gaussian_filter(pres, sigma=3.0)
   
6. mxdata = imagelib.minimum_filter(pres, size=25)
   
7. mxy, mxx = np.where(mxdata == pres)
   
8. midata = imagelib.maximum_filter(pres, size=25)
   
9. miy, mix = np.where(midata == pres)
   
10. 
    
11. #Plot
    
12. axesm()
    
13. geoshow('country')
    
14. layer = contourm(lon, lat, pres, 10, color=(0,153,204))
    
15. clabel(layer, fontsize=10)
    
16. 
    
17. for j in range(len(mxy)):
    
18.     geoshow(lat[mxy[j]], lon[mxx[j]], color='b', edgecolor=None, size=4)
    
19.     text(lon[mxx[j]], lat[mxy[j]], 'L', color='r', xalign='center',
    
20.         yalign='center', fontsize=16)
    
21.     print(mxdata[mxy[j],mxx[j]])
    
22. 
    
23. for j in range(len(miy)):
    
24.     geoshow(lat[miy[j]], lon[mix[j]], color='r', edgecolor=None, size=4)
    
25.     text(lon[mix[j]], lat[miy[j]], 'H', color='b', xalign='center',
    
26.         yalign='center', fontsize=16)
    
27.     print(midata[miy[j],mix[j]])
    
28. 
    
29. title('Plot low and high pressure centers')

f = addfile('D:/Temp/GrADS/model.ctl')<br /> pres = f['PS'][0,'-60:70','50:280']<br /> lat = pres.dimvalue(0)<br /> lon = pres.dimvalue(1)<br /> pres = imagelib.gaussian_filter(pres, sigma=3.0)<br /> mxdata = imagelib.minimum_filter(pres, size=25)<br /> mxy, mxx = np.where(mxdata == pres)<br /> midata = imagelib.maximum_filter(pres, size=25)<br /> miy, mix = np.where(midata == pres)<br /> <br /> #Plot<br /> axesm()<br /> geoshow('country')<br /> layer = contourm(lon, lat, pres, 10, color=(0,153,204))<br /> clabel(layer, fontsize=10)<br /> <br /> for j in range(len(mxy)):<br />     geoshow(lat[mxy[j]], lon[mxx[j]], color='b', edgecolor=None, size=4)<br />     text(lon[mxx[j]], lat[mxy[j]], 'L', color='r', xalign='center',<br />         yalign='center', fontsize=16)<br />     print(mxdata[mxy[j],mxx[j]])<br /> <br /> for j in range(len(miy)):<br />     geoshow(lat[miy[j]], lon[mix[j]], color='r', edgecolor=None, size=4)<br />     text(lon[mix[j]], lat[miy[j]], 'H', color='b', xalign='center',<br />         yalign='center', fontsize=16)<br />     print(midata[miy[j],mix[j]])<br /> <br /> title('Plot low and high pressure centers')

kmdqlb

王老师威武

river

这个厉害了，很实用的功能

kmdqlb

谢谢王老师

随缘

{:eb502:}

karenlk

老师您好，最近meteoinfo下载不起来啊

MeteoInfo

karenlk 发表于 2020-2-27 15:34
老师您好，最近meteoinfo下载不起来啊

可以在百度网盘上下载：https://pan.baidu.com/s/1jHFAxoy

karenlk

MeteoInfo 发表于 2020-2-27 16:24
可以在百度网盘上下载：https://pan.baidu.com/s/1jHFAxoy

谢谢老师，太棒了