实验数据是一个35,35的数组,采用线性插值得到375,375的数组。
构造interp2d
对象时,要求输入原始数据的行列坐标和数值,并指定插值方式。
然后向构造好的对象中输入想要得到的数组坐标值。
import matplotlib.pyplot as plt
import numpy as np
from scipy.interpolate import interp2d
residual = np.load(r'F:\Py_project\Spatial_RF\residual.npy').reshape(35, 35)
f = interp2d(np.arange(0, 375, 11), np.arange(0, 375, 11),
residual, kind='linear')
residual_hs = f(np.arange(375), np.arange(375))
结果及与ENVI双线性插值结果对比