[论文解读] Least squares fitting of circles and lines

N. Chernov, C. Lesort|ArXiv.org|Jan 1, 2003

Image and Object Detection Techniques参考文献 32被引用 41

一句话总结

本文对圆及圆弧的最小二乘拟合进行了全面分析，提出了一种新型迭代算法（LMA），在可靠性与效率方面优于现有方法，尤其适用于小段圆弧。研究评估了代数拟合与几何拟合技术，表明将先进的代数预拟合（AF2/AF3）与Levenberg-Marquardt算法结合，可在各种数据配置下实现更优的收敛性与鲁棒性。

ABSTRACT

We study theoretical and computational aspects of the least squares fit (LSF) of circles and circular arcs. First we discuss the existence and uniqueness of LSF and various parametrization schemes. Then we evaluate several popular circle fitting algorithms and propose a new one that surpasses the existing methods in reliability. We also discuss and compare direct (algebraic) circle fits.

研究动机与目标

• 建立最小二乘法拟合圆的理论基础，包括解的存在性与唯一性。
• 评估并比较现有几何拟合算法（Levenberg-Marquardt、Landau、Späth）在鲁棒性与计算成本方面的表现。
• 开发并验证一种新型迭代算法（LMA），以提升收敛的可靠性与效率。
• 评估代数拟合方法（AF1、AF2、AF3、TRI、CEN）作为迭代优化预处理步骤的性能。
• 为涉及噪声或稀疏数据的应用（如计算机视觉与图像处理）提供一种实用且高性能的圆拟合解决方案。

提出的方法

• 提出一种基于Levenberg-Marquardt优化的新迭代算法（LMA），采用改进的参数更新策略以增强收敛稳定性。
• 使用代数预拟合（AF2与AF3）为迭代优化过程生成高质量初始估计。
• 采用几何目标函数 F = Σd_i²，其中 d_i 为每个数据点到拟合圆的欧氏距离。
• 分析圆拟合的参数空间（圆心 (a,b)，半径 R），研究 F 的行为以确保收敛至全局最小值。
• 通过每数据点的浮点运算次数（flops）比较计算成本，以评估算法效率。
• 通过在不同弧长、噪声水平与点分布的合成数据上进行大量数值实验，验证性能。

实验结果

研究问题

• RQ1 在何种条件下，一组点的最小二乘圆拟合存在，且解唯一？
• RQ2 不同迭代算法（Levenberg-Marquardt、Landau、Späth）在小段圆弧上的收敛可靠性与计算成本表现如何？
• RQ3 代数拟合方法（AF2、AF3）能否作为更优的预拟合步骤，以提升迭代几何拟合算法的收敛性？
• RQ4 数据配置（如共线点或对称点集）对最小二乘拟合的存在性与唯一性有何影响？
• RQ5 在不同弧长与噪声水平下，所提出的LMA算法在收敛概率与计算成本方面相较于现有方法表现如何？

主要发现

• 所提出的LMA算法在小段圆弧（5°–15°）上实现了最高的收敛概率（接近100%），优于LMC及其他迭代方案。
• AF2或AF3代数预拟合与LMA算法结合，实现了更优的可靠性与效率，所有测试配置下的收敛概率均稳定高于95%。
• 对于弧长低于30°的情况，Späth与Landau方法频繁出现不收敛或超过2000次flops成本阈值，因此在小弧段上不实用。
• 基于Levenberg-Marquardt的算法（LMC与LMA）即使在5°弧段上也保持稳定收敛，平均每次拟合迭代19–20次，尽管难度逐渐增加。
• AF2与AF3预拟合方法显著提升了迭代算法的收敛速度与成功率，在可靠性与效率方面优于AF1、TRI与CEN。
• LMA的计算成本保持适中（每数据点低于1000次flops），同时保持高可靠性，适用于实时或大规模应用场景。