留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

坐标参数化的平面基准转换方法

王利朋 刘成龙 刘胜 何林烜

王利朋, 刘成龙, 刘胜, 何林烜. 坐标参数化的平面基准转换方法[J]. 武汉大学学报 ● 信息科学版, 2016, 41(10): 1409-1413. doi: 10.13203/j.whugis20140416
引用本文: 王利朋, 刘成龙, 刘胜, 何林烜. 坐标参数化的平面基准转换方法[J]. 武汉大学学报 ● 信息科学版, 2016, 41(10): 1409-1413. doi: 10.13203/j.whugis20140416
WANG Lipeng, LIU Chenglong, LIU Sheng, HE Linxuan. A Method of Plane Datum Transformation by Coordinate Parameterization[J]. Geomatics and Information Science of Wuhan University, 2016, 41(10): 1409-1413. doi: 10.13203/j.whugis20140416
Citation: WANG Lipeng, LIU Chenglong, LIU Sheng, HE Linxuan. A Method of Plane Datum Transformation by Coordinate Parameterization[J]. Geomatics and Information Science of Wuhan University, 2016, 41(10): 1409-1413. doi: 10.13203/j.whugis20140416

坐标参数化的平面基准转换方法

doi: 10.13203/j.whugis20140416
基金项目: 

长江学者和创新团队发展计划 No.IRT13092

详细信息
    作者简介:

    王利朋,硕士生,主要从事精密工程测量与变形监测研究。Wanglpswjtu@sina.com

  • 中图分类号: P207.2

A Method of Plane Datum Transformation by Coordinate Parameterization

Funds: 

The Program for Changjiang Scholars and Innovative Research Teams in Universities No.IRT13092

计量
  • 文章访问数:  1054
  • HTML全文浏览量:  43
  • PDF下载量:  383
  • 被引次数: 0
出版历程
  • 收稿日期:  2014-09-18
  • 刊出日期:  2016-10-05

坐标参数化的平面基准转换方法

doi: 10.13203/j.whugis20140416
    基金项目:

    长江学者和创新团队发展计划 No.IRT13092

    作者简介:

    王利朋,硕士生,主要从事精密工程测量与变形监测研究。Wanglpswjtu@sina.com

  • 中图分类号: P207.2

摘要: 工程测量中经常需要实现不同坐标系下成果的相互转换,而高精度的转换参数是完成这一工作的基础。获取基准转换参数的实质就是利用公共点在两套坐标系下的坐标,根据相似变换原理建立误差方程求解。传统的最小二乘(LS)相似变换法只考虑了公共点在一套坐标系下的误差,与实际情况不符。基于此,探讨了坐标参数化的平面基准转换方法,解决了考虑公共点在两套坐标系下坐标都含有误差时高斯-马尔科夫(Gauss-Markov,G-M)模型不成立的问题,以相似变换原理为基础,采取通用的最小二乘方法解算基准转换参数。

English Abstract

王利朋, 刘成龙, 刘胜, 何林烜. 坐标参数化的平面基准转换方法[J]. 武汉大学学报 ● 信息科学版, 2016, 41(10): 1409-1413. doi: 10.13203/j.whugis20140416
引用本文: 王利朋, 刘成龙, 刘胜, 何林烜. 坐标参数化的平面基准转换方法[J]. 武汉大学学报 ● 信息科学版, 2016, 41(10): 1409-1413. doi: 10.13203/j.whugis20140416
WANG Lipeng, LIU Chenglong, LIU Sheng, HE Linxuan. A Method of Plane Datum Transformation by Coordinate Parameterization[J]. Geomatics and Information Science of Wuhan University, 2016, 41(10): 1409-1413. doi: 10.13203/j.whugis20140416
Citation: WANG Lipeng, LIU Chenglong, LIU Sheng, HE Linxuan. A Method of Plane Datum Transformation by Coordinate Parameterization[J]. Geomatics and Information Science of Wuhan University, 2016, 41(10): 1409-1413. doi: 10.13203/j.whugis20140416
  • 二维基准转换在实际工程中使用频繁,不同平面坐标系统之间的转换就显得非常重要,而常用的二维基准转换方法为相似变换或正形投影变换[1]。为了求解转换参数,必须已知足够的公共点,建立高斯-马尔科夫(Gauss-Markov,G-M)平差模型[2]。实际上采用的公共点坐标本质上也是通过测量手段获取的,公共点在不同坐标系下的坐标均含有误差。若公共点两套坐标都含有误差,则平面坐标转换中误差方程的系数矩阵B不再是完全准确的,这种情况下经典的G-M模型将不再适用[3]

    为了解决坐标转换中公共点在两套坐标系下的坐标均含有误差这一问题,近年来,众多研究提出了各种不同的坐标转换模型或理论,无论是在二维平面还是三维空间坐标转换中都具有良好的改善效果。例如无缝三维基准转换模型l[4]、非线性坐标转换的稳健估计[5]、整体最小二乘法[6]、加权总体最小二乘法[7]、改进的加权总体最小二乘法[8]、约束总体最小二乘法以及约束加权混合整体最小二乘法[3, 9]等。以上所有模型都是考虑系数矩阵B含有误差eB,基于目标函数elTQl-1el+eBTQB-1eB=min,或同时考虑限制条件和定权方法而建立的,其模型复杂,求解过程较为繁琐。为了高效、简便地获取基准转换参数,本文在二维平面坐标系统下,基于相似变换原理,探讨了坐标参数化的最小二乘(coordinates parameterization least squares,CPLS)平面基准四参数转换方法。基本思想是:考虑公共点在两套平面坐标系下的坐标均含有误差,将公共点在 Ⅰ 套坐标系下的坐标(X,Y)和转换参数作为未知参数,按照目标函数VTPV=min求解转换参数并评定精度,V为公共点在Ⅱ套坐标系中的坐标改正数。

    • 平面基准转换的相似变换原理为:

      (1)

      其矩阵形式为:

      (2)

      式中,公共点i=1,2,…,n;(Xi,Yi)、(Xi,Yi)分别为第i个公共点在Ⅱ套和Ⅰ套坐标系下的坐标;[X0Y0]T为平移参数;θ为旋转角;m为尺度因子。

      Z1=mcosθ、Z2=msinθ,以[X0 Y0 Z1 Z2]T为参数,按照G-M模型建立误差方程:

      (3)

      式中,[δX0δY0δZ1δZ2]iT为参数改正数;(X0,Yi0)为第i个公共点在Ⅱ套坐标系下的近似坐标,可以通过参数近似值和公共点的Ⅰ套坐标推算得出;[vX1vY1…vXnvYn]T为各点在Ⅱ套坐标系下坐标的改正数。

      传统的LS平面基准转换不考虑Ⅰ套坐标系下[XiYi]T的误差,只考虑Ⅱ套坐标系下[XiYi]T的误差,按照VTPV=min的准则进行平差计算,其中V=[VXVY]T。不考虑公共点在Ⅰ套坐标系下坐标的误差,显然不符合实际情况。为此衍生出了TLS等多种改进的求解模型和方法[2-9]。文献[2-9]中所提方法均以系数矩阵B含有误差eB为基础,对整体最小二乘坐标转换进行推导和分析。此类方法求解过程复杂,所需模型较多而且步骤繁琐。因此,本文探讨了更加简便的坐标参数化的平面基准转换方法,该方法的核心是将公共点在Ⅰ套坐标系下的坐标作为未知参数进行平差,从而解决公共点两套坐标均含有误差这一问题,使得转换参数更加准确。

    • 对式(2) 进行变形得:

      (4)

      将尺度参数m与Ⅰ套坐标(X,Y)相乘构造缩放后的Ⅰ套坐标(a,b)。本节内容根据公共点的Ⅰ套缩放坐标和Ⅱ套坐标求解平移参数和旋转角,尺度参数m在§1.3节中解算。

      ai=mXi,bi=mYi,c=cosθ,d=sinθ,则式(1) 可转化为:

      (5)

      Xi、X0、Y0、ai、bi、c、d的近似值分别为Xi0、X00、Y00、ai0、bi0、c0、d0,相应的改正数为vXi、δX0、δYo、δai、δbi、δc、δd,忽略微小量δai、δbi、δc、δd之间的乘积,则式(5) 可转化为:

      (6)

      式中,δai、δbi为第i个公共点Ⅰ套缩放坐标的改正数。

      考虑到公共点在两套坐标系下的坐标都含有误差,则平面基准转换的观测方程为:

      (7)

      vai、vbi为第i个公共点Ⅰ套缩放坐标的残差,vai=ai-ai0,vbi=bi-bi0,将vai、vbi和式(6) 代入到式(7) 中得误差方程为:

      (8)

      方程中观测值个数为4n,未知数个数为2n+4,误差方程的矩阵形式为:

      (9)

      式中,I为2n阶单位阵;A为2n阶对角阵;B为2n×4阶矩阵。

      通过以上步骤,考虑公共点在两套坐标系下的坐标误差,建立了G-M平差模型,按照间接平差法进行平差计算,假设观测值的权为单位权,由式(9) 得最小二乘解为:

      (10)
    • 根据式(10) 可得:

      (11)

      由于公共点Ⅰ套缩放坐标ai=mXi,bi=mYi,为了求解尺度参数m,由全微分公式得:

      (12)

      式中,m0为尺度近似值;δm为尺度改正数;(Xi0,Yi0)为公共点在Ⅰ套坐标系下的坐标。实例证明了在尺度接近1而且坐标转换范围为小区域的情况下,Ⅰ套坐标系下各点坐标改正数的微小差异对尺度的精度影响可以忽略不计。因此,假设各公共点在Ⅰ套坐标系下的坐标改正数相同,采用虚拟的坐标改正数(δX,δY)与尺度改正数δm构成误差方程(13) 的解向量。

      式(12) 的矩阵形式为:

      (13)

      式(13) 为超定方程组,采用超定方程组的最小二乘解法可得:

      (14)

      从而得到尺度参数:

      (15)

      将式(10) 代入式(9) 可以得到CPLS平面基准转换的残差:

      (16)

      按照目标函数F=VTV=min,可得验后单位权中误差:

      (17)
    • 本文按照平移参数X0=5.000 0 m、Y0=8.000 0 m、尺度参数m=1.0和旋转角θ=50.0°构造了6个公共点在两套坐标系下的坐标(${\hat{X}}$,${\hat{Y}}$)和(${\hat{X}}$,${\hat{Y}}$),然后对所有的坐标加0~2 mm的随机误差,采用含有误差的坐标作为算例分析的样本数据,公共点在两套坐标系下的坐标如表 1所示。对表 1中含扰动误差的数据分别采用LS和CPLS两种方法解算参数并进行精度评定,对比结果如表 2所示。为了比较不同方法求解参数的可靠性,分别采用两种方法解算的转换参数,根据公共点在Ⅰ套坐标系下的坐标计算公共点在Ⅱ套坐标系下的坐标(XLS,YLS)和(XCPLS,YCPLS),并与X、Y比较计算坐标差ΔXLS、ΔYLS、ΔXCPLS和ΔYCPLS,如表 3所示。将ΔXCPLS、ΔYCPLS的绝对值与ΔXLS、ΔYLS的绝对值坐标差得ΔXCPLS-LS、ΔYCPLS-LS,如表 4所示。

      表 1  平面基准转换样本点/m

      Table 1.  Sample Points of Plane Datum Transformation/m

      点号 无误差 加扰动误差
      ${\hat{X}}$Ⅰ ${\hat{Y}}$Ⅰ ${\hat{X}}$Ⅱ ${\hat{Y}}$Ⅱ X Y X Y
      1 1.000 00 8.000 00 11.771 14 12.376 26 1.000 88 8.000 27 11.772 68 12.376 43
      2 1.000 00 9.000 00 12.537 19 13.019 04 1.001 77 9.001 94 12.538 21 13.020 22
      3 1.000 00 10.000 00 13.303 23 13.661 83 1.001 06 10.001 81 13.304 65 13.663 61
      4 3.000 00 8.000 00 13.056 72 10.844 17 3.001 42 8.000 09 13.058 37 10.845 02
      5 3.000 00 9.000 00 13.822 76 11.486 96 3.000 51 9.000 20 13.823 73 11.487 82
      6 3.000 00 10.000 00 14.588 81 12.129 74 3.001 99 10.000 83 14.590 30 12.130 81

      表 2  平面基准转换不同方法解算的精度和参数估值

      Table 2.  Estimation of Precision and Parameters from Different Plane Datum Transformation Methods

      方法σ0 /mmX0 /mY0 /mθ/ (°)尺度参数
      LS 0.60 4.999 68 8.000 45 49.995 938 1.000 054 99
      CPLS 0.42 5.002 47 7.996 23 49.996 494 1.000 000 00

      表 3  平面基准转换不同方法计算的Ⅱ套坐标及坐标差

      Table 3.  Coordinates and Difference in Coordinate System Ⅱ from Different Plane Datum Transformation Methods

      点号 XLS/m YLS/m XCPLS/m YCPLS/m ΔXLS/mm ΔYLS/mm ΔXCPLS/mm ΔYCPLS/mm
      1 11.771 66 12.376 92 11.771 66 12.376 15 1.02 -0.49 1.02 0.28
      2 12.539 55 13.020 20 12.539 16 13.019 80 -1.34 0.02 -0.95 0.42
      3 13.305 04 13.663 53 13.304 30 13.663 51 -0.39 0.08 0.35 0.09
      4 13.057 61 10.844 31 13.058 38 10.844 32 0.75 0.71 -0.01 0.70
      5 13.823 16 11.487 96 13.823 53 11.488 35 0.58 -0.14 0.20 -0.53
      6 14.590 63 12.130 11 14.590 65 12.130 87 -0.33 0.70 -0.35 -0.07

      表 4  平面基准转换不同方法计算的Ⅱ套坐标差之差/mm

      Table 4.  Difference of Coordinate Difference from Different Plane Datum Transformation Methods/mm

      点号 |ΔXLS| |ΔYLS| |ΔXCPLS| |ΔYCPLS| ΔXCPLS-LS ΔYCPLS-LS
      1 1.02 0.49 1.02 0.28 0.00 -0.21
      2 1.34 0.02 0.95 0.42 -0.39 0.39
      3 0.39 0.08 0.35 0.09 -0.04 0.02
      4 0.75 0.71 0.01 0.70 -0.74 -0.01
      5 0.58 0.14 0.20 0.53 -0.37 0.39
      6 0.33 0.70 0.35 0.07 0.02 -0.63

      表 2中可以看出,两种方法的验后单位权中误差分别为σ0LS=0.60 mm、σ0CPLS=0.42 mm。这说明基于坐标参数化的平面基准转换精度高于普通的平面基准转换方法。结合LS、CPLS解算的参数估值与表 3中两种方法计算的坐标差可以看出,CPLS解算的坐标较差总体上更小,证明该方法解算的参数更加接近真值。为了比较两种方法的差异,采取了不同方法计算的Ⅱ套坐标差之差这一指标,如表 4所示,取绝对值来表示反算的坐标与样本点坐标的偏离程度,坐标较差的绝对值越大,偏离程度越大,反之就越小。用ΔXCPLS-LS、ΔYCPLS-LS可以比较出两种方法反算坐标偏离程度的大小,根据大量仿真实例计算的统计结果得出,CPLS反算坐标偏离程度小于LS反算坐标偏离程度的比例为70%左右,即CPLS有效地提高了平面基准转换的精度,表 2~表 4中的精度指标和坐标差等数据说明CPLS较LS可以求得更高精度的平面基准转换参数。

    • 本文从理论推导和实例计算上验证了坐标参数化的平面基准转换方法要优于普通的平面基准转换方法,因为基于CPLS的方法同时考虑了公共点在两套坐标系下的坐标误差,在选择合适公共点的前提下,CPLS获得的基准转换参数更符合实际情况。

      仿真实例的统计结果表明,CPLS验后单位权中误差小于LS验后单位权中误差的比例为100%,CPLS反算公共点在Ⅱ套坐标系下的坐标偏离程度小于LS反算坐标偏离程度的比例约为70%。CPLS验后单位权方差较LS有所改进,前者反算公共点在Ⅱ套坐标系下的坐标偏离程度较后者有所降低。

参考文献 (12)

目录

    /

    返回文章
    返回