python – 重叠区域的卫星图像插值

我试图通过使用以下代码执行插值：

import scipy as sc
import scipy.interpolate as spi
import numpy as np

def interpolate_data(img,lat,lon):

p_x,p_y = sc.mgrid[-1:1:lat.shape[0]*1j,-1:1:lon.shape[1]*1j] 
n_x,n_y = sc.mgrid[-1:1:img.shape[0]*1j,-1:1:img.shape[1]*1j] 
n_lat = spi.griddata((p_x.ravel(),p_y.ravel()),lat.ravel(),(n_x,n_y),method='linear') 
n_lon = spi.griddata((p_x.ravel(),lon.ravel(),method='linear')
print "n_lat:",n_lat.shape
print "n_lon:",n_lon.shape
print "img:",img.shape
return n_lat,n_lon

但如果我以这种方式进行插值,我将得到不正确的数据.

在我的情况下,卫星将拍摄9次地球.在曲线A和曲线B之间是卫星在第一次拍摄中拍摄的区域.在曲线C和曲线D之间是卫星在第二次拍摄中所占的面积. E和F是重叠的区域.

如果我使用上面的代码进行插值,我会将重叠区域计算两次并给出错误的数据.所以我想问一下是否有办法只计算一次这些重叠区域.谢谢.

如果我使用上面的代码进行插值,我会将重叠区域计算两次并给出错误的数据.所以我想问一下是否有办法只计算一次这些重叠区域.谢谢.