随着测绘科技的发展,我国的测绘基准从过去长期使用参心坐标系( 主要有: 1954 年北京坐标系、1980 西安坐标系等) ,已全面进入使用地心坐标系( 2000 国家大地坐标系) 时代。无论是1954 年北京坐标系还是1980 西安坐标系都在我国长期的国民经济建设中发挥了不可替代的作用。测绘工作者利用这些基准测绘了大量成果,现在使用这些成果时需要利用一定的方法和技术将这些成果进行坐标系统变换,统一到同一个测绘基准下才能正确使用。根据多年的工作实践,就常见几种格式测绘地理信息数据利用几款常见GIS 软件进行坐标转换的方法做一说明,以期对此类工作有所借鉴。
1、我国常用的大地坐标系
1954 年北京坐标系属于参心坐标系,采用了克拉索夫斯基椭球参数,并与前苏联1942 坐标系进行联测,原点位于前苏联的普尔科沃,基准面为1956年青岛验潮站求出的黄海平均海水面。
1980 西安坐标系采用1975 年国际大地测量与地球物理联合会( IUGG) 第十六届大会推荐的地球椭球基本参数。属于参心坐标系,以JYD1968. O 系统为椭球定向基准,其原点位于陕两西安,基准面采用1985 国家高程基准。
2000 国家大地坐标系是地心坐标系,其原点为包括海洋和大气的整个地球的质量中心; 2000 国家大地坐标系的z 轴由原点指向历元2000. 0 的地球参考极的方向,该历元的指向由国际时问局给定的历元为1984. O 的初始指向推算,定向的时间演化保证相对于地壳不产生残余的全球旋转,x 轴由原点指向格林尼治参考子午线与地球赤道面( 历元2000. 0)的交点,Y 轴与Z 轴、X 轴构成右手正交坐标系。采用广义相对论意义下的尺度[1]。
2、坐标转换理论和数学模型
对坐标转换方法的研究,目前有很多理论和成果,归纳下来主要涉及两类情况: 一类是同一基准下不同的坐标形式之间的转换,如同一坐标系下空间直角坐标与大地坐标的转换; 第二类是不同基准下坐标之间的转换,是不同基准之间的转换,需要确定基准之间的转换参数,也称基准转换,如1980 坐标系到2000 国家大地坐标系的转换[2]。
2. 1 同一基准下坐标转换
2. 1. 1 大地坐标与空间直角坐标间的转换
将同一坐标参照系下的大地坐标( B、L、H) 转换为空间直角坐标( X、Y、Z) 的公式如下:
式中,N 为卯酉圈的半径,
,a 为参考椭球的长半轴,e 为参考椭球的第一偏心率。
空间直角坐标( X、Y、Z) 转换为大地坐标( B、L、H) 的公式如下:
式中,e' 为参考椭球的第二偏心率,θ = arctan