我需要对已经通过同态滤波器的图像应用阈值。
我的阈值必须是图像强度的平均值+标准差。
我使用Jan Motl 的阈值代码如下:
function J = bernsen_thres(I)
T = thres_val(I);
J = bernsen(I, [T T],20,'replicate');
end
function T = thres_val(I)
mn = mean(I(:));
sd = std(double(I(:)));
thres = round((mn+sd));
if(is_odd(thres))
T = thres;
else
T = thres+1;
end
function ret = is_odd(val)
if(mod(val,2) == 0);
ret = 0;
else
ret = 1;
end
我使用了来自 Steve Eddins 的同态过滤器代码,如下所示,
clear_all();
I = gray_imread('cameraman.png');
I = steve_homo_filter(I);
Ithres = bernsen_thres(I);
imshowpair(I, Ithres, 'montage')
但是输出完全是黑色的,
我应该怎么做才能解决这个问题?