keeplocalmax.m

image thresholding otsu method

function im2 = keeplocalmax(im)

ok = 0*im;

for v = -1:1
    for u = -1:1
        if ((u ~= 0) | (v ~= 0))
            ok = ok - (im<immove(im,v,u));
        end
    end
end

im2 = im .* (ok==0);