坐标反算里程
关于坐标反算里程的超常规方法
关于坐标反算里程的超常规方法
(仅适用于隧道工程)
对于曲线隧道超开挖线放样测量,大家都觉得用坐标反算里程时计算的时间长,有的可能在三秒以内,有的可能会超过五秒钟,那让人期待的滋味的确让人难以忍受,如果不慎把坐标数据输入错误,更会让人抓狂(呵呵,因为必须得输正确坐标数据,错了要返工.),这样的无奈我也深有感触啊.有没有什么好方法让计算时间缩短呢?肯定有,下面我会详细介绍.
不管什么样式的曲线,它在某些情况下是可以看成直线的,(比如说 1.5米,2米这样的曲线线段),我说的这种方法就是把曲线当成直线来计算,这样就少去了很多计算步骤,而且能直接计算出结果.在现实施工中,我们的开挖面都是倾斜的,这个倾斜面就存在一个里程差,我们就可以把这个里程差当做直线来计算,我们只需取两个点就可以了,一个是凹进去最深的点,一个是最突出的点,我们架设好仪器后就先测算出这两点的里程,然后就把我们计算范围界定在这两个里程之间,把这两个里程点之间的当着一条直线,而且这条直线一般都很短,很少有超过2米的,然后我们在放开挖线的时候测出来的坐标都通过这条直线来反算它对应的真实里程,计算出测量点的偏移值.
接下来我们再说说精度的问题,我曾在CAD中画过这样的图,里程差取2米,半径取值为50米.隧道半幅宽度按8米计算,垂直于这条直线的距离和半径方向的距离误差也就只有0.0077米,半径50米2米长的弧线外矢距也就0.01米(现实中设计的曲线半径很少有小于50米的),完全能满足要求.且我们实际的开挖尺寸也是按直线掘进的,进尺完全大于我们所取的这两个里程差.所以这个误差完全是可以忽略不记的.
以上这种方法,我已经在实际当中用过,真的不错哟.希望同行们留意,也许能对你有所帮助,本人QQ468076885,欢迎同行灌水.。
浅谈用切线迭代法反算坐标里程
强烈 。 同时 ， 论在压 密 路基 、 不 上部 路基 、 下部 路基 底部 节
【ｌ 波 ， 英 ， 东 生 ． 一路 垂 向耦 合 系 统 的动 力 分 析 【 ． 道 ５梁 蔡 朱 车 Ｊ铁 】
２ ０ ２ （） — ． ６ 点 处 的加 速度 时 程 分析 曲线 波 动和 线 性 变化 趋 势 基 本一 学报 ， ０ ０， ２ １：５ ７１
【】 强 ． 速 铁 路 路 桥 过 渡 段 动 力 学特 性 分析 及 工 程 试 验 研 究 ６罗 高 致 ， 们 的曲线形 式具 有较 好的 吻合性。 他 【 】 都 ： 南 交 通 大 学 ，０ ３ Ｄ． 成 西 ２０ ． ３ 结语 【】 谦 ． 速 铁 路 路 基 空 闻 时 变耦 合 系统 动 力 分 析 模 璎 及 其 应 ７苏 高 通过 前 述 分析 ， 以得 到 如 下 认 识 ： 速 度 时 程 曲线 可 加 用研 究【 ． 都： 南 交通 大学 ， ０ １ Ｄ１ 成 西 ２０． 对于 越 处于路 基上 层 的土体 波动 越 大 ，波 动 幅度越 强 烈。 ［】 谦 ， 英 ． 速 铁 路 路 基 结 构 空 间 时 变 系 统 耦 合 动 力 分 析 ８苏 蔡 高 加速 度 时程 分析 曲线波 动和 线性 变化 趋势 基本 一致 ， 有 ［］ 具 Ｊ． 西南 交 通 大 学 学报 ，０ １ ３ （） ０ — １ ． ２ ０ ，６５： ９ ５ ３ ５ 较好 的 吻合性。 【】 志 红 ． 速 铁 路 轨 道 路 基 竖 向 动 力 响应 研 究 Ｉ ］ 沙 ： 南 ９聂 高 Ｄ． 长 中
【２边 学 成 ， 云 敏 于 ２ ５维 有 限 元 方 法 分析 列 车 荷 载 产 生 １】 陈 基 ．
ｇｏｕ ｉｒｔ ｎ ｎ ｕ ｅ ｙ ｈｇｈ ｐ ｅ ｒｉｓ ａ ｅ ｓ ａ ｒ ［］ 分析 【 ＿ 动 与 冲 击 ， ０ ７ ２ （） ４ ６ ｒ ｄ ｖｂａ ｉ ｓ ｉｄ ｃ ｄ ｂ ｉ －ｓ ｅ ｄ ｔａｎ ｔ Ｌ ｄ ｇ ａｄＪ ． ｏ Ｊ振 】 ２ ０ ，６３： — ７ ６
线路上点坐标反算其里程桩号及中桩距
一、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、VPol(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、WPol(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、TAbs U≤180°⇒Prog“V”△A bs U＞180°⇒U=T+360°-F：U＞G⇒Prog“W”：≠⇒Prog“X”△4、UAbsU＞180°⇒Prog“V”：≠⇒AbsU≤G⇒Prog“X”⊿AbsU＞G⇒Prog“W”⊿5、XE“LI”=H-Q（2πR÷360°×U）◢K“B”=Q（P-R）◢说明：一、本程序用于计算直线与圆曲线组合的线路上由点坐标计算相应桩号及其从中桩左、右偏移量。
由坐标反算出路线里程桩号在工程测量中的应用
由坐标反算出路线里程桩号在工程测量中的应用摘要：文章系统地介绍了由已知坐标，反向推算出该坐标位于公路路线的里程桩号及距中距离的基本概念和方法，重点介绍利用卡西欧fx-4800计算器编程，由已知坐标反向推算出公路路线中常见的直线段、圆曲线段的里程桩号方法及其在工程测量中的应用，对缩短工程测量时间，提高测量效率起到重要作用。
关键词：坐标反向推算；公路路线里程；编程；工程测量；引言在公路施工测量中，我们常用的测设方法是通过设计文件、图纸中的直线、曲线及转角表，来计算出路线的中桩、边桩坐标，再利用全站仪的放样功能，将中桩及边桩坐标在实地位置测定出来。
那么我们能不能由某一已知坐标，通过卡西欧fx-4800计算器编程计算，就能知道该坐标的准确的路线里程桩号以及其距中线的距离，甚至其设计高程，这对我们提高施工测量的效率有很大的帮助。
本文将阐述由已知坐标反向推算出公路路线中较为常见的直线段、圆曲线段里程桩号的基本概念和方法，以及介绍利用卡西欧fx-4800计算器编程，快速反算出已知坐标的里程桩号及距中距离。
1 直线段坐标反向推算路线里程桩号如图1，设直线段路线起始点O的里程桩号为K1，中桩平面坐标为（X0，Y0），路线方位角为α，P点平面坐标为（X，Y），求出P点位于该直线段的里程桩号，以及距中距离，计算步骤为：①计算直线起点O至P点的直线距离，计算公式为:②反算直线起点O至P点的坐标方位角，计算公式为:③计算直线0P与路线的夹角，计算公式为:④计算P点垂直投影于路线的垂直距离，即P点到路线中心的距离，计算公式为⑤计算P点的里程桩号，计算公式为：图1从图1中，可以看出P点的里程桩号由P点垂直投影于路线的P1点与直线起始点O之间的直线距离，与O点里程K0+000相加后，即为P点的里程桩号，L=Dsinβ即为P点位于路线右侧的距中距离，数值为正值；如果P点位于路线左侧，L=Dsinβ数值为负值。
例1：一直线段，起始桩号为K176+600，计算方位角α=299°06′58″，起点的平面坐标为X=94342.979m，Y=10235.344m，P点平面坐标为X=94387.488m，Y=10145.669m，试计算P点位于该路线的里程桩号及距中距离。
坐标正反算程序(4850)
坐标正反算程序（4850）正算主程序：ZS1.Prog “SJ”2.Lbl 0: {ZS}3.Z<A=>Rec(sqrt(S2+(Z-Q)2),F+180-tan-1(S÷(Q-Z))):X=N+I:Pause 0:Y=E+J◢4.Goto 0◣第一直线5.Z＜B=>L=Z-A:W=F+90GL2÷πRU+90:H=U:O=Z[5]:P=Z[6]:T=F+90G:Goto 1◣第一回旋线6.Z<C=>L=Z-B7.W=F+G((90U+180L)÷πR+270)8.Rec(R-GS,W):X=Z[9]+I: Pause 0:Y=Z[10]+J◢9.Goto 0◣圆曲线10.Z<D=>L=D-Z:W=F+G(K-90L2÷πR V)+90:H=V:O=Z[7]:P=Z[8]:T=F+G(K+90):Goto 1◣第二回旋线11.Lbl 1:X=L-L5÷40R2V2+L9÷3456R4V412.Y=L3÷6RV-L7÷336R3V3+L11÷42240R5V513.Rec(sqrt(X2+Y2),T):X=O+I:Y=P+J14.Rec(S,W): X=X+I: Pause 0:Y=Y+J◢15.Goto 0◣16.Z>D=>Rec(sqrt(S2+(Z-D+Z[4])2),F+G(K+tan-1(S÷(Z-D+Z[4])))):X=N+I: Pause 0:Y=E+I◢17.Goto 0◣第二直线反算主程序：FS1.Prog"SJ"2.H=90U÷πR第一回旋线所对圆心角β013.T=90V÷πR第二回旋线所对圆心角β024.Rec((Z[3] +Z[4]÷cosK)/tanK,F+90G):Z[11]=Z[5]+I:Z[12]=Z[6]+J ZH、HZ点垂线交点坐标5.Z[13]=F+90G+180 ZH点与ZH、HZ点垂线交点连线方位角6.Z[14]=Z[13]+GH HY点与圆心连线方位角7.Z[15]=Z[13]+G(K-H) YH点与圆心连线方位角8.Z[16]=Z[15]+GT HZ点与ZH、HZ点垂线交点连线方位角9.Lbl 0:{XY}10.Pol(X-Z[11],Y-Z[12]):J<0=>J=J+360◣11.GJ<GZ[13]=> Pol(N-X,E-Y):Rec(I,F-J): “Z”:Z=Q-I:Pause 0:“S”:J◢12.Goto 0◣第一直线13.GJ>GZ[16]=> Pol(X-N,Y-E):Rec(I,J-F-GK): “Z”:Z=I-Z[4]+D:Pause 0:“S”:J◢14.Goto 0◣第二直线15.GJ>GZ[13]=>Pol(X-Z[9],Y-Z[10]):J<0=>J=J+360◣16.GJ<GZ[14]=>P=Z[13]:H=U:M=A:T=1:Z[17]=Z[5]:Z[18]=Z[6]:Z[19]=F:Goto 1◣第一缓和曲线17.GJ<GZ[15]=> “Z”:Z=B+RG(J-Z[14]):Pause 0:“S”:S=G(R-I)◢18.Goto 0◣圆曲线19.GJ>GZ[15]=>P=Z[16]:H=V:M=D:T=-1:Z[17]=Z[7]:Z[18]=Z[8]:Z[19]=F+GK+180:Goto 1◣第二缓和曲线20.Lbl 1:Pol(X-Z[11],Y-Z[12]):J<0=>J=J+360:L= Abs(J-P) πR÷180 “L=H×Abs(J-P)÷2÷(90H÷πR)”21.Lbl 2:O=L-L5÷40R2H2+L9÷3456R4H4-L13÷599040R6H622.P=L3÷6RH-L7÷336R3H3+L11÷42240R5H5-L15÷9676800R7H723.Rec(sqrt(O2+P2),Z[19]+TGtan-1(P÷O)):O=Z[17]+I:P=Z[18]+J24.Pol(X-O,Y-P):Rec(I,J-(Z[19]+90TGL2÷πRH)):AbsI>0.001=>L=L+I:Goto 2:≠>“Z”:Z=M+TL: Pause 0: “S”:TJ◢25.Goto 0◣子程序(曲线要素数据)：SJ1.Defm 19:M:M=1=>F=*:K“ZJ”=*:G=*:U(Ls1)=*:V(Ls2)=*:R=*:Q=*:N“XJ”=*:E“YJ”=*:Goto 1◣2.M=2=>F=*:K“ZJ”=*:G=*:U(Ls1)=*:V(Ls2)=*:R=*:Q=*:N“XJ”=*:E“YJ”=*:Goto 1◣3.………………4.Lbl 1:Z[1]=U2÷24R-U4÷2688R3+U6÷506880R5Z[1]第一回旋线内移值P15.Z[2]=V2÷24R-V4÷2688R3+V6÷506880R5Z[2]第二回旋线内移值P26.M=(R+Z[2]-(R+Z[1])cosK)÷sinK7.P=(R+Z[1]-(R+Z[2])cosK)÷sinK8.Z[3]= M+U÷2-U3÷240R2Z[3]第一回旋线切线长T19.Z[4]= P+V÷2-V3÷240R2Z[4]第二回旋线切线长T210.L=RKπ÷180+(U+V)÷2曲线长(Ls1+圆+Ls2)11.A=Q-Z[3]:B=A+U:D=A+L:C=D-V ZH,HY,HZ,YH里程桩号12.I=0:J=0:Rec(Z[3],F+180):Z[5]=N+I:Z[6]=E+J ZH点坐标增量计算，ZH点坐标13.Rec(Z[4],F+GK):Z[7]=N+I:Z[8]=E+J HZ点坐标增量计算，HZ点坐标14.Rec(sqrt(M2+(R+Z[1])2),F+G(K+tan-1((R+Z[2])÷P))):Z[9]=N+I:Z[10]=E+J 圆心坐标注：F－起始边方位角K－转角G－线路左转为－，右转为+U－第一回旋线长V－第二回旋线长R－半径Q－交点桩号N,E－交点坐标Z－待求点桩号S－偏距(左偏为-,右偏为+)sqrt为根号。
线路测量坐标正、反算计算原理及卡西欧5800程序说明
线路测量坐标正、反算计算原理及卡西欧fx-5800P程序说明一、计算原理在各测量书中对于坐标正算的相关计算式均有说明，故在此不做详解，仅对正算过程中需要用到的原理及公式做一汇总。
对于坐标反算，虽然都采用无限趋近原理进行计算，但计算方式各有差别，本文仅对其中一种自认为相对简单易懂并便于操作的原理进行详解。
1.1 坐标转换[1]如图1，设X P、Y P为P点在国家控制网坐标系中的坐标；x P、y P为P点在工程独立控制网坐标系中的坐标。
X O、Y O为工程独立坐标系原点o在国家坐标系中的坐标，Δα为两坐标系纵坐标轴的交角。
如果一条边在国家坐标系中的坐标方位角为A，而在工程独立坐标系中的坐标方位角为α，则：∆α=A−α(1-1)当由工程独立坐标系坐标换算至国家坐标系坐标时，换算公式为：X=x cos∆α−y sin∆α+X O(1-2)Y=x sin∆α+y cos∆α+Y O}当由国家坐标系坐标换算至工程独立坐标系坐标时，也可使用式(1-2)，此时应将X、Y与x、y互换，且∆α=α−A。
1.2 坐标方位角关系计算1.2.1 正、反坐标方位角[2]一条直线的坐标方位角与直线的前进方向有关，沿直线前进方向的坐标方位角称为正坐标方位角，与其相反方向的坐标方位角称为反坐标方位角。
如图2，由于轴子午线之间是互相平行的，因此同一直线的正、反坐标方位角相差180°，即：α正=α反±180°(2-1)当α反<180°时，取“+”号；当α反>180°时，取“-”号。
1.2.2 坐标方位角的推算[3] 1.2.2.1 转折角为右角如图3(a)，α12为已知边坐标方位角，α23为推算边的坐标方位角，β右为该两边所夹的右角，则：α23=α12±180°−β右=α21−β右 (2-2)1.2.2.2 转折角为左角如图3(b)，α12为已知边坐标方位角，α23为推算边的坐标方位角，β左为该两边所夹的右角，则：α23=α12+β左±180°=α21+β左(2-3) 无论用右角还是左角推算，如遇出现负数的情形，应加上360°。
坐标反算里程偏距公式
坐标反算里程偏距公式引言在土木工程、道路建设以及测绘等领域中，经常需要进行坐标反算和测量。
坐标反算是指根据已知的坐标点，通过一定的计算方法来确定点之间的距离和方位角。
其中，里程偏距公式是一种常用的计算方法，用于确定某一点在车辆行驶里程上的位置与其与公路中心线垂直方向的距离。
坐标反算和里程偏距公式坐标反算是通过已知的坐标点和其对应的里程数，来确定另一点在公路行驶里程上的位置以及其与公路中心线的垂直距离。
里程偏距公式是用于计算该垂直距离的数学模型。
坐标反算在坐标反算中，我们需要知道基准点的坐标和里程数，以及目标点的坐标。
假设我们有两个已知的坐标点A和B，对应的里程数分别为a和b。
我们希望确定目标点C的里程数和垂直距离。
首先，我们可以通过坐标点A和B的坐标差来计算它们之间的水平距离d：d = sqrt((xb - xa)^2 + (yb - ya)^2)接下来，我们可以使用三角形相似性原理来计算目标点C的里程数x和垂直距离y。
假设AC和BC分别为三角形ABC的两条边，根据三角形相似性原理，我们可以得到以下关系：x / (b - a) = y / d通过变换，我们可以得到目标点C的里程数公式：x = a + (b - a) * (xb - xa) / d里程偏距公式已知目标点C的里程数x，我们可以通过里程偏距公式来计算其与公路中心线的垂直距离y。
里程偏距公式通常是一个多项式函数，其形式可以根据实际情况进行调整。
一个常用的里程偏距公式是二次函数形式：y = A * x^2 + B * x + C其中，A、B和C是公式的系数，可以根据实际测量数据进行拟合。
根据目标点C的里程数x，我们可以计算出其对应的垂直距离y。
使用示例现在我们来看一个使用坐标反算和里程偏距公式的实际示例。
假设我们有两个已知坐标点A和B的信息如下： - 坐标点A：(3, 5)，里程数：50km - 坐标点B：(6, 8)，里程数：100km我们希望确定目标点C的里程数和垂直距离。
坐标反算公式、示意图及相应5800程序案例
反算原理在图中，A点为已知坐标而待求对应中桩桩号及边距的点。
B点为假定的A点对应中桩桩号点。
显然，B点并不对应于A点。
做出B点的切线，过A点做辅助线垂直于B点的切线，相交于C点。
设：B点的切线方位角为α，B点的桩号为K B,B点的坐标分别为X B、Y B，A点的桩号为K A,A点的坐标分别为X A、Y A，“B-A”的方位角为β，“B-A”的距离为N，“B-C”的距离为L，“C-A”的距离为Z。
根据前面的坐标正算的公式可以得到α，X B、Y B值。
根据计算器内置的Pol（X A-X B，Y A-Y B）公式（直角坐标转换为极坐标）能得到β，N值。
（1）（2）当L=0时，B点是对应于A点的，K B=K A，Z即为A点的距中桩的距离。
当L≠0时，则采用K B=K B+L,对B点进行新的假定，进而再次对L进行解算，直至L=0，或则L值在容许误差范围之内。
坐标正反算通用程序(终极篇)1. 坐标正算主程序(命名为KP-XY)第1行：Lbl 0“DK=”?K:“PJ=”? P:“α=”?W第2行：Prog “GP”(GP为线路名称)第3行：“X=”:X+Pcos(F+W)→X◢第4行：“Y=”:Y+Psin(F+W)→Y◢第5行：“FWJ=”:F►DMS◢第6行： Goto 0K——计算点的里程PJ——计算点到中桩的距离（左负右正）α——取线路前进方向右夹角为正X、Y--为计算坐标、F为方位角GP--为线路平曲线名称2. 坐标反算桩号和偏距主程序(命名为XY-KP)第1行：“DK=”?K:“X:=”? U:“Y:=”?V第2行：Lbl 0Prog “GP”第3行：Pol(U-X,V-Y):Icos(F-J)→Z[1]:K+Z[1]→K 第4行：Abs(Z[1])>0.0001=>Goto 0第5行：“DK=”:K◢第6行：“PJ=”:Isin(J-F)→P◢X:= —取样点的X坐标Y:= —取样点的Y坐标DK= 输入时为计算起始点(在线路内即可)，输出时为反算点的桩号P—偏距（左负右正）3. 计算坐标子程序(命名为XYF)为了简洁，本程序由数据库直接调用，上述中的正反算主程序不直接调用此程序第1行：K-G→S：(D-C)÷L→I第2行：X+∫(cos(F+X(2C+XI)×90÷π),0,S)→X第3行：Y+∫(sin(F+X(2C+XI)×90÷π),0,S)→Y第4行：F+S(2C+SI)×90÷π→F第5行：F<0=>F+360→F: F>360=>F-360→F4. 数据库（命名为A）第1行：0→G:3925525.975→X:502796.176→Y:3°47′30.2″→F:1/996→C:1/180→K≤G+L =>GoTo1（第一缓和曲线）第2行：245.607→G: 3925566.736→X:502801.006→Y:11°30′39.4″→F: 1/180→C:1/180→D:91.027→L:K≤G+L =>Goto1（圆曲线）第3行：318.522→G: 3925647.682→X:502804.484→X:40°29′08.7″→F: 1/180→C: 1/295.54→D: 70.37→L:K≤G+L =>Goto1（第二缓和曲线）第4行：Lbl 1:Prog “XYF”G——曲线段起点的里程X——曲线段起点的x坐标Y——曲线段起点的y坐标F——曲线段起点的坐标方位角C——曲线段起点的曲率(半径倒数，直线为0，左负右正)D——曲线段终点的曲率(半径倒数，直线为0，左负右正)L——曲线段长度（尽量使用长度，为计算断链方便）说明：（1）适用于任意线形：直线（0→C、0→D）、圆曲线（圆半径倒数→C、圆半径倒数→D）、缓和曲线（0或圆半径倒数→C、圆半径倒数或0→D）、卵形曲线（接起点圆的半径倒数→C、接终点圆的半径倒数→D），曲线左转多加一负号。
5800计算器全线坐标计算放样正反算程序
5800计算器全线坐标计算放样正反算程序（定稿）彭赐明主程序坐标正反算程序名称：ZBZFSLB1 A↙Cls : Fix 4 : 30→Dim Z :"0=ZS,1=FS"?Z ↙If Z=0:Then Goto B:IfEnd↙(Z=0进入里程点坐标正算）If Z=1:Then Prog"FSLCBZ ": Goto A:IfEnd↙(Z=1进入反算里程边桩）Lb1 B ↙“K=”?K ：（计算里程）Prog"SJK1 ":Prog”ZBFY” :Goto A↙子程序数据库 SJK1IF K<本段曲线终点里程 AND K≥上段曲线终点里程：THEN 本段终点里程→Z[1] : 上段曲线终点里程→Z[2] ：1→O（注：左偏曲线输入-1→O,右偏曲线输入1→O）: 半径→R : 曲线偏角→A:第一缓和曲线→Z[6] : 第二缓和曲线→Z[7] : 交点X→B :交点Y→C : 小里程向交点方位角→E : 交点向大里程方位角→F : Prog”XLZBJSCX”:Return: IFEND↙IF…………Prog”XLZBJSCX”:Return:IFEND（曲线段分段输入）↙补充直线段输入如下(单独直线段）IF K<本段直线终点里程AND K≥本段直线起点里程:THEN 本段直线终点里程→Z[3]:终点X→Z[16]:终点Y→Z[17]:方位角→E: Z[16]+ (K- Z[3])*COS E→Z[18]:Z[17]+ (K- Z[3])*SIN E→Z[19] : Return:IFEND子程序坐标正算，名称：XLZBJSCXLB1 2 ↙（曲线要素计算）Z[6]/2- Z[6]^3/(240*R^2)+ Z[6]^5/(34560*R^4) →Z[8] ↙（M1）Z[7]/2- Z[7]^3/(240*R^2)+ Z[7]^5/(34560*R^4) →Z[9] ↙（M2）Z[6]^2/(24*R)- Z[6]^4/(2688*R^3) →Z[10] ↙（P1）Z[7]^2/(24*R)- Z[7]^4/(2688*R^3) →Z[11] ↙（P2）π*A*R/180+0.5*( Z[6]+ Z[7])→Z[25] ↙（曲线总长）90* Z[6]/(R*π) →Z[14] ↙（第一缓和曲线总偏角）90* Z[7]/(R*π) →Z[15] ↙（第二缓和曲线总偏角,可以省略）Z[8]＋（R+Z[10])TAN(A/2)-(Z[10]-Z[11] )/SIN A→Z[12]↙(切线T1)Z[9]＋（R+Z[11])TAN(A/2)+(Z[10]-Z[11] )/SIN A→Z[13]↙(切线T2)B+ Z[12]*COS (E+180)→ Z[16] ↙（ZH点X）C+ Z[12]*SIN(E+180)→ Z[17] ↙（ZH点Y）Z[1]- Z[25]→Z[3] ↙(ZH点里程)Z[3]+ Z[6]→Z[4] ↙(HY点里程)Z[1]- Z[7]→Z[5] ↙(YH点里程)GOTO 3 ↙LB1 3 ↙(判断里程点与曲线关系)IF K≤Z[3] AND K> Z[2] : THEN GOTO 4 : IFEND ↙IF K≤Z[4] AND K> Z[3] : THEN GOTO 5 : IFEND ↙IF K≤Z[5] AND K> Z[4] : THEN GOTO 6 : IFEND ↙IF K≤Z[1] AND K> Z[5] : THEN GOTO 7 : IFEND ↙LB1 4 ↙（里程小于直缓点直线独立坐标）K- Z[3] →Z[23] : 0→Z[24] : E→T : GOTO 8↙LB1 5 ↙（第一缓和曲线独立坐标）K- Z[3] →H ↙H-H^5/(40*R^2* Z[6]^2)+H^9/(3456*R^4* Z[6]^4) →Z[23] ↙H^3/(6*R* Z[6])-H^7/(336*R^3* Z[6]^3) →Z[24] ↙90*H^2/( R*π* Z[6]) →T ↙IF O>0 ：THEN T +E→T : ELSE E-T →T : T<0=>360+T→T : IFEND ↙GOTO 8 ↙LB1 6 ↙（圆曲线独立坐标）K- Z[4] →H ↙H*180/( R*π)+ Z[14]→T ↙R*SIN T+ Z[8]→Z[23] ↙R*(1-COS T)+ Z[10]→Z[24] ↙IF O>0 ：THEN T +E→T : ELSE E-T →T : T<0=>360+T→T : IFEND ↙GOTO 8 ↙LB1 7 ↙（第二缓和曲线独立坐标）Z[1] -K →H ↙H-H^5/(40*R^2* Z[7]^2)+H^9/(3456*R^4* Z[7]^4) →U↙H^3/(6*R* Z[7])-H^7/(336*R^3* Z[7]^3) →V ↙90*H^2/( R*π* Z[7]) →T ↙Z[13]COS A+ Z[12]-U*COS A-V*S IN A→Z[23] ↙Z[13]*SIN A-U*SIN A+V*COS A→Z[24]↙IF O>0 ：THEN F-T→T : T<0=>360+T→T : ELSE F+T →T : IFEND ↙GOTO 8 ↙LLb1 8IF O<0 : THEN - Z[24]→Z[24] : IFEND ↙Z[16]+Z[23]*COS E-Z[24]*SIN E→Z[18] ↙Z[17]+Z[23]*SIN E+Z[24]*COS E→Z[19] ↙ReTurn↙子程序反算里程边桩，名称：FSLCBZLb1 1↙"XK="?X:"YK="?Y↙（输入任意测点的XY坐标）:“K=”?K （试算里程，任意输入标段内里程点即可，也可以输入大致的估算里程加快速度）↙Lb1 2↙Prog"SJK1 ":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos W-(X-Z[18])*Sin W)→S↙If S<0.0001:Then Goto 4:Else Goto 3:Ifend↙Lb1 3↙K+S→K: Prog"S JK1":T-90 →W:W<0=>360+W→W:Abs((Y-Z[19])*Cos W-(X-Z[18])*Sin W)→Q↙If Q<0.0001 :Then Q→S: Goto 4:Else IF Q<S: THEN K+Q→K: Goto 2 :Else IF Q>S :THEN K-Q→K:G oto 2:Ifend:Ifend:Ifend↙Lb1 4↙Pol (X-Z[18],Y-Z[19]:"DP(-Z+Y)=":I◢(偏距）（由于该程序不能准确判断边桩左右方向，暂作修改，取消左右边判定）"K=":K+S→K◢（里程）Return↙子程序坐标放样：ZBFYLB1 0 ↙“XHS="?G(后视点X):"YHS="?L(后视点Y):"XZJ="?M(置镜点X):"YZJ="?N(置镜点Y):Pol(G-M,L-N):"DH=":I(后视距)◢J<0=>J+360→J:"FH=":J→DMS◢(后视方位角) “QXJ=” :T◢（计算里程点切线方位角，可以不显示）“XI=” : Z[18] ◢（中线X）“YI=” : Z[19] ◢（中线Y）Pol(Z[18]-M,Z[19]-N):"DI=":I（中桩放样距）◢J<0=>J+360→J:"FI=":J→DMS◢（中桩放样方位角）“PJ=”?P◢(输入边桩与线路夹角，左-右+)“PD=”?D◢（输入边桩距）Z[18]+D*COS(T+P) →Z[20] ↙Z[19]+D*SI N(T+P) →Z[21] ↙“XP=”: Z[20] ◢（边桩X）“YP=”: Z[21] ◢（边桩Y）Pol(Z[20]-M,Z[21]-N):"DP=":I◢（边桩放样距）J<0=>J+360→J:"FP=":J→DMS◢（边桩放样方位角）Return↙。
铁路线路任意里程坐标正反算程序使用说明
铁路线路任意里程坐标正反算程序(有需要程序的可联系陈工，QQ：285242895)1、程序开发背景在铁路线路测量中,在曲线要素已定的情况下,已知某点的里程及距中线的距离,计算该点的坐标,我们称之为线路坐标正算。
相反地,已知某点的坐标,确定该点在已定线路中的里程及距中线距离的过程,我们称之为线路坐标反算。
对于一条完整的曲线，它包括直线、第一缓和曲线、圆直线及第二缓和曲线。
而一条完整的铁路线路，通常都包含不止一条曲线，如果我们根据铁路线路多个曲线的曲线要素，构建一个线路模型，然后给出任意里程点，自动计算出对应的线路坐标，也可以给出任意坐标，计算出对应的线路里程和偏距，这将在测量和放样工作有着较为实际的应用。
比如用于逐桩坐标计算、隧道开挖及土石方开挖、线路征地界坐标计算、线路测量中线质量的检查、地质钻孔位置、桥梁桩基坐标计算等方面。
2、程序界面3、程序功能1、可以根据点的里程及距中线的距离，计算出该点的坐标，显示数据文件导入结果及计算结果，最后以csv格式文件保存计算的里程数据成果及曲线要素。
2、可以根据任意点的坐标，计算出点在已定线路中的里程及距中线距离，同时显示数据计算结果，最后以.zb格式文件保存计算的坐标成果。
4、程序特色3.1 本程序采用易于交互操作的对话框模板和MSFlexGrid控件，在MFC开发环境下利用VC语言进行编写，整个程序的计算过程及结果均可在图表中直接呈现，便于数据的检查，整个程序的界面简洁直观，功能清晰、易学易用。
3.2 结合铁路测量的实际情况，在导入曲线要素时，不需要输入曲线的五大桩要素以及曲线偏向，只需要曲线数据文件中包含曲线半径、缓长及曲线两侧各两个直线点坐标，就可以计算出其他曲线要素，进而构建完整的线路模型。
3.3 在线路正算时，里程数据既可以从文件中导入，也可以在程序界面上获取。
当采用从文件导入时，里程数据可以是乱序排列的。
当从界面获取时，程序可以自动计算出连续里程数据，3.4 在线路反算中，当我们给定任意点的坐标时，程序不仅可以计算出对应线路中的里程、距离及垂点坐标，还可以计算出此点是否在线路对应范围内以及位于曲线上的具体位置。
9860线元法隧道测量(坐标正反计算、高程)计算程序
CASIO 9860G SD线元法隧道三维（坐标正反计算、高程）计算程序1. A（此为主程序）Lbl 0：〝1.LC=>XY〝：〝2.XY=>LC〝：〝3.ZHZL=>GC〝：〝PB=>V=1,2,3〞？→V：If V=1：Then GOTO 1 ：IfEnd ：If V=2：Then GOTO 2 ：IfEnd ：If V=3：Then GOTO 3 ：Else GOTO 0 ：IfEnd：Lbl 3：〝ZH=H〝？→H ：〝SDZF=Z〝？→Z：Prog〝ZGCZCX〝：GOTO 0：Lbl 1 ：〝ZH=L〝？→L：If L>173000 And L<174661.96：Then GOTO 4 ：Else GOTO 0 ：IfEnd ：Lbl 4：L→L：〝SDZF=Q〝？→Q：〝XLZJ,-Z+Y=Q〝：Q+0.125→Q：Prog 〝ZBQXYS〝：〝JSJD=J〝：90→J▲Prog 〝ZSZB〝：〝X=〝：X ▲ 〝Y=〝：Y▲〝FWJ=O〝：O▼DMS▲L→H：Q-0.125→Z：Prog 〝ZGCZCX〝：GOTO 0：Lbl 2：〝XO=M〝?→M：〝YO=R〝?→R：173300→L：If M>3845505.273 And M<3846506.099And R>499371.832 And R<500352.224 ：Then GOTO 5：Else GOTO 2：IfEnd ：Lbl 5：0→Q：0→J：Prog 〝ZBFS〝：〝LC=L〝：L ▲〝JL=Q〝：Q▲ 〝SDZJ，-Z+Y=Q〝：Q-0.125→Q ▲L→H：Q→Z：Prog 〝ZGCZCX〝：GOTO 02.正算坐标ZBZS( L-S ) / 4→H：90/π→F：HHF（1/T-1/I）/（K-S）→U：2HF/ I→D：C+4D+16 U→Ｏ：O+J→P ：C+ D+ U→E：C+2D+4U→W：C+3D+9U→G：A+AbsH/3*(cosC+4(cosG+cosE)+2cosW+cosO)+Qcos P→X ：B+AbsH/3*(sinC+4(sinG+sinE)+2sinW+sinO)+Qsin P→Y3.反算坐标：ZBFSLbl 0：Prog 〝ZBQXYS〝：Prog 〝ZBZS〝：O-90→Z：(R-Y)cosZ-(M-X) sinZ→P ：L+P→L：If Abs P≥0.001：Then GOTO 0 ：Else GOTO 1：IfEnd ：Lbl 1： (R-Y)cosO-(M-X) sinO→Q4. 曲线元要素数据库：ZBQXYSIf L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：If L≥S And L＜K ：Th en **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：……………………………If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ┘(注：如有多个曲线元要素继续添加入数据库ZBQXYS中)5 高程计算主程序ZGCZCX (后有修改说明)Lbi1 ：〝SCGC=X〝？→X：〝R=M〝：5.98→M：〝CS=N〝：1.603→N：Prog＂GC SJK＂：C-D→F：Abs(RF÷2) →T：R AbsF÷F→R：If H≤B-T ：Then 0→K：GOTO 2：IfEnd ：If H>B-T And H<B+T ：Then H-B+T→K ：GOTO 2：IfEnd ：If H≥B+T ：Then 0→K：D→C：GOTO 2：IfEnd ：Lbi 2 ：〝XLZG=G〝：A-（B-H）C-K^2÷2R→G▲If Z≥0：Then If X>100 Then 〝YO1XGC=J〝：G+N→J▲〝YGCFSKD=J〝：Abs√(M^2-(X-（G+N）)^2) →J▲〝YKDCQ,+C,-Q=J〝：J- Abs (Z+0.000) →J▲〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J▲〝YGCCQ,C+,Q-=J〝：X-J→J▲Else〝YO1XGC=J〝：G+N→J▲〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J▲ IfEnd ：Else 〝ZO1XGC=J〝：G+N→J▲If X>100 ：Then〝ZGCFSKD=J〝：Abs√(M^2-(X-（G+N）)^2) →J▲〝ZKDCQ,+C,-Q=J〝：J- Abs (Z+0.000) →J▲〝ZSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J▲〝ZGCCQ,C+,Q-=J〝：X-J→J▲Else〝ZSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J ▲IfEnd ：IfEnd6高程计算主程序子程序：GCSJKIf H>起点桩号 And H≤第一个竖曲线起点桩号Then 第一竖曲线交点高程→A：第一竖曲线交点桩号→B：第一竖曲线前坡度→C：第一竖曲线后坡度→D：第一竖曲线半径→R：IfEnd：If H>第一竖曲线止点桩号 And H≤第二竖曲线起点桩号Then 第一竖曲线交点高程→A：第一竖曲线交点桩号→B：第一竖曲线前坡度→C：第一竖曲线后坡度→D：第一竖曲线半径→R：IfEnd：If H>第一竖曲线止点桩号 And H≤第二竖曲线起点桩号Then 第一竖曲线交点高程→A：第一竖曲线交点桩号→B：第一竖曲线前坡度→C：第一竖曲线后坡度→D：第一竖曲线半径→R：IfEnd………………………(继续添加要素)说明：第一部分坐标部分（1、2、3、4）V=1进入坐标正算 V=2进入坐标反算V=3进入单独的高程计算当V不等于1、2、3时，则返回程序，要求再次输入V值。
坐标正反算
坐标正反算一、坐标正算与坐标反算1、坐标正算已知点的坐标、边的方位角、两点间的水平距离，计算待定点的坐标，称为坐标正算。
如图6-6所示，点的坐标可由下式计算：式中、为两导线点坐标之差，称为坐标增量，即：=1000、=1000、方位角【例题6-1】已知点A坐标，=35°17＇36.5\，2、坐标反算已知两点的坐标，计算两点水平距离=200.416，计算点的坐标？35o17＇36.5\35o17＇36.5\两点的水平距离与坐标方位角，称为坐标反算。
如图6-6可知，由下式计算水平距离与坐标方位角。
（6-3）（6-4）式中反正切函数的值域是-90°～+90°，而坐标方位角为0°～360°，因此坐标方位角的值，可根据正切角值换算为坐标方位角。
、的正负号所在象限，将反【例题6-2】=3712227.860、、水平距离=3712232.528、=523620.436、=523611.598，计算坐标方位角计算坐标方位角。
=62°09＇29.4\29.4\注意：一直线有两个方向，存在两个方位角，式中：、的计算是过A点坐标纵轴至直线的坐标方位角，若所求坐标方位角为，则应是A点坐标减点坐标。
坐标正算与反算，可以利用普通科学电子计算器的极坐标和直角坐标相互转换功能计算，普通科学电子计算器的类型比较多，操作方法不相同，下面介绍一种方法。
【例题6-3】坐标反算，已知=1771.03、=2365.16、=1181.77、、水平距离。
=1719.24，试计算坐标方位角键入1771.03-2365.16按等号键［=］等于纵坐标增量，按储存键［］，键入1719.24-1181.77按等号键［=］等于横坐标增量，按［］键输入，按［］显示横坐标增量，按［］键输入，按第二功能键［2ndF］，再按［］键，屏显为距离，再按［］键，屏显为方位角。
【例题6-4】坐标正算，已知坐标方位角=200.40，试计算纵坐标增量横坐标增量=294°42＇51\，。
纬地道路坐标反算里程偏距
纬地道路坐标反算里程偏距以纬地道路坐标反算里程偏距为标题，我们来探讨一下这个话题。
在许多导航系统和地图应用中，我们经常会看到里程偏距这个概念。
那么，什么是纬地道路坐标反算里程偏距呢？为什么要进行反算呢？让我们一起来了解一下。
纬地道路坐标反算里程偏距是一种根据地理坐标反算相应里程的方法。
在现代导航系统中，我们可以通过GPS定位来获取当前位置的纬度和经度信息。
而在地图上，道路通常是以经纬度的形式表示的。
因此，如果我们知道了起点和终点的经纬度，就可以通过纬地道路坐标反算里程偏距的方法来计算出两点之间的里程。
在进行纬地道路坐标反算里程偏距之前，我们需要先了解一下地球的形状。
地球并不是一个完全的球体，而是一个稍微扁平的椭球体。
这就意味着，地球上不同纬度的距离并不相同。
因此，在进行纬地道路坐标反算里程偏距时，我们需要考虑这个因素。
我们需要将起点和终点的经纬度转换为弧度表示。
然后，我们可以利用球面三角学中的一些公式来计算出两点之间的弧长。
最后，我们将弧长转换为里程，就可以得到纬地道路坐标反算里程偏距了。
纬地道路坐标反算里程偏距的计算方法相对复杂，需要使用一些数学公式和算法。
在这里，我们就不再深入讨论了。
如果你对此感兴趣，可以搜索相关的资料进行学习。
那么，纬地道路坐标反算里程偏距有什么应用呢？首先，它可以用于导航系统中的路径规划。
通过纬地道路坐标反算里程偏距，我们可以根据起点和终点的经纬度信息，计算出最短路径的里程。
这对于司机来说非常实用，可以帮助他们选择最优的行驶路线。
纬地道路坐标反算里程偏距还可以用于交通流量分析。
通过分析不同道路之间的里程偏距，我们可以了解到不同路段的交通状况。
比如，如果两个道路的里程偏距较小，说明这两个道路之间的交通流量较大，可能存在拥堵情况。
这对于交通管理部门来说非常重要，可以帮助他们进行交通拥堵预测和调度。
除了导航系统和交通流量分析，纬地道路坐标反算里程偏距还可以应用于地理信息系统（GIS）和地图制作等领域。
卡西欧5800公路坐标正反算程序
卡西欧5800公路坐标正反算程序(总19页)-CAL-FENGHAI.-(YICAI)-Company One1-CAL-本页仅作为文档封面，使用请直接删除目录一、坐标正算基本公式 (02)二、坐标反算原理 (04)三、高程数据库录入变换 (05)四、计算器程序 (07)01、ZBZS(坐标正算) (07)02、ZBFS（坐标反算） (08)03、GCJF(高程积分) (09)04、PJFY（坡脚放样） (10)05、JFCX（积分程序） (11)06、ZBFY（坐标放样） (11)07、DT（递推） (12)08、HP（横坡） (13)09、LK（路宽） (14)10、SJK1（平面数据库） (14)11、SJK2（纵面数据库） (14)12、SJK3（左路宽度数据库） (15)13、SJK4（右路宽度数据库） (15)14、SJK5（横坡数据库） (16)15、SJK6（下边坡数据库） (16)16、SJK7（左上边坡数据库） (17)17、SJK8（右上边坡数据库） (18)五、后记 (19)CASIO 5800计算器公路工程测量程序一、正算所涉及的计算公式 X R i d XαβBd Y d l d βI图表 1在图1中，A 点为回旋曲线起点，B 点为回旋曲线止点，I 点为所求坐标点。
设：A 点的X 坐标为X A ，Y 坐标为Y A ，A 点的切线方位角为α，A 点的曲率为ρA ,A 点的里程为L A ,B 点的曲率为ρB ,B 点的里程为L B ,I 点的曲率为ρI ,I 点的里程为L I 。
I 点的切线角为β。
由于回旋线上各点曲率半径R i 和该点至曲线起点的距离L 成反比。
故此任意点的曲率为;CL R i i ==1ρ (c 为常数) （1） 由式(1)可知，回旋曲线任意点的曲率按线性变化，由此回旋曲线上里程为L i 点的曲率为;AB A i A B A i L L L L --⨯-+=)(ρρρρ （2） 当曲线右偏时ρB 、ρA 取正值，反之取负值。
fx-5800p全线坐标正反算带高程计算程序(线元法
曲线任意里程中边桩坐标正反算(CASIO fx-5800P计算器)程序一、程序功能及原理1.功能说明：本程序由一个主程序(TYQXJS)和五个子程——正算子程序(SUB-ZS)、反算子程序(SUB-FS)等构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。
本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。
本程序在CASIO fx-5800P计算器运行。
2．计算原理：利用Gauss-Legendre 5点通用公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb<0=>该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。
二、源程序1.主程序(TYQXJS)Deg:fix 3119→DimZ“INPUT(0) Or DATA(Else)”?ILbl 0:“1.SZ=>XY,2.XY=>SZ，3.TF=>CK,4.SD=>FY，5.TW=>FY”?NIf N=1 Or N=5:Then Goto 1Else If N=2 Or N=3 Or N=4:Then Goto 2Else Goto 3IfEnd:IfEndLbl 1:“K(m)=”?SIf S<0:Then Goto 0:IfEnd“JL(m)=”?ZIf Z≠0:Then “ANGLE→R(Deg)=”?M:IfEndIf I=0:Then Prog “DAT1”:Else Prog “DAT2”:IfEndS-O→W:If W<0:Then Goto 0:Else If W>H:Then Goto 0:IfEnd:IfEndProg “SUB-ZS”:Prog “SUB-GC”If Z<0:Then“XL(m)=”:X◢“YL(m)=”:Y◢If N=5:Then Prog “SUB-TW”:IfEndElse If Z>0:Then “XR(m)=”:X◢“YR(m)=”:Y◢If N=5:Then Prog “SUB-TW”:IfEndElse “X(m)=”:X◢“Y(m)=”:Y◢“Hs(m)=”:L◢“FWJ=”: F►DMS◢IfEnd:IfEndGoto 1Lbl 2:“X(m)=”?X:If X<0:Then Goto 0:IfEnd“Y(m)=”?YIf N=3 Or N=4:Then “H(m)=”?→Z[2]:IfEndX→Z[4]:Y→Z[5]:90→MIf I=0:Then Prog “DAT1”G+QEH(C+HD)→F:1→J:U→Z[8J]:V→Z[8J+1]G-M→A:F-M→B:Prog “SUB-ZX1”:Goto BG+M→A:F+M→B:Prog “SUB-ZX1”:Goto B:IfEnd0→J:Lbl A: Isz JZ[8J+3]-M→A:Z[8(J+1)+3]-M→B:Prog “SUB-ZX1”If Z[6]Z[7]<0:Then J→J:Prog“DAT1” :Goto B:IfEndZ[8J+3]+M→A:Z[8(J+1)+3]+M→B:Prog “SUB-ZX1”If Z[6]Z[7]<0:Then J→J:Prog “DAT1” :Goto B:Else Goto A:IfEndLbl B:Prog“SUB-FS”O+W→S:Prog “SUB-GC”“K(m)=”:S◢“Hs(m)=”:L◢“JL(m)=”:Z◢If N=3:Then Prog “WIDE”:IfEndIf N=4:Then Prog “SDFY”:IfEndGoto 2Lbl 3:“TYQXJS→END”2. 正算子程序(SUB-ZS)0.1184634425→A: 0.2393143352→B:0.2844444444→Z[1]:0.0469100770→K:0.2307653449→L:0.5→Z[3]U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Z[1]cos(G+QEZ[3]W(C+Z[3]WD))+Bcos(G+QE(1 -L)W(C+(1-L)WD))+Acos(G+QE(1-K)W(C+(1-K)WD)))→XV+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Z[1]sin(G+QEZ[3]W(C+Z[3]WD))+Bsin(G+QE(1-L)W(C+(1-L)WD))+Asin(G+QE(1-K)W(C+(1-K)WD)))→YG+QEW(C+WD)+M→FX+Zcos(F)→X:Y+Zsin(F)→YF-M→F:If F<0:Then F+360→F:IfEndReturn3. 反算子程序(SUB-FS)G-M→T:Abs((Y-V)cos(T)-(X-U)sin(T))→W:0→ZLbl C:Prog “SUB-ZS”T+QEW(C+WD)→L:(Z[5]-Y)cos(L)-(Z[4]-X)sin(L)→ZIf Abs（Z）<10-6:Then Goto D:Else W+Z→W：Goto C:IfEndLbl D:0→Z:Prog “SUB-ZS”:(Z[5]-Y)÷cos(F)→Z:Return4．高程计算子程序（SUB-GC）Prog “DAT3”0．5RAbs(0.01(E-D))→TIf E>D:Then 1→C:Else -1→C:IfEndIf S<A:Then -0.01D→W:Else 0.01E→W:IfEndAbs(A-S)→LIf L>T:Then B+WL→L:Else B+WL+C(T-L)2÷2÷R→L:IfEndReturn5. 垂距计算子程序(SUB-ZX1)(Z[5]-Z[8J+1])cos(A)-(Z[4]-Z[8J])sin(A)→Z[6](Z[5]-Z[8(J+1)+1])cos(B)-(Z[4]-Z[8(J+1)])sin(B)→Z[7]Return6.曲线元要素数据库：DAT1(两条线路，可扩充为多条线路)If I=0:Then Goto E: Else If I=2:Then Goto F:IfEnd:IfEndJ=1=>***→Z[8]:***→Z[9]:***→Z[10]:***→Z[11]:***→Z[12]:***→Z[13]:***→Z[14]:*** →Z[15]：IfEndJ=2=>***→Z[16]:***→Z[17]:***→Z[18]:***→Z[19]:***→Z[20]:***→Z[21]:***→Z[22]:***→Z[23]:IfEndJ=3=>***→Z[24]***→Z[25]:***→Z[26]:***→Z[27]:***→Z[28]:***→Z[29]:***→Z[30]:***→[31]:IfEnd………………………………………………..Goto HLbl F: J=1=>***→Z[8]:***→Z[9]:***→Z[10]:***→Z[11]:***→Z[12]:***→Z[13]:***→Z[14]:*** →Z[15]：IfEndJ=2=>***→Z[16]:***→Z[17]:***→Z[18]:***→Z[19]:***→Z[20]:***→Z[21]:***→Z[22]:***→Z[23]:IfEndJ=3=>***→Z[24]***→Z[25]:***→Z[26]:***→Z[27]:***→Z[28]:***→Z[29]:***→Z[30]:***→[31]:IfEnd(注：如有多个曲线元要素继续添加入数据库DAT1中)Goto GLbl E:“X0=”?U:“Y0=”?V:“S0=”?O:“F0=”?G:“LS=”?H:“R0=”?P:“RN=”?R:“Q=”?Q Lbl G:1÷P→C:(P-R)÷(2HPR)→D:180÷π→E:Return7．隧道放样子程序：SDFY5.5→A:8.5→B:L+1.5→CIf I=1:Then -1→K:Else 1→K:IfEnd 给左右隧道符号赋值1.25K→F:4.25K→L:7.25K→K 给各圆心距路线的偏距赋值Z[2]-C→D 计算测点至圆心的高差D，If D>A:Then Goto 1:Else If D<0:Then Goto 2:IfEnd:IfEndZ-L→E:Abs(E)→E:√(D2+E2)→W:W-A→T: 计算测点至圆心的平距E，实测半径W“R(m)=”:W◢“DR(m)=”:T◢If E≤A:Then D-√(A2-E2)→C:E-√(A2-D2)→L 计算测点水平偏移量L、垂直偏移量C “Dh(m)=”:C◢显示测点垂直偏移量C“DL(m)=”:L◢显示测点水平偏移量LElse E-√(A2-D2)→L“DL(m)=”:L◢IfEnd:ReturnLbl 1Z-L→E:Abs(E)→E:√(D2+E2)→W:W-A→T: 计算测点至圆心的平距E，实测半径W“R(m)=”:W◢“DR(m)=”:T◢If E≤A:Then D-√(A2-E2)→C“Dh(m)=”:C◢IfEnd:ReturnLbl 2:Abs(D)→D:Abs(Z)→Z:Abs(K)→K:Abs(F)→FIf Z>K:Then Z-F→E:Else Z+K→E:IfEnd√(D2+E2)→W:W-B→T“R(m)=”:W◢“Rc(m)=”:T◢If E≤B:Then D-√(B2-E2)→C:E-√(B2-D2)→L 计算测点水平偏移量L、垂直偏移量C “Dh(m)=”:C◢显示测点垂直偏移量C“DL(m)=”:L◢显示测点水平偏移量LElse E-√(B2-D2)→L“DL(m)=”:L◢显示测点水平偏移量LIfEnd:Return8．曲线元要素判断数据库：DAT2(两条线路，可扩充为多条线路)If I=2:Then Goto A:IfEndIf S<第一线元终点里程:Then 1→J:Prog “DAT1”:ReturnElse If S<第二线元终点里程:Then 2→J:Prog “DAT1”:ReturnElse If S<第三线元终点里程:Then 3→J:Prog “DAT1”:ReturnElse If S<第四线元终点里程:Then 4→J:Prog “DAT1”:ReturnIfEnd:IfEnd:IfEnd:IfEndIf S≥第四线元起点里程And S<第五线元终点里程：Then 5→J：Prog “DAT1”：Return .................................Else If S<第n线元里程：Then n→J：Prog “DAT1”：ReturnIfEnd:IfEnd:IfEnd:IfEndLbl AIf S<第一线元终点里程:Then 1→J:Prog “DAT1”:ReturnElse If S<第二线元终点里程:Then 2→J:Prog “DAT1”:ReturnElse If S<第三线元终点里程:Then 3→J:Prog “DAT1”:ReturnElse If S<第四线元终点里程:Then 4→J:Prog “DAT1”:ReturnIfEnd:IfEnd:IfEnd:IfEndIf S≥第四线元起点里程And S<第五线元终点里程：Then 5→J：Prog “DAT1”：Return Else If S<第n线元里程：Then n→J：Prog “DAT1”：ReturnIfEnd:IfEnd:IfEnd:IfEnd .................................9．高程数据库子程序：DAT3（两条线路，可扩充为多条线路）If I=0:Then Goto A:Else If I=2:Then Goto B:IfEnd:IfEndIf S≤46500:Then 46150→A:361.26→B:40000→R:-0.3→D:-1.2→E:ReturnElse If S≤47200:Then 46800→A:353.46→B:20000→R:-1.2→D:1→E:ReturnElse If S≤48400:Then 47500→A:360.46→B:20000→R:1→D:-1.6→E:ReturnElse If S≤49600:Then 48700→A:341.26→B:10000→R:-1.6→D:2.45→E:Return IfEnd:IfEnd:IfEnd:IfEndIf S>49600 And S≤50300:Then 49840→A:369.19→B:100000→R:2.45→D:2.8→E:Return Else If I=2:Then49800→A:368.21→B:150000→R:2.45→D:2.7→E:ReturnIfEnd:IfEndLbl A:“C PEG(m)=”?→A“C H(m)=”?→B“R(m)=”?→R“I1=”?→D“I2=”?→EReturnLbl B49800→A:368.21→B:150000→K:2.45→D:2.7→E:Return10.填方路基宽度检查子程序（WIDE）“WIDTH CHECK”◢“H→MID(m)=”?→A 起坡点与中桩高差常数“DIST→MID(m)=”?→B 起坡点与中桩距离常数“i=”?→C 路基填筑坡度CAbs(Z[2]-L+A)→DAbs(Z)-B-D→E“DL(m)=”:E◢显示实有宽度与应有宽度之差值Return11.坡口坡脚线坐标计算子程序:SUB-TW“TIANWA LINE”◢“DIST→MID(m)=”?A 起坡点与中桩距离常数“H→MID(m)=”?B 起坡点与中桩高差常数“T(1) Or W(Else)=”?C 填方边坡为1，其它为挖方边坡“i=”?E 填方或挖方边坡坡度Abs(Z)-A→D: DE→KIf C=1:Then L-B-K→W:Else L+B+K→W:IfEnd“JS H(m)=”:W◢计算出应有高程“SC H(m)=”？→D 输入实测高程D-W→T“GC(m)=”:T◢显示实测高程与计算高程的差异Return三、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
CASIO fx-5800P线元法坐标正反算程序
CASIO fx-5800P线元法坐标正反算程序说明：本程序适用于卡西欧计算器 CASIO fx-5800P,可对全线贯通坐标正反算、竖曲线高程计算。
该程序可计算任意线型，包含（直线、圆曲线、缓和曲线、卵形曲线）等,还可以能通过坐标反推该点里程和距中线距离,适用测量员专用。
主程序名：ABCYT第1步Deg：Fix 3：10→DimZ第2步Lbl 3：＂1.DK=>XY＂：＂2.XY=>DK＂：＂Q＂?W：＂DK＂?S：Prog＂ABCYTSJ＂：If P=0：Then 10^(45)→P：IfEnd：If R=0：Then 10^(45)→R：IfEnd第3步1÷P→C：(P-R)÷(2HPR)→D：180÷π→E：If W=1：Then Goto 1：Else Goto2：IfEnd第4步Lbl 1：＂W＂?Z：＂α＂?N：Abs(S-O)→W：Prog ＂ABCYTZ＂第5步Cls：＂F=＂：Locate 3,1,F°：＂X=＂：Locate 3,2,X：＂Y=＂：Locate 3,3,Y◢第6步Prog＂ABCYTSQX＂：Cls：＂H=＂：Locate 3,1,H◢第7步1→W：90→N：Goto 3第8步Lbl 2：?X：?Y：X→I：Y→J：Prog＂ABCYTF＂：O+W→S第9步Cls：＂K=＂：Locate 3,1,S：＂S=＂：Locate 3,2,Z◢第10步2→W：Goto 3正算子程序名：ABCYTZ第1步0.1739274226→A：0.3260725774→B：0.0694318442→K：0.3300094782→L第2步1-L→F：1-K→M第3步U+W×(A×cos(G+Q×E×K×W×(C+K×W×D))+B×cos(G+Q×E×L×W×(C+L×W×D))+B×cos(G+Q×E×F×W×(C+F×W×D))+A×cos(G+Q×E×M×W×(C+M×W×D)))→X第4步V+W×(A×sin(G+Q×E×K×W×(C+K×W×D))+B×sin(G+Q×E×L×W×(C+L×W×D))+B×sin(G+Q×E×F×W×(C+F×W×D))+A×sin(G+Q×E×M×W×(C+M×W×D)))→Y第5步G+Q×E×W×(C+W×D)→F：F+ N→Z[1]第6步X+Z×cos（Z[1]）→X：Y+Z×sin（Z[1]）→Y反算子程序名：ABCYTF第1步Lbl 2：（S-O）→W：0→Z：Prog ＂ABCYTZ＂：F-90→Z[9]：(J-Y)×cos(Z[9])-(I-X)×sin(Z[9])→Z[10]第2步If Abs(Z[10])>0.001：Then S+Z[10]→S：Goto 2：Else Goto 1：IfEnd第3步Lbl 1：(Y-J)÷sin(Z[9])→Z第4步Pol(X-I，Y-J)：If Z<0：Then -1×I→Z：Else 1×I→Z：IfEnd数据库名：ABCYTSJ第1步Goto1第2步Lbl 1：If S<7586.707 Or S>13346.96：Then Cls：Locate 2,2,＂PQX＂：Locate 4,3,＂CHAOXIAN＂：Locate 10,4,＂→Stop＂◢第3步Stop：IfEnd第4步Lbl 1：IF S<7946.707：Then 98°56′56″→G：7586.707→o：3378605.445→U：453648.704→V：0→P：4500→R：360→H：1→Q：Return：IfEnd第5步Lbl 1：IF S<11766.03：Then 101°14′26″→G：7946.707→o：3378544.714→U：454003.518→V：4500→P：4500→R：3819.323→H：1→Q：Return：IfEnd第6步Lbl 1：IF S<12126.03：Then 149°52′11″→G：11766.03→o：3376389.890→U：457018.324→V：4500→P：0→R：360→H：1→Q：Return：IfEnd第7步Lbl 1：IF S<13346.96：Then 152°09′41.68″→G：12126.03→o：3376073.846→U：457190.654→V：0→P：0→R：1220.93→H：0→Q：Return：IfEnd第n步……………………………………………………数据输入说明：第1步Goto1第2步Lbl 1：If S<本条线路起点里程S>本条线路止点里程Cls：Locate 2,2,＂PQX＂：Locate 4,3,＂CHAOXIAN＂：Locate 10,4,＂→Stop＂◢第3步Stop：IfEnd第4步Lbl 1：If S<本线元止点里程：Then线元起点切线方位角→G：线元起点桩号→O：线元起点坐标X→U：线元起点坐标Y→V：线元起点半径（直线为0、曲线为半径）→P：线元止点半径（直线为0、曲线为半径）→R：线元长度→H：线元转向（左转为-1、右转为1、直线为0）→Q：Return：IfEnd第n步Lbl 1：每增加一行则为增加一个线元要素。
曲线(含直线)任意里程中边桩坐标正反算
曲线（含直线）任意里程中边桩坐标正反算（CASIO fx-5800）J-PQX(平曲线数据输入,自动切换到J-JSMS)“JD”?A:“JDX”?B:“JDY”?C:“FJ”?N:“ZJ：Z-，Y+”?O:"R"?R:“LS1”?D:“LS2”?K：DD÷(24R)-D^4÷(2688RRR)->Z[1]:D÷2-D^3÷(240RR)->Z[2]：(DD-KK)÷(24R)÷sin(Abs(O))->Z[37]：“T1=”:(R+Z[1])tan((AbsO÷2))+Z[2]-Z[37]->Z[3]▲“T2=”:（(R+KK÷(24R))-K^4÷(2688RRR)）tan(Abs(O÷2))+K÷2-KKK÷(240RR)+Z[37]->Z[4]▲"L=":Abs(O)πR÷180+(D+K)÷2->L▲tan-1((R+Z[1])÷(Z[3]-Z[2]))->J：“E=”:(R+Z[1])÷sin(J)-R->Z[34]▲A-Z[3]->Z[35]：X+E->Z[36]：“ZH=”:A-Z[3]->Z[35]▲“HY=”:Z[35]+D->Z[36]▲“QZ=”:Z[36]+(L-K-D)÷2->Z[31]▲Z[35]+L-K->Z[32]：X+L-> Z[31]:“YH=”:Z[35]+L-K->Z[31]▲“HZ=”:Z[35]+L->Z[32] ▲Prog“J-JSMS”J-JSMS(放样模式主程序)“1-ZS,2-FS,4-DMFY”Lb1 0：“MS=”?Z:If Z≤1:Then Goto 1Else if Z ≤ 4 :Then Goto 2Lb1 1："P"?P:“BZ”?V:“BJ”?W：Prog“JP”：“X=”:Z[15]->X▲“Y=”:Z[16]->Y ▲Goto 0Lb1 2：“X”?X:“Y”?Y：“BJ”?W: Prog“JF”“P=”P->P▲“BZ=”:V-2.42->V▲“G1=”:1.225.353+（P-O）*1.5/1000-> [58] (高程计算公式）▲“S=”?SProng“J”Goto 0JP(平曲线正算子程序)FixmLb1 1：N->J：B-Z[3]cosN->Z[21]：C-Z[3]sinN->Z[22]If P≤A-Z[3]+D:Then P-A+Z[3]->I：90II÷(RDπ)->H：O＜0=＞-H->H:H+W+N->HElse If P≤A-Z[3]+L-K:Then P-(A-Z[3]+D)->I： Goto 4Else Goto 2Lb1 2：B+Z[4]cos（J+O）->Z[21]：C+Z[4]sin（J+O）->Z[22]：J+O+180->J: A-Z[3]+L-P->I：90II÷(RDπ)->H：O＞0=＞-H->H:H+J+W+180->H： Goto 3Lb1 3： I-I^5÷(40RRKK)->M：III÷(6KR)-I^7÷(336RRRKKK)->T： Goto 5Lb1 4：(D+2I)×90÷(πR)->H：RsinH+Z[2]->M：R(1-cosH)+Z[1]->T：O＜0=＞-H->H:N+W+H->H： Goto 5Lb1 5：If P≤A-Z[3]+L-K:Then O＜0=＞-T->T:Goto 6Else 0>0=＞-T->TGoto 6Lb1 6：Z[21]+Mcos J-Tsin J+Vcos H->Z[15]：Z[22]+Msin J+T cosJ+Vsin H-> Z[16]JF(平曲线反算子程序)0->V:N-90->JA+（Y-C)cosJ-(X-B)sin J->PLb11:Prog"JP":H-180->Z[10]:(Y-Z[16])cos(Z[10])-(X-Z[15])sin(Z[10])->I:If Abs(I)＜1*-4:Then Goto 210Else P+I->P:Goto1Lb1 2:(Y-Z[16])/sinH->VJ(超欠挖算子程序)If S>G-1.601:Then “YY=”:∫(2.022-(G+1.108-S)2)+3.3-∫(B2)▲Else “YG=”:∫(13.912- B2)-(G+11.54-S▲If End一、程序简介1、本套程序共有2个主程序，5个子程序。
卡西欧fx5800p计算器 坐标正反算程序
(以下程序是专业人士编写，本店铺不对程序负责，仅供您参考使用。
)卡西欧fx5800p计算器坐标正反算程序一、程序功能本程序由 6 个主程序、 5 个次子程序及 5 个参数子程序组成。
主要用于公路测量中坐标正反算，设计任意点高程及横坡计算 , 桥涵放样，路基开挖口及填方坡脚线放样。
程序坐标计算适应于任何线型 .二、源程序1. 主程序 1 ：一般放样反算程序(① 正算坐标、放样点至置仪点方位角及距离；② 反算桩号及距中距离 )程序名 :1ZD-XYLb1 0:Norm 2F=1 ： ( 正反算判别， F=1 正算， F=2 反算 , 也可以改 F 前加？，改 F 为变量 )Z[1]=90 （与路线右边夹角）Prog ＂ THB ＂： F=1=>Goto 1:F=2=>Goto 2Lb1 1: Ｆ ix 3: ＂Ｘ = ＂： Locate 6,4, Ｘ◢＂Ｙ＝＂： Locate 6,4, Ｙ◢Ｐ rog ＂3JS”:Goto 0:Lb1 2:Fix 3: ＂ＫＭ＝＂： Locate 6,4, Ｚ◢＂Ｄ＝＂： Locate 6,4, Ｄ◢Ｇ oto 0２．主程序２：高程序横坡程序 ( 设计任意点高程及横坡 )程序名： 2GCLbI 0:Norm 2“KM”?Z:?D:Prog”H”:Fix 3:” H=”:Locate 6,4,H◢“ I=”: Locate 6,4,I◢Goto 03. 主程序 3 ：极坐放样计算程序 ( 计算放样点至置仪点方位角及距离 )程序名： 3JSX ： Y ：1268 ．123→K( 置仪点 X 坐标 )2243 ．545→L （置仪点 Y 坐标，都是手工输入 , 也可以建导线点数据库子程序 , 个人认为太麻烦）Y-L→E ： X-K→F ： Pol(F,E):IF J<0:ThenJ+360→J:Int(J)+0.01Int(60Frac(J))+0.006Frac(60Frac(J)) →J:( 不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢ 来直接显示) Fix 4:” FWJ=”: Locate 6,4,J◢( 不习惯小数点后四位为角度显示的，也可以用命令J◢DMS◢ 来直接显示 )Fix 3:” S=”:Locate 6,4,I◢4 ．主程序 4 ：涵洞放样程序（由涵中心桩号计算出各涵角坐标、在主程序 3 中输入置仪点坐标后计算放样点至置仪点方位角及距离 )程序名： 4JH-XYLbI 0:Norm 290→Z[1]( 涵洞中心桩与右边夹角，手工输入，也可以修改成前面加？后变为变量 )1→F:Prog”THB”:?L:Z[2]-Z[1] →E:X+Lcos(E) →X:Y+Lsin(E) →Y:Fix 3: ＂Ｘ = ＂： Locate 6,4, Ｘ◢＂Ｙ＝＂： Locate 6,4, Ｙ◢Ｐ rog ＂3JS”:Goto 0:5 ．主程序 5: 路基开挖边线及填方坡脚线放样程序（输入大概桩号及测量坐标、地面标高计算出偏移距离、桩号、距中距离、填挖高度）程序名： 5FBXLbI 0:Norm 2: 18→DimZ:2→F:90→Z[1]:Prog“THB”:Z:D:”M0”?M:M→Z[4]:D→Z[3]:Prog”6GD”:L→Z[6]:If D<0:Then 0.75-L→D:Goto H:Else L-0.75→D:Goto H:IfEndLbI H:Prog”H”:H-0.03-Z[4] →Z[5]:Z[6] →L:If Z[5]<0:Then –z[5]→G:Goto W:Else Z[5] →G:Goto T:LbI W:Prog “W0”:Z[10]+Z[11] →A: If G>A:Then Goto 1:Else IfG>Z[10]:Then Goto 2:Else Goto 3:IfEnd:LbI 1:L+Z[12]+Z[13]+Z[14]+(G-A)×Z[9]+Z[11]×Z[8]+Z[10]×Z[7]:Goto Z: LbI 2:L+Z[12]+Z[13]+(G-Z[10])×Z[8]+Z[10]×Z[7]:G oto Z:LbI 3:L+Z[12]+G×Z[7]:Goto z:LbI T:L+0.5→N:If G>Z[17]:Then(N+Z[18]+(G-Z[17])×Z[16]+Z[17]×Z[15])→S:Goto Z:Else(N+G×Z[15])→S:Goto z:LbI Z:Z[3]→D:Fix 2:Abs(D)-S→T:” L0=”:L Locate 6,4,T◢＂ＫＭ＝＂： Locate 6,4, Ｚ◢＂Ｄ＝＂： Locate 6,4, Ｄ◢“ TW=”: Locate 6,4,Z[5]◢Goto 06 ．主程序 6 ：路基标准半幅宽度计算程序 ( 对于设计有加宽渐变的有用，如路基宽度无变化，则把此程序直接输入半幅宽度值至 L)程序名： 6GDProg “G0”Z-C→E:(B-A)×E/S+A→L:L:7 ．坐标计算次程序（ＴＨＢ）程序名： THB18→DimZ ：＂ＫＭ＂ ?Z ： Prog ＂Ｘ０＂1÷P→ Ｃ： (P-R)÷(2HPR) → Ｓ：180÷π→ Ｅ：Ｆ =1=>Goto 1 ：Ｆ=2=>Goto 2←┘Lbl 1 ： ? Ｄ： Abs( Ｚ－O) →W ： Prog " Ａ " ： X ：Ｙ：Ｇ oto 3 LbI 2: Ｘ：Ｙ：Ｘ→I ：Ｙ→J ： Prog ＂Ｂ＂：Ｏ＋Ｗ→Z ：Ｄ→ Ｄ：Ｇ oto 3LbI 3:IF Ｆ＝１ Then Ｘ：Ｙ： Else Ｚ：Ｄ8. 正算子程序 ( Ａ )程序名： A0.1184634425→A ：0.2393143352→B ：0.2844444444→N 0.046910077→K ：0.2307653449→L ：0.5→M ： U+W(Acos(G+QEKW(C+KWＳ ))+Bcos(G+QELW(C+LW Ｓ ))+Ncos(G+QEMW(C+MWＳ ))+Bcos(G+QE(1-L)W(C+(1-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS))) →X ：V+W(Asin(G+QEKW(C+KW Ｓ ))+Bsin(G+QELW(C+LWＳ ))+Nsin(G+QEMW(C+MWS))+Bsin(G+QE(1-L)W(C+(1-L)WS))+Asin(G+QE(1-K)W (C+(1-K)W Ｓ))) →Y ：G+QEW(C+W Ｓ )+ Ｚ［１］→ Ｚ［２］： X+ Ｄ cos （Ｚ［２］）→ Ｘ： Y+ Ｄ sin （Ｚ［２］）→ Ｙ9. 反算子程序 ( Ｂ )程序名： BG-90→ Ｔ： Abs((Y-V)cos(T)-(X-U)sin(T)) → Ｗ：0→ Ｄ： Lbl 0 ： Prog " Ａ " ： T+QEW(C+W Ｓ) → Ｌ： (J-Y)cos(L)-(I-X)sin(L)→ Ｄ： IF Abs( Ｄ )<0.01:Then Goto1 ： Else W+ Ｄ→W:Goto 0←┘L bl 1:0→D ： Prog " Ａ " ： (J-Y)÷sin( Ｚ［２］) →D:10 ．高程计算子程序（ H ）程序名： HProg “S0”:R:T:C:G:I:C-T→F:Z-F→L:C+T→E:G-TI→Q:If T=O:ThenQ+LI→H:Goto 0:Else If Z<F:Then Q+LI→H:Goto 0:Else If Z≤E:ThenQ+LI+L2÷2÷R→H:Goto 0:LbI 0:H:If D=0:Then Goto I:Else Prog “I”:H+V→H:Goto I:LbI I:H:I: 11. 高程超高计算程序（ I ）程序名： IProg”I0”:W=1=> Goto 0:W=2=>Goto 1:LbI 0:If L=0:Then Abs(D)×M→V:Goto 2:ElseAbs(D)×((N-M)×(Z-C)÷L+M)→V:Goto 2:IfEnd:LbI 1:If L=0:Then Abs(D)×M→V:Goto 2:ElseAbs(D)×(((3((Z-C)÷L)2-2((Z-C)÷L)∧(3))×(N-M))+M)→V:Goto 2:IfEnd: LbI 2:Abs(D)→E:V÷E→I:I(E-K)→V:12 ．数据子程序 ( 附后示例 )① 程序名： X0 （坐标计算要素程序）If Z≥25900 And Z≤26615.555:Then25900→O:11587.421→U:1847.983→V:101 。
坐标反算里程偏距RTk
坐标反算里程偏距RTk
今天就来看看RTK单点校正的操作步骤：单点校正：利用一个点的WGS84坐标和当地坐标可以求出3个平移参数，旋转为零，比例因子为1。
在不知道当地坐标系统的旋转、比例因子的情况下，单点校正的精度无法保障，控制范围更无法确定。
因此建议尽量不要使用这种方式。
一．基准站架设基准站架设的好坏，将影响移动站工作的速度，并对移动站测量质量有着深远的影响，因此用户注意使观测站位置具有以下条件：在15度截止高度角以上的空间部应没有障碍物；邻近不应有强电磁辐射源，比如电视发射塔、雷达电视发射天线等，以免对RTK电信号造成干扰，离其距离不得小于200m；基准站最好选在地势相对高的地方以利于电台的作用距离；地面稳固，易于点的保存。