坐标系转换教材文档格式.docx
《坐标系转换教材文档格式.docx》由会员分享,可在线阅读,更多相关《坐标系转换教材文档格式.docx(13页珍藏版)》请在冰豆网上搜索。
另外,鉴于该坐标系是按局部平差逐步提供大地点成果的,不是整体平差结果,因而不可避免地出现一些矛盾和不够合理的地方。
随着我国测绘事业的发展,现在已经具备条件,可以利用我国测量资料和其他有关资料,建立起适合我国情况的新的坐标系。
1.2新1954北京坐标系
新1954年北京坐标系,是由1980年国家大地坐标系转换得来的,简称;
原1954年北京坐标系又称为旧1954年北京坐标系。
由于在全国的以GDZ80为基准的测绘成果建立之前,的测绘成果仍将存在较长时间,而与GDZ80两者之间差距较大,给成果的使用带来不便,所以又建立了作为过渡坐标系。
经过渡坐标系的转换,和的控制点的高斯平面坐标,其差值在全国80%地区内小于5m,局部地区最大达12.9m。
是在GDZ80基础上,改变GDZ80相对应的IUGG1975椭球几何参数为克拉索夫斯基椭球参数,并将坐标原点(椭球中心)平移,使坐标轴保持平行而建立起来的。
的特点是:
(1)采用克拉索夫斯基椭球参数。
(2)是综合GDZ80和建立起来的参心坐标系。
(3)采用多点定位,但椭球面与大地水准面在我国境内不是最佳拟合。
(4)定向明确,坐标轴与GDZ80相平行,椭球短轴平行于地球质心指向地极原点方向,起始大地子午面平行于我国起始天文子午面,。
(5)大地原点与GDZ80相同,但大地起算数据不同。
(6)大地高程基准采用1956年黄海高程系。
(7)与相比,所采用的椭球参数相同,其定位相近,但定向不同。
的坐标是局部平差结果,而是GDZ80整体平差结果的转换值,两者之间无全国统一的转换参数,只能进行局部转换。
1.31980年国家大地坐标系
为了适应我国大地测量发展的需要,在1978年4月于西安召开的“全国天文大地网整体平差会议”上,参加会议的专家对建立我国比1954年北京坐标系更精确的新大地坐标系进行了讨论和研究。
到会专家普遍认为1954年北京坐标系相对应的椭球参数不够精确,其椭球面与我国大地水准面差距较大,在东部经济发达地区差距高达60余米,因而建立我国新的大地坐标系是必要的。
该次会议关于建立新大地坐标系提出了如下原则:
(1)全国天文大地网整体平差要在新的坐标系的参考椭球面上进行。
为此,首先要建立一个新的大地坐标系,并命名为1980国家大地坐标系;
(2)1980国家大地坐标系的大地原点固定在我国中部,具体选址是陕西省泾阳县永乐镇;
(3)采用国际大地测量和地球物理联合会1975年推荐的四个地球椭球基本参数(),并根据这四个参数求解椭球扁率和其他参数。
(4)1980国家大地坐标系的椭球短轴平行于地球质心指向我国地极原点方向,大地起始子午面平行于格林尼治平均天文台的子午面;
(5)椭球定位参数以我国范围内高程异常值平方和等于最小为条件求解。
1980国家大地坐标系就是根据以上原则在1954年北京坐标系基础上建立起来的。
1980国家大地坐标系的特点是:
(6)采用1975年国际大地测量与地球物理联合会(IUGG)第16届大会上推荐的4个椭球基本参数。
地球椭球长半径a=6378140m,
地心引力参数GM=3.986005×
,
地球重力场二阶带球谐系数=1.08263×
地球自转角速度。
根据物理大地测量学中的有关公式,可由上述4个参数算得
地球椭球扁率a=1/298.257,
赤道的正常重力值。
(7)参心大地坐标系是在1954年北京坐标系基础上建立起来的。
(8)椭球面同似大地水准面在我国境内最为密合,是多点定位。
(9)定向明确。
椭球短轴平行于地球质心指向地极原点方向,起始大地子午面平行于我国起始天文子午面,。
(10)大地原点地处我国中部,位于西安市以北60km处的泾阳县永乐镇,简称西安原点。
(11)大地高程基准采用1956年黄海高程系。
精度高,点位中误差±
0.84m。
2.坐标的换带计算
为了限制高斯投影长度变形,将椭球面按一定经度的子午线划分成不同的投影带;
或者为了抵偿长度变形,选择某一经度的子午线作为测区的中央子午线。
由于中央子午线的经度不同,使得椭球面上统一的大地坐标系,变成了各自独立的平面直角坐标系,就需要将一个投影带的平面直角坐标系,换算成另外一个投影带的平面直角坐标,称为坐标换带。
2.1坐标换带的方法
坐标换带有直接换带计算法和间接换带计算法两种。
目前采用间接换带计算法,因此下面仅就此方法作一介绍。
如将第一带(东带或西带)的平面坐标换算为第二带(西带或东带)的平面坐标,方法是先根据第一带的平面坐标x,y和中央子午线的经度L。
按高斯投影坐标反算公式求得大地坐标B,L然后根据B,L和第二带的中央子午线经度按高斯投影坐标正算公式求得在第二带中的平面坐标,。
由于在换带计算中,把椭球面上的大地坐标作为过渡坐标,因而称为间接换带法。
这种方法理论上是严密的,精度高,而且通用性强,他适用于6°
带与6°
带,3°
带与3°
带,6°
带之间的坐标换带。
如图所示,图(a)是P点在Ⅰ带平面直角坐标系的投影,它的平面直角坐标是;
图(b)是该P点在Ⅱ带平面直角坐标系的投影,它的平面直角坐标是。
假若已知要求,或已知要求,这就是所谓的邻带坐标换算。
利用高斯投影正反算公式进行邻带坐标换算的实质是把椭球面上的大地坐标作为过渡坐标。
其解法是,首先利用高斯投影坐标反算公式,根据换算成椭球面大地坐标,进而得到。
然后再由大地坐标,利用高斯投影坐标正算公式,根据计算该点在Ⅱ度带的平面直角坐标,但在这一步计算时,要根据第Ⅱ带的中央子午线的经度计算P点在第Ⅱ带的经差。
为了检核计算的正确性,每步都需要往返计算。
2.2坐标换带的实际应用
在生产实践中通常有以下两种情况需要换带计算
(1)控制网中的已知点位于相邻的两个投影带中。
如图一
(图一:
坐标换带示意图)
中的附合导线,A,B,C,D为已知高级点。
A,B两点位于西带内,具有西带的高斯平面直角坐标值;
C,D两点位于东带内,具有东带的高斯平面直角坐标值。
在坐标平差计算时,就必须将它们的坐标系统统一起来,或是将A,B点的西带坐标值换算至东带,或是将C,D点的东带坐标值换算至西带。
(2)国家控制点的坐标通常是6°
带的坐标,而在工程测量中往往需要采用带或1.5°
带,这就产生了6°
带与带或1.5°
带之间的坐标换算问题。
我们知道,带的中央子午线中,有半数与6°
带的中央子午线重合。
所以,由6°
带到3°
带的换算区分为2种情况:
①3°
带的中央子午线重合如图二所示,3°
带第
(图二:
41带与6°
第21带的中央子午线重合。
既然中央子午线一致,坐标系统也就一致。
所以,图中P1点在6°
带第21带的坐标,也就是该点在3°
带第41带的坐标。
在这种情况下,6°
带之间,不存在换带计算问题。
②3°
带中央子午线与6°
带分带子午线不重合如图所示,若已知P2点在6°
带第21带的坐标,求它在3°
带第42带的坐标。
由于这2个投影带的中央子午线不同,坐标系统不一致,必须进行换带计算。
不过P2点在6°
带第21带的坐标与它在3°
第41带的坐标相同,所以6°
带坐标换算,也可看作是3°
带的邻带坐标换算。
换带计算目前广泛采用高斯投影坐标正反算方法,他适用于任何情况下的换带计算工作。
这种方法的程序是:
首先将某投影带的已知平面坐标,按高斯投影坐标反算公式求得其大地坐标(B,L);
然后根据纬度B和对于所选定的中央子午线的经差,按高斯投影坐标正算公式求其在选定的投影带的平面坐标。
例如,某点A在新54坐标系6°
带的平面坐标为
=3589644.287=20679136.439
求A点在3°
带的平面直角坐标(x2,y2).
首先确定A点所在投影带中央子午线经度。
由横坐标的规定值可以直接判定,A点位于6°
带第20带,其中央子午线经度L。
=117°
,横坐标的自然值为y1=679136.493-500000=+179136.439m;
该坐标等同于3°
带第39带的平面坐标。
其次将已知的6°
带坐标反算为大地坐标。
为此,可直接应用坐标反算公式进行计算,其结果为
B=32°
2457.6522"L=118°
5415.2206"
由大地经度L可判断,A点位于3°
第40带,中央子午线为L。
=120°
。
最后根据高斯投影坐标正算公式,由已知的纬度B和经度计算A点在3°
带第40带的平面直角坐标,得
=3588576.591=40396922.874
其中横坐标为规定值。
3.坐标系转换
首先,我们要弄清楚几种坐标表示方法。
大致有三种坐标表示方法:
经纬度和高程,空间直角坐标,平面坐标和高程。
我们通常说的西安80坐标是经纬度和高程这一种,北京54坐标是平面坐标和高程着一种。
现在,再搞清楚转换的严密性问题,在同一个椭球里的转换都是严密的,而在不同的椭球之间的转换这时不严密的。
举个例子,在西安坐标和北京54坐标之间是不存在一套转换参数可以全国通用的,在每个地方会不一样,因为它们是两个不同的椭球基准。
那么,两个椭球间的坐标转换应该是怎样的呢?
一般而言比较严密的是用七参数法,即X平移,Y平移,Z平移,X旋转,Y旋转,Z旋转,尺度变化K。
要求得七参数就需要在一个地区需要3个以上的已知点,如果区域范围不大,最远点间的距离不大于30Km(经验值),这可以用三参数,即X平移,Y平移,Z平移,而将X旋转,Y旋转,Z旋转,尺度变化K视为0,所以三参数只是七参数的一种特例。
在本软件中提供了计算三参数、七参数的功能。
在一个椭球的不同坐标系中转换需要用到四参数转换,举个例子,在深圳既有北京54坐标又有深圳坐标,在这两种坐标之间转换就用到四参数,计算四参数需要两个已知点。
下面举例子说明:
在湖北有一个测区,需要完成北京54坐标系到西安80坐标系的转换,整个转换过程是:
3.1四参数转换
如果区域范围不大,可采用4参数方法计算坐标转换,计算4参数法需要两个已知点的80和54两套坐标。
该方法计算简单且在一定范围内计算结果精度较高,有一定的实用性。
四参数即两个平移参数,一个旋转参数,一个尺度参数,四参数没有转换模型,是强制转换。
实例:
湖北用户在一个测区内有一些点的北京54的坐标,现在希望将其转换为国家80坐标。
用户有测区的七个GPSC级点(112.5度带)控制点,这些控制点既有北京54坐标也有国家80坐标,转换步骤为:
3.1.24参数计算