基于VBA平台七参数法坐标系统转换程序设计.docVIP

基于VBA平台七参数法坐标系统转换程序设计 　　【摘 要】 本文研究坐标系之间的坐标转换方法，特别是空间直角坐标系之间的坐标转换。通过详细研究1954北京坐标系和1980国家大地坐标系之间的相互转换方法，例如欧勒角、三参数法、七参数法及多项式回归模型等，通过Excel2007 VBA平台来实现坐标系转换过程中的参数求取、精度估算，坐标转换等工作。 　　【关键词】 1954北京坐标系 1980国家大地坐标系 坐标系间转换 VBA 　　1 坐标转换概述 　　坐标转换是空间实体的位置描述，是从一种坐标系统变换到另一种坐标系统的过程。通过建立两个坐标系统之间一一对应关系来实现。是各种比例尺地图测量和编绘中建立地图数学基础必不可少的步骤。在测量中，有两种意义的坐标转换，一是地图投影变换，即从一种地图投影转换到另一种地图投影，地图上各点坐标均发生变化；另一是量测系统坐标转换，即从大地坐标系??地图坐标系、数字化仪坐标系、绘图仪坐标系或显示器坐标系之间的坐标转换。 　　2 坐标系统转换方法概述 　　2.1 不同空间直角坐标系间的转换 　　2.1.1 欧勒角 　　高等数学解析几何中，关于空间直角坐标系的转换，包括坐标轴的平移和坐标轴的旋转，坐标轴的三个旋转角叫欧勒角。对于不同的空间直角坐标系与，设它们的原点一致，相应的坐标轴互不平行，其欧勒角分别为。按以下步骤，将转换为。 　　2.1.2 三参数法 　　设两个空间直角坐标系分别为O1-X1Y1Z1与O2-X2Y2Z2，它们的原点不一致，相应的坐标轴相互平行，则有： 　　三参数坐标转换公式是在假设两坐标系间各坐标轴相互平行，轴系间不存在欧勒角的条件下得出的，实际应用中，因为欧勒角不大，可以用三参数公式近似地进行空间直角坐标系统的转换。 　　2.1.3 七参数法 　　进行两个不同空间直角坐标系统之间的坐标转换，需要求出坐标系统之间的转换参数。转换参数一般是利用重合点的两套坐标值通过一定的数学模型进行计算。当重合点数为三个以上时，可以采用七参数法进行转换[1]。 　　七参数公式比较三参数公式能获得较高精度的转换结果。实际应用中，也可以舍弃不显著的参数，例如个别欧勒角，选择四、五或六个参数进行不同空间直角坐标系统的转换。 　　转换参数可以通过联测一些公共点获得。 　　2.1.4 坐标转换多项式回归模型 　　坐标转换七参数公式模型属于相似变换模型。大地控制网中的系统误差一般呈区域性，当区域较小时，区域性的系统误差被相似变换参数拟合，故局部区域的坐标转换采用七参数公式模型是比较适合的。但对全国或是一个省区范围内的坐标转换，可以采用多项式回归模型，将各区域的系统偏差拟合到回归参数中，从而提高坐标转换精度。 　　2.2 不同平面坐标系的转换 　　对于局部地区不同的平面坐标系统的转换采用平面坐标系统相似变换模型会更容易些。研究表明，在局部范围内，例如在每十万分之一图幅内，同一点在1980年国家坐标系和原1954年北京坐标系的高斯平面直角坐标，在米级精度上都只相差一个常数。对于米级以下的微小差异，可以看成是这个局部区域内的两个平面坐标系统之间存在着某种旋转和尺度伸缩造成的。这样，就可以用平面相似变换公式来模拟两个平面坐标系统的转换关系。 　　3 基于VBA for Microsoft Excel2007平台坐标系统转换程序设计 　　3.1 VBA编程简介 　　Visual Basic for Applications（简称VBA）是新一代标准宏语言，是基于Visual Basic for Windows 发展而来的。它与传统的宏语言不同，传统的宏语言不具有高级语言的特征，没有面向对象的程序设计概念和方法。而VBA 提供了面向对象的程序设计方法，提供了相当完整的程序设计语言。VBA 易于学习掌握，可以使用宏记录器记录用户的各种操作并将其转换为VBA 程序代码。这样用户可以容易地将日常工作转换为VBA 程序代码，使工作自动化。 　　3.2 七参数法程序设计 　　3.2.1 七参数法转换模型 　　定义I点在54和80坐标系的坐标分别为和，则根据七参数布尔莎转换模型，有如下关系： 　　定义系数矩阵Bi和权因数阵P，由下式确定： 　　组成常数项矩阵，由下式确定： 　　求取法方程矩阵，如下： 　　转换参数矩阵： 　　精度评定部分： 　　改正数矩阵， 　　单位权中误差， 　　方差协方差矩阵为：，其主对角线上元素即为各参数的精度。 　　坐标转换部分： 　　将求得的七参数带入如下矩阵，即可求出转换后坐标[2]： 　　3.2.2 七参数法转换程序部分 　　运用VBA进行七参数程序转化，七参数求取及坐标转换程序主界面如图1。 　　程序共含两个Excel文件和