interpolate_grid.py 文件源码

def interpolate_in_depth(self,z,zi,Var):
        shp = Var.shape
        ndepths = z.size

        shp1= [shp[0],ndepths]

        im, zm = np.meshgrid(np.arange(shp[0]),z)
        imi, zmi = np.meshgrid(np.arange(shp[0]),zi)

        im = im.T.ravel()
        zm= zm.T.ravel()

        imi = imi.T.ravel()
        zmi= zmi.T.ravel()

        Vari = interpolate.griddata((im,zm),Var.ravel(),(imi,zmi)).reshape(shp[0],zi.size)

        return Vari