5800通用線路測量程序及編程原理

1、5800通用線路測量程序及編程原理 關(guān)于坐標(biāo)計算內(nèi)核1、程序采用數(shù)據(jù)保存在擴(kuò)充存儲的方式，更換或補(bǔ)充數(shù)據(jù)不用修改程序。2、線路采用線路號來區(qū)分，程序內(nèi)部和使用時都用自然里程參加計算。3、平曲線不用區(qū)分匝道和主線，所有平曲線都轉(zhuǎn)化成交點(diǎn)資料，此交點(diǎn)可以緩和曲線不對稱、不完整，主線通常由交點(diǎn)起算、匝道均由曲線起點(diǎn)起算。每條平曲線交點(diǎn)格式為9個數(shù)據(jù)。4、直線、圓曲線、緩和曲線分別用數(shù)學(xué)模型不同的計算子程序，保證了計算速度。直線、圓曲線子程序為精確公式，緩和曲線計算子程序為原創(chuàng)的復(fù)化段數(shù)可調(diào)的復(fù)化高斯2節(jié)點(diǎn)公式，計算精度可自由控制。（現(xiàn)暫改為固定4次復(fù)化）*數(shù)據(jù)庫編制規(guī)則基本規(guī)則：圓曲線半徑

2、右轉(zhuǎn)為正，坡度上坡為正，角度單位為度，邊樁轉(zhuǎn)角右轉(zhuǎn)為正,線路號為非負(fù)整數(shù),數(shù)據(jù)庫中的直圓里程應(yīng)該在起點(diǎn)和緩圓各輸入一個互差0.0001的里程。圓直點(diǎn)同樣處理。一、目錄部分：格式為：起點(diǎn)里程、終點(diǎn)里程、平曲線資料起始位置、豎曲線資料起始位置、坡度資料起始位置起點(diǎn)里程、終點(diǎn)里程、平曲線資料起始位置、豎曲線資料起始位置、坡度資料起始位置各條線路依次連續(xù)排列。二、平曲線部分：格式為： 圓半徑起點(diǎn)里程、緩圓里程、圓緩里程、終點(diǎn)里程、切線長、交點(diǎn)X、交點(diǎn)Y、切線方位角、圓半徑起點(diǎn)里程、緩圓里程、圓緩里程、終點(diǎn)里程、切線長、交點(diǎn)X、交點(diǎn)Y、切線方位角、圓半徑終點(diǎn)里程、圓半徑 各交點(diǎn)依次連續(xù)排列。主線按圖紙設(shè)

3、計數(shù)據(jù)整理，匝道的線元數(shù)據(jù)組合成交點(diǎn)組合時，每條曲線由直線、緩和曲線、圓曲線、緩和曲線、直線組成，對不完整緩曲要對數(shù)據(jù)進(jìn)行簡單處理，不完整緩曲可以劃歸任一側(cè)，不完整緩曲所在的起（終）點(diǎn)里程改為負(fù)，交點(diǎn)坐標(biāo)輸入為起點(diǎn)坐標(biāo)，切線長輸入0即可。三、豎曲線部分：格式為:前坡度、變坡點(diǎn)里程、變坡點(diǎn)高程、豎曲線半徑前坡度、變坡點(diǎn)里程、變坡點(diǎn)高程、豎曲線半徑后坡度、（0）各豎曲線依次連續(xù)排列。注意本程序計算單位為直坡+豎曲線，如線路最后為直坡，必須在最后一個變坡點(diǎn)里程位置填一個0，最后無直坡段可不加0*C-測量程序使用說明運(yùn)行C程序XL0？選線路號JX?輸入儀器X、JY?輸入儀器Y0-END 1-ZS 1-

4、F? 選1正算，選2反算，選0退出一、正算計算流程顯示起點(diǎn)里程、終點(diǎn)里程K?輸入K中樁里程，輸出中樁坐標(biāo)、高、放樣數(shù)據(jù)L?輸入相對中樁或上點(diǎn)偏距，按EXE輸入0到F?輸入相對線路前進(jìn)方向右角，輸出坐標(biāo)、放樣數(shù)據(jù)，到如正算里程輸入超出范圍重新選正反算時執(zhí)行反算，則反算剛計算的正算角樁坐標(biāo)     二、反算計算流程X?Y?輸入測點(diǎn)坐標(biāo), 顯示起點(diǎn)里程、終點(diǎn)里程K?輸入近似里程（默認(rèn)是線路中心里程） 輸出里程、外移距、中樁高到反算過程中隨時顯示里程逼近誤差計算過程中右下角顯示的數(shù)字為緩曲計算時復(fù)化高斯公式分段數(shù)。正反算計算過程輸入里程或任何時候中間結(jié)

5、果里程超出對應(yīng)線路里程會回到程序起點(diǎn)！*C-測量主程序Do:Deg:Norm 1:Cls:"XL 0-?"?H: /輸入線路號"JX"?V:"JY"?W:"0-END 1-Z 2-F"?U:Fix 3:/輸入置儀坐標(biāo)，正反算選項While U=1:Prog "Z":WhileEnd/循環(huán)正算流程（除非里程超限）While U=2:Prog "F":WhileEnd/循環(huán)反算流程（除非里程超限）LpWhile U0:Cls:"END"/里程超限時返回程序起點(diǎn)

6、，選項為0則退出程序D-讀ZGL:ZG+1M:ZG+2N:ZG+3O:ZG+4P/從ZG起讀5個擴(kuò)充存儲到L、M、N、O、PF-三維反算XV:YW:"X"?V:"X"?W:1.5U: /輸入測點(diǎn)坐標(biāo)Prog "K":U=3=>Retrun:2U:/近似里程，里程超限則退出Do:Prog "P":V-XI:W-YJ:I=0=>0.0001I:Pol(I,J): /求到假定點(diǎn)斜距及方位角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/計算直線段坐標(biāo)H緩和曲線計算L(P-B)÷2÷DD:/求測點(diǎn)到起點(diǎn)曲率變化D÷8M:L÷8L:1÷(3R:OA:/求拆分后各段曲率差、曲線長、高斯2節(jié)點(diǎn)系數(shù)For 1C To 7 Step 2:/循環(huán)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/求測點(diǎn)方位角I圓曲線計算LPr÷2A:2÷P*sin(A)L:O+AO:Prog "G":O+AO/用弦切角公式計算圓曲線坐標(biāo)K里程判斷5H+2G:Prog "D":U=1.5=>(L+M)÷2K:NG:OZ:/提取平、縱、橫坡要素要素指針I(yè)f Int(U)=1:Then Cls:"QD=":"ZD=":Locate 4

9、,1,L:Locate 4,2,M:?K:IfEnd:/提示起終點(diǎn)里程Abs(2K-L-M)>M-L=>3U/里程超限令正反算選項為3P-平面While K>ZG+9:G+9G:WhileEnd:/查詢符合要求的交點(diǎn) G+4G:Prog "D":-LL:MX:NY:Prog "G":1÷PP:/讀取切線長、交點(diǎn)坐標(biāo)、方位角、半徑G-5G:ZG+1Q:K-Abs QL:/讀取起點(diǎn)里程、求測點(diǎn)到起點(diǎn)長If L0:Then Prog "G":Return:IfEnd:/小于起點(diǎn)計算直線坐標(biāo)，并退出子程序0B:Q&

10、lt;0=>1÷ZGB:/計算起終點(diǎn)曲率，起點(diǎn)里程小于0起點(diǎn)曲率為上交點(diǎn)半徑倒數(shù)Abs(QD:ZG+2Q:Q-DD:/ 讀取緩圓里程、求第一緩曲全長If LD:Then Prog "H":Return:IfEnd:/如果測點(diǎn)到起點(diǎn)長小于緩一長計算緩曲坐標(biāo)，并退出子程序DL:Prog "H":K-QL:/令計算長等于緩曲長計算緩圓坐標(biāo)、求測點(diǎn)到緩圓長QD:ZG+3Q:Q-DD:/求圓曲全長，讀取圓緩里程If LD:Then Prog "I":Return:IfEnd:/如果測點(diǎn)到緩圓長小于圓曲長計算圓曲坐標(biāo)，并退出子程序

11、DL:Prog "I":K-QL:/令計算長等于圓曲長計算圓緩坐標(biāo)、求測點(diǎn)到圓緩長QD:ZG+4Q:Abs(Q)-DD:/ 求第二緩曲全長，讀取緩直里程PB:0P:Q<0=>1÷ZG+18P:-Q=ZG+10=>1÷ZG+11P/緩二起終點(diǎn)曲率If LD:Then Prog "H":Return:IfEnd:/如測點(diǎn)到圓緩小于緩二全長計算緩曲坐標(biāo)，并退出子程序DL:Prog "H":K-Abs(QL:/令計算長等于緩二長計算緩直坐標(biāo)，求測點(diǎn)到緩直長Prog "G"/計算第二直線

12、坐標(biāo)S-豎曲線ZG:G=0=>Retrun:/讀取豎曲線指針Do:Prog "D":P-LD:Oabs(D)÷DR:/讀前坡、變坡里程、變坡高、半徑、后坡，算坡差、半徑凹凸RD÷2T:G+4G:/計算切線長，完成豎曲要素準(zhǔn)備LpWhile K>M+T And ZG+10:/里程大于本豎曲終點(diǎn)且下一變坡里程不為0讀下一曲線K-MM:N+LMZ:M+TM:/計算測點(diǎn)到變坡里程長，計算前坡高，求測點(diǎn)到豎曲起點(diǎn)長If M>0:Then Z+M2÷2÷RZ:L+MD÷2÷TL:IfEnd:/測點(diǎn)在豎曲范圍，

13、計算豎曲線高、測點(diǎn)坡度M-TM:IF M>T:Then N+PMZ:PL:IfEnd/里程大于豎曲線終點(diǎn)，計算后坡高，后坡坡度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-顯示坐標(biāo)Cls:"X=":"Y=":Locate 3,1,X:Locate 3,2,Y:360Frac(360+O)/360DMS/顯示坐標(biāo)，方位角 Y-高程顯示Locate 1,4,"Z=":Locate 3,4,Z/顯示高程Z-正算Prog "K":U=3=>Retrun:/ 輸入里程，里程超限則退出Prog "P":Prog "X"/顯