matlab示例：从点云重构曲线

1、命令模块一：

% Create a set of points representing the point cloud

numpts = 192;

t = linspace( -pi, pi, numpts+1 )';

t(end) = [];

r = 0.1 + 5*sqrt( cos( 6*t ).^2 + (0.7).^2 );

x = r.*cos(t);

y = r.*sin(t);

ri = randperm(numpts);

x = x(ri);

y = y(ri);

2、命令模块二：

% Construct a Delaunay Triangulation of the point set.

dt = delaunayTriangulation(x,y);

tri = dt(:,:);

3、命令模块三：

% Insert the location of the Voronoi vertices into the existing

% triangulation

V = dt.voronoiDiagram();

% Remove the infinite vertex

V(1,:) = [];

numv = size(V,1);

dt.Points(end+(1:numv),:) = V;

4、警告：已检测到重复的数据点并将其删除。

相对于delaunayTriangulation属性X中的唯一点集定义了Triangulation索引。

5、命令模块四：

% The Delaunay edges that connect pairs of sample points represent the

% boundary.

delEdges = dt.edges();

validx = delEdges(:,1) <= numpts;

validy = delEdges(:,2) <= numpts;

boundaryEdges = delEdges((validx & validy), :)';

xb = x(boundaryEdges);

yb = y(boundaryEdges);

clf;

triplot(tri,x,y);

axis equal;

hold on;

plot(x,y,'*r');

plot(xb,yb, '-r');

xlabel('Curve reconstruction from a point cloud', 'fontweight','b');

hold off;

6、按“Enter”键。

如图1所示。