高さの等高線を含む地形図があり、scipy.interpolate.griddataを使用して、そこから高さマップ(各セルの高さを値として持つグリッド)を生成したいと考えています。私の考えは、等高線のポイントを非構造化データポイントとして内挿アルゴリズムに追加し、グリッド全体の内挿値(高さ)を尋ねることです。
何らかの理由で、griddata()は、ほとんどが2つの等高線間で補間されないグリッドを返し、2つの等高線間のグリッドポイントは、2つの等高線値間の値ではなく、低い値を取得します。
以下は、等高線(左)と結果の補間グリッド(1000x1000)を含む私の例の20x20画像です。等高線からのゼロ以外の値(黒以外のピクセル)のみが、補間への入力として使用されます。同じ問題は、「最も近い」、「線形」、および「立方体」の内挿法でも発生します。期待しているハイマップを取得できない理由はありますか?
ソースコードはここにあります-http://pastebin.com/umNmSrE2
更新:問題は解決しました。実際に悪い入力画像が問題でした(以下のコメントを参照)。新しい結果は次のとおりです。