GPS测量中坐标系统、坐标系地转换过程问题详解详解

GPS测量中坐标系统、坐标系的转换过程详解

摘要：GPS在测量领域得到了广泛的应用，本文介绍将GPS所采集到的WGS-84坐标转换成工程所需的坐标的过程。

关键词：GPS 坐标系统坐标系转换

一、概述GPS及其应用

GPS即全球定位系统(Global Positioning System)是美国从本世纪70年代开始研制，历时20年，耗资200亿美元，于1994年全面建成的卫星导航定位系统。作为新一代的卫星导航定位系统经过二十多年的发展，已成为在航空、航天、军事、交通运输、资源勘探、通信气象等所有的领域中一种被广泛采用的系统。我国测绘部门使用GPS也近十年了，它最初主要用于高精度大地测量和控制测量，建立各种类型和等级的测量控制网，现在它除了继续在这些领域发挥着重要作用外还在测量领域的其它方面得到充分的应用，如用于各种类型的工程测量、变形观测、航空摄影测量、海洋测量和地理信息系统中地理数据的采集等。G PS以测量精度高; 操作简便，仪器体积小，便于携带; 全天候操作;观测点之间无须通视;测量结果统一在WGS84坐标下，信息自动接收、存储，减少繁琐的中间处理环节、高效益等显著特点,赢得广大测绘工作者的信赖。

二、GPS测量常用的坐标系统

1．WGS-84坐标系

WGS-84坐标系是目前GPS所采用的坐标系统，GPS所发布的星历参数就是基于此坐标系统的。WGS-84坐标系统的全称是World Geodical System-84（世界大地坐标系-84），它是一个地心地固坐标系统。WGS-84坐标系统由美国国防部制图局建立，于1987年取代

了当时GPS所采用的坐标系统―WGS-72坐标系统而成为GPS的所使用的坐标系统。WGS -84坐标系的坐标原点位于地球的质心，Z轴指向BIH1984.0定义的协议地球极方向，X轴指向BIH1984.0的启始子午面和赤道的交点，Y轴与X轴和Z轴构成右手系。采用椭球参数为：a = 6378137m f = 1/298.257223563

2．1954年北京坐标系

1954年北京坐标系是我国目前广泛采用的大地测量坐标系，是一种参心坐标系统。该坐标系源自于原苏联采用过的1942年普尔科夫坐标系。该坐标系采用的参考椭球是克拉索夫斯基椭球，该椭球的参数为：a = 6378245m f = 1/298.3。我国地形图上的平面坐标位置都是以这个数据为基准推算的。

3．地方坐标系（任意独立坐标系）

在我们测量过程中时常会遇到的如一些某城市坐标系、某城建坐标系、某港口坐标系等，或我们自己为了测量方便而临时建立的独立坐标系。

三、坐标系统的转换

在工程应用中使用GPS卫星定位系统采集到的数据是WGS-84坐标系数据,而目前我们测量成果普遍使用的是以1954年北京坐标系或是地方（任意）独立坐标系为基础的坐标数据。因此必须将WGS-84坐标转换到BJ-54坐标系或地方（任意）独立坐标系。

目前一般采用布尔莎公式（七参数法）完成WGS-84坐标系到北京54坐标系的转换，得到北京54坐标数据。

XBJ54=XWGS84+ KXWGS84+Δx+YWGS84ξZ＂/ρ＂－ＺＷＧＳ８４ξＹ＂/ρ＂

YBJ54=YWGS84+ KYWGS84+ΔY-XWGS84ξZ＂/ρ＂+ＺＷＧＳ８４ξX＂/ρ＂

ZBJ54=ZWGS84+ KZWGS84+ΔZ+XWGS84ξY＂/ρ＂－ＺＷＧＳ８４ξX＂/ρ＂

四、坐标系的变换

同一坐标系统下坐标有多种不同的表现形式，一种形式实际上就是一种坐标系。如空间直角坐标系（X，Y，Z）、大地坐标系（B，L）、平面直角坐标（x,y）等。通过坐标统的转换我们得到了BJ54坐标系统下的空间直角坐标，我们还须在BJ54坐标系统下再进行各种坐标系的转换，直至得到工程所需的坐标。

1．将空间直角坐标系转换成大地坐标系，得到大地坐标(B,L)：

L=arctan(Y/X)

B=arctan {(Z+Ne2sinB)/(X2+Y2)0.5}

H=(X2+Y2)0.5sinB-N

用2．将大地坐标系转换成高斯坐标系，得到高斯坐标（x,y）

按高斯投影的方法求得高斯坐标，x=F1(B,L),y=F2(B,L)

3．将高斯坐标系转换成任意独立坐标系，得到独立坐标（x’,y’）

在小范围内测量，我们可以将地面当作平面，用简单的旋转、平移便可将高斯坐标换成工程中所采用坐标系的坐标(x’,y’)，

x’=xcosα+ysinα

y’=ycosα-xsinα

五、小结

由于GPS测量的种种优点，GPS 定位技术现已基本上取代了常规测量手段成为了主要的技术手段，市面上出现了许多转换软件和不同型号的GPS数据处理配套软件（包含了怎样将GPS测量中所得到的WGS-84转换成工程中所须坐标的功能），万变不离其宗，只要我们明白了WGS-84转换到独立坐标系的转换过程，便可很容易的使用该软件了，甚至可以自己编写程序，将WGS-84坐标转换成独立坐标系坐标。

本文主要是介绍坐标系统、坐标系的转换过程，文中提及的符号及具体转换方法请参阅相关文献。

上式采用迭代法求出大地坐标（B，L）

坐标转换工具说明书-1208

§10.2坐标转换工具 HGO 数据处理软件包提供了坐标转换程序，可以进行地方坐标与WGS-84坐标的相互转换，同时具备参数求解功能。 下面对这个工具进行介绍： 10.2.1概述 首先，介绍一下常见的三种坐标表示方法：经纬度和椭球高（BLH），空间直角坐标（XYZ），平面坐标和水准高程（xyh/NEU）。注意：椭球高是一个几何量，而水准高是一个物理量。 我们通常说的WGS-84坐标是经纬度和椭球这一种，北京54坐标是平面坐标和水准高程这一种，实质是有平面基准和高程基准组成的。 此外，再注意一下坐标转换的严密性问题，在同一个椭球里的纯几何转换都是严密的（BLH<->XYZ），而在不同的基准之间的转换是不严密的。举个例子，在WGS-84坐标和北京54坐标之间是不存在一套转换参数可以全国通用的，因为前者是一个地心坐标系，后者是一个参心坐标系。高程转换是由几何高向物理高转换。因此在每个地方必须用椭球进行局部拟合，通常用7参数模型来拟合。 那么，两个椭球间的坐标转换应该是怎样的呢？一般而言比较严密的是用七参数法（或称布尔莎模型），即X平移，Y平移，Z平移，X旋转，Y旋转，Z旋转，尺度变化K。要求得七参数就需要在一个地区需要3个以上的已知点（7个参数至少7个方程可解，所以需要三个点列出9个方程），如果区域范围不大、最远点间的距离不大于30Km(经验值)的情况可以用三参数，即X平移，Y平移，Z平移，而将X旋转，Y旋转，Z旋转，尺度变化K视为0，所以三参数只是七参数的一种特例。 七参数模型的实质是用一个局部椭球去拟合地方坐标系的形态；所以转换后获得的地方椭球高就是水准高。当然我们也可以把平面和高程两个方向分别进行拟合。例如平面用四参数模型拟合，高程方向则用二次曲面等模型来拟合。这样分开处理的模式相对七参数模型自由度更高。但是由于四参数模型参数较少，表达能力较弱，通常只用于小区域坐标转换。 综上所述，从实用的角度出发，坐标转换程序提供了两种转换策略供给客户选择使用： 1.七参数模型，一步得到地方平面和水准数据。 2.四参数加高程拟合模型，分两步得到地方平面和水准数据。 由于各厂家的模型和流程定义可能是不一样的，这里就我们公司的转换流程描述如下：七参数的转换过程是这样的：

不同类型地图使用的投影与坐标系

不同类型地图使用的投影与坐标系 (2016-08-12 15:29:29) 不同类型地图使用的投影与坐标系 1.概念辨析 地图投影跟大地坐标系是完全两个东西，尽管具有相关性。地球椭球体则是另一个东西。实际上地图编绘涉及三个基本的东西：椭球体、地图投影、大地坐标系。三者密切关联。（百科知识） 要绘制地图，首先考虑用什么椭球体，这是投影和坐标系的基础——我国三代坐标系使用三种椭球体。 三者之间的关系：先有个椭球体，然后是投影到承影面，然后是添加经纬网。椭球体是基础，投影是转换函数，是数学关系，大地坐标系是参照系。因此，同一椭球体可以用不同的投影；而同一投影，也可以用不同的大地坐标系。 但是一般三者是协调一致的，如我国的三代坐标系，有对应的椭球体、投影类型、基准面（坐标系）。 从地图反映地球表面来看，整个过程涉及五个环节：地球～椭球体～投影～坐标系～地图。而地球是球面的，是一个曲面，而地图是平面的，二者的结构性矛盾，导致我们不得不采用一系列转换，这个转换中不可避免地产生扭曲、变形和误差。具体关系：总结：地球（地球表面，存在高低起伏）→椭球体（光滑球面，相关参数）→投影（投影方式：几何投影与解析投影）→坐标系（地理坐标系与平面直角坐标系）→地图。 2. 我国三代坐标系 我们经常给影像投影时用到的北京54、西安80和2000坐标系是投影直角坐标系，如下表所示为国内坐标系采用的主要参数。从中可以看到我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国的大地基准面。 表：北京54、西安80和2000坐标系参数列表 坐标名称投影类型椭球体基准面 北京54Gauss Kruger （Transverse Mercator） Krasovsky D_Beijing_1954 西安80Gauss Kruger （Transverse Mercator） IAG75D_Xian_1980 CGCS2000Gauss Kruger CGCS2000D_China_2000

施工坐标系与测量坐标系之间的相互转换关系

施工坐标系与测量坐标系之间的相互转换 一、用Microsoft Excel 编辑转换 如图（1-1）所示：设Y O X -- 为测量坐标系，y o x -'- 为施工坐标，如果知道了施工坐标系的原点o '的测量坐标为（'0X ,'0Y ）、定向点I 的测量坐标为（XI,YI ），定向坐标方位角 x -'0α （即纵轴的旋转角，因为0=-X o α为正北方向，则x -'0α=X o -α+α）。则所求P 点由施工坐标P （p p y x ,）换算成为测量坐标P （p p Y X ,）的公式则为： α αsin *cos *0p p p y x X X -+=' ααcos *sin *0p p p y x Y Y ++=' 上面两式在Excel 中编辑公式为： [][]180/()*sin *180/()*cos *0Pi y Pi x X X p p p αα-+=' [][]180/()*cos *180/()*sin *0Pi y Pi x Y Y p p p αα++=' 而如果知道了施工坐标系（第二坐标系）的原点的测量坐标 o '为

（'0X 、'0Y ）、坐标方位角 x -'0α （即纵轴的旋转角，因为0=-X o α为正北方向，则x -'0α=X o -α+α）。则所求P 点由测量坐标P （p p Y X ,）转换算为施工坐标P （p p y x ,）其公式为： ααsin *)(cos *)(00''-+-=Y Y X X x p p p ααcos *)(sin *)(00''-+--=Y Y X X y p p p 上面两式在Excel 中编辑公式为： [][]180/()*sin *)(180/()*cos *)(00Pi Y Y Pi X X x p p p αα''-+-= [][]180/()*cos *)(180/()*sin *)(00Pi Y Y Pi X X y p p p αα''-+--= 以上各式中施工坐标系原点o ' 的测量坐标（'0X ，'0Y ）与方位角α ，可在设计资料中查找或用图解法得出。 附： 如（图1-2）直线AB 的坐标方位角 ? ?? ? ??--=-A B A B AB x x y y 1tan α B ( x ,y ) β B B C ( x ,y ) C C A ( x ,y ) A A α A B α A C 图（1-2） 如（图1-2）直线AB 与直线AC 的夹角 β ???? ??---???? ??--=-=--A B A B A C A C A B A C x x y y x x y y 11tan tan ααβ

坐标变换就是两种坐标类型

坐标变换就是两种坐标类型、不同参照体系之间的变换 坐标变换因不同的坐标类型、体系变换方法不一样，没有固定的公式 比方说测量地球，就有多种坐标体系： 1。以地心为原点的空间直角坐标 2。经纬度坐标 3。把地球表面分成很多格子，对于一个小格子区，球面接近平面，在这个平面上设一个平面直角坐标系，就是北京54坐标等坐标形式 这些坐标来回转换，比较复杂，甚至是学术性的问题，一般根据不同的观点和精度，有一些小程序，做转换工作 工程施工过程中，常常会遇到不同坐标系统间，坐标转换的问题。目前国内常见的转换有以下几种：1，大地坐标（BLH）对平面直角坐标（XYZ）；2，北京54全国80及WGS84坐标系的相互转换；3，任意两空间坐标系的转换。其中第2类可归入第三类中。所谓坐标转换的过程就是转换参数的求解过程。常用的方法有三参数法、四参数法和七参数法。以下对上述三种情况作详细描述如下： 1，大地坐标（BLH）对平面直角坐标（XYZ） 常规的转换应先确定转换参数，即椭球参数、分带标准（3度，6度）和中央子午线的经度。椭球参数就是指平面直角坐标系采用什么样的椭球基准，对应有不同的长短轴及扁率。一般的工程中3度带应用较为广泛。对于中央子午线的确定有两种方法，一是取平面直角坐标系中Y坐标的前两位*3，即可得到对应的中央子午线的经度。如x=3250212m， y=395121123m，则中央子午线的经度=39*3=117度。另一种方法是根据大地坐标经度，如果经度是在155.5~185.5度之间，那么对应的中央子午线的经度=（155.5+185.5）/2=117度，其他情况可以据此3度类推。 另外一些工程采用自身特殊的分带标准，则对应的参数确定不在上述之列。 确定参数之后，可以用软件进行转换，以下提供坐标转换的程序下载。 2，北京54全国80及WGS84坐标系的相互转换 这三个坐标系统是当前国内较为常用的，它们均采用不同的椭球基准。

ARCGIS中坐标转换

ArcGIS 坐标转换 1.坐标分析 问题：对于某地A中心点坐标为455299.845,3223622.525的CAD矩形，CAD施工图。将其转换为WGS-84坐标，如何转换？ 分析：分析455299.845为6位，则为东向Y坐标，省去了带号，加上了5000000加常数，其最大为为4，说名在中央子午线的左侧（左侧为负值，加上500万后肯定小于500万，首位为4。若在中央子午线右侧，则最大位数为5）；3223622.525为7位，为北向X坐标。 查看“某地A”的经度为92.5度，因为为CAD施工图，比例尺肯定大于1:5万，所以为3度带，所以此点的中央子午线为93E，带号为Beijing_54_Zone_31。 2.CAD转为shp格式并设定坐标系： ArcTool box-Convesion Tools->To Geodatabse->CAD to Geodatabase: 其中空间参考坐标系选择Beijing_1954_3_Degree_GK_CM_93E。 具体原因：选择投影坐标系-Gauss Kruger-Bei Jing54，此时3度带有两种：Beijing_1954_3_Degree_GK_CM_93E和Beijing_54_Zone_31，前者表示中央子午线为93E的3度带，后者表示北京54 31度带，二者意义一样，但选择哪种呢？因为点坐标东向为455299.845为6位，不带带号，因此选择Beijing_1954_3_Degree_GK_CM_93E（若东向坐标

为31455299.845，则选择Beijing_54_Zone_31）， 3.北京54到WGS84坐标的转换 1.1加载图层： 打开ArcTool box-Data Management Tools->Project and transformation->feature->Project，加载shp图层，弹出下列窗口： 出现红色“X”号，说明原始图层坐标系没有识别出，则需要首先设定其坐标系后再转换。具体设坐标系参考“9 设置或改变Shp文件坐标系” 1.2选择输出图层地址和名称： 在Out Put Dataset or Feature处输入输出图层名：

坐标系间的转换

坐标系间的转换 针对西安80坐标系和北京54坐标系之间椭球参数的转换，采用七参数布尔莎模型，进行不同坐标系之间的坐标转换。 标签：七参数布尔莎模型参考椭球MAPGIS平台 0 引言 我们现在改用的西安80坐标系与以前的北京54坐标系的参考椭球体参数是不相同的。54坐标系转换成80坐标系由于椭球参数、定位和定向的变化，必然引起地形图的图廓线、方里线位置以及地形图内地形、地物相关位置的改变。为此，若同时使用根据两种坐标系测制的地形图的情况下，一定要涉及到54坐标系向80坐标系转换问题。转换的原理和方法：大地坐标系变更后，国家基本系列地形图的变更和处理，必须在高斯平面内进行。由于新旧椭球参数不同，参心所在位置也不同，反映在高斯平面上，在同一个投影带里，它们的纵横坐标轴不重合，因此，地面上某一点经过不同椭球面而投影到高斯平面上，它距两系统坐标轴之距离是不等的，在X轴和Y轴上必定都有一个差值。我们按照一定的数学法则将地球面上的经纬网转换到平面上，使地面的地理坐标与平面直角坐标建立起函数关系，实现由曲面向平面的转化。常用的投影大概有二三十种，投影的选取要考虑地图的用途，投影的形变大小等众多因素。 1 北京54坐标系与西安80坐标系 1.1 54国家坐标系：是我国建国初期，为了迅速开展我国的测绘事业，鉴于当时的实际情况，将我国一等锁与原苏联远东一等锁相连接，然后以连接处呼玛、吉拉宁、东宁基线网扩大边端点的原苏联1942年普尔科沃坐标系的坐标为起算数据，平差我国东北及东部区一等锁，这样传算过来的坐标系就定名为1954年北京坐标系。因此，P54可归结为：①属参心大地坐标系；②采用克拉索夫斯基椭球的两个几何参数；③大地原点在原苏联的普尔科沃；④采用多点定位法进行椭球定位；⑤高程基准为1956年青岛验潮站求出的黄海平均海水面；⑥高程异常以原苏联1955年大地水准面重新平差结果为起算数据。按我国天文水准路线推算而得。 自P54建立以来，在该坐标系内进行了许多地区的局部平差，其成果得到了广泛的应用。 1954北京坐标系参考椭球基本几何参数 长半轴a＝6378245m 短半轴b＝6356863.0188m

参考系坐标系及转换汇总

1 天球坐标系、地球坐标系和卫星测量中常用的坐标系的建立方法。天球直角坐标系 天球坐标系 天球球面坐标系 坐标系 地球直角坐标系 地球坐标系 地球大地坐标系 常用的天球坐标系：天球赤道坐标系、天球地平坐标系和天文坐标系。在天球坐标系中，天体的空间位置可用天球空间直角坐标系或天球球面坐标系两种方式来描述。 1 天球空间直角坐标系的定义 地球质心O为坐标原点，Z轴指向天球北极，X轴指向春分点，Y轴垂直于XOZ平面，与X轴和Z轴构成右手坐标系。则在此坐标系下，空间点的位置由坐标（X，Y，Z）来描述。 春分点：当太阳在地球的黄道上由天球南半球进入北半球，黄道与赤道的交 点）．

2 天球球面坐标系的定义 地球质心O为坐标原点，春分点轴与天轴（天轴：地球自转的轴）所在平面为天球经度（赤经）测量基准——基准子午面，赤道为天球纬度测量基准而建立球面坐标。空间点的位置在天球坐标系下的表述为（r，α，δ）。

表示：2-1天球空间直角坐标系与天球球面坐标系的关系可用图

岁差和章动的影响 岁差：地球实际上不是一个理想的球体，地球自转轴方向不再保持不变，这 使春分点在黄道上产生缓慢的西移，这种现象在天文学中称为岁差。章动：在日月引力等因素的影响下，瞬时北天极将绕瞬时平北天极旋转，大致呈椭圆，这种现象称为章动。 极移：地球自转轴相对地球体的位置并不是固定的，因而，地极点在地球表面上的位置，是随时间而变化的，这种现象称为极移。地球的自转轴不仅受日、月引力作用而使其在空间变化，而且还受地球内部质量不均匀影响在地球内部运动。前者导致岁差和章动，后者导致极移。 协议天球坐标系：为了建立一个与惯性坐标系统相接近的坐标系，人们通常选择某一时刻，作为标准历元，并将此刻地球的瞬时自转轴（指向北极）和地心至瞬时春分点的方向，经过瞬时的岁差和章动改正后，分别作为X轴和Z轴的指向，。协议天球坐标系由此建立的坐标系称为 3 地球坐标系

常用坐标系与高程系简介

常用坐标系与高程系简介 2009-09-27 10:06:45| 分类：GIS技术| 标签：|字号大中小订阅 坐标系的概念 1.坐标系的定义： 如果空间上任意一点P的位置，可以用一组基于某一时间系统时刻t的空间结构的数学描述来确定，则这个空间结构可以称为坐标系，数学描述称为P点在该坐标系中的坐标。牛顿运动学原理要求坐标系是惯性的，惯性是每个物体所固有的当没有外力作用时保持静止或匀速直线运动的属性，基于这个特性，惯性坐标系的定义需与时间无关，通常这样的坐标系需要三个属性来描述(这应该是三维空间的本性吧)，首先一个是原点(O)，就是坐标系的中心点，第二个是过原点的任意直线(这里称为Z轴)，第三个是过原点且与Z轴不重合的任意直线(这里称为X轴)，如果X轴与Z轴垂直，会带来较优美的数学描述，我们称这样的坐标系是笛卡尔坐标系。P点的位置可以用P到原点的距离r，OP与Z轴的夹角，OP与X 轴的夹角来描述(当然也可以有其它等价描述)，可以证明这个描述确定的P点是唯一的。 2.GPS领域常用坐标系模型： 在GPS测量中，最常用的坐标系模型是协议地球坐标系，该坐标系随同地球一起旋转，讨论随地球一起自转的目标位置，用这类坐标系方便；另外一类是协议天球坐标系，这个坐标系随同太阳系一同旋转，与地球自转无关，讨论卫星轨道运动时，用这类坐标系方便。 天球坐标系的定义是这样的，原点是地球质心(O)，Z轴指向地球自转轴(天极，向北为正)，X轴指向春分点，根据春分点的定义可以证明X轴与Z轴互相垂直，且X轴在赤道面上，同时为数学描述方便，引入与XOZ成右手旋转关系的Y轴。因为地球自转轴受其它天体影响(日、月)在空间产生进动，使得春分点变化(章动和岁差)，导致用“瞬时天极”定义的坐标系不断旋转，而旋转的坐标系表现出非惯性的特性，不能直接应用牛顿定律。我们可以用某一历元时刻的天极和春分点(协议天极和协议春分点)定义一个三轴指向不变的天球 坐标系，称为固定极天球坐标系。 地球坐标系的定义是这样的，原点为地球质心(O)，Z轴为地球自转轴，X轴指向地球上赤道的某一固定“刚性”点，所谓“刚性”是指其自转速度与地球一致，同时也为数学描述方便，引入与XOZ成右手旋转关系的Y轴。地球不是一个严格刚性的球体，Z轴在地球上随时间而变，称为极移，同天球坐标系一样，需要指定一个固定极为Z轴，这样的地球坐标系称为固定极地球坐标系。可以证明当观察地球上的物体时，该坐标系是惯性的。如果一个坐标系OXYZ，O不是地球质心，Z轴与地球自转轴平行，则这个坐标系具有与地球相同的自转角速度，我们也把此类坐标系称为地球坐标系。 3.协议坐标系统： 那么，什么是“协议”坐标系呢？通常，理论上坐标系由定义的坐标原点和坐标轴指向来确定。坐标系一经定义，任意几何点都具有唯一一组在该坐标系内的坐标值，反之，一组该坐标系内的坐标值就唯一定义了一个几何点。实际应用中，在已知若干参考点的坐标值后，通过观测又可反过来定义该坐标系。可以将前一种方式称为坐标系的理论定义。而由一系列已知点所定义的坐标系称为协议坐标系，这些已知参考点构成所谓的坐标框架。在点位坐标值不存在误差的情况下，这两种方式对坐标系的定义是一致的。事实上点位的坐标值通常是通过一定的测量手段得到，它们总是有误差的，由它们定义的协议坐标系与原来的理论定义的坐标系会有所不同，凡依据这些点测定的其它点位坐标值均属于这一协议坐标系而不属于理论定义的坐标系。由坐标框架定义的固定极天球坐标系和固定极地球坐标系，称为协议天 球坐标系和协议地球坐标系。

大地测量坐标系统及其转换

大地测量坐标系统及其转换 雷伟伟 河南理工大学测绘学院 wwlei@https://www.wendangku.net/doc/8811770729.html,

基本坐标系 1、大地坐标系 坐标表示形式：(, ,)L B H 大地经度L ：地面一点P 地的大地子午面N P S 与起始大地子午面所构成的二面角； 大地纬度B ：P 地点对椭球面的法线P P K 地与赤道面所夹的锐角； 大地高H ：P 地点沿法线到椭球面的距离。 赤道面 S W 2、空间直角坐标系 坐标表示形式：(,,)X Y Z 以椭球中心O 为坐标原点，起始子午面N G S 与赤道面的交线为X 轴，椭球的短轴为Z 轴（向北为正），在赤道面上与X 轴正交的方向为Y 轴，构成右手直角坐标系O X YZ 。

Y W 3、子午平面坐标系 坐标表示形式：(,,) L x y 设P点的大地经度为L，在过P点的子午面上，以椭圆的中心为原点，建立x、y平 面直角坐标系。则点P的位置用(,,) L x y表示。 x

坐标表示形式：(,,)L u H 设椭球面上的点P 的大地经度为L 。在此子午面，以椭球中心O 为圆心，以椭球长半径a 为半径，做一个辅助圆。过P 点做一纵轴的平行线，交横轴于1P 点，交辅助圆于2P 点，连结2P 、O 点，则21P O P 称为P 点的归化纬度，用u 来表示。P 点的位置用(,)L u 表示。 当P 点不在椭球面上时，则应将P 沿法线投影到椭球面上，得到点0P ，0PP 即为P 点的大地高，0P 点的归化纬度，就是P 点的归化纬度。P 点的位置用(,,)L u H 表示。 x y P u 点在椭球面上时的 P u 点不在椭球面上时的x

参考系坐标系及转换

1天球坐标系、地球坐标系和卫星测量中常用的坐标系的建立方法。 L天球直角坐标系 厂天球坐标系 天球球面坐标系 地球直角坐标系地球大地坐标系 常用的天球坐标系：天球赤道坐标系、天球地平坐标系和天文坐标系。 在天球坐标系中，天体的空间位置可用天球空间直角坐标系或天球球面坐标系两种方式来描述。 1天球空间直角坐标系的定义 地球质心0为坐标原点，Z轴指向天球北极，X轴指向春分点，丫轴垂直于XOZ 平面，与X轴和Z轴构成右手坐标系。则在此坐标系下，空间点的位置由坐标（X，丫Z）来描述。 春分点：当太阳在地球的黄道上由天球南半球进入北半球，黄道与赤道的交点）

A ＜空闵直笥坐瑟厂K V ： z 丿的楚辽” 2天球球面坐标系的定义 地球质心0为坐标原点，春分点轴与天轴（天轴：地球自转的轴）所在平面为天 球经度（赤经）测量基准一一基准子午面，赤道为天球纬度测量基准而建立球面 坐标。空间点的位置在天球坐标系下的表述为（r ,a,S ）。 天欢申诗与地球质?M 重合T 赤礙刊为舍天黏 和感分点的天球子牛面 与过天体$的天球子牛面 之间的夾角，未纬 S 为 原点Mi 天体￡的连規与 天球击道面之间的夹角, 旬題丫为展点Mi 天体S 球球】?坐抚1就，S 1 r ）的C 义： 天球空间直角坐标系与天球球面坐标系的关系可用图 2-1表示: 感鼻—地I 球质心M 一孑塾一指向天球北奴Pn 、 ￥菇'一垂直于XMZ 平面, 与X 抽和Z 抽枸成右 手坐 标系统。 Pn A Z y X 1 \y X 奋 My\5 Ps / /

对同一空间点，直角坐标糸与其著效的球面坐标糸参教间有如下转换关务: C X - /cos a cos S < Y= / sin cos -Z = ysin 5 Y V a = arctan —— L Xz d -arctail . 岁差和章动的影响 岁差：地球实际上不是一个理想的球体，地球自转轴方向不再保持不变，这使春分点在黄道上产生缓慢的西移，这种现象在天文学中称为岁差。 章动：在日月引力等因素的影响下，瞬时北天极将绕瞬时平北天极旋转，大致呈椭圆，这种现象称为章动。 极移：地球自转轴相对地球体的位置并不是固定的，因而，地极点在地球表面上的位置，是随时间而变化的，这种现象称为极移。地球的自转轴不仅受日、月引力作用而使其在空间变化，而且还受地球内部质量不均匀影响在地球内部运动。 前者导致岁差和章动，后者导致极移。 协议天球坐标系：为了建立一个与惯性坐标系统相接近的坐标系，人们通常选择某一时刻，作为标准历元，并将此刻地球的瞬时自转轴（指向北极）和地心至瞬 时春分点的方向，经过瞬时的岁差和章动改正后，分别作为 X轴和Z轴的指向, 由此建立的坐标系称为协议天球坐标系。天味奋 5 y X X Ps

北京54坐标系转换工具

北京54坐标系转换工具 利用ARCGIS进行自定义坐标系和投影转换 ARCGIS种通过三参数和其参数进行精确投影转换 注意：投影转换成54坐标系需要下载无偏移卫星图像进行转换，有偏移的转换将导致转换后的卫星图像扭曲，坐标错误，无法配准。 第一步：选择无偏移地图源，下载你所需要的卫星图像。 第二步：选择BIGEMAP软件右边工具栏，选择【投影转换】，如下图所示： 2.1 选择说明： 1. 源文件：选择下载好的卫星图像文件（下载目录中后缀为tiff的文件） 2. 源坐标系：打开的源文件的投影坐标系（自动读取，不需要手动填写） 3. 输出文件：选择转换后你要保持文件的文件路径和文件名 4. 目标坐标系：选择你要转换成的目标坐标系，如下图：

选择上图的更多，如下图所示： 1：选择 -Beijing 1954 2：选择地区3：选择分度带对应的带号（一般默认，也可以手动修改）

选择对应的分度带或者中央子午线（请参看：如何选择分度带？），点击【确定】 5. 重采样算法：投影转换需要将影像的像素重新排列，一次每种算法的效率不一样，一般选择【立方卷积采样】，以达到最好的效果。如下图： 6. 指定变换参数：在不知道的情况下，可以不用填此处信息，如果√上，则如下图：

此参数为【三参数】或者【七参数】，均为国家保密参数，需要到当地的测绘部门或者国土部门，以单位名义签保密协议进行购买，此参数各地都不一样，是严格保密的，请不要随便流通。 第三步：点击【确定】，开始转换，如下图：

第四步：完成后，打开你刚才选择的输出文件夹，里面就是转换后的卫星图像。 第五步：如果你需要套合你手里已经有的矢量文件，请参看：【BIGEMAP无偏移影像叠加配准】

各种坐标系的关系

WGS84经纬度坐标与北京54坐标或者西安80坐标的关系一般来讲，GPS直接提供的坐标（B,L,H）是1984年世界大地坐标系(Word Geodetic System1984即WGS-84)的坐标，其中B为纬度，L为经度，H为大地高即是到WGS-84椭球面的高度。而在实际应用中，我国地图采用的是1954北京坐标系或者1980西安坐标系下的高斯投影坐标（x,y,），不过也有一些电子地图采用1954北京坐标系或者1980西安坐标系下的经纬度坐标（B,L），高程一般为海拔高度h。 GPS的测量结果与我国的54系或80系坐标相差几十米至一百多米，随区域不同，差别也不同，经粗落统计，我国西部相差70米左右，东北部140米左右，南部75米左右，中部45米左右。 1、1984世界大地坐标系 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。 2、1954北京坐标系 1954北京坐标系是将我国大地控制网与前苏联1942年普尔科沃大地坐标系相联结后建立的我国大地坐标系。属于参心大地坐标系，采用了前苏联的克拉索夫斯基椭球体。其长半轴 a=6378245，扁率

f=1/298.3。1954年北京坐标系虽然是苏联1942年坐标系的延伸,但不能说它们完全相同。 3、1980西安坐标系 1978年，我国决定建立新的国家大地坐标系统,并且在新的大地坐标系统中进行全国天文大地网的整体平差,这个坐标系统定名为1980年西安坐标系。属参心大地坐标系。1980年西安坐标系Xi'an Geodetic Coordinate System 1980 采用1975国际椭球,以JYD 1968.0系统为椭球定向基准,大地原点设在陕西省泾阳县永乐镇,采用多点定位所建立的大地坐标系.其椭球参数采用1975年国际大地测量与地球物理联合会推荐值,它们为:其长半轴a=6378140m; 扁率f=1/298.257。 4 高斯平面直角坐标系和UTM 一般的地图均为平面图,其对应的也是平面坐标.因此,需要将椭球面上各点的大地坐标,按照一定的数学规律投影到平面上成为平面直角坐标.目前世界各国采用最广泛的高斯- 克吕格投影和墨卡托投影(UTM)均是正形投影（等角投影），即该投影在小区域范围内使平面图形与椭球面上的图形保持相似。为了限制长度变形，根据国际测量协会规定,将全球按一定经差分成若干带。我国采用6度带或3度带，6度带是自零度子午线起每隔经度。 高斯平面直角坐标系一般以中央经线（L0）投影为纵轴X, 赤道投影为横轴Y，两轴交点即为各带的坐标原点。为了避免横坐标出现负值，在投影中规定将坐标纵轴西移500公里当作起始轴。为了区

大地测量坐标系统及其转换(精)

大地测量坐标系统及其转换 基本坐标系 1、大地坐标系 坐标表示形式:(, ,L B H 大地经度L :地面一点P 地的大地子午面N P S 与起始大地子午面所构成的二面角; 大地纬度B :P 地点对椭球面的法线P P K 地与赤道面所夹的锐角; 大地高 H :P 地点沿法线到椭球面的距离。 赤道面 S W 2、空间直角坐标系

坐标表示形式:(,,X Y Z 以椭球中心O 为坐标原点,起始子午面N G S 与赤道面的交线为X 轴,椭球的短轴为Z 轴(向北为正,在赤道面上与X 轴正交的方向为Y 轴,构成右手直角坐标系O X YZ 。 Y W 3、子午平面坐标系 坐标表示形式:(,, L x y 设P点的大地经度为L,在过P点的子午面上,以椭圆的中心为原点,建立x、y 平

面直角坐标系。则点P的位置用(,, L x y表示。 x 坐标表示形式:(,,L u H 设椭球面上的点P 的大地经度为L 。在此子午面,以椭球中心O 为圆心,以椭球长半径a 为半径,做一个辅助圆。过P 点做一纵轴的平行线,交横轴于1P 点,交辅助圆于2P 点,连结2P 、O 点,则21P O P 称为P 点的归化纬度,用u 来表示。P 点的位置用(,L u 表示。 当P 点不在椭球面上时,则应将P 沿法线投影到椭球面上,得到点0P ,0PP 即为P 点的大地高,0P 点的归化纬度,就是P 点的归化纬度。P 点的位置用(,,L u H 表示。

x y P u 点在椭球面上时的 P u 点不在椭球面上时的x

坐标表示形式:(,, L φρ 设P 点的大地经度为L ,连结O P ,则POx φ∠=,称为球心纬度,OP ρ=,称为P 点的向径。P 点的位置用(,,L φρ表示。 x 6、大地极坐标系 坐标表示形式:(,S A 以椭球面上某点0P 为极点,以0P 的子午线为极轴,从0P 出发,作一族A =常数的大地线和S =常数的大地圆。它们构成相互正交的坐标系曲线,即椭球面上的大地极坐标系,简称地极坐标系。在大地极坐标系中,点的位置用(,S A 来表示。 P A =常数 S =常数 坐标表示形式:1(,,P X Y Z -

坐标转换器使用说明

大地坐标（BLH） 平面直角坐标（XYZ） 四参数：X 平移、Y 平移、旋转角和比例 七参数：X平移，Y平移，Z 平移，X 轴旋转，Y 轴旋转，Z 轴旋转，缩放比例（尺度比） GPS控制网是由相对定位所求的的基线向量而构成的空间基线基线向量网，在GPS控制网的平差中，是以基线向量及协方差为基本观测量。 图3-1表示为HDS2003数据处理软件进行网平差的基本步骤，从图中可以看到，网平差实际上可以分为三个过程： l、前期的准备工作，这部分是用户进行的。即在网平差之前，需要进行坐标系的设置、并输入已知点的经纬度、平面坐标、高程等。 2、网平差的实际进行，这部分是软件自动完成的； 3、对处理结果的质量分析与控制，这部分也是需要用户分析处理的过程。 图3-1 平差过程 坐标系选择 针对不同的平差，要相应选择不同的坐标系，是否输入相应信息。在笔者接触过的项目中，平差时先通过三维无约束平差后，再进行二维约束平差。由于先进行的时三维无约束平差，是在WGS84坐标系统下进行的。 首先更改项目的坐标系统。在菜单“项目”->“坐标系统”或在工具栏“坐标系统”，则弹出“坐标

系统”对话框，选择WGS-84坐标。 图3-2 坐标系统 这里注意的是，在“投影”下见图，中央子午线是114°。很多情况下这里需要进行修改。 图3-3 WGS84投影 软件中自带的“中国-WGS 84”是允许修改的，我们换种方法：就是新建一个坐标文件，其他参数都和“中国-WGS84”一致，仅仅将中央子午线修改下。 在上图中，点击“新建”，得到“COORD GM”对话框，在“文件”->“新建”，如图

图3-4 新建坐标系统 然后在“设置”->“地图投影”，直接修改中央子午线，这里以81°为例，点击确定后，返回“COORD GM”对话框。 图3-5 投影设置 将输入源坐标和输入目标坐标的椭球，均改为WGS84。在“文件”->“保存”，输入名称和国家（中国），退出操作。

54坐标系、80坐标系、84坐标系之间的转换关系

工程施工过程中，常常会遇到不同坐标系统间，坐标转换的问题。目前国内常见的转换有以下几种：1，大地坐标（BLH）对平面直角坐标（XYZ）；2，北京54全国80及WGS84坐标系的相互转换；3，任意两空间坐标系的转换。其中第2类可归入第三类中。所谓坐标转换的过程就是转换参数的求解过程。常用的方法有三参数法、四参数法和七参数法。以下对上述三种情况作详细描述如下： 1，大地坐标（BLH）对平面直角坐标（XYZ） 常规的转换应先确定转换参数，即椭球参数、分带标准（3度，6度）和中央子午线的经度。椭球参数就是指平面直角坐标系采用什么样的椭球基准，对应有不同的长短轴及扁率。一般的工程中3度带应用较为广泛。对于中央子午线的确定有两种方法，一是取平面直角坐标系中Y坐标的前两位*3，即可得到对应的中央子午线的经度。如x=3250212m，y=395121123m，则中央子午线的经度=39*3=117度。另一种方法是根据大地坐标经度，如果经度是在155.5~185.5度之间，那么对应的中央子午线的经度=（155.5+185.5）/2=117度，其他情况可以据此3度类推。 另外一些工程采用自身特殊的分带标准，则对应的参数确定不在上述之列。 确定参数之后，可以用软件进行转换，以下提供坐标转换的程序下载。 2，北京54全国80及WGS84坐标系的相互转换 这三个坐标系统是当前国内较为常用的，它们均采用不同的椭球基准。 其中北京54坐标系，属三心坐标系，大地原点在苏联的普而科沃，长轴6378245m，短轴6356863，扁率1/298.3；西安80坐标系，属三心坐标系，大地原点在陕西省径阳县永乐镇，长轴6378140m，短轴6356755，扁率1/298.25722101；WGS84坐标系，长轴6378137.000m，短轴6356752.314，扁率1/298.257223563。由于采用的椭球基准不一样，并且由于投影的局限性，使的全国各地并不存在一至的转换参数。对于这种转换由于量较大，有条件的话，一般都采用GPS联测已知点，应用GPS软件自动完成坐标的转换。当然若条件不许可，且有足够的重合点，也可以进行人工解算。详细方法见第三类。 3，任意两空间坐标系的转换 由于测量坐标系和施工坐标系采用不同的标准，要进行精确转换，必须知道至少3个重合点（即为在两坐标系中坐标均为已知的点。采用布尔莎模型进行求解。布尔莎公式: 对该公式进行变换等价得到： 解算这七个参数，至少要用到三个已知点（2个坐标系统的坐标都知道），采用间接平差模型进行解算： 其中：V 为残差矩阵; X 为未知七参数; A 为系数矩阵; 解之：L 为闭合差 解得七参数后，利用布尔莎公式就可以进行未知点的坐标转换了，每输入一组坐标值，就能求出它在新坐标系中的坐标。但是要想GPS观测成果用于工程或者测绘，还需要将地方直

CORS坐标转换软件使用说明

坐标转换软件使用说明　 1、功能介绍　 在南京进行测量的同行一直受到坐标系统和已知控制点的困扰， 所以往往许多测量成果因坐标系统问题得不到承认，浪费了大量的人 力物力。基于此：本公司集全部精干技术力量，研发本款坐标转换软 件，可以说：它是全体测量工作者的福音。　 南京CORS因为其免费，应用十分广泛，但是使用南京CORS在 很多情况下，因为已知控制点原因无法实地取得平面坐标而限制了 CORS优势的发挥。本软件可以实现基于南京CORS测量的WGS84 坐标与92南京地方坐标双向自由转换，转换精度与权威部门转换成 果比较（在南京市6800平方公里范围内，包括高淳、溧水、六合、 浦口）：平面残差中误差优于±5mm、高程残差中误差均优于±1cm。精度完全具有保障，免去到处寻找控制点带来的人力、财力和时间浪费。按照最新城市规范规定，这种模式可以实现城市E级GPS控制 点的平面测量。　 本软件是一款后处理软件，即：内业处理软件，它不能在实地计 算坐标，通过事后（采集）或事前（放样）数据处理，同样可以让你 在野外无忧无障碍开展工作。　 适用平台：Windows 32位所有系统平台。　 2、外业采集数据转换操作介绍　 外业测量数据从RTK手簿中以WGS84坐标格式导出，导出以后 将文件复制到计算机，假设文件名为0513.dat。在电脑中启动软件，

界面如下：　 图一：程序启动界面　 首先选择转换方向下拉列表框，此时选择“WGS84—>NJ92”，表示将WGS84坐标转向92南京地方坐标，此时软件会出现一个按钮 键读入数据并转换，点击该按钮，在弹出的文件对话框中选择从手簿 导出的外业坐标文件。如：0513.dat，点击打开按钮即可完成转换。如图二：　 图二：选择原始数据文件　 记得一定要选择你的原始数据文件格式在点击打开按钮。转换完 成以后又会在对话框中再出现一个按钮导出转换成果，点击它即可将

坐标系之间的转换

大地坐标（BLH经纬度高程）和北京54等坐标系之间的转换 2008-12-11 16:25:23| 分类：默认分类| 标签：|字号大中小订阅 工程施工过程中，常常会遇到不同坐标系统间，坐标转换的问题。目前国内常见的转换有以下几种：1，大地坐标（BLH）对平面直角坐标（XYZ）；2，北京54全国80及WGS84坐标系的相互转换；3，任意两空间坐标系的转换。其中第2类可归入第三类中。所谓坐标转换的过程就是转换参数的求解过程。常用的方法有三参数法、四参数法和七参数法。以下对上述三种情况作详细描述如下： 1，大地坐标（BLH）对平面直角坐标（XYZ） 常规的转换应先确定转换参数，即椭球参数、分带标准（3度，6度）和中央子午线的经度。椭球参数就是指平面直角坐标系采用什么样的椭球基准，对应有不同的长短轴及扁率。一般的工程中3度带应用较为广泛。对于中央子午线的确定有两种方法，一是取平面直角坐标系中Y坐标的前两位*3，即可得到对应的中央子午线的经度。如x=3250212m，y=395121123m，则中央子午线的经度=39*3=117度。另一种方法是根据大地坐标经度，如果经度是在155.5~185.5度之间，那么对应的中央子午线的经度=（155.5+185.5）/2=117度，其他情况可以据此3度类推。 另外一些工程采用自身特殊的分带标准，则对应的参数确定不在上述之列。 确定参数之后，可以用软件进行转换，以下提供坐标转换的程序下载。 2，北京54全国80及WGS84坐标系的相互转换 这三个坐标系统是当前国内较为常用的，它们均采用不同的椭球基准。 其中北京54坐标系，属三心坐标系，大地原点在苏联的普而科沃，长轴6378245m，短轴6356863，扁率1/298.3；西安80坐标系，属三心坐标系，大地原点在陕西省径阳县永乐镇，长轴6378140m，短轴6356755，扁率1/298.25722101；WGS84坐标系，长轴6378137.000m，短轴6356752.314，扁率1/298.257223563。由于采用的椭球基准不一样，并且由于投影的局限性，使的全国各地并不存在一至的转换参数。对于这种转换由于量较大，有条件的话，一般都采用GPS联测已知点，应用GPS软件自动完成坐标的转换。当然若条件不许可，且有足够的重合点，也可以进行人工解算。详细方法见第三类。 3，任意两空间坐标系的转换 由于测量坐标系和施工坐标系采用不同的标准，要进行精确转换，必须知道至少3个重合点（即为在两坐标系中坐标均为已知的点。采用布尔莎模型进行求解。布尔莎公式: 对该公式进行变换等价得到： 解算这七个参数，至少要用到三个已知点（2个坐标系统的坐标都知道），采用间接平差模型进行解算： 其中：V 为残差矩阵; X 为未知七参数; A 为系数矩阵; 解之：L 为闭合差

施工坐标系与测量坐标系之间的相互转换关系

施工坐标系与测量坐标系之间的相互转换一、用Microsoft Excel 编辑转换

如图（1-1）所示：设Y O X -- 为测量坐标系，y o x -'- 为施工坐标，如果知道了施工坐标系的原点o '的测量坐标为（'0X ,'0Y ）、定向点I 的测量坐标为（XI,YI ），定向坐标方位角 x -'0α （即纵轴的旋转角，因为0=-X o α为正北方向，则x -'0α=X o -α+α）。则所求P 点由施工坐标P （p p y x ,）换算成为测量坐标P （p p Y X ,）的公式则为： ααsin *cos *0p p p y x X X -+=' ααcos *sin *0p p p y x Y Y ++=' 上面两式在Excel 中编辑公式为： [][]180/()*sin *180/()*cos *0Pi y Pi x X X p p p αα-+=' [][]180/()*cos *180/()*sin *0Pi y Pi x Y Y p p p αα++=' 而如果知道了施工坐标系（第二坐标系）的原点的测量坐标 o '为（'0X 、'0Y ）、坐标方位角 x -'0α （即纵轴的旋转角，因为0=-X o α为正北方向，则x -'0α=X o -α+α）。则所求P 点由测量坐标P （p p Y X ,）转换算为施工坐标P （p p y x ,）其公式为： ααsin *)(cos *)(00''-+-=Y Y X X x p p p ααcos *)(sin *)(00''-+--=Y Y X X y p p p 上面两式在Excel 中编辑公式为：

参考坐标与动坐标系之间的旋转变换

坐标系之间的坐标变换 取一参考坐标系Z Y X O ''''，该坐标系为船舶平衡位置上，不随船舶摇荡。 取一动坐标系OXYZ ，该坐标系与船体固结，随船舶一起做摇荡运动，OX 轴位于纵中剖面内，并指向船首，OY 垂直向上，OZ 轴指向船舶右舷。 再取一坐标系Z Y X O ???，它与参考坐标系平行，原点与动坐标系重合，且仅随船体作振荡运动。这三个坐标系之间的相对位置如图所示： 角位移用欧拉角来定义。我们假设动坐标系OXYZ 的原始位置为Z Y X O ???，经三次转动转到目前的位置。 首先将坐标系Z Y X O ???绕X O ?轴转动α角，使其转到OZ 和X O ?所确定的平面，然后绕Y O ?轴旋转β角使Z O ?与OZ 重合，此时平面Y X O ''??和平面OXY 重合，最后将得到的Z Y X O ''??绕OZ 轴转动γ角，这样，坐标系OXYZ 和坐标系Z Y X O ???就完全重合。 第一次旋转可以写为： ααααcos ?sin ??sin ?cos ????Z Y Z Z Y Y X X '+'='-'== 写为矩阵形式为 ????? ? ??''????? ??-=?????? ??Z Y X Z Y X ???cos sin 0sin cos 000 1???αα αα

同理，第二次旋转得 ?????? ??''????? ??-=?????? ??''Z Y X Z Y X ??cos 0sin 010sin 0cos ???ββ ββ 第三次旋转得， ???? ? ??????? ??-=?????? ??''Z Y X Z Y X 10 0cos sin 0sin cos ??γγγ γ 综合上面三式，得 ???? ? ????? ? ? ??++--+-+-=?????? ??Z Y X Z Y X βαγ αγβαγ αγβαβαγαγβαγαγβαβγ βγβcos cos cos sin sin sin cos sin sin cos sin cos cos sin cos cos sin sin sin sin cos cos sin sin sin sin cos cos cos ???则 [][][]X r X O '+='