坐标转换快捷法——用EXCEL计算

坐标转换快捷法——用EXCEL计算
坐标转换快捷法——用EXCEL计算

1 概述

GPS 测量作为现代高新技术手段,能方便快捷的测定点位坐标,在操作上比全站仪等其

他常规测量设备有明显的优越性。随着我国各地GPS 连续运行参考站(CORS)的不断建设,

使得单机作业模式得到大范围推广,在CORS 的支持下定位的精度大大提高,已经达到了亚

米级、厘米级精度,能够满足国土资源调查、土地利用更新、道路工程测量、遥感监测等工作的应用。

一般情况,我们使用的是北京54 坐标系统或西安80 坐标系统,而GPS 测定的坐标是WGS-84 坐标系坐标,需要进行坐标系转换。对于非测量专业的工作人员来说,虽然GPS 定

位操作非常容易,但坐标转换则难以掌握,EXCEL 是比较普及的电子表格软件,能够处理较

复杂的数学运算,用它来进行GPS 坐标转换会非常轻松自如。

要进行坐标系转换,离不开高斯投影换算,下面分别介绍用EXCEL 进行换算的方法将GPS 坐标转换为平面坐标格式。

2 坐标转换

从经纬度BL 换算到高斯平面直角坐标XY(高斯投影正算),或从XY 换算成BL(高斯

投影反算),一般需要专用计算机软件完成,在目前流行的换算软件中,存在一个共同的不足之处,就是灵活性较差,大都需要一个点一个点地进行,不能成批量地完成,而具有批量转换功能的软件大多是收费的,价格不菲,这样给实际的工作带来许多不便。

2.1 用EXCEL 表格进行坐标转换

EXCEL 可以处理复杂的数学运算,可以很直观、方便地完成坐标换算工作,不需要写代码,不需要编制任何软件,只需要在EXCEL 的相应单元格中输入相应的公式即可。下面以北

京54 坐标系统为例,介绍具体的计算方法。

本方法完成经纬度坐标BL 到平面直角坐标XY 的换算,在EXCEL 表格中大约需要占用21 列,也可以通过简化计算公式或考虑直观性,适当减少或增加所占列数。在EXCEL 中,

输入公式的起始单元格不同,则反映出来的公式不同,以公式从第2 行第1 列(A2 格)为起始单元格为例,各单元格的公式如下:

表1 单元格公式说明

按上面表格中的公式输入到相应单元格后,就可方便地由经纬度求得平面直角坐标。当

输入完所有的经纬度后,用鼠标下拉一拖即可得到所有的计算结果。表中的许多单元格公式为中间过程,可以用EXCEL 的列隐藏功能把这些没有必要显示的列隐藏起来,表面上形成

准的计算报表,使整个计算表简单明了。

2.2 用例

用此表格可方便GPS 定位的GPGGA 信息离散度分析,我们知道,GPS 定位坐标的离散度从一定程度上代表着接收机定位结果的精度,通过此数据常用来实际分析某款GPS 接收机平面坐标的可靠性。

分析GPGGA 定位信息,将保存的数据导入到EXCEL 表格中,只保留B L H,在EXCEL 中

填写对应单元格的内容,即可计算出北京54 坐标系下的平面坐标X Y

图1 计算X Y 坐标

这样通过平面坐标的方式就更加直观的体现GPS 接收机定位结果的精度了,可利用EXCEL 继续分析坐标的离散度,在此不做赘述。

3 结论

现代PC 机的性能越来越强劲,在EXCEL 的支持下,从理论上讲,可计算的数据量是无限的,当第一次输入公式后,相当于自己完成了一软件的编制,可另存起来供今后重复使用,非常方便。

空间大地坐标系与平面直角坐标系转换公式

§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右侧所示。 图2-5 高斯投影 x 方向指北,y 方向指东。 可见,高斯投影存在长度变形,为使其在测图和用图时影响很小,应相隔一定的地区,另立中央子午线,采取分带投影的办法。我国国家测量规定采用六度带和三度带两种分带方法。六度带和三度带与中央子午线存在如下关系: 366 N L =中; n L 33=中 其中,N 、n 分别为6度带和3度带的带号。

线元法线路坐标正反算程序

经苦心钻研,奋战多日,终于编写出了代码短,速度快,精度高,功能全的线路坐标正反算程序,欢迎试用并提出宝贵意见。 功能简介及特点: 1、选用高斯-勒让德公式作计算内核,保证精度,模块化设计,便于扩充功能。 2、线元数据可自动从数据库调用,也可手工输入。 3、可管理多条线路,如里程不在线路或线元范围,将警告里程偏大、偏小。 4、边桩计算设计为导线式递推方式,可用于由一个中桩推出结构物所有角点坐标。 5、反算实现了智能化操作,只需输入线路号(或手工输线元资料)、坐标,不需近似里程,即可自动从起点向后开始试算出里程、位置,如对算出里程、位置表示怀疑,还可以让计算器从终点起再向前试算下一个可能的位置(匝道、回头曲线同一坐标可能会有一个以上结果)。第三次及以后试算才要求输入近似里程。 6、程序代码规范简洁,便于阅读、理解。 完整程序清单: ZFS %正反算主程序 B=.1739274226:C=.5-B: Lbl 1:U"0 ZS 1 FS"=0=>Prog "ZS": ≠>U=1=>Prog"FS":≠>Goto 1

ZS %正算子程序 {K}:Prog"ZZ":I=0:{I}:I"L"≠0=>"Prog"WY":≠>Prog"ZB" FS %反算子程序 {KVW}:V"XC"W"YC":Lbl 2:Prog "ZZ":I=V-S:J=W-T:Pol(I,J: J=J-F:K=K+Rec(I,J:AbsI<1m=>Prog"WZ":≠>Goto 2Δ M=0:{M}:M"0 NEXT"=0=>U=U+1:Goto 2:≠>U=1 ZZ %高斯法中桩子程序(4节点) Prog"XL":M=K-L:O=(P-R)÷2PQR: D=.0694318442:E=.3300094782:F=1:G=1-E:H=1-D: I=5:Lbl 1:C[I]=A+MrC[I](1÷P+OMC[I]:Dsz I:Goto 1: S=X+M(BcosD+CcosE+CcosG+BcosH: T=Y+M(BsinD+CsinE+CsinG+BsinH WY %外移点计算子程序 Lbl 1:J=90:{J}:J=F+J"<":F=J:S=S+Rec(I,J:T=T+J: Prog"ZB":I=0:{I}:I"L"≠0=>Goto 1 WZ %位置显示子程序 "KJ":K:Pause 1:J◢ ZB %坐标显示子程序 "XY":S:Pause 1:T◢ YC %异常处理子程序 U=1=>K=L:U=2Δ U=3=>K=M:U=4Δ

利用EXCEL表格计算线路中边桩坐标

利用EXCEL表格计算线路中边桩坐标 高速公路施工中,桥梁、隧道施工以及路沿石施工时对测量要求相当严格,要求总体宽度误差在1cm以内。这就要求测量工程师必须把线路边桩都计算出来,每隔10米(曲线段)或者20米(直线段)放样线路边桩,以便指导施工。如果采用常规的计算器计算,不仅繁琐、费力而且容易出现差错。本人在计算中边桩坐标时,试着利用EXCEL表格功能,编辑函数大批量计算线路中边桩坐标,达到了高效、准确的目的。由于直线部分相对简单,复曲线又特别繁琐。现就圆曲线举例计算如下: 例:在浙江省龙丽一级公路施工中,左线圆曲线起点里程为K86+966.6,曲线半径为1500m,曲线左偏,起始方位角为147°43 ′58.2″,线路中心和隧道中心偏差15cm。计算隧道中心线和衬砌边线(半径5.1m)。 1.在B4方格中输入: =$B$3+3000*SIN((A4-$A$3)/3000)*COS(2.5784244-(A4-$A$3)/3000)+ 0.15*COS(4.149220679-2*(A4-$A$3)/3000) 按回车键确认即可以计算出该里程隧道中心线X轴坐标。 2.在C4方格中输入: =$C$3+3000*SIN((A4-$A$3)/3000)*SIN(2.5784244-(A4-$A$3)/3000)+

0.15*SIN(4.149220679-2*(A4-$A$3)/3000) 按回车键确认即可以计算出该里程隧道中心线Y轴坐标。 同理在 3.D4方格中输入: =$B$3+3000*SIN((A4-$A$3)/3000)*COS(2.5784244-(A4-$A$3)/3000)+ COS(4.149220679-2*(A4-$A$3)/3000) *(0.15+5.1) 4.在E4方格中输入: =$C$3+3000*SIN((A4-$A$3)/3000)*SIN(2.5784244-(A4-$A$3)/3000)+ SIN(4.149220679-2*(A4-$A$3)/3000) *(0.15+5.1) 其中 2.5784244为方位角147°43 ′58.2″的弧度表达方式,4.149220679为起始方位角147°43 ′58.2″加90°后的弧度值。由于EXCEL自身问题无法计算度,只能以弧度方式计算,所以计算时,必须把角度转换成弧度计算。 5.使用EXCEL自动复制功能,采用鼠标拖动,就可以迅速计算出其他桩位点的隧道中心和右边桩坐标。当需要不同里程点的坐标时只需改动第一列中线路中心里程即可。由于有些版本的EXCEL无法识别K87+367这种格式,所以所有里程必须采用正常的阿拉伯数字形式。 本方法避免了手工计算,输入输出的很多麻烦,为大批量的计算公路曲线要素提供了一种有效途径,有着比计算器和测量软件更快速、准确出计算成果的优点。在本人担负测量的很多公路施工中采用此方法,达到了节省人力、准确、系统计算的功能。

坐标转换之计算公式

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

线元法万能坐标计算程序

线元法万能坐标计算程序(适用于CASIO fx-9750GⅡ计算器) 论文https://www.360docs.net/doc/1d13128381.html,/:本论文仅供学习交流使用,本站仅作合理转载,原作者可来邮要求删除论 文。 摘要:我国公路建设事业正处于一个高速发展的时期,在公路工程施工过程中,施工技术人员经常要使用全站仪、水准仪进行施工放样、高程测量,在测量过程中,手工计算速度慢,失误率高,工作效率极低。利用CASIO fx-9750GⅡ编程函数计算器强大的内存(可诸存63000个字符)和编程功能,编写各种计算程序,能够在2秒钟内计算出施工放样、桩点坐标等施工过程中的各项数据资料,同时也使我们有更多的时间去挑战更富有创造性的工作。 关键词:坐标放线线元测量程序 1、前言 本程序采用Gauss-Legendre(高斯-勒让德)五节点公式作内核,计算速度(太约2秒)适中,计算精度很高。在此之前,本人曾用过以下公式作内核:①积分公式simpson法②双重循环复化高斯2节点③高斯-勒让德3节点④求和公式复化simpson法⑤双重循环复化simpson法⑥高斯-勒让德4节点,⑦高斯-勒让德5节点,经过测试③计算最快,⑦代码稍长但计算速度只比③⑥稍慢,精度最高,可满足线元长小于1/2πD 的所有线形的精度要求。⑦作内核分别计算圆曲线长1/4πD、1/2πD、3/4πD、πD处的精度,1/4πD时偏差为0.001mm,1/2πD时偏差为0.55m m,3/4πD时偏差为31.63mm,πD时偏差为968mm,偏差按半径倍数增大,如线元长大于1/2πD(1/2圆周长)时,可将其拆分二个或多个线元单位,以确计算保精度。 2、程序特点 事先将所有的平曲线交点的线元要素诸存到计算器内,测量时只输桩号、边距等程序会自动寻找各类要素,一气呵成地完成施工测量任务,中途不需人工转换各类要素数据,本程序可诸存几百条线路的要素数据,计算时可按需选择线路编号进行测量。测量时不需查阅及携带图纸,仅一台CASIO fx-9750GⅡ编程函数计算器即可。 本程序含一个主程序:3XYF,五个子程序:GL(公式内核)、QD(线路选择)、XL(线路要素判断)、GF(坐标反算)、File 1 (要素存放的串列工作簿)。可以根据曲线段——直线、圆曲线、缓和曲线(完整或非完整型)的线元要素(起点坐标、起点里程、起点切线方位角、终点里程、起点曲率半径、止点曲率半径)及里程边距或坐标,对该线元段范围内任意里程中边桩坐标进行正反算。 3、计算公式及原理 如图:BC 间为一曲线元,曲线元上任一点的曲率随至B 点的弧长作线性变化。设起点B 的曲率为KA ,终点C 的曲率为KB ,R 为曲线半径。±表示曲线元的偏向,当曲线元左偏时取负号,当曲线元右偏时取正号,直线段以1的45次方代替(即半径无穷大)。 式中:αΑ=起始方位角l =p 点到B的距离lS=曲线总长αp=p 点切线方位角 R1=R5=0.118463442528095 ,R2 = R4 = 0.239314335249683 , R3 = 0.28444444444444 V1=1-V5= 0.046910070 ,V 2= 1-V4 = 1 0.2307653449 V3= 0.5 利用上面公式及CASIO fx-9750GⅡ编程函数计算器可编写下列计算程序。 4、程序清单 (1)、3XYF(主程序) "1→XY2→FS"?→V:V=1=>Goto 1:V=2=>Goto 2↙(选择计算功能) Lbl 1:File 1:”XLn”?→S:Prog “QD”↙(选择线路)

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

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

交点法线元法坐标计算

3、交点法、线元法坐标计算 坐标计算是根据图纸中“直线及曲线转角一览表”提供的数据计算道路中桩坐标,然后和图纸提供的“逐桩坐标表”比对,如果一样则说明输入平曲线参数输入正确,可以计算边桩坐标和其他结构物坐标了;如果中桩坐标不一样,一般是平曲线参数输入有误,需要重新检查输入,另一种结果是图纸有错,这种情况少见,但不代表没有。“直线及曲线转角一览表”和“逐桩坐标表”见附件1、附件2。 线元法是以路线的起点坐标、方位角、起终点桩号等节点元素来计算出要求的坐标;交点法是以路线的交点要素和路线的主要要素来求得坐标。 ①交点法 交点:路线的转折点,路线改变方向是相邻两直线的延长线相交的点。用JD表示, 有些图 纸上用 IP表示。 看下图: 交 点是针对曲线的(包含圆曲线和缓和曲线),一段曲线就有一个交点。交点参数有:坐标(X,Y)、交点桩号、转角值、圆曲线半径R、缓和曲线长度。 教学提供软件(轻松测量、双心软件、测量工具)交点法曲线要素输入说明: 1、QD起点坐标: 起点坐标必须在直线段上,或填写前一交点的坐标。

2、JD交点曲线要素: (1)交点桩号 (2)交点坐标(X,Y) (3)曲线半径R 始点的话,起始里程有时候需要校正,当然,并不是每个图纸给出的起点里程都需要校正,大多数图纸的起点里程已经被设计院校正过,我们输入平曲线的时候需要验证一下。如果我们按照图纸给出的起点里程输入,发现后面的交点里程都和图纸相差一个相同的值,这就表明我们输入的起点里程需要校正。 起始点里程正常输入,第二、三个交点输入完成后,检查第二个交点的切线长和交点

里程是否和图纸一样,如果切线长正确,交点里程不正确,说明起点里程需要校正,将第二个交点的里程与正确里程的差值,应用到起点里程中,从而使第二个交点里程和后面交点的里程与图纸吻合。 注意:交点法计算坐标适用的平曲线为对称或不对称缓和曲线、圆曲线。对于非普通的三单元曲线,交点法不适用。非普通的三单曲线例如下页的JD18及JD19处的平曲线, 的输入是否正确,有的图纸给的方位角数据较少,需要每隔几个线元才能检验方位角。

电子表格常用函数公式

电子表格常用函数公式 1、自动排序函数: =RANK(第1数坐标,$第1数纵坐标$横坐标:$最后数纵坐标$横坐标,升降序号1降0升) 例如:=RANK(X3,$X$3:$X$155,0) 说明:从X3 到X 155自动排序 2、多位数中间取部分连续数值: =MID(该多位数所在位置坐标,所取多位数的第一个数字的排列位数,所取数值的总个数) 例如:612730************在B4坐标位置,取中间出生年月日,共8位数 =MID(B4,7,8) =19820711 说明:B4指该数据的位置坐标,7指从第7位开始取值,8指一共取8个数字 3、若在所取的数值中间添加其他字样, 例如:612730************在B4坐标位置,取中间出生年、月、日,要求****年**月**日格式 =MID(B4,7,4)&〝年〞&MID(B4,11,2) &〝月〞& MID(B4,13,2) &〝月〞&

=1982年07月11日 说明:B4指该数据的位置坐标,7、11指开始取值的第一位数排序号,4、2指所取数值个数,引号必须是英文引号。 4、批量打印奖状。 第一步建立奖状模板:首先利用Word制作一个奖状模板并保存为“奖状.doc”,将其中班级、姓名、获奖类别先空出,确保打印输出后的格式与奖状纸相符(如图1所示)。 第二步用Excel建立获奖数据库:在Excel表格中输入获奖人以及获几等奖等相关信息并保存为“奖状数据.xls”,格式如图2所示。 第三步关联数据库与奖状:打开“奖状.doc”,依次选择视图→工具栏→邮件合并,在新出现的工具栏中选择“打开数据源”,并选择“奖状数据.xls”,打开后选择相应的工作簿,默认为sheet1,并按确定。将鼠标定位到需要插入班级的地方,单击“插入域”,在弹出的对话框中选择“班级”,并按“插入”。同样的方法完成姓名、项目、等第的插入。 第四步预览并打印:选择“查看合并数据”,然后用前后箭头就可以浏览合并数据后的效果,选择“合并到新文档”可以生成一个包含所有奖状的Word文档,这时就可以批量打印了。

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

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

坐标自动计算表格

K18+000-K20+934.86(1075m) 桩号偏距X(m)Y(m)方位角高程1808002961525.733487467.350137.72801651 1808502961529.688487470.409737.72801651 1809002961533.642487473.469237.72801651 1809502961537.597487476.528837.72801651 1810002961541.551487479.588437.72801651 1810502961545.506487482.64837.72801651 1811002961549.461487485.707537.72801651 1811502961553.415487488.767137.72801651 1812002961557.37487491.826737.72801651 1812502961561.325487494.886237.72801651 1813002961565.279487497.945837.72801651 1813502961569.234487501.005437.72801651 1814002961573.188487504.064937.72801651 1814502961577.143487507.124537.72801651 1815002961581.098487510.184137.72801651 1815502961585.052487513.243637.72801651 1816002961589.007487516.303237.72801651 1816502961592.962487519.362837.72801651 1817002961596.916487522.422437.72801651 1817502961600.871487525.481937.72801651 1818002961604.825487528.541537.72801651 1818502961608.78487531.601137.72801651 1819002961612.735487534.660637.72801651 1819502961616.689487537.720237.72801651 1820002961620.644487540.779837.72801651 1820502961624.599487543.839337.72801651 1821002961628.553487546.898937.72801651 1821502961632.508487549.958537.72801651 1822002961636.462487553.01837.72801651 1822502961640.417487556.077637.72801651 1823002961644.372487559.137237.72801651 1823502961648.326487562.196837.72801651 第 1 页,共 23 页

坐标转换模型

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

2020年excel表格计算公式大全、使用技巧

AND “与”运算,返回逻辑值,仅当有参数的结果均为逻辑“真(TRUE)”时返回逻辑“真(TRUE)”,反之返回逻辑“假(FALSE)”。条件判断 AVERAGE 求出所有参数的算术平均值。数据计算 COLUMN 显示所引用单元格的列标号值。显示位置 将多个字符文本或单元格中的数据连接在一起,显示在一个单元格中。字符CONCATENATE 合并 COUNTIF 统计某个单元格区域中符合指定条件的单元格数目。条件统计 DATE 给出指定数值的日期。显示日期 DATEDIF 计算返回两个日期参数的差值。计算天数 DAY 计算参数中指定日期或引用单元格中的日期天数。计算天数 DCOUNT 返回数据库或列表的列中满足指定条件并且包含数字的单元格数目。条件统计 FREQUENCY 以一列垂直数组返回某个区域中数据的频率分布。概率计算 IF 根据对指定条件的逻辑判断的真假结果,返回相对应条件触发的计算结果。条件计算 INDEX 返回列表或数组中的元素值,此元素由行序号和列序号的索引值进行确定。数据定位 INT 将数值向下取整为最接近的整数。数据计算 ISERROR 用于测试函数式返回的数值是否有错。如果有错,该函数返回TRUE,反之返回FALSE。逻辑判断 LEFT 从一个文本字符串的第一个字符开始,截取指定数目的字符。截取数据 LEN 统计文本字符串中字符数目。字符统计 MATCH 返回在指定方式下与指定数值匹配的数组中元素的相应位置。匹配位置 MAX 求出一组数中的最大值。数据计算 MID 从一个文本字符串的指定位置开始,截取指定数目的字符。字符截取 MIN 求出一组数中的最小值。数据计算

不同坐标系之间的变换

不同坐标系之间的变换 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;

Excel表格自动计算技巧

Excel表格自动计算技巧 一、显示单元格例有计算式的结果的设置方法 首先:插入-名称-定义在弹出的对话框“当前工作薄的名称”中输入: X 或“结果”的自定义名称,再在“引用位置”处粘贴=EVALUATE(SUBSTITUT E(SUBSTITUTE(SUBSTITUTE(SUBSTITUTE(计算!$C$1,"[","("),"]",")"),"×","*"),"÷","/"))公式(注意要有“=”号。再对公式中“计算!$C$1”选择上,然后再到需要做公式的单元格中点击即可。如要相对引用,则要删除$字符。 已经OK,你在 C1输入表达式比如 15+5×3 ,在D1中输入 =x 看看(应该是30)。 二、如何在Excel中输入计算式后另一单元显示计算结果 菜单--插入--名称--定义:输入AA(任意取名),在下面输入公式: =EVALUATE($A$1) 然后在B1单元格输入公式:=AA 但本式不能识别如:[、×、÷等符号进行计算,使用第一种较好。 三、如何在Excel中编写自定义函数,象在表格中调用SUM()一样?

第一种情况: 单元格 A1 = 2;单元格 B1 = 2;单元格 C1 = 2 单元格 D1=(A1+B1)×C1 显示结果为 8 ; 那么如何才能在单元格 E1中显示(2+2)×2的计算表达式,并且建立关联,当单元格D1变成 =(A1+BA)^C1 计算式后,显示结果为16;那么E1也就自动显示为 (2+2)^2的计算表达式,也就是说随着单元格D1的计算公式变化,单元格E1显示的计算表达式也随之变化。 第二种情况: 是当计算的单元格任意变化时,怎么办?比如说计算式E1=A1+B1+C1+D1,也有D2=A2+B2+C2,还有F3=A1+B1+C1+D1+E1时。如何将计算式变为计算表达式。即计算式可以在任意一个单元格,计算公式所引用数据的单元格可以任意变化,需要显示计算表达式的单元格也是任意的。 针对第一种问题,分两步做: 1,自定义一个函数 Public Function K(x As Range) K = Right(x.Formula, Len(x.Formula) - 1) End Function 自定义函数操作过程如下:

EXCEL角度与坐标计算

主题2 EXCEL角度与坐标计算 1.3 Excel角度与坐标计算 电子表格Excel是微软的Office办公软件的重要组件,可用于编制表格、进行各种数值的批量计算,还可以利用内置的VBA程序进行二次开发自动计算,是工程中使用最广泛的软件之一。 Excel在工程计算应用中的难点是角度和三角函数的应用。 1.3.1 Excel三角函数计算 在Excel中,三角函数中的角度以及反三角函数计算出的角度结果,格式均为弧度。由于工程实际中使用度(以及60进制的度分秒),而Excel三角函数中使用弧度,两者之间必须进行转换。 除三角函数和反三角函数外,Excel中关于角度的函数,有三个比较重要: (1)圆周率常数:PI(); (2)角度转换函数(弧度转换成十进制度):DEGREES(); (3)角度转换函数(十进制度转换成弧度):RADIANS()。 Excel中角度和三角函数的基本计算见表1-1,使用时需要头脑清醒,灵活运用。 表1-1 Excel中角度和三角函数计算 输入公式计算结果说明 =PI() 3.141592654 圆周率常数,半圆对应的弧度,函数后面括号内无参数,但括号不能省略 =PI()/6 0.523598776 =RADIANS(30) 0.523598776 十进制度转换为弧度 =DEGREES(PI()/6) 30 弧度转换为十进制度 =SIN(PI()/6) 0.5 三角函数使用角度参数为弧度 =SIN(RADIANS(30)) 0.5 =ASIN(0.5) 0.523598776 反三角函数计算结果为弧度 =DEGREES(ASIN(0.5)) 30 计算结果弧度转换为十进制度 【操作1-8】某条支导线如图1-24所示,已知B点坐标及AB边的坐标方位角 AB a,观

坐标系转换步骤以及公式

一、各坐标系下椭球参数 二、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 。

相关文档
最新文档