工程施工过程中,常常会遇到不同坐标系统间,坐标转换的问题。目前国内常见的转换有以下几种:1,大地坐标(BLH)对平面直角坐标(XYZ);2,北京54全国80及WGS84坐标系的相互转换;3,任意两空间坐标系的转换。其中第2类可归入第三类中。所谓坐标转换的过程就是转换参数的求解过程。常用的方法有三参数法、四参数法和七参数法。以下对上述三种情况作详细描述如下:
1,大地坐标(BLH)对平面直角坐标(XYZ)
常规的转换应先确定转换参数,即椭球参数、分带标准(3度,6度)和中央子午线的经度。椭球参数就是指平面直角坐标系采用什么样的椭球基准,对应有不同的长短轴及扁率。一般的工程中3度带应用较为广泛。对于中央子午线的确定有两种方法,一是取平面直角坐标系中Y坐标的前两位3,即可得到对应的中央子午线的经度。如x=3250212m,y=395121123m,则中央子午线的经度=393=117度。另一种方法是根据大地坐标经度,如果经度是在1555~1855度之间,那么对应的中央子午线的经度=(1555+1855)/2=117度,其他情况可以据此3度类推。
另外一些工程采用自身特殊的分带标准,则对应的参数确定不在上述之列。
确定参数之后,可以用软件进行转换,以下提供坐标转换的程序下载。
2,北京54全国80及WGS84坐标系的相互转换
这三个坐标系统是当前国内较为常用的,它们均采用不同的椭球基准。
其中北京54坐标系,属三心坐标系,大地原点在苏联的普而科沃,长轴6378245m,短轴6356863,扁率1/2983;西安80坐标系,属三心坐标系,大地原点在陕西省径阳县永乐镇,长轴6378140m,短轴6356755,扁率1/29825722101;WGS84坐标系,长轴6378137000m,短轴6356752314,扁率1/298257223563。由于采用的椭球基准不一样,并且由于投影的局限性,使的全国各地并不存在一至的转换参数。对于这种转换由于量较大,有条件的话,一般都采用GPS联测已知点,应用GPS软件自动完成坐标的转换。当然若条件不许可,且有足够的重合点,也可以进行人工解算。详细方法见第三类。
3,任意两空间坐标系的转换
由于测量坐标系和施工坐标系采用不同的标准,要进行精确转换,必须知道至少3个重合点(即为在两坐标系中坐标均为已知的点。采用布尔莎模型进行求解。布尔莎公式:
对该公式进行变换等价得到:
解算这七个参数,至少要用到三个已知点(2个坐标系统的坐标都知道),采用间接平差模型进行解算:
其中: V 为残差矩阵;
X 为未知七参数;
A 为系数矩阵;
解之:L 为闭合差
解得七参数后,利用布尔莎公式就可以进行未知点的坐标转换了,每输入一组坐标值,就能求出它在新坐标系中的坐标。 但是要想GPS观测成果用于工程或者测绘,还需要将地方直角坐标转换为大地坐标,最后还要转换为平面高斯坐标。
上述方法类同于我们的间接平差,解算起来较复杂,以下提供坐标转换程序,只需输入三个已知点的坐标即可求解出坐标转换的七个参数。如果已知点的数量较多,可以进行参数间的平差运算,则精度更高。
当已知点的数量只有两个时,我们可以采用简单变换法,此法较为方便易行,适于手算,只是精度受到一定的限制。
详细解算方程如下:
式中调x,y和x\'、y\'分别为新旧(或;旧新)网重合点的坐标,a、b、、k为变换参数,显然要解算出a、b、、k,必须至少有两个重合点,列出四个方程。
即可进行通常的参数平差,解求a、x、b、c、d各参数值。将之代人(3)式,可得各拟合点的残差(改正数)代人(2)式,可得待换点的坐标。
求出解算参数之后,可在Excel中,进行其余坐标的转换。
上次笔者用此法进行过80和54坐标的转换,由于当时没有多余的点可供验证和平差,所以转换精度不得而知,但转换之后各点的相对位置不变。估计,实际的转换误差应该是10m量级的。
还有一些情况是先将大地坐标转换 为直角坐标,然后进行相关转换。
对于坐标系的转换,给很多GPS的使用者造成一些迷惑,尤其是对于刚刚接触的人,搞不明白到底是怎么一回事。我对坐标系的转换问题,也是一知半解,对于没学过测量专业的人来说,各种参数的搞来搞去实在让人迷糊。在我有限的理解范围内,我想在这里简单介绍一下,主要是抛砖引玉,希望能引出更多的高手来指点迷津。
我们常见的坐标转换问题,多数为WGS84转换成北京54或西安80坐标系。其中WGS84坐标系属于大地坐标,就是我们常说的经纬度坐标,而北京54或者西安80属于平面直角坐标。对于什么是大地坐标,什么是平面直角坐标,以及他们如何建立,我们可以另外讨论。这里不多啰嗦。
那么,为什么要做这样的坐标转换呢?
因为GPS卫星星历是以WGS84坐标系为根据而建立的,我国目前应用的地形图却属于1954年北京坐标系或1980年国家大地坐标系;因为不同坐标系之间存在着平移和旋转关系(WGS84坐标系与我国应用的坐标系之间的误差约为80),所以在我国应用GPS进行绝对定位必须进行坐标转换,转换后的绝对定位精度可由80提高到5-10米。简单的来说,就一句话,减小误差,提高精度。
下面要说到的,才是我们要讨论的根本问题:如何在WGS84坐标系和北京54坐标系之间进行转换。
说到坐标系转换,还要罗嗦两句,就是上面提到过的椭球模型。我们都知道,地球是一个近似的椭球体。因此为了研究方便,科学家们根据各自的理论建立了不同的椭球模型来模拟地球的形状。而且我们刚才讨论了半天的各种坐标系也是建立在这些椭球基准之上的。比如北京54坐标系采用的就是克拉索夫斯基椭球模型。而对应于 WGS84坐标系有一个WGS84椭球,其常数采用 IUGG第17届大会大地测量常数的推荐值。WGS84椭球两个最常用的几何常数:长半轴:6378137±2(m);扁率:1:298257223563
之所以说到半长轴和扁率倒数是因为要在不同的坐标系之间转换,就需要转换不同的椭球基准。这就需要两个很重要的转换参数dA、dF。
dA的含义是两个椭球基准之间半长轴的差;dF的含义是两个椭球基准之间扁率倒数的差。
在进行坐标转换时,这两个转换参数是固定的,这里,我们给出在进行84—〉54,84—〉80坐标转换时候的这两个参数如下:
WGS84>北京54:DA:-108;DF:00000005
WGS84>西安80:DA: -3 ;DF: 0
椭球的基准转换过来了,那么由于建立椭球的原点还是不一致的,还需要在dXdYdZ这三个空间平移参量,来将两个不同的椭球原点重合,这样一来才能使两个坐标系的椭球完全转换过来。而由于各地的地理位置不同,所以在各个地方的这三个坐标轴平移参量也是不同的,因此需要用当地的已知点来计算这三个参数。具体的计算方法是:
第一步:搜集应用区域内GPS“B”级网三个以上网点WGS84坐标系B、L、H值及我国坐标系(BJ54或西安80)B、L、h、x值。(注:B、L、H分别为大地坐标系中的大地纬度、大地经度及大地高,h、x分别为大地坐标系中的高程及高程异常。各参数可以通过各省级测绘局或测绘院具有“A”级、“B”级网的单位获得。)
第二步:计算不同坐标系三维直角坐标值。计算公式如下:
X=(N+H)cosBcosL
Y=(N+H)cosBsinL
Z=[N(1-e2)+H]sinB
1、设置四参数在求转换参数界面首先点击右上角的设置按钮,将“坐标转换方法”改为“一步法”,点击“确定”,则可以开始四参数的设置。
2、添加坐标参数“一步法”设置完成后,在求转换参数界面点击“添加”14输入坐标参数平面坐标输入已知点坐标,大地坐标输入RTK测的已知点的坐标。
几何线型变化的话(因为你提到了坐标),因此我建议你采用开源的计算几何库,矩阵变化,平移,向量运算都有
看看这个
>
以上就是关于高分悬赏!! 大地坐标与平面坐标的换算问题,地图学的专家帮忙解惑!全部的内容,包括:高分悬赏!! 大地坐标与平面坐标的换算问题,地图学的专家帮忙解惑!、如何wgs84与当地坐标之间的转换参数、2000坐标系不用求转换参数怎么 *** 作等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)