function [res, cpt] = squelette(I)
    res=soustraction(I, ouverture(I, 1));
    ero=erosion(I, 1);
    cpt=1;
    while(max(max(ero))>0)
        res=addition(res, soustraction(ero, ouverture(ero, 1)));
        ero=erosion(ero, 1);
        cpt=cpt+1;
    end
end