thinning.m

matlab代码

function out = thinning(in);


%update on 2/12

%thinning process obtain a one pixel wide image skeleton
% use function n_sum and t_sum

[w,h] = size(in);

out = in;

for i= 3:h-2
for j= 3:w-2

        if in(i,j) == 1

                if 1 < n_sum(i,j,in) &  n_sum(i,j,in) < 7 
                        if  t_sum(i,j,in) == 2
                                if or ( and_157(i,j,in) == 0  , t_sum(i,j-1,in) ~= 2 )
                                        if or (and_357(i,j,in) == 0  , t_sum(i+1,j,in) ~=2 )                              
                                                out(i,j) = 0;
                                        else
                                                out(i,j) = 1;
                                        end;
                                else
                                        out(i,j) = 1;
                                end;
                        else
                                out(i,j) = 1;

                        end;
                end;

        end;                

end;
end;