文章快速检索     高级检索
  大地测量与地球动力学  2018, Vol. 38 Issue (11): 1191-1195  DOI: 10.14075/j.jgg.2018.11.018

引用本文  

邱德超, 鲁铁定, 邓小渊. 一种基于Partial EIV模型的圆曲线拟合解法[J]. 大地测量与地球动力学, 2018, 38(11): 1191-1195.
QIU Dechao, LU Tieding, DENG Xiaoyuan. A Circular Curve Fitting Solution Based on Partial EIV Model[J]. Journal of Geodesy and Geodynamics, 2018, 38(11): 1191-1195.

项目来源

国家自然科学基金(41374007,41464001);江西省科技落地计划(KJLD12077);江西省教育厅科技项目(GJJ13457);江西省自然科学基金(2017BAB203032);国家重点研发计划(2016YFB0501405,2016YFB0502601-04)。

Foundation support

National Natural Science Foundation of China, No. 41374007, 41464001; Jiangxi Province Science and Technology Landing Project, No. KJLD12077;Jiangxi Provincial Department of Education Science and Technology Project, No. GJJ13457; Natural Science Foundation of Jiangxi Province, No. 2017BAB203032; National Key Research and Development Program of China, No. 2016YFB0501405, 2016YFB0502601-04.

通讯作者

鲁铁定,教授,主要从事测绘数据处理研究,E-mail:tdlu@whu.edu.cn

第一作者简介

邱德超,硕士生,主要从事现代大地测量数据处理研究,E-mail:765730354@qq.com

About the first author

QIU Dechao, postgraduate, majors in modern geodetic data processing, E-mail: 765730354@qq.com.

文章历史

收稿日期:2017-11-23
一种基于Partial EIV模型的圆曲线拟合解法
邱德超1,2,3     鲁铁定1,2     邓小渊4     
1. 东华理工大学测绘工程学院,南昌市广兰大道418号, 330013;
2. 流域生态与地理环境监测国家测绘地理信息局重点实验室,南昌市广兰大道418号, 33001;
3. 3 云浮市国土资源和城乡规划管理局,广东省云浮市云城区府前路11号,527300;
4. 浙江省地理信息中心,杭州市保俶北路83号,310012
摘要:针对圆曲线拟合问题,以圆曲线的参数方程为基础建立圆曲线拟合的EIV模型,根据系数矩阵的特点将模型转化为更合理的Partial EIV模型,通过公式变形为最小二乘形式,采用两步迭代法求解模型参数,保证系数矩阵中相同元素的改正数一致,常数元素的改正数为零。算例数据结果表明,所提算法的可行性、拟合精度相对较优。
关键词圆曲线拟合参数方程总体最小二乘Partial EIV模型

针对圆曲线拟合问题,许多学者进行了各种算法研究。文献[1]以圆曲线一般方程为基础,采用最小二乘原理求解出圆曲线参数;文献[2]同样以圆曲线一般方程为基础,以原半径与拟合半径的几何距离差值平方和最小为准则,使用泰勒公式展开成最小二乘形式迭代求解;文献[3]以圆曲线正交距离残差平方和最小为拟合准则,基于参数方程提出采用正交最小二乘法拟合圆曲线;文献[4]采用总体最小二乘(total least square,TLS)算法求解待定参数;文献[5]提出利用协方差传播定律定权的加权总体最小二乘(weight total least square, WTLS)算法解算圆拟合参数的思想,但缺少实例;文献[6]在参数方程的拟合模型中采用总体最小二乘算法求解拟合参数,但平差计算量较大;文献[7]以圆曲线参数方程为基础构建拟合模型,采用最小二乘(least square,LS)算法取得较好的效果,但中间参数初始值的确定较为复杂。

针对总体最小二乘算法,文献[8]根据系数矩阵既含非随机元素又含重复的随机元素的特点提出Partial EIV(partial errors-in-variables)模型,该模型将系数矩阵中非随机元素与重复的随机元素分离并进行平差计算,使平差模型更具一般性;文献[9]对Partial EIV模型进行移项变化,转换为间接平差形式,采用两步迭代法求解,简化了参数的表达形式且提高了计算效率;文献[10]概述了整体最小二乘估计的发展历史及研究进展,侧重分析了各种算法的本质特征。

本文首先分析现有的圆曲线模型的构建及其参数求解方法,然后以圆曲线的参数方程为基础,针对模型的系数矩阵特点,建立圆曲线拟合的Partial EIV模型,采用两步迭代法求解模型参数。结合实例数据和模拟数据,对比分析现有的以一般方程构建的圆曲线拟合模型对应的最小二乘算法、EIV模型的总体最小二乘算法。实验结果表明,本文提出的基于Partial EIV模型的总体最小二乘法求解由参数方程构建的拟合模型是可行的。

1 拟合模型构建 1.1 圆曲线一般方程形式

圆曲线一般方程形式为[11]

$ {x^2} - 2ax + {y^2} - 2by + c = 0 $ (1)

式中,(a, b)为圆心,半径$ r = \sqrt {{a^2} + {b^2} - c} $

1.2 圆曲线参数方程的EIV模型

平面上确定圆心(a, b)和半径r,就可以唯一确定圆。假设圆上一点m(x, y), 采用圆的参数表达形式为:

$ \left\{ \begin{array}{l} x = a + r\cos \theta \\ y = b + r\sin \theta \end{array} \right. $ (2)

式中,θ∈[0 2π]是旋转角参数。式(2)采用平差值表示为:

$ \left\{ \begin{array}{l} \hat x = \hat a + \hat r\cos \hat \theta \\ \hat y = \hat b + \hat r\sin \hat \theta \end{array} \right. $ (3)

式中,$\left( {\hat a, \hat b} \right)$为圆心坐标平差值,$ {\hat r}$${\hat \theta }$为半径和旋转角的平差值。令${\hat x}$=x+vx, ${\hat y}$=y+vy, ${\hat r}$=r0+δr, ${\hat \theta }$=θ0+δθ, ${\hat a}$ =a0+δa, ${\hat b}$ =b0+δb,将式(3)线性化后误差方程为:

$ \left\{ \begin{array}{l} {v_x} = \delta a + \cos {\theta ^0}\delta r - {r^0}\sin {\theta ^0}\frac{{\delta \theta }}{\rho } - {L_x}\\ {v_y} = \delta b + \sin {\theta ^0}\delta r + {r^0}\cos {\theta ^0}\frac{{\delta \theta }}{\rho } - {L_y} \end{array} \right. $ (4)

式中,$\left\{ \begin{array}{l} {L_x} = x - ({a^0} + {r^0}{\rm{cos}}{\theta ^0}) = x - {x^0}\\ {L_y} = y - ({b^0} + {r^0}{\rm{sin}}{\theta ^0}) = y - {y^0} \end{array} \right. $θ0x0y0计算如下:

$ \left\{ \begin{array}{l} \cos {\theta ^0} = \frac{{x - a}}{{\sqrt {{{\left( {x - a} \right)}^2} + {{\left( {y - b} \right)}^2}} }}\\ \sin {\theta ^0} = \frac{{y - b}}{{\sqrt {{{\left( {x - a} \right)}^2} + {{\left( {y - b} \right)}^2}} }}\\ {x^0} = a + {r^0}\cos {\theta ^0}\\ {y^0} = b + {r^0}\sin {\theta ^0} \end{array} \right. $ (5)

式(4)写成矩阵形式为:

$ \begin{array}{*{20}{c}} {\left[ {\begin{array}{*{20}{c}} {{v_x}}\\ {{v_y}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} 1&0&{\cos {\theta ^0}}&{ - {r^0}\sin {\theta ^0}/\rho }\\ 0&1&{\sin {\theta ^0}}&{{r^0}\cos {\theta ^0}/\rho } \end{array}} \right] \cdot }\\ {\left[ {\begin{array}{*{20}{c}} {\delta a}\\ {\delta b}\\ {\delta r}\\ {\delta \theta } \end{array}} \right] - \left[ {\begin{array}{*{20}{c}} {x - {x^0}}\\ {y - {y^0}} \end{array}} \right]} \end{array} $ (6)
$ {\rm{令}}\;\;\;\;\;\;\;\mathit{\boldsymbol{l}} = \left[ {\begin{array}{*{20}{c}} {x - {x^0}}\\ {y - {y^0}} \end{array}} \right],\;\;\;\;{\mathit{\boldsymbol{e}}_l} = \left[ {\begin{array}{*{20}{c}} {{v_x}}\\ {{v_y}} \end{array}} \right] $
$ \mathit{\boldsymbol{A}} = \left[ {\begin{array}{*{20}{c}} 1&0&{\cos {\theta ^0}}&{ - {r^0}\sin {\theta ^0}/\rho }\\ 0&1&{\sin {\theta ^0}}&{{r^0}\cos {\theta ^0}/\rho } \end{array}} \right],\mathit{\boldsymbol{\xi }} = \left[ {\begin{array}{*{20}{c}} {\delta a}\\ {\delta b}\\ {\delta r}\\ {\delta \theta } \end{array}} \right] $

则式(6)简化为:

$ {\mathit{\boldsymbol{e}}_l} = \mathit{\boldsymbol{A\xi }} - \mathit{\boldsymbol{l}} $ (7)

式中,l为观测向量,el为对应的观测向量误差;A为系数矩阵;ξ为参数改正数;ρ= $\left( {\frac{{180}}{\pi }} \right)$×60×60;系数矩阵 A中元素求解时需采用带有误差的观测值(x, y),故A含有随机误差,由此可以构建圆曲线拟合的EIV模型。

假设有圆曲线的一组数据为pi(xi, yi),需要由这组点拟合出圆曲线的参数,相应的误差方程可以根据式(5)表示出来。其中,

$ \begin{array}{l} \mathit{\boldsymbol{l = }}\left[ {\begin{array}{*{20}{c}} {{x_1} - x_1^0}\\ {{y_1} - y_1^0}\\ \vdots \\ {{x_i} - x_i^0}\\ {{y_i} - y_i^0} \end{array}} \right] \in {R^{2n \times 1}},\mathit{\boldsymbol{A}} = \\ \left[ {\begin{array}{*{20}{c}} 1&0&{\cos \theta _1^0}&{ - {r^0}\sin \theta _1^0/\rho }&0& \cdots &0\\ 0&1&{\sin \theta _1^0}&{ - {r^0}\cos \theta _1^0/\rho }&0& \cdots &0\\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\ 1&0&{\cos \theta _n^0}&0& \cdots &0&{ - {r^0}\sin \theta _i^0/\rho }\\ 0&1&{\sin \theta _n^0}&0& \cdots &0&{ - {r^0}\cos \theta _i^0/\rho } \end{array}} \right] \in \end{array} $

R2n×(n+3)ξ =[δa δb δc δθ1δθn]TR(n+3)×1。采用总体最小二乘解出 ξ,即确定了圆曲线方程。

2 基于Partial EIV模型的参数估计

分析圆曲线参数方程拟合模型式(7)可发现,系数矩阵 A既含有随机元素又含有重复的非随机元素。为了保证重复元素的改正数一致,且减少平差计算量,采用基于Partial EIV模型的总体最小二乘求解拟合参数。Partial EIV模型是将系数矩阵中的随机元素提取出来,所对应的数学模型如下[8]

1) 函数模型:

$\mathit{\boldsymbol{l}}\mathit{ = }\left( {{\mathit{\boldsymbol{\xi }}^{\rm{T}}} \otimes {\mathit{\boldsymbol{I}}_{2n}}} \right)\left( {\mathit{\boldsymbol{h}} + \mathit{\boldsymbol{B\bar t}}} \right) - {\mathit{\boldsymbol{e}}_l}\\ \mathit{\boldsymbol{t}} = \mathit{\boldsymbol{\bar t}} + {\mathit{\boldsymbol{e}}_\mathit{\boldsymbol{t}}} $ (8)

2) 随机模型:

$ \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{e}}_\mathit{\boldsymbol{l}}}}\\ {{\mathit{\boldsymbol{e}}_\mathit{\boldsymbol{t}}}} \end{array}} \right] \sim {N_{6n}}\left( {\left[ \begin{array}{l} 0\\ 0 \end{array} \right],\sigma _0^2\left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{Q}}_\mathit{\boldsymbol{l}}}}&0\\ 0&{{\mathit{\boldsymbol{Q}}_\mathit{\boldsymbol{t}}}} \end{array}} \right]} \right) $ (9)

式中,l 为2n×1的观测向量;h 是由系数矩阵 A按列拉直后的零元素和常数元素所组成的2n(n+3)×1的列向量;QlQt分别是观测向量 y和列向量 a的协因数阵;B是构造出的2n(n+3)×4n的已知矩阵;I2n为2n×2n的单位矩阵;σ02为单位权方差;⊗为Kronecker积;t 是由 A中不同随机元素所组成的4n×1向量,t表示 t 的真值,ett 的随机误差。圆拟合的参数方程对应的向量t =[cosθ10 sinθ10 cosθ20 sinθ20 … cosθn0 sinθn0 -r0sinθ10/ρ r0cosθ10/ρ … -r0 sinθn0/ρ r0cosθn0/ρ]T,采用Partial EIV模型进行参数计算时,所对应的已知向量 h 构造如下:

$ \mathit{\boldsymbol{h}} = {\left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{h}}_1}}&{{\mathit{\boldsymbol{h}}_2}}&{{\mathit{\boldsymbol{h}}_3}} \end{array}} \right]^{\rm{T}}} $ (10)

式中,h1n个[1 0]T组成的列向量,h2n个[0 1]T组成的列向量,h3为2n(n+1)×1的零向量。

以式(8)为基础,将第一个式子展开并移项[9]为:

$ \left. \begin{array}{l} \mathit{\boldsymbol{l}} - \left( {{\mathit{\boldsymbol{\xi }}^{\rm{T}}} \otimes {\mathit{\boldsymbol{I}}_n}} \right)\mathit{\boldsymbol{h}} = \left( {{\mathit{\boldsymbol{\xi }}^{\rm{T}}} \otimes {\mathit{\boldsymbol{I}}_n}} \right)\mathit{\boldsymbol{B\bar t}} - {\mathit{\boldsymbol{e}}_l}\\ \mathit{\boldsymbol{t}} = {\mathit{\boldsymbol{I}}_{n + 3}}\mathit{\boldsymbol{\bar t}} + {\mathit{\boldsymbol{e}}_\mathit{\boldsymbol{t}}} \end{array} \right\} $ (11)

$\mathit{\boldsymbol{L = }}\left[{\begin{array}{*{20}{c}} {\mathit{\boldsymbol{l}}-\left( {{\mathit{\boldsymbol{\xi }}^{\rm{T}}} \otimes {\mathit{\boldsymbol{I}}_{2n}}} \right)\mathit{\boldsymbol{h}}}\\ \mathit{\boldsymbol{t}} \end{array}} \right]$, $\mathit{\boldsymbol{C = }}\left[{\begin{array}{*{20}{c}} {\left( {{\mathit{\boldsymbol{\xi }}^{\rm{T}}} \otimes {\mathit{\boldsymbol{I}}_{2n}}} \right)\mathit{\boldsymbol{B}}}\\ {{\mathit{\boldsymbol{I}}_{n + 3}}} \end{array}} \right]$$\mathit{\boldsymbol{ \boldsymbol{\varDelta} = }}\left[{\begin{array}{*{20}{c}} {-{\mathit{\boldsymbol{e}}_l}}\\ {{\mathit{\boldsymbol{e}}_t}} \end{array}} \right]$$\mathit{\boldsymbol{P}} = \left[{\begin{array}{*{20}{c}} {\mathit{\boldsymbol{Q}}_l^{-1}}&0\\ 0&{\mathit{\boldsymbol{Q}}_t^{-1}} \end{array}} \right]$,则式(11)等同于:

$ \mathit{\boldsymbol{L}} = \mathit{\boldsymbol{C\bar t}} + \mathit{\boldsymbol{ \boldsymbol{\varDelta} }} $ (12)

式中,C 为新构建的系数矩阵,并把 t作为参数平差。

Δ的估值为V,并在准则VTPV =min下按间接平差原理解算平差参数 t的估值${\mathit{\boldsymbol{\hat t}}}$

$ \mathit{\boldsymbol{\hat t}} = {\left( {{\mathit{\boldsymbol{C}}^{\rm{T}}}\mathit{\boldsymbol{PC}}} \right)^{ - 1}}{\mathit{\boldsymbol{C}}^{\rm{T}}}\mathit{\boldsymbol{PL}} $ (13)

进而由${\mathit{\boldsymbol{\hat t}}}$重新构造成新的系数矩阵${\mathit{\boldsymbol{\hat A}}}$,再次按间接平差原理求解模型参数${\mathit{\boldsymbol{\hat \xi}}}$

$ \mathit{\boldsymbol{\hat \xi }} = {\left( {{{\mathit{\boldsymbol{\hat A}}}^{\rm{T}}}{\mathit{\boldsymbol{P}}_l}\mathit{\boldsymbol{\hat A}}} \right)^{ - 1}}{{\mathit{\boldsymbol{\hat A}}}^{\rm{T}}}{\mathit{\boldsymbol{P}}_l}\mathit{\boldsymbol{l}} $ (14)

最后将计算出的${\mathit{\boldsymbol{\hat \xi }}}$前3项加上初始值(${{{\hat a}^0}, {{\hat b}^0}}$)、r0得到最终的圆心坐标和半径。

基于Partial EIV模型的总体最小二乘方法求解参数方程的具体计算步骤如下。

1) 初始值:采用一般方程形式的最小二乘解作为圆心初始值(${{{\hat a}^0}, {{\hat b}^0}}$)和半径的初始值r0,由(${{{\hat a}^0}, {{\hat b}^0}}$)、r0根据式(5)、(6)、(7)得到参数方程的系数矩阵初始值 A(0)和观测向量初始值 l(0)。再用式(7)求得的${\mathit{\boldsymbol{\hat \xi }}}$作为拟合参数的初值${{\mathit{\boldsymbol{\hat \xi }}}^{\left( 0 \right)}}$,并令${{\mathit{\boldsymbol{\hat t}}}^{\left( 0 \right)}} = \mathit{\boldsymbol{t}}$

2) 由 ${{\mathit{\boldsymbol{\hat \xi }}}^{\left( i \right)}}$${{\mathit{\boldsymbol{\hat t }}}^{\left( i \right)}}$构造新的观测向量L(i)和系数矩阵C(i)

3) 将步骤2)计算的 L(i)C(i)代入式(13)计算${{\mathit{\boldsymbol{\hat t}}}^{\left( {i + 1} \right)}}$

4) 把上一步求得的${{\mathit{\boldsymbol{\hat t}}}^{\left( {i + 1} \right)}}$重新构造成改正后的系数矩阵${{\mathit{\boldsymbol{\hat A}}}^{\left( {i + 1} \right)}}$

5) 根据式(14)的迭代公式计算新的参数估值${{\mathit{\boldsymbol{\hat \xi}}}^{\left( {i + 1} \right)}}$

6) 重复步骤2)~5),直到前后两次计算的未知参数估值${{\mathit{\boldsymbol{\hat \xi}}}^{\left( {i } \right)}}$满足‖${{\mathit{\boldsymbol{\hat \xi}}}^{\left( {i + 1} \right)}}$-${{\mathit{\boldsymbol{\hat \xi}}}^{\left( {i } \right)}}$‖<ε且‖${{\mathit{\boldsymbol{\hat t}}}^{\left( {i + 1} \right)}}$-${{\mathit{\boldsymbol{\hat t}}}^{\left( {i } \right)}}$‖<δ时,迭代终止,输出所求参数的最终估值。

7) 将输出的${{\mathit{\boldsymbol{\hat \xi}}}^{\left( {i } \right)}}$前3个参数加上初始值$\left( {{{\hat a}^0}, {{\hat b}^0}} \right)$r0,得到最终的圆心坐标和半径。

一般总体最小二乘系数矩阵误差改正量为2n(n+3),而本文Partial EIV模型的总体最小二乘系数矩阵误差改正量为4n,较前者减少了2n2+2n个误差改正量的计算。特别是当拟合点个数n较多时,模型的待估量大大减少。

3 圆曲线拟合的精度评定

拟合点与圆心之间的几何距离(或称为拟合前的原半径)与拟合后的半径之差di可以作为圆拟合精度评定的一个指标[2, 11]di的计算公式为:

$ d_i^2 = {\sum\limits_{i = 1}^n {\left[ {\sqrt {{{\left( {{x_i} - a} \right)}^2} + {{\left( {{y_i} - b} \right)}^2}} - r} \right]} ^2} $ (15)

对于模拟的圆曲线,可以把已知圆曲线的方程表示成一般方程形式或式(2)的参数表达形式。通过比较模拟点解算的参数平差值与真实值分别计算的残差值,可以评价拟合效果。

4 算例 4.1 实例数据

采用文献[7]提供的一组实例数据,坐标值如表 1,分别采用以一般方程构建的高斯-马尔科夫模型的LS算法,以参数方程构建的EIV模型、Partial EIV模型的TLS算法3种算法计算,结果如表 2所示。

表 1 实例数据 Tab. 1 Measured data

表 2 3种模型下参数平差结果 Tab. 2 Comparison of parameters under three models

表 2可以看出,本文算法求解的拟合参数与EIV模型的TLS解基本一致,且两种算法的距离残差和均为1.227 599 1,这说明本文算法在圆拟合中的有效性。这主要是因为Partial EIV模型是EIV模型更为一般的表达形式,涵盖了一般总体最小二乘算法需要特殊处理的各种情况,例如部分元素为非随机量导致的 Q 奇异,A呈现结构性特征的结构总体最小二乘等。本文算法的系数矩阵改正量个数为24,较EIV模型的108个减少77.78%,可以看出本文算法大幅减少了平差计算量。

4.2 模拟数据

为进一步验证算法的可行性,采用仿真数据实验。假定拟合的圆方程为:

$ {\left( {x + 1} \right)^2} + {\left( {y + 2} \right)^2} = {3^2} $ (16)

利用Matlab在圆上取6个点,坐标如表 3所示。在xy中分别加入中误差为0.01的随机误差进行圆拟合,同样采用上面3种算法求解拟合参数。模拟100次,取参数计算的平均值,结果如表 4所示。‖(Δa, Δb, Δr)‖2为参数值与真值残差值的二范数。

表 3 未加入误差的坐标值 Tab. 3 Simulated data without putting in errors

表 4 3种模型计算100次参数平均值 Tab. 4 The three models calculate the 100 parameter mean values

表 4可以看出,以参数方程建立的拟合模型的TLS解较一般方程的LS解更接近真值,说明顾及系数矩阵误差的TLS算法拟合效果好于LS算法,而以参数方程构建的Partial EIV拟合模型总体最小二乘解与EIV模型解基本一致,且均接近真值,这说明本文解法的有效性。本文方法的拟合距离残差和为0.242 239×10-3,参数值与真值残差值的二范数为2.799 848×10-3,在3种方法中均为最小,这说明本文方法在圆曲线拟合中有一定的优势。这是因为Partial EIV模型的总体最小二乘算法能保证系数矩阵相同元素的改正数一致,常数元素的改正数为零,较EIV模型的总体最小二乘算法更为合理严谨。同算例1,本算例中系数矩阵改正量数为24个,较EIV模型的108个减少77.78%。

5 结语

针对圆曲线拟合参数求解问题,提出一种基于Partial EIV模型的总体最小二乘算法。通过模拟数据和实例数据对算法的正确性、解算精度、适用性进行了验证,得到以下结论:

1) 不同形式构建的圆曲线拟合模型均可用于求解拟合参数。以一般方程构建的最小二乘法求解简单,但拟合精度较低。

2) 以圆曲线参数形式为基础构建的拟合模型,当观测坐标有误差导致系数矩阵有误差时,需采用总体最小二乘算法求解。

3) 本文给出的基于Partial EIV模型的总体最小二乘算法适用于圆曲线拟合,为圆拟合提供了一种新方法,且在拟合精度上具有一定优势。

4) 本文方法较EIV模型的总体最小二乘法,具有系数矩阵改正量少的优点。

参考文献
[1]
Coope I D. Circle Fitting by Linear and Nonlinear Least Squares[J]. Journal of Optimization Theory and Applications, 1993, 76(2): 381-388 DOI:10.1007/BF00939613 (0)
[2]
陈明晶, 方源敏, 陈杰. 最小二乘法和迭代法圆曲线拟合[J]. 测绘科学, 2016, 41(1): 194-197 (Chen Mingjing, Fang Yuanmin, Chen Jie. Least Square Method and Iterative Method of Circular Curve Fitting[J]. Science of Surveying and Mapping, 2016, 41(1): 194-197) (0)
[3]
丁克良, 刘全利, 陈翔. 正交距离圆曲线拟合方法[J]. 测绘科学, 2008(增3): 74-75 (Ding Keliang, Liu Quanli, Chen Xiang. Orthogonal Distance Curve Fitting Method[J]. Science of Surveying and Mapping, 2008(S3): 74-75) (0)
[4]
Davis T G. Total Least Squares Spiral Circle Fitting[J]. Journal of Surveying Engineering, 1999, l25(4): 159-173 (0)
[5]
袁豹, 岳东杰, 许义, 等. 加权总体最小二乘在曲面拟合中的应用[J]. 勘察科学技术, 2013(3): 47-50 (Yuan Bao, Yue Dongjie, Xu Yi, et al. Application of Weighted Total Least Squares in Surface Fitting[J]. Journal of Surveying Science and Technology, 2013(3): 47-50 DOI:10.3969/j.issn.1001-3946.2013.03.013) (0)
[6]
鲁铁定.总体最小二乘平差理论及其在测绘数据处理中的应用[D].武汉: 武汉大学, 2010 (Lu Tieding. General Least Squares Adjustment Theory and Its Application in Surveying and Mapping Data Processing[D]. Wuhan: Wuhan University, 2010) http://www.cnki.com.cn/Article/CJFDTotal-CHXB201304028.htm (0)
[7]
Gander W, Golub G H, Strebel R. Least-Squares Fitting of Circles and Ellipses[J]. BIT Numerical Mathematics, 1994, 34(4): 558-578 DOI:10.1007/BF01934268 (0)
[8]
Xu P, Liu J, Shi C. Total Least Squares Adjustment in Partial Errors-In-Variables Models: Algorithm and Statistical Analysis[J]. Journal of Geodesy, 2012, 86(8): 661-675 DOI:10.1007/s00190-012-0552-9 (0)
[9]
王乐洋, 余航, 陈晓勇. Partial EIV模型的解法[J]. 测绘学报, 2016, 45(1): 22-29 (Wang Leyang, Yu Hang, Chen Xiaoyong. Partial EIV Model Solution[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(1): 22-29) (0)
[10]
刘经南, 曾文宪, 徐培亮. 整体最小二乘估计的研究进展[J]. 武汉大学学报:信息科学版, 2013, 38(5): 505-512 (Liu Jingnan, Zeng Wenxian, Xu Peiliang. Overview of Total Least Squares Methods[J]. Geomatics and Information Science of Wuhan University, 2013, 38(5): 505-512) (0)
[11]
同济大学数学系. 高等数学[M]. 北京: 高等教育出版社, 2010 (Department of Mathematics, Tongji University. Higher Mathematics[M]. Beijing: Higher Education Press, 2010) (0)
A Circular Curve Fitting Solution Based on Partial EIV Model
QIU Dechao1,2,3     LU Tieding1,2     DENG Xiaoyuan4     
1. Faculty of Geomatics, East China University of Technology, 418 Guanglan Road, Nanchang 330013, China;
2. Key Laboratory of Watershed Ecology and Geographical Environment Monitoring, NASMG, 418 Guanglan Road, Nanchang 330013, China;
3. Yunfu City Land and Resources Planning Administration, 11 Fuqian Road, Yunfu 527300, China;
4. Geomatics Center of Zhejiang Province, 83 North-Baochu Road, Hangzhou 310012, China
Abstract: Aiming at the problem of circular curve fitting, based on the parametric equation of circular curve, firstly this paper establishes the EIV model of circular curve fitting, and transforms the model into a more reasonable Partial EIV model for the characteristics of the coefficient matrix. Then, the formula is transformed into the least squares form, and the two-step iterative method is used to solve parameters of the model, ensuring that the corrections of the same elements in the coefficient matrix are same, and the corrections of the constant elements equals zero. Finally, combing case data shows the feasibility of this algorithm, and fitting accuracy is relatively superior.
Key words: circle fitting; parametric equation; total least squares; Partial EIV model