lb = bwlabel( inputBW, 4 );
 st = regionprops( lb, 'Centroid', 'Area' );
 [bgArea toRemove] = max( [st(:).Area] );
 st(toRemove) = []; %// discard background
 cnt = vertcat( st.Centroid ); %// all centers

figure;
 imshow( inputBW ); hold on;
 scatter( cnt(:,1), cnt(:,2), 50, 'xg' );