空间大地坐标系与平面直角坐标系转换公式详细版.doc
§2.3.1 坐标系的分类
正如前面所提及的,所谓坐标系指的是描述空间位置的表达形式,即采用什么方法来表示空间位置。人们为了描述空间位置，采用了多种方法，从而也产生了不同的坐标系，如直角坐标系、极坐标系等。
在测量中常用的坐标系有以下几种：
一、空间直角坐标系
空间直角坐标系的坐标系原点位于参考椭球的中心，Z 轴指向参考椭球的北极，X 轴指向起始子午面与赤道的交点，Y 轴位于赤道面上且按右手系与X 轴呈90°夹角。某点在空间中的坐标可用该点在此坐标系的各个坐标轴上的投影来表示。空间直角坐标系可用图2-3来表示：
图2-3 空间直角坐标系
二、空间大地坐标系
空间大地坐标系是采用大地经、纬度和大地高来描述空间位置的。纬度是空间的点与参考椭球面的法线与赤道面的夹角；经度是空间中的点与参考椭球的自转轴所在的面与参考椭球的起始子午面的夹角；大地高是空间点沿参考椭球的法线方向到参考椭球面的距离。空间大地坐标系可用图2-4来表示：
图2-4空间大地坐标系
三、平面直角坐标系
平面直角坐标系是利用投影变换，将空间坐标空间直角坐标或空间大地坐标通过某种数学变换映射到平面上，这种变换又称为投影变换。投影变换的方法有很多，如横轴墨卡托投影、UTM 投影、兰勃特投影等。在我国采用的是高斯-克吕格投影也称为高斯投影。UTM 投影和高斯投影都是横轴墨卡托投影的特例，只是投影的个别参数不同而已。
高斯投影是一种横轴、椭圆柱面、等角投影。从几何意义上讲，是一种横轴椭圆柱正切投影。如图左侧所示，设想有一个椭圆柱面横套在椭球外面，并与某一子午线相切（此子午线称为中央子午线或轴子午线），椭球轴的中心轴CC ’通过椭球中心而与地轴垂直。
高斯投影满足以下两个条件：
1、 它是正形投影；
2、 中央子午线投影后应为x 轴，且长度保持不变。
将中央子午线东西各一定经差（一般为6度或3度）范围内的地区投影到椭圆柱面上，再将此柱面沿某一棱线展开，便构成了高斯平面直角坐标系，如下图２－５右侧所示。
图2-5 高斯投影
x 方向指北，y 方向指东。
可见，高斯投影存在长度变形，为使其在测图和用图时影响很小，应相隔一定的地区，另立中央子午线，采取分带投影的办法。我国国家测量规定采用六度带和三度带两种分带方法。六度带和三度带与中央子午线存在如下关系：
366 N L ＝中； n L 33＝中
其中，N 、n 分别为6度带和3度带的带号。
另外，为了避免y 出现负号，规定y 值认为地加上500000m ；又为了区别不同投影带，前面还要冠以带号，如第20号六度带中，y=-200.25m ，则成果表中写为y 假定＝20499799.75m 。x 值在北半球总显正值，就无需改变其观测值了。
1、空间直角坐标系与空间大地坐标系间的转换
图２－６表示了空间直角坐标系与空间大地坐标系之间的关系。
图2-6 地球空间直角坐标系与大地坐标系
在相同的基准下空间大地坐标系向空间直角坐标系的转换公式为：
??
???+-=+=+=B H e N Z L B H N Y L B H N X sin ])1([sin cos )(cos cos )(2 (2-1)
式中， W
a N =，a 为椭球的长半轴，N 为椭球的卯酉圈曲率半径 a =6378.137km
B e W 22sin 1-=
22
22
a b a e -=，e 为椭球的第一偏心率，b 为椭球的短半轴 b =6356.7523141km
在相同的基准下空间直角坐标系向空间大地坐标系的转换公式为
?????
??????-Φ=??? ??=?????????? ??+Φ=N B R H X Y arctg L W B Z ae tg arctg B cos cos sin 12 (2-2) 式中
??????+=Φ22Y
X Z arctg 222Z Y X R ++=
2、空间坐标系与平面直角坐标系间的转换
空间坐标系与平面直角坐标系间的转换采用的是投影变换的方法。在我国一般采用的是高斯投影。因为高斯投影和UTM 投影都是横轴墨卡托的特例，因此，高斯投影和UTM 投影都可以套用横轴墨卡托投影的投影公式。
横轴墨卡托投影的投影的正反算公式可参见有关资料，它们的区别在于轴子午线投影到平面上后，其长度的系数，对于高斯投影，系数为1，对于UTM 投影，其系数为0.9996。 3、变动高程归化面的影响
用户在建立地方独立坐标系时，有时变动高程归化面，这将产生一个新椭球，这就必须计算新常数，新椭球常数按下列方法和步骤进行：
1) 新椭球是在国家坐标系的参考椭球上扩大形成的，它的扁率应与国家坐标系参考椭球的扁率相等，即a a ='。
2) 计算该坐标系中央地区的新椭球平均曲率半径和新椭球长半轴。
新椭球平均曲率半径为：
m m
m m m m H B e e a H W a W e a H MN H R R +--=+-=+=+=222
32sin 11)1(' (2.10) 式中
m H ───该地区平均大地高；
m B ───该地区的平均纬度。
新椭球的长半轴按下式计算：
2221sin 1''e B e R a m
--= (2.11)
将新的椭球参数代入，就可以进行投影的正反计算了。
二、坐标系统的转换方法
不同坐标系统的转换本质上是不同基准间的转换，不同基准间的转换方法有很多，其中最为常用的有布尔沙模型，又称为七参数转换法。
七参数转换法是：
设两空间直角坐标系间有七个转换参数：3 个平移参数()z y x ???、3 个旋转参数()z y x εεε和1 个尺度参数k 。比如，由空间直角坐标系A 转换到空间直角坐标系B 可
采用下面的公式： A
x y x z y z A B Z Y X Z Y X k z y x Z Y X ????? ??????? ??--+????? ??++????? ?????=????? ??000)1(εεεεεε §2.3.4 GPS 测量中常用的坐标系统
一、世界大地坐标系WGS-84
WGS-84 坐标系是目前GPS 所采用的坐标系统，GPS 所发布的星历参数和历书参数等都是基于此坐标系统的。
WGS-84 坐标系统的全称是World Geodical System-84 (世界大地坐标系-84), 它是一个地心地固坐标系统。WGS-84 坐标系统由美国国防部制图局建立，于1987 年取代了当时GPS 所采用的坐标系统WGS-72 坐标系统而成为现在GPS 所使用的坐标系统。
WGS-84 坐标系的坐标原点位于地球的质心，Z 轴指向BIH 1984.0 定义的协议地球极方向，X 轴指向BIH 1984.0 的启始子午面和赤道的交点，Y 轴与X 轴和Z 轴构成右手系。
WGS-84 系所采用椭球参数为见表2.1。
二、1954 年北京坐标系
1954 年北京坐标系是我国目前广泛采用的大地测量坐标系。该坐标系源自于原苏联采用过的1942 年普尔科夫坐标系。该坐标系采用的参考椭球是克拉索夫斯基椭球。该椭球的参数见表2.1。
遗憾的是该椭球并未依据当时我国的天文观测资料进行重新定位，而是由前苏联西伯利亚地区的一等锁经我国的东北地区传算过来的，该坐标系的高程异常是以前苏联1955 年大地水准面重新平差的结果为起算值，按我国天文水准路线推算出来的，而高程又是以1956 年青岛验潮站的黄海平均海水面为基准。
由于当时条件的限制1954 年北京坐标系存在着很多缺点主要表现在以下几个方面：
1. 克拉索夫斯基椭球参数同现代精确的椭球参数的差异较大，并且不包含表示地球物理特性的参数，因而给理论和实际工作带来了许多不便。
2. 椭球定向不十分明确，椭球的短半轴既不指向国际通用的CIO 极，也不指向目前我国使用的JYD极。参考椭球面与我国大地水准面呈西高东低的系统性倾斜，东部高程异常达60余米，最大达67 米。
3. 该坐标系统的大地点坐标是经过局部分区平差得到的。因此全国的天文大地控制点实际上不能形成一个整体，区与区之间有较大的隙距，如在有的接合部中同一点在不同区的坐标值相差1-2 米，不同分区的尺度差异也很大，而且坐标传递是从东北到西北和西南，后一区是以前一区的最弱部作为坐标起算点，因而一等锁具有明显的坐标积累误差。
三、1980 年西安大地坐标系
1978 年我国决定重新对全国天文大地网施行整体平差，并且建立新的国家大地坐标系统。整体平差在新大地坐标系统中进行，这个坐标系统就是1980 年西安大地坐标系统。1980 年西安大地坐标系统所采用的地球椭球参数的四个几何和物理参数采用了IAG 1975 年的推荐值，见表2.1中的西安80。
椭球的短轴平行于地球的自转轴（由地球质心指向1968.0 JYD 地极原点方向），起始子午面平行于格林尼治平均天文子午面，椭球面同似大地水准面在我国境内符合最好，高程系统以1956 年黄海平均海水面为高程起算基准。
四、几种常用的坐标系统的几何和物理参数
下表列出了几种常用的坐标系统的几何和物理参数，用户需要时可以查阅：
表 2.1 GPS 测量中常用的坐标系统的几何和物理参数
§2.4 GPS 高程系统
在测量中常用的高程系统有大地高系统、正高系统和正常高系统。
§2.4.1 大地高系统
大地高系统是以参考椭球面为基准面的高程系统，某点的大地高是该点到通过该点的参考椭球的法线与参考椭球面的交点间的距离。大地高也称为椭球高。大地高一般用符号H 表示。
大地高是一个纯几何量，不具有物理意义，同一个点在不同的基准下具有不同的大地高。 通常，GPS 接收机单点定位得到的高程为WGS-84下的大地高。
§2.4.2 正高系统
正高系统是以大地水准面为基准面的高程系统，某点的正高是该点到通过该点的铅垂线与大地水准面的交点之间的距离。正高用符号 H g 表示。
§2.4.3 正常高
正常高系统是以似大地水准面为基准的高程系统，某点的正常高是该点到通过该点的铅垂线与似大地水准面的交点之间的距离，正常高用 H γ 表示。
§2.4.4 高程系统之间的转换关系
大地水准面到参考椭球面的距离称为大地水准面差距，记为 h g ，大地高与正高之间的关系可以表示为：
正 高：g g h H H -=
似大地水准面到参考椭球面的距离，称为高程异常，记为ζ。大地高与正常高之间的关系可以表示为：
正常高：ζγ-=H H
高程之间的相互关系可以用下图2-7来表示：
图2-7 高程系统间的相互关系
空间大地坐标系与平面直角坐标系转换公式
§2.3.1 坐标系的分类 正如前面所提及的,所谓坐标系指的是描述空间位置的表达形式,即采用什么方法来表示空间位置。人们为了描述空间位置，采用了多种方法，从而也产生了不同的坐标系，如直角坐标系、极坐标系等。 在测量中常用的坐标系有以下几种： 一、空间直角坐标系 空间直角坐标系的坐标系原点位于参考椭球的中心，Z 轴指向参考椭球的北极，X 轴指向起始子午面与赤道的交点，Y 轴位于赤道面上且按右手系与X 轴呈90°夹角。某点在空间中的坐标可用该点在此坐标系的各个坐标轴上的投影来表示。空间直角坐标系可用图2-3来表示： 图2-3 空间直角坐标系 二、空间大地坐标系 空间大地坐标系是采用大地经、纬度和大地高来描述空间位置的。纬度是空间的点与参考椭球面的法线与赤道面的夹角；经度是空间中的点与参考椭球的自转轴所在的面与参考椭球的起始子午面的夹角；大地高是空间点沿参考椭球的法线方向到参考椭球面的距离。空间大地坐标系可用图2-4来表示：
图2-4空间大地坐标系 三、平面直角坐标系 平面直角坐标系是利用投影变换，将空间坐标空间直角坐标或空间大地坐标通过某种数学变换映射到平面上，这种变换又称为投影变换。投影变换的方法有很多，如横轴墨卡托投影、UTM 投影、兰勃特投影等。在我国采用的是高斯-克吕格投影也称为高斯投影。UTM 投影和高斯投影都是横轴墨卡托投影的特例，只是投影的个别参数不同而已。 高斯投影是一种横轴、椭圆柱面、等角投影。从几何意义上讲，是一种横轴椭圆柱正切投影。如图左侧所示，设想有一个椭圆柱面横套在椭球外面，并与某一子午线相切（此子午线称为中央子午线或轴子午线），椭球轴的中心轴CC ’通过椭球中心而与地轴垂直。 高斯投影满足以下两个条件： 1、 它是正形投影； 2、 中央子午线投影后应为x 轴，且长度保持不变。 将中央子午线东西各一定经差（一般为6度或3度）范围内的地区投影到椭圆柱面上，再将此柱面沿某一棱线展开，便构成了高斯平面直角坐标系，如下图２－５右侧所示。 图2-5 高斯投影 x 方向指北，y 方向指东。 可见，高斯投影存在长度变形，为使其在测图和用图时影响很小，应相隔一定的地区，另立中央子午线，采取分带投影的办法。我国国家测量规定采用六度带和三度带两种分带方法。六度带和三度带与中央子午线存在如下关系： 366 N L ＝中； n L 33＝中 其中，N 、n 分别为6度带和3度带的带号。
不同空间直角坐标系的转换
不同空间直角坐标系的转换 欧勒角 不同空间直角坐标系的转换，包括三个坐标轴的平移和坐标轴的旋转，以及两个坐标系的尺度比参数，坐标轴之间的三个旋转角叫欧勒角。 三参数法 三参数坐标转换公式是在假设两坐标系间各坐标轴相互平行，轴系间不存在欧勒角的条件下得出的。实际应用中，因为欧勒角不大，可以用三参数公式近似地进行空间直角坐标系统的转换。公共点只有一个时,采用三参数公式进行转换。
七参数法 用七参数进行空间直角坐标转换有布尔莎公式，莫洛琴斯基公式和范氏公式等。下面给出布尔莎七参数公式： 坐标转换多项式回归模型 坐标转换七参数公式属于相似变换模型。大地控制网中的系统误差一般呈区域性，当区域较小时，区域性的系统误差被相似变换参数拟合，故局部区域的坐标转换采用七参数公式模型是比较适宜的。但对全国或一个省区范围内的坐标转换，可以采用多项式回归模型，将各区域的系统偏差拟合到回归参数中，从而提高坐标转换精度。 两种不同空间直角坐标系转换时，坐标转换的精度取决于坐标转换的数学模型和求解转换系数的公共点坐标精度，此外，还与公共点的分布有关。鉴于地面控制网系统误差在???? ??????+??????????=??????????000111222Z Y X Z Y X Z Y X ???? ??????+????????????????????---+??????????+=??????????000111111222000)1(Z Y X Z Y X Z Y X m Z Y X X Y X Z Y Z εεεεεε
不同区域并非是一个常数，所以采用分区进行坐标转换能更好地反映实际情况，提高坐标转换的精度。
空间直角坐标系与大地坐标系转换程序
空间直角坐标系与大地坐标系转换程序 #include
坐标转换之计算公式
坐标转换之计算公式 一、参心大地坐标与参心空间直角坐标转换 1名词解释： A ：参心空间直角坐标系： a) 以参心0为坐标原点； b) Z 轴与参考椭球的短轴（旋转轴）相重合； c) X 轴与起始子午面和赤道的交线重合； d) Y 轴在赤道面上与X 轴垂直，构成右手直角坐标系0-XYZ ； e) 地面点P 的点位用（X ，Y ，Z ）表示； B ：参心大地坐标系： a) 以参考椭球的中心为坐标原点，椭球的短轴与参考椭球旋转轴重合； b) 大地纬度B ：以过地面点的椭球法线与椭球赤道面的夹角为大地纬度B ； c) 大地经度L ：以过地面点的椭球子午面与起始子午面之间的夹角为大地经度L ； d) 大地高H ：地面点沿椭球法线至椭球面的距离为大地高H ； e) 地面点的点位用（B ，L ，H ）表示。 2 参心大地坐标转换为参心空间直角坐标： ?? ???+-=+=+=B H e N Z L B H N Y L B H N X sin *])1(*[sin *cos *)(cos *cos *)(2 公式中，N 为椭球面卯酉圈的曲率半径，e 为椭球的第一偏心率，a 、b 椭球的长短半 径，f 椭球扁率，W 为第一辅助系数 a b a e 2 2-= 或 f f e 1*2-= W a N B W e =-=22sin *1( 3 参心空间直角坐标转换参心大地坐标
[]N B Y X H H e N Y X H N Z B X Y L -+=+-++==cos ))1(**)()(*arctan( )arctan(2 2222 二 高斯投影及高斯直角坐标系 1、高斯投影概述 高斯-克吕格投影的条件：1. 是正形投影；2. 中央子午线不变形 高斯投影的性质：1. 投影后角度不变；2. 长度比与点位有关，与方向无关； 3. 离中央子午线越远变形越大 为控制投影后的长度变形，采用分带投影的方法。常用3度带或6度带分带，城市或工 程控制网坐标可采用不按3度带中央子午线的任意带。 2、高斯投影正算公式： 5 2224253 2236 4254 42232)5814185(cos 120 )1(cos 6 cos )5861(cos sin 720 495(cos sin 24 cos sin 2l t t t B N l t B N Bl N y l t t B B N l t B B N Bl B N X x ηηηηη-++-++-+=+-+++-++=） 3、高斯投影反算公式：
空间坐标转换说明
空间坐标转换说明 集团文件版本号：（M928-T898-M248-WU2669-I2896-DQ586-M1988）
坐标转换说明 GPS 接收机接收到GPS （大地坐标：经度、纬度和高度值）信号后，并不利于显示，需要将大地坐标进行转换，现选用东北天坐标系（也叫站心坐标系）作为显示的依据。 GPS 接收机接收到的第一个信号L （经度）、B （纬度）和H （高度），作为东北天坐标系的原点。当接收到第二个信号时L 1、B 1和H 1，应用坐标转换公式，转换到东北天坐标系下进行显示。依次类推，凡是接收到的GPS 信号都转换到东北天坐标系下进行显示，在东北天坐标系下预测出来的坐标值通过坐标转换公式在显示屏上显示大地坐标（经度、纬度和高度）。 1.大地坐标与直角坐标的相互转化 对空间某一点，大地坐标系(L ,B ,H )到直角坐标系(X ,Y ,Z )的转换关系如下： ?? ? ?? +-=+=+=B H e N Z L B H N Y L B H N X sin ])1([sin cos )(cos cos )(2（1） 由直角坐标系(X ,Y ,Z )转化到大地坐标系(L ,B ,H )的公式如下： ??? ? ??? --=+-++==)1(sin /]})1((/[)(arctan{) /arctan(2222e N B Z H H e N Y X H N Z B X Y L （2） 式中：B e a N 22sin 1/-=，N 为该点的卯酉圈曲率半径；2222/)(a b a e -=，a 、 b 、e 分别为该大地坐标系对应参考椭球的长半轴、短半轴和第一偏心率。长半轴 a =6378137±2m ，短半轴 b =6356.7523142km ，90130066943799.02=e 。 从公式（2）看出，经度比较容易求得，纬度和高度必须通过迭代计算获直接计算得到。迭代计算的次序为：N H B →→，通常迭代四次可以达到H 优于0.001m ，B 优于0.00001''的计算精度；教科书中给出的直接法计算公式比较繁琐，有的计算公式的应用条件受到一定限制，例如要求大地高度小于10000m 时，才能使B 、H 达到上述计算精度，有的直接计算公式精度较低。 根据[张华海]提供的方法，本文建议采用该方法将直角坐标（X ,Y ,Z ）转变成大地坐标（L ,B ,H ）。该方法的公式形式比较简便，B 、H 的计算精度高；用计算出
大地坐标转换成施工坐标公式
大地(高斯平面)坐标系工程坐标系转换大地坐标系--->工程坐标系 ======================== 待转换点为P，大地坐标为：Xp、Yp 工程坐标系原点o: 大地坐标：Xo、Yo 工程坐标：xo、yo 工程坐标系x轴之大地方位角：a dX=Xp-Xo dY=Yp-Yo P点转换后之工程坐标为xp、yp: xp=dX*COS(a)+dY*SIN(a)+xo yp=-dX*SIN(a)+dY*COS(a)+yo 工程坐标系--->大地坐标系 ======================== 待转换点为P，工程坐标为：xp、yp 工程坐标系原点o: 大地坐标：Xo、Yo 工程坐标：xo、yo 工程坐标系x轴之大地方位角：a dx=xp-xo dy=yp-yo P点转换后之工程坐标为xp、yp: xp=Xo+dx*COS(a)-dy*SIN(a)
yp=Yo+dx*SIN(a)+dy*COS(a) 坐标方位角计算程序 置镜点坐标：ZX ZY 后视点坐标：HX HY 方位角：W 两点间距离: S Lb1 0← ｛A, B, C, D｝← A〝ZX=〞:B〝ZY=〞:C〝HX=〞:D 〝HY=〞:W=tg1((D-B)÷(C-A)):(D-B)>0=>(C-A)>0=>W=W:∟∟(D-B)>0=>(C-A)<0=>W=W+180:∟∟(D-B)<0=>(C-A)<0=>W=W+180:∟∟(D-B)<0=>(C-A)>0=>W=360+W∟∟W=W◢ S=√((D-B)2+(C-A)2) ◢ Goto 0← CASIO fx－4500p坐标计算程序 根据坐标计算方位角 W＝W＋360△W：“ALF(1~2)＝”L1 A“X1＝”：B“Y1＝”：Pol(C“X2”－A，D“Y2”－B：“S＝”▲W＜0 直线段坐标计算 L1 X“X(0)”：Y“Y(0)”：S“S(0)”：A“ALF” L2 Lb1 2 L3 {L}：L“LX”
直角坐标系下的画图及其转换公式
直角坐标系下的画图及其转换公式 在直角坐标系下我们的圆方程是： 222()()x a y b R -+-= 其中，a 和b 是圆心，R 是半径。但在画圆的时候，你就会发现如果按该公式画圆，多半是不成功的，或者画了一半，所以在matlab 中画圆，一半采用极坐标形式 圆对应的极坐标转换公式为： cos sin x R y R θ θ =?? =?（公式1） 这个很容易理解，你画个单位圆来看看就知道了。 那么上面那个黑色的点的x 坐标和y 坐标用半径和连线与坐标轴x 的夹角来表示，就得到了公式1。 观察这个公式，我们发现，在极坐标系下，圆的半径没变，夹角是在不断变化的，所以，在matlab 中极坐标系下画单位圆的问题可以这样来考虑： 首先将夹角360等分，也就是每一个步长为360度/360; 但需要指出的是，matlab 中正弦预先函数的变量其实是弧度，并不是度。这个你在matlab 命令窗里就可以试： 比如你要得到30度的正弦值，一般是sin （pi/6）,而不是sin(30)。这里的pi 是3.1415926的在matlab 中的表示。 所以我们的步长应该是弧度制的，我们知道，1度对应的弧度为360/(2*pi)。也即180/pi; 所以我们的夹角应该是： Theta=0:180/pi:2*pi-180/pi; 注意，由于是从零开始画图的，所以最后一个应该是2*pi-180/pi;而不是2*pi ； 这个时候我们可以开始画图了 X=R*cos(Theta); Y=R*sin(Theta); Plot(x,y,’r.’) axis square %保证画出来的圆是圆的。
空间直角坐标系坐标转换方法
坐标转换方法 空间直角坐标系如果其原点不动，绕着某一个轴旋转而构成的新的坐标系，这个过程就叫做坐标旋转。在旧坐标系中的坐标与在旋转后新坐标系中的坐标有一定的转换关系，这种转换关系可以用转换矩阵来表示。 如图5.7，直角坐标系XYZ，P点的坐标为（x, y, z），其相应的在XY 平面，XZ平面，YZ平面分别为M（x, y,0），Q（x,0, z）和N（0, y, z）。 图5.7直角坐标系XYZ 设?表示第j 轴的旋转角度，R j (?) 表示绕第j 轴的旋转，其正方向是沿坐标轴向原点看去的逆时针方向。很明显当j 轴为旋转轴时，它对应的坐标中的j 分量是不变的。由于直角坐标系是对称的，下面我们以绕Z轴旋转为例推导其旋转变换矩阵，其它两个轴推导和它是一样的。 设图5.7的坐标绕Z轴逆时针旋转θ角度，新坐标为X 'Y'Z'，如图5.8所示： 图5.8 坐标绕Z 轴逆时针旋转θ角度 由于坐标中的z 分量不变，我们可以简化地在XY 平面进行分分析，如图
5.9所示： 图5.9坐标绕Z 轴逆时针旋转θ 角度的XY 平面示意图 点 M X 和点M X ' 分别是M 点在X 轴和X '轴的投影。如图5.9 cos cos() sin sin() X X X X x OM OM MOM OM y MM OM MOM OM ?θ?θ==∠=-??==∠=-? （5-1） cos cos sin sin X X X X x OM OM MOM OM y MM OM MOM OM ? ?'''''==∠=??'==∠=? （5-2） 把（5-1）式按照三角函数展开得： cos cos sin sin sin cos cos sin x OM OM y OM OM ?θ?θ ?θ?θ=+??=+? （5-3） 把（5-2）式代入（5-3）式得： cos sin sin cos x x y y x y θθ θθ''=+??''=-+? （5-4） 坐标中的z 分量不变，即z = z'这样整个三维坐标变换就可以写成（用新坐标表 示旧坐标） cos sin sin cos x x y y x y z z θθ θθ''=+? ?''=-+??' =? （5-5） 把式（5-5）用一个坐标旋转变换矩阵R Z (θ) 表示可以写成：
大地坐标与直角空间坐标转换计算公式
大地坐标与直角空间坐标转换计算公式 一、参心大地坐标与参心空间直角坐标转换 1名词解释： A ：参心空间直角坐标系： a) 以参心0为坐标原点； b) Z 轴与参考椭球的短轴（旋转轴）相重合； c) X 轴与起始子午面和赤道的交线重合； d) Y 轴在赤道面上与X 轴垂直，构成右手直角坐标系0-XYZ ； e) 地面点P 的点位用（X ，Y ，Z ）表示； B ：参心大地坐标系： a) 以参考椭球的中心为坐标原点，椭球的短轴与参考椭球旋转轴重合； b) 大地纬度B ：以过地面点的椭球法线与椭球赤道面的夹角为大地纬度B ； c) 大地经度L ：以过地面点的椭球子午面与起始子午面之间的夹角为大地经度L ； d) 大地高H ：地面点沿椭球法线至椭球面的距离为大地高H ； e) 地面点的点位用（B ，L ，H ）表示。 2 参心大地坐标转换为参心空间直角坐标： ?? ? ?? +-=+=+=B H e N Z L B H N Y L B H N X sin *])1(*[sin *cos *)(cos *cos *)(2 公式中，N 为椭球面卯酉圈的曲率半径，e 为椭球的第一偏心率，a 、b 椭球的长短半径，f 椭球扁率，W 为第一辅助系数 a b a e 2 2-= 或 f f e 1 *2-= W a N B W e = -=22 sin *1( 西安80椭球参数： 长半轴a=6378140±5（m ）
短半轴b=6356755.2882m 扁 率α=1/298.257 3 参心空间直角坐标转换参心大地坐标 [ ] N B Y X H H e N Y X H N Z B X Y L -+= +-++==cos ))1(**)() (*arctan() arctan(2 22 2 2 二 高斯投影及高斯直角坐标系 1、高斯投影概述 高斯-克吕格投影的条件：1. 是正形投影；2. 中央子午线不变形 高斯投影的性质：1. 投影后角度不变；2. 长度比与点位有关，与方向无关； 3. 离中央子午线越远变形越大 为控制投影后的长度变形，采用分带投影的方法。常用3度带或6度带分带，城市或工程控制网坐标可采用不按3度带中央子午线的任意带。 2、高斯投影正算公式： 52224253 2236 425442232)5814185(cos 120 )1(cos 6 cos )5861(cos sin 720 495(cos sin 24cos sin 2l t t t B N l t B N Bl N y l t t B B N l t B B N Bl B N X x ηηηηη-++-++-+=+-+++-++ =） 3、高斯投影反算公式：
平面直角坐标变换
§5.7 平面直角坐标变换 为了考虑同一图形在不同的坐标系下的方程之间的关系，我们首先需要建立同一个点在不同的坐标系下的坐标之间的关系，这就是坐标变换的问题，因为我们研究的图形是点的轨迹． 我们仅考虑平面直角坐标变换． 设在平面上给出了由两个标架 {O ；i , j } 和 {O'；i', j' } 所决定的右手直角坐标系，这里i 和j 以及i' 和j' 是两组坐标基向量，它们是平面上的两个标准正交基，我们依次称这两个坐标系为旧坐标系和新坐标系． 由于坐标系的位置完全由原点和坐标基向量所决定，所以新坐标系与旧坐标系之间的关系，就由O' 在 {O ；i , j } 中的坐标以及i' 和j' 在 {O ；i , j } 中的分量所决定． 任一直角坐标变换总可以分解成移轴（也叫坐标平移）和转轴（也叫坐标旋转）两个步骤． 1．移轴 如果两个标架 {O ；i , j } 和 {O'；i , j' } 的原点O 与O' 不同，O' 在{O ；i , j }中的坐标为 (x 0，y 0)，但两标架的坐标基向量相同，即有 i' = i ， j' = j 那么标架 {O'；i', j'} 可以看成是由标架 {O ；i , j } 将原点平移到O'点而得来的（图5.7.1）．这种坐标变换叫做移轴（坐标平移）． 设P 是平面内任意一点，它对标架 {O ；i , j } 和 {O'；i', j'} 的坐标分别为 (x ，y ) 与 (y x '',)，则有 P O O O OP '+= 但 j i y x +=， j i y x O '+'='， j i 00y x O +=' 于是有 j i j i )()(00y y x x y x +'++'=+ 故 {x ，y } = {x 0，y 0} ＋ {x'，y' } 根据向量相等的定义得移轴公式为 图5.7.1 ? ? ?+'=+'=00 y y y x x x (5.7－1) 从中解出x' 和y'，就得逆变换公式为 ? ? ?-='-='00 y y y x x x (5.7－2) 2．转轴 若两个标架 {O ；i , j } 和 {O'；i', j'} 的原点相同，即O = O'，但坐标基向量不同，且有∠(i ，i' ) = α，则标架 {O'；i'，j'} 可以看成是由标架 {O ；i ，j } 绕O 点旋转α 角而得
北京54坐标转换为地理坐标的简易方法
北京54坐标转换为地理坐标的简易方法 1. 椭球体、基准面及地图投影 GIS中的坐标系定义是GIS系统的基础，正确定义GIS系统的坐标系非常重要。GIS中的坐标系定义由基准面和地图投影两组参数确定，而基准面的定义则由特定椭球体及其对应的转换参数确定，因此欲正确定义GIS系统坐标系，首先必须弄清地球椭球体(Ellipsoid)、大地基准面(Datum)及地图投影(Projection)三者的基本概念及它们之间的关系。 基准面是利用特定椭球体对特定地区地球表面的逼近，因此每个国家或地区均有各自的基准面，我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国的两个大地基准面。我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐标系，1978年采用国际大地测量协会推荐的1975地球椭球体建立了我国新的大地坐标系--西安80坐标系，目前大地测量基本上仍以北京54坐标系作为参照，北京54与西安80坐标之间的转换可查阅国家测绘局公布的对照表。 WGS1984基准面采用WGS84椭球体，它是一地心坐标系，即以地心作为椭球体中心，目前GPS测量数据多以WGS1984为基准。
上述3个椭球体参数如下： 椭球体与基准面之间的关系是一对多的关系，也就是基准面是在椭球体基础上建立的，但椭球体不能代表基准面，同样的椭球体能定义不同的基准面，如前苏联的Pulkovo 1942、非洲索马里的Afgooye基准面都采用了Krassovsky椭球体，但它们的基准面显然是不同的。 地图投影是将地图从球面转换到平面的数学变换，如果有人说：该点北京54坐标值为X=4231898,Y=21655933，实际上指的是北京54基准面下的投影坐标，也就是北京54基准面下的经纬度坐标在直角平面坐标上的投影结果。 2. GIS中基准面的定义与转换 虽然现有GIS平台中都预定义有上百个基准面供用户选用，但均没有我们国家的基准面定义。假如精度要求不高，可利用前苏联的Pulkovo 1942基准面(Mapinfo中代号为1001)代替北京54坐标系；假如精度要求较高，如土地利用、海域使用、城市基建等GIS系统，则需要自定义基准面。 GIS系统中的基准面通过当地基准面向WGS1984的转换7参数来定义，转换通过相似变换方法实现，具体算法可参考科学出版社1999年出
坐标转换模型
坐标转换模型 1.空间直角坐标系间的转换模型（七参数模型） ①公式（布尔莎模型）： ②分析： (1)将O-XYZ中的长度单位缩放l+m倍，使其与O'-X'Y'Z'的长度单位一致； (2)从X反向看向原点O，以O为旋转点，让O-XYZ绕X轴顺时针旋转Wx角，使经过旋转后的Y轴与O'-X'Y'Z’平面平行； (3)从Y反向看向原点O，以O为旋转点，让O-XYZ绕Y轴顺时针旋转Wy角，使经过旋转后的X轴与O'-X'Y'Z'平面平行。显然，此时Z轴也与Z'轴平行； (4)从Z反向看向原点O，以O点为旋转点，O-XYZ绕Z轴顺时针旋转Wz角，使经过旋转后的X轴与X’轴平行。显然，此时O-XYZ的三个坐标轴己与O'-X'Y'Z’中相应的坐标轴平行； 原坐标为O-XYZ，转换到新坐标O-X’Y’Z’.(两坐标系都为空间直角坐标系）其中（dX dY dZ）为坐标原点的平移参数，即将坐标O-XYZ的原点分别沿三个坐标轴平移-dX,-dY,-dZ，使原坐标轴与O-X’Y’Z’的点重合。m为尺度参数，（w1 w2 w3）分别为坐标轴的旋转参量（角度），构成的旋转矩阵分别为： 分别将R1 R2 R3代入上式，可得：
当旋转角度w1 w2 w3很小时（<=10),cos(w)=1,sin(w)=0;在误差允许范围内可以将模型简化为：(同样七参数模型） 四参数模型是在七参数模型的特例，没有考虑坐标轴的旋转量，只考虑坐标轴的平移。 总结： 类似布尔莎模型（以坐标原点为参考点），还有莫洛金斯基坐标模型（以目标点为变换中心）、武测转换模型和范士转换模型（以控制网参考点的站心地平坐标系的三个坐标轴为旋转轴），这些坐标转换模型很容易实现相关坐标在不同坐标系的转换，但是参考位置的偏移向量的相关参数，在实际运用中这些参量是很难测定的，并且受地球重力等物理因素的影响，两个坐标系统即使经过相似变换，仍可能存在较大的残差，所以这些模型适用于简单且规则模型中。 ④程序： clc clear all dX=input('please input value of dX=');
坐标转换之计算公式
创作编号： GB8878185555334563BT9125XW 创作者： 凤呜大王* 坐标转换之计算公式 一、参心大地坐标与参心空间直角坐标转换 1名词解释： A ：参心空间直角坐标系： a) 以参心0为坐标原点； b) Z 轴与参考椭球的短轴（旋转轴）相重合； c) X 轴与起始子午面和赤道的交线重合； d) Y 轴在赤道面上与X 轴垂直，构成右手直角坐标系0-XYZ ； e) 地面点P 的点位用（X ，Y ，Z ）表示； B ：参心大地坐标系： a) 以参考椭球的中心为坐标原点，椭球的短轴与参考椭球旋转轴重合； b) 大地纬度B ：以过地面点的椭球法线与椭球赤道面的夹角为大地纬度B ； c) 大地经度L ：以过地面点的椭球子午面与起始子午面之间的夹角为大地经度 L ； d) 大地高H ：地面点沿椭球法线至椭球面的距离为大地高H ； e) 地面点的点位用（B ，L ，H ）表示。 2 参心大地坐标转换为参心空间直角坐标： ?? ? ?? +-=+=+=B H e N Z L B H N Y L B H N X sin *])1(*[sin *cos *)(cos *cos *)(2 公式中，N 为椭球面卯酉圈的曲率半径，e 为椭球的第一偏心率，a 、b 椭球的长短半径，f 椭球扁率，W 为第一辅助系数
a b a e 2 2-= 或 f f e 1 *2-= W a N B W e = -=22 sin *1( 3 参心空间直角坐标转换参心大地坐标 [ ] N B Y X H H e N Y X H N Z B X Y L -+= +-++==cos ))1(**)() (*arctan() arctan(2 22 2 2 二 高斯投影及高斯直角坐标系 1、高斯投影概述 高斯-克吕格投影的条件：1. 是正形投影；2. 中央子午线不变形 高斯投影的性质：1. 投影后角度不变；2. 长度比与点位有关，与方向无关； 3. 离中央子午线越远变形越大 为控制投影后的长度变形，采用分带投影的方法。常用3度带或6度带分带，城市或工程控制网坐标可采用不按3度带中央子午线的任意带。 2、高斯投影正算公式：
直角坐标与极坐标的区别与转换
直角坐标 直角坐标系在数学中应用广泛，是数学大厦最重要的根基之一。 在平面内画两条 直角坐标 直角坐标 互相垂直，并且有公共原点的数轴。其中横轴为X轴，纵轴为Y轴。这样我们就说在平面上建立了平面直角坐标系，简称直角坐标系。 直角坐标中的点 直角坐标中的点 坐标：对于平面内任意一点C，过点分C别向X轴、Y轴作垂线，垂足在X 轴、Y轴上的对应点a,b分别叫做点C的横坐标、纵坐标，有序数对（a,b）叫做点C的坐标。坐标平面：坐标系所在平面。 坐标原点：两坐标轴的公共原点。 象限：X轴和Y轴把坐标平面分成四个象限，右上面的叫做第一象限，其他三个部分按逆时针方向依次叫做第二象限、第三象限和第四象限。象限以数轴为界，横轴、纵轴上的点不属于任何象限。
极坐标 极坐标系 polar coordinates 在平面内由极点、极轴和极径组成的坐标系。在平面上取定一点O，称为极点。从O出发引一条射线Ox，称为极轴。再取定一个长度单位，通常规定角度取逆时针方向为正。这样，平面上任一点P的位置就可以用线段OP的长度ρ以及从Ox到OP 的角度θ来确定，有序数对（ρ，θ）就称为P点的极坐标，记为P（ρ，θ）；ρ称为P 点的极径，θ称为P点的极角。当限制ρ≥0，0≤θ<2π时，平面上除极点Ο以外，其他每一点都有唯一的一个极坐标。极点的极径为零，极角任意。若除去上述限制，平面上每一点都有无数多组极坐标，一般地，如果（ρ，θ）是一个点的极坐标，那么（ρ，θ+2nπ），（－ρ，θ+（2n+1）π），都可作为它的极坐标，这里n 是任意整数。平面上有些曲线，采用极坐标时，方程比较简单。例如以原点为中心，r为半径的圆的极坐标方程为ρ=r 等速螺线的极坐标方程为ρ=aθ 。此外，椭圆、双曲线和抛物线这3种不同的圆锥曲线，可以用一个统一的极坐标方程表示。 极坐标系到直角坐标系的转化： 在极坐标系与平面直角坐标系（笛卡尔坐标系）间转换极坐标系中的两个坐标ρ和θ可以由下面的公式转换为直角坐标系下的坐标值 x=ρcosθ y=ρsinθ 由上述二公式，可得到从直角坐标系中x和y两坐标如何计算出极坐标下的坐标θ=arctany/x ( x不等于0) 在x= 0的情况下：若y为正数θ= 90° (π/2 radians)；若y为负，则θ= 270° (3π/2 radians). 极坐标的方程 用极坐标系描述的曲线方程称作极坐标方程，通常表示为r为自变量θ的函数。 极坐标方程经常会表现出不同的对称形式，如果r(?θ) = r(θ)，则曲线关于极点
大地坐标与直角空间坐标转换计算公式
坐标与直角空间坐标转换计算公式 一、参心坐标与参心空间直角坐标转换 1名词解释： A ：参心空间直角坐标系： a) 以参心0为坐标原点； b) Z 轴与参考椭球的短轴（旋转轴）相重合； c) X 轴与起始子午面和赤道的交线重合； d) Y 轴在赤道面上与X 轴垂直，构成右手直角坐标系0-XYZ ； e) 地面点P 的点位用（X ，Y ，Z ）表示； B ：参心坐标系： a) 以参考椭球的中心为坐标原点，椭球的短轴与参考椭球旋转轴重合； b) 纬度B ：以过地面点的椭球法线与椭球赤道面的夹角为纬度B ； c) 经度L ：以过地面点的椭球子午面与起始子午面之间的夹角为经度L ； d) 高H ：地面点沿椭球法线至椭球面的距离为高H ； e) 地面点的点位用（B ，L ，H ）表示。 2 参心坐标转换为参心空间直角坐标： ?? ? ?? +-=+=+=B H e N Z L B H N Y L B H N X sin *])1(*[sin *cos *)(cos *cos *)(2 公式中，N 为椭球面卯酉圈的曲率半径，e 为椭球的第一偏心率，a 、b 椭球的长短半径，f 椭球扁率，W 为第一辅助系数 a b a e 2 2-= 或 f f e 1*2-= W a N B W e = -=22 sin *1( 80椭球参数： 长半轴a=6378140±5（m ）
短半轴b=6356755.2882m 扁 率α=1/298.257 3 参心空间直角坐标转换参心坐标 [ ] N B Y X H H e N Y X H N Z B X Y L -+= +-++==cos ))1(**)() (*arctan() arctan(2 22 2 2 二 高斯投影及高斯直角坐标系 1、高斯投影概述 高斯-克吕格投影的条件：1. 是正形投影；2. 中央子午线不变形 高斯投影的性质：1. 投影后角度不变；2. 长度比与点位有关，与方向无关； 3. 离中央子午线越远变形越大 为控制投影后的长度变形，采用分带投影的方法。常用3度带或6度带分带，城市或工程控制网坐标可采用不按3度带中央子午线的任意带。 2、高斯投影正算公式： 52224253 2236 425442232)5814185(cos 120 )1(cos 6 cos )5861(cos sin 720 495(cos sin 24cos sin 2l t t t B N l t B N Bl N y l t t B B N l t B B N Bl B N X x ηηηηη-++-++-+=+-+++-++ =） 3、高斯投影反算公式：
#地理信息中各种坐标系区别和转换总结
地理信息中各种坐标系区别和转换总结 一、北京54坐标到西安80坐标转换小结 1、北京54和西安80是两种不同的大地基准面，不同的参考椭球体，因而两种地图下，同一个点的坐标是不同的，无论是三度带六度带坐标还是经纬度坐标都是不同的。 2、数字化后的得到的坐标其实不是WGS84的经纬度坐标，因为54和80的转换参数至今没有公布，一般的软件中都没有54或80投影系的选项，往往会选择WGS84投影。 3、WGS8 4、北京54、西安80之间，没有现成的公式来完成转换。 4、对于54或80坐标，从经纬度到平面坐标（三度带或六度带）的相互转换可以借助软件完成。 5、54和80间的转换，必须借助现有的点和两种坐标，推算出变换参数，再对待转换坐标进行转换。（均靠软件实现） 6、在选择参考点时，注意不能选取河流、等高线、地名、高程点，公路尽量不选。这些在两幅地图上变化很大，不能用作参考。而应该选择固定物，如电站，桥梁等。 二、西安80坐标系和北京54坐标系转换 西安80坐标系和北京54坐标系其实是一种椭球参数的转换作为这种转换在同一个椭球里的转换都是严密的，而在不同的椭球之间的转换是不严密，因此不存在一套转换参数可以全国通用的，在每个地方会不一样，因为它们是两个不同的椭球基准。那么，两个椭球间的坐标转换，一般而言比较严密的是用七参数布尔莎模型，即 X 平移， Y 平移， Z 平移， X 旋转（WX）， Y 旋转（WY）， Z 旋转（W Z），尺度变化（DM ）。要求得七参数就需要在一个地区需要 3 个以上的已知点。如果区域范围不大，最远点间的距离不大于 30Km（经验值），这可以用三参数，即 X 平移， Y 平移， Z 平移，而将 X 旋转， Y 旋转， Z 旋转，尺度变化面DM视为 0 。 在MAPGIS平台中实现步骤： 第一步：向地方测绘局（或其它地方）找本区域三个公共点坐标对（即54坐标x，y，z和80坐标x，y，z）； 第二步：将三个点的坐标对全部转换以弧度为单位。（菜单：投影转换/输入单点投影转换，计算出这三个点的弧度值并记录下来） 第三步：求公共点求操作系数（菜单：投影转换/坐标系转换）。如果求出转换系数后，记录下来。 第四步：编辑坐标转换系数。（菜单：投影转换/编辑坐标转换系数。）最后进行投影变换，“当前投影”输入80坐标系参数，“目的投影”输入54坐标系参数。进行转换时系统会自动调用曾编辑过的坐标转换系数。 三、地理坐标系和投影坐标系的区别 1、首先理解地理坐标系（Geographic coordinate system），Geographic coordinate system直译为地理坐标系统，是以经纬度为地图的存储单位的。很明显，Geographic coordinate system是球面坐标系统。我们要将地球上的数字化信息存放到球面坐标系统上，如何进行操作呢？地球是一个不规则的椭球，如何将数据信息以科学的方法存放到椭球上？这必然要求我们找到这样的一个椭球体。这样的椭球体具有特点：可以量化计算的。具有长半轴，短 半轴，偏心率。以下几行便是Krasovsky_1940椭球及其相应参数。
不同坐标系之间的变换
不同坐标系之间的变换 SANY GROUP system office room 【SANYUA16H-
§10.6不同坐标系之间的变换 10.6.1欧勒角与旋转矩阵 对于二维直角坐标，如图所示，有： ?? ? ?????????-=??????1122cos sin sin cos y x y x θθθθ(10-8) 在三维空间直角坐标系中，具有相同原点的两坐标系间的变换一般需要在三个坐标平面上，通过三次旋转才能完成。如图所示，设旋转次序为： ①绕1OZ 旋转Z ε角，11,OY OX 旋 转至0 0,OY OX ； ②绕0 OY 旋转Y ε角 10 ,OZ OX 旋转至0 2 ,OZ OX ； ③绕2OX 旋转X ε角， 0,OZ OY 旋转至22,OZ OY 。 Z Y X εεε,,为三维空间直角坐标变换的三个旋转角，也称欧勒角，与 它相对应的旋转矩阵分别为： ???? ? ?????-=X X X X X R εεεεεcos sin 0sin cos 00 01 )(1 (10-10)
????? ?????-=Y Y Y Y Y R εεεεεcos 0sin 010sin 0cos )(2 (10-11) ???? ? ?????-=10 0cos sin 0sin cos )(3Z Z Z Z Z R εεεεε (10-12) 令 )()()(3210Z Y X R R R R εεε= (10- 13) 则有： ???? ? ?????=??????????=??????????1110111321222)()()(Z Y X R Z Y X R R R Z Y X Z Y X εεε (10-14) 代入： ???? ??? ??? +-+++--=Y X Z Y X Z X Z Y X Z X Y X Z Y X Z X Z Y X Z X Y Z Y Z Y R εεεεεεεεεεεεεεεεεεεεεεεεεεεεεcos cos sin sin cos cos sin cos sin cos sin sin cos sin sin sin sin cos cos cos sin sin sin cos sin sin cos cos cos 0一般Z Y X εεε,,为微小转角，可取： sin sin sin sin sin sin sin ,sin ,sin 1cos cos cos =========Z Y Z X Y X Z Z Y Y X X Z Y X εεεεεεεεεεεεεεε 于是可化简
平面直角坐标变换
平面直角坐标变换 【摘要】对利用EXCEL电子表格进行高斯投影换算的方法进行了较详细的介绍，对如何进行GPS坐标系转换进行了分析，提出了一种简单实用的坐标改正转换方法，介绍了用EXCEL完成转换的思路。 [关键字] 电子表格；GPS；坐标转换 作为尖端技术GPS，能方便快捷性地测定出点位坐标，无论是操作上还是精度上，比全站仪等其他常规测量设备有明显的优越性。随着我国各地GPS差分台站的不断建立以及美国SA政策的取消，使得单机定位的精度大大提高，有的已经达到了亚米级精度，能够满足国土资源调查、土地利用更新、遥感监测、海域使用权清查等工作的应用。在一般情况下，我们使用的是1954年北京坐标系或1980年西安坐标系（以下分别简称54系和80系），而GPS测定的坐标是WGS-84坐标系坐标，需要进行坐标系转换。对于非测量专业的工作人员来说，虽然GPS定位操作非常容易，但坐标转换则难以掌握，EXCEL是比较普及的电子表格软件，能够处理较复杂的数学运算，用它来进行GPS坐标转换、面积计算会非常轻松自如。要进行坐标系转换，离不开高斯投影换算，下面分别介绍用EXCEL进行换算的方法和GPS 坐标转换方法。 一、用EXCEL进行高斯投影换算 从经纬度BL换算到高斯平面直角坐标XY（高斯投影正算），或从XY换算成BL（高斯投影反算），一般需要专用计算机软件完成，在目前流行的换算软件中，存在一个共同的不足之处，就是灵活性较差，大都需要一个点一个点地进行，不能成批量地完成，给实际工作带来许多不便。笔者发现，用EXCEL可以很直观、方便地完成坐标换算工作，不需要编制任何软件，只需要在EX CEL的相应单元格中输入相应的公式即可。下面以54系为例，介绍具体的计算方法。 完成经纬度BL到平面直角坐标XY的换算，在EXCEL中大约需要占用21列，当然读者可以通过简化计算公式或考虑直观性，适当增加或减少所占列数。在EXCEL中，输入公式的起始单元格不同，则反映出来的公式不同，以公式从第2行第1列（A2格）为起始单元格为例，各单元格的公式如下： 单元格 单元格内容 说明A2 输入中央子午线，以度.分秒形式输入，如115度30分则输入1 15.30 起算数据L0 B2 =INT(A2)+(INT(A2*100)-INT(A2)*100)/60+(A2*10000-INT(A2* 100)*100)/3600 把L0化成度 C2 以度小数形式输入纬度值，如38°14′20″则输入38.1420 起算数据B D2 以度小数形式输入经度值 起算数据L E2 =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-B2 L-L0 H2 =G2/57.2957795130823 化作弧度 I2 =TAN(RADIANS(E2)) Tan(B) J2 =COS(RADIANS(E2)) COS(B)