bw = imread('circles.png');                        % 讀入圖像
L=bwlabel(bw);                                     % 圖像標記
s=regionprops(L, 'centroid');                      % 求區域屬性   
a=regionprops(L, 'area')                           % 求區域屬性
centroids = cat(1, s.Centroid);                    
imshow(bw)                                         % 顯示原始圖像 
hold on                                            % 繼續繪圖
plot(centroids(:,1), centroids(:,2), 'b*')         % 質心標識，藍色星號
hold off                                           % 停止繼續繪圖