function imo = HON_cap(im)
%

[M,N,P] = size(im);

imo = zeros(M,N);

for r = 1:floor(M/4)
    for c = 1:N
        if im(r,c,1)>2.6*im(r,c,2)&im(r,c,1)>2.6*im(r,c,3)
            imo(r,c) = 1;
        end
    end
end

tch = 0;
