Kuvan osien numerointi

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ään24-bittinen RGB-kuva, jossa eri osat omilla väreillään

%cat(värien määrä,punainen, vihreä ,sininen)

%kerrotaan255:llä, jotta saadaan oikea alue0.0->1.0 to0->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');