CASIO程序与全站仪在油田路基坡脚线、坡面放样中的应用

史上最全的casio测量放样、坐标放点、导线放1 XY边长及方位角反算Prog“A":L=J:Prog"B"Fix 3:I"D="Lbl 1:{S}:V"X"=X+Rec(S, L U"Y"=Y+J Goto 1输入：X1,Y1,X2,Y2.输出: L°=,D=.2 XY-1输入里程桩号计算中间点坐标Prog"A":L=J:Prog"B"Fix 3:I"D="{Z}: Z:Lbl 1:{G}:V"X"=X+Rec(Abs(G-Z, LU"Y"=Y+J Z-起始点里程桩号Goto 1Ｇ－任一点里程桩号测量程序.xls 第1 页3 JD两直线交点坐标计算Prog"A": Prog"C":L"L1":Prog"F":E=W+K:{L}:L"L2":Prog"F":C=V+K:S=tan E:T=tan C:G=S-T:A=A+0.0001 U"X"=(SA-MT-B+N)/G Z"Y"=S(U-A)+BPol(U-M,Z-N:I"D="输入X1,Y1,X2,Y2,X3,Y3. X4,Y4,L1,L2,输出X=,Y=,D=,4 ZBFD坐标放点Prog"A":Fix 3:{MN}:Pol(M"X3"-A,N"Y3"-B: I"D="J=J-W+180:J0=J=J+360L=J:Prog"B"Goto 3输入:X1,Y1,X2,Y2,X3,Y3. 输出:D,L°测量程序.xls第2 页5 Z自由导线推算E=0:1F=0:Prog"A":Fixm:C=A:D=B:J[F]=J:{LH}:L"L0":Prog"F":I[F]=J[F]+K-180:H[E+F]=H"D":Prog"R":C"X="D"Y="J[F]=I[F]:Goto 46 ZDX支导线计算X1Y1,X2Y2,为已知点坐标Prog"A":L°-角度用小数输入Lbl 4:输一组L°、D.得一个点的坐标{LS}:L:Prog"F":M"X3"=A+Rec(S,W+K+180N"Y3"=B+JGoto 4输入:X1,Y1,X2,Y2,L,S.输出:X3,Y3.L°―观测角度以下按小数点输入测量程序.xls第3 页7 DX单一导线简易平差计算Prog"A":Prog"C":J=W:Ｉ=V:C=A-M:D=B-N:F=E"N":I[F]=J:Lbl 1:F=F-1:{GH}:I[F]=J[F]+G"L°"-180:I[F]0=I[F]=I[F]+360:≠I[F]360=I[F]=I[F]-360H[E+F]=H"D":F1=Goto 1{G}:G=I[F]+G"L°"-180:G0=G=G+360:≠G360=G=G-360I=G-I:W=12″[S]=620.790G"W"=602FracI 相对中误差:fS/[S]=1/*****fX=-0.007fY=0.000G=0:Lbl 2:I[F]=I[F]-I(E-F)/E:Prog"R":说明:F=F+1:N-观测角数,L°-左角按度分秒输入G=H+G:W-角度闭合差,fX,fY坐标闭合差F Goto 2若用外(右)角计算闭合差程序中之I=G-I-360G"[S]="若用内(左)角计算闭合差程序中之I=G-II"fX"=C 本程序只能计算八个未知点,多于八个时每增加一个C=A:未知点,必需扩展两个储存器.A"fY"=D D=B:B"fS/[S]"=Int(G/ √ (I I+AA Lbl 3: F=F －1:Prog"R":C"X"=C-I H/G D"Y"=D-AH/G F1=Goto 3测量程序.xls 第4 页8 WDX 无定向导线Prog"A":D=I:C=0:M=0:N=0:P:Lbl 0:P=P-1:{LS}:L"B°":Prog"F":F=K:C=C+F:C=0=Goto 1C=C+180:Goto 1X1Y1,X2Y2,-已知点坐标Lbl 1:P-边数,B°-观测角,第一个角输入0,Rec(S,C:S-边长,Xf,Yf-增量累加数,记录侍用Fix 4:输入一组Xf,Yf,求得一组X Y,M"Xf"=M+I 用左角,度以下按小数输入.N"Yf"=N+JP0=Goto 0Pol(-M,-N:R=W+180-J:K=D/I:V=Cos R:点号观测角边长Xf Yf X Y E=Sin R:Lbl 2:{HO}:H"Xf":O"Yf":Fix 3:T"X"=X+K(HV-OE Z"Y"=Y+K(HE+OV Goto 2B 为观测角,第一个角输入04119.***-*****.***-*****.***-*****.***-*****.***-*****.630142.*****.8800-66.*****.*****.0000-211.5145-303.*****0°00′00″ 42 28 06 277 33 42 A 1164.*****.*****.***-*****.952 测量程序.xls 第5 页9 QJD前交点(正弦公式)Prog"A":Prog"J":S=I Sin E/Sin(E+K:Prog"O"A°、B°-观测角,X1、Y1、X2、Y2-正切公式Prog"A":Prog"J":S=tan E:T=tan K:G=S －1＋T -1:M"X3"=(X/T+A/S-Y+B)/G N"Y3"=(Y/T+B/S+X-A)/G 输入:X1,Y1.X2,Y2,A°,B°,输出:X3Y310 FX 方向交会Prog"A":Prog"J":S=tan E:T=tan K:G=S-T:M"X3"=(SX-AT-Y+B)/G N"Y3"=S(M-X)+YX1Y1,X2Y2,为起算点坐标.A°,B°-方位角.X,Y--交点坐标11 CJ 侧方交会点Prog"A":Prog"J":S=I sin(E+K)/Sin K:{K}:K=0=M"X"=A+Rec(S,W+E-180≠N"Y"=B+JX1Y1,X1'Y1',X2Y2--已知点坐标,A°B°-观测角,左右分开.左K=0,右K=1A=63° 33′06″B=92° 20′51″从左右两组算出坐标取中数为XYA ′=109° 59′39″B ′=50° 57′06″测量程序.xls 第6 页12 HJ后交点计算Prog"A":M"X3":X1=992.692X2=1012.499N"Y3":Prog"J":S=tan E:T=tan K:H=(Y-B)/S+X-A:F=(X-A)/S-Y+B:Q=(N-B)/T-M+A:V=(M-A)/T+N-B:K=(H+Q)/(F+V:X3=1063.248D=(H-FK)/(1+K2:Y4=999.997Y3=1000.000 O"X4"=A+DX1Y1,X2Y2,X3Y3---已知点坐标P"Y4"=B+KD A°、B°----观测角X4Y4--测点坐标12HJ―1ProgAUD=IU（另一种后Pol（A－MX3，B－NY3U交公式）S=IUG=JUProgJUZ=180－（W－G＋E＋K）÷2UQ=tan－1（D Sin C÷S Sin EUH=tan－1（tan Z÷tan（45＋QUR=Z＋HUT=D Sin（E＋R）÷Sin EU OX=X＋T Cos（W＋RPY=Y＋T Sin（J＋R测量程序.xls第7 页13 LDHJ两点后交Prog"A": L"A°": Prog"F":C=K:{L}:X2=857.025X1=944.202 L"B°":Prog"F":O=K:{L}:L"E°":Prog"F":E=K:{L}:L"F°":Prog"F":F=K:Y3=139.053N=Sin O/Sin(O+F:P=√(S2+M2-2SM Cos(C-O:观测角A°B°E°F°度以下按小数点输入.T=√(D2+N2-2DN Cos(F-E:K"D1"=2I"D0"/(P+TS"S1"=KSD"S2"=KDM"S3"=KMN"S4"=KNH=Cos-1((I2+S2-M2)/2IS:L=H:Prog"B":Z"H°="G=Cos-1((I2+N2-D2)/2IN:L=G:Prog"B":Z"G°="Fix 3:Q"X3"=A+Rec(S,W+H+180R"Y3"=B+JQ"X4"=A+Rec(D,W+E-F-LR"Y4"=B+J测量程序.xls第8 页Prog"A":Lbl 0:Prog"W":Fix 2:S=D(Sin K)2H=S/tan K+G{L}:Prog"F":输入:X1,Y1,X2,Y2,H0-起点高,T仪高, Prog"O"D-视距,S-平距,Z°―天顶距,P标高Goto 0输出:X3,Y3,S平距15 JZB极坐标法测碎部点坐标、高程Prog"A":G=E"H0"+T:Lbl 8:{LSV}:Prog"F":Prog"O"H=G+V"h"-PGoto 8输入:X1,Y1,X2,Y2, H0,T,L°，S，h求:X,Y,H.测量程序.xls第9 页X3=273.***** SD 垂距,垂点坐标计算Y3=208.502Prog"A":O=Sin J:P=CosJ:E=M"X3"-X:F=N"Y3"-Y:S=EP+FO D=EO-FPY=252.157V"X"=X+Rec(S,W U"Y"=Y+J注:D 为正时,3点在1-2之左侧;S 为正,3在1-2之线上;为负则在2--117 MM 坐标法计算多边形面积N: P=0: V=0:A"X":B"Y":E=A: F=B:Fix 3:Lbl 0:{XY}:X: Y:N=4(边数)Lbl 1:Pol(X-A,Y-B: I"S"P=P+(Y-B)(X+A)/2:A=X:B=Y:V=V+1:V ≠N=Goto 0P"m 2"=Abs P Q"mu"=1.5 m P测量程序.xls 第10 页18断面图封闭图形面积计算(高程相当于X,距离相当于Y,用相对坐标计算,)N:P=0:V=0:Y=0:A"H":B"D":E=A:69.8F=B: Fix 3:Lbl 0:N=V+1=H=E: Y=F: Goto 1｛HD ｝:H:Y=Y+D:Lbl 1:m =32.795A=H:B=Y:N=多边形边数,V=V+1:H=转折点高程V ≠N= Goto 0D=相邻两点距离P"m 2"=Abs P第一点D 输入0。

卡西欧5800边坡测量程序U″XO″：V″YO″：A″X1″：B″Y1″：POL（G″X2″-A,D″Y2″-B：W=J：I：LBI1：{LKM}：X″XA="=A+LCOSW+KCOS（W+M)◢Y"YA="=B+LSinW+Ksin（W+M）◢POL（S-U),Y-U：J <0=>J=J+36O⊿X″DIX=″=I◢FiX3：O″FW=″=J◢Goto15800边坡测量程序″XO=″?U：″YO=″?V：″X1=″?A：″Y1″?B：″X2″?G:″Y2″?D↙POL（G-A,D-B)：J→W:Fix 3:CLS↙Lbi 1 ↙″L=″?L:″K=″?K:″M=″?M↙A+LCOSW+KCOS（W+M)→X↙B+LSinW+Ksin（W+M)→Y↙Cls↙"XA=":Locate 5,1,X:"YA=":Locate 5,2,Y◢Pol（X-U,Y-V)：If J <0:Then J+36O→J :Else J→J↙Cls↙″DIX=″:Locate 5,1,I:″FW=″:J▲DMS◢Goto 1程序特点：全线标高数据一次性程序化输入，扩充变量数据库，无需修改程序内容；全线贯通计算，标高计算时将全线数据分为4等分，首先判断待求点在哪个区间，再循环选择，提高运算时间；傻瓜操作，适用初级用户。

一、程序：M CW“1.BG 2.SZ ”：W=1=> Goto 1ΔO“KOU LING”：O≠1 23456=>O=0：Goto ５ΔO=0：Z[1]=0: Goto 0←┘Lbi 0←┘”No.”：Z[1]+1◢{ZHPR}：Z“D0” ：H“BG” ：P“PD” ：R“R0” ：Z[Z[1]×4+ 3]=Z：Z[Z[1]×4+4]=H：Z[Z[1]×4+5]=P：Z[Z[1]×4+6]=R：“NEXT”◢Isz Z[1]：Goto0←┘Lbi 1←┘Z[2]= Z[1]←┘{D}：D≤Z[Int(Z [2]÷4+1)×4+3]=>Z[2]= Int(Z [2]÷4) ：Goto2ΔD≤Z [Int (Z[2] ÷2+1)×4+3]=>Z[2]=Int(Z[2]÷2)：Goto2ΔD≤Z [Int (Z[2]×.75+1)×4+3]=>Z[2]=Int (Z[2]×.75)：Goto2ΔGoto 2←┘Lbi 2←┘D<Z[Z[2]×4+3]=>Dsz Z[2]：Goto 2ΔGoto 3←┘Lbi 3←┘Z[Z[2]×4+6]≠0=>G=Z[Z[2]×4+4]+(D-Z[Z[2]×4+3])×Z[Z[2]×4+ 5]÷100+(D-Z[Z[2]×4+3]) ^2÷2÷Z[Z[2]×4+6]：Goto 4ΔG=Z[Z [2]×4+4]+(D-Z[Z[2]×4+3])×Z[Z[2]×4+5]÷100：Goto 4←┘Lbi 4←┘{Z}：Z：G=G-0.02Abs(Z):“G”: G◢（4850输出）{Z}：Z：G=G-0.02Abs(Z): G“G”◢（4800输出）Goto 1←┘Lbi 5←┘“OUT”◢程序操作说明程序分为两部分：1.BG为标高计算，2.SZ为参数设置。

CASIOFX—4800P计算器在公路测量中的应用中铁七局三公司邢汾高速公路L13标项目部胡会杰二0一一年三月二十日CASIO Fx-4800p计算器在公路测量中的应用随着高速公路工程建设的深入发展，建设中对测量施工放样精度要求越来越高，公路工程施工线形具体有平曲线和竖曲线。

平曲线有直线段、圆曲线段和缓和曲线段。

在一般情况计算，计算慢和计算时产生累计误差，在精度方面已达不到施工要求，怎样使计算快速、精度得以保证，成为测量中一项重要工作，针对并结合本人多年来的工作实践经验，简述在Fx-4800p计算器中编制的程序在公路建设测量中的应用。

并与各位同行共同探讨。

敬请各位同仁提出宝贵意见。

一：计算原理1-1：平曲线计算1-1-1：直线段坐标计算通过坐标正算求得，根据已知点坐标、边长和方位角，求未知点坐标，如图(1-1)可知ΔX ab=S ab×COS αabΔY ab=S ab×Sin αabX b=X a+Δx abY b=Y a+Δy abαab ab方位角S ab ab两点之间长度图(1-1) 1-1-2：圆曲线段坐标计算圆曲线坐标计算有好几种计算方法，在这使用圆心角计算法，则圆曲线上相邻点间的弧距为L，其所对的圆心角为θ、半径和弧距的关系式为：θ=L/R×180/π通过圆心与切点的方位角，可求得圆心点到任一圆曲线点的方位角，应注意圆曲线走向，右转θ为正，左转θ为负，公式如下，如图(1-2)可知αi=αzy±θi zy式中：αzy切点与圆心点的方位角JDθθθo θyz图(1-2) 该方法在计算边桩时只变动半径长度，方位角不变，计算方便，在计算时首先求圆心点到任一圆曲线点的方位角。

以图(1-2)可以看出，坐标正算求得，未知点坐标。

1-1-3：缓和曲线段坐标计算缓和曲线段计算一般使用近似方法，具体内容见«工程测量学»，公式如下：X i=L i-L i 5 / (40R²L o²)Y i=L i³ / (6RL o)-L i7 / (336R³L o³)L i缓和曲线任一段长L o缓和曲线总长缓和曲线坐标公式为施工坐标系，要转换测量坐标系，计算公式如下：X p=X o+X i×COSα-Y i×SinαY p=Y o+Y i×COSα-X i×SinαX o,Y oX i,Y i缓和曲线相对坐标α转换方位角1-1-4：放样数据计算利用坐标反算，根据已知两点坐标、求出两点间方位角和距离，公式如下：tgαab=ΔY ab/ΔX ab=(Y b-Y a) / (X b-X a)S ab=Δy ab / Sinα=ΔX ab / COSαab程序编制方法：根据线路的已知曲线类型,把已知曲线条件，曲线划分界线、判断条件，输进主程序,根据桩号判断采用的数据，通过共用子程序，进行任意点的坐标计算。

CASIO fx-4800P(4500PA)计算器程序清单程序说明：本工程可以计算公路平面曲线中桩、边桩坐标（只限于缓和曲线长度相同之曲线）。

本工程中行首数字为行号，无须输入计算器。

程序清单如下：GLZBJS 公路坐标计算主程序（文件名）10 Defm 5 扩展内存20 F “JSFWJ”：R:S“LS”：A“ZJ”：输入曲线计算要素：F-计算方位角；I“L-,R+”：K“JDZH”：M“XJD”：R-曲线半径：N“YJD”S-缓和曲线长度；A-转角值；I-转角判30 C=28.6479S/R 定，当左转时输入+1，右转时输入-1；40 D=S2/(24R)-S^4/(2384R^3) K-JD桩号；M\N-JD坐标50 Q=S/2-S^3/(240R2) C-计算缓和曲线角；D-曲线内移植；60 T=(R+D)tan(A/2)+Q◢ Q-切线增值70 L“LQX”=（A-2C）Rπ/180+2S◢ T-切线长； L“LQX”-曲线长度80 E=（R+D）/cos(A/2)-R◢ E-曲线外距。

90 O“LY”=Rπ(A-2C)/180 O“LY”-圆曲线长度100 Z[1]“ZH”=K-T◢ ZH点桩号（如S=0,则为ZY点桩号）110 Z[2]“HY”=(K-T)+S◢ HY点桩号（同上）120 Z[3]“QZ”=K-T+L/2◢ QZ点桩号130 Z[4]“YH”=K-T+S+O◢ YH点桩号（如S=0,则为ZY点桩号）140 Z[5]“HZ”=K-T+L◢ HZ点桩号（同上）150 Lbl A 循环标识160 {W，B，C}输入W-1计算点桩号；B-宽度；C与170 W﹤K-T=﹥Prog“ZX SUB”：前进方向（轴线）的夹角。

≠=﹥W﹤(K-T)+S=﹥Prog“HHQX1SUB”：≠=﹥W﹤K-T+S+O=﹥Prog 判别输入的桩号在直线段、第一缓“YQX SUB”：≠和曲线段、圆曲线段还是第二缓和曲线=﹥W﹤K-T+O-2S=﹥Prog“HHQX2 段，然后分别用各子程序进行计算。