Tristan Ursell
Distance Transform for Arbitrary Point Cloud
March 2012
[Xvec,Yvec,Dist]=clouddist(Xin,Yin,xlim,ylim,res);
The function creates a Euclidean distance transform from an arbitrary set of points, which is the inplane distance to to closest point. This does not require the Image Processing Toolbox, and in fact, can be used as a more versatile version of 'bwdist' transform.
(Xin, Yin) are a list of input points from which to calculate the distance transform.
'xlim' and 'ylim' are each two component vectors that specify the extent of the distance transform calculation.
'res' sets the spatial resolution of the distance transform. For instance is the image goes from 1,1 in both X and Y, you might set res = 0.01, giving you 201 points along each dimension. Generally, good practice is the have 'res' set to be much less than the closest distance between any two points.
The output vectors 'Xvec' and 'Yvec' specify the exact points at which the distance transform 'Dist' has been calculated.
see also: bwdist
Example: general points in 2D
Xin=rand(1,50);
Yin=rand(1,50);
xlim=[1/2,3/2];
ylim=[1/2,3/2];
res=0.01;
[Xvec,Yvec,Dist]=clouddist(Xin,Yin,xlim,ylim,res);
figure;
hold on
imagesc(Xvec,Yvec,Dist)
plot(Xin,Yin,'wx')
title('Point Cloud Distance Transform')
xlabel('X')
ylabel('Y')
axis equal tight
box on
Example: mimicing the bwdist transform
Xin=ceil(249*rand(1,50));
Yin=ceil(199*rand(1,50));
xlim=[1,250];
ylim=[1,200];
res=1;
[Xvec,Yvec,Dist]=clouddist(Xin,Yin,xlim,ylim,res);
figure;
hold on
imagesc(Xvec,Yvec,Dist)
plot(Xin,Yin,'wx')
title('Point Cloud Distance Transform')
xlabel('X')
ylabel('Y')
axis equal tight
box on
