画像処理のプログラムについて

Question

function I_rec = notSoBaselineAnymoreFast(I, mask)
%set not set pixel to zero, just to be sure
I(~mask) = 0;

% calculate inpainting with a 3x3 kernel
souroundedPatchCount3 = imfilter(double(mask),ones(3));
I_rec = imfilter(I,ones(3)) ./ souroundedPatchCount3;
% calculate inpainting with a 5x5 kernel
souroundedPatchCount5 = imfilter(double(mask),ones(5));
I_rec5 = imfilter(I,ones(5)) ./ souroundedPatchCount5;
% remove dev by zero by overwriting it with a overall mean, just in case.
I_rec5(souroundedPatchCount5 <= 0) = mean( I(mask == 1) );

% set in the missing pixel in the smaller aproximation
I_rec(souroundedPatchCount3 <= 0) = I_rec5(souroundedPatchCount3 <= 0);

%%% post processing
I_rec(mask) = I(mask);
h = [0.01, 0.08, 0.01; 0.08, 0.64, 0.08; 0.01, 0.08, 0.01];
I_rec = imfilter(I_rec,h,'symmetric');

%smoothing with box filter
I_rec(mask) = I(mask);
h = [0, 0.125, 0; 0.125, 0.5000, 0.1250; 0, 0.125, 0];
I_rec = imfilter(I_rec,h,'symmetric');

I_rec(mask) = I(mask);

end

---

I=imread('lena_destroyed.png');
mask=7
re=notSoBaselineAnymoreFast(I, mask);
imshow(re)

というプログラムを作成しました。

画像を表示した結果真っ黒な画像が出てきました。
maskをどのようにしたら画像が修復されますか。
ご回答よろしくお願い致します

Answer 1

filter処理する際に、画像をdouble型にして渡していますが、imshowはuint8なら輝度表示範囲がデフォルトで[0 255]でああるのに対し、それ以外ではデフォルトで[0 1]となってしまい、真っ黒に表示されます。

imshow(re,[])とすれば表示上はreの最小値を真っ黒、reの最大値を真っ白に表示するので、真っ黒に表示されるのを直せると思います。

もしもこれをuint8に直し保存したいなら、reを[0 255]に収まるように正規化し手からuint8に直せば修復できます。
そのフィルタ係数ですと、単にuint8にするだけでは真っ黒になるので注意が必要です。

Answer 2

imshow(re)を
imagesc(re)に置き換える