利用激光扫描点云测量墙体夹角的方法及Matlab实现咨询

完全可以通过点云处理实现墙体夹角核验，适配你正在使用的Matlab环境，核心逻辑是先拟合两面墙体的平面，再通过平面法向量计算夹角，具体操作步骤如下：

1. 点云预处理 #

首先需要从原始点云中分离出待检测的两面墙体点云：

• 读入点云：调用pcread()函数读取扫描仪导出的pcd/ply格式点云文件
• 去噪处理：调用pcdenoise()移除扫描产生的离群噪声点，避免干扰后续平面拟合
• 墙体分割：优先用pcfitplane()或者pcsegdist函数自动分割两面墙体的点云，分别存储为ptCloud1和ptCloud2两个点云对象

提示：如果点云中包含地面、家具等无关元素，可以先拟合地面平面后调用pcremove()移除干扰点，再做墙体分割精度更高。如果自动分割效果不符合预期，也可以直接使用Matlab自带的Point Cloud Labeler工具手动框选两面墙体的点云，导出为单独的点云对象即可，更适合新手操作。

2. 拟合墙体平面 #

对两个墙体点云子集分别拟合平面，获取平面法向量，示例代码如下：

% maxDistance为点到拟合平面的最大允许距离，室内扫描通常设为0.01~0.05m，可根据扫描精度调整
maxDistance = 0.02;
% 拟合第一面墙体的平面
[model1, ~] = pcfitplane(ptCloud1, maxDistance);
% 拟合第二面墙体的平面
[model2, ~] = pcfitplane(ptCloud2, maxDistance);

返回结果中的model1.Normal、model2.Normal即为两个墙体平面的单位法向量。

3. 计算墙体夹角 #

两个平面的夹角等于其法向量夹角或其补角，取和实际墙体角度一致的数值即可，计算代码如下：

% 计算两个法向量的点积
dot_val = dot(model1.Normal, model2.Normal);
% 计算夹角（弧度制）
angle_rad = acos(abs(dot_val));
% 转换为角度制
angle_deg = rad2deg(angle_rad);

如果需要核验是否为90度，只需判断angle_deg与90的差值是否在你论文设定的误差允许范围内即可。

结果验证 #

你可以通过Matlab可视化功能验证拟合和计算结果是否准确：调用pcshow()显示原始点云，用hold on叠加显示两个拟合平面，确认平面贴合墙体点云，也可将计算得到的夹角数值标注在图中，直接作为论文结果素材使用。

内容的提问来源于stack exchange，提问作者Kzzt