坐标拟合的双向解算与矩阵系数的生成.pdf

第5卷第1集 V01．5，No．1 2005年12月 现代工程测量技术发展与应用研讨交流会论文集 Dec．2005 坐标拟合的双向解算与矩阵系数的生成 周 卫 (南京师范大学地理科学学院江苏南京210097) 摘要：在工程测绘等领域，四参数与七参数拟合是坐标转换的重要方法。采取重合点坐标可以获得正 反算两套拟合参数，从而实现正反算求解。本文阐述了基于四参数与七参数这两种转换模型，通过重组矩阵 模型，推求反算数学模型，采用一套拟合参数进行正反算的方法，并给出了相应系数矩阵的自动生成程序和 解算过程。 关键词：四参数、七参数、拟合模型、正算、反算 1 引言 四参数与七参数拟合坐标转换是实现坐标转换的重要方法。在工程测量等领域有着广泛的 应用。目前，GPSRTK技术也是采用参数拟合的方法进行实时的坐标转换。采取重合点坐标 可以获得正反算两套拟合参数，从而实现正反算求解。本文阐述了基于这两种转换模型，采 用一套拟合参数实现坐标正反算的求解方法。推求出四参数与七参数反算数学模型，给出了 相应参数矩阵的生成程序。根据本文阐述的方法和提供的有关模型可以开发出功能实用的坐 标换算程序或坐标换算的动态连接库，提供工程测绘等坐标转换模块的集成与开发。 2拟合模型表达与解算 2．1四参数拟合模型的矩阵表达(赫尔默特法)： 阱卧叶m僻=捌 (1) 将上述矩阵展开，其线性方程组如下： X0+Xl*k*cosA—y1*k*sinA=x2 i (2) yo+Xl*k*snA十yl水k木cosA=y2 (上式中，k=l+m) 令：kc=k*cosA，ks=k*sinA，则有线性方程组： l*xo+O*yo+xl木kc+(一y1)*ks=x2 0爿cxo+1：lcyo+yL木kc+xt*ks=yz (3) ．182· 现代工程测量技术发展与应用研讨交流会论文集 第5卷 上述线性方程组的矩阵表达如下： AS=B (4) 其中： ％ o yo A：『1 x-咄]，s： Lo 1 Y1 zlJ kc廿圈 ks 以上述方程系数组成矩阵解算，获得Xo、Y。、kc、ks。将该四参数代入(3)式即获得四参 数的正算方程。 。 如有必要，可继续解算A和K。由于sin2A+COS2A=I，令ksc=ks／kc，则有： A=arcsin(ksc／sqr(1+ksc2)) (5) 由此继续解出k，注意在旋转角很小情况下应采用下式解算k： ’ k=ks／sinA (6) 上述式中：xly。一源(旧)坐标，x2yz一目标(新)坐标；Xo、yo一坐标平移，A一旋转角， k一尺度比(尺度系数m=k-1)。 2．2七参数拟合模型的矩阵表达(布尔沙一沃尔夫Bursa．Wolf模型)： X2 Xo 0 —Zl ‘