function imo = HON_label(im)
%

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

imo = zeros(M,N);

for r = floor(M/2):M
    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
imo = bwmorph(bwmorph(bwmorph(imo,'erode',3),'dilate',21),'erode',12);

tch = 0;
