excel转换坐标文件的公式

在excel表格中用公式计算测量数据

用坐标计算方位角公式(格式为：“度.分秒”)：

=180+N(TRUNC(ATAN((I5-I4)/(H5-H4))*180/PI()))+N(TRUNC((ATAN((I5-I4)/(H 5-H4))*180/PI()-TRUNC(ATAN((I5-I4)/(H5-H4))*180/PI()))*60))/100+ROUND(((A TAN((I5-I4)/(H5-H4))*180/PI()-TRUNC(ATAN((I5-I4)/(H5-H4))*180/PI()))*60-TR UNC((ATAN((I5-I4)/(H5-H4))*180/PI()-TRUNC(ATAN((I5-I4)/(H5-H4))*180/PI())) *60))*60,0)/10000

“度.分秒”转换为“度”的计算公式：

=INT(B6)+(INT(MOD(B6,1)*100)/60)+(MOD(B6,0.01)*10000/3600)

“度”转换为“度.分秒”的计算公式：

=INT(C6)+INT((C6-INT(C6))*60)/100+((C6-INT(C6))*60- INT((C6-INT(C6))*60))*60/10000

坐标增量计算：

△x=F6*COS(D6*PI()/180)

△y=F6*SIN(D6*PI()/180)

两坐标点间长度计算：

=SQRT((C5-C4)*(C5-C4)+(B5-B4)*(B5-B4))

“度”转换为“弧度”的计算公式：

=RADIANS(d5) or =d5*pi()/180

“弧度”转换为“度”的计算公式：

=d5*180/pi()

求方位角公式

度格式：

=(PI()*(1 - SIGN(B3-$B$1) / 2) - ATAN((A3-$A$1) /(B3-$B$1)))*180/PI()

Excel 中求方位角公式：a1,b1放起始点坐标 a3,b3放终点坐标。

度分秒格式：

=INT((PI()*(1 - SIGN(B3-$b$1) / 2) - ATAN((A3-$a$1) /(B3-$b$1)))*180/PI()) &"-"& INT( ((PI()*(1 - SIGN(B3-$b$1) / 2) - ATAN((A3-$a$1) /(B3-$b$1)))*180

/PI()-INT((PI()*(1 - SIGN(B3-$b$1) / 2) - ATAN((A3-$a$1) /(B3-$b$1)))*180/

PI()))*60)&"-"&INT( (((PI()*(1 - SIGN(B3-$b$1) / 2) - ATAN((A3-$a$1) /(B3-

$b$1)))*180/PI()-INT((PI()*(1 - SIGN(B3-$b$1) / 2) - ATAN((A3-$a$1) /(B3-

$b$1)))*180/PI()))*60-INT(((PI()*(1 - SIGN(B3-$b$1) / 2) - ATAN((A3-$a$1) /

(B3-$b$1)))*180/PI()-INT((PI()*(1 - SIGN(B3-$b$1) / 2) - ATAN((A3-$a$1) /

(B3-$b$1)))*180/PI()))*60))*600)/10

其中：A1,B1中存放测站坐标,a3,b3放终点坐标。

上面的计算出来的是度分秒格式，也就是字符串格式，不能用来计算，只是用来看的哟！下面这个简单一点：

=(PI()*(1 - SIGN(B3-B1) / 2) - ATAN((A3-A1) /(B3-B1)))*180/PI()

Excel 中求方位角公式：a1,b1放起始点坐标 a3,b3放终点坐标。

求距离公式:

=Round(SQRT(POWER((A3-$A$1),2)+POWER((B3-$B$1),2)),3)

【最后一个3表示取的小数点位数，最大为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度带的带号。

坐标转换之计算公式

坐标转换之计算公式一、参心大地坐标与参心空间直角坐标转换 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、高斯投影反算公式：

大地坐标转换成施工坐标公式

大地(高斯平面)坐标系工程坐标系转换大地坐标系--->工程坐标系 ======================== 待转换点为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”

平面直角坐标变换

平面直角坐标变换【摘要】对利用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)

大地坐标与直角空间坐标转换计算公式

大地坐标与直角空间坐标转换计算公式一、参心大地坐标与参心空间直角坐标转换 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、高斯投影反算公式：

经纬度转换公式

经纬度转换公式一、用EXCEL进行高斯投影换算从经纬度BL换算到高斯平面直角坐标XY（高斯投影正算），或从XY换算成BL（高斯投影反算），一般需要专用计算机软件完成，在目前流行的换算软件中，存在一个共同的不足之处，就是灵活性较差，大都需要一个点一个点地进行，不能成批量地完成，给实际工作带来许多不便。笔者发现，用EXCEL可以很直观、方便地完成坐标换算工作，不需要编制任何软件，只需要在EXCEL的相应单元格中输入相应的公式即可。下面以54系为例，介绍具体的计算方法。完成经纬度BL到平面直角坐标XY的换算，在EXCEL中大约需要占用21列，当然读者可以通过简化计算公式或考虑直观性，适当增加或减少所占列数。在EXCEL中，输入公式的起始单元格不同，则反映出来的公式不同，以公式从第2行第1列（A2格）为起始单元格为例，各单元格的公式如下：单元格单元格内容说明 A2 输入中央子午线，以度.分秒形式输入，如115度30分则输入115.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

坐标转换模型

坐标转换模型 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、高斯投影正算公式：

大地坐标和经纬度之间的换算完整版

大地坐标和经纬度之间的换算 Document serial number【NL89WT-NY98YT-NC8CB-NNUUT-NUT108】

大地坐标和经纬度之间的换算地质工作中常要对进行大地坐标转经纬度和经纬度换大地坐标,我写一下一般的过程,希望对大家有点帮助. 大地座标-----→经纬度(地理坐标) 1,输入大地坐标数据，格式为 Y空格X，输入到文本就行如下，原始的大地坐标由一8位的Y和一个7位的X组成，这组坐标数据中的Y的前两位为31，是分带号，一般使用的分带有三分带，六分带，这里的坐标是三分带的，记下Y前的这两位数，在原始数据中去除掉，现在数据变为，Y，6位，X，7位保存这个TXT的文本文件。 2，打开MAPGIS，启动坐标投影变形程序接下来选择投影转换>>>>用户文件投影转换点打开文件，打开刚才的大地坐标的文本文件，设置输入数据的格式，点击用户投影参数，并完成设置。这里我们的大地座标为3度带的第31带，注意填好，坐标单位为米好了以为设置输出的格式，我们要求输出的是经纬度，点结果转换参数，完成以下设置我们输出的经纬度的单位应该是DDDMMMSS。SS注意点写到文件，保存就大功告成了，注意保存的文件要写上.TXT的后缀下面是计算出的结果文件 XP为经度，1234234。357就是123度42分34。357秒， YP为纬度，403950。225就是40度39分50。255秒（纬度没有最多90，所以没有三位数） ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 经纬度——→大地座标输入文件格式如下，这里面的数据前面的为经度，格式为ＤＤＤＭＭＳＳ，后面的为纬度，格式为ＤＤＭＭＳＳ接下来的转换过程和大地坐标转换一样，只要将刚才的用户转换参数和结果转换参数交换即可，要注意分带号的确定，如果你不知道分带号，就应该先计算分带号，算法是经度／３得到的整数为三度带的分带号经度／６得到的整数为六度带的分带号其中的ＸＰ为地图上的Ｙ坐标，记得在前面加上带号，其中的ＹＰ为地图上的Ｘ坐标

不同坐标系之间的变换

不同坐标系之间的变换 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 εεεεεεεεεεεεεεε 于是可化简

推导坐标旋转公式

推导坐标旋转公式数学知识2010-09-12 21:03:53 阅读151 评论0 字号：大中小订阅在《Flash actionScript 3.0 动画教程》一书中有一个旋转公式： x1=cos(angle)*x-sin(angle)*y; y1=cos(angle)*y+sin(angle)*x; 其中x，y表示物体相对于旋转点旋转angle的角度之前的坐标，x1，y1表示物体旋转angle 后相对于旋转点的坐标从数学上来说，此公式可以用来计算某个点绕另外一点旋转一定角度后的坐标，例如：A（x，y）绕B（a，b)旋转β度后的位置为C（c，d），则x，y，a，b，β，c，d有如下关系式： 1。设A点旋转前的角度为δ，则旋转(逆时针)到C点后角度为δ+β 2。求A，B两点的距离：dist1=|AB|=y/sin(δ)=x/cos(δ) 3。求C，B两点的距离：dist2=|CB|=d/sin(δ+β)=c/cos(δ+β) 4。显然dist1=dist2，设dist1=r所以： r=x/cos(δ)=y/sin(δ)=d/sin(δ+β)=c/cos(δ+β) 5。由三角函数两角和差公式知： sin(δ+β)=sin(δ)cos(β)+cos(δ)sin(β) cos(δ+β)=cos(δ)cos(β)-sin(δ)sin(β) 所以得出：

c=r*cos(δ+β)=r*cos(δ)cos(β)-r*sin(δ)sin(β)=xcos(β)-ysin(β) d=r*sin(δ+β)=r*sin(δ)cos(β)+r*cos(δ)sin(β)=ycos(β)+xsin(β) 即旋转后的坐标c，d只与旋转前的坐标x，y及旋转的角度β有关从图中可以很容易理解出A点旋转后的C点总是在圆周上运动，圆周的半径为|AB|，利用这点就可以使物体绕圆周运动，即旋转物体。上面公式是相对于B点坐标来的，也就是假如B点位（0,0）可以这么做。现在给出可以适合任意情况的公式： x0 = dx * cos(a) - dy * sin(a) y0 = dy * cos(a) + dx * sin(a) 参数解释： x0,y0是旋转后相对于中心点的坐标，也就是原点的坐标，但不是之前点旋转后的实际坐标，还要计算一步，a旋转角度，可以是顺时针或者逆时针。 dx是旋转前的x坐标-旋转后的x坐标 dy是旋转前的y坐标-旋转后的y坐标 x1=b+x0; y1=c+y0; 上面才是旋转后的实际坐标，其中b，c是原点坐标下面是上面图的公式解答： x0=(x-b)*cos(a)-(y-c)*sin(a); y0=(y-c)*cos(a)+(x-b)*sin(a); x1=x0+b; y1=y0+c;

GPS经纬度及小数点如何换算

GPS经纬度及小数点如何换算经度最大是180度,纬度最大是90度.一度等于60分,一分等于60秒.这个关系有了什么都能算.。例：如果测量出120度18分30秒秒换算成分用30秒/60=0.5分就换算成120度18.5分想要熟悉GPS中的经纬度，就必须先了解GPS，知道经纬度的来源： 1. GPS系统组成 GPS是Gloabal Positioning System 的简称，意为全球定位系统，主要由地面的控制站、天上飞的卫星、咱们手里拿的接收机三大块组成，我们所使用的GPS 包括手持机和车载导航机本质上都是GPS接受机。 2. GPS接收机接收机大大小小，千姿百态，有袖珍式、背负式、车载、船载、机载什么的。一般常见的手持机接收L1信号，还有双频的接收机，做精密定位用的。 3. 坐标系地形图坐标系：我国的地形图采用高斯－克吕格平面直角坐标系。在该坐标系中，横轴：赤道，用Ｙ表示；纵轴：中心经线，用Ｘ表示；坐标原点：中心经线与赤道的交点，用Ｏ表示。赤道以南为负，以北为正；中心经线以东为正，以西为

负。我国位于北半球，故纵坐标均为正值，但为避免中心经度线以西为负值的情况，将坐标纵轴西移５００公里。北京５４坐标系：１９５４年我国在北京设立了大地坐标原点，采用克拉索夫斯基椭球体，依此计算出来的各大地控制点的坐标，称为北京５４坐标系。ＧＳ８４坐标系：即世界通用的经纬度坐标系。６度带、３度带、中心经线。 https://www.360docs.net/doc/861465434.html, 我国采用６度分带和３度分带：１∶２．５万及１∶５万的地形图采用６度分带投影，即经差为６度，从零度子午线开始，自西向东每个经差６度为一投影带，全球共分６０个带，用１，２，３，４，５，……表示。１∶１万的地形图采用３度分带，从东经１．５度的经线开始，每隔３度为一带，用１，２，３，……表示，全球共划分１２０个投影带 4. 经纬度的来源为了精确地表明各地在地球上的位置，人们给地球表面假设了一个坐标系，这就是经纬度线。那么，最初的经纬度线是怎么产生？又是如何测定的呢？公元344年，亚历山大渡海南侵，继而东征，随军地理学家尼尔库斯沿途搜索资料，预备绘一幅“世界地图”。他发现沿着亚历山大东征的路线，由西向东，无论季节变换与日照长短都很相仿。于是做出了一个重要贡献——第一次在地球上划出了一条纬线，

坐标系转换步骤以及公式

一、各坐标系下椭球参数二、WGS84转北京54一般步骤（转80一样，只是椭球参数不同）前期工作：收集测区高等级控制点资料。在应用手持GPS 接收机观测的区域内找出三个以上分布均匀的等级点（精度越高越好）或GPS “ B ”级网网点，点位最好是周围无电磁波干扰，视野开阔，卫星信号强。并到测绘管理部门抄取这些点的54北京坐标系的高斯平面直角坐标(x 、y)，大地经纬度（B 、L ），高程h ，高程异常值ξ和WGS-84坐标系的大地经纬度（B 、L ），大地高H 。如果没有收集到WGS-84下的大地坐标，则直接用手持GPS 测定已知点B 、L 、H 值。转换步骤： 1、把从GPS 中接收到84坐标系下的大地坐标（经纬度高程B 、L, H ，其中B 为纬度，L 为经度，H 为高程），使用84坐标系的椭球参数转换为84坐标系下的地心直角坐标（空间坐标）：式中，N 为法线长度，为椭球长半径，b 为椭球短半径，为第一偏心率。 2、使用七参数转换为54坐标系下的地心直角坐标（x ，y ，z ）： x = △x + k*X- β*Z + γ*Y+ X y = △y + k*Y + α*Z - γ*X + Y z = △z + k*Z - α*Y + β*X + Z

其中，△x，△y，△z为三个坐标方向的平移参数；α，β，γ为三个方向的旋转角参数；k为尺度参数。（采用收集到的控制点计算转换参数，并需要验证参数）在小范围内可使用七参数的特殊形式即三参数，即k、α、β、γ都等于0，变成： x = △x+ X y = △y+ Y z = △z + Z 3、根据54下的椭球参数，将第二步得到的地心坐标转换为大地坐标（B54,L54,H54）计算B时要采用迭代，推荐迭代算法为： 4、根据工程需要以及各种投影（如高斯克吕格）规则进行投影得到对应的投影坐标，即平面直角坐标。（投影正算）三、北京54转WGS84一般步骤（80转84一样，只是椭球参数不同） 1、将所有点的BJ54高斯平面直角坐标（x，y）化算为大地坐标(B,L )。（投影反算） 2、顾及水准高h后将三维大地坐标（B，L，h），按54椭球参数化算为地心直角坐标（X，Y，Z ）。（公式同上面第一步） 3、根据公共点求转换七参数或多项式拟合系数并将54下的（X,Y,Z）转为84 下的（X,Y,Z）。（公式同上面第二步）. 4、将转换后的三维直角坐标WGS-84XYZ化算为大地坐标WGS-84(BLH) 。（公式同上面第三步） 5 、引入基于WGS-84椭球的高程异常值由水准高求得基于WGS-84椭球的大地高H 。

经纬度转换XYZ

我们用经纬度来表示地球上的点，用三维体系来表示空间中的点，那么在地球上构造一个三维体系，是否会对我们有很大的帮助？我们要想让三维体系使我们最方便，那应该把三维体系构造在地球上的哪个位置？三维体系坐标与地球经纬度又有什么关系？本文将会对这些问题进行探索。【关键字】：经纬线三维体系坐标关系 1 前言由于地球的极半径与赤道半径长度相差不多，所以我们通常地把地球看成是一个规则球体，本文将会把三维体系建立在地球地心处，让z轴与地轴重合，y轴与0—180度经线圈于0度纬线的直径重合，则x轴与90E—90W经线圈于0度纬线的直径重合。如图① 图① 然后通过一定的计算，让地球的点（纬度，经度，海拔）表示为三维体系中的点（x,y,z）。 2 地球经纬度与三维体系坐标的关系是什么设纬度α，经度为β，海拔为h（米）。若地球上任意一点为（α,β,h）表示为三维体系中的点（x,y,z），则：东经：x≈(6372.8+h/1000)(cos α)(sin β)西经：x=(6372.8+h/1000)(cos α)sin (180+β) y≈(6372.8+h/1000)(cos α)(cos β) 北纬：z≈(6372.8+h/1000)(sin α)南纬：z=(6372.8+h/1000)sin (90+α) 海拔计算时单位是米，x,y,z单位是千米。 3 如何证明“地球经纬度与三维体系坐标的关系” 切换仰视视角：

点A是南半球的一点，作与A同一海拔，同一经度，纬度为0度的点B，如图③ 海拔=h米=h/1000千米 ∴B到地心的距离BO=地球平均半径+h/1000 ≈6372.8+h/1000千米切换回立体的视角图④ 作AC⊥BO，C为垂直，r为平均半径

不同坐标系之间的变换

???? ??????-=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 εεεεεεεεεεεεεεε 于是可化简 ???? ? ?????---=111 0X Y X Z Y Z R εεεεεε (10-16) 上式称微分旋转矩阵。

经纬度和平面坐标的相互转换

经纬度和平面坐标的相互转换首先，RTK中测量的坐标，想要再转换成经纬度，很简单，其实不用转。直接打开所测量的工程，在坐标管理库中有数据导出的功能，可以直接导出来你所测的每一个点的经纬度或者是把工程文件中的.RTK文件复制出来，用EXCEL表格打开，直接提取经纬度即可。如果是别人提供的平面坐标，或者以前测量的原始文件删除了，那就需要通过软件进行转换来获得大地坐标。 GPStool GPS工具箱是常用的坐标转换软件，已上传到百度网盘。第一步，打开软件新建作业——起名，保存

第二步，设置转换参数源椭球，肯定选择WGS84。目标椭球，以西安80为例

投影参数设置，中心经度，也就是中央子午线，输入进去，在这里中央子午线的输入格式是“度.分秒”格式，（举例126度02分03.55秒就输126.020355 注意千万不要输成126.2355 中间的0不能少），其他一般不用。四参数，校正参数，拟合参数，七参数这四个参数，有哪些，输哪些。一般都是四参数+高程拟合参数+校正参数或者七参数+校正参数。这里以校正参数为例，记得使用校正参数一定打勾！

转换前坐标设置为大地坐标，格式有多种选择，一般选度或者无格式。度就是度的格式，举个例子，125度30分，度的格式下，就应该输125.5。无格式的情况下，就输入125.30或者125.3（末位的0可以不用输）说到这，告诉大家怎么区分“度”和“度.分秒”的区别，当你拿到一些经纬度时，出现43.6579 125.7484这种情况，肯定就是“度”的格式了，因为度分秒中，分秒不会大于60的。如果所有的经纬度中，小数点后第一位和第三位都小于6，那么基本就可以判断，这是“度.分秒”的格式。转换后类型，根据需要，自己选择。经纬度转平面，转换前椭球选择WGS84，类型选择大地坐标，转换后椭球选择80，类型选择平面坐标。平面转经纬度，转换前椭球选择80，类型选择平面坐标，转换后椭球选择WGS84，类型选择大地坐标。这是单点转换的例子，很简单，输入经纬度点箭头，就出平面。

坐标换算公式

关于A u t o C A D在变电所工程测量定位中的一些应用在变电所土建总平面图上，所内各建（构）筑物的平面位置系用施工坐标系和测量坐标系分别表示的。变电所的施工坐标系的原点一般虚设在变电所围墙的西南角上，从而使所内所有建（构）筑物的坐标皆为正值；而整个变电所的整体平面位置则用测量坐标系来表示，测量坐标系统系平面直角坐标，一般有国家坐标系统、城市坐标系统等。所以在我们进行变电所工程定位及所内施工控制网的布设时就需要将施工坐标系统与测量坐标系统进行转换计算。在我们以往的工程施工中，较为常用的是采用施工坐标系统与测量坐标系统的转换公式进行换算，但是这种较为繁琐的公式计算，包括距离计算、角度计算，工作量大，且很容易出现计算错误或计算精度达不到施工要求的问题。现在我向大家介绍一种无需进行公式计算，仅使用Auto CAD进行变电所两种坐标系统自由转换的方法，我们以110kV宁阳变电所工程的测量定位为例，分别采用上述两种方法进行计算、比较。根据《总平面布置图》宁阳变电所的测量坐标系统是1954年的北京坐标系。所区建（构）筑物采用施工坐标系统，取变电所西南角围墙轴线交点为施工坐标系的原点（A=0.00 B=0.00），其中A=0.00相当于X0=3386346.750m，B=0.00相当于Y0=496024.938m，A轴与指北针的夹角为北偏西18°。设计院交桩记录中给出的三个城市测量坐标控制点分别为：A1（X1=3386375.145m Y1=496019.325m）、A2（X2=3386418.782m Y2=496011.617m）、A3（X3=3386462.756m Y3=495977.459m）。在实际工程施工中需要根据施工场地的面积、建（构）筑物的位置及实际施工需要，布设四个控制桩点作为施工控制网：K1（A=48.00 B=10.00）、K2（A=48.00 B=38.00）、K3（A=77.00 B=10.00）、K4（A=77.00 B=38.00）。为了测定上述四个主控网控制桩点，则需要根据设计院提供的测量控制桩点进行坐标换算，从而在同一坐标系中进行距离及角度计算。由于两个坐标系的旋向不同，则施工坐标系与测量坐标系之间的换算关系式为： X＝X0＋A*COSθ－B*SINθ Y＝Y0＋A*SINθ＋B*COSθ 式中：θ＝-18° 代入坐标数据 K1： X1＝3386346.750＋48.00×COS（-18°）－10.00×SIN（-18°）＝3386395.491 Y1＝496024.938＋48.00×SIN（-18°）＋10.00×COS（-18°）＝496019.616

excel经纬度转距离公式文库

excel经纬度转距离公式知道两个点的经纬度，怎么用excel转换成距离？例如：A(118°19'20",35°4'4");B(118°19'56”,35°4’46”) ======================函数分割线========================= A1 : 第一点经度 B1 ：第一点纬度 A2 : 第二点经度 B2 ：第二点纬度经纬度格式：118°19'20" （度分秒的字符不要搞错）如： 118°19'20" 35°4'4" 118°19'56" 35°4'46" 计算结果是：1708.610943 米。当然，将地球视作标准圆球 =6371000*ACOS(COS(RADIANS(SUM(1*LEFT(A2,FIND("°",A2)-1),MID(A2,FIND("°",A2)+1,FIND("'",A2)-FIND("°",A2)-1)/60,RIGHT(LEFT(A2,LEN(A2)-1),LEN(A2)-FIND("'",A2)-1)/3600)-SUM(1*LEFT(A1 ,FIND("°",A1)-1),MID(A1,FIND("°",A1)+1,FIND("'",A1)-FIND("°",A1)-1)/60,RIGHT(LEFT(A1,LEN(A1)-1),LEN(A1)-FIND("'",A1)-1)/3600)))*COS(RADIANS (SUM(1*LEFT(B2,FIND("°",B2)-1),MID(B2,FIND("°",B2)+1,FIND("'",B2)-FIND("°",B2)-1)/60,RIGHT(LEFT(B2,LEN(B2)-1),LEN(B2)-FIND("'",B2)-1)/3600)-SUM(1*LEFT(B1 ,FIND("°",B1)-1),MID(B1,FIND("°",B1)+1,FIND("'",B1)-FIND("°",B1)-1)/60,RIGHT(LEFT(B1,LEN(B1)-1),LEN(B1)-FIND("'",B1)-1)/3600)))) ==========================算法分割线========================= 假设A点经纬度坐标为（a0，a1），B点经纬度坐标为（b0，b1），地球半径为R。则理论上AB两点间的弧长为 R * arccos[cos(b0-a0) * cos(b1-a1)] =======================详细算法分割线======================== 用EXCEL进行高斯投影换算从经纬度BL换算到高斯平面直角坐标XY（高斯投影正算），或从XY换算成BL（高斯投影反算），一般需要专用计算机软件完成，在目前流行的换算软件中，存在一个共同的不足之处，就是灵活性较差，大都需要一个点一个点地进行，不能成批量地完成，给实际工作带来许多不便。其实用EXCEL就可以很直观、方便地完成坐标换算工作，不需要编制任何软件，只需要在EXCEL的相应单元格中输入相应的公式即可。下面以54系为例，介绍具体的计算方法。完成经纬度BL到平面直角坐标XY的换算，在EXCEL中大约需要占用21列，当然读者可以通过简化计算公式或考虑直观性，适当增加或减少所占列数。在EXCEL中，输入公式的起始单元格不同，则反映出来的公式不同，以公式从第2行第1列（A2格）为起始单元格为例，各单元格的公式如下：