interpolate_grid.py 文件源码

def interpola(self,x,y,xi,yi,var,ndepths):
        Vari = np.zeros((xi.shape[0],ndepths))

        for k in range(ndepths):
            Vari[:,k] = interpolate.griddata((x,y),var[k,:,:].ravel(),(xi,yi),method='nearest')
        Vari[Vari==Vari.min()]=Vari.max()
        return Vari