平面点集下三条共交点直线的拟合方法咨询

2D平面场景求解逻辑 #

问题参数化 #

我们先把约束条件转化为可优化的参数形式：

• 设三条直线的公共交点为 (x0, y0)，任意过该点的直线可表示为 a_i(x - x0) + b_i(y - y0) = 0，其中i=1,2,3分别对应红、蓝、黄三个点集
• 单个点 (x_{i,j}, y_{i,j}) 到对应直线的平方距离为 [a_i(x_{i,j}-x0) + b_i(y_{i,j}-y0)]^2 / (a_i^2 + b_i^2)
• 优化目标为最小化所有点到对应直线的平方距离总和

分阶段求解 #

该问题可以拆分为两个子问题迭代求解，效率很高且保证全局最优：

1. 固定公共交点求解最优直线
   当(x0, y0)固定时，将所有点做坐标偏移得到(u_{i,j}, v_{i,j}) = (x_{i,j}-x0, y_{i,j}-y0)，此时问题转化为对每个偏移后的点集求过原点的最小二乘拟合直线，直接对每个点集的偏移坐标矩阵做SVD分解，最小奇异值对应的右奇异向量就是直线的法向量(a_i, b_i)，同时可直接算出当前(x0, y0)对应的总损失值。
2. 优化公共交点坐标
   此时问题退化为仅含两个变量(x0, y0)的无约束凸优化问题，用Nelder-Mead、梯度下降等常用无约束优化算法就能快速收敛到全局最优解。

3D场景迁移方法 #

该方案可以完全无缝迁移到你提到的三维空间共交线平面拟合场景：

• 先将公共交线参数化，做坐标变换把公共交线对齐到Z轴
• 对每个点集做坐标变换后，求过Z轴的最优拟合平面，同样可以用SVD直接求解
• 再优化公共交线的参数即可，整体逻辑和2D场景完全一致

轻量落地替代方案 #

如果不想实现迭代优化，也可以用两步闭式方案实现，精度损失极低且实现难度极小：

1. 先去掉共交点约束，分别对三个点集做标准最小二乘直线拟合，得到三条初始直线
2. 求解到三条初始直线平方距离和最小的点，作为近似公共交点(x0, y0)，该问题有闭式解无需迭代
3. 以(x0, y0)为公共交点，分别对三个点集求过该点的最小二乘拟合直线即可

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