MeteoInfoLab脚本示例：数据插值

MeteoInfo

MeteoInfoLab的numeric包中有interpolate模块，其中有两个类（class）用于一维和二维插值：interp1d, RectBivariateSpline。利用已知数据创建类的对象后可以像函数一样使用进行数据插值计算。

一维插值：

1. from mipylib.numeric import interpolate
   
2. 
   
3. x = linspace(0, 10, num=11, endpoint=True)
   
4. y = cos(-x**2/9.0)
   
5. f = interpolate.interp1d(x, y)
   
6. f2 = interpolate.interp1d(x, y, kind='cubic')
   
7. 
   
8. xnew = linspace(0, 10, num=100, endpoint=True)
   
9. plot(x, y, 'bo', xnew, f(xnew), 'g-', xnew, f2(xnew), 'r--')
   
10. ylim(-1.5, 1.2)
    
11. legend(['data','linear','cubic'], loc='lower left')
    
12. title('Interpolation example')

from mipylib.numeric import interpolate<br /> <br /> x = linspace(0, 10, num=11, endpoint=True)<br /> y = cos(-x**2/9.0)<br /> f = interpolate.interp1d(x, y)<br /> f2 = interpolate.interp1d(x, y, kind='cubic')<br /> <br /> xnew = linspace(0, 10, num=100, endpoint=True)<br /> plot(x, y, 'bo', xnew, f(xnew), 'g-', xnew, f2(xnew), 'r--')<br /> ylim(-1.5, 1.2)<br /> legend(['data','linear','cubic'], loc='lower left')<br /> title('Interpolation example')

二维插值（矩形格点数据）：

1. from mipylib.numeric import interpolate
   
2. 
   
3. x, y = meshgrid(linspace(-1, 1, 20), linspace(-1, 1, 20))
   
4. z = (x+y) * exp(-6.0*(x*x+y*y))
   
5. f = interpolate.RectBivariateSpline(x[0,:], y[:,0], z)
   
6. 
   
7. subplot(1,2,1)
   
8. im = imshow(x[0,:], y[:,0], z, 40)
   
9. colorbar(im)
   
10. 
    
11. xnew, ynew = meshgrid(linspace(-1, 1, 70), linspace(-1, 1, 70))
    
12. znew = f(xnew, ynew)
    
13. subplot(1,2,2)
    
14. im = imshow(xnew[0,:], ynew[:,0], znew, 40)
    
15. colorbar(im)
    
16. suptitle('2-D grid data interpolation example')

from mipylib.numeric import interpolate<br /> <br /> x, y = meshgrid(linspace(-1, 1, 20), linspace(-1, 1, 20))<br /> z = (x+y) * exp(-6.0*(x*x+y*y))<br /> f = interpolate.RectBivariateSpline(x[0,:], y[:,0], z)<br /> <br /> subplot(1,2,1)<br /> im = imshow(x[0,:], y[:,0], z, 40)<br /> colorbar(im)<br /> <br /> xnew, ynew = meshgrid(linspace(-1, 1, 70), linspace(-1, 1, 70))<br /> znew = f(xnew, ynew)<br /> subplot(1,2,2)<br /> im = imshow(xnew[0,:], ynew[:,0], znew, 40)<br /> colorbar(im)<br /> suptitle('2-D grid data interpolation example')

kmdqlb

王老师威武

ddd212

学习了。。。

liuxiunuli

学习了。。。谢谢楼主{:5_213:}