Otra posibilidad es utilizar la función BWBOUNDARIES, que:
traza los contornos exteriores de los objetos, así como límites de agujeros dentro de estos objetos
Esa información está contenida en el cuarto resultado A
, una matriz de adyacencia que representa las dependencias padre-hijo-agujero.
%# read binary image
bw = imread('SUvif.png');
%# find all boundaries
[B,L,N,A] = bwboundaries(bw, 8, 'holes');
%# exclude inner holes
[r,~] = find(A(:,N+1:end)); %# find inner boundaries that enclose stuff
[rr,~] = find(A(:,r)); %# stuff they enclose
idx = setdiff(1:numel(B), [r(:);rr(:)]); %# exclude both
bw2 = ismember(L,idx); %# filled image
%# compare results
subplot(311), imshow(bw), title('original')
subplot(312), imshow(imfill(bw,'holes')), title('imfill')
subplot(313), imshow(bw2), title('bwboundaries')
solución impresionante. Muchas gracias! ¡Lo único que cambié fue eliminar objetos cuyo espacio es menor que el 2% de la imagen y funciona muy bien! – user1240792