python 拟合圆心_圆的散点拟合, 已知圆的采样点, 求圆的圆心和半径.

UTF8gbsn

本文介绍一个圆的拟合方法. 参照论文为 Coope, I. D. (1993). Circle fitting

by linear and nonlinear least squares. Journal of Optimization Theory

and Applications, 76(2), 381–388. https://doi.org/10.1007/BF00939613

问题描述

首先还是先来描述一下问题,

假如我们有一个采样数据A={a1,a2,⋯ ,aN}A=\{a_1,a_2, \cdots, a_N\}A={a1​,a2​,⋯,aN​}, 其中aia_iai​,

是n维空间中点的坐标向量. 我们希望去拟合一个圆心xxx, 和一个半径rrr,

注意xxx是一个n维向量而rrr是一个标量. 注意重申一下我们的已知集合AAA,

未知x,rx, rx,r.

求解过程

首先, 我们需要写出我们的最优化化方程

arg⁡min⁡x,r=∑i=1N[∥x−ai∥22−r2]2\arg\min_{x,r}=\sum_{i=1}^{N}[\|x-a_i\|_2^2-r^2]^2argx,rmin​=i=1∑N​[∥x−ai​∥22​−r2]2

接下来我们来分析下中括号内部的这个东西怎么来写.

∥x−ai∥22−r2=(x−ai)T(x−ai)−r2=(xT−aiT)(x−ai)−r2=xTx−2aiTx+aiTai−r2\left.

\begin{aligned}

\|x-a_i\|_2^2-r^2&=(x-a_i)^T(x-a_i)-r^2\\

&=(x^T-a_i^T)(x-a_i)-r^2\\

&=x^