Euler ja keskipiste

function f = HaePrikatJaMerkitseKeskipiste()
I = imread('PrikkojaJaLevyja.bmp');
imshow(I), title('alkuperainen kuva');
L = bwlabel(I);%numeroi alueet
Ominaisuus = regionprops(L,'all')
idx = find([Ominaisuus.EulerNumber] == 0 ); % 0=prikka 1=levy
%tee uusi kuva, jossa vain prikat tai levyt
I2 = ismember(L,idx);
figure,imshow(I2), title('Haetaan Eulerin numerolla prikat tai levyt ja merkitään massakeskipisteet');
%numeroi uudet alueet
NumeroidutAlueet=bwlabel(I2);
OminaisuusI2 = regionprops(NumeroidutAlueet,'all');
hold on %pidä kuva "auki" lisäystä varten
Maara=max(NumeroidutAlueet(:))%hae kpl määrä kuvasta
%haetaan kaikkien kuvan objektien keskipisteet
for i=1:Maara
    Keskipiste=OminaisuusI2(i).Centroid
    plot(Keskipiste(1),Keskipiste(2),'X')
end
hold off