由坐标反算出路线里程桩号在工程测量中的应用
由坐标反算出路线里程桩号在工程测量中的应用
摘要：文章系统地介绍了由已知坐标，反向推算出该坐标位于公路路线的里程桩号及距中距离的基本概念和方法，重点介绍利用卡西欧fx-4800计算器编程，由已知坐标反向推算出公路路线中常见的直线段、圆曲线段的里程桩号方法及其在工程测量中的应用，对缩短工程测量时间，提高测量效率起到重要作用。
关键词：坐标反向推算；公路路线里程；编程；工程测量；
引言
在公路施工测量中，我们常用的测设方法是通过设计文件、图纸中的直线、曲线及转角表，来计算出路线的中桩、边桩坐标，再利用全站仪的放样功能，将中桩及边桩坐标在实地位置测定出来。那么我们能不能由某一已知坐标，通过卡西欧fx-4800计算器编程计算，就能知道该坐标的准确的路线里程桩号以及其距中线的距离，甚至其设计高程，这对我们提高施工测量的效率有很大的帮助。本文将阐述由已知坐标反向推算出公路路线中较为常见的直线段、圆曲线段里程桩号的基本概念和方法，以及介绍利用卡西欧fx-4800计算器编程，快速反算出已知坐标的里程桩号及距中距离。
1 直线段坐标反向推算路线里程桩号
如图1，设直线段路线起始点O的里程桩号为K1，中桩平面坐标为（X0，Y0），路线方位角为α，P点平面坐标为（X，Y），求出P点位于该直线段的里程桩号，以及距中距离，计算步骤为：
①计算直线起点O至P点的直线距离，计算公式为:
②反算直线起点O至P点的坐标方位角，计算公式为:
③计算直线0P与路线的夹角，计算公式为:
④计算P点垂直投影于路线的垂直距离，即P点到路线中心的距离，计算公式为
坐标正算反算公式讲解(借鉴材料)
一 方位角： 在高斯直角坐标系中，由坐标纵轴方向的北端起，顺时针量到直线间的夹角，称为该直线的坐标方位角，常简称方位角，用a 表示。 1、第一象限的方位角 Y X 第一象限第二象限 第三象限 第四象限 o A a 图1 2、第二象限的方位角 Y X 第一象限 第二象限第三象限 第四象限 o A a 图2
3、第三象限的方位角 Y X 第一象限 第二象限 第三象限 第四象限 o A a 图3 4、第四象限的方位角 Y X 第一象限 第二象限 第三象限 第四象限 o A a 图4 方位角计算公式：
x =a -1 tan A Y O Y -A X O X - 方位角的计算器计算程序：Pol(X A -X O ,Y A -Y O ) 直线OA 方位角度值赋予给计算器的字母J ，0≤J ＜360。 直线段OA 的距离值赋予给计算器的字母I,I ＞0 直线OA 与直线AO 的方位角关系： 1、 当直线OA 的方位角≤180°时，其反方位角等于a+180°。 2、 当直线 OA 的方位角＞180°时，其反方位角等于a-180°。 二 方位角的推算 （一）几个基本公式 1、坐标方位角的推算
或： 注意：若计算出的方位角>360°，则减去360°；若为负值，则加上360°。 例题:方位角的推算 已知：α12=30°，各观测角β如图，求各边坐标方位角α23、α34、α 45 、α51。 13 图5
解：α23= α12-β2+180°=30°-130°+180°=80° α34= α23-β3+180°=80°-65°+180°=195° α45=α34-β4+180°=195°-128°+180°=247° α51=α45-β5+180°=247°-122°+180°=305° α12=α51-β1+180°=305°-95°+180°=30°（检查） 三坐标正算 一、直线段的坐标计算 o B D A C E a a p 图6 设起点O的坐标（X O,Y O）,直线OP的方位角为F op，求A、C、E点的坐标 1、设直线段OA长度为L,则A点坐标为 X A=X O+L×Cos(F op)
高斯投影正反算公式 新
高斯投影坐标正反算 一、相关概念 大地坐标系由大地基准面和地图投影确定，由地图投影到特定椭圆柱面后在南北两极剪开展开而成，是对地球表面的逼近，各国或地区有各自的大地基准面，我国目前主要采用的基准面为：基准面，为GPS基准面，17届国际大地测量协会上推荐，椭圆柱长半轴a=6378137m，短半轴b=； 2.西安80坐标系，1975年国际大地测量协会上推荐，椭圆柱长半轴a=6378140m，短半轴b=; 3.北京54坐标系，参照前苏联克拉索夫斯基椭球体建立,椭圆柱长半轴a=6378245m, 短半轴b=; 通常所说的高斯投影有三种，即投影后： a)角度不变（正角投影），投影后经线和纬线仍然垂直； b)长度不变； c)面积不变； 大地坐标一般采用高斯正角投影，即在地球球心放一点光源，地图投影到过与中央经线相切的椭圆柱面上而成；可分带投影，按中央经线经度值分带，有每6度一带或每3度一带两种(起始带中央经线经度为均为3度，即：6度带1带位置0-6度，3度带1带位置度)，即所谓的高斯-克吕格投影。
图表11高斯投影和分带 地球某点经度（L）为过该点和地球自转轴的半圆与子午线所在半圆夹角，东半球为东经，西半球为西经；地球某点纬度（B）为所在水平面法线与赤道圆面的线面角。 正算是已知大地坐标（L，B），求解高斯平面坐标（X，Y）,为确保Y值为正，Y增加500公里；反算则是由高斯平面坐标（X，Y）求解大地坐标（L，B）。 二、计算模型： 地球椭球面由椭圆绕地球自转轴旋转180度而成。 图表 1 椭圆 椭圆长半轴a，椭圆短半轴b, 椭圆方程：
(1) 图表2椭球面 椭球面方程： y2 a2+ x2 b2 + z2 a2 =1 /*************************************** 与网上充斥的将函数关系先展开为泰勒级数，再依据投影规则确定各参数不同，本文直接依据空间立体三角函数关系得出结果。 *****/ （一）正算 由图表1，
关于坐标反算里程的超常规方法
关于坐标反算里程的超常规方法 (仅适用于隧道工程) 对于曲线隧道超开挖线放样测量,大家都觉得用坐标反算里程时计算的时间长,有的可能在三秒以内,有的可能会超过五秒钟,那让人期待的滋味的确让人难以忍受,如果不慎把坐标数据输入错误,更会让人抓狂(呵呵,因为必须得输正确坐标数据,错了要返工.),这样的无奈我也深有感触啊.有没有什么好方法让计算时间缩短呢?肯定有,下面我会详细介绍. 不管什么样式的曲线,它在某些情况下是可以看成直线的,(比如说 1.5米,2米这样的曲线线段),我说的这种方法就是把曲线当成直线来计算,这样就少去了很多计算步骤,而且能直接计算出结果.在现实施工中,我们的开挖面都是倾斜的,这个倾斜面就存在一个里程差,我们就可以把这个里程差当做直线来计算,我们只需取两个点就可以了,一个是凹进去最深的点,一个是最突出的点,我们架设好仪器后就先测算出这两点的里程,然后就把我们计算范围界定在这两个里程之间,把这两个里程点之间的当着一条直线,而且这条直线一般都很短,很少有超过2米的,然后我们在放开挖线的时候测出来的坐标都通过这条直线来反算它对应的真实里程,计算出测量点的偏移值. 接下来我们再说说精度的问题,我曾在CAD中画过这样的图,里程差取2米,半径取值为50米.隧道半幅宽度按8米计算,垂直于这条直线的距离和半径方向的距离误差也就只有0.0077米,半径50米2米长的弧线外矢距也就0.01米(现实中设计的曲线半径很少有小于50米的),完全能满足要求.且我们实际的开挖尺寸也是按直线掘进的,进尺完全大于我们所取的这两个里程差.所以这个误差完全是可以忽略不记的. 以上这种方法,我已经在实际当中用过,真的不错哟.希望同行们留意,也许能对你有所帮助,本人QQ468076885,欢迎同行灌水.
5800简单全线坐标计算程序
5800全线任意坐标计算程序 1. 正算主程序(ZHCX) （不运行） 8→DimZ 1÷P→Z[4 ]：(P-R)÷(2HPR)→D: 180÷π→E “Z=”?Z:”YJJ=”?A:Abs(S-O)→W 0.26→Z[1 ]: 0.74→B: 0.02→K: 0.82→Z[3 ]: 1-Z[3 ]→F:1-K→Z[2 ] U+W(Z[1 ]cos(G+QEKW(Z[4 ]+KWD))+Bcos(G+Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+Bcos(G+QEFW (Z[4 ]+FWD))+ Z[1 ]cos(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→X: V+W(Z[1 ] sin (G+QEKW(Z[4 ]+KWD))+B sin(G+ Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+B sin(G+QEFW (Z[4 ]+FWD))+ Z[1 ] sin(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→Y： G+QEW(Z[4 ]+WD)→F：X+Zcos(F+A)→X：Y+Zsin(F+A)→Y:If F≧360:Then F-360→F:IfEnd ”X=”:X→X◢ ”Y=”:Y→Y◢ If F﹤0:Then F+360→F:IfEnd ”QX FWJ=”:F▼DMS◢ “C=1=>XX: C=2=>XZ”: ”C=”?C: ”QHJU=”?L: If C=1：Then Goto 1:Else Goto 2: IfEnd 可以计算斜交斜做或斜交正做的桥涵坐标 Lbi 1 X+L cos(F)→X：Y+Lsin(F)→Y: Goto 3 Lbi 2 X+L cos(F+A-90)→X：Y+Lsin(F+A-90)→Y: Goto 3 Lbi 3 “QH-X=”: X →X◢ “QH-Y=”: Y →Y◢ Prog “FY” 2 . 参数子程序（直接运行） M(主线) 一条线路一个名称 “S=”?S If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd … … If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd Lbi 1 Prog “ZBJS” 3. 放样程序（FY）（不运行） “X0=”?M:“Y0=”？N Pol((X-M, Y-N)
坐标反算
在现场工作中，以往我们都是已知某点的里程及边距，来计算出该点的坐标，但有时我们如果能在测得某点坐标后，计算出该点的里程和距线路中心的距离(在这里我姑且称之为坐标反算)的话，将会帮我们大大减轻野外工作量，提高我们的工作效率。例如：路基填了几层后要精确检查一下路基是否够宽，那么按照我们以往的做法，就是要先将线路中心线放出来，然后用尺拉一下路基宽度，与其在此高程的设计宽度作比较，这样做对高填方而言极不方便。或者是先按所测高程,计算的宽度放出路基边桩，再与所填边线作比较。以上两种方法现场工作量都比较大。较为简便的方法是，我们可以测一下已填路基边线上任一点的三维坐标，然后将其反算求出该点的里程，及其距中线的距离(即所填宽度),由计算出的里程，可算出该里程的路面设计高程，再有所测高程，可计算出该点的设计宽度，两宽度作比较即可。同样在桥面铺装施工时，我们也无须再像以往那样，先放出某点再测其高程，然后与设计高程比较计，算出该点铺装厚度，而可以沿桥面外边线随意布点，测其三维坐标，计算出其里程及到中心线的距离，便可由其里程及距中心距离，计算出该点的设计高程，与其测得高程作比较得出应铺厚度。这样便大大减轻了外业工作强度(由放出点后再测其高程，变为测任意点高程)，而内业计算量与常规相当。另外在临时增加桥涵时，也常用到此方法来计算变更桥涵的中心里程(斜交或正交均可).如目前我标段就存在很多临时变更涵洞,按以往我们的方法是先估计该处大概里程,然后放出所估计里程的中心桩,再用皮尺量出所要增加涵洞处与该中心桩的距离,以此来推算出涵洞的中心里程,这一过程即繁琐又不准确。而目前我们采用的方法是用全站仪测得跨路基现有水沟两端的沟底坐标,计算出其与路基的夹角,按所测坐标及此夹角就可以准确、快速地反算出水沟中心所对应的线路中心里程了。我们在日常测量工作中的很多方面，也会用到这一方法来减轻野外工作量。在目前我标段的S334分离式立交桥的架设过程中,也同样用到了此方法.支座安装好后,对支座中心位置检及高程查无误后开始梁板架设,但是尽管测量控制放样符合规范要求,可是因为其它方面的各种原因可能会使梁板出现偏位高程也可能会出现偏差,那么对现在这种问题该如何检查呢?其实方法是一样的,首先我们可以用全站仪测得架设好后梁的边板外边缘任一点的三维坐标,由此坐标反算出该点所对应的中心里程和距中心的距离,就可以和设计图纸上的距离作比较来检查其是否存在偏位,该点的设计高程也可以由反算所得的中心里程和距中心的边距算出,与所测得的实际高程作一下比较也就可以了.那么通过以上讨论问题归结到了一点，那就是如何在测得任一点坐标后，计算出其所对应的线路中心里程，及其到线路中心的距离（或是斜交的长度）呢？解决此类问题，对目前一些测量软件来说早已不成问题，但是在现场工作中我们用的更多、更方便的还是计算器，那么能否用我们常用的4800或5800计算器编程，来计算此类问
线元法万能坐标计算程序
线元法万能坐标计算程序(适用于CASIO fx-9750GⅡ计算器) 论文https://www.360docs.net/doc/118585277.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”↙（选择线路）
坐标正反算计算公式
坐标正反算公式
一、GPS数据处理相关术语 1、三维无约束平差 三维无约束平差是以基线解算所得到的三维静态基线向量为观测值，待定参数主要为GPS 网中点的坐标；同时，利用基线解算时随基线向量一同输出的基线向量的方差阵，形成平差的随机模型，最终形成平差完整的数学模型。随后对所形成的数学模型进行求解，根据平差结果来确定观测值中是否存在粗差，数学模型是否有需要改进的部分，若存在问题，则采用相应的方法进行处理并重新进行求解；若未发现问题，则输出最终结果，并进行后续的数据处理。 2、三维约束平差 三维约束平差是以基线解算所得到的三维静态基线向量为观测值，在平差过程中引入会使GPS 网的尺度、方向和位置发生变化的外部起算数据，从而实现GPS 网成果由基线解算时GPS 卫星星历所采用的参照系（WGS84 ）到特定参照系的转换，得到在特定参照系下的经过用户约束条件约束的点三维空间坐标。 二、南方GPS数据处理软件的平差方式
三维约束平差是指在基线解算后，WGS84坐标系下的三维平差，在三维平差中是不需要当地平面直角坐标系下的已知点坐标，当需要用到WGS84经纬度或空间直角坐标的用户可加载已知点的WGS84空间坐标（如果只有经纬度时，可采用COORD4.1软件进行转换，本站免费提供）进行三维约束平差，即可得到与已知点相匹配的WGS84坐标。 一般情况下，在“已知点坐标录入”窗口中，我们都没有输入WGS8坐标，而只输入当地坐标系下的已知坐标，此时GPS处理软件会自动识取一个坐标点的WGS84坐标进行约束平差。如下图：
如果在某些控制测量中,需要得到精确的WGS84经纬度或空间坐标时,让系统自动识取显然是不行的,此时我们只要为参与平差的已知点的WGS84空间坐标输入后再进行三维平差即可 在这里,我们加入了两个已知点的WGS84空间坐标,三维平差后,列表中会显示两个"固定"字样的点,说明,在进行三维平差中,我们把这两个点做为起算点,进行平差别的未知点。
线元法线路坐标正反算程序
经苦心钻研，奋战多日，终于编写出了代码短，速度快，精度高，功能全的线路坐标正反算程序，欢迎试用并提出宝贵意见。 功能简介及特点： 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Δ
坐标正反算定义及公式
坐标正反算定义及公式 Corporation standardization office #QS8QHH-HHGX8Q8-GNHHJ8
第六章→第三节→导线测量内业计算 导线计算的目的是要计算出导线点的坐标，计算导线测量的精度是否满足要求。首先要查实起算点的坐标、起始边的方位角，校核外业观测资料，确保外业资料的计算正确、合格无误。 一、坐标正算与坐标反算 1、坐标正算 已知点的坐标、边的方位角、两点间的水平距离，计算待定点的坐标，称为坐标正算。如图6-6 所示，点的坐标可由下式计算： 式中、为两导线点坐标之差，称为坐标增量，即： 【例题6-1】已知点A坐标，=1000、=1000、方位角=35°17＇36.5"，两点水平距离=200.416，计算点的坐标？
35o17＇36.5"=1163.580 35o17＇36.5"=1115.793 2、坐标反算 已知两点的坐标，计算两点的水平距离与坐标方位角，称为坐标反算。可知，由下式计算水平距离与坐标方位角。 （6-3） （6-4） 式中反正切函数的值域是-90°～+90°，而坐标方位角为0°～360°，因此坐标方位角的值，可根据、的正负号所在象限，将反正切角值换算为坐标方位角。 【例题6-2】=3712232.528、=523620.436、=3712227.860、=523611.598，计算坐标方位角计算坐标方位角 、水平距离。
=62°09＇29.4"+180°=242°09＇29.4" 注意：一直线有两个方向，存在两个方位角，式中：、的计算是过A点坐标纵轴至直线的坐标方位角，若所求坐标方位角为，则应是A点坐标减点坐标。 坐标正算与反算，可以利用普通科学电子计算器的极坐标和直角坐标相互转换功能计算，普通科学电子计算器的类型比较多，操作方法不相同，下面介绍一种方法。 【例题6-3】坐标反算，已知=2365.16、=1181.77、 =1771.03、=1719.24，试计算坐标方位角、水平距离。 键入1771.03-2365.16按等号键［=］等于纵坐标增量，按储存键［］， 键入1719.24-1181.77按等号键［=］等于横坐标增量，按［］键输入，按［］显示横坐标增量，按［］键输入，按第二功能键［2ndF］，再按［］键，屏显为距离，再按［］键，屏显为方位角。 【例题6-4】坐标正算，已知坐标方位角=294°42＇51"， =200.40，试计算纵坐标增量横坐标增量。
坐标正反算程序计算器
一、 Lbl 3:"1→ZS,2→FS"?Q Q=1=>Goto 1:Q=2=>Goto 2 Lbl 1:"CE:X"?M:"CE:Y"?F:"JL"?L:"FWJ"?A:Rec(L,A):M+I→C:F+J→D Cls "X=":Locate 3,1,C:"Y=":Locate 3,2,D◢ Goto 3 Lbl 2::"CE:X"?G:"CE:Y"?H:"(HOU)FY:X"?N:"(HOU)FY:Y"?E Pol(N-G,E-H) If J<0:Then J+360→Y:Else J→Y:IfEnd Cls "FY JL=":Locate 10,1,I:"FY FWJ=":Y◆DMS◢ Goto 3 进入程序运行如下: 1→ZS,2→FS? 输入1为正算,2为反算. 以输入1为例: CE:X? 测站点X(5796.717) CE:Y? 测站点Y(5212.569) JL? 仪器测得的距离(321.889) FWJ? 仪器测得的方位角(193-41-07) 得到:X=5483.966 Y=5136.414 再按EXE,输2为例: CE:X? 测站点X(5796.717) CE:Y? 测站点Y(5212.569) (HOU)FY:X? 后视或放样的X(5483.966) (HOU)FY:Y? 后视或放样的Y(5136.414) 得到:FY JL=321.889 FY FWJ=193-41-6.79 二、 Deg : Fix 3 : “XZ→0:YZ→1”？A : If A = 1: Then Goto 1 : IfEnd ↙ If A = 0 : Then “BS→0:XY→1:AND→2:DK→3:L(I)→4 ”？O : IfEnd ↙ If O = 4: Then Goto 1 : IfEnd ↙ If O = 3: Then Prog “F.2 ”: If X= 0 : Then Goto 1 : IfEnd : IfEnd ↙ If O≠1: Then “X1 ”？X : “Y1”？Y : X→Z［11］: Y→Z［12］: “X2 ”？P : “Y2”？Q : Pol( P－X , Q－Y) : If J﹤0 : Then J + 360→J : IfEnd : Cls : “S12= ”: Locate 6 ,1, I : “B12= ”: J ?DMS◣
fx-4800P缓和曲线和圆曲线坐标正反算程序(正确版)
缓和曲线和圆曲线坐标正反算程序主程序“TYQXJS” Lb1 0↙→（EXE） {NUVOGHPRQ}:“1.SZ=>XY”: “2.xy=>SZ”: N:U“QDX”:V“QDY”:O“QDLC”:G“FWJ”: H“LS”:P“RO”:R“RN”:Q“ZP=-1,YP=+1,ZZ=0” :C=1÷P:D=(P-R)÷(2HPR):E=180÷π:N=1=>Goto 1:≠=>Goto 2◣↙Lb1 1:{SZ}:SZ:W= Abs(S-O): Prog“1”: X“XS”=X◢ Y“YS”=Y◢ F“FS”=F-90◢ Goto 3↙ Lb1 2:{XY}:XY:I=X:J=Y: Prog“2”:S“S”=O+W◢ Z“Z”=Z◢ Goto 3↙ Lb1 3↙ {DE}:E“QX-JJ,Z-1,Y+1”:D“BZ-JL”↙ F=F+E↙ X=X+D Cos F◢ Y=Y+D Sin F◢ {DE}:D“BZ-JJ”:E“JJ,Z-1,Y+1”↙ F=F+E↙
X=X+D Cos F◢ Y=Y+D Sin F◢ Goto 0 子程序1：“1” A=0.1739274226:B=0.3260725774:K=0.0694318442:L=0.3300094782:F=1 -L:M=1-K:X=U+W(ACos(G+QEKW(C+KWD))+BCos(G+QELW(C+LWD))+BCos(G+Q EFW(C+FWD))+ACos(G+QEMW(C+MWD))):Y=V+W(ASin(G+QEKW(C+KWD))+BSin (G+QELW(C+LWD))+BSin(G+QEFW(C+FWD))+ASin(G+QEMW(C+MWD))):F=G+QE W(C+WD)+90:X=X+Z Cos F:Y=Y+Z Sin F 子程序2: “2” T=G-90:W=Abs((Y-V)Cos T-(X-U)Sin T):Z=0: Lb1 0:Prog“1”:L=T+QEW(C+WD):Z=(J-Y)Cos L-(I-X)Sin L:Abs Z<1E-6=>Goto 1: ≠=>W=W+Z:Goto 0◣↙(E为：4800P键盘的EXE键) Lb1 1:Z=0Prog“1”:Z=(J-Y)÷Sin F◣ 注：第一缓和曲线起点半径输入无穷大（10 45），终点输入圆曲线半径；第二缓和曲线起点半径输入圆曲线半径，终点半径输入无穷大（10 45）；圆曲线输入给出的起点和终点半径；直线段则都输入无穷大（10 45）。 其中起点切线方位角用此程序计算。 ◣(代替空心) 坐标正算 程序“ZBZS” Lb1 0↙ X“HSX”:Y“HSY”:U“CZX”:V“CZY”↙ X-U≥0=>Goto 1: ≠=>Goto 2◣↙
坐标反算正算计算公式
坐标反算正算计算公式 一、坐标正算 根据A点的坐标X A、Y A和直线AB的水平距离D AB与坐标方位角O AB，推算B点的坐标X B、Y B，为坐标正算，其计算公式为： X B = X A + AX AB Y B = X A + AY AB(1-18 ) 二式中，AX AB与AY AB分别称为A?B的纵、横坐标增量，其计算公式为： AX AB = X B—X A = D AB COS O AB AY AB = Y B—Y A = D AB sin O AB(1-19 ) 注意，AX AB和AY AB均有正、负，其符号取决于直线AB的坐标方位角所在的象限。 二、坐标反算 根据A、B两点的坐标X A、Y A和X B、Y B，推算直线AB的水平距离D AB与坐标方位角 OCAB ， 为坐标反算。其计算公式为： (1-20 ) 注意，由(1-20 )式计算OCAB时往往得到的是象限角的数值，必须先根据AX AB、AY AB的正、负号，确定直线AB所在的象限，再将象限角换算为坐标方位角。 三角函数内容规律 三角函数看似很多，很复杂，但只要掌握了三角函数的本质及内部规律就会发现 三角函数各个公式之间有强大的联系。而掌握三角函数的内部规律及本质也是学好三 角函数的关键所在. 1、三角函数本质： 三角函数的本质来源于定义，如右图： 根据右图，有 sin 0 =y/ R; cos 0 =x/R; tan 0 =y/x; cot 0 =x/y。 深刻理解了这一点，下面所有的三角公式都可以从这里出发推导出来，比如以推导 si n( A+B) = si nAcosB+cosAs inB 为例： 推导： 首先画单位圆交X轴于C，D，在单位圆上有任意A，B点。角AOD为a，BO D为B，旋转AOB使0B与0D重合，形成新A'OD。 A(cos a ,sin a ),B(cos 3 ,sin 3 ),A'(cos( - BM,sin( 诩)) OA'=OA=OB=OD=1,D(1,0) [cos( a- 3 >1]A2+[sin( a- 3 )]A2=(cos a cos 3 )A2+(sin a-sin 3 )A2 和差化积及积化和差用还原法结合上面公式可推出(换(a+b)/2与(a-b)/2 ) [1] (1-21 )
隧道坐标正反算
隧道测量---快速坐标正反算程序(4800-4850计算器） 2008-11-23 20:47:00 点击：212 ZB （坐标正算） 1. Lb1 1 2. ｛L｝ 3.SMNARCL 4. K=L-S 5. R=0=>X“X”=M+KcosA▲ 6. Y“Y”=N+KsinA▲ 7. Z“FW”=A▲ 8.Goto 4△R≠0=>Abs K≤C=>Goto 2△R≠0=>K＞C=>Goto 3△ R≠0=>K<0=>Abs K＞C=>Goto 1△ 9. Lb1 2 10. O=90K2÷（πRC） 11. G=K3÷（6RC）-K?7÷（336（RC）3＋K?11÷（42240（RC）?5-K?15÷9676800（RC）?7）＋K?19÷（3530096640（RC）?9） 12. F=K－K?5÷（（40（RC）2）－K?9÷（3456（RC）?4）＋K?13÷（599040（RC）?6）－K?17÷（40320×4532（RC）?8） 13. J=0=>X“X”=M＋GsinA＋cosA▲ 14. Y“Y”=N＋FsinA- cosA▲ 15. Z“FW”=A-O▲ Goto 4△ 16．J≠0=>X“X”=M＋F cosA－GsinA▲ 17．Y“Y”=N＋G cosA＋F sinA▲ 18．Z“FW”=A＋O▲ Goto 4△ 19．LbI 3 20．B=90（2K-C）÷（лR） 21.U=C÷2－C3÷（240R2）＋R sinB＋C?5÷（34560R?4）-C?7÷（599040 R?6） 22．V=C2÷（24R）＋R（1－cosB）－C?4÷（2688R3）＋C?6÷（42240 R?5） 23．J=0=>X“X”=M＋UcosA＋VsinA▲ Y“Y”=N＋U sinA－V cosA▲ Z“FW”=A-B▲ Goto 4△ 24．J≠0=>X“X”=M＋U cosA－VsinA▲ Y“Y”=N＋U sinA＋V cosA▲ Z“FW”=A＋B▲ Goto 4△ 25．Lb1 4 26．｛D｝：｛T｝ 27．D=0 => Goto 1△ 28．D≠0=>X“LX”=X＋Dcos（T＋Z）▲
线路上点坐标反算其里程桩号及中桩距
一、ZUOBIAO—ZH（由线路上任意点坐标计算相应的桩号、中桩距[偏左或偏右的距离]） 主程序： AbsJ“Fixmc?”=7?Fixm：Goto1△{ABFHGRQ}：A“X-JD”：B“Y-JD”：G“A0”：G＜0?Q=1△G＞0?Q=-1△F=AbsG:Prog“S3”：G=F：R“R”：O“Gra0(ZY-JD)”:F=O：Prog“S3”：F=F：H“L-ZY”：S=R÷COS（G÷2）：M“X（YXIAN）”=A+Scos（F+180°+Q×（90°-G÷2））◢N“Y（YIAN）”=B+Ssin（F+180°+Q×（90°-G÷2））◢L=2πR÷360°×G：F=F- Q×90°+180°：F＜0?F=F+360°△F＞360° ?F=F-360°：≠?F=F△Lbl1:{CD}:C“X”：D“Y”：POL（M-C，N-D）：T=J： T=Ans+180°：P=I：U=T-F：-Q×U≦0?Prog“T”：Goto1△-Q×U＞0?Prog“U”：Goto1△ (主程序完) 子程序： 1、V Pol(A—C，B—D)：T=J：T=Ans+180°：P=I：U=T-（F+90°Q）：E“LI”=H+Rtan(G ÷2)-PcosU◢K“B”=-Psin(u)◢ 2、W Pol(A—C，B—D)：T=J：T=Ans+180°：P=I：U=T-（F-GQ-90°Q）：E“LI”=H+L+（PcosU-Rtan(G÷2)）◢K“B”=Psin(u)◢ 3、T Abs U≤180°?Prog“V”△A bs U＞180°?U=T+360°-F：U＞G?Prog“W”：≠?Prog“X”△ 4、U AbsU＞180°?Prog“V”：≠?AbsU≤G?Prog“X”⊿AbsU＞G?Prog“W”⊿5、X E“LI”=H-Q（2πR÷360°×U）◢K“B”=Q（P-R）◢ 说明： 一、本程序用于计算直线与圆曲线组合的线路上由点坐标计算相应桩号及其从中桩左、右偏移量。 二、输入：X-JD、Y-JD：交点坐标 A0：转向角（以十进制的形式，即13°25′30″按13.2530 的形式输入,左拐曲线取“-”值，右拐曲线取“+”值） R：半径 Gra0(ZY-JD):ZY点至JD方向坐标方位角,输入形式同A0(转向角) L-ZY:ZY点桩号 X，Y：所要计算点的坐标 三、显示：X（YXIAN）、Y（YXIAN）：圆心坐标 LI：计算所得桩号 B：所计算点至中桩的距离，负值为左偏量，正值为右偏量，0为中桩
坐标计算方法
已知：①缓和曲线上任一点离ZH点的长度：l ②圆曲线的半径：R ③缓和曲线的长度：l0 ④转向角系数：K(1或－1) ⑤过ZH点的切线方位角：α ⑥点ZH的坐标：xZ，yZ 计算过程： 说明：当曲线为左转向时，K=1，为右转向时，K=-1，公式中n的取值如下： 当计算第二缓和曲线上的点坐标时，则： l为到点HZ的长度 α为过点HZ的切线方位角再加上180° K值与计算第一缓和曲线时相反 xZ，yZ为点HZ的坐标 切线角计算公式：
已知：①圆曲线上任一点离ZH点的长度：l ②圆曲线的半径：R ③缓和曲线的长度：l0 ④转向角系数：K(1或－1) ⑤过ZH点的切线方位角：α ⑥点ZH的坐标：xZ，yZ 计算过程： 说明：当曲线为左转向时，K=1，为右转向时，K=-1，公式中n的取值如下： 当只知道HZ点的坐标时，则： l为到点HZ的长度 α为过点HZ的切线方位角再加上180° K值与知道ZH点坐标时相反
xZ，yZ为点HZ的坐标 三、曲线要素计算公式 公式中各符号说明： l——任意点到起点的曲线长度（或缓曲上任意点到缓曲起点的长度）l1——第一缓和曲线长度 l2——第二缓和曲线长度 l0——对应的缓和曲线长度 R——圆曲线半径 R1——曲线起点处的半径 R2——曲线终点处的半径
P1——曲线起点处的曲率 P2——曲线终点处的曲率 α——曲线转角值 四、竖曲线上高程计算 已知：①第一坡度：i1(上坡为“＋”，下坡为“－”) ②第二坡度：i2(上坡为“＋”，下坡为“－”) ③变坡点桩号：SZ ④变坡点高程：HZ ⑤竖曲线的切线长度：T ⑥待求点桩号：S 计算过程： 五、超高缓和过渡段的横坡计算
计算坐标与坐标方位角基本公式
二 计算坐标与坐标方位角的基本公式 控制测量的主要目的是通过测量和计算求出控制点的坐标，控制点的坐标是根据边长及方位角计算出来的。下面介绍计算坐标与坐标方位角的基本公式，这些公式是矿山测量工中最基本最常用的公式。 一、坐标正算和坐标反算公式 1．坐标正算 根据已知点的坐标和已知点到待定点的坐标方位角、边长计算待定点的坐标，这种计算在测量中称为坐标正算。 如图5—5所示，已知A 点的坐标为A x 、A y ，A 到B 的边长和坐标方位角分别为AB S 和AB α，则待定点B 的坐标为 AB A B AB A B y y y x x x ?+=?+= ｝ （5—1） 式中 AB x ? 、AB y ?——坐标增量。 由图5—5可知 AB AB AB AB AB AB S y S x ααsin cos =?=? ｝ （5—2） 式中 AB S ——水平边长； AB α——坐标方位角。 将式（5-2）代入式（5-1），则有 AB AB A B AB AB A B S y y S x x ααsin cos +=+= ｝
（5—3） 当A 点的坐标A x 、A y 和边长AB S 及其坐标方位角AB α为已知时，就可以用上述公式计算出待定点B 的坐标。式（5—2）是计算坐标增量的基本公式，式（5—3）是计算坐标的基本公式，称为坐标正算公式。 从图5—5可以看出AB x ?是边长AB S 在x 轴上的投影长度， AB y ?是边长AB S 在 y 轴上的投影长度，边长是有向线段，是在 实地由A 量到B 得到的正值。而公式中的坐标方位角可以从0°到360°变化，根据三角函数定义，坐标方位角的正弦值和余弦值就有正负两种 情况，其正负符号取决于坐标方位角所在的象限，如图5—6所示。从式（5—2）知，由于三角函数值的正负决定了坐标增量的正负，其符号归纳成表5—3。
缓和曲线上任意点坐标计算程序
第一缓和曲线加圆曲线上任意点坐标计算程序：L1：U=U"Ｘ0"：V=V"Y0"：F"FANG"=F：E=E"LEFT-1"：LbI 0 L2：｛B｝：｛D｝：｛P｝ L3：L=AbS（B-A"ZHD"） L 4： L5：X=L-LX Y5/（40R2S2） L6：Y= LX Y3/（6RS）- LX Y7/（336RX Y3SX Y3）：G=90L2/（∏RS） L GOtO 2 7： L8：LbI 1 L9：L=L-S L10：O=90S/（∏R）+90L/（∏R） L11：M=2（Rsin（90L/∏/R）） L12：X=S-SX Y3/（40R2）+Mcos O L13：Y=S2/（6R）+MsinO:G=90S/（∏R）+180L/（∏R） L GOtO 2 14： L15：LbI 2 L16：W=tan-1（Y/X）:Q=√（X2+Y2） L 17： L E=1=>G=-G 18： L19：X[1]=U+Qcos（F+W）+Dcos（F+G+P）◢ L20：Y[1]=V+Qsin（W+F）+Dsin（F+G+P）◢ L21：GOtO 0
注、○1、XO—为起点X坐标 EXE ○2、YO—为起点Y坐标 EXE ○3、F？—方位角 EXE ○4、LEFT-1？—左偏取1右偏取0 EXE ○5、B？—所求坐标点里程（起点输0时为到起点长度）EXE ○6、ZHD？—为直缓点里程或直圆点里程（起点可以输0）EXE ○7、S？—缓和曲线长、圆曲线时输为0 EXE ○8、R？—半径EXE ○9、D？—中桩到边桩长度EXE ○10、P？—左右方向与中线切线交角、法线方向时左-90右+90 EXE ○11、上述每一步输完后必须确认、结果显示字后转到B进行循环操作。
卡西欧fx-5800P直线圆曲线坐标正反算程序
[精] fx 5800 直线圆曲线坐标正反算程序 (2010-07-01 21:50:11) 标签： 杂谈 直线计算程序 0→I:0→J:”X 0”?D:”Y0”?E:”X1”?B:”Y1”?C:Pol(B-D,C-E):J→A:If A<0: Then A+360→A:Else A→A:IfEnd:Lbl 0:?O:?S:If O≠0:Then Goto 1:IfEnd: D+Scos(A) →X:”X=”:X◢ E+Ssin(A) →Y:”Y=”:Y◢ Goto 0:Lbl 1:D+Scos(A)+Ocos(A+90) →X:”X=”:X◢ E+Ssin(A)+Osin(A+90) →Y:”Y=”:Y◢ Goto 0 输入程序时注意区别字母O与数字0 程序运行时符号说明 X0? Y0?分别输入直线起点的XY坐标值 X1? Y1?分别输入直线终点的XY坐标值 O? 输入边桩与中桩的距离(左边桩为负值,右边桩为正值),如计算中桩坐标输入0 S? 输入所求点到直线起点的距离 圆曲线计算程序 0→I:0→J: ”X0”?C:”Y0”?D:”X1”?E:”Y1”?F:?R:”L:-1 R:1”?N:”ZY”?W:Lbl 0:”LN”?T: T-W→O:Pol(E-C,F-D):I→S:J→A: If A<0:Then A+360→A:Else A→A:IfEnd: sin-1(S÷(2R)) →K:2∏RK÷180→L:180O÷(2∏R) →G:(2R)sin(G) →H:C+Hcos(A-KN+GN) →X:”X=”:X◢ D+Hsin(A-KN+GN) →Y:”Y=”:Y◢ R-0.5√(4R2-H2) →Q:”Q=”:Q◢ “S=”?V:If V=0:Then Goto 0:IfEnd:”L:-90 R:90”?U:X+Vcos(A-KN+2GN+U) →X:”X=”:X◢ Y+Vsin(A-KN+2GN+U) →Y:”Y=”:Y◢ Goto 0 输入程序时注意区别字母O与数字0 程序运行时符号说明 X0? Y0?分别输入直线起点的XY坐标值 X1? Y1?分别输入直线终点的XY坐标值 R? 输入圆曲线半径 L:-1 R:1?圆曲线向左转弯时输入-1,向右转弯时输入1 ZY? 输入起点桩号
坐标转换计算方式
72绝对坐标转换为相对坐标在直线段施工测量中，可以把绝对坐标转换为相对坐标进行放线测量，此方法比较快捷实用。 如，已知直线段线路中线A点的里程与绝对坐标X1，Y1.和其直线A点至线路前进方向的方位角a。同样已知附近的控制点Q的绝对坐标QX1，QY1.那么现在为了使用方便，要将其Q点的绝对坐标转换为相对于直线段的相对坐标，计算方法如下： 根据以上所知，根据坐标发算可以得出点A至控制点Q 的距离为L，以及点A至控制点Q方向的方位角简称R。已知线路中心线前进方向的方位角a，那么由点A至线路前进方向，和点A至控制点Q方向就形成一个夹角r，r=R-a。现在做控制点到线路中线的垂直线Y,（也就是所谓的Y坐标数据）。根据直角三角形计算方式得出Y=SIN r×L(L,是点A至点Q的距离)那么相对于线路X的坐标计算方式（X坐标表示里程）。X=COSr×L+A点里程。 即得出控制点Q相对于直线的相对坐标。 例题：例如，ZDK400至ZDK700为直线段，已知里程400的线路中心线坐标X=22580.40165 Y=27356.42893 里程700的线路中心线坐标X=22558.58105 Y=27655.63522 欲求J2点X=22562.1789 Y=27510.4874相对于400至700的相对坐标，图示如下：
解：根据已知，经过坐标反算可以求得点A至点B的坐标方位角为94 10 16 AB距离为300。 A 至D的坐标方位角为96 44 45.26 距离为155.132 那么可求得角FAD=2 34 29.26 因现已知AD=155.132 角FAD=2 24 29.26 根据三角函数可计算DF=sinfa d×AD=0.045×155.132=6.969 AF=cosfad×AD=0.999×155.132=154.975