オーバーラップ領域を含む衛星画像の補間

Question

衛星画像の補間に問題があります。私は理解してどのように地球の衛星撮影画像：

と私は、次のコードを使用して補間を行うためにしようとしています：

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は重なる領域です。

補間に上記のコードを使用すると、重複領域が2回カウントされ、誤ったデータが返されます。だから私は、それらの重複領域を数えて1回だけカウントする方法があるかどうか尋ねたい。ありがとう。

私は補間する上記のコードを使用している場合、私は二回重複領域をカウントし、私に間違ったデータが得られます。だから私は、それらの重複領域を数えて1回だけカウントする方法があるかどうか尋ねたい。ありがとう。

Answer 1

イメージに地理的エクステントデータ（GeoTIFFファイルなど）が添付されている場合は、GDAL libraryの使用を検討することができますが、このライブラリはC/C++で書かれていますが、python packageもあります。

これは、数多くの異なるフォーマットで動作し、画像のワーピングおよび変換のためのいくつかの異なる画像処理方法を備えています。

Answer 2

9ショットすべてのデータを連結して一度に補間することはできませんか？私はこれらの9つのデータセットを別々に扱う理由は見当たりません。あなたが持っている場合、私はあなたがより良い9つの異なるイメージを作成すると思います。