scipy - Pythonで不良ピクセルをマークするための画像の補間

Question

pyfits.getdataこんにちは、画像の隅に大きなケラレがある画像を読み込んでいます。元の画像と同じ次元の 2 番目の配列を作成しnumpy、不良ピクセルの値を 1、使用可能なピクセルの値を 0 にpysexします。元の画像で使用する点光源を検出します。このpysexルーチンは、ケラレ領域の変動を点源として検出します。不良ピクセルをマークする配列を補間するにはどうすればよいですか。試してみscipy.interpolate.griddataましscipy.interpolate.interp2dたが、必要な引数が明確ではありません。最も簡単な解決策は何ですか?

Answer 1

map_coordinatesあなたが望むものをあなたに与えることができるはずだと思います：

index_vignetted = scipy.ndimage.map_coordinates(vignette, point_sources)

これにより、各点源の位置の下にあるピクセルのビネットの値を持つベクトルが返されます。

Answer 2

(x0,y0) = 5,10 に不良ピクセルがあり、(x1,y1) = 50,100 に別のピクセルがあります。

output = Scipy.interpolate.map_coordinat(array,np.array([[5,50],[10,100]]))
1bad_value = output[0]
2bad_value = output[1]

それは [x0,x1,x2,x3],[y0,y1,y2,y3] のようなもので、あまり直感的ではありません。

あなたはもっと欲しいです, これは本当に良い投稿です : グリッドデータの高速補間

しかし、実際には、悪いピクセルには scipy.filter.uniform_filter を使用します。ネイブールの平均値を示します。正確な aa キュービック補間ではありませんが、より堅牢です。2 つの不良ピクセルが互いに近くにある場合は、さらに優れた scipy.filter.median_filter を使用してください。