python - 整列されたポイントにぶら下がらないscipy.interpolate.griddataの代替

Question

グリッド上で内挿しようとしているポイントデータセットがあります。これらのポイントはグリッド形式で配置され、一部のポイントが欠落しています。以下を参照してください。

複雑なことに、他の入力ポイントセットがグリッド上で整列していない可能性があるため、scipy.interpolate.griddataこれらの値を通常のグリッドに補間するために使用しようとしています。ただし、基になるグリッドが入力ポイントデータセットのサンプリングレートと完全に一致し、griddataハングすることがあります。この質問 によるとscipy.interpolate.griddata、サンプルが完全に整列している3つのポイントにまたがって発生する場合、これは時々発生しますが、うまく機能しません。

Pythonの通常のグリッドにポイントデータを補間するための高性能な代替手段はありますか？

Answer 1

私の解決策は、各入力点のx座標とy座標を、の間のランダムな係数で確率的に摂動させること1e-6 * grid_cell_sizeでした。1e-6経験則として、32ビット浮動小数点の最小デルタ差よりも数桁大きいが、それによって発生するエラーが補間スキームによって発生するエラーによって隠されるほど小さいものを選択します。

私はまだ他のアイデアを受け入れていますが、このようにしscipyて、非常に小さな数値誤差を導入することによって提供される株式補間関数を使用できます。