(总第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
2
(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).