reikien taytto

 I = imread('PrikkojaJaNappeja.bmp');
%imfinfo('PrikkojaJaNappeja.bmp')
%täyteään kuassa mahdollisesti olevat reiät
%kokeile suurentaa "diskin" kokoa esim. 10, mikä on tulos ?
se = strel('disk',3);
I=imdilate(I,se);%suljetaan roskia eli pieniä reikiä, kokeile kuvasta myös imclose funktiota!
%Numeroidaan kuvan osat
L = bwlabel(I);
RGB1 = label2rgb(L,'jet');%tehdään muuten vaan värikuva
%Haetaan osien ominaisuudet EulerNumber = kohde-reiät
ominaisuus = regionprops(L,'EulerNumber');
reikia_1 = find([ominaisuus.EulerNumber]==0);% yksireikä löytyi
reikia_2 = find([ominaisuus.EulerNumber]==-1);
reikia_3 = find([ominaisuus.EulerNumber]==-2);
%sijoitetaan osat reikien maäärän mukaan uusiin kuviin
I1 = ismember(L,reikia_1);
I2 = ismember(L,reikia_2);
I3 = ismember(L,reikia_3);
%Tehdään 24-bittinen RGB-kuva, jossa eri osat omilla väreillään
%cat(värien määrä,punainen, vihreä ,sininen)
%kerrotaan 255:llä, jotta saadaan oikea alue 0.0->1.0 to 0->255
RGB3=uint8(cat(3,I1,I2,I3))*255;
subplot(3,2,1),imshow(I), title('alkuperäinen');
subplot(3,2,2),imshow(I1), title('1-reikäiset');
subplot(3,2,3),imshow(I2), title('2-reikäiset');
subplot(3,2,4),imshow(I3), title('3-reikäiset');
subplot(3,2,5),imshow(RGB1), title('jet');
subplot(3,2,6),imshow(RGB3), title('eri värit');