公路坐标计算程序1.0

公路测量坐标计算公式引言公路测量是一项基础工作，用于确定公路建设或维护所需的各个节点坐标。

在公路工程领域，测量坐标计算公式是至关重要的工具，用于测算和确定公路各个位置的坐标信息。

本文将介绍公路测量坐标计算公式的原理和应用。

坐标系统在公路测量中，使用的坐标系统通常是平面直角坐标系。

这个坐标系统由X轴和Y轴组成，其中X轴表示东西方向，Y轴表示南北方向。

公路测量中，测量标准一般会规定一个起始点作为基准点，所有的测量点都以该基准点为原点建立坐标。

公路测量坐标计算公式坐标计算原理坐标计算公式的原理是通过已知的测量数据和几何关系，推导出待求点的坐标。

常用的坐标计算方法包括三角测量法、坐标转换法和横断面测量法等。

三角测量法三角测量法是基于三角形几何学原理的坐标计算方法。

它通过测量已知点与待求点之间的夹角和距离，利用三角函数关系计算出待求点的坐标。

三角测量法适用于平面内的测量，并具有较高的精度。

坐标转换法坐标转换法是将已知点的坐标转换到待求点坐标的计算方法。

它利用已知点和待求点在同一坐标系统中的相对位置关系，通过坐标转换公式计算待求点的坐标。

坐标转换法适用于已知点坐标较为准确的情况。

横断面测量法横断面测量法适用于公路等线性工程的坐标计算。

它通过测量已知点的高程和于待求点的高程差，利用高程差和水平距离之间的关系计算待求点的坐标。

横断面测量法适用于公路中断面的测量和计算。

应用示例公路测量坐标计算公式在实际工程中具有广泛的应用。

下面以一个简单的示例来说明其应用过程：假设有一段公路，已知起点的坐标为（0,0），终点的坐标为（1000,0）。

现在需要测算该公路上距离起点100米处的坐标。

根据三角测量法，可以通过测量起点和待求点之间的夹角和距离来计算待求点的坐标。

假设测量得到的夹角为45度，距离为100米。

根据三角函数的性质，可以计算出待求点的坐标为（100,100）。

总结公路测量坐标计算公式是公路工程中不可或缺的工具。

通过三角测量法、坐标转换法和横断面测量法等方法，可以准确计算公路上各个位置的坐标信息。

高速公路测量计算CASIO程序全套坐标主程序XY(坐标程序)LbI 2｛K｝PROG “SUB3”X=X◢Y=Y◢｛v｝v=√2 PROG “SUB 2”：⊿GOTO 2FANS｛MNWT｝M“X0=” ：N“Y0=”W“X-=” ：T“Y-=”Z=tan-1((T-N)/(W-M+0.19)LbI 7｛UV｝U“JIAO=”V“JU=”W＞M U=U+Z：GOTO 3：⊿U=U+Z-180LbI 3P=M+VcosU◢Q=N+VsinU◢｛K｝K“DGLC=”LbI 6Prog “SUB 3”O=（Q-Y+（P-X）/（tanE+0.19））×sinEK=K+OO=Abs OO≥0.17GOTO 6：⊿K=K◢V=（Q-Y+（X-P）×tanE）×cosE ◢GOTO 7子程序ZHIX 1H= K-80100 （80+100为起点里程，直线段）A=3026441.522 （A、B为K80+100坐标）B=515965.236E=270.265 为K80+100方位角X=A+HcosE Y=B+HsinEQUX 1 （K80+100~ K80+948.114曲线右偏）R=900：L=170H=K-80310.927 （80+310.927为ZH点里程）I=K-80480.927 （80+480.927为HY点里程）I≤0A=3026442.499B=515754.311 （A、B为ZH点坐标）PROG “SUB 1”E=270.265+3×G （270.265为ZH点方位角）X=A+Fcos（E-2G）Y=B+Fsin（E-2G）Goto 2 ⊿J=K-80778.114J≤0A=3026448.633 ：B=515584.489 （A、B为HY点坐标）C=I×90÷π÷RD=2×R×sinCE=275.6763+2C （275.6763为HY点方位角）X=A+D×cos（E-C）Y=B+D×sin（E-C）Goto 2 ⊿S=K-80948.114S≤0A=3026606.193 ：B=515149.157 （A、B为HZ点坐标）H=-SPROG “SUB 1”E=300.0072-3×G （300.0072为HZ点方位角）X=A+Fcos（E+180+2G）Y=B+Fsin（E+180+2G）Goto 2 ⊿LbI 2备注：曲线左偏与右偏，只有中间几个“＋－”号的区别。

高速公路坐标高程计算程序本软件简要说明：一、平曲线计算（主程序）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）。

坐标计算程序及坐标计算公式一、坐标计算公式直线段：中桩公式：x=x1+(z-c)*cosay=y1+(z-c)*sina边桩公式：x=x1+(z-c)*cosa±d*cos(a-90)y=y1+(z-c)*sina±d*sin(a-90)说明： x1——起点x坐标，y1——起点y坐标，z——计算点桩号， c——起点桩号，a——方位角，d——距中桩距离。

“±”左边桩为“+”，右边桩为“-”。

二、方位角计算公式1、直线段方位角图纸提供。

2、若为单一圆曲线时，起点用直线段方位角图纸提供。

3、若为缓和曲线时：第一段缓和曲线方位角为直线段方位角图纸提供；第二段圆曲线起点方位角a1=a0±Ls*180/2∏r。

第三段缓和曲线方位角a2=a1±L*360/2∏r。

（a0为直线段方位角，Ls为缓和曲线长度，L为圆曲线长度，r为圆曲线半径，“±”右偏角为“+”、左偏角为“-”。

）三、5800计算器坐标程序坐标程序由1个主程序、5个子程序和1个数据库组成进入计算器编程模式（5：PROG—1：NEW新建程序名—3：EDIT），输写程序。

1、主程序adminFix 3（回车换行）Lb1 0:150→DimZ:“ZX→0,A→1,B1→2,B2→3,C→4,D→5,CR→6”?N:N→Z[149]:Prog “DAT2”：“ZS=>1，FS=>2”?N:If N=1:Then Goto 1:Else N=2=>Goto 2:IfEnd:Goto 0:Lb1 1:?S:“ANG=”?M:?Z:S=0=>Goto 0:0→N（回车换行）Lb1 5:Isz N:If S≤Z[8N+2]+Z[8N+4]:Then N→J:Prog“DAT1”:Else Goto5:IfEnd（回车换行）Abs(S-0)→W:Prog“SUB1”：“XS=”：X◢“YS=” ：Y◢“FWI=”：F-M→F：If F≤360:Then F→F:Else F-360→F:IfEnd:F►DMS◢Goto4（回车换行）Lb1 2:?X:?Y:“ANG=”？M：M→Z[148]：If M＜90:Then 180-M→M:IfEnd（回车换行）X-Z[4]:Y-Z[5]:X=0=>Goto 0:0→N（回车换行）Lb1 A:Isz N:N→Z[150]:Z[8N+3]-M→A:Z[8(N+1)+3]-M→B:Prog “ZX1”:If Z[6]×Z[7]≤0:Then N→J:Prog “DAT1”:Goto B:IfEnd（回车换行）Z[8N+3]+M→A:Z[8(N+1)+3]+M→B:Prog “ZX1”:If Z[6]×Z[7] ≤0:Then N→J:Prog “DAT1”:Goto B:Else Goto A: IfEnd（回车换行）Lb1 B:Prog “SUB2” （回车换行）Z[150]→N:0+W→S:If S>Z[8N+2]+Z[8N+4]+.001:Then Goto A: IfEnd（回车换行）If N>13:Then 0→N: Goto A: IfEnd（回车换行）If Z[148]>90: Then S+2Zsin(M-90) →S: IfEnd:“S=”:S◢“Z=”: Z◢“OK→2,NO→1”?N:If N=1: Then Z[150]→N: Goto A:Else Goto 2: IfEnd（回车换行）Lb1 4 :0→J:“DIST=”：Pol （X－Z[1]，Y－Z[2]）→I◢ J-F：If F<0：Then F+360→F：IfEnd：“FW=”： F►DMS◢ Goto 12、子程序DAT1Z[8J]→U：Z[8J+1]→V：Z[8J+2]→O：Z[8J+3]→G：Z[8J+4]→G：Z[8J+4]→H：Z[8J+5]→P：Z[8J+6]→R：Z[8J+7]→Q：1÷P→C：（P－R）÷（2HPR）→D：（180÷∏）→E：Return3、子程序ss“A”？→A：If A<0:Then 10A◢Else 9A◢IfEnd:Ans×1.05（回车）4、子程序SUB1.1184634425→A：.2393143352→B：.28444444444→N：.046910077→K：.2307653449→L：.5→Z[3]:U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Ncos(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)))→X:V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Nsin(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))) →Y：G+QEW(C+WD)+M→F：X+Zcos（F）→X：Y+Zsin（F）→Y: Return5、子程序SUB2G-M→T：Abs（（Y-V）cos（T）-（X-U）sin（T））→W：0→Z：Lb1 0：Pros “SUB1” （回车换行）T+QEW（C+WD）→L：（Z[5]-Y）cos（L）-（Z[4]-X）sin（L）→Z：If Abs（Z）<1×10∧(-4):Then Goto 1:Else W+Z→W: Goto 0: IfEnd（回车换行）Lb1 1:0→Z:Prog“SUB1”：（Z[5]-Y）÷sin（F-2M+180）→Z：Return6、子程序ZX1（Z[5]-Z[8N+1]）cos（A）-（Z[4]-z[8N]）sin（A）→Z[6]：（Z[5]-Z[8（N+1）+1]）C5]-Z[8（N+1）+1]）cos（B）-（Z[4]-Z[8（N+1）]sin（B）→Z[7]：Return7、数据库DAT2If N=0:Then起点X坐标→Z[8]：起点Y坐标→Z[9]：起点桩号→Z[10]：起点坐标方位角→Z[11]：曲线长度→Z[12]：起点半径→Z[13]：终点半径→Z[14]：曲线转向（左转为“-1”右转为“+1”直线为“0”）→Z[15]：Return IfEnd（依次把所有平曲线要素输完）If N=1:Then起点X坐标→Z[8]：起点Y坐标→Z[9]：起点桩号→Z[10]：起点坐标方位角→Z[11]：曲线长度→Z[12]：起点半径→Z[13]：终点半径→Z[14]：曲线转向（左转为“-1”右转为“+1”直线为“0”）→Z[15]：Return IfEnd（依次把所有平曲线要素输完）说明：1、所有路线都是从Z[8]开始。

5800计算器公路全线坐标正、反算计算程序FX5800全线贯通万能正、反算程序（一体化、超好用、短小、易懂）FX5800计算器的积分程序（正反算、全线贯通、新线路）终极版ZHUCHENGXU主程序"1.ZS,2.FS" ?→Q←┘输入1正算，输入2反算“NEW=0,OLD≠0”?Z←┘IfZ=0:Then“X0=”?A:“Y0=”?B:“C0=”?C:“1/R0=”?D:“1/RI=”?E:“SP=”?F:“EP=”?G:Ifend:Q=2=>Goto 2←┘Lbl1 :“KM=,<0 Stop”?H:H<0=>Stop:“PJ=”?O:“PY=”?L←┘LblZ:Z=1=> Prog“01”：Z=2=> Prog“02”←┘选择数据库文件，可增加H- F→X:0.5(E-D)÷(G-F)→N←┘C+(XD+NX2)*180÷π→P：P<0=>P+360→P：P>360=>P-360→P←┘-A+∫(cos(C+(XD+NX2)*180÷π),0,X)+Lcos(P+O)→U←┘B+∫(sin(C+(XD+NX2)*180÷π),0,X)+Lsin(P+O)→V←┘Q=2=>Goto4：Cls：Fix 3←┘"Xn=":Locate4，1，U："Yn=": Locate5，2，V：“FWJ=”：PDMS◢Norm 2：Cls：Goto1←┘Lbl2:“XD=,<0,STOP”？R：R<0=>Stop:“YD=”？S←┘“KMDG=”？H ：90→O:0→L：GotoZ←┘（H线路范围内的任意桩号）Lbl4:Pol（R-U，S-V）：J<0 => J+360→J←┘Whileabs（Icos（J-P））≤0.001:P-J>180=> J+360→J： P-J<-180=> P+360→P:If P-J>0:then -I→L:else I→L :Ifend:Goto3: Whileend：H+Icos（J-P）→H:GotoZ←┘Lbl3:Cls：Fix 3←┘“KM=”: Locate4，1，H：“PY=”: Locate4，2，L◢Norm 2：Cls：Goto2←┘01（数据库子程序）If H<=第一曲线终点桩号：then第一曲线起点X→A：第一曲线起点Y→B：第一曲线起点方位角→C：起点曲率→D：终点曲率→E：起点桩号→F：终点桩号→G：return：ifend ……………程序说明：1、该程序可以计算任意线形（直线、圆曲线、缓和曲线、不完整曲线）任意桩号的坐标（正算，输入1），也可根据坐标计算该点到线路的距离及垂足桩号（反算，输入2）；2、（NEW=0，OLD≠0）？如果要计算的点为数据库线路中的点，则输入数据库编号（以整数1、2、3…代替输入）；如果在数据库中没有要计算线路的数据，则输入曲线要素X0：曲线起点X坐标；Y0：曲线起点Y坐标；C0:曲线起点方位角；R0-1、 RI-1:曲线起点、终点曲率，直线为0，曲线左偏输入负值，右偏输入正值；SP、 EP：曲线起点桩号，终点桩号；KM：待求点桩号；PJ：正斜交的设定；PY：偏中距离，线路上的点输入0，右偏输入+值，左偏输入-值；3、正算显示坐标及切线方位角；反算输入线路的任意桩号（此桩号越接近真实值计算速度越快）、待求点坐标，显示待求点桩号及偏中距离；4、正算子程序为积分公式编写而成；反算子程序为角度趋近的方法编写，计算速度有点慢。

ZBZS★坐标正算主程序Deg:FixmLbl 0:｛O｝:O“1,47264-50167.162,50150-50673.503,50700-54700”Lbl 1:｛ZKB｝:Z:K:B:O=1=＞Goto 2△O=2=＞Goto 3△O=3=＞Goto 4△Lbl 2:Z＜47264=＞Goto 5△Z＞50167.16=＞Goto 5△Z＜48708.52=＞Prog“55”: Goto6△Z＜49310.50=＞Prog“56”: Goto6△Z＜49874.21=＞Prog“57”: Goto6△Prog“57-1”: Goto6△Lbl 3:Z＜50150=＞Goto 5△Z＞50673.5=＞Goto 5△Prog“58”: Goto6△Lbl 4:Z＜50700=＞Goto 5△Z＞54700=＞Goto 5△Z＜51265.12=＞Prog“59”: Goto6△Z＜52266.78=＞Prog“60”: Goto6△Z＜53341.82=＞Prog“61”: Goto6△Z＜53868.03=＞Prog“62”: Goto6△Prog“63”:Goto 6△Lbl 5:“ZH ERROR”▲Goto 0Lbl 6:Prog“DQXY”:“X=”:X:Pause 0:“Y=”:Y▲Pol（X-Q,Y-U）:J＜0=＞J=360+J△:I=Intg（1000I+.5）÷1000“L=”:I:Pause 0:“T=”:J→DMS▲Goto 1 说明：本程序适用于casio4850计算器。

把线路从断链处分为三段，O的取值，待求点在第一段输1，第二段输2，第三段输3。

Z、待求点的桩号。

K、待求点在左幅输0右幅输1。

B、待求点距中距离。

Goto 5前的桩号为起止点或断链点桩号。

Prog“xx”: Goto 6前的桩号为曲线组合的HZ点桩号。

Q、U为测站坐标，运行程序前先存入。

公路逐桩坐标计算程序(可以计算对称、不对称缓和曲线)Lb1 0Z=？V=？W=V+2：Fixm｛K｝Lb1 1K>Z[W+5Z+4]=>W=W+1:Goto 1⊿(判断桩号在哪个交点范围，就是该交点曲线起点至下一交点曲线起点) S=K-Z[W+5Z+3] （计算该桩号与曲线起点的距离）R=Z[W+2Z+2]：L=Z[W+3Z+2]:E=Z[W+4Z+2] （读取该交点曲线要素R、Ls1 、Ls2）Pol（Z[W]-Z[W-1],Z[W+Z+2]-Z[W+Z+1]）（计算该交点与下一交点直线方位角）J<0=>J=J+360⊿A=JPol（Z[W-1]-Z[W-2],Z[W+Z+1]-Z[W+Z]）（计算该交点与上一交点直线方位角）J<0=>J=J+360⊿C=A-J:A=J （计算偏角）W=V+2=>Goto2⊿（如果桩号在起点与第一交点曲线起点之间，则转Lb1 2 ）I=Abs（tan（c÷2））M=L÷2-L^3÷240R^2:N=E÷2-E^3÷240R^2P=L^2÷6R-L^4÷336R^3-R（1-cos（90L÷πR））Q=E^2÷6R-E^4÷336R^3-R（1-cos（90E÷πR））D=（P-Q）I÷2 : F=（P+Q+2R）I÷2M=F+M-D:Q=F+N+DN=πRAbsC÷180+（L+E）÷2X=Z[W-1]-McosAY=Z[W+Z+1]-MsinAM=Z[W-1]+Qcos（A+C）V=Z[W+Z+1]+Qsin（A+C）Q=AbsC÷CS≤L=>P=0:Goto3⊿（如果桩号在第一缓和曲线内，则转Lb1 3）S≤N-E=>S=S-L:Goto4⊿（如果桩号在圆曲线内，则转Lb1 4）S≤N=>S=N-SQ=-Q:A=A+C-180:X=M:Y=V:L=E：P=180:Goto3 ⊿（如果桩号在第二缓和曲线内，则转Lb1 3）P=A+C:S=S-N:D=M+ScosP:F=V+SsinPGoto6 （如果桩号在直线内，则转Lb1 6）Lb1 2P=A+CD=Z[W-1]+ScosPF=Z[W+Z+1]+SsinP:Goto6Lb1 3I=S-S^5÷40R^2÷L^2+S^9÷3456R^4÷L^4J=Q（S^3÷6RL-S^7÷336R^3÷L^3）P=P+A+90Q S^2÷πRL:Goto5Lb1 4M=90（2S+L）÷πRI=RsinM+L÷2-L^3÷240R^2J=Q（L^2÷24R+R（1-cosM））P=A+QMLb1 5D=X+IcosA-jsinA:F=Y+JcosA+IsinALb1 6D″X=″◢（结果显示X坐标）F″Y=″◢（结果显示Y坐标）P″AT=″◢（结果显示该桩号方位角）{BO}：B″S″O″⊿″ （输入边桩距离，交角）P=P+OL″XB″=D+BcosP◢（结果显示边桩X坐标）M″YB″=F+Bs inP◢（结果显示边桩Y坐标）以上是坐标计算程序，括号内是程序计算的大致原理及说明，中间部分为直线、圆曲线、缓和曲线计算的各种公式，大家也知道，书上也有。