您好,欢迎来到化拓教育网。
搜索
您的当前位置:首页03 WGS_84和BJ54坐标转换问题的分析

03 WGS_84和BJ54坐标转换问题的分析

来源:化拓教育网
(总第161期)2008年第4期安徽建筑

169

WGS-84和BJ54坐标转换问题的分析

TheAnalysisofCoordinateTransformationBetweenWGS-84andBJ54

(安徽省建设工程勘察设计院,安徽合肥230001)

要:文章介绍了从GPS定位结果至平面坐标的两种转换模型,即平面转换模型和空间转换模型。通过实例,比较了它们的转换精度。并得出

了一些有益的结论。

关键词:WGS-84;BJ54;平面转换;空间转换中图分类号:TU198

文献标识码:B

文章编号:1007-7359(2008)04-0169-02

ZhangCan(AnhuiGeotechnicalInvestigationandArchitecturalDesignInstitute,Hefei230001,China)

Abstract:ThepaperintroducestwoalgorithmsfortransferringtheGPSresulttolocalgridframe,namelytheplanetransformmodelandthespacetransformmodel.Thetransformprecisioniscomparedinexamplesandsomebeneficialresultsareconcluded.Keywords:WGS-84;BJ54;planetransformmodel;spacetransformmodel

0引言其中:(x'g.y'g)T为假定的地方平面坐标;

随着GPS定位精度的发展,GPS技术在测量中的应用也越来越广泛。但是GPS定位的结果是位于WGS-84坐标系统中,

(xg,yg)T为地方平面坐标;(x0,y0)T为坐标平移量;

r为缩放尺度;

)sin(ψ)cos(ψ

为旋转矩阵;R(ψ)=

-sin(ψ)cos(ψ)

WGS-84坐标是一个协议地球参考系,其坐标原点位于地球的

质心上。而我国测绘成果普遍表示在BJ54和1980年国家大地

坐标系下,它是参心坐标系,BJ54生标系是以克拉索夫斯基椭球为参考椭球,采用高斯投影方式进行投影。因此,需要将GPS测量得到的位于WGS-84坐标系的成果转换到BJ54坐标系下。在进行WGS-84坐标系和BJ54坐标系转换时有两种转换思想和模型,即平面转换模型和空间转换模型。

&(ψ为旋转角。

为反算式缩放尺度和旋转参数,至少需要两(3)中的平移、个已知WGS-84坐标和BJ54坐标的点,如多于两个点,可按最小二乘法求解。

对所有GPS测定的点需经过以上4个过程求得当地平面坐标。而正常高需要另外计算,可由大地高扣除高程异常求得。

平面转换模型原理简单,数值稳定可靠,但由于(3)式是一个线性变换公式,而高斯投影是非线性的,因此平面转换模型只适合范围较小的工程使用,对于大范围的GPS测量应使用空间转换模型。

1平面转换模型

假设BJ54椭球的中心和坐标轴与WGS-84椭球一致,可

通过平面转换模型,将GPS定位得到的大地经纬度和大地高通过以下过程转换为平面坐标。

①在WGS-84椭球参数约束下将WGS-84大地坐标转换为WGS-84空间直角坐标;

!###\"###$

X=(N+H)cosBcosLY=(N+H)cosBsinLZ=[N(1-e2)+H]sinB

(1)

空间转换模型

在进行空间转换时,首先必须假定WGS-84坐标测定点中

至少有3个已知的BJ54坐标系的平面坐标。根据公共点的坐标首先进行7参数的测定,然后进行WGS-84坐标和BJ54坐标的转换。

(1)式将WGS-84大地①在WGS-84椭球参数约束下通过坐标转换成WGS-84空间直角坐标;

(2)

(4)式,将WGS-84空间直角坐标转②通过空间转换模型换成地方空间直角坐标;

②将WGS-84空间直角坐标等同于BJ54坐标系下的空间

直角坐标,将其在BJ54椭球参数约束下转换为假定的大地坐标;

!

#####\"#####$

L=arctan(Y)

XB=arctan[(Z+Ne2sinB)/%X2+Y2]H=%X2+Y2secB-N

③在地方椭球参数约束下,将地方空间直角坐标转换成地

方大地坐标;

③取当地子午线,将假定的地方大地坐标通过高斯

投影转换成假定地方平面坐标;

④取当地子午线,将地方大地坐标高斯投影成地方平

面坐标。

空间转换模型见式(4)。

)

*******+

④通过平面转换模型将假定的地方平面坐标转换成地方

平面坐标。平面转换模型见式(3)。

&’&(xgyg

=x0y0

+(1+r)R(ψ)

&(x'g

y'g

(3)

X地方,--Y地方-=----Z地方.

)

*******+

)

Dx,X84*-*-*-Dy-+(1+k)R(α)R(β)R(γ)*Y84

*-*-*-Dz.Z84+

,

-------.

(4)

收稿日期:2008-04-16作者简介:张灿(1974-),男,安徽六安人,毕业于同济大学,硕士,工程师。

其中,[Dx,Dy,Dz]T是空间转换坐标平移量;

k为缩放尺度参数;

170

张灿WGS-84和BJ54坐标转换问题的分析第15卷第4期

为旋转参数。,β,γα

为反算出式缩放尺度和旋转参数,至少需要(4)中的平移、

点名

检测点已知WGS-84坐标和已知BJ54坐标

WGS-84

(度分秒)B

(度分秒)L

(m)H

(m)x

表1

BJ54

(m)y

3个已知WGS-84坐标和地方坐标的点。如多于3个点,可按

最小二乘进行拟合求解。

空间转换模型适合于大范围GPS测量,但在实际施工过程中,根据施工精度的要求,又分为3种情况:在式(4)中,需求解7个参数,故称为七参数转换模型;若不需求尺度参数,则称为六参数转换;若尺度参数和旋转参数均不求,则称为三参数转换。对于7参数模型的求解,至少需要3个公共点,6参数模型的求解也至少需要3个公共点,因为尽管两个公共点有6个坐标分量,按式(4)可以列出6个观测方程,但这6个坐标分量中,由于两点间的距离是固定的,只有5个是的,3参数模型可在只有1个公共点的情况下求解。

通过空间转换模型,可以同时求得平面坐标和高程,不过按上面步骤求得的是地方坐标系下的大地高。而通常情况下在求解空间转换参数时,使用正常高代替地方坐标系下的大地高来换算地方空间直角坐标,这样GPS测定的WGS-84坐标就直接转换为地方坐标系下的平面坐标和正常高。用正常高代替大地高时,尽管转换模型对高程并不敏感,转换出来的平面坐标不会有多大的影响,但大地水准面的不规则性将体现在高程转换残差中,这种方法在大地水准面不复杂的地区使用比较有效。但在大地水准面复杂的地区转换出来的正常高精度较低,实际应用时可以通过拟合模型,对高程转换残差进行拟合,在转换出来的高程上加一个修正量以获得更精确的正常高。

C1C2C3C4

点名点名

C1C2C3C4

315359.82226315336.69010315222.88286315135.26930

1171432.248541171620.458631171630.187951171821.99622

12.53814.68249.2215.640

3530834.3633530128.5863527855.8803529671.753

522853.259525698.415525959.830522045.889

空间转换模型转换值与已知值比较

已知BJ54

转换BJ54

差值

表2

x(m)

C1C2C3C4

3530834.3633530128.5863527855.8803529671.753

y(m)522853.259525698.415525959.830522045.889

x(m)3530834.3613530128.5833527855.8763529671.758

y(m)522853.261525698.417525959.830522045.906

x(cm)0.200.300.42-0.49

y(cm)-0.20-0.200.01-1.68

平面转换模型转换值与已知值比较

已知BJ54

转换BJ54

差值

表3

x(m)3530834.3633530128.5863527855.8803529671.753

y(m)522853.259525698.415525959.830522045.889

x(m)3530834.3693530128.5793527855.8773529671.762

y(m)522853.264525698.413525959.824522045.901

x(cm)-0.600.700.30-0.90

y(cm)-0.500.200.60-1.20

3实例分析

现已知7个公共点A1、A2、A3、B1、B2、B3、B4的WGS-84

可见,在用来求解转换参数的公共点坐标精度得到保障的情况下,通过空间转换模型和平面转换模型,都可以得到较高的转换精度,都能适合于实际的生产测绘作业。

坐标和BJ54平面坐标,而且这7个公共点的等级相同。测区范围内有4个检测点C1、C2、C3、C4,这4个点的WGS-84坐标和BJ54坐标已知,首先将这11个点的BJ54平面坐标展在一副工作图上(如图1所示),为了检验平面转换模型和空间转换模型对转换结果的影响,选择A1、A2、A3三个控制范围较大的点建立空间转换模型,B1、B2、B3、B4四个点建立平面转换模型。通过空间和平面转换模型,比较4个检测点C1、C2、C3、C4转换值和已知值的差异,从而评价空间和平面转换模型的转换精度。

A1

4结束语

通过对WGS-84与BJ54坐标转换问题的分析,我们可以

获得如下一些有用的结论:

①GPS测定点可以通过先投影,再用平面转换模型转换到

当地平面坐标,高程可以用高程拟合的方法获取,即平面与高程分别转换得到。这种转换模型数值上稳定,但含有高斯投影变形的影响,只适合测区范围较小的情况;

②GPS测定点通过空间转换,可同时得到平面坐标和高

程,通常在测区范围较大情况下使用。在大地水准面不复杂的地区,高程转换的精度也能得到保障;

B1

B4

C1C2

图1C3已知坐标展点图C4B2

③在用来求解转换参数的公共点坐标精度得到保障的情

况下,平面转换模型和空间转换模型都能得到很高的精度;

④转换的精度除了取决于坐标转换的思想和转换模型外,还

与求解转换参数的公共点坐标精度、公共点的点数和分布有关。

B3

A3

A2

图1

已知坐标展点图

参考文献

[1][2][3][4]

李连伟,荣艳妮.WGS84和BJ54坐标转换问题探讨[J].测绘与空间地理信息,2004(2).蒋

辉,等.北京54坐标系坐标转换的探讨[J].南京工业大学学报,2007(7).

丁士俊,张忠明.几种不同坐标变换方法问题的研究[J].四川测绘,

鉴于54坐标点的高程异常通常不能精确得知,故对这两种转换方法得到的平面坐标的精度进行了比较。

方法一:通过A1、A2、A3三个控制范围较大的点建立空间

7参数转换模型,转换值与已知值的比较列于表2。

方法二:通过B1、B2、B3、B4四个控制范围较小的点建立

平面4参数转换模型,转换值与已知值的比较列于表3。

2005(1).

徐仕祺,张晓帆,等.关于利用七参数法进行WGS-84和BJ-54坐标转换问题的探讨[J].测绘与空间地理信息,2007(10).

因篇幅问题不能全部显示,请点此查看更多更全内容

Copyright © 2019- huatuo9.cn 版权所有 赣ICP备2023008801号-1

违法及侵权请联系:TEL:199 18 7713 E-MAIL:2724546146@qq.com

本站由北京市万商天勤律师事务所王兴未律师提供法律服务