如何对格点场进行升/降尺度？

zyfws

有没有函数直接对获得的格点场进行升降尺度，比如把0.25*0.25的格点场变成0.125*0.125？

MeteoInfo

可以用 interp2d 函数来实现（缺省是双线性插值），脚本示例程序：

1. f = addfile('D:/Temp/GrADS/model.ctl')
   
2. vname = 'PS'
   
3. psv = f[vname]
   
4. ps = psv[0,'10:60','60:140']
   
5. 
   
6. #Plot
   
7. subplot(1,2,1)
   
8. axesm(newaxes=False)
   
9. geoshow('country', edgecolor=(0,0,255))
   
10. layer = imshowm(ps, 20)
    
11. #scatterm(ps)
    
12. title('Pressure - low resolution')
    
13. colorbar(layer, orientation='horizontal', extendrect=False, shrink=0.8, aspect=12)
    
14. 
    
15. subplot(1,2,2)
    
16. axesm(newaxes=False)
    
17. geoshow('country', edgecolor=(0,0,255))
    
18. x = linspace(60., 140., 34)
    
19. y = linspace(10., 62., 28)
    
20. xx, yy = meshgrid(x, y)
    
21. ps1 = interp2d(ps, xx, yy)
    
22. layer = imshowm(x, y, ps1, 20)
    
23. #scatterm(x, y, ps1)
    
24. title('Pressure - high resolution')
    
25. colorbar(layer, orientation='horizontal', extendrect=False, shrink=0.8, aspect=12)

f = addfile('D:/Temp/GrADS/model.ctl')<br /> vname = 'PS'<br /> psv = f[vname]<br /> ps = psv[0,'10:60','60:140']<br /> <br /> #Plot<br /> subplot(1,2,1)<br /> axesm(newaxes=False)<br /> geoshow('country', edgecolor=(0,0,255))<br /> layer = imshowm(ps, 20)<br /> #scatterm(ps)<br /> title('Pressure - low resolution')<br /> colorbar(layer, orientation='horizontal', extendrect=False, shrink=0.8, aspect=12)<br /> <br /> subplot(1,2,2)<br /> axesm(newaxes=False)<br /> geoshow('country', edgecolor=(0,0,255))<br /> x = linspace(60., 140., 34)<br /> y = linspace(10., 62., 28)<br /> xx, yy = meshgrid(x, y)<br /> ps1 = interp2d(ps, xx, yy)<br /> layer = imshowm(x, y, ps1, 20)<br /> #scatterm(x, y, ps1)<br /> title('Pressure - high resolution')<br /> colorbar(layer, orientation='horizontal', extendrect=False, shrink=0.8, aspect=12)

御风者

ncl 里的esmf功能也可以实现。

zyfws

MeteoInfo 发表于 2018-1-23 23:43
可以用 interp2d 函数来实现（缺省是双线性插值），脚本示例程序：

谢谢王老师！