道路坐标计算程序

坐标计算程序及坐标计算公式一、坐标计算公式直线段：中桩公式：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]开始。

公路道路上有关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分别为第一圆曲线和第一缓和曲线长度，左转取“—”，右转取“+”。

道路曲线坐标计算步骤已知量：转角α=23°3’38”缓和曲线长L 0=100m圆曲线半径R=1000m圆曲线长L=00180)2(πβα⨯-⨯R 曲线长 L h ==000180)2(2πβα⨯-⨯+R L起始边方位角A zh-jd =tan -1(Y jd -Y zh /X jd -X zh )= tan -1((750-500)/ (750-500))=45°切线加长 q=23002402R L L -圆曲线相对切线内移量 p=L 02/(24R)切线长T h = q ＋(R ＋p)•tan(α/2)第一步计算ZH 坐标：方位角A jd-zh =A zh-jd +180°=225°T= q ＋(R ＋p)•tan(α/2) ___这里α为转角Xzh=Xjd+TcosA jd-zhYzh=Yjd+TsinA jd-zh第二步计算HY 坐标：缓和曲线切线角 βi =π001802⨯R L缓和曲线偏角(i)： δi=βi /3=L i 2/6RL s *180°/π缓和曲线方位角：аi 缓=A jd-zh+δi====== A jd-zh 已知（45°）用方位角(线路向左转A jd-zh -δi ) 缓和曲线坐标(i)： X i =230040R L L -Y I = L 02/6R缓和曲线ZH 的弦长坐标(i)：C i =22Y X +所以缓和曲线最后要求坐标为： X HY =X zh + C i cos аi 缓Y HY =Y zh + C i sin аi 缓 第三步计算QZ 的坐标：外矢距 E h = (R ＋p)/cos(α/2)－R 这里α为转角αjd-qz = A jd-zh +90°+α/2 这里α为转角线路向左转（A jd-zh -90°-α/2） QZ 的坐标X qz =X jd +Ecos αjd-qz Y qz =Y jd +Esin αjd-qz第四步DK1+100坐标计算：ZH 点里程=JD 里程-T=1300-254.08=1045.92=DK1+45.92L i = DK1+100-ZH 点里程=1100-1045.92=54.08缓和曲线切线角 βi =π0021802⨯RL L i缓和曲线偏角(i)： δi=βi /3缓和曲线方位角：аi =A jd-zh+δi====== A jd-zh 已知（45°）用方位角(线路向左转A jd-zh -δi ) 缓和曲线坐标(i)： X i =02540L R L L i i -Y i =L i 3/6RL 0缓和曲线ZH 的弦长坐标(i)：C i =22i i Y X +所以缓和曲线最后要求坐标为： X HY =X zh + C i cos аiY HY =Y zh + C i sin аi 第五步DK1+280坐标计算：HY 点到DK1+280的弧长L i =DK1+280- DK1+45.92-L 0=134.08 HY 点到DK1+280的弧长对应的圆心角Ψi =π0180⨯R L i 圆曲线HY 点到任意一点DK1+280的偏角△i =Ψi /2 圆曲线HY 点到任意一点DK1+280的弦长C i =2Rsin(△i ) 圆曲线HY 点到任意一点DK1+280的方位角аi =аi 缓+△i 圆曲线上任意一点DK1+280的坐标：X i =X hy + C i cos аi Y i =Y hy + C i sin аi。

道路桩号算中边桩坐标高程计算程序道路桩号是指道路上的标志桩，用于表示道路上的位置和距离。

在道路规划、设计和施工中，需要根据桩号来确定道路的线形和纵断面，并计算出桩号对应的坐标和高程。

道路桩号的计算程序可以分为以下几个步骤：1.确定基准点：选择一个具备准确坐标和高程的点作为道路的起点，确定其坐标和高程。

2.确定桩号起点：确定一个参考点作为桩号的起点，通常选择道路的起点或其他规定的地点。

为了方便计算，可以选择一个整数作为起点桩号，如0、100等。

3.桩号计算：根据道路设计和实际情况，确定桩号的计数方式和间隔。

通常情况下，桩号以米为单位，从起点开始递增或递减。

4.桩号与坐标的关系：桩号与坐标之间存在一定的数学关系，可以根据道路的几何特征和设计参数进行计算。

例如，对于一条平直无坡道路，可以使用线性插值法计算桩号对应的坐标。

5.桩号与高程的关系：桩号与高程之间也存在一定的数学关系，可以根据道路的纵断面和地形特征进行计算。

例如，对于一条按规定坡度设计的道路，可以使用坡比法计算桩号对应的高程。

6.精度控制：在桩号计算过程中，需要考虑测量误差和计算方法的精度。

为了提高计算结果的准确性，可以采用较精确的测量方法和计算算法，并进行误差修正。

7.应用场景：道路桩号的计算程序可以应用于道路工程中的位置控制、导线布设、测量定位、横断面绘制等方面，为道路规划、建设和维护提供准确的空间位置和高程信息。

总结起来，道路桩号的计算程序是根据道路的设计和实际情况，通过选择基准点和起点桩号，确定桩号计算方式和间隔，以及桩号与坐标、高程之间的关系，计算出桩号对应的坐标和高程。

这个程序可以应用于道路工程中的各个环节，为道路的设计、施工和维护提供准确的空间位置和高程信息，提高工程质量和效率。

小白自学ExcelVBA道路坐标高程计算程序流程1声明1.1 本人是个工作时间不长测量施工员。

1.2听说想轻松愉快的干好测量得学门编程，又听说ExcelVBA是个简单的可视化的常用的……1.3 本人很懒，面对厚厚的全是“昏天暗地”代码的ExcelVBA教科书籍，实在是没心没力。

其实俺还是懂点儿VB编程的（学校学过“=”是赋值，if、while、end的意思，嘿嘿，还有画窗体图框，其他的就都还给老师了）。

1.4 谷歌、度娘，是咱的好帮手。

1.5 仰望的存在，神一般的存在，无私的化身！道路中边桩坐标计算程序120424.xls（王中伟”教授”QQ：595077）曲线坐标计算程序VBA 4.6.xls（陈超”中铁”QQ：295188316）陈师傅的VBA代码乍一看感觉比王老师的让人头痛，怎么说咱们也是个有理想有目标的人，要啃就啃硬骨头！就选~~~嘿嘿~~~王老师的程序吧！1.6 道路中边桩坐标计算程序120424.xls代码获取。

王老师、陈师傅担心不懂编程的人，不小心修改代码导致程序错误，设置了vba密码（个人认为），其次我跟他们不熟，没有正大光明的要的密码，于是求助谷歌度娘，这个关于VBA 密码破解，于是很不厚道的得到了密码，为了方便本次学习，又很不厚道的把道路中边桩坐标计算程序120424.xls的密码贴出来（df750726），神人勿怪呀！2 代码理解ing…我们通过密码打开vba代码窗口，大家可以先把每个表格打开浏览一下代码，是不是感觉很短呀，感觉应该是超级链接的意思，嘿嘿，主要的代码在模块“lx”里，打开下，代码多吧，这就是咱们要肯的骨头。

现在咱们冲呀！2.1 第一段代码理解貌似第一段代码是定义的意思“定义π值”、'定义桩号点”、“'定义坡段”、“'定义基本曲线”①Public Const PI As Double = 3.14159265358979谷歌翻译：Public Const PI As Double公共常量PI 作为双精度浮点型百度百科：双精度浮点数双精度浮点数（Double）用来表示带有小数部分的实数，一般用于科学计算，用8个字节（64位）存储空间，其数值范围为1.7E-308～1.7E+308，双精度浮点数最多有15或16位十进制有效数字，双精度浮点数的指数用“D”或“d”表示。

CASIO Fx-5800p计算器道路坐标计算程序直线程序（ZX）“FWJ=”？Q 」（“」”回车按EXE即可）“QIDZH=”？Z 」“QIDX=”？A 」“QIDY=”？B 」Lb1 0“KP=”？H 」A+(H-Z)Cos(Q)→X 」B+( H-Z)Sin(Q)→Y 」Lb1 1」“XP=”:X▲“YP=”:Y▲Lb1 2」“BIAN JU(m)=”?D」D=0＝＞Goto 0」Lb1 3」“JIA JIAO L(-),R(+)=”?E」X+D*Cos(Q+E) →M」Y+D*Sin(Q+E) →N」Lb1 4」“M［I］=”：M▲“N［I］=”：N▲Goto 2卵形曲线（LXQX）Lb1 0」“XA=”？A：“YB=”？B：“CA=”？C：“1/RA=”？D：“1/RB=”？E：“DKA=”？F：“DKB=”？G 」Lb1 1 」“DKI=”？H：“DL=”？O：“DR=”？R：H＞G ＝＞Goto 3 」（E-D）/Abs（G-F）→P：Abs（H-F）→Q：P×Q→ I：D+I →T 」C+（I+2D）×Q×90/∏→J 」“J=”：J▲C+（I/4+2D）×Q×45/（2∏）→M 」C+（3I/4+2D）×Q×135/（2∏）→N 」C+(I/2+2D) ×Q×45/∏→K 」A+Q/12×(Cos（C）+4(Cos（M）+Cos（N）)+2Cos（K）+Cos（J）) →X」B+Q/12×(Sin（C）+4(Sin（M）+Sin（N）)+2Sin（K）+Sin(J)) →Y」“X=”：X▲“Y=”：Y▲X+OCos(J-90) →U 」Y+OSin(J-90) →V 」X+R Cos(J+90) →W 」Y+RSin(J+90) →Z 」“U=”：U▲“V=”：V▲“W=”：W▲“Z=”：Z▲Goto 1说明：XA-起点X坐标；XB-起点Y坐标；CA-起点方位角；1/RA-1除以起点半径（左转为负，右转为正）；1/RB-1除以终点半径（左转为负，右转为正）；DKA-起点桩号；DKB-终点桩号；DKI-待求点的桩号；DL-左边距；DR-右边距；J-显示待求点切线方位角；X-显示待求点的中桩X坐标；Y-显示待求点的中桩Y 坐标；U-显示待求点的左边桩X坐标；V-显示待求点的左边桩Y坐标；W-显示待求点的右边桩X坐标；Z-显示待求点的右边桩Y坐标。

Fx-5800道路三维坐标计算器程序QXJS-000 主程序（文件名）“1.SZ=>NE”:“2.NE=>SZ”:?Q↙Lbl 4:?S:“PJ=”?K:Prog“QXJS-SUB0”↙Lbl 0:Q=1 => Goto1:Q=2 => Goto2:↙Lbl 1:?Z:Prog“QXJS-SUB1”: “F=”:F◢DMS◢“N=”:N◢“E=”:E◢Prog“GCJS”: Goto4↙Lbl 2: “N=”:?B: “E=”:?C:B→N: C→E:Prog“QXJS-SUB2”: “S=”:S◢“Z=”:Z◢ Goto4↙QXJS-SUB0 数据库子程序（文件名）Goto1↙同时保存多个曲线时的指针Lbl 1IF S<***（线元终点里程）:Then***→A（线元起点方位角）:***→O（线元起点里程）:***→U（线元起点X）:***→V（线元起点Y）:***→P（线元起点曲率半径）:***→R （线元终点曲率半径）: ***→L（线元起点至终点长度）: Return:IfEnd↙IF S<***:Then***→A:***→O:***→U:***→V:***→P:***→R: ***→L: Return:IfEnd↙………………………..为了便于解读，每增加一个线元增加一行语句，每增加一条曲线增加一个Lbl，每增加一个工程增加一个文件。

QXJS-SUB1 正算子程序（文件名）0.5（1÷R-1÷P）÷L→D:S-O→X↙U+∫(cos(A+(X÷P+DX2)×180÷π,0,X)→N↙V+∫(sin(A+(X÷P+DX2)×180÷π,0,X)→E↙A+(X÷P+DX2)×180÷π→F↙N+Zcos(F+K) →N:E+Zsin(F+K) →EQXJS-SUB2 反算子程序（文件名）Lbl 1:0→Z：1→Q：Prog“QXJS-SUB0”: Prog“QXJS-SUB1”↙Pol(N-B+10^(-46), E-C+10^(-46)):Isin(F-K-J) →W:S+W→S↙Abs(W)>0.0001 => Goto1↙Lbl 2: 0→Z：Prog“QXJS-SUB1”:(C-E) ÷sin(F+K) →ZGCJS(高程计算—文件名)Lbl 0：Prog“SJK”：0.005R Abs(I -L) →T：T2÷2÷R→E:Z[1]-T→A：Z[1]+T→B：S-A→C：B-S→D：S-Z[1]→E←┘If L >I：Then 1→J：Else -1→J：IfEnd←┘If S<A：Then 0→C：Else S>B =>0→D：IfEnd←┘If S≦Z[1]：Then I→Y：C→G：Else S >Z[1]=>L→Y：D→G：IfEnd←┘Lbl 1：Fix 3：“ZP=”：H+EY÷100+JG2÷2÷R◢Goto 0←┘SJK(高程计算子程序—数据库)If S≦下一竖曲线起点And S >上一竖曲线终点：Then ……→Z[1](变坡桩号): ……→H （变坡高程）: ……→R（竖曲线半径）:……→I（第一纵坡坡度，带正负号）: ……→L(第二纵坡坡度，带正负号): IfEnd←┘坐标正反算输入与使用使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当曲线半径在左时，P、R 取负值，当曲线半径在右时，P、R取正值，当曲线半径为无穷大（即直线）时，P、R以10的45次代替。