algoritmopavlidis.m

algoritmo pavilis, para el adelgazamiento de imagenes

function y = algoritmoPavlidis(imagenBinaria)
% Funcion que se encarga de esqueletizar una imagen por adelgazamientos
% iterativos

% en las plantillas definimos a cada valor A con -3, a cada valor B con -2
% y a cada valor de p con -1; para poder representarlo en la matriz
% el elemento p es 1 siempre
Plantillas(:,:,1) = [-3 -3 -3;  0  1  0; -2 -2 -2];
Plantillas(:,:,2) = [-3  0 -2; -3  1 -2; -3  0 -2];
Plantillas(:,:,3) = [-2 -2 -2;  0  1  0; -3 -3 -3];
Plantillas(:,:,4) = [-2  0 -3; -2  1 -3; -2  0 -3];
Plantillas(:,:,5) = [-3 -3 -3; -3  1  0; -3  0  2];
Plantillas(:,:,6) = [-3  0  2; -3  1  0; -3 -3 -3];
Plantillas(:,:,7) = [ 2  0 -3;  0  1 -3; -3 -3 -3];
Plantillas(:,:,8) = [-3 -3 -3;  0  1 -3;  2  0 -3];
Plantillas(:,:,9) = [-3 -3 -3; -3  1  0; -3  0  1];
Plantillas(:,:,10) = [-3  0  1; -3  1  0; -3 -3 -3];
Plantillas(:,:,11) = [ 1  0 -3;  0  1 -3; -3 -3 -3];
Plantillas(:,:,12) = [-3 -3 -3;  0  1 -3;  1  0 -3];

REMAIN = 1; % se usa para denotar la existencia de pixeles esqueletales
SKEL = 0;   % se fija a cierta cuando alguna vecindad de un pixel es similar a
% alguna de las plantillas existentes en P
imagen = convertirCeroUno(imagenBinaria);
img = redimensionarImagen(imagen,3);
[nf nc] = size(img);
d = 1;% (3-1)/2
while(REMAIN == 1)
    REMAIN = 0;
    for k=1:4   % los indices de los vecinos de un determinado punto
        % se recorre toda la imagen
        for i=1+d:nf-d
            for j=1+d:nc-d
                mat = extraerMatrizImagen(img,3,i,j);
                vecino = obtenerVecino(k,mat);
                if vecino == 0 && img(i,j)==1                    
                    SKEL = compararMatrizPlantillas(mat,Plantillas);%compara si se asemeja a la plantilla
                    if SKEL == 1% si es cierto y la plantilla concuerda con la matriz
                        img(i,j) = 2;
                    else
                        img(i,j) = 3;
                        REMAIN = 1;
                    end
                end
                %%%%%%%%%%%%%%%%%%%%%%%%%%%%
            end
        end
        % se reemplazan los 3 por 0
        for i=1:nf
            for j=1:nc
                if img(i,j)==3
                    img(i,j)=0;
                end
            end
        end
    end
end % fin del algoritmo
imgAux = redimInvImage(img,1,1,1,1);
imagenEsqueletizada = desconvertirImagen(imgAux);
subplot(1,2,1),imshow(uint8(imagenBinaria)),title('Imagen Original');
subplot(1,2,2),imshow(uint8(imagenEsqueletizada)),title('Imagen Esqueletizada');% fin del algoritmo
y = imagenEsqueletizada;