公路坐标计算程序
公路任意点坐标计算公式,基本上包含公路中所有的线形,满足你的好奇心!
公路道路上有关P 点的坐标计算：在公路施工过程中，需要进行放样的点位，不外乎两种情况：一种是该点位于公路中心线上，即公路中桩，另一类则是点位在中线以外，位于某个中桩的横断方向上。
这样无论哪种情况，需要放样的点的桩号首先是已知的。
以下就这两种情况，对公路点位放样计算进行一下阐述，讲述一下坐标计算方法。
一：P 点位于直线段上，各桩坐标计算： 1、 P 点在直线上各中桩坐标计算当需要放样的P 点位于直线上时，有两种情况：位于YZ 到ZY 或者HZ 到ZH 之间， 或者位于公路QZ 和ZH （ZY ）之间，其计算方法相同，公式如下：[公式（1）]X p =X 0+l cosA i-1,i Y p =Y 0+l sin A i-1,I式中， （X 0 ,Y 0） 为该段直线的起点（可以是YZ ，HZ 或QZ ）坐标 l 为要求计算的P 点与该直线段起点的桩号差（距离）。
2、 P 点位于横断面上，其所对应的中桩位于直线上时：X p =X z +Dcos （A i-1,i ±90） Y p =Y z +Dsin （A i-1,i ±90）式中， （X z ，Y z ）为P 点对应的中桩的坐标 P 点位于左幅时，取“—”反之取“+” D 为P 点到直线上的法线距离二、P 点位于单圆曲线上，各桩坐标计算：1、当需要放样的P 点位于单圆曲线上，其中桩坐标计算如下：[公式（2）]式中， （X 0，Y 0）为ZY 点坐标，R 为圆曲线半径 l 为P 点与ZY 点的桩号差（弧长） 当路线左转时，取“—”，反之取“+”2、P 点位于横断面上，其所对应的中桩位于单圆曲线上时：式中，第一个“”号，路线左转取“—”，右转取“+”第二个“”号，P 点位于左幅时，取“—”，反之取“+”三、P 点位于带缓和曲线的圆曲线上，各桩坐标计算：当P 点位于带缓和曲线的圆曲线时，分为以下三种情况： 第一种情况，ZH 到HY 段，中桩和边桩计算： 1、ZH 到HY 段，中桩坐标计算：[公式（3）]式中，c = l -(X 0，Y O )为ZH 点坐标l 为P 点与ZH 点桩号差，L s 为缓和曲线长 当路线左转时，取“—”，反之取“+”2、ZH 到HY 段，P 点对应的中桩位于带缓和曲线的圆曲线上：式中，（Xz ，Yz ）为P 点对应的中桩坐标l 为P 点对应的中桩与ZH 点桩号差，Ls 为缓和曲线长 第一个“”号，路线左转取“—”，右转取“+” 第二个“”号，P 点位于左幅时，取“—”，反之取“+” 第二种情况，HY 到YH 段，中桩和边桩计算：1、HY 到YH 段，中桩坐标计算：[公式（4）]式中，(X 0，Y O )为HY 点坐标l 为P 点与HY 点桩号差，Ls 为缓和曲线长 当路线左转时，取“—”，反之取“+”2、HY 到YH 段，P 点对应的中桩位于带缓和曲线的圆曲线上：式中，（Xz ，Yz ）为HY 点坐标l 为P 点对应的中桩与HY 点桩号差，Ls 为缓和曲线长 前两个“”号，路线左转取“”，右转取“” 第三个“”号，P 点位于左幅时，取“—”，反之取“+”第三种情况，YH 到HZ 段，中桩和边桩计算： 1、YH 到HZ 段，中桩坐标计算：[公式（5）]c = l -(X 0，Y O )为HZ 点坐标l 为HZ 点与P 点桩号差，L s 为缓和曲线长 当路线左转时，取“+”，反之取“—”2、YH 到HZ 段，P 点对应的中桩位于带缓和曲线的圆曲线上：式中，（Xz ，Yz ）为P 点对应中桩坐标l 为HZ 点桩号与P 点对应的中桩桩号差，Ls 为缓和曲线长 第一个“”号，路线左转取“+”，右转取“—” 第二个“”号，P 点位于左幅时，取“—”，反之取“+”四、复曲线上各点的坐标计算：1、 当复曲线中间不设缓和曲线时，采用以下方法进行计算：对于第一缓和曲线、第一段圆曲线以及第二缓和曲线，分别用公式（3）、公式（4）和公式（5）计算；对于第二段圆曲线，用公式（2）计算，计算时将公式（2）中的换成，l 1，L s1分别为第一圆曲线和第一缓和曲线长度，左转取“—”，右转取“+”。
高速MC公路测量CASIO4800&4850万能坐标计算程序(完整版)
高速公路测量CASIO4800&4850万能坐标计算程序(完整版）程序特点：真正的全线贯通坐标正反计算、任意斜角计算！！！程序中加入测站点，真正的实现了“坐标法”与“极坐标法”两种放样方法的同时显示的功能，使得放样操作方法选择时更加灵活！！！在曲线元要素输入时仅需要输入第一段全部曲线元要素，后面曲线元要素除起点半径、终点半径、曲线长、转向需输入外其他要素均从前一曲线按辛普森8等分计算得出，解决了主线坐标计算无法获得第二段及其以后曲线元起点参数的问题；辛普森公式任意等分，满足所有精度要求；全线曲线元数据一次性程序化输入，参数存储采用扩充变量数据库，无需修改程序内容；多功能采用单程序编程，避免频繁调用子程序，提高运算速度。
一、程序：ZBJSW“1.ZS 2.FS 3.SZ”:W=1=>Z[2]=0:V=0：Goto 1 ΔW=2=> Goto 4ΔW=3=> O “KOU LING”:O≠123456=>O=0: “OUT”◢Goto CΔO=0: V=0:Z[1]=0:Goto 0←┘Lbi 0←┘”N0.”：Z[1]+1 ◢Z[1]=0=>{ABCREFGUKO}:A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G “LS”:U“G”:K“X（00）”: O“Y（00）”: Z[Z[1]×8+3]=A:Z[Z[1]×8+4]=B:Z[Z[1]×8+5]=C:Z[Z[1]×8+6]= R-1:Z[Z[1]×8+7]= E-1:Z[Z[1]×8+8]=F: Z[Z[1]×8+9]=F+G: Z[Z[1]×8+10]=U:“NEXT”◢Isz Z[1]: Goto 0ΔZ[1]=1=>D=Z[9]:Z=0:Z[2]=0:GOTO 2ΔD=Z[(Z[1]-1)×8+9]:Z=0:Z[2]=Z[1]-1:GOTO 2←┘Lbi A←┘Z[Z[1]×8+3]=X:Z[Z[1]×8+4]=Y:Z[Z[1]×8+5]=J: Z[Z[1]×8+8]=D: {REGU}:R“R0”：E “RN”: G“LS”:U“G”: Z[Z[1]×8+6]=R-1：Z[Z[1]×8+7]=E-1: Z[Z[1]×8+9]=D+G: Z[Z[1]×8+10]=U:“NEXT”◢Isz Z[1]: Goto 0←┘Lbi 1←┘{DZT }:D：Z：T“RJ”：Z[2]=0:Goto 2←┘Lbi 2←┘V≠1=>Z[2]>Z[1] =>GoToCΔΔD≤Z[Z[2]×8+9]=> A=Z[Z[2]×8+3]:B=Z[Z[2]×8+4]: C =Z[Z[2]×8+5]:R=Z[Z[2]×8+6]: E=Z[Z[2]×8+7]: F=Z[Z[2]×8+8]: G=Z[Z[2]×8+9]: U=Z[Z[2]×8+10]: Goto3ΔIsz Z[2]:Goto 2←┘Lbi 3←┘W=3 =>N=8：≠P＝U（E-R）÷Abs（G-F）：Q＝Abs（D-F）÷N：S=90Q÷π：J＝C+(NPQ+2UR)NS：L=1←┘X＝A+Q÷6×(Cos C+Cos J +4∑（Cos (C+((L+0.5)PQ+2UR)×(L+0.5)S),L,0,(N-1)）+2∑（Cos (C+((LPQ+2UR)LS,L,1,(N-1)))+ZCos（J+ T）←┘Y＝B+Q÷6×(Sin C+Sin J +4∑（Sin (C+((L+0.5)PQ+2UR)×(L+0.5)S),L,0,(N-1)）+2∑（Sin (C+((LPQ+2UR)LS,L,1,(N-1)))+Z Sin（J+T）:V=1=>Goto6ΔV=2=>Goto9ΔV=3=> GOTO CΔW=3=>GOTO AΔZ=0=>“X（Z）=”:X:Pause 0: “Y（Z）=”:Y◢Pol（（X－K），（Y－O））←┘“S（Z）=”:I ◢J<0=> J=J+360Δ“F（Z）=”: J→DMS◢Goto 1ΔZ<0=>“X（L）=”:X:Pause 0: “Y（L）=”:Y◢Pol（（X－K），（Y－O））←┘fx4850①“S（L）=”:I ◢J<0=> J=J+360Δ“F（L）=”: J→DMS◢Goto 1ΔZ>0=>“X（R）=”:X:Pause 0: “Y（R）=”:Y ◢Pol（（X－K），（Y－O））←┘“S（R）=”:J ◢J<0=> J=J+360Δ“F（R）=”: J→DMS◢Goto 1 ←┘Z=0=> X “X（Z）=”◢Y “Y（Z）=”◢Pol（（X－K），（Y－O））←┘I“S（Z）=”◢J<0=> J=J+360ΔJ“F（Z）=”◢Goto 1ΔZ<0=> X “X（L）=”◢Y “Y（L）”◢Pol（（X－K），（Y－O））←┘fx4800②I“S（L）=”◢J<0=> J=J+360ΔJ“F（L）=”◢Goto 1ΔZ>0=> X “X（R）=”◢Y “Y（R）=”◢Pol（（X－K），（Y－O））←┘I“S（R）=”◢J<0=> J=J+360ΔJ“F（R）=”◢Goto 1 ←┘Lbi 4←┘{MH} ：M“X”：H“Y”：Z[2]=0：GOTO 5←┘Lbi 5←┘V=1：D= Z[Z[2]×8+9]：Z=0：T=90：GOTO 2←┘Lbi 6←┘K=（（H -B）Cos（C-90）-（M-A）Sin（C-90））×（（H -Y）Cos（J-90）-（M-X）Sin（J-90））：K≤0=> Goto 7ΔIsz Z[2]：Goto5←┘Lbi 7←┘D=F+Abs((H -B)Cos(C-90)-(M-A)Sin（C-90）)：D>G=> Isz Z[2]: Goto5ΔGoto 8←┘Lbi 8←┘V=2 ：GOTO 3←┘Lbi 9 ←┘K=(H -Y)Cos(J-90)-(M-X)Sin(J-90)：Abs K<（1÷E）^3=>Goto BΔD=D+K ：GOTO 8←┘Lbi B←┘V=3 ：Z=0：Goto 3←┘Lbi C←┘Z=(H-Y) ÷Sin(J+90)：“D”:D:Pause 0: “Z”: Z◢4850输出(Z=(H-Y) ÷Sin(J+90)：D“D”◢Z “Z”◢4800输出)GOTO 4←┘Lbi C←┘二、说明a、编制说明本程序是运用复化辛普生公式根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距，对该曲线段范围内任意里程中边桩坐标进行计算，以及对卡西欧扩充变量的灵活应用，实现了真正意义上的的全线贯通及曲线要素输入程序化（在不修改程序内容的情况下可通过运行程序输入任意多段曲线元要素）。
公路测量坐标计算步骤
复曲线可分解为ZH~HY、ZH~HY~YH、HY~YH~HZ、YH~HZ等可在缓和程序、圆曲线程序进行计算。
6、其它复合型曲线
复合型曲线主要由多个圆曲线、缓和曲线组合而成，分解后可分别在“圆曲线”、“完整缓和曲线”、“非完整缓和曲线”。
“完整缓和曲线”是指起点直缓点或缓直点已知的缓和曲线，“非完整缓和曲线”是指起点直缓点或缓直点未知，仅知道缓和曲线上的任意一点或两点的缓和曲线，计算时需推算缓和曲线起点。
1、C型曲线
C型曲线主要由两个同向的单交点基本型曲线组成，可分解为两个单交点基本型曲线分别进行计算；
2、S型曲线
S型曲线主要由两个反向的单交点基本型曲线组成，可分解为两个单交点基本型曲线分别进行计算；
3、卵形曲线
卵形曲线主要组合为ZH~ห้องสมุดไป่ตู้Y1~Y1H~HY2~Y2H~HZ，可分解为ZH~HY1、HY1~Y1H、Y1H-HY2、HY2~Y2H、Y2H~HZ五段进行计算，分别为第一缓和线段、圆曲线段、卵刑曲线段、圆曲线段、第二缓和线段进行计算。
主要组合为ZH~HY~YH~HZ，A1=A2(或Ls1=Ls2)，其中A为缓和曲线回旋参数，A*A=Ls*R。
4、单交点非对称型曲线
主要组合为ZH~HY~YH~HZ，A1≠A2(或Ls1≠Ls2)，其中A为缓和曲线回旋参数，A*A=Ls*R。
(三)复杂曲线
本程序所能计算的复杂曲线主要为C型曲线、S型曲线、卵型曲线、凸型曲线、复曲线及其它复合型曲线：
公路测量坐标计算说明
一、软件主要功能：
本程序适用于公路主线、立交匝道及铁路的施工放样工作，可根据设计给定参数精确进行公路、铁路的线路坐标计算。
公路逐桩坐标计算程序
2 )sinαAB +(
P 点为顺时针方向时，其方位角为 αAB+900 P 点为逆时针方向时，其方位角为 αAB−900
第一段缓和曲线时：以直缓点（ZH）为起点计算，αAB 为 ZH 点的坐标方位角，L 为 P 点距 ZH 点的距离。 第二段缓和曲线时：以缓直点（HZ）为起点计算，αAB 为 HZ 点坐标方位角的反 方向即 HZ 点方位角加 180 度，L 为 P 点距 HZ 点的距离，加减 90 度刚好与第一 段缓和曲线相反。图如下：
缓和曲线转角公式：β
=
L2
2RLs
（2） 边桩坐标计算公式
左侧
XA=XP+T1cos(αAB± β-900) YA=YP+T1sin(αAB± β-900)
右侧
XB=XP+T2cos(αAB± β+900) YB=YP+T2sin(αAB± β+900)
*第二段缓和曲线计算边桩坐标时， 注意加减 90 度时， 与第一段缓和曲线相反。
2、 缓和曲线段 （1） 中桩坐标计算公式
Xp=X1+(L − Yp=Y1+(L −
L5 40R 2 Ls L5 40R 2 Ls
2 )cosαAB+(
L3 6RLs L3 6RLs
− −
L7 336R 3 ������s 3 L7 336R 3 ������s 3
)sin(αAB±900) )sin(αAB±900)
右侧
XB=XP+T2cos(αAB± β+900) YB=YP+T2sin(αAB± β+900)
第二章
公路导线测量计算
1
公路导线测量为附合导线测量，按路线前进方向测量右角。β
CASio 9860公路三维坐标计算程序帮助
程序使用说明--CASIO 9860G系列公路三维坐标计算程序程序说明：本程序是一款运行用CASIO 9860G 系列上的公路三维坐标计算程序，程序约8KB大小，使用9860中的表文件作为程序数据，由于受9860主存大小限制，一次约能进行200个交点（交点法）或是约600个线元（线元法）的公路平面坐标和中桩高程，如使用内部后备存储器（1.5M）在需要时手动调入进行计算，或是使用SD卡为数据载体则计算只受这两个容量的限制。
程序需运行在OS2.0以上的系统上，如不是OS2.0的请自行升级成OS2.0。
程序功能:程序能选择交点法（一般适用用主线）或是线元法（一般适用于岔道）计算，他们有相同的功能：1、全线三维坐标正算,计算全线路上任一桩号偏距点三维坐标（支持斜交）2、全线三维坐标反算，实测一点坐标，计算该点所在的桩号及偏距、高程计算可选3、批量计算三维坐标，从指定的桩号向结束处或开始处按间距连续计算点的三维坐标，所有点存在表文件中，倒入EXCEL中即可生成点位坐标表4、结构物相对位置正反算，计算结构物前后左右的点在线路上的坐标或是实测一点后反算实测点到结构特轴线的前后左右（该程序需单独起动不在统一界面选择下）5、路面挂线测量，支持全站仪和水准仪操作，实测一点计算该点到指定的结构层的高度，用于指导结构层抄平工作6、任意断面隧道超欠挖程序，实测断面上任一点即可得到该点的超或是欠挖值，支持单独断面，或是将断面加入到线路中去,加入线路后，实测一点三维坐标，即可得出该点到断面轮廓线的距离，快速的指导施工放样。
（该程序需单独起动不在统一界面选择下）7、线路放样功能，放任一点设站后只需输入放样桩号和偏距即可得出放样数据8、所有使用子程序预留标准接口，方便自己程序进行调用，以增加现有功能9、放样功能，设站后，只需输入需要放点的桩号和偏距(斜交），即可得到方位和平距10、主点元素计算及导出功能（交点法），计算所线路上所有交点五大桩的坐标，桩号，可导出到EXCEL，即生成主点元素表11、测量中常用的坐标正、反算程序12、不需输入点数的多边形面积计算程序13、单线元坐标正反算程序以及常用功能示例程序14、锥坡放样功能，实测一点反算该点在水平方向和高程方向和设计值的偏差，支持斜交锥坡15、任意多边形面积解析法计算程序16、路线任意级数边坡放样程序，只需测定边坡上一点三维坐标，指示该点向设计坡面的移动值。
4800公路测量放样计算程序
4800公路测量放样计算程序CASIO fx-4800P计算器程序DA DI——N E (公路弯道坐标计算)K“JD”：R：S“LS”：F“PJ——L－R＋”：A“ZH——N”：B“ZH——E”：C“JD——N”：D“JD——E”M“ZH”=K－（R＋S2÷（24R））tan（Abs F÷2）－S÷2＋S^3÷（240R2）▲N“HY”=M＋S▲O“QZ”=M＋（πRAbs F÷180＋S）÷2▲Z[1]“YH”=M＋πRAbs F÷180▲Q“HZ”= Z[1]＋S▲T=K－M▲L=Q－M▲E=（R＋S2÷（24R））÷cos（F÷2）－R▲Lb1 0：{G}：G“ZHUANG HAO”：G＜O＝=＞Z=G－M≠=＞Z=Q－G：△{U}：U“I——B OUT＋IN－”Z＜S＋0.005＝=＞H=tan-1（24Z2S2R2－Z^6）÷（48Z^3R^3－6Z^4SR））：V=Z－Z^5÷（40S2R2）＋Usin H：W=Z^3÷（6SR）－Z ^7÷（336Z^3R^3）－Ucos H：≠=＞H=180（Z－S÷2）÷（πR）：V=Rsin H＋S÷2－S^3÷（240R2）＋Usin H：W=R－Rcos H＋S2÷（24R）－Ucos H：△F＜0＝=＞G＜O＝=＞X=V：Y=W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Tsin Abs F －Vsin Abs F＋Wcos F：△≠=＞G＜O＝=＞X=V：Y=－W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Vsin Abs F－Tsin Abs F－Wcos F：△△C=A＝=＞D≥B＝=＞P=90：≠=＞P=270：△≠=＞C＜A＝=＞P= tan-1（（D－B）÷（C－A））＋180：≠=＞P= tan-1（（D－B）÷（C－A））：△△I“COM——N”=A＋√（X2＋Y2）cos（P－tan-1（Y÷X））▲J“COM——E”=B＋√（X2＋Y2）sin（P－tan-1（Y÷X））▲G=G＋20Goto 0程序运行输入：JD？——交点桩号R？——圆曲线半径LS？——缓和曲线长度PJ——L－R＋？——偏角右偏为正ZH——N？——ZH点N坐标ZH——E？——ZH点E坐标JD——N？——JD点N坐标JD——E？——JD点E坐标输出：ZH= ——ZH点桩号HY= ——HY点桩号QZ= ——QZ点桩号YH= ——YH点桩号HZ= ——HZ点桩号T= ——切线长L= ——曲线长E= ——外矢距再输入：ZHUANG HAO？——输入任一点桩号I——B OUT＋IN－？——距中线距离（外侧为正）再输出：COM——N= ——计算出N坐标COM——E= ——计算出E坐标以此循环——再输入——再输出程序中各字母所代表的含义：A——ZH点N坐标B——ZH点E坐标C——JD点N坐标D——JD点E坐标E——外矢距F——偏角（右偏为正）G——弯道内任一点桩号H——弯道内任一点切线角I——计算出的N坐标J——计算出的E坐标K——交点桩号L——曲线长M——ZH点桩号N——HY点桩号O——QZ点桩号P——ZH--JD方位角Q——HZ点桩号R——圆曲线半径S——缓和曲线长度T——切线长U——距中线距离（外侧为正）V——支距W——支距X——弯道内坐标Y——弯道内坐标Z——任一点曲线长度Z[1]——YH点桩号。
5800计算器公路全线坐标正、反算计算程序
5800计算器公路全线坐标正、反算计算程序正算主程序(ZS)Lb1 0FIX4：?S：?Z：Prog “PM-SJ”：Abs(S-O) → W：Prog "SUB1"："XS="：X◢"YS="：Y◢F-90→F(需要时可以让他显示，不要时不必输入)：Prog“SQX-SJ”: Prog“SQX”：“H=”：H◢Goto 0反算主程序(FS)正算子程序(SUB1)：(P-R)÷(2HPR) → D：180÷π→ E：0.1739274226→ A：0.3260725774→ B：0.0694318442→ K：0.3300094782→ L：1-L→ F：1-K→ M：U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+F WD))+Acos(G+QEMW(C+MWD))) → X：V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+F WD))+Asin(G+QEMW(C+MWD))) → Y：G+QEW(C+WD)+90→ F：X+ZcosF→ X：Y+ZsinF→ Y反算子程序(SUB2)：Abs((Y-V)cosT-(X-U)sin（T）) → W：0→ Z：Lbl 0：Prog "SUB1"：T+QEW(C+WD) →L：(J-Y)cosL-(I-X)sinL → Z：ifAbsZ<1E-6：thenGoto1：ElssW+Z → W：Goto 0：IfEnd←┘：0→ Z：Prog "SUB1"：(J-Y)÷sinF→ Z子程序（平面线形数据库）PM-SJ1÷P→ CG-90→TLbl 1ifS ≥500（线元起点里程）AndS＜769.256（线元止点里程）：Then19942.837 →U （线元起点X坐标）：28343.561 →V（线元起点Y坐标）：500 →O（线元起点里程）：125°16°31°→G（线元起点方位角）：269.265 →H（线元长度）：10^（45）→P（线元起点曲率半径）：10^（45）→R（线元终点曲率半径）：0 →Q （线元左右偏标志：左负右正）：IfEnd←┘ifS ≥769.256（线元起点里程）AndS＜806.748（线元止点里程）：Then19787.34 →U （线元起点X坐标）：28563.378→V（线元起点Y坐标）：769.256→O（线元起点里程）：125°16°31°→G（线元起点方位角）：37.492→H（线元长度）：10^（45）→P（线元起点曲率半径）：221.75→R（线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘ifS ≥806.748（线元起点里程）AndS＜919.527（线元止点里程）：Then19766.566→U （线元起点X坐标）：28594.574→V（线元起点Y坐标）：806.748→O（线元起点里程）：120°25°54.07°→G（线元起点方位角）：112.779→H（线元长度）：221.75→P（线元起点曲率半径）：221.75→R（线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘ifS ≥919.527（线元起点里程）AndS＜999.812（线元止点里程）：Then19736.072 →U（线元起点X坐标）：28701.893 →V（线元起点Y坐标）：919.527→O（线元起点里程）：91°17°30.63°→G（线元起点方位角）：80.285→H（线元长度）：221.75→P（线元起点曲率半径）：9579.228→R（线元终点曲率半径）：-1→Q （线元左右偏标志：左负右正）：IfEnd←┘ifS ≥999.812（线元起点里程）：Then19744.038 →U（线元起点X坐标）：28781.659→V（线元起点Y坐标）：999.812→O（线元起点里程）：80°40°50°→G （线元起点方位角）：100→H（线元长度）：10^（45）→P（线元起点曲率半径）：10^（45）→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左负右正）：IfEnd←┘子程序（竖曲线数据库）SQX-SJifS≤999.812（竖曲线终点里程）：then0.0357（前坡度：上坡正下坡负）→E：0.0227（后坡度：上坡正下坡负）→F：600→R（竖曲线半径）：780→G（变坡点里程）：110→C（变坡点高程）：IfEnd←┘子程序（竖曲线计算公式）SQXE-F→J：ifJ＞0：then-R→R：ElssR→R：Abs（R*J÷2）→T：IfEnd←┘ifS≤G-T：thenG-S→L：C-LE→H：IfEnd←┘ifS＞G-T And S≤G：thenG-S→L：（S+T-G）²÷2÷R→N：C-LE+N → H：IfEnd←┘ifS＞G And S≤G+T：then S-G→L：（G +T- S）²÷2÷R→N：C+LF+N → H：IfEnd←┘ifS＞G+T ：thenS-G→L：C+LF → H：IfEnd←┘三、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
计算公路的方法
计算公路的方法一、程序功能本程序由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:Then J+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:IfEnd LbI 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 If G>Z[10]:Then Goto 2:Else Goto 3:IfEnd:LbI1: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]:Goto 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:F ix 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))+Asi n(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←┘Lbl 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:Then Q+LI→H:Goto 0:Else If Z<F:Then Q+LI→H:Goto 0:Else If Z≤E:Then Q+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:Else Abs(D)×((N-M)×(Z-C)÷L+M)→V:Goto 2:IfEnd:LbI 1:If L=0:Then Abs(D)×M→V:Goto 2:Else Abs(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:Then 25900→O:11587.421→U:1847.983→V:101。
卡西欧计算器坐标计算程序
CASIO fx-4800P计算器程序DA DI——N E (公路弯道坐标计算)K“JD”：R：S“LS”：F“PJ——L－R＋”：A“ZH——N”：B“ZH——E”：C“JD——N”：D“JD——E”M“ZH”=K－（R＋S2÷（24R））tan（Abs F÷2）－S÷2＋S^3÷（2 40R2）▲N“HY”=M＋S▲O“QZ”=M＋（πRAbs F÷180＋S）÷2▲Z[1]“YH”=M＋πRAbs F÷180▲Q“HZ”= Z[1]＋S▲T=K－M▲L=Q－M▲E=（R＋S2÷（24R））÷cos（F÷2）－R▲Lb1 0：{G}：G“ZHUANG HAO”：G＜O＝=＞Z=G－M≠=＞Z=Q－G：△{U}：U“I——B OUT＋IN－”Z＜S＋0.005＝=＞H=tan-1（24Z2S2R2－Z^6）÷（48Z^3R^3－6Z^4 SR））：V=Z－Z^5÷（40S2R2）＋Usin H：W=Z^3÷（6SR）－Z ^ 7÷（336Z^3R^3）－Ucos H：≠=＞H=180（Z－S÷2）÷（πR）：V=Rsin H＋S÷2－S^3÷（240R2）＋Usin H：W=R－Rcos H＋S2÷（24R）－Ucos H：△F＜0＝=＞G＜O＝=＞X=V：Y=W：≠=＞X=T＋Tcos F－Vcos F－Ws in Abs F：Y=Tsin Abs F－Vsin Abs F＋Wcos F：△≠=＞G＜O＝=＞X=V：Y=－W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Vsi n Abs F－Tsin Abs F－Wcos F：△△C=A＝=＞D≥B＝=＞P=90：≠=＞P=270：△≠=＞C＜A＝=＞P= tan-1（（D－B）÷（C－A））＋180：≠=＞P= tan-1（（D－B）÷（C－A））：△△I“COM——N”=A＋√（X2＋Y2）cos（P－tan-1（Y÷X））▲J“COM——E”=B＋√（X2＋Y2）sin（P－tan-1（Y÷X））▲G=G＋20Goto 0程序运行输入：JD？——交点桩号R？——圆曲线半径LS？——缓和曲线长度PJ——L－R＋？——偏角右偏为正ZH——N？——ZH点N坐标ZH——E？——ZH点E坐标JD——N？——JD点N坐标JD——E？——JD点E坐标输出：ZH= ——ZH点桩号HY= ——HY点桩号QZ= ——QZ点桩号YH= ——YH点桩号HZ= ——HZ点桩号T= ——切线长L= ——曲线长E= ——外矢距再输入：ZHUANG HAO？——输入任一点桩号I——B OUT＋IN－？——距中线距离（外侧为正）再输出：COM——N= ——计算出N坐标COM——E= ——计算出E坐标以此循环——再输入——再输出程序中各字母所代表的含义：A——ZH点N坐标B——ZH点E坐标C——JD点N坐标D——JD点E坐标E——外矢距F——偏角（右偏为正）G——弯道内任一点桩号H——弯道内任一点切线角I——计算出的N坐标J——计算出的E坐标K——交点桩号L——曲线长M——ZH点桩号N——HY点桩号O——QZ点桩号P——ZH--JD方位角Q——HZ点桩号R——圆曲线半径S——缓和曲线长度T——切线长U——距中线距离（外侧为正）V——支距W——支距X——弯道内坐标Y——弯道内坐标Z——任一点曲线长度Z[1]——YH点桩号。
公路路线坐标计算说明书_1.3.0
智璟安卓工程计算系统公路路线坐标计算说明书一、计算原理1、公路路线坐标计算是用线元法方式进行坐标、方位角（均指中桩切线方位角）的计算的，其中缓和曲线是用复化辛普森公式进行计算的。
线元法计算方式通用性强，不用区分直线、圆曲线、缓和曲线的组合方式，计算对象为每一个独立的线元，而非组合线型。
计算方式为通过逐个计算每个线元的终点坐标、方位角来计算整条路线上任意一点的坐标和方位角。
2、公路路线坐标计算所指的线元仅为直线、圆曲线和缓和曲线（回旋线）这三种，要求两线元交接点的曲率半径和方位角相同，即前一线元终点的曲率半径、方位角等于后一线元起点的曲率半径、方位角；如果两线元交接点的曲率半径不同，需要采用加入虚线元的方式进行处理，但仍需保证两线元交接点的方位角相同。
3、计算前先要建立整条路线的基本数据及线元终点桩号、终点曲率半径的数据库，计算时选择路线文件后只要输入所在路线上的任何一个桩号，即可进行此桩号坐标和方位角的计算。
4、建立路线数据库可以在程序界面提示下逐个手工输入数据，也可以通过外部TXT 文本文件导入数据，毕竟用电脑在TXT文件中输入数据更快、更方便。
5、由于计算任意桩号的坐标、方位角都是从整条路线的起点坐标、方位角开始逐个线元计算出来的，因此为了减少累积误差，均采用较高的计算精度，在利用辛普森公式迭代计算时，曲线分段数（迭代次数）均给予了偏高的取值，保证每次计算X、Y坐标同时精确到0.000001；同时角度（8位小数）、半径（6位小数）数据均要求较高的输入精度。
由于计算方式、计算精度、数据输入精度（角度、半径）的不同，计算结果可能与设计给出的结果存在微小的偏差，属正常现象。
6、程序未提供断链处理功能，如果整条路线有断链情况，请从断链处分段建立不同的路线文件。
7、程序只适合于一般手机竖屏显示状态下的正常使用。
二、功能介绍1、对整条路线上任意桩号的中桩、正交边桩、斜交边桩的坐标进行计算，包括中桩切线方位角。
CAISO5800公路坐标、高程计算程序
“E=”(R+(Z[6]+Z[7])/2)/cos(A/2)-R→E◢
“K(ZH)=”:K-S→Z[1]
“K(HY)=”:Z[1]+B→Z[2]◢
“K(QZ)=”:Z[1]+L/2+(B-C)/4→Z[3]◢
“K(YH)=”:Z[1]+L-C→Z[4]◢
“K(HZ)=”:Z[4]+C→Z[5]◢
Deg:Fix 3↓
“X0”?X:”Y0”?Y:”R”?R”FWJ”?:C:”KS”?S↓
Lbl 0↓
“KP”?G:G-S→L↓
90L/R/π→M↓
Abs(2RSin(M))→O↓
“XP=”:X+Ocos(C+M)→I◢
“YP=”:Y+Osin(C+M)→J◢
C+2M→F↓
“BP=”:F▲DMS◢
“RE”?T:If T≠0:Then 1/T→A:Else 0→A:IfEnd↓
“X0”?U:”Y0”?V↓
Lbl 0↓
“KP”?P↓
P-S→L↓
(B-A)/2/C→M↓
πQ/180→K↓
“XP=”:W+∫(cos(K+AX+MX^2),0,L)→F◢
“YP=”:Y+∫(sin(K+AX+MX^2),0,L)→G◢
“D1=”I◢
Lbl 2↓
“D ANGEL”?H:If H=0:Then Goto 0:IfEnd:”D”?D↓
“XB=”:F+Dcos(Z+H)→D◢
“YB=”:G+Dsin(Z+H)→H◢
高速公路坐标计算方法
高速公路坐标高程计算程序本软件简要说明：一、平曲线计算（主程序）1、J为起算点里程，C、D为起算点的X、Y坐标，F为起算点的切线方位角，R为圆曲线半径（左偏取负，右偏取正），A、B为第一、第二缓和曲线回旋参数，O为圆曲线长度，Ki为该分段的终点里程；2、对于直线段或圆曲线段，起算点可取直线或圆曲线上的任意一点；3、对于带第一、第二缓和曲线的平曲线段，起算点应取HY点；4、K为所求点的里程，T、P为第一偏距、偏角，S、Z为第二偏距、偏角，偏角取从该点的切线顺时针旋转的夹角；5、分段法则：直线单独分段；单一的圆曲线单独分段；缓和曲线1+圆曲线+缓和曲线2为一个整体单独分段，若不存在第一或第二缓和曲线（即不完全缓和曲线）仍然可以计算，A或B可取任意不为零的值；若不存在圆曲线，则O取零；6、无论任何时候A、B不能取零，否则可能导致被零除的错误；7、F、Q切线方位角输入输出均为度.分秒的格式，例如153°24′05.24″=153.240524。
Q改变时，可按照新方位角为基准，结合第一第二偏距、偏角重新计算所求点；8、输入平曲线参数后，默认为计算全线坐标，可修改来计算某段曲线，默认间距也可修改；9、可参考CAD图《平曲线计算图例》；10、生成的中桩CAD脚本设置成在世界坐标系下生成，注意的是世界坐标系与大地测量坐标系的区别是XY坐标是互换的，否则画出的图形与实际相反。
先打开CAD，设置好图层名称、颜色，并设置为当前层，然后单击CAD的工具==>运行脚本==>选中生成的脚本文件即可。
11、输出的坐标结果可以导入到EXCEL中，操作办法为：打开EXCEL，然后把坐标数据复制到单元格里，然后单击数据==>分列==>选中分隔符号==>下一步==>选中TAB键和逗号==>下一步==>完成即可。
下一次可直接在此表中粘贴，数据自动分列。
二、缓和曲线计算（辅助程序）1、本程序为辅助程序，用来从ZH点或HZ点计算整条完全的缓和曲线，若不知道HY点X、Y、Q参数，可用此程序计算出来，然后输入平曲线参数；2、参数设置参考平曲线计算；3、导出到EXCEL的办法同平曲线计算；三、直线计算（辅助程序）1、本程序为辅助程序，若已知P1（X1，Y1），P1-->P2的距离I及方位角J（度.分秒格式），可计算坐标P2（X2，Y2）。
公路中线上点的坐标计算
1. 公路中线上点的坐标计算当需放样的点位于公路中线上时，如图1，各JDi的坐标(Xi，Yi)在控制测量阶段就已经测定(或由施工图文件中《直线、曲线及转角表》中查出)，相邻JD连线的坐标方位角Ai -1，i可由同样方法查出，或利用JD坐标反算推出。
各曲线主点坐标可由《直线、曲线及转角表》查出，或由曲线要素值及，计算得到。
1.1直线上各中桩坐标计算当需要放样的P点位于直线上时，有两种情况：位于YZ（HZ）之间和ZY（ZH）之间，或者位于公路QD和ZH（ZY）之间，其计算方法相同，公式如下：式中为该段直线的起点（可以是YZ，HZ，或QD）坐标为要求的P点与该段直线起点的桩号差（距离）1.2 单圆曲线上各中桩坐标计算当需要放样P点位于单圆曲线上时，其坐标计算如下：式中为ZY点坐标，为圆曲线半径为P点与ZY 点的桩号差（弧长）当路线左转时，取“-”，反之取“+”1.3 带缓和曲线的圆曲线上各中桩坐标计算当P点位于带缓和曲线的圆曲线时，又分为以下三种情况：ZH到HY段式中为ZH点坐标为P点与ZH点桩号差，为缓和曲线长当路线左转时，取“-”，反之取“+”HY到YH段式中为HY点坐标为P点与ZH点桩号差，为缓和曲线长当路线左转时，取“-”，反之取“+”YH到HZ段式中，为HZ点坐标，为HZ点与P点的桩号差当路线左转时，取“+”，反之取“-”1.4 复曲线上各中桩坐标计算1.4.1 当复曲线中间不设缓和曲线时，采用以下方法进行计算：对于第一缓和曲线、第一段圆曲线以及第二缓和曲线，分别用公式（3）、公式（4）和公式（5）计算；对于第二段圆曲线，用公式（2）计算，计算时将公式（2）中的换成，分别为第一圆曲线和第一缓和曲线长度，左转取“-”，右转取“+”。
1.4.2 当复曲线中间有缓和曲线时，即构成卵型曲线。
如图2，缓和曲线AB的长度为，A、B点的曲率半径分别为和，为缓和曲线上曲率为零的点，AB段内任意点的坐标从点推算。
公路工程卡西欧5800坐标及高程计算程序(带超高计算程序)
本套程序包含的主程序及子程序：
1、主程序：0---ZXJS
2、坐标计算判断数据库：ZXDATA
3、坐标正算子程序：SUB1
4、坐标反算子程序：SUB2
5、竖曲线判断数据库程序：ZXGC
6、超高段判断数据库程序：ZXGC
7、竖曲线高程计算子程序：SUB3
8、竖曲线超高计算子程序：SUB4
H-Z((J-I)(3((S-O)÷A)^2-2((S-O)÷A)^3)+I)→H
Return
If Abs(Z)<1e^(-4):Then Goto 1:Else W+Z→W
Goto 0
Lbl 1
0→Z:Prog “SUB1”:(J-Y)÷sin(F)→Z
5：竖曲线判断数据库程序：ZXGC
参数含义：
B：变坡点桩号
G：变坡点高程
I：变坡点前坡度(区分正负，带符号)
J：变坡点后坡度(区分正负，带符号)
V+W(Z[1]sin(G+QKZ[3]W(C+Z[3]WD))+Z[2]sin(G+QKZ[4]W(C+Z[4]WD))+Z[2]sin(G+QKFW(C+FWD))+Z[1]sin(G+QKMW(C+MWD))→Y
G+QKW(C+WD)→Z[5]
Z[5]+E→F
X+Zcos(F)→X
Y+Zsin(F)→Y
Goto B
2、坐标计算判断数据库：ZXDATA
遵循 直线-缓和曲线-圆曲线-缓和曲线-直线 法则 无缓和曲线可以省略缓和曲线参数
公路中线坐标计算程序
公路中线坐标计算程序[CASIO4500]公路中线坐标计算程序 [CASIO 4500][转帖]公路中线坐标计算程序 [casio 4500]ZXCLCXDefm16Lbl1X=1:{X}:X“ZJF-1 ZBF-2 TT-3 SJD-4 P2L-5 B2A-6”X=1＝>Prog"ZJF"△X=2＝>Prog“ZBF”△X=3＝>Prog“TT"△X=4＝>Prog“SJD"△X=5＝>Prog“P2L"△X=6＝>Prog“B2A"△Goto1ZBF“-----ZXZBJS-----”:{YD}:Y“I1”“I2”Y≠D=>Z=0rog“FP”:A=B:B=0:Goto3△X=0:{TPXVWU}:T“HJD(X+Yi)”Lbl1X“XJA=1”“JD(X+Yi)”:U“JDLC”:Y=Arg(P-T)Y<0=>Y=Y+360△Y“I”◢C“L”=Abs(P-T)◢X=0=>V=P:Goto2△V“JDB(X+Yi)”=Arg(V-P)D<0=>D=D+360△D“I”◢Z“L”=Abs(V-P)◢Prog"FP":A=BLbl2W“QJD(X+Yi)”=Arg(W-V)D<0=>D=D+360△D“I”◢C“L”=Abs(W-V)◢Lbl4Prog“FP”X=1=>B=B-A△X=0=>A=B:B=0:Z=0:E=Rep P:F=ImpP:Goto3△M=Rep(P-T):N=Imp(P-T):I=Rep(W-V):J=Imp(W-V):F“Y”=(M Imp T÷N-I Imp V÷J+Rep V-Rep T)÷(M÷N-I÷J)◢E“X”=M(F-Imp T)÷N+Rep T◢C“L1”=√((E-Rep P)2+(F-Imp P)2)◢U=U+C:C“L2”=√((E-Rep V)2+(F-Imp V)2)◢Lbl3X=0:{X}:X“ZBF(0)ZJF(1)”X=1=>O=Zrog"ZJF"△A=A+Brog"ZBJS"ZBJSI=0:J=0:B=0:{RGH}:R:G“LS1”:H“LS2”rog“TWE”Z=0=>E“JD(X)”:F“JD(Y)”:U“JDLC”△Rec(-Z[7],Y):Z[15]=I+E:Z[16]=J+F:Rec(Z[8],D):Z[13]=I+E:Z[14]=J+F:C=U-Z[7]:Z[6]“ZH”=C◢Prog“ZDZ”C=Z[6]Lbl1{C}:C“JIAZHUANG”:Z=C-Z[6]Z≤0=>Rec(Z,Y):Z[11]“X”=I+Z[15]◢Z[12]“Y”=J+Z[16]◢Prog“BZ”:Goto1△Z≤G=>S=G=Z[1]:Q=Z[3]rog“H”ol(Z[11],Z[12]):Rec(I,Y-JO):Z[11]“X”=I+Z[15]◢Z[12]“Y”=J+Z[16]◢Prog“BZ”:Goto1△Z≤W-H=>S=G=Z[1]:Q=Z[3]rog“Y”ol(Z[11],Z[12]):Rec(I,Y-JO):Z[11]“X”=I+Z[15]◢Z[12]“Y”=J+Z[16]◢Prog"BZ":Goto1△Z≤W=>Z=W-Z=H=Z[2]:Q=Z[4]rog"H"ol(Z[11],Z[12]):Rec(I,D+180+JO):Z[11]“X"=I+Z[13]◢Z[12]“Y"=J+Z[14]◢Z=W-Zrog“BZ”:Goto1△Z＞W=>Z=Z-W:Rec(Z,D):Z[11]“X"=I+Z[13]◢Z[12]“Y"=J+Z[14]◢Z=W+Zrog “BZ”:Goto1△FPB=Y-DB>180=>B=B-360△B<-180=>B=B+360△O=B÷Abs B:B=Abs BTWES=Grog“PQ”:Z[1]=P:Z[3]=Q=Hrog“PQ”:Z[2]=P:Z[4]=QZ[7]“T1”=(R+Z[2]-(R+Z[1])cos(A+B))÷sin(A+B)+Z[3]◢Z[8]“T2”=(R+Z[1]-(R+Z[2])cos(A+B))÷sin(A+B)+Z[4]◢W“L”=πR(A+B)÷180+.5G+.5H◢X“E”=(R+Z[1])÷cos(tan-1(((R+Z[2])÷(R+Z[1])-cos(A+B))÷sin(A+B)))-R◢Z[9]=Z[7]-Osin B÷sin(A+B):Z[10]=Z[8]-Osin A÷sin(A+B)B≠0=>Z[9]“TA”◢Z[10]“TB”◢△。
卡西欧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 取正值，反之取负值。