中心确定方法、点云滤波与截面类型对激光扫描树干体积估算精度的影响

《Sensors》：Effect of Center Determination Method, Point Cloud Filtering and Cross-Section Type on the Accuracy of Tree Stem Volume Estimation Using Laser Scanning

　　

静态激光扫描是一种先进的木材体积估算方法，其精度现已达到或超过传统林业体积估算方法，使其能够在木材加工管理中实际应用。本研究评估了点云估算树干体积的精度如何受到树干表面状况（带树皮和无树皮）、截面面积计算方法、点云滤波以及树干中心确定的影响。结果表明，树干表面状况从根本上控制了系统误差的方向和大小。对于带树皮的树干，使用直径平方计算的截面面积导致系统性的体积低估，滤波后低估加剧，而使用均方根直径则导致高估，滤波后高估大幅降低并接近计算机断层扫描（CT）获得的参考值。对于去皮树干，平方直径法导致系统性的高估，滤波后有效缓解，而均方根直径则导致持续的高估。体积法本身对精度只有次要影响。体积估算精度主要由截面面积公式和树干表面状况决定。对于带树皮的原木，算术平均半径平方产生系统性的体积低估，滤波后增加约1%，而均方根半径滤波后高估约2%。对于去皮原木，使用算术平均半径平方时，滤波后体积高估从4.05–4.72%降至1.57–1.80%，而使用均方根半径时高估仍然很高（6.55–6.92%）。相比之下，体积法和中心确定的选择对体积估算的影响小于1%。这些发现证实了静态激光扫描的巨大潜力，并强调了需要适应点云的计算模型。

**论文解读：点云滤波与截面公式对静态激光扫描树干体积估算精度的影响**

**研究背景与问题**
森林生态系统在提供自然资源和生态服务中具有基础性作用，精确评估林分特征是森林管理规划、火灾建模和碳储量估算的关键前提。传统林业测量方法虽然可靠，但成本高、劳动密集。现代技术如地面激光扫描（TLS）和计算机断层扫描（CT）能快速获取数据，但TLS受遮挡、噪声和复杂工作流程限制，CT则因设备昂贵在实际应用中受限。树干形状无法用单一经典几何体建模，原木也常呈现不规则锥度，传统林业测量方法固有估计误差。因此，研究人员评估了点云中心确定方法、点云滤波和截面类型对静态激光扫描树干体积估算精度的影响，并对比有皮与去皮原木，以期找到适用性强、精度高的体积估算流程。该研究发表在《Sensors》期刊。

**研究内容与结论**
研究人员从木材场选取13根原木（含6个树种：泡桐、栎树、鹅耳枥、挪威云杉、欧洲赤松、欧洲冷杉），分别进行CT扫描获取参考体积，并用Gocator激光轮廓仪获取点云数据。通过比较两个中心确定方法（扫描设备直接定义中心、基于主成分分析PCA估算中心）、两个截面面积计算方式（算术平均半径平方、均方根半径）、两种体积积分法（Huber法、Frustum法），以及点云滤波（Open3D库的统计离群值去除）的影响，得出以下结论：截面面积公式主导系统误差方向，滤波和中心确定起校正作用；有皮原木用算术平均半径平方导致低估，用均方根半径导致高估，高估经滤波后接近CT参考值；去皮原木用算术平均半径平方导致高估，滤波后高估值显著降低（从4.05–4.72%降至1.57–1.80%），而均方根半径导致持续高估（6.55–6.92%）；体积积分法和中心确定对体积估算影响小于1%。该研究证实静态激光扫描在木材体积估算中的高精度潜力，强调需开发适应点云特性的计算模型。

**关键技术方法（不超过250字）**
本研究使用的主要技术方法包括：（1）数据采集：采用Microtec三维CT扫描仪（意大利）获取参考体积（纵向分辨率10 mm，横向2×2 mm），采用Gocator 2670激光轮廓仪（三个传感器同步捕获）获取点云数据（空间分辨率1 mm）；（2）中心确定：基于扫描设备三束激光交点定义直接中心，以及基于主成分分析（PCA）估算质心作为中心；（3）截面面积计算：采用算术平均半径平方公式（S = π(R?)2）及均方根半径公式（S = π(R²)?）两种简化方法；（4）点云滤波：利用Python环境中的Open3D库执行统计离群值去除（邻域点数30，标准差比率2）；（5）体积计算：采用Huber法（基于中点截面面积）和Frustum法（线性锥度模型），将原木分割为厚度2 cm的切片，汇总切片体积得到总体积。样本来源为木材场存储的原木（13根，6个树种，有皮与去皮）。

**研究结果**
**带树皮原木的体积估算**：
通过比较Huber法与Frustum法，以及滤波前后的体积估算与CT参考体积，发现使用算术平均半径平方计算截面面积时，所有树种的体积均系统性地低估，滤波后低估增加约1%；使用均方根半径时，体积系统性地高估，滤波后高估降至约2%，更接近CT参考值。无论滤波与否或中心确定方法如何，各树种的相对排名保持不变，Huber法与Frustum法的差异小于1%，中心确定（PCA）导致体积一致减小约0.4%。

**去皮原木的体积估算**：
在使用算术平均半径平方时，滤波前体积系统高估4.05–4.72%，滤波后高估降至1.57–1.80%；使用均方根半径时，滤波前高估6.64–7.31%，滤波后仍高达6.55–6.92%。Huber法与Frustum法差异小于1%，PCA中心确定导致体积一致减小约0.3%。

**体积积分法与中心确定的影响**：
Huber法与Frustum法对最终体积估算的影响最小（差异<1%），因为原木被分为100个切片，几何失真效应被衰减。中心确定方法（PCA与设备直接中心）影响较小（约0.3–0.4%），但方向保持一致，表明PCA轻微削弱局部偏心，导致有效直径减小。

**讨论与结论总结**
讨论部分指出，截面面积公式是系统误差的主导因素，而点云滤波和中心确定为校正因素。有皮与去皮的相反误差方向源于树皮空间变厚导致的半径分布差异。对于点云滤波，统计离群值去除可能优先移除远离轴心点，导致几何收缩、体积低估。最佳组合（均方根半径结合Open3D滤波）使带皮原木体积估算偏差最小。未来研究应扩大样本量、分析更多树种及点云密度影响。
结论翻译：静态激光扫描是一种快速、精确且成本低的锯木厂木材体积测量现代方法，其精度已超越传统林业方法。本研究表明，树皮存在显著影响系统误差行为。对于带皮原木，算术平均半径平方导致低估（滤波后加重），均方根半径导致高估（可被滤波大幅降低）；对于去皮原木，算术平均半径平方引起的高估可通过滤波有效缓解，而均方根半径的高估由模型误差主导且持续存在。体积积分法（Huber或Frustum）和中心确定对精度影响次要，截面面积计算方法和树干表面状况是主要影响因素。未来应基于点云推导新体积方程，并进一步研究点云质量与计算模型对精度的交互作用。