我遇到了衛星圖像插值問題。我瞭解地球的衛星拍攝的圖像:重疊區域的衛星圖像插值
,我嘗試使用下面的代碼來進行內插:
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(), p_y.ravel()), lon.ravel(), (n_x, n_y), 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是重疊的區域。
如果我使用上面的代碼進行插值,我將計算兩次重疊區域並給我錯誤的數據。所以我想問一下,是否有一種方法來統計這些重疊區域。謝謝。
如果我用上面的代碼進行插值,我會算兩次重疊區域,並給我錯誤的數據。所以我想問一下,是否有一種方法來統計這些重疊區域。謝謝。