一、常用坐标系1、北京坐标系北京54坐标系为参心大地坐标系,大地上的一点可用经度L54、纬度M54和大地高H54定位,它是以克拉索夫斯基椭球为基础,经局部平差后产生的坐标系。
1954年北京坐标系的历史:新中国成立以后,我国大地测量进入了全面发展时期,再全国范围内开展了正规的,全面的大地测量和测图工作,迫切需要建立一个参心大地坐标系。
由于当时的“一边倒”政治趋向,故我国采用了前苏联的克拉索夫斯基椭球参数,并与前苏联1942年坐标系进行联测,通过计算建立了我国大地坐标系,定名为1954年北京坐标系。
因此,1954年北京坐标系可以认为是前苏联1942年坐标系的延伸。
它的原点不在北京而是在前苏联的普尔科沃。
北京54坐标系,属三心坐标系,长轴6378245m,短轴6356863,扁率1/298.3;2、西安80坐标系1978年4月在西安召开全国天文大地网平差会议,确定重新定位,建立我国新的坐标系。
为此有了1980年国家大地坐标系。
1980年国家大地坐标系采用地球椭球基本参数为1975年国际大地测量与地球物理联合会第十六届大会推荐的数据,即IAG75地球椭球体。
该坐标系的大地原点设在我国中部的陕西省泾阳县永乐镇,位于西安市西北方向约60公里,故称1980年西安坐标系,又简称西安大地原点。
基准面采用青岛大港验潮站1952-1979年确定的黄海平均海水面(即1985国家高程基准)。
西安80坐标系,属三心坐标系,长轴6378140m,短轴6356755,扁率1/298.257221013、2000国家大地坐标系的定义国家大地坐标系的定义包括坐标系的原点、三个坐标轴的指向、尺度以及地球椭球的4个基本参数的定义。
2000国家大地坐标系的原点为包括海洋和大气的整个地球的质量中心;2000国家大地坐标系的Z轴由原点指向历元2000.0的地球参考极的方向,该历元的指向由国际时间局给定的历元为1984.0的初始指向推算,定向的时间演化保证相对于地壳不产生残余的全球旋转,X轴由原点指向格林尼治参考子午线与地球赤道面(历元2000.0)的交点,Y轴与Z轴、X轴构成右手正交坐标系。
采用广义相对论意义下的尺度。
2000国家大地坐标系,长半轴6378137m,扁率f=1/298.257222101,地心引力常数GM =3.986004418×1014m3s-2,自转角速度ω=7.292l15×10-5rads-1。
4、1984世界大地坐标系(WGS84坐标系WorldGeodeticSystem)wgs-84坐标系是美国国防部研制确定的大地坐标系,是一种协议地球坐标系。
wgs-84坐标系的定义是:原点是地球的质心,空间直角坐标系的z轴指向bih(1984.0)定义的地极(ctp)方向,即国际协议原点cio,它由iau和iugg共同推荐。
x轴指向bih定义的零度子午面和ctp 赤道的交点,y轴和z,x轴构成右手坐标系。
wgs-84椭球采用国际大地测量与地球物理联合会第17届大会测量常数推荐值,采用的两个常用基本几何参数:长半轴a=6378137m;扁率f=1:298.257223563。
GPS广播星历是以WGS-84坐标系为根据的。
5 高斯平面直角坐标系和utm一般的地图均为平面图,其对应的也是平面坐标.因此,需要将椭球面上各点的大地坐标,按照一定的数学规律投影到平面上成为平面直角坐标.目前世界各国采用最广泛的高斯- 克吕格投影和墨卡托投影(utm)均是正形投影(等角投影),即该投影在小区域范围内使平面图形与椭球面上的图形保持相似。
为了限制长度变形,,根据国际测量协会规定,将全球按一定经差分成若干带。
我国采用6度带或3度带,6度带是自零度子午线起每隔经度。
高斯平面直角坐标系一般以中央经线(l0)投影为纵轴x, 赤道投影为横轴y,两轴交点即为各带的坐标原点。
为了避免横坐标出现负值,在投影中规定将坐标纵轴西移500公里当作起始轴。
为了区别某一坐标系统属于哪一带,通常在横轴坐标前加上带号,如(4231898m,21655933m),其中21即为带号。
城建坐标多采用三度带的高斯-克吕格投影。
同一坐标系下的大地坐标(即经纬度坐标b,l)与其对应的高斯平面直角坐标(x,y)有严格的转换关系。
现行的测绘的教科书的一般都有。
6、地方独立坐标系在我国许多城市测量与工程测量中,若直接采用国家坐标系下的高斯平面直角坐标,则可能会由于远离中央子午线,或由于测区平均高程较大,而导致长度投影变形较大,难以满足工程上或实用上的精度要求。
另一方面,对于一些特殊的测量,如大桥施工测量,水利水坝测量,滑坡变形监测等,采用国家坐标系在实用中也会很不方便。
因此,基于限制变形,以及方便实用,科学的目的,在许多城市和工程测量中,常常会建立适合本地区的地方独立坐标系。
建立地方独立坐标系,实际上就是通过一些元素的确定来决定地方参考椭球与投影面.地方参考椭球一般选择与当地平均高程相对应的参考椭球,该椭球的中心,轴向和扁率与国家参考椭球相同。
其椭球半径α1增大为:α1=α+δα1,δα1=hm+ζ0式中:hm为当地平均海拔高程,ζ0为该地区的平均高程异常。
而地方投影面的确定中,选取过测区中心的经线或某个起算点的经线作为独立中央子午线.以某个特定方便使用的点和方位为地方独立坐标系的起算原点和方位,并选取当地平均高程面hm为投影面。
6、坐标系统间的转换几种常用坐标系和坐标转换介绍GPS直接提供的坐标(B,L,H)是1984年世界大地坐标系(Word Geodetic System 1984即WGS-84)的坐标,其中B为纬度,L为经度,H为大地高即是到WGS-84椭球面的高度。
而在实际应用中,我国地图采用的是1954北京坐标系或者1980西安坐标系下的高斯投影坐标(x,y,),不过也有一些电子地图采用1954北京坐标系或者1980西安坐标系下的经纬度坐标(B,L),高程一般为海拔高度h。
GPS的测量结果与我国的54系或80系坐标相差几十米至一百多米,随区域不同,差别也不同,经粗落统计,我国西部相差70米左右,东北部140米左右,南部75米左右,中部45米左右。
现就上述几种坐标系进行简单介绍,供大家参阅,并提供各坐标系的基本参数,以便大家在使用过程中自定义坐标系。
既然说到了不同的坐标系,就存在坐标转换的问题。
关于坐标转换,首先要搞清楚转换的严密性问题,即在同一个椭球里的坐标转换都是严密的,而在不同的椭球之间的转换这时不严密的。
例如,由1954北京坐标系的大地坐标转换到954北京坐标系的高斯平面直角坐标是在同一参考椭球体范畴内的坐标转换,其转换过程是严密的。
由1954北京坐标系的大地坐标转换到wgs-84的大地坐标,就属于不同椭球体间的转换。
不同椭球体间的坐标转换在局部地区的采用的常用办法是相似变换法,即利用部分分布相对合理高等级公共点求出相应的转换参数。
一般而言,比较严密的是用七参数的相似变换法,即x平移,y平移,z平移,x旋转,y旋转,z旋转,尺度变化k。
要求得七参数就需要在一个地区需要3个以上的已知点,如果区域范围不大,最远点间的距离不大于30km(经验值),这可以用三参数,即x平移,y平移,z平移,而将x旋转,y旋转,z旋转,尺度变化k视为0,所以三参数只是七参数的一种特例。
如果不考虑高程的影响,对于不同椭球体下的高斯平面直角坐标可采用四参数的相似变换法,即四参数(x平移,y平移,尺度变化m,旋转角度α)。
如果用户要求的精度低于20米,在一定范围(2'*2')内,就直接可以用二参数法(δb,δl)或(δx,δy)修正。
但在实际操作中,这也取决于选取的公共点是否合理,并保证其足够的精度。
由于采用的椭球基准不一样,并且由于投影的局限性,使的全国各地并不存在一致的转换参数。
对于这种转换由于量较大,有条件的话,一般都采用GPS联测已知点,应用GPS 软件自动完成坐标的转换。
当然若条件不许可,且有足够的重合点,也可以进行人工解算。
在实际应用中需要将GPS观测成果点位的WGS-84坐标转换为地面网的坐标,首先要把点位的WGS-84坐标转换成国家(或地区)的大地坐标,然后再把大地坐标转换成高斯平面直角坐标。
7、EXCEL完成GPS坐标转换的简易方法作为尖端技术GPS,能方便快捷性地测定出点位坐标,无论是操作上还是精度上,比全站仪等其他常规测量设备有明显的优越性。
随着我国各地GPS差分台站的不断建立以及美国SA政策的取消,使得单机定位的精度大大提高,有的已经达到了亚米级精度,能够满足国土资源调查、土地利用更新、遥感监测、海域使用权清查等工作的应用。
在一般情况下,我们使用的是1954年北京坐标系或1980年西安坐标系(以下分别简称54系和80系),而GPS测定的坐标是WGS-84坐标系坐标,需要进行坐标系转换。
对于非测量专业的工作人员来说,虽然GPS定位操作非常容易,但坐标转换则难以掌握,EXCEL是比较普及的电子表格软件,能够处理较复杂的数学运算,用它来进行GPS坐标转换、面积计算会非常轻松自如。
要进行坐标系转换,离不开高斯投影换算,下面分别介绍用EXCEL进行换算的方法和GPS坐标转换方法。
1)、用EXCEL进行高斯投影换算从经纬度BL换算到高斯平面直角坐标XY(高斯投影正算),或从XY换算成BL(高斯投影反算),一般需要专用计算机软件完成,在目前流行的换算软件中,存在一个共同的不足之处,就是灵活性较差,大都需要一个点一个点地进行,不能成批量地完成,给实际工作带来许多不便。
笔者发现,用EXCEL可以很直观、方便地完成坐标换算工作,不需要编制任何软件,只需要在EXCEL的相应单元格中输入相应的公式即可。
下面以54系为例,介绍具体的计算方法。
完成经纬度BL到平面直角坐标XY的换算,在EXCEL中大约需要占用21列,当然读者可以通过简化计算公式或考虑直观性,适当增加或减少所占列数。
在EXCEL中,输入公式的起始单元格不同,则反映出来的公式不同,以公式从第2行第1列(A2格)为起始单元格为例,各单元格的公式如下:单元格单元格内容说明A2输入中央子午线,以度.分秒形式输入,如115度30分则输入115.30起算数据L0B2=INT(A2)+(INT(A2*100)-INT(A2)*100)/60+(A2*10000-INT(A2*100)*100)/3600 把L0化成度C2以度小数形式输入纬度值,如38°14′20″则输入38.1420起算数据BD2以度小数形式输入经度值起算数据LE2=INT(C2)+(INT(C2*100)-INT(C2)*100)/60+(C2*10000-INT(C2*100)*100)/3600 把B化成度F2=INT(D2)+(INT(D2*100)-INT(D2)*100)/60+(D2*10000-INT(D2*100)*100)/3600 把L化成度G2=F2-B2L-L0H2=G2/57.2957795130823化作弧度I2=TAN(RADIANS(E2))Tan(B)J2=COS(RADIANS(E2))COS(B)K2=0.006738525415*J2*J2L2=I2*I2M2=1+K2N2=6399698.9018/SQRT(M2)O2=H2*H2*J2*J2P2=I2*J2Q2=P2*P2R2=(32005.78006+Q2*(133.92133+Q2*0.7031))S2=6367558.49686*E2/57.29577951308-P2*J2*R2+((((L2-58)*L2+61)*O2/30+(4*K2+5)*M2-L2)*O2/12+1)*N2*I2*O2/2计算结果XT2=((((L2-18)*L2-(58*L2-14)*K2+5)*O2/20+M2-L2)*O2/6+1)*N2*(H2*J2)计算结果Y表中公式的来源及EXCEL软件的操作方法,请参阅有关资料,这里不再赘述。