5800通用线路测量程序及编程原理

1、5800通用线路测量程序及编程原理 关于坐标计算内核1、程序采用数据保存在扩充存储的方式，更换或补充数据不用修改程序。2、线路采用线路号来区分，程序内部和使用时都用自然里程参加计算。3、平曲线不用区分匝道和主线，所有平曲线都转化成交点资料，此交点可以缓和曲线不对称、不完整，主线通常由交点起算、匝道均由曲线起点起算。每条平曲线交点格式为9个数据。4、直线、圆曲线、缓和曲线分别用数学模型不同的计算子程序，保证了计算速度。直线、圆曲线子程序为精确公式，缓和曲线计算子程序为原创的复化段数可调的复化高斯2节点公式，计算精度可自由控制。（现暂改为固定4次复化）*数据库编制规则基本规则：圆曲线半径

2、右转为正，坡度上坡为正，角度单位为度，边桩转角右转为正,线路号为非负整数,数据库中的直圆里程应该在起点和缓圆各输入一个互差0.0001的里程。圆直点同样处理。一、目录部分：格式为：起点里程、终点里程、平曲线资料起始位置、竖曲线资料起始位置、坡度资料起始位置起点里程、终点里程、平曲线资料起始位置、竖曲线资料起始位置、坡度资料起始位置各条线路依次连续排列。二、平曲线部分：格式为： 圆半径起点里程、缓圆里程、圆缓里程、终点里程、切线长、交点X、交点Y、切线方位角、圆半径起点里程、缓圆里程、圆缓里程、终点里程、切线长、交点X、交点Y、切线方位角、圆半径终点里程、圆半径 各交点依次连续排列。主线按图纸设

3、计数据整理，匝道的线元数据组合成交点组合时，每条曲线由直线、缓和曲线、圆曲线、缓和曲线、直线组成，对不完整缓曲要对数据进行简单处理，不完整缓曲可以划归任一侧，不完整缓曲所在的起（终）点里程改为负，交点坐标输入为起点坐标，切线长输入0即可。三、竖曲线部分：格式为:前坡度、变坡点里程、变坡点高程、竖曲线半径前坡度、变坡点里程、变坡点高程、竖曲线半径后坡度、（0）各竖曲线依次连续排列。注意本程序计算单位为直坡+竖曲线，如线路最后为直坡，必须在最后一个变坡点里程位置填一个0，最后无直坡段可不加0*C-测量程序使用说明运行C程序XL0？选线路号JX?输入仪器X、JY?输入仪器Y0-END 1-ZS 1-

4、F? 选1正算，选2反算，选0退出一、正算计算流程显示起点里程、终点里程K?输入K中桩里程，输出中桩坐标、高、放样数据L?输入相对中桩或上点偏距，按EXE输入0到F?输入相对线路前进方向右角，输出坐标、放样数据，到如正算里程输入超出范围重新选正反算时执行反算，则反算刚计算的正算角桩坐标     二、反算计算流程X?Y?输入测点坐标, 显示起点里程、终点里程K?输入近似里程（默认是线路中心里程） 输出里程、外移距、中桩高到反算过程中随时显示里程逼近误差计算过程中右下角显示的数字为缓曲计算时复化高斯公式分段数。正反算计算过程输入里程或任何时候中间结

5、果里程超出对应线路里程会回到程序起点！*C-测量主程序Do:Deg:Norm 1:Cls:"XL 0-?"?H: /输入线路号"JX"?V:"JY"?W:"0-END 1-Z 2-F"?U:Fix 3:/输入置仪坐标，正反算选项While U=1:Prog "Z":WhileEnd/循环正算流程（除非里程超限）While U=2:Prog "F":WhileEnd/循环反算流程（除非里程超限）LpWhile U0:Cls:"END"/里程超限时返回程序起点

6、，选项为0则退出程序D-读ZGL:ZG+1M:ZG+2N:ZG+3O:ZG+4P/从ZG起读5个扩充存储到L、M、N、O、PF-三维反算XV:YW:"X"?V:"X"?W:1.5U: /输入测点坐标Prog "K":U=3=>Retrun:2U:/近似里程，里程超限则退出Do:Prog "P":V-XI:W-YJ:I=0=>0.0001I:Pol(I,J): /求到假定点斜距及方位角J-OJ:Rec(I,J):K+IK:/得出偏距和修正后里程Prog "K":U=3=>Retr

7、un:LpWhile Abs（I）>1m/里程超限则退出、修正值大于1mm则重新试算Prog "W":Prog "S":Prog "Y"/计算里程、位置、求中桩高程G直线计算X+Lcos(OX:Y+Lsin(OY/计算直线段坐标H缓和曲线计算L(P-B)÷2÷DD:/求测点到起点曲率变化D÷8M:L÷8L:1÷(3R:OA:/求拆分后各段曲率差、曲线长、高斯2节点系数For 1C To 7 Step 2:/循环4次A+Lr(C+R)(M(C+R)+BO:Prog "G&

8、quot;:/计算每段的第一部分A+Lr(C-R)(M(C-R)+BO:Prog "G":Next:/计算每段的第二部分A+8Lr(D+BO/求测点方位角I圆曲线计算LPr÷2A:2÷P*sin(A)L:O+AO:Prog "G":O+AO/用弦切角公式计算圆曲线坐标K里程判断5H+2G:Prog "D":U=1.5=>(L+M)÷2K:NG:OZ:/提取平、纵、横坡要素要素指针If Int(U)=1:Then Cls:"QD=":"ZD=":Locate 4

9、,1,L:Locate 4,2,M:?K:IfEnd:/提示起终点里程Abs(2K-L-M)>M-L=>3U/里程超限令正反算选项为3P-平面While K>ZG+9:G+9G:WhileEnd:/查询符合要求的交点 G+4G:Prog "D":-LL:MX:NY:Prog "G":1÷PP:/读取切线长、交点坐标、方位角、半径G-5G:ZG+1Q:K-Abs QL:/读取起点里程、求测点到起点长If L0:Then Prog "G":Return:IfEnd:/小于起点计算直线坐标，并退出子程序0B:Q&

10、lt;0=>1÷ZGB:/计算起终点曲率，起点里程小于0起点曲率为上交点半径倒数Abs(QD:ZG+2Q:Q-DD:/ 读取缓圆里程、求第一缓曲全长If LD:Then Prog "H":Return:IfEnd:/如果测点到起点长小于缓一长计算缓曲坐标，并退出子程序DL:Prog "H":K-QL:/令计算长等于缓曲长计算缓圆坐标、求测点到缓圆长QD:ZG+3Q:Q-DD:/求圆曲全长，读取圆缓里程If LD:Then Prog "I":Return:IfEnd:/如果测点到缓圆长小于圆曲长计算圆曲坐标，并退出子程序

11、DL:Prog "I":K-QL:/令计算长等于圆曲长计算圆缓坐标、求测点到圆缓长QD:ZG+4Q:Abs(Q)-DD:/ 求第二缓曲全长，读取缓直里程PB:0P:Q<0=>1÷ZG+18P:-Q=ZG+10=>1÷ZG+11P/缓二起终点曲率If LD:Then Prog "H":Return:IfEnd:/如测点到圆缓小于缓二全长计算缓曲坐标，并退出子程序DL:Prog "H":K-Abs(QL:/令计算长等于缓二长计算缓直坐标，求测点到缓直长Prog "G"/计算第二直线

12、坐标S-竖曲线ZG:G=0=>Retrun:/读取竖曲线指针Do:Prog "D":P-LD:Oabs(D)÷DR:/读前坡、变坡里程、变坡高、半径、后坡，算坡差、半径凹凸RD÷2T:G+4G:/计算切线长，完成竖曲要素准备LpWhile K>M+T And ZG+10:/里程大于本竖曲终点且下一变坡里程不为0读下一曲线K-MM:N+LMZ:M+TM:/计算测点到变坡里程长，计算前坡高，求测点到竖曲起点长If M>0:Then Z+M2÷2÷RZ:L+MD÷2÷TL:IfEnd:/测点在竖曲范围，

13、计算竖曲线高、测点坡度M-TM:IF M>T:Then N+PMZ:PL:IfEnd/里程大于竖曲线终点，计算后坡高，后坡坡度V-放样距离角度Cls:"L=":"<=":Pol(X-V,Y-W:Locate 3,1,I:/计算并显示放样距离J<0=>J+360J:JDMS/显示放样方位角W-位置显示Cls:Fix 3:"K=":"J=":Locate 3,1,K:Locate 3,2,J/显示里程、偏移值X-显示坐标Cls:"X=":"Y=":Locate 3,1,X:Locate 3,2,Y:360Frac(360+O)/360DMS/显示坐标，方位角 Y-高程显示Locate 1,4,"Z=":Locate 3,4,Z/显示高程Z-正算Prog "K":U=3=>Retrun:/ 输入里程，里程超限则退出Prog "P":Prog "X"/显