deconvlucy
ルーシー・リチャードソン アルゴリズムを使用したイメージのブレ除去
構文
説明
例
ブレやノイズのない初期状態のイメージを読み取って表示します。この例ではオプションで、左上 (x,y) 座標が (2,50) にある 256 x 256 のサイズにイメージをトリミングします。
I = imread("board.tif"); I = imcrop(I,[2 50 255 255]); imshow(I) title("Original Image")
標準偏差 5 とサイズ 5 x 5 のフィルターを使用して、ガウスぼかしを表す PSF を作成します。
PSF = fspecial("gaussian",5,5);イメージのブレをシミュレートします。
blurred = imfilter(I,PSF,"symmetric","conv");
シミュレートしたゼロ平均のガウス ノイズを加えます。
V = 0.002; blurred_noisy = imnoise(blurred,"gaussian",0,V); imshow(blurred_noisy) title("Blurred and Noisy Image")
deconvlucy を使用して、ブレとノイズを含むイメージを復元します。ブレを作成するときに使用した PSF を指定し、反復回数を 5 に減らします。
luc1 = deconvlucy(blurred_noisy,PSF,5);
imshow(luc1)
title("Restored Image")
サンプル イメージを作成し、ブレを加えます。
I = checkerboard(8); PSF = fspecial('gaussian',7,10); V = .0001; BlurredNoisy = imnoise(imfilter(I,PSF),'gaussian',0,V);
重み配列を作成し、各オプションのパラメーターを使用して deconvlucy を呼び出します。
WT = zeros(size(I)); WT(5:end-4,5:end-4) = 1; J1 = deconvlucy(BlurredNoisy,PSF); J2 = deconvlucy(BlurredNoisy,PSF,20,sqrt(V)); J3 = deconvlucy(BlurredNoisy,PSF,20,sqrt(V),WT);
結果を表示します。
subplot(221);imshow(BlurredNoisy); title('A = Blurred and Noisy'); subplot(222);imshow(J1); title('deconvlucy(A,PSF)'); subplot(223);imshow(J2); title('deconvlucy(A,PSF,NI,DP)'); subplot(224);imshow(J3); title('deconvlucy(A,PSF,NI,DP,WT)');
入力引数
出力引数
ヒント
deconvlucyを使用すると、逆畳み込みを前回停止した箇所から再開できます。この機能を使用するには、入力イメージIを cell 配列{I}として渡します。これを行う場合、関数deconvlucyは、出力イメージJを cell 配列として返します。これを次に呼び出すdeconvlucyに入力配列として渡すことができます。出力 cell 配列Jには次の 4 つの要素が含まれています。J{1}は元のイメージIを格納します。J{2}は最後の反復の結果を格納します。J{3}は最後から 2 番目の反復の結果を格納します。J{4}は反復アルゴリズムによって生成される配列です。出力イメージ
Jは、アルゴリズムで使用されている離散フーリエ変換に起因するリンギングを示すことがあります。リンギングを軽減するには、deconvlucyを呼び出す前にI = edgetaper(I,psf)を使用してください。deconvlucyは、PSF を正規化しないでdoubleに変換します。deconvlucyは、出力イメージが入力イメージの範囲外の領域に入る場合にも値を返すことができます。
参照
[1] D. S. C. Biggs and M. Andrews, Acceleration of iterative image restoration algorithms, Applied Optics, Vol. 36, No. 8, 1997.
[2] R. J. Hanisch, R. L. White, and R. L. Gilliland, Deconvolutions of Hubble Space Telescope Images and Spectra, Deconvolution of Images and Spectra, Ed. P.A. Jansson, 2nd ed., Academic Press, CA, 1997.
