用CASIO fx4800p编高斯投影坐标转换程序
卡西欧fx-4800P计算器程序集(小刚)
计算器对应程序名称：{ZDZB}（转点坐标）A“XZ”B“YZ”C“XH”D“YH”：{RL}I=1：J=1：Pol（（C－A），（D－B））：F=R＋JRec（L，F）X“X”=I＋A◢Y“Y”=J＋B◢注：用于已知两点坐标，及待测点与测站点的连线与已知边（两已知点的连线）的水平夹角和待测点至测站点的水平距离求待测点坐标。
XZ、YZ—测站点坐标；XH、YH—后视点坐标；R—水平角（顺时针为正，逆时针为负）；L—待测点至测站点水平距离。
计算器对应程序名称：{FWJ}(方位角)X“X1”：Y “Y1”：U“X2”：V“Y2”G=U-X:H=V-YPol(G,H):I◢J<0==>J=+360⊿J“J0=”注：用于已知两点坐标，求两点之间的方位角和距离。
X1、Y1—前视（后视）点坐标；X2、Y2—测站点坐标；L—测站点至前视（后视）点距离；J0—方位角(按SHIFT键再按°′″键转换为60制角度)。
任意多边形面积计算程序fx-4800p（三）计算器对应程序名称：{MJ}Lbl 0:M“N”：N“E”：E=M：F=N：S=0Lbl 1:{XY}:S=S-(XN-YM)/2:X=E==>Y=F==>“㎡”:S ◢Goto O ⊿⊿M=X:N=Y:Goto 1注：用于计算任意多边形的面积。
N-起始点X坐标（闭合点X坐标）E-起始点Y坐标（闭合点Y坐标）X-前进点X坐标Y-前进点Y坐标S-多边形面积（取绝对值）应顺序依次输入点坐标，输完所有点坐标后应再次输入起点坐标，这样就形成了多边形的闭合，程序将自动计算出多边形的面积。
极坐标测设建筑物（放样）计算程序fx-4800p（四）计算器对应程序名称：{ZBCS}E:F:K:LG=K-E:H=L-FPOL(G,H):J<0==>B=J+360：≠==>B=J⊿Lbl 0{XY}:X=0==>G o t o E⊿G=X-E:H=Y-FPol(G,H):I◢J<0==>J=J+360⊿C=J-B:C<0==>C=C+360⊿C◢Goto 0Lbl E注：用于计算已知方向边AB（已知点A、B两点的连线）与A（测站点）点至任意放样点i方向边（测站点与待放样点的连线）的水平夹角Bi和水平距离。
卡西欧计算器坐标计算程序4800
CASIO fx-4800P计算器程序DA DI——N E (公路弯道坐标计算)K“JD”：R：S“LS”：F“PJ——L－R＋”：A“ZH——N”：B“ZH——E”：C“JD——N”：D“JD——E”M“ZH”=K－（R＋S2÷（24R））tan（Abs F÷2）－S÷2＋S^3÷（2 40R2）▲N“HY”=M＋S▲O“QZ”=M＋（πRAbs F÷180＋S）÷2▲Z[1]“YH”=M＋πRAbs F÷180▲Q“HZ”= Z[1]＋S▲T=K－M▲L=Q－M▲E=（R＋S2÷（24R））÷cos（F÷2）－R▲Lb1 0：{G}：G“ZHUANG HAO”：G＜O＝=＞Z=G－M≠=＞Z=Q－G：△{U}：U“I——B OUT＋IN－”Z＜S＋0.005＝=＞H=tan-1（24Z2S2R2－Z^6）÷（48Z^3R^3－6Z^4 SR））：V=Z－Z^5÷（40S2R2）＋Usin H：W=Z^3÷（6SR）－Z ^ 7÷（336Z^3R^3）－Ucos H：≠=＞H=180（Z－S÷2）÷（πR）：V=Rsin H＋S÷2－S^3÷（240R2）＋Usin H：W=R－Rcos H＋S2÷（24R）－Ucos H：△F＜0＝=＞G＜O＝=＞X=V：Y=W：≠=＞X=T＋Tcos F－Vcos F－Ws in Abs F：Y=Tsin Abs F－Vsin Abs F＋Wcos F：△≠=＞G＜O＝=＞X=V：Y=－W：≠=＞X=T＋Tcos F－Vcos F－Wsin Abs F：Y=Vsi n Abs F－Tsin Abs F－Wcos F：△△C=A＝=＞D≥B＝=＞P=90：≠=＞P=270：△≠=＞C＜A＝=＞P= tan-1（（D－B）÷（C－A））＋180：≠=＞P= tan-1（（D－B）÷（C－A））：△△I“COM——N”=A＋√（X2＋Y2）cos（P－tan-1（Y÷X））▲J“COM——E”=B＋√（X2＋Y2）sin（P－tan-1（Y÷X））▲G=G＋20Goto 0程序运行输入：JD？——交点桩号R？——圆曲线半径LS？——缓和曲线长度PJ——L－R＋？——偏角右偏为正ZH——N？——ZH点N坐标ZH——E？——ZH点E坐标JD——N？——JD点N坐标JD——E？——JD点E坐标输出：ZH= ——ZH点桩号HY= ——HY点桩号QZ= ——QZ点桩号YH= ——YH点桩号HZ= ——HZ点桩号T= ——切线长L= ——曲线长E= ——外矢距再输入：ZHUANG HAO？——输入任一点桩号I——B OUT＋IN－？——距中线距离（外侧为正）再输出：COM——N= ——计算出N坐标COM——E= ——计算出E坐标以此循环——再输入——再输出程序中各字母所代表的含义：A——ZH点N坐标B——ZH点E坐标C——JD点N坐标D——JD点E坐标E——外矢距F——偏角（右偏为正）G——弯道内任一点桩号H——弯道内任一点切线角I——计算出的N坐标J——计算出的E坐标K——交点桩号L——曲线长M——ZH点桩号N——HY点桩号O——QZ点桩号P——ZH--JD方位角Q——HZ点桩号R——圆曲线半径S——缓和曲线长度T——切线长U——距中线距离（外侧为正）V——支距W——支距X——弯道内坐标Y——弯道内坐标Z——任一点曲线长度Z[1]——YH点桩号。
CASIO fx-4500P、4800P计算器程序
CASIO fx－4500P与CASIO fx－4800P计算器部分按键及程序语句不同之处比较【坐标反算程序】文件名：ZBFS（调用子程序文件名：Z.FS）L1 Lbl 0L2 {XY}:X″X1″:Y″Y1″:U″X0″:Z″Y0″L3 Prog ″Z.FS″L4 ″DIST=″:D:Pause 0L5 D=0=>Goto 0:◣A″ANGLE=″◢Goto 0使用说明屏幕提示执行操作屏幕显示计算结果X1? 输入视点坐标XDIST= 置镜点与视点的距离Y1? 输入视点坐标YX0? 输入置镜点坐标XANGLE= 置镜点→视点方位角Y0? 输入置镜点坐标Y注意fx-4500P型计算器将L3行改成：Prog Z.FS【坐标正算程序】文件名：ZBZSL1 Lbl 0L2 {DA}:D″DIST″:A″ANGLE″:U″X0″:Z″Y0″L3 ″X=″:X=U+Dcos A:Pause 0L4 Y=Z+Dsin A◢Goto 0使用说明屏幕提示执行操作屏幕显示计算结果DIST? 输入置镜点与视点的距离X= 视点坐标X ANGLE? 输入置镜点→视点方位角X0? 输入置镜点坐标XY= 视点坐标Y Y0? 输入置镜点坐标Y 【转角放样程序】文件名：ZJFY（调用子程序文件名：Z.FS）L1 Lbl 0L2 {STUZ}:X=S″HSX″:Y=T″HSY″:U″X0″:Z″Y0″L3 Prog ″Z.FS″L4 D=0=>Goto 0:◣J=AL5 Lbl 1L6 {XY}:X″QSX″:Y″QSY″L7 Prog ″Z.FS″L8 ″DIST=″:D:Pause 0L9 D=0=>Goto 1:◣A≥J=>A″R-ANGLE″=A-J◢≠=>A″R-ANGLE″=A-J+360◢◣Goto 1使用说明屏幕提示执行操作屏幕显示计算结果HSX? 输入后视点坐标XDIST= 置镜点与前视点距离HSY? 输入后视点坐标YX0? 输入置镜点坐标XY0? 输入置镜点坐标YR-ANGLE=置镜点→后视点到置镜点→前视点的转角QSX? 输入前视点坐标XQSY? 输入前视点坐标Y注意fx-4500P型计算器将L3、L7行改成：Prog Z.FS【单交点基本型（对称、非对称）平曲线要素、中桩、边桩坐标计算程序】文件名：PQXJS（调用子程序文件名：PQXZB）L1 R:H″LS1″:H=0=>I=0:◣I″LS2″:Z:D″JD″:S=Abs Z L2 N=H2÷24R-H^4÷2384R^3L3 P=I2÷24R-I^4÷2384R^3L4 K=H÷2-H^3÷240R2L5 M=I÷2-I^3÷240R2L6 B=H÷2R:C=I÷2RL7 T″T1″=(R+P-(R+N)cos S)÷sin S+K◢U″T2″=(R+N-(R+P)cos S)÷sin S+M◢L=H+I+(Sπ÷180-B-C)R◢H=I=>E=(R+N)÷cos (S÷2)-R◢J=2T-L◢◣N″LY″=L-H-I◢B″ZH″=D-T◢C″HY″=B+H◢E″YH″=C+L-H-I◢F″HZ″=E+I◢G″QZ″=C+N÷2◢Prog ″PQXZB″使用说明屏幕提示执行操作屏幕显示计算结果R? 输入圆曲线半径T1= 小里程段切线长LS1?输入小里程段缓和曲线长（不设缓和曲线时，输入0）T2= 大里程段切线长LS2?输入大里程段缓和曲线长（不设缓和曲线时，不提示）L= 平曲线总长Z? 输入转向角(左转负、右转正) E= 外距(T1=T2时显示)JD? 输入交点桩号J= 校正数(T1=T2时显示)LY= 圆曲线长ZH= 直缓点桩号不设缓和曲线时，ZH=HY，为直圆点桩号HY= 缓圆点桩号YH= 圆缓点桩号不设缓和曲线时，YH=HZ，为圆直点桩号HZ= 缓直点桩号QZ= 曲中点桩号XJD? 输入交点的X坐标XZH= 直缓点X坐标（不设缓和曲线时，为直圆点X坐标）YJD? 输入交点的Y坐标YZH= 直缓点Y坐标（不设缓和曲线时，为直圆点Y坐标）F1? 输入小里程段导线的方位角XHZ= 缓直点X坐标（不设缓和曲线时，为圆直点X坐标）YHZ= 缓直点Y坐标（不设缓和曲线时，为圆直点Y坐标）XHY= 缓圆点X坐标不设缓和曲线时，不显示YHY= 缓圆点Y坐标XYH= 圆缓点X坐标YYH= 圆缓点Y坐标DQD? 输入待求点桩号FWJ= 待求点切线方位角X= 待求点中桩X坐标Y= 待求点中桩Y坐标ZB? 输入左幅宽度XZ= 待求点左幅边桩X坐标YZ= 待求点左幅边桩Y坐标YB? 输入右幅宽度XY= 待求点右幅边桩X坐标YY= 待求点右幅边桩Y坐标注意fx-4500P型计算器将L7行最后一句改成：Prog PQXZB；子程序中Prog ″PQX.F″改成Prog PQX.F，Prog ″PQX.BZ″改成Prog PQX.BZ【单交点基本型（对称、非对称）平曲线中桩、边桩坐标计算子程序】文件名：PQXZB（可单独执行，调用子程序PQX.F，PQX.BZ）L1 R:H″LS1″:H=0=>I=0:◣I″LS2″:Z:{XY}:X″XJD″:Y″YJD″:T″T1″:H=0=>U=T:◣U″T2″:B″ZH″:H=0=>C=B:◣C″HY″:E″YH″:H=0=>F=E:◣F″HZ″:{A}:A″F1″:Q=A+Z:Z>0=>S=1:≠=>S=-1:◣ J″XZH″=X+Tcos (A+180)◢K″YZH″=Y+Tsin (A+180)◢P″XHZ″=X+Ucos Q◢O″YHZ″=Y+Usin Q◢H≠0=>V=C-B-(C-B)^5÷40R2H2:W=(C-B)^3÷6RH-(C-B)^7÷336R^3H^3:M″XHY″=J+Vcos A+Wcos (A+90S)◢N″YHY″=K+Vsin A+Wsin (A+90S)◢V=F-E-(F-E)^5÷40R2I2:W=(F-E)^3÷6RI-(F-E)^7÷336R^3I^3:Z ″XYH″=P+Vcos (Q+180)+Wcos (Q+90S)◢Z ″YYH″=O+Vsin (Q+180)+Wsin (Q+90S)◢◣Lbl 0L2 {Z}:Z″DQD″L3 Z<B=>Goto 0:◣Z>F=>Goto 0:◣H=0=>G=A+180(Z-B)S÷πR:Prog ″PQX.F″:M=J:N=K:Goto 1:◣Z≤C=>G=A+90(Z-B)2S÷πRH:Prog ″PQX.F″:V=Z-B-(Z-B)^5÷40R2H2:W=(Z-B)^3÷6RH-(Z-B)^7÷336R^3H^3:X=J+Vcos A+Wcos (A+90S)◢Y=K+Vsin A+Wsin (A+90S)◢Prog ″PQX.BZ″:Goto 0:◣Lbl 1L4 Z≤E=>G=A+90(C-B)2S÷πRH+180(Z-C)S÷πR:Prog ″PQX.F″:V=90H÷πR:W=90(Z-C)÷πR:X=M+2Rsin Wcos (A+SV+SW)◢Y=N+2Rsin Wsin (A+SV+SW)◢Prog ″PQX.BZ″:Goto 0:◣Z>E=>G=Q-90(F-Z)2S÷πRI:Prog ″PQX.F″:V=F-Z-(F-Z)^5÷40R2I2:W=(F-Z)^3÷6RI-(F-Z)^7÷336R^3I^3:X=P+Vcos (Q+180)+Wcos (Q+90S)◢Y=O+Vsin (Q+180)+Wsin (Q+90S)◢Prog ″PQX.BZ″:◣Goto 0【平曲线切线方位角计算结果处理子程序】文件名：PQX.FL1 G>0=>G<360=>G″FWJ=″◢≠>G″FWJ″=G-360◢◣≠>G″FWJ″=G+360◢◣【平曲线边桩坐标计算子程序】文件名：PQX.BZL1 {L}:V=G-90:W″XZ″=X+L″ZB″cos V◢W″YZ″=Y+Lsin V◢{L}:V=G+90:W″XY″=X+L″YB″cos V◢W″YY″=Y+Lsin V◢【契形体体积计算程序】文件名：XXTTJL1 V=H÷6×(A″A1″×B″B1″+C″A2″×D″B2″+(A+C)(B+D)) 【竖曲线要素及高程计算程序】文件名：SQXJSL1 M″BZH″:G″BGC″:A″I1″:B″I2″:RL2 W=B-A◢L=RAbs W◢T=L÷2◢E=T2÷2R◢F″QZH″=M-T◢J″ZZH″=M+T◢Lbl 0L3 {Z}:Z″ZH″L4 Z≤F=>H=G-(M-Z)A◢Goto 0:◣Z≥J=>H=G+(Z-M)B◢Goto 0:◣Y=(Z-F)2÷2RL5 H=G-(M-Z)AL6 W≥0=>H=H+Y◢≠=>H=H-Y◢◣Goto 0 使用说明屏幕提示执行操作屏幕显示计算结果BZH? 输入变坡点桩号W= 前后纵坡坡度差BGC? 输入变坡顶点高程L= 竖曲线长度I1? 输入小里程段纵坡(前进方向上坡为+，下坡为-) T= 竖曲线切线长I2? 输入大里程段纵坡(前进方向上坡为+，下坡为-) E= 竖曲线外距R? 输入竖曲线半径QZH= 竖曲线起点桩号ZH? 输入计算点的桩号ZZH= 竖曲线终点桩号H= 计算点的中桩高程【无竖曲线段纵坡高程连续计算程序】文件名：ZPGC L1 F″ZH0″:G″GC0″:IL2 Lbl 0:{N}:N=0=>Goto 0:◣Z=FL3 Lbl 1L4 ″ZH=″:Z=Z+N:Pause 0L5 ″H=″:H=G+(Z-F)I:Pause 2L6 Goto 1使用说明屏幕提示执行操作屏幕显示计算结果ZH0? 输入已知桩号ZH= 计算点的桩号GC0? 输入已知桩号的高程I? 输入纵坡(前进方向下坡为负，上坡为正)H= 计算点的高程N? 输入桩号增量【余弦定理的坐标运用】文件名：YXDL-ZBL1 A″AX″:B″AY″:C″BX″:D″BY″:E″CX″:F″CY″L2 G″AB″=√((C-A)2+(D-B)2)◢H″BC″=√((E-C)2+(F-D)2)◢I″CA″=√((E-A)2+(F-B)2)◢J″<A″=cos-1 ((G2+I2-H2)÷2GI)◢J″360°-<A″=360-J◢J″<B″=cos-1 ((G2+H2-I2)÷2GH)◢J″360°-<B″=360-J◢J″<C″=cos-1 ((H2+I2-G2)÷2HI)◢J″360°-<C″=360-J使用说明屏幕提示执行操作屏幕显示计算结果AX? 输入A点坐标X AB= AB距离AY? 输入A点坐标Y BC= BC距离BX? 输入B点坐标X CA= CA距离BY? 输入B点坐标Y <A= <A角度CX? 输入C点坐标X <B= <B角度CY? 输入C点坐标Y <C= <C角度【施工坐标与大地坐标相互转换程序】文件名：ZBZH（调用子程序文件名：Z.FS，Z.ZH）L1 Lbl 0L2 {XYFG}:X″AX1″:Y″AY1″:U=F″OX1″:Z=G″OY1″L3 Prog ″Z.FS″L4 D=0=>Goto 0:◣E=AL5 Lbl 1L6 {XYHI}:X″AX2″:Y″AY2″:U=H″OX2″:Z=I″OY2″L7 Prog ″Z.FS″L8 D=0=>Goto 1:◣E≥A=>E″ANGLE1->2″=E-A◢≠=>E″ANGLE1->2″=E-A+360◢◣Prog ″Z.ZH″使用说明屏幕提示执行操作屏幕显示计算结果AX1? 输入A点在坐标系1中的坐标XANGLE1->2= 坐标系1到坐标系2的转角AY1? 输入A点在坐标系1中的坐标Y OX1? 输入O点在坐标系1中的坐标X OY1? 输入O点在坐标系1中的坐标Y AX2? 输入A点在坐标系2中的坐标X AY2? 输入A点在坐标系2中的坐标YOX2? 输入O点在坐标系2中的坐标XX2= 计算点在坐标系2中的坐标XOY2? 输入O点在坐标系2中的坐标YX1? 输入计算点在坐标系1中的坐标XY2= 计算点在坐标系2中的坐标YY1? 输入计算点在坐标系1中的坐标Y注意fx-4500P型计算器将L3、L7行改成：Prog Z.FSL8行最后一句改成：Prog Z.ZH【坐标反算子程序】文件名：Z.FSL1 B=X-U:C=Y-ZL2 D=√(B2+C2)L3 D=0=>Goto 0:◣C≥0=>A=cos-1(B÷D):≠=>A=360-cos-1(B÷D):◣Lbl 0【坐标转换子程序】文件名：Z.ZH（可单独执行，调用子程序：Z.FS）L1 U=F″OX1″:Z=G″0Y1″:H″0X2″:I″OY2″:E″ANGLE1->2″L2 Lbl 0L3 {XY}:X″X1″:Y″Y1″L4 Prog ″Z.FS″L5 D=0=>Goto 0:◣A=A-EL6 ″X2=″:X=H+Dcos A:Pause 0L7 Y″Y2″=I+Dsin A◢Goto 0使用说明屏幕提示执行操作屏幕显示计算结果OX1? 输入O点在坐标系1中的坐标XX2=计算点在坐标系2中的坐标XOY1? 输入O点在坐标系1中的坐标YOX2? 输入O点在坐标系2中的坐标XY2=计算点在坐标系2中的坐标YOY2? 输入O点在坐标系2中的坐标YANGLE1->2? 输入坐标系1到坐标系2的转角X1? 输入计算点在坐标系1中的坐标XY1? 输入计算点在坐标系1中的坐标Y注意fx-4500P型计算器将L4行改成：Prog Z.FS【仁宗海交通桥拱圈悬链线坐标计算程序】文件名：XLXL1 M=2.514:F=12.195:L=30.484L2 Lbl 0L3 {X}:J=30.484-XL4 P=-0.00011837194X2+0.0072169X◢A=Jln (M+√(M2-1))÷LL5 D=tan-1(FAsinh A÷J(M-1))◢Y″YZ″=F-F(cosh A-1)÷(M-1)+2363.821+P◢Z″YS″=Y+0.7÷cosD◢{H}:U=Z-H◢W″YX″=Y-0.7÷cosD◢Q=W-H◢Goto 0使用说明屏幕提示执行操作屏幕显示计算结果X? 输入拱圈X坐标(以拱圈左下角为坐标原点)P= 计算点的设计预拱度YZ= 中轴设计高程（已加了设计预拱度P）YS= 上缘设计高程（已加了设计预拱度P）H? 输入实测高程U=上缘设计高程与实测高程的高差：负数表明实测高程高于设计高程正数表明实测高程低于设计高程YX= 下缘设计高程（已加了设计预拱度P）Q=下缘设计高程与实测高程的高差：负数表明实测高程高于设计高程正数表明实测高程低于设计高程注意fx-4500P型计算器将L5行中sinh A改成(e A-e-A)÷2；将coshA改成(e A+e-A)÷2【高斯坐标转换程序】文件名：GSZHL1 A″H0″:Lbl 1:{BCH}:E=B″H1″+C″H2″:F″HD0″=(1+(A-.5E)÷6378245)H″HD″◢Goto 1【高斯坐标转换程序1】文件名：GSZH1L1 A″H0″:R=6378245:Lbl 1:{BCH}:E=B″H1″+C″H2″:F″HD0″=(R+A)H″HD″÷(R+.5E)◢Goto 1使用说明屏幕提示执行操作屏幕显示计算结果H0? 输入投影面高程HD0=高程投影后的平距H1?、H2? 输入该边两端点的高程HD? 输入平距【40米T型梁梁底模反拱度高程控制程序】文件名：FGD-Q（使用全站仪）L1 Lbl 0L2 {XQ}L3 P=.003X-.04+(20-X)2÷10000◢Z=P+2717.90358◢B=Z-Q◢Goto 0使用说明屏幕提示执行操作屏幕显示计算结果X? 输入X坐标（以0＃台梁端为(0,0)点，输入范围0～40）P= 计算点与0,0点的设计高差Z= 计算点的设计高程Q? 输入实测高程B= 设计高程－实测高程（正值表示该点偏低，负值表示该点偏高）注意fx-4500P型计算器将L3行中“÷”改为“/”，0＃台梁端高程为2717.90358【40米T型梁梁底模反拱度高程控制程序】文件名：FGD-S（使用水准仪）L1 Lbl 0L2 {XQ}L3 P=.003X-.04+(20-X)2÷10000◢Z=P+2717.90358◢B=Z-S+Q◢Goto 0使用说明屏幕提示执行操作屏幕显示计算结果X?输入X坐标（以0＃台梁端为(0,0)点，输入范围0～40）P= 计算点与0,0点的设计高差Z= 计算点的设计高程S?输入视线高程（单位：m，只输入一次）B=设计高程－实测高程（正值表示该点偏低，负值表示该点偏高）Q? 输入前视读数（单位：m）注意fx-4500P型计算器将L3行中“÷”改为“/”，0＃台梁端高程为2717.90358。
FX4800P测量程序集
CASIO fx 系列编程计算器测量程序集编写: 王劲松第一部分：CASIO fx-4800P 编程计算器测量程序集一．SHI JU CE TU （视距测图）1． 用途：该程序用于“经纬仪视距法碎部测量”时，计算测站点至碎部点间的平距及碎部点的高程。
2． 程序清单：C “H0”：I ：A “UP ”：B “DOWN ”：V “MIDDLE ”：L ：E=90-L ↙D=100Abs （A-B ）（cosE ）2◢ H=（C+I ）+DtanE-V3． 程序说明：H0——测站点高程；I ——测站仪器高；UP ——塔尺上丝读数；DOWN ——塔尺下丝读数；MIDDLE ——塔尺中丝读数；L ——经纬仪竖盘读数盘左读数；D ——测站点至碎部点间平距；H ——碎部点的高程。
5．该程序若在fx-3950中的程序清单为：？→M ：？→A ：？→B ：？→C ：？→D ：100（A-B ）（cos(90-D)）2→X ◢ M+Xtan(90-D)-C →Y说明：M ——将测站高程H0+仪器高I ；A ——上丝读数；B ——下丝读数；C ——中丝读数；D ——竖盘盘左读数L 。
X ——测站点至碎部点的平距；Y ——碎部点的高程。
二．FAN SUAN （反算）1．用途：该程序是“坐标反算公式”程序，用于根据两点1、2的平面坐标（X1，Y1）、（X2，Y2），计算两点间的距离D12及坐标方位角A12。
2． 程序清单：G “X1”：B “Y1”：C “X2”：H “Y2”：E=C-G ：F=H-B ↙Pol （E ，F ）：I “D12”◢ ∆+=⇒<3600J J J J “A12”3． 程序说明：X1，Y1——点1的平面坐标；X2，Y2——点2的平面坐标； D12——点1至点2间平距；A12——点1至点2坐标方位角。
由于编程时用到了极坐标函数POL （），故在程序运行中，会出现I=？XXX （数字）及J=？XXX （数字）这是函数POL计算出的两个量，不能改变，直接回车默认即可。
CASIO fx-4800P程序
1、方位角计算程序:【FWJ】 (2)2、转换坐标与高程程序:【0】 (2)3、坐标计算主程序:【XY】 (2)4、高程计算主程序:【H】 (5)5、坐标计算子程序:【XY-1】 (5)6、高程计算子程序:【H-1】 (6)7、水平距离归算到高斯平面计算程序:【JLGF】 (8)8、任意多边形面积计算程序:【MJ】 (8)9、四边形平差计算程序:【DDSBX】 (8)10、坐标换带计算程序:【ZBFD】 (11)11、两直线交点计算程序:【L-XY-L】 (12)12、高斯投影坐标计算程序:【GSBL-XY】 (12)13、竖曲线计算程序:【SQX】 (13)14、复曲线计算程序:【FQX】 (13)15、匝道计算子程序:【FQX】 (14)1、方位角计算程序:【FWJ】置镜点A＂X【Z】＂:B＂Y【Z】＂:后视点C＂X【H】＂:D＂Y【H】＂:前视点E＂X【Q】＂:F＂Y【Q】＂↵M=C-A:N=D-B:V=E-A:W=F-B:L=√（M2+N2）:S=√（V2+W2）:I=AbsTan-1（N÷M）:Q= AbsTan-1（W÷V）↵M＞0＝＞Goto 1:≠＝＞Goto 2↵Lbl 1↵N＞0＝＞O=I:Goto 3:≠＝＞O=360-I:Goto 3↵Lbl 2↵N＞0＝＞O=180-I:Goto 3:≠＝＞O=180+I:Goto 3↵Lbl 3↵V＞0＝＞Goto 4:≠＝＞Goto 5↵Lbl 4↵W＞0＝＞H=Q:Goto 6:≠＝＞H=360-Q:Goto 6↵Lbl 5↵W＞0＝＞H=180-Q:Goto 6:≠＝＞H=180+Q:Goto 6↵Lbl 6↵O＞H＝＞X=O-H:Goto 9:≠＝＞X=O-H+360:Goto 9↵Lbl 9↵＂F【Z-H】=＂:O ◢＂L【Z-H】=＂:L ◢＂F【Z-Q】=＂:H ◢＂L【Z-Q】=＂:S ◢＂f【H-Z-Q】=＂:360-X（左角）◢＂f【Q-Z-H】=＂: X（右角）计算示例:2、转换坐标与高程程序:【0】计算选择:1.坐标计算;2.高程计算。
坐标正算与反算、施工与大地坐标换算程序(4800P)
Goto 1
算例：1点，X=10.000，Y=20.000
2点，X=500.000，Y=600.000
3点，X=700.000，Y=400.000
专业的3S站
DB子程序：V=V+1:B=Z[V]+T:Prog“A”:D=Z[V+N]:Prog“XY”
DEG子程序：P=Frac B*100:B=INTB+INTP/60+Frac P/36
三、 支导线的计算
支导线就是从一个已知点出发，既不符合到另一个已知点也不回到原来的点上。支导线没有检核条件，不易发现错误，过不宜采用。
X？ 输入5点X坐标
Y？ 输入5点Y坐标
S=218500 （1.2.3.4.5点围成的面积）
X？ 输入6点X坐标
Y？ 输入6点Y坐标
S=230500 （1.2.3.4.5.6点围成的面积）
.........
.........
Lbl 3：｛XY｝：X“X=”：Y“Y=”：A=（X－C）cos R+(Y－D) sin R:“A=” ：Pause 0:B=（Y－D）cos R－(X－C) sin R:B“B=” ▲ Goto 0 △
Lbl 4：｛AB｝：A“A=”：B“B=”：X=C+A cos R－B sin R：“X=”：Pause 0:Y=D+A sin R+B cos R：Y“Y=” ▲ Goto 0 △
坐标正算与反算、施工与大地坐标换算程序（4800P）
Deg：Fix4：J=0：C“X=”：D“Y=”：R“R=“：
CASIO-fx4800P型计算器连续计算线路高程、坐标Word文档
CASIO-fx4800P型计算器连续计算线路高程、坐标及放样程序【简述】为了适应现场快速、准确、灵活的放样要求，作者根据线路施工测量特点，运用计算器的编程功能，使整条线路的施工测量计算数据能连续快速的计算，只要把整条线路各个交点处的计算要素输入子程序中，你便拥有整条线路的“数据库”。
在忙碌的施工现场，你不必携带大量的施工图和资料来查找平面、纵断面的设计要素，更不必担心设计要素的输入错误，从而使得计算和放样轻松、方便、快捷。
【计算范围】坐标部分包括：各等级公路和高速路的直线段、圆曲线段、加对称缓和曲线的圆曲线段中、边桩坐标及放样计算。
高程部分包括：各等级公路和高速路的直线段、竖曲线段、缓和超高段中、边桩高程及放样计算。
一、坐标及放样程序【起算数据】：交点里程桩号、交点坐标、前直线方位角、交点转角、圆曲线半径、缓和曲线长（一）主程序：XYFY2Fix4M"Xc=":P"Yc=" ＇如须放样输入测站坐标值,否则输入0Lbl0{HS}:H"K0+0":S"B=(-B,0,B)":Prog"A" ＇变量输入和声明S=0=>E=0:≠>E=90 ＇求边桩坐标时的偏角值判定L=πRN÷180+V ＇含有缓曲线的曲线总长T=(V÷2-V^3÷(240R^2))+(R+(V^2÷(24R)-V^4÷(2688R^3)))tan(N÷2)＇切线长A=Q-T:B=A+V:D=A+L:C=D-V ＇A-ZH(ZY)点桩号,B-HY点桩号,C-YH点桩号,D-HZ(YZ)点桩号Rec(T,F+180) ＇求ZH(ZY)点的坐标增量Z[1]=W+I:Z[2]=K+J ＇ZH(ZY)点坐标值Rec(T,F+GN) ＇求HZ(YZ)点的坐标增量Z[3]=W+I:Z[4]=K+J ＇HZ(YZ)点的坐标值Lbl1U=S÷(Abs(S+10E-19)):S=AbsS ＇左右偏角每件判定H<A=>Goto2: ≠> ＇第一直线上坐标计算判定H<B=>Goto3: ≠> ＇第一缓和曲线上坐标计算判定H<C=>Goto4: ≠> ＇净圆曲线上坐标计算判定H<D=>Goto5: ≠> Goto6⊿⊿⊿⊿＇第二缓和曲线和第二直线上坐标计算判定Lbl2Rec(Q-H,F+180) ＇第一直线上中桩坐标增量计算X=W+I:Y=K+J ＇中桩坐标值计算Rec(S,F+180-(180-E)U) ＇第一直线上与中桩同一断面的边桩增量计算X=X+I:Y=Y+J ＇边桩坐标值计算S=SU ＇边长条件判定Goto8Lbl3＇进入第一缓曲线段计算Z=H-A ＇待求点至ZH点的距离O=90Z^2÷(πRV) ＇所求点缓曲线对应的圆心角X=Z-Z^5÷(40R^2V^2)+Z^9÷(3456R^4V^4) ＇缓曲线上ZH点至待求的支距X坐标Z=Z^3÷(6RV)-Z^7÷(336R^3V^3)+Z^11÷(42240R^5V^5) ＇缓曲线上ZH点至待求点的支距Y坐标Lbl7Rec(X,F)X=Z[1]+I:Y=Z[2]+JRec(Z,F+90G) ＇把中桩支距坐标转换成测量统一坐标X=X+I:Y=Y+JRec(S,F+O G+EU) ＇曲线边桩坐标增量计算X=X+I:Y=Y+J ＇坐标计算结果S=SU ＇边长条件Goto8Lbl4＇进入主圆计算部分Z=H-A-V÷2 ＇曲线上待求点到HY点的里程O=180Z÷(Rπ) ＇待求点方位角所对应的圆心角X=Rsin O+(V÷2-V^3÷(240R^2))Z=R(1-cos O)+(V^2÷(24R)-V^4÷(2688R^3))＇主圆部分待求点支距坐标XY Goto7Lbl5＇进入第二缓和曲线计算Z=D-H ＇曲线上待求点到HZ点的里程O=90Z^2÷(πRV) ＇待求点方位角所对应的圆心角X=Z-Z^5÷(40Z R^2V^2)+^9÷(3456R^4V^4)Z=Z^3÷(6RV)-Z^7÷(336R^3V^3)+Z^11÷(42240R^5V^5) ＇待求点支距坐标XY Rec(X,F+GN+180)X=Z[3]+I:Y=Z[4]+J ＇把中桩支距坐标转换成测量统一坐标Rec(Z,F+GN+180-90G)X=X+I:Y=Y+JRec(S,F+GN+180-O G-(180-E)U) ＇曲线边桩坐标增量计算X=X+I:Y=Y+J ＇坐标计算结果S=SUGoto8Lbl6＇进入第二直线段计算部分Rec(H-D+T,F+GN) ＇中桩坐标增量计算X=W+I:Y=K+J ＇中桩坐标计算结果Rec(S,F+GN+EU) ＇边桩坐标增量计算X=X+I:Y=Y+J ＇边桩坐标计算结果S=SUGoto8Lbl8X=X◢＇显示坐标计算结果Y=Y◢M≠0=>Goto9: ≠>Goto0 ⊿Lbl9 ＇进入放样部分计算Pol(X-M,Y-P) ＇增量计算J<0=>J=J+360: ≠>J=J ⊿J"A"=Int J+0.01Int(60Frac J)+0.006Frac(60Frac J) ◢＇方位角计算结果（此显示值为度分秒格式）I"L"=I◢＇极距计算结果Goto0(二)子程序:AH≤ZY2(ZH2)=>Q=JD1:W=X1:K=Y1:R=R1:F=F0:N=N1:V=Ls1:G=1(-1):=>H≤ZY3(ZH3)=>Q=JD2:W=X2:K=Y2:R=R2:F=F1:N=N2:V=Ls2:G=1(-1):=>H≤ZY4(ZH4)=>Q=JD3:W=X3:K=Y3:R=R3:F=F2:N=N3:V=Ls3:G=1(-1):=>H≤ZY5(ZH5)=>Q=JD4:W=X4:K=Y4:R=R4:F=F3:N=N4:V=Ls4:G=1(-1):=>……………………………H≤ZY n+1(ZH n+1)=>Q=JD n:W=X n:K=Y n:R=R n:F=F n-1:N=N n:V=Ls n:G=1(-1)注:H—待求点桩号ZY n+1(ZH n+1)—后曲线起点桩号 Q—本曲线交点桩号 W—本曲线交点X坐标K—本曲线交点Y坐标 R—本曲线内圆曲线半径 F—前直线方位角 N—本交点转角V—本曲线内缓和曲线长 G—线路转向符,左转角-1,右转角为+1二、高程及放样程序【起算数据】：变坡点桩号A、变坡点高程B、前后纵坡C D、竖曲半径R、横坡度F、最大超高横坡度G、ZH点桩号、YH点桩号、HZ点桩号、施工层厚度N、视线高W （一）主程序：GCFYW"SXG="：N＇常量输入（视线高程、施工层厚度）Lbl 1｛LО｝:L"K0+0"：О"B=" :Prog"B"＇变量输入及声明K=C-D ＇坡度差及竖曲线偏角K>0=>J=-1:≠>J=1⊿＇凹凸竖曲线判断条件T=Abs(RK)÷2 ＇竖曲线切线长,如需显示则在后面加上◢H=A-T ＇竖曲线起点ZY桩号,如需显示则在后面加上◢X=A+T ＇竖曲线终点YZ桩号,如需显示则在后面加上◢M=L-A ＇待求点至变坡点的距离L≤H=>P=B+MC-N-W: ≠> ＇前纵坡直线段中桩高程计算L≤A=>P=B+MC+J(T-AbsM)2÷(2R)-N-W: ≠>＇前纵坡竖曲线中桩高程计算L≤X=>P=B+MD+J(T-AbsM)2÷(2R)-N-W: ≠>＇后纵坡竖曲线中桩高程计算L>X=>P=B+MD-N-W ⊿⊿⊿⊿⊿＇后纵坡直线段中桩高程计算E=0=>Goto 2 ＇直线部分中、边桩计算判断条件≠>L≤E=> Goto 2 ＇直线部分中、边桩计算判断条件≠>L>V=> Goto 2 ＇直线部分中、边桩计算判断条件≠>L>E=> Goto 3 ＇前缓曲超高部分中、边桩计算判断条件≠> L>U=> Goto 4⊿⊿⊿⊿＇后缓曲超高部分中、边桩计算判断条件Lbl 2O=0=>Z“H”=AbsP◢＇显示直线部分中、边桩高程计算结果≠>Y“HB”=Abs(P-OF)◢不需放样显示高程,否则显示放样水准尺计数⊿Goto 1Lbl 3S=（L-E）（F+G）÷（V-U）-F＇前缓曲线超高横坡度计算S≤F=> Goto 5：≠>S≤G=> Goto 6 ＇超高起点(终点)至零界面、零界面至全超高及全超高≠> Goto 7 ⊿⊿部分的计算判断条件Lbl 4S=（V-L）（F+G）÷（V-U）-F＇后缓曲线超高横坡度计算S≤F=> Goto 5：≠>S≤G=> Goto 6 ＇超高起点(终点)至超高变坡零界面、零界面至全≠> Goto 7 ⊿⊿超高、全超高部分计算判断条件Lbl 5O=0=>Z“H1”=AbsP◢＇超高起点(终点)至超高变坡零界面中、边桩高程≠>Y“H-D1”=Abs(P-OF)◢计算结果显示Q“H-G1”=Abs(P+OS)◢⊿Goto 1Lbl 6O=0=>Z“H2”=AbsP◢＇超高变坡零界面至全超高中、边桩高程≠>Y“H-D2”=Abs(P-OS)◢计算结果显示Q“H-G2”=Abs(P+OS)◢⊿Goto 1Lbl 7O=0=>Z“H3”=AbsP◢＇全超高中、边桩高程计算结果显示≠>Y“H-D3”=Abs(P-OG)◢Q“H-G3”=Abs(P+OG)◢⊿Goto 1注：W“SXG”-视线高N-施工层厚度L"K0+0"-待求点桩号О"B="-中桩至边桩距离（二）子程序：BL≤ZY2=>A=JD1:B=H1:C=I0:D=I1: R=R1:F=I h:E=ZH1:U=YH1:V=HZ1:G=I max1 :≠>L≤ZY3=>A=JD2:B=H2:C=I1:D=I2: R=R2:F=I h:E=ZH2:U=YH2:V=HZ2:G=I max2 :≠>L≤ZY4=>A=JD3:B=H3:C=I2:D=I3: R=R3:F=I h:E=ZH3:U=YH3:V=HZ3:G=I max3 :≠>L≤ZY4=>A=JD3:B=H3:C=I2:D=I3: R=R3:F=I h:E=ZH3:U=YH3:V=HZ3:G=I max3 :≠> ……………………L≤ZY n+1=>A=JD n:B=H n:C=I n-1:D=I n: R=R n:F=I h:E=ZH n:U=YH n:V=HZ n:G=I maxn注:ZY N+1:-下一个竖曲线起点桩号（如：变坡点2的计算范围为前竖曲线终点YZ1至后竖曲线起点ZY3，则在程序中“L≤ZY3”就应输入ZY3的里程桩。
CASIO4800P计算器坐标计算程序
坐标放样程序Defm 60：A”JX”:B”JY”:Q:R:L”L1”:M”L2”:T”T1”: C”T2”:E”ZH”:D=E+L:O”HZ”:I=O-M:F:U :Fix3：回车Lbi 0:{K}:K≤E=〉G0t0 1：≠=〉G0t0 2⊿回车Lbi1：{VJZ}:Z[1]=A-T CosF:Z[2]=B-TSinF: Z[3]=E-K: Z[4]=F+VJ: Z[20]= Z[1]- Z[3] CosF:Z[21]=Z[2]-Z[3]SinF：X=Z[20]+ZCos （Z[4]）：Y=Z[21] +ZSin（Z[4]）：X“X=”▲Y“Y=”▲{N}：N＜1=〉Prog“FMH”:≠=> Goto 0回车Lbi 2：K≤D=> Goto 3: ≠=> Goto 4⊿回车Lbi3:{VJZ}:Z[5]=A-TCosF:Z[6]=B-TSinF: H=K-E: Z[7]=H-（(H＾5)÷（40R＾2L＾2)）-（(H＾13)÷（599040R＾6L＾6)）+（(H＾9)÷（3456R＾4L＾4)）:W=F+U(90H＾2÷（RL ∏）)+VJ:Z[22]=Z[5]+Z[7] ÷Cos(（30H＾2）÷（RL∏)） *Cos(F+U(30H^2)÷(RL∏)): Z[23]= Z[6]+ Z[7] ÷Cos((30H^2)÷(RL∏)) *Sin(F+U(30H^2)÷(RL∏)):回车X= Z[22]+Z CosW:Y=Z[23]+Z SinW: 回车X“X=”▲Y“Y=”▲Prog“AA”:回车{N}：N＜1=〉Prog“FMH”:≠=> Goto 0回车Lbi 4:K≤I=> Goto 5:≠=> Goto 6⊿回车Lbi5:{VJZ}:Z[8]=A+TCos(F+180):Z[9]=B+ TSin(F+180):Z[10]=L-(L^3)÷(40R^2)-(L^7)÷(599040R^6)+(L^5)÷(3456R^4): Z[11]=(30L)÷(R∏): Z[12]=F+U(30L)÷(R∏):P=K-D: Z[13]=F+90U（L+P）÷(R∏): Z[14]=(90P)÷(R∏):W=F+U(90L)÷(R∏)+(180UP)÷(R ∏)+VJ: 回车Z[24]= Z[8]+ Z[10]÷Cos(Z[11]) *Cos(Z[12])+2RSin(Z[14]) *Cos(Z[13]): Z[25]= Z[9]+ Z[10] ÷Cos(Z[11]) *Sin(Z[12])+2RSin(Z[14]) Sin(Z[13]):X=Z[24]+Z CosW:Y= Z[25]+Z SinW: X“X=”▲Y“Y=”▲Prog“BB”:{N}：回车N〈1=>Prog“FMH”:≠=> Goto 0⊿回车Lbi 6:K≤O=> Goto 7: ≠=> Goto 8⊿回车Lbi 7:{VJZ}:P=F+Q:Z[15]=A+C CosP: Z[16]=B+C SinP:H=O-K: Z[17]=H-(H∧5)÷(40R∧2M∧2)-(H∧13)÷(599040R∧6M∧6)+(H∧9)÷(3456R∧4M∧4): Z[18]=(30H∧2) ÷(RM ∏): Z[19]=P+180-U(30H^2)÷(RM ∏):W=P+180-U((90H∧2) ÷(RM∏))-VJ:Z[26]= Z[15]+ Z[17] ÷Cos (Z[18]) *Cos(Z[19]): Z[27]= Z[16]+ Z[17] ÷Cos( Z[18])* Sin (Z[19]):X= Z[26]+ Z CosW:Y= Z[27]+Z SinW:X“X=”▲Y“Y=”▲Prog“CC”:回车{N}：N〈1=>Prog“FMH”:≠=> Goto 0回车Lbi 8:K≤O=> Goto7: ≠=> Goto9⊿回车Lbi9:{VJZ}:P=F+Q:W=P+VJ:Z[28]=A+(C+K-O) CosP: Z[29]=B+(C+K-O) SinP:X= Z[28]+Z CosW:Y= Z[29]+Z SinW:X“X=”▲Y“X=”▲Prog“DD”:{N}：N〈1=>Prog“FM H”:≠=> Goto 0⊿2、文件名“ FMH”Lbi1:{S}:S<0=> Goto2: ≠=> Goto3⊿回车Lbi2{J}:Z[30]=X- Z[20]: Z[31]=Y- Z[21]: Z[32]= Z[31] ÷Z: Z[33]= Z[30] ÷Z:V= Z[20]+JZ[33]-SZ[32]:K= Z[21]+ JZ[32]+S Z[33]:V“X=”▲K“Y=”▲Goto13、文件名“AA”Z[20]= Z[22]: Z[21]= Z[23]4、文件名“BB”Z[20]= Z[24]: Z[21]= Z[25]5、文件名“CC”Z[20]= Z[26]: Z[21]= Z[27]6、文件名“DD”Z[20]= Z[28]: Z[21]= Z[29]道路前进方向操作：A:B:坐标:Q转角（左负，右正）：R：半径：L1：第一缓和曲线长度（L）：L2:第二缓和曲线长度(M):T1:第一起线长度（T）:T2: 第二起线长度(C):ZH:直缓点里程(E):HZ:缓直点里程(O):F：方位角：U：转向系数（左转负，U=-1,右正，U=1）:K:待求点里程：V:方向系数（左取－1，右取1）：J：边桩与道路中线的夹角:Z:边距：在左边上时S取正J取0左边下时S取负J 为距中宽度，在右边上时S取负J取距宽度中，在右边下时S取正J取0。
fx_4800P型计算器坐标变换程序
f x 24800P 型计算器坐标变换程序杜荣学1,段泓娟2,段计平3(1.中国冶金地质勘查工程总局第三地质勘查院,太原030002;2.中国冶金地质勘查工程总局三局312队,山西忻州034000;3.中国冶金地质勘查工程总局三局深蓝公司,太原030002)摘 要: 用CASIO f x 24800P 型可编程计算器进行经纬度、3°带和6°带坐标转换应用程序实现了用使用率较高的普通编程计算器进行一般精度要求的高斯投影面坐标转换的计算,文章介绍了该应用程序及其计算精度和使用方法。
关键词: fx 24800P 型可编程计算器;坐标转换应用程序;高斯投影面;精度可达百分位(厘米)中图分类号: TP317;P282.2 文献标识码: A 文章编号:100121412(2005)S020194202 经纬度、3°带坐标及6°带坐标三者间的相互转换在地质工作、测量工作及其它作业中应用愈来愈广泛,电脑、手持GPS 机中坐标转换及相关软件的开发也较普遍,然而对于目前工作中利用率很高的普通可编程计算器,该类应用程序的开发仍然较少,近年来笔者就CASIO f x 24800P 型可编程计算器编制了一套坐标转换应用程序,经工作中应用,操作简单,显示明了,可适用于一般精度的坐标转换工作,具有一定的实用价值。
1 程序内容Lb11↓Fix 2:{Z }Z =1]“-BEGIN -XY-BL ”:G oto 2:≠]△↓Z =0]“-BEGIN -X Y -X Y ”:G oto 2:≠]△↓Z =-1]“-BEGIN -BL -X Y ”:G oto 4:≠]G oto 1↓Lb12↓{X YA }↓Fix 2:X YA “L0”↓G =180×0.157046064172×10-6×X ÷π:H =cos G :I =sin G :J =G +HI (5.051773759×10-3+I 2(-2.9837302×10-5+I 2×2.38189×10-7))×180÷π:K =tan J :L =cos J :N =6.73852541468×10-3×L 2+1:O =6399698.9018÷N 0.5:P =(Y -500000)÷O:Q =P 2:R =Q ÷30×(61+(90+45K 2)K 2):↓B =J -QN K ÷2×(1-Q ÷12×(5+(N -1)+(3-9(N -1))K 2-R ))×180÷π:T =1-Q ÷6×(N +2K 2-Q ÷20×(5+(6+8K 2)(N -1)+(28+24K 2)K 2)):C =T P ÷L ×180÷π+A :↓Lb13↓Z =0]G oto 6:≠]Δ↓“------”Fix 6:V “L ”=Int C +Int (Frac C ×60)÷100+Frac (Frac C ×60)×60÷10000▲W “B ”=Int B +Int (Frac B ×60)÷100+Frac (Frac B ×60)×60÷10000▲“------”G oto 2Lb14↓{D }↓Z =-1]{US}:≠]Δ↓Fix 4:S “L ”U “B ”:↓C =Int S +Int (Frac S ×100)÷60+Frac (Frac S ×100)÷36:B =Int U +Int (Frac U ×100)÷60+Frac (Frac U ×100)÷36↓L b16↓{D }↓Fix 2:D “L 1”↓E =tan B :F =cos B :G =sin B :H =6.73852541468×10-3×F 2:I =1+H :J =6399698.9018÷I 0.5:收稿日期:2004203220作者简介:杜荣学(19702),男,山西襄汾人,工程师,1991年毕业于沈阳黄金学院,现从事地质勘查工作。
CASIO_fx-4800P常用测量程序
一、QXFY 辛甫森公式放样程序1. U“X0”：V“Y0”2. A“XA”：B“Y A”：C“CA”：D“1÷RA”：E“1÷RB”：F“KA”：G“KB”3. Lb1 1：{HLR}：H“KI”：L“JJ”：R“DZ+Y-”4. H>G Goto 15. H<F Goto 16. P=(E-D) ÷Abs(G-F)：Q=Abs(H –F)：Z=P×Q：T=D+Z7. W=C+（Z+2D）Q×90÷π8. M=C+(Z÷4+2D)Q×22.5÷π：N=C+（3Z÷4+2D）Q×67.5÷π：K=C+(Z÷2+2D)Q×45÷π9. X=A+Q（cosC+4(cosM+cosN)+2cosK+cosW）÷12+Rcos(W-L)10. Y=B+Q(sinC+4(sinM+sinN)+2sinK+sinW) ÷12+ Rsin (W-L)11. Pol((X-U)，（Y-V）)J≤0J=J+360J“FWJ”12. I“S”13. Goto 1程序显示说明：须输入参数：X0:置镜点X (对应变量U)Y0:置镜点Y (对应变量V)XA:曲线起算点X (对应变量A)Y A:曲线起算点Y (对应变量B)CA：曲线起算点切线方位角(对应变量C) 1÷RA：1÷半径，即起算点曲率，右偏为正左偏为负(对应变量D)1÷RB：1÷半径，即曲终点曲率，右偏为正左偏为负(对应变量E)KA:起算点里程(对应变量F)KB:曲终点里程(对应变量G)KI:待求点里程(对应变量H)JJ:夹角(与前进方向切线方位角夹角，向右为正)(对应变量L)D“Z+Y－”：偏距，左正右负(对应变量R)计算结果: W: 待求点切线方位角(对应变量W) X:计算点x (对应变量X)Y:计算点y (对应变量Y)FWJ:放样方位角(对应变量J)S:放样距离(对应变量I)K BX 0Y 0本程序依据复化辛甫森公式计算原理改进而成，特点是把曲线按曲率变化点分成若干计算单元单独计算，不论直线（曲率为0）、圆曲线（曲率为1/R ）、卵形曲线，只要按曲率变化点分解弄清变化点曲率半径，右偏左偏曲线由曲率的正负号分别，既可求得该计算单元内任意里程中桩及斜交或正交的边桩。
CASIO fx-4800P程序编制(里程中边桩计算坐标及反算)
一、源程序(SZ-XY-SZ)"1.SZ => XY"："2. XY => SZ"：NU"X0"V"Y0"O"S0"G"F0"Q：Q=0=>R=1E-45：≠>R△J=0：M=AbsQ：K=1-M：N=1=>Goto 1：≠>Goto 2△←┙Lbl 1：{SZ}：SZ：L=S-O：B=90QL÷π÷R：F=G+B：E=F+B+90：C=LK+2MRsinQB：X"XS" = U+Rec(C,F)+ZcosE◢Y"YS"=V+J+ZsinE◢Goto 1←┙Lbl 2：W=G+90Q：E=U+Rec(R,W)：F=V+J←┙Lbl 3：{XY}：XY：D=Pol(X-E,Y-F):T=J：J<0=>T=J+360△B=T-W-180M：P=Abs B：S"S"=O+ DKcosP+ MπRP÷180◢Z"Z"=DKsinB+MQ(R-D)◢Goto 3←┙二、使用说明1.由里程、中边桩计算坐标启动程序后按以下提示输入数据：N ？1 （N=1，由里程中边桩计算坐标；N≠1，由坐标反算里程及中边桩)X0 ？142.353 (直线或圆曲线起点的X坐标）Y0 ？368.166 (直线或圆曲线起点的Y坐标）S0 ？155.252 (直线或圆曲线起点的里程）F0 ？60°54′03″(直线或圆曲线起点的切线方位角）Q ？-1 (当线元为直线时Q=0；当线元为圆曲线且左偏时Q=-1；当线元为圆曲线且右偏时Q=1)R ？100 (圆曲线的半径)S ？194.016 (中线里程)Z ？-5.251 (中线左侧5.251m的边桩；在中线上Z=0，在左侧<0，在右侧Z>0)显示计算结果XS=172.7467767YS=395.23168072.由坐标反算里程、中边桩启动程序后按以下提示输入数据：N ？2 （N=1，由里程中边桩计算坐标；N≠1，由坐标反算里程及中边桩)X0 ？142.353 (直线或圆曲线起点的X坐标）Y0 ？368.166 (直线或圆曲线起点的Y坐标）S0 ？155.252 (直线或圆曲线起点的里程）F0 ？60°54′03″(直线或圆曲线起点的切线方位角）Q ？-1 (当线元为直线时Q=0；当线元为圆曲线且左偏时Q=-1；当线元为圆曲线且右偏时Q=1)R ？100 (圆曲线的半径)X ？172.7467767 (待求里程及中边桩的点的X坐标)Y ？395.2316807 (待求里程及中边桩的点的Y坐标)显示计算结果S=194.016Z=-5.251三、注意事项1.本程序主要是针对市政道路工程中的直线和圆曲线而编制的，未考虑缓和曲线等曲线的计算；2.程序中“R=1E-45”的意思为10的-45次方，输入程序时切莫输错。
卡西欧4800P计算器大地坐标放样程序
计算点K
LX
置仪点
置方位角 预
后视点
eh曲线外距p输入1时计算测点的左边点输入0时计算中点
主程序 WTO
(卡西欧 4800P 计算器用坐标放样程序)
L“L1”E“L2”A“ZH”：Z：F：U“JX”V“JY”M“X0”N“Y0”：J=180/πR：B“HZ”=A+Z/J+L/2+E/2▲ W“-1，1”：Prog“K”↙ LIb0：｛K｝：｛H｝：｛P｝：K≤A ＞Goto1△K≤A+L ＞Goto2△K＜B-E ＞Goto3△ K≤B ＞Goto4△Goto5△↙ LIb1：S=K-A：X=U-TcosF+ScosF+Hcos（F+90P）▲Y=V-TsinF+SsinF+Hsin（F+90P）▲Goto6↙ LIb2：S=K-A：Prog“A”：Goto6↙ LIb3：S=K-A-L：Prog“B”：Goto6↙ LIb4：S=B-K：Prog“C”：Goto6↙ LIb5：S=K-B：X=U+(I+S)cos(F+WZ)+Hcos(F+WZ+90P)▲ Y=V+(I+S)sin(F+WZ)+Hsin(F+WZ+90P)▲Goto6↙ LIb6：X≠M ＞Prog“G”↙ Goto0 子程序G：
L1－第一缓和曲线总长；L2－第二缓和曲线总长；Z－曲线总转角；F－ZH点至JD点的方位角：
-1 2 缓点桩号；HZ－缓直点桩号；JX－交点的横坐标(N)；JY－交点的纵坐标(E)； X0－置仪点的横坐标(N)；ZY－置仪点的纵坐标(E)；R－曲线半径：HZ－缓直点桩号； “-1，1?"—曲线左转时输如“-1”，曲线右转时输入“+1”；Eh——曲线外距 K－测点的桩号；H－测点边桩与中线之间的距离； P－输入“-1”时计算测点的左边点，输入“0”时计算中点；输入“+1”时计算右边点 X－测点的点的横坐标(N)；Y－测点的纵坐标(E)； LX－测距；C－置仪点与测点间的方位角(即拨角)。 注：应用此程序时须预置置仪点至后视点间的方位角。
fx4800编程计算器程序二
6. 圆曲线段 文件名：YQX (COMP) 程式： L”K0”:O”X0”:P”Y0”:W”ALF”: Lbl 1:{K}: J=(-1)^N (K-L)÷R×180÷π:D=2Rsin((-1)^N J÷2): X”XZ”=O+Dcos(W+J÷2) ◢ Y”YZ”=P+Dsin(W+J÷2) ◢ {B}:S”XL”=X-B”LL”cos(W+J+Q) ◢ T”YL”=Y-Bsin(W+J+Q) ◢ {C}:U”XR”=X+C”LR”cos(W+J+Q) ◢ V”YR”=Y+Csin(W+J+Q) ◢ Goto 1 注：若没有直接 HY 点方位角，则 ALFHY=ALFZH±Ls/2/R×180/π,(左转-，右转 +)。
3.缓和曲线曲线要素： 程序步骤： B=0°1718.87′﹡L/R◢ X=L-L＾3/40/ R2◢ Y=L2/6/ R◢ P=Y-R(1-cos B◢ Q= X-Rsin B◢ T=(R+P)tan(A/2)+Q◢ E=(R+P)(cos(A/2))-1-R◢ Z= R(A-2B)л/180+2L◢ J=2T-Z◢ D=X-Y/tan B 操作过程： HQQXYS→EXE→输入 L 值（即缓和曲线总长）→EXE→输入 R 值（即圆曲线半径） →EXE→得β角度→EXE→得 Xh 值→EXE→得 Yh 值→EXE→得 P 值→EXE→得 Q 值 →EXE→输入 A 角（例 125°31′23.25″）→EXE→得 T 值→EXE→得 E 值→EXE →得 Z 值→EXE→得 J 值→EXE→得 D 值 注：此程序可循环计算。
CASIO-fx4800P型计算器连续计算线路高程、坐标
CASIO-fx4800P型计算器连续计算线路高程、坐标及放样程序【简述】为了适应现场快速、准确、灵活的放样要求，作者根据线路施工测量特点，运用计算器的编程功能，使整条线路的施工测量计算数据能连续快速的计算，只要把整条线路各个交点处的计算要素输入子程序中，你便拥有整条线路的“数据库”。
在忙碌的施工现场，你不必携带大量的施工图和资料来查找平面、纵断面的设计要素，更不必担心设计要素的输入错误，从而使得计算和放样轻松、方便、快捷。
【计算范围】坐标部分包括：各等级公路和高速路的直线段、圆曲线段、加对称缓和曲线的圆曲线段中、边桩坐标及放样计算。
高程部分包括：各等级公路和高速路的直线段、竖曲线段、缓和超高段中、边桩高程及放样计算。
一、坐标及放样程序【起算数据】：交点里程桩号、交点坐标、前直线方位角、交点转角、圆曲线半径、缓和曲线长（一）主程序：XYFY2Fix4M"Xc=":P"Yc=" ＇如须放样输入测站坐标值,否则输入0Lbl0{HS}:H"K0+0":S"B=(-B,0,B)":Prog"A" ＇变量输入和声明S=0=>E=0:≠>E=90 ＇求边桩坐标时的偏角值判定L=πRN÷180+V ＇含有缓曲线的曲线总长T=(V÷2-V^3÷(240R^2))+(R+(V^2÷(24R)-V^4÷(2688R^3)))tan(N÷2) ＇切线长A=Q-T:B=A+V:D=A+L:C=D-V ＇A-ZH(ZY)点桩号,B-HY点桩号,C-YH点桩号,D-HZ(YZ)点桩号Rec(T,F+180) ＇求ZH(ZY)点的坐标增量Z[1]=W+I:Z[2]=K+J ＇ZH(ZY)点坐标值Rec(T,F+GN) ＇求HZ(YZ)点的坐标增量Z[3]=W+I:Z[4]=K+J ＇HZ(YZ)点的坐标值Lbl1U=S÷(Abs(S+10E-19)):S=AbsS ＇左右偏角每件判定H<A=>Goto2: ≠> ＇第一直线上坐标计算判定H<B=>Goto3: ≠> ＇第一缓和曲线上坐标计算判定H<C=>Goto4: ≠> ＇净圆曲线上坐标计算判定H<D=>Goto5: ≠> Goto6⊿⊿⊿⊿ ＇第二缓和曲线和第二直线上坐标计算判定Lbl2Rec(Q-H,F+180) ＇第一直线上中桩坐标增量计算X=W+I:Y=K+J ＇中桩坐标值计算Rec(S,F+180-(180-E)U) ＇第一直线上与中桩同一断面的边桩增量计算X=X+I:Y=Y+J ＇边桩坐标值计算S=SU ＇边长条件判定Goto8Lbl3 ＇进入第一缓曲线段计算Z=H-A ＇待求点至ZH点的距离O=90Z^2÷(πRV) ＇所求点缓曲线对应的圆心角X=Z-Z^5÷(40R^2V^2)+Z^9÷(3456R^4V^4) ＇缓曲线上ZH点至待求的支距X坐标Z=Z^3÷(6RV)-Z^7÷(336R^3V^3)+Z^11÷(42240R^5V^5) ＇缓曲线上ZH点至待求点的支距Y坐标Lbl7Rec(X,F)X=Z[1]+I:Y=Z[2]+JRec(Z,F+90G) ＇把中桩支距坐标转换成测量统一坐标X=X+I:Y=Y+JRec(S,F+O G+EU) ＇曲线边桩坐标增量计算X=X+I:Y=Y+J ＇坐标计算结果S=SU ＇边长条件Goto8Lbl4 ＇进入主圆计算部分Z=H-A-V÷2 ＇曲线上待求点到HY点的里程O=180Z÷(Rπ) ＇待求点方位角所对应的圆心角X=Rsin O+(V÷2-V^3÷(240R^2))Z=R(1-cos O)+(V^2÷(24R)-V^4÷(2688R^3)) ＇主圆部分待求点支距坐标XY Goto7Lbl5 ＇进入第二缓和曲线计算Z=D-H ＇曲线上待求点到HZ点的里程O=90Z^2÷(πRV) ＇待求点方位角所对应的圆心角X=Z-Z^5÷(40Z R^2V^2)+^9÷(3456R^4V^4)Z=Z^3÷(6RV)-Z^7÷(336R^3V^3)+Z^11÷(42240R^5V^5) ＇待求点支距坐标XY Rec(X,F+GN+180)X=Z[3]+I:Y=Z[4]+J ＇把中桩支距坐标转换成测量统一坐标Rec(Z,F+GN+180-90G)X=X+I:Y=Y+JRec(S,F+GN+180-O G-(180-E)U) ＇曲线边桩坐标增量计算X=X+I:Y=Y+J ＇坐标计算结果S=SUGoto8Lbl6 ＇进入第二直线段计算部分Rec(H-D+T,F+GN) ＇中桩坐标增量计算X=W+I:Y=K+J ＇中桩坐标计算结果Rec(S,F+GN+EU) ＇边桩坐标增量计算X=X+I:Y=Y+J ＇边桩坐标计算结果S=SUGoto8Lbl8X=X◢ ＇显示坐标计算结果Y=Y◢M≠0=>Goto9: ≠>Goto0 ⊿Lbl9 ＇进入放样部分计算Pol(X-M,Y-P) ＇增量计算J<0=>J=J+360: ≠>J=J ⊿J"A"=Int J+0.01Int(60Frac J)+0.006Frac(60Frac J) ◢ ＇方位角计算结果（此显示值为度分秒格式）I"L"=I◢ ＇极距计算结果Goto0∙(二)子程序:AH≤ZY2(ZH2)=>Q=JD1:W=X1:K=Y1:R=R1:F=F0:N=N1:V=Ls1:G=1(-1):=>H≤ZY3(ZH3)=>Q=JD2:W=X2:K=Y2:R=R2:F=F1:N=N2:V=Ls2:G=1(-1):=>H≤ZY4(ZH4)=>Q=JD3:W=X3:K=Y3:R=R3:F=F2:N=N3:V=Ls3:G=1(-1):=>H≤ZY5(ZH5)=>Q=JD4:W=X4:K=Y4:R=R4:F=F3:N=N4:V=Ls4:G=1(-1):=>……………………………H≤ZY n+1(ZH n+1)=>Q=JD n:W=X n:K=Y n:R=R n:F=F n-1:N=N n:V=Ls n:G=1(-1)注:H—待求点桩号ZY n+1(ZH n+1)—后曲线起点桩号 Q—本曲线交点桩号 W—本曲线交点X坐标K—本曲线交点Y坐标 R—本曲线内圆曲线半径 F—前直线方位角 N—本交点转角V—本曲线内缓和曲线长 G—线路转向符,左转角-1,右转角为+1二、高程及放样程序【起算数据】：变坡点桩号A、变坡点高程B、前后纵坡C D、竖曲半径R、横坡度F、最大超高横坡度G、ZH点桩号、YH点桩号、HZ点桩号、施工层厚度N、视线高W （一）主程序：GCFYW"SXG="：N ＇常量输入（视线高程、施工层厚度）Lbl 1｛LО｝:L"K0+0"：О"B=" :Prog"B" ＇变量输入及声明K=C-D ＇坡度差及竖曲线偏角K>0=>J=-1:≠>J=1⊿ ＇凹凸竖曲线判断条件T=Abs(RK)÷2 ＇竖曲线切线长,如需显示则在后面加上◢H=A-T ＇竖曲线起点ZY桩号,如需显示则在后面加上◢X=A+T ＇竖曲线终点YZ桩号,如需显示则在后面加上◢M=L-A ＇待求点至变坡点的距离L≤H=>P=B+MC-N-W: ≠> ＇前纵坡直线段中桩高程计算L≤A=>P=B+MC+J(T-AbsM)2÷(2R)-N-W: ≠> ＇前纵坡竖曲线中桩高程计算L≤X=>P=B+MD+J(T-AbsM)2÷(2R)-N-W: ≠> ＇后纵坡竖曲线中桩高程计算L>X=>P=B+MD-N-W⊿⊿⊿⊿⊿ ＇后纵坡直线段中桩高程计算E=0=>Goto 2 ＇直线部分中、边桩计算判断条件≠>L≤E=> Goto 2 ＇直线部分中、边桩计算判断条件≠>L>V=> Goto 2 ＇直线部分中、边桩计算判断条件≠>L>E=> Goto 3 ＇前缓曲超高部分中、边桩计算判断条件≠> L>U=> Goto 4⊿⊿⊿⊿ ＇后缓曲超高部分中、边桩计算判断条件Lbl 2O=0=>Z“H”=AbsP◢＇显示直线部分中、边桩高程计算结果≠>Y“HB”=Abs(P-OF)◢不需放样显示高程,否则显示放样水准尺计数⊿Goto 1Lbl 3S=（L-E）（F+G）÷（V-U）-F ＇前缓曲线超高横坡度计算S≤F=> Goto 5：≠>S≤G=> Goto 6 ＇超高起点(终点)至零界面、零界面至全超高及全超高≠> Goto 7 ⊿⊿ 部分的计算判断条件Lbl 4S=（V-L）（F+G）÷（V-U）-F ＇后缓曲线超高横坡度计算S≤F=> Goto 5：≠>S≤G=> Goto 6 ＇超高起点(终点)至超高变坡零界面、零界面至全≠> Goto 7 ⊿⊿ 超高、全超高部分计算判断条件Lbl 5O=0=>Z“H1”=AbsP◢＇超高起点(终点)至超高变坡零界面中、边桩高程≠>Y“H-D1”=Abs(P-OF)◢计算结果显示Q“H-G1”=Abs(P+OS)◢⊿Goto 1Lbl 6O=0=>Z“H2”=AbsP◢＇超高变坡零界面至全超高中、边桩高程≠>Y“H-D2”=Abs(P-OS)◢计算结果显示Q“H-G2”=Abs(P+OS)◢⊿Goto 1Lbl 7O=0=>Z“H3”=AbsP◢＇全超高中、边桩高程计算结果显示≠>Y“H-D3”=Abs(P-OG)◢Q“H-G3”=Abs(P+OG)◢⊿Goto 1∙∙注：W“SXG”-视线高N-施工层厚度L"K0+0"-待求点桩号О"B="-中桩至边桩距离（二）子程序：BL≤ZY2=>A=JD1:B=H1:C=I0:D=I1: R=R1:F=I h:E=ZH1:U=YH1:V=HZ1:G=I max1 :≠>L≤ZY3=>A=JD2:B=H2:C=I1:D=I2: R=R2:F=I h:E=ZH2:U=YH2:V=HZ2:G=I max2 :≠>L≤ZY4=>A=JD3:B=H3:C=I2:D=I3: R=R3:F=I h:E=ZH3:U=YH3:V=HZ3:G=I max3 :≠>L≤ZY4=>A=JD3:B=H3:C=I2:D=I3: R=R3:F=I h:E=ZH3:U=YH3:V=HZ3:G=I max3 :≠> ……………………L≤ZY n+1=>A=JD n:B=H n:C=I n-1:D=I n: R=R n:F=I h:E=ZH n:U=YH n:V=HZ n:G=I maxn注:ZY N+1:-下一个竖曲线起点桩号（如：变坡点2的计算范围为前竖曲线终点YZ1至后竖曲线起点ZY3，则在程序中“L≤ZY3”就应输入ZY3的里程桩。
CASIO4800综合曲线坐标计算程序
[正确] CASIO fx-4800P综合曲线坐标计算程序说明：该程序适用于计算器 CASIO fx-4800P，可计算线路中心的缓和曲线、圆曲线、直线段，中、边桩坐标及切线方位角。
1、A? 输入转角:左转为负，右转为正2、R? 输入圆曲线半径3、LS? 输入缓和曲线长度4、JD(DK)? 输入交点里程桩号5、X(JD)? 输入本交点X坐标6、Y(JD)? 输入本交点Y坐标7、FWJ? 输入待求点切线方位角9、J? 输入0程序计算中桩，输入1程序计算边桩10、Z? 输入里程桩号1 A：R：C“LS”：D“JD(DK)”2 P=C∧2/24/R-C∧4/2688/R∧33 Q=C/2-C∧3/240/R∧24 B=90C/兀/R5 T=(R+P)tan(AbsA/2)+Q◢6 W=(R+P)/cos(A/2)-R◢7 L=((AbsA)-2B)兀R/180+2C◢8 G“ZH”=D-T◢9 H“HY”=G+C◢10 I“QZ”=G+L/2◢11 K“YH”=G+L-C◢12 M“HZ”=G+L◢13 N＂X(JD)＂：E＂Y(JD)＂：F＂FWJ＂: J14 A＜0=>S=-1:≠=>S=1⊿ (提示：0为数字“0”)15 U=F+A/2+90S16 V=W+R17 B=N+VcosU18 O=E+VsinU (提示：O为字母“O”)19 Lbl 120 {Z}21 Z≤G=>L=T+G-Z22 V=F+18023 U=F25 Z≤H=>L=Z-G26 V=L-L∧5/(90R∧2C∧2)27 L=30L∧2S/(兀RC)28 P=F+18029 Q=F+L30 U“FWJ”=F+3L◢31 Goto 4⊿32 Z≤K=>L=F+A/2+90S+180+180(Z-I)S/R/兀33 U“FWJ”=L+90S◢34 Goto 5⊿35 Z≤M=>L=M-Z36 V=L-L∧5/(90R∧2C∧2)37 L=30SL∧2/(兀RC)38 P=F+A39 Q=F+A+180-L40 U=F-3L+A◢41 Goto4⊿42 Z>M=>L=Z-M+T43 U=F+A44 V=U45 Goto 246 Lbl 247 X=N+LcosV◢48 Y=E+LsinV◢49 Goto 6⊿50 Lbl 351 {W}52 P“XL”=X+Wcos(U-90) ◢53 Q“YL”=Y+Wsin(U-90) ◢54 P“XR”=X+Wcos(U+90) ◢55 Q“YR”=Y+Wsin (U+90) ◢56 Goto 157 Lbl 458 X=N+TcosP+VcosQ◢59 Y=E+TsinP+VsinQ◢60 Goto 661 Lbl 562 X=B+R cosL◢63 Y=O+RsinL◢ (提示：O为字母“O”)65 Lbl 666 J=1=> Goto 3⊿67 Goto 1注:1、◢为输出指令，若在后面加上，即可显示前面的计算结果输出在屏幕上。
Fx4800p在高斯投影中的应用
‘’'/3Xm ’一 t2 由公式: ‘’ ) 5 “" } 式中子午弧长
J = A
_
一 一B s m i f Co s B - G P
丑
_
. _
s ln ' 五COS B
(32005 . 7799sin B+ 133 .9238(sinB)?3+0. 6976(sin B) ?5+0. 0039(sin B) ?7)cos B: Fix m : X=F+NT(M2 二2+(5- T2+9E+4E) M?4 -. 24+ (6 1- 58T2+T)M?6 二720)Y=N (M+ ( 1一 2+ E)M ?3 - 6+ ( 5一 T 18T 2 +T ? 4+ 14E- 58ET 2)M : 120)
91 1 5 0 +45i)n6 7 4+0.2 (61+9 tf 1 1t1)1
to 二
令则 。=csB 1 ， ‘ o 命0 得
80n- 30(1+2 +'7f)ri3+ tf
(c o s G )2 : N 二 / ( l + E ) : 6 399 69 8 . Y 90178271:B=G一 +E) = 。x T(90N2- 7. 5 (l ( 5 + 3 T 2 + E 一9 E T 2 ) N ? 4 + 0 . 2 5 (61+90T2+45T?4)N?6)U= Frac(Frac B*60) *60L= ( 180N 一 1+ 2T 2+ E)N ?3+ 1 . 5 30( (5+28T2+24T?4)N?5) 二令cos GV= Frac (Frac L*60)*60 注; 在输人的估计纬度 B 进行迭代计算 时， 一定要大于实际纬度。 G 为代表底点纬度Bf ,B 为取纬度的度 和分值， 为纬度B 的秒值;L 为取经度的度 U 和分值; V 为经度的秒值。
fx-4800P型计算器坐标变换程序
fx-4800P型计算器坐标变换程序
杜荣学;段泓娟;段计平
【期刊名称】《地质找矿论丛》
【年(卷),期】2005(020)B08
【摘要】：用CASIO fx-4800P型可编程计算器进行经纬度、3°带和6°带坐标转换应用程序实现了用使用率较高的普通编程汁算器进行一般精度要求的高斯投影面坐标转换的计算，文章介绍了该应用程序及其汁算精度和使用力法。
【总页数】2页(P194-195)
【作者】杜荣学;段泓娟;段计平
【作者单位】中国冶金地质勘查工程总局第三地质勘查院,太原030002;中国冶金地质勘查工程总局三局312队,山西忻州034000;中国冶金地质勘查工程总局三局深蓝公司,太原030002
【正文语种】中文
【中图分类】TP317
【相关文献】
1.CASIO FX-4800P计算器公路放线程序 [J], 兰明
2.fx-4800P型计算器坐标变换程序 [J], 杜荣学;段泓娟;段计平
3.fx-4800P型计算器钻孔内地质点定位程序 [J], 杜荣学;邢福林;段泓娟
4.fx-4800P型计算器钻孔内地质点定位程序 [J], 杜荣学;邢福林;段泓娟
5.罗盘仪导线测量面积Casio fx-4800P型计算器内业计算程序 [J], 陆景星
因版权原因，仅展示原文概要，查看原文内容请购买。
用CASIO fx4800p编高斯投影坐标转换程序
用CASIO fx4800p编高斯投影坐标转换程序
申传明;熊忠招;周志丹
【期刊名称】《资源环境与工程》
【年(卷),期】2007(21)6
【摘要】应用相关的转换公式,采用CASIO fx4800p编程以实现任何参考椭球体下大地坐标与高斯平面坐标之间的相互转换,达到高斯平面坐标转换的目的.【总页数】3页(P739-740,748)
【作者】申传明;熊忠招;周志丹
【作者单位】湖北省国土测绘院,湖北,武汉,430010;湖北省国土测绘院,湖北,武汉,430010;湖北省国土测绘院,湖北,武汉,430010
【正文语种】中文
【中图分类】P226
【相关文献】
1.用CASIO fx4800P编空间直角坐标与大地坐标之间相互转换程序 [J], 申传明
2.测绘工作中自编坐标转换程序和其他坐标转换程序对比 [J], 孙慧梅;李昊;高源;彭秀忠;孟建国
3.Fx4800p在高斯投影中的应用 [J], 刘克福
4.用CASIO fx—4500P进行高斯投影带坐标换算 [J], 张崇安
5.基于VB的高斯投影换带转换程序的实现 [J], 司举香;周小祥;王冬梅
因版权原因，仅展示原文概要，查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
种被称 为“ 投影 ” 的过 程实 现 。 目前 世界 上 有些 国家
采用通用 的横 轴 墨卡托 投 影 的投影 方 法 （ 称 Ｕ Ｍ 投 简 Ｔ 影 ） 在 中国根据 有关 测绘 方 面的 法规规 定 进行 测 量时 ， 统一采用高斯投影 。球 面坐标与平面坐标之 间的相互 Ｊ
３ Ｍ＾ ４ ５ ３÷ ８＋１５ ５÷ ５ ： ＝３ ５ ０ Ｍ＾ ２ ６ Ｊ １ Ｍ ４÷５ ２ Ｋ＝ １：
Ｂ Ｐｏ “ ” Ｐ＝Ｅ 兀 ÷１０ ： ｒｇ Ｐ ： Ｐ ８
Ｘ＝Ｐ＋ Ｎ（ ｃｓ ÷ Ｔ Ｌ ｏＢ） ２×（ Ｉ Ｌ ｏＢ） （ １－（ ｃｓ ５一Ｔ － － Ｉ －
（ ８ １３ ５＋３６ ３ ２ ９ Ｔ４＋１５ ５ ＾ ３ Ｔ ＋４０ ５ ＇ ７ Ｔ ６）÷２ ６ Ｎ ６ ）× ０ １０ ＾ ）
１０÷兀： ８ Ｌ＝Ｙ÷Ｎ ｏＰ×（ Ｉ ｃｓ １－ （一１ 一Ｏ）÷ Ｎ － Ｙ —２ １ Ｉ ２ ６ ＋Ｙ ４ ５－ ８ ２ Ｉ４＋６Ｏ＋８＇ ｆ Ｉ Ｔ ＋ ４＂ － ２ Ｉ 。Ｏ）÷１０ ＾ ２ Ｎ ４＋Ｙ ６ ＾
．
Ｌ ｌ ： ＝Ｙ 一５ ５ Ｋ ＝１０ ÷兀 Ｅ： ｂ Ｙ ３ Ｅ： ８Ｘ Ｆ＝３ ÷２— Ｍ ２ Ｍ＾÷ ２ ２９ ５＋ １ ： ＝ １ ６— ５ ４÷３： ： ７ ３ ３ ＋ ６ Ｍ＇ ５２ Ｈ ２ Ｍ ÷１ ５ Ｍ＾ ２ Ｉ １１ ３＋ ６— １ ５ ２ ： １０７Ｍ＾ ５２ Ｐｏ “ ， ５Ｍ＇ ９ ４７Ｍ＇÷１８ Ｊ＝ ９ ４＋ １： ｒ Ｐ’ ｇ ．
Ｌ ｌ ： Ｂ ： ： ＝ Ｂ： ｒｇ Ｄ” Ｂ ＝Ｗ ： ＝Ｌ： ｂ ｛ Ｌ｝ ＢＬ Ｗ ０ Ｐｏ “ ： Ｗ
Ｐｒ ｇ ｏ “Ｄ’： ＝Ｗ ’Ｌ
Ｌ ｌ： ｒｇ Ｑ ： ｂ２ Ｐｏ “ ” Ｌ＝（ Ｌ—Ｖ）兀 ÷１０ Ｆ＝ 一 Ｍ ÷ ８： ３ ２
＋９Ｍ ＾ ÷ １ — ３ ＾ ÷ ３ Ｈ ： １ Ｍ ÷ １ — １ Ｍ ÷ ３ Ｉ ３ ６ Ｍ ５ ２： ５ ６ ５ ４ ２：
＝ 一
转换 问题 ， 是在一定坐标参照 系下大地坐标 （ ， ） Ｂ Ｌ 与高斯
坐标（ Ｙ） ｘ， 的相互转换 。通过这 样的转换 以实现在一定 坐标参照系下坐标 的正反算 和换带计算 。本文采用 Ｃ — Ａ ＳＯｆ ８０ Ｉ ４０ ｐ的功能和有关 公式 编制它们 的转换 程序 ， ｘ 使
一
Ｂ＝ ：ｒ “ ”Ｂ＝ Ｐ Ｐｏ Ｑ ： Ｐ＋（ ｙ ＋ Ｎ ×（ ｇ Ｔ ２ ２ ２ （一１ Ｏ）Ｉ ５－ 一 －Ｙ （ Ｉ － －
３＂ ＋６Ｏ 一 １ １ ２ ６ ＂Ｏ一３ Ｏ 一９ ＂Ｏ ）÷１ Ｎ ２ １ ２ ２ ＋Ｙ＇（一６ ４ １— ０ ９
Ｉ 一 ５ ＇ ０ ＋１２ａ０－ ５ ＾ ＇ ４Ｔ 。 ４—１７Ｏ ６ ｒ Ｉ Ｔ４Ｏ１÷ ６ Ｎ ４＋Ｙ ６ ｌ ０ ￣ ６ ２ ２ Ｔ４—７ ０ ￣ １ ２ Ｔ６）÷５０ ０ ￣ Ｎ ６） ４
Ｗ ＝ Ｂ： ｒ ｇ Ｐ ｏ “ＤＭ Ｓ”： ： Ｗ “Ｂ ＝ ”： ｕ ｅ Ｐａ ｓ Ｏ：Ｌ ＝ １ ０Ｌ ÷ ８
兀 ＋Ｕ： ＝Ｌ Ｐｏ “ Ｍ ” Ｗ “ Ｗ ： ｒｇ Ｄ Ｓ ： Ｌ＝” ▲
１ 概 况
１ １ 关 于 ｆ ８０ ． ｘ ０ｐ ４
Ｇ＝０ ｛ ： Ｇｏｏ０ ｊ Ｖ｝ Ｖ： ｔ
ｆ ８ ０ 是一种运行速度快 、 带方便 、 ｘ ０ｐ ４ 携 计算 功能 强、 适用性广的计算工具 ， 广泛应用于测绘领域。用它可 随时
Ｇ＝１ ｛ ： Ｇ ｔ １ ｊ Ｕ｝ Ｕ： ｏｏ ｛ Ｖ Ｙ｝Ｕ Ｘ ： ｏｏ３ Ｕ Ｘ ：Ｖ Ｙ Ｇｔ Ｌｌ ： Ｘ ： Ｙ ｂ ｛ Ｙ｝Ｘ １
用 Ｃ ＳＯ ｆ ８０ Ａ Ｉ ４ ０ ｐ编 高 斯 投 影坐 标 转 换 程序 ｘ
申传 明，熊忠招 ，周 志丹
（ 北 省 国 土ｌ 绘 院 ， 湖 测 湖北 武 汉 ４０１） ３０ ０
摘
要 ：应用相关的转换公式 ，采用 Ｃ ＳＯ ｆ ８０ Ａ Ｉ ４０ ｐ编程 以实现任何参 考椭球体 下大地 坐标与 高斯 平面坐标之 间 ｘ
随地实现任何参考椭球体下 高斯 坐标转换与大地 坐标之
间的正反算 和换带 计算 ， 在土地 、 地灾 、 地质 找矿 、 区域调
查等方面的测绘工作 中具有一定的应 用价值 。
１ ２ 坐标 转换 在测 绘工作 中的重 要意 义 ． 坐标 转换包 括 同一 基 准下不 同坐 标形 式 之 问的转 换 、 同基准下 测量坐 标之 间的转 换 、 面 坐标 与 平面 不 球 坐 标之 间 的转 换 等 … ， 测绘 的 目的就 是要 将 不 同形 式 的坐标转 换为 我们 所 需要 的 坐标 ， 后 为 我们 的工 程 然 服务 。因此 坐标 转换在 测绘 工作 中必 不可 少 。球 面坐 标 与平 面坐标之 问 的转 换就 是其 中 的一 种 。 为 了进行球 面坐 标与平面坐标 问 的转 换 ， 以通过 可
维普资讯
第 ２ 卷 ６期 １ 第
２ ０７ 皋 １ 月 ０ ２
资 源 环 境 与 工 程
Ｒｅ ｏ ｒ ｅ ｖｒ ｎ ｅ ｔ＆ Ｅｎ ｉ ｅ ｒｎ ｓ ｕ ｃ ｓＥｎ ｉｏ ｍ ｎ ｇ ｎ ｅ ｉｇ
Ｖｏ ． １ ２１， ． Ｎｏ ６ Ｄｅ ， ０ ｃ．２ ０７
的相 互 转 换 ，达 到 高 斯 平 面 坐标 转换 的 目的 。 关 键 词 ：ｆ ８０ ； 高斯 投 影 坐 标 转 换 ；程序 ｘ ０ｐ ４ 中 图分 类 号 ：Ｐ２ ２６ 文献 标 识 码 ：Ａ 文 章 编 号 ：１７ ６ １—１１ （０ ７ ０ ０ ３ ０ ２ １２ ０ ）６— ７ ９— ２