[size=large]
clc;
clear all;
close all;
workspace;
fontSize = 24;
originalImage = peaks(256);
[rows columns] = size(originalImage);
subplot(2,2,1);
imshow(originalImage, []);
title('Original "Peaks" Image', 'FontSize', fontSize);
set(gcf, 'Position', get(0, 'ScreenSize')); % Maximize figure.
% Have user locate points on the image.
binaryImage = originalImage > 5;
subplot(2,2,2);
imshow(binaryImage);
title('Binary Image', 'FontSize', fontSize);


% Label each blob so we can make measurements of it
[labeledImage numberOfBlobs] = bwlabel(binaryImage, 8);
blobMeasurements = regionprops(labeledImage, originalImage,
'Centroid');
xCentroid = round(blobMeasurements(1).Centroid(1))
yCentroid = round(blobMeasurements(1).Centroid(2))
% Make image to hold centroid.
centroidImage = false(rows, columns);
centroidImage(yCentroid, xCentroid) = 1;
subplot(2,2,3);
imshow(centroidImage);
[size=xx-large]title('Centroid Image', 'FontSize', fontSize);

% Calculate the distance of every pixel in the image
% from the centroid of the blob.
edtImage = bwdist(centroidImage);
subplot(2,2,4);
imshow(edtImage, []);
title('Distance Transform', 'FontSize', fontSize);

1, 1/2*3, 1/2*6, 1/2*9, 1*2*12, ...