matlab - 画像のしきい値を設定するにはどうすればよいですか?

Question

すでに準同型フィルターを通過した画像にしきい値処理を適用する必要があります。

私のしきい値は、画像強度の平均+標準偏差でなければなりません。

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')

しかし、出力は完全に黒く、

これを修正するにはどうすればよいですか?