MeteoInfoLab脚本示例：格点插值到格点

MeteoInfo

格点数据有时需要插值到另一个范围和分辨率的网格上，可以用 linint2 函数来插值。

1. fn = os.path.join(migl.get_sample_folder(), 'GrADS', 'model.ctl')
   
2. f = addfile(fn)
   
3. ps = f['PS'][:]
   
4. ps = ps[:,'10:60','60:140']
   
5. lon = arange(50, 141, 2.5)
   
6. lat = arange(5, 66, 2.5)
   
7. 
   
8. #Interpolate
   
9. nps = interpolate.linint2(ps, lon, lat)
   
10. 
    
11. #Plot
    
12. levs = arange(500, 1021, 20)
    
13. subplot(2,1,1,axestype='map')
    
14. geoshow('continent', edgecolor=(0,0,255))
    
15. imshow(ps[1,:,:], levs)
    
16. title('Pressure - origin')
    
17. colorbar()
    
18. 
    
19. subplot(2,1,2,axestype='map')
    
20. geoshow('continent', edgecolor=(0,0,255))
    
21. imshow(lon, lat, nps[1,:,:], levs)
    
22. title('Pressure - linint2')
    
23. colorbar()

fn = os.path.join(migl.get_sample_folder(), 'GrADS', 'model.ctl')<br /> f = addfile(fn)<br /> ps = f['PS'][:]<br /> ps = ps[:,'10:60','60:140']<br /> lon = arange(50, 141, 2.5)<br /> lat = arange(5, 66, 2.5)<br /> <br /> #Interpolate<br /> nps = interpolate.linint2(ps, lon, lat)<br /> <br /> #Plot<br /> levs = arange(500, 1021, 20)<br /> subplot(2,1,1,axestype='map')<br /> geoshow('continent', edgecolor=(0,0,255))<br /> imshow(ps[1,:,:], levs)<br /> title('Pressure - origin')<br /> colorbar()<br /> <br /> subplot(2,1,2,axestype='map')<br /> geoshow('continent', edgecolor=(0,0,255))<br /> imshow(lon, lat, nps[1,:,:], levs)<br /> title('Pressure - linint2')<br /> colorbar()

随缘

{:eb502:}{:eb502:}