小区域矿区测绘中七参数坐标转换模型研究及实现.docVIP

小区域矿区测绘中七参数坐标转换模型研究及实现 　　【摘 要】目前国家矿区基础测绘目普遍采用西安80坐标系，在进行控制测量时必将涉及54坐标系与80坐标系之间的转换问题。本文在分析七参数Bursa Wolf转换模型的基础上，利用最小二乘分解法解决了矩阵求逆过程中出现的数值不稳定问题，并基于VC++平台实现了两种坐标系统间的转换问题。根据内蒙古包头哈不沁铁矿测量点，对部分GPS网控制点进行54坐标与80坐标的计算转换试验，该模型简单、方便、可行，具有一定的实用意义。 　　【关键词】坐标转换算法；七参数；GPS；最小二乘分解法 　　随着GPS技术的广泛应用，以WGS84坐标系统成果越来越多。我国法定的国家大地坐标系为北京54坐标系和西安80坐标系，北京54坐标系采用的参考椭球为克拉索夫斯基椭球体，而西安80坐标系采用的参考椭球为IAG75椭球体[1]。目前国家基础测绘及一些工程项目普遍采用西安80坐标系，在进行控制测量时必将涉及54坐标系与80坐标系之间的转换问题。宗刚军等[2]对坐标转换算法做了相关研究。本文在分析七参数Bursa Wolf转换模型的基础上，利用最小二乘分解法解决了矩阵求逆过程中出现的数值不稳定问题，并基于VC++实现了两种坐标系统间的转换问题。根据内蒙古包头哈不沁铁矿测量点，对部分GPS网控制点进行54坐标与80坐标的计算转换试验，该模型简单、方便、可行，具有一定的实用意义。 　　1 坐标转换求解模型 　　在大地测量坐标转换中，广泛使用七参数Bursa Wolf转换模型[3-4]解决54坐标系至80坐标系的转换问题，其模型为： 　　式中：X0、Y0、Z0是平移参数； 、 、 是旋转参数；X54、Y54、Z54是大地点在旧坐标系中的坐标；X80、Y80、Z80是大地点在新坐标系中的坐标。 　　上述矩阵模型的线性方程组表达式如下： 　　在具体转换时先选用最好选用4组重合点（X54，Y54，Z54）和（X80，Y80，Z80），采用最小二乘法解算出方程中的转换参数[5-7]，然后利用上述方程式实现其他数据的统一转换。 　　2最小二乘分解法 　　在用最小二乘法求解转换参数时涉及矩阵的求逆运算，以式（3）为例，其误差方程的系数阵A为N*7阶，其中N为公共点数目。利用最小二乘法时需要求ATBA的逆矩阵，也即求一个7*7阶矩阵的逆矩阵，并且对于坐标转换而言，坐标数据一般较大，因此很容易导致求逆的数值不稳定。一种合理的解决方法是采用最小二乘分解法[8]，将系数矩阵A分解为： 　　这样通过上式无需通过矩阵求逆即可得到转换参数S的最小二乘解。 　　3算法实现 　　3.1定义坐标系结构体及转换参数结构体定义坐标系结构体，用CArray地址指针存储和读取数据。 　　Typedef struct COORD struct 　　{double x；double y；double z；}COORD； 　　typedefCArray CCoordinateArray； //定义坐标数组 　　定义转换参数结构体： 　　Struct CTLSPARA 　　{doubleX0；doubleY0；doubleZ0；doubleY0；double x；double y；doublez；double k；} 　　3.2七参数坐标转换算法在已知某一矿区北京54坐标系大地测量点坐标情况下，根据坐标转换参数，可计算相应的西安80坐标系坐标数据。针对矩阵表达式（3）用最小二乘法解算出其中的转换参数，具体实现过程如下： 　　1）输入一组54坐标到坐标数组COORD54和一组80坐标到坐标数组COORD80； 　　2）把矩阵A分解成Q矩阵和R矩阵； 　　3）根据式（7）计算出坐标换算参数（七参数）：3个平移参数X0、Y0、Z0，3个旋转参数 、 、 和1个尺度参数k； 　　4）把计算后的参数代入式（3），计算相应的80坐标系的空间间直角坐标X80、Y80、Z80。在用最小二乘法求解未知的测量点数据时，利用QR矩阵分解可避免矩阵求逆过程中出现的不稳定性，使求得的测量点数据与实际数据之间误差的绝对值之和最小，从而提高54坐标向80坐标的转换精度。 　　4实例分析 　　利用内蒙古包头哈不沁铁矿的实际测量点坐标数据对程序进行验证，从10个测量点中选取3个同名点参与坐标转换参数的计算。相应的坐标数据见表1。 　　根据表1计算各坐标轴上的坐标误差XAvg、YAvg、ZAvg： 　　XAvg = dx/10=-0m） 　　YAvg= dy/10=0.0005652328（m） 　　ZAvg= dz/10=0.001953834（m） 　　把计算得到的坐标转换参数代入式（2），输入该矿区测量点的