BwlabelRegionpropsSubplot

function JajittelePrikat()
%kommentoi käskyt Matlab helpin avulla. kokeile erilaisia diskin halkaisija
%muuta subplot käsky siten, että tulostetaan 2 riviä ja  kolme saraketta
%
I = imread('PrikkojaJaNappeja.bmp');
%imfinfo('PrikkojaJaNappeja.bmp')
%Täytetään kuvasta pienet reiät, koeta kommentoida kaksi seuraavaa riviä
%mikä on tulos ?
%kokeile myös suurentaa "diskin" kokoa esim. 10, mikä on tulos ?
se = strel('disk',3);
I=imdilate(I,se);%suljetaan roskia eli pieniä reikiä kuvasta myös imclose toimii !
%Haetaan osat kuvasta
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);
reikia_2 = find([ominaisuus.EulerNumber]==-1);
reikia_3 = find([ominaisuus.EulerNumber]==-2);
%sijoitetaan osat reikien mää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');