最大角度法检测点云边界的原理，以及其对应的matlab代码，代码中全为matlab自带的内置函数

最大角度法是一种常见的点云边界检测方法，其基本原理是通过计算每个点的法向量与邻域点法向量的夹角，找到最大的夹角作为该点的边界特征。

具体的实现步骤如下：

1. 选定一个点作为起始点，并计算其法向量。

2. 在该点周围选取一定数量的邻域点（可以使用matlab自带的kdtree或balltree等算法进行快速查询）。

3. 计算每个邻域点的法向量。

4. 计算该点的法向量与每个邻域点法向量的夹角，找到夹角最大的邻域点。

5. 如果夹角大于一定阈值，则该点为边界点，否则将其标记为非边界点。

6. 将被标记为非边界点的邻域点和该点从点云中移除，重复以上步骤，直到所有点都被标记为边界点或非边界点。

下面是使用matlab实现最大角度法检测点云边界的代码示例：

% 读取点云数据
ptCloud = pcread('pointCloud.ply');

% 定义邻域搜索半径和夹角阈值
searchRadius = 0.1;
angleThreshold = 30;

% 构建kdtree
model = createns(ptCloud.Location,'NSMethod','kdtree');

% 初始化点云标记矩阵
isBoundary = zeros(ptCloud.Count,1);

% 选取起始点
startIndex = 1;
isBoundary(startIndex) = true;

while true
    % 获取起始点
    startPoint = ptCloud.Location(startIndex,:);
    
    % 搜索邻域点
    [neighborIndices,~] = rangesearch(model,startPoint,searchRadius);
    neighborIndices = neighborIndices{1};
    neighborIndices(neighborIndices == startIndex) = [];
    
    % 计算法向量
    startPointNormal = pcnormals(ptCloud, startIndex);
    neighborNormals = pcnormals(ptCloud, neighborIndices);
    
    % 计算夹角
    angles = acosd(abs(dot(startPointNormal, neighborNormals, 2)));
    
    % 找到最大夹角
    [maxAngle, maxIndex] = max(angles);
    
    % 判断是否为边界点
    if maxAngle > angleThreshold
        isBoundary(neighborIndices(maxIndex)) = true;
        startIndex = neighborIndices(maxIndex);
    else
        break;
    end
end

% 可视化结果
color = uint8([255, 0, 0]);
color(isBoundary,:) = repmat(uint8([0, 255, 0]), sum(isBoundary), 1);
ptCloud.Color = color;
pcshow(ptCloud);

在以上代码中，我们首先读取点云数据，并定义了邻域搜索半径和夹角阈值。然后使用matlab内置的`createns`函数构建kdtree用于快速搜索邻域点。接着，我们初始化点云标记矩阵，并选取起始点开始搜索边界点。在每次搜索中，我们首先获取起始点，并使用`rangesearch`函数查询其邻域点。然后，我们计算起始点和邻域点的法向量，并计算夹角。找到夹角最大的邻域点后，我们判断是否为边界点，并更新起始点。如果夹角小于阈值，则搜索结束。

最后，我们将边界点标记为绿色，非边界点标记为红色，并可视化结果。