K0+540577.493 缓和曲线 1357.277 1364.654 2655.881 2692.600 18.24 18.24 1374.734 1382.687 2650.593 2689.862 18.20 18.20 1339.859 1346.660 2661.157 2695.332
K0+600+640 直线段 1368.032 1374.027 2714.852 2754.399 18.24 18.24 1386.065 1392.070 2712.114 2751.661 18.20 18.20 1350.038 1356.043 2717.584 2757.131 5.结论
(1)本文程序是采用n=2时的复化辛甫生公式； (2)本程序适用于线路的各种线型；
(3)本程序既可按里程增加方向计算也可按里程减少方向计算，计算时要注意曲线的偏向。
4800线路放样程序（中边桩）
最近有人还是问这些程序！再提供一次，希望以后没有人再我4800线路程序了 4800P极坐标放样输入步骤(中边桩放样) A
XYD“X1”E“Y1”F“X2”G“Y2”H“X3”V“Y3”RL“LS”W“JDK”
I=D-F：J=E-G：Porg“C”：O=J：I=H-F：J=V-G：Porg“C”：U=J
M=O-J：M<0=>M=M+360⊿M“A”=180-M◢M<0=>C=-1：≠> C=1⊿M=AbsM：
P=L2÷(24R)-L^4÷(2688R^3)：Q=L÷2- L^3÷(240R2)： (R+P)tan(M÷2)+Q：Prog“Z”： T=Ans
◢: I“E”=(R+P)÷cos(M÷2)-R◢RπM÷180+L：Prog“Z”： S=Ans◢I“ZH”=W-T◢I“HY”=W-T+L◢I“QZ”=W-T+S÷2◢I“YH”=W-T+S-L◢I“HZ”=W-T+S◢:Porg“B” B
Lbl 0：Z“Z=90”：Fixm：{KB}：I=K-W+T：J=Z：B：I
I≤0=>\ZX2\:M=I-B cosZ：N=-BsinZ：≠>\LS2\:C=-1=>J=180-Z⊿Porg“D”：N=-N⊿J=U：Porg“G”：Goto0 C
11
Pol(I，J)：J<0=>J=J+360⊿ D
A=RL：M=I-I^5÷(40A^2)+I^9÷(3456A^4)-I^13÷599040A^6+I^17÷175472640A^8：
N=I^3÷(6A)-I^7÷(336A^3)+I^11÷42240A^5-I^15÷9676800A^7+I^19÷3530096640A^9：A=90I^2÷(Aπ)：Porg“E” E
A=A-J：M=M-CBcosA：N=CN-BsinA F
A=90÷(πR)×(2I-L)：M=RsinA+Q：N=P+R(1-cosA)：Porg“E” G
A=cosJ：J=sinJ：I=F+TA-MA+NJ：J=G+TJ-MJ-NA Fix3:I“X”=I◢ J“Y”=J◢ I=I-X:J=J-Y
Porg“C”: I“JU LI(m)”=I◢Norw:J“FWJ(°°°)”=J◢ Z
Fix3：Rnd：Norm
以下程序是把所有的交点坐标全部存入计算机里 文件名：JJ Lbl O
J=1=>D交点坐标E=: F=: G= :H=: V= :R=: L=; W=JDK; Porg“A”:GotoO
用此程序可直接将所有的交点坐标及半径、LS等全部存入计算器，放桩号只要交点点号，它自己会判断。 R: 半径 LS: 缓和曲线长度 A: 偏角 T: 切线长 E: 外距 JDK: 交点桩号 K: 放样中桩桩号 FWJ: 方位角
JU LI: 距离 B: 距中桩宽度(放边桩左+右-)
CASIO4500缓和曲线、圆曲线、直线段坐标放样计算程序 CASIO4500缓和曲线、圆曲线、直线段坐标放样计算程序
1、 该程序适用于计算器 CASVO fx-4500，可计算与线路中心成任意夹角的缓和曲线、圆曲线、直线段中、边桩坐标、切线方位角和放样计算程序。
2、 PJ？ 输入转角：左转为负，右转为正 3、 R？ 输入圆曲线半径 4、 LS？ 输入缓和曲线长度 5、 A=？ 输入起始方位角 6、 XJ 、YJ？ 输入交点X、Y坐标 7、 X0、Y0？ 输入测站点X、Y坐标 8、 Z？ 输入交点桩号 7、 B？ 输入左右边距
8、 U 待求点切线方位角(大于360-360) 9、 W 放样点方位角 10、XL、YL 放样点X、Y坐标
12
11、如计算交角为60°时，则公式改为（U-60），（U+120）
1 A“PJ”：R：C“LS”：D“JD” 2 P=C∧2/24/R-C∧4/2688/R∧3 3 Q=C/2-C∧3/240/R∧2 4 B=90C/兀/R
5 T=(R+P)tan(AbsA/2)+Q◢ 6 W=(R+P)/cos(A/2)-R 7 L=(AbsA-2B)兀R/180+2C◢ 8 G=D-T◢ 9 H=G+C◢ 10 I=G+L/2◢ 11 K=G+L-C◢ 12 M=G+L◢
13 N “A=” ： E “XJ”：F“YJ”: X“X0=” : Y“Y0=” 14 A＜0=>S=-1:≠=>S=1⊿◢ 15 U=N+A/2+90S 16 V=W+R 17 J=N+VcosU 18 O=E+VsinU 19 Lbl 1 20 {Z}
21 Z≤G=>L=T+G-Z 22 V=F+180 23 U=F 24 Goto 2⊿ 25 Z≤H=>L=Z-G
26 V=L-L∧5/(90R∧2C∧2) 27 L=30L∧2S/(兀RC) 28 W=N+180 29 L=N+L 30 U=N+3L◢ 31 Goto 4⊿
32 Z≤K=>L=N+A/2+90S+180+180(Z-I)S/R/兀 33 U=L+90S◢ 34 Goto 5⊿ 35 Z≤M=>L=M-Z
36 V=L-L∧5/(90R∧2C∧2) 37 L=30SL∧2/(兀RC) 38 W=N+A 39 L=N+A+180-L 40 U=N-3L+A◢ 41 Goto 4⊿
42 Z>M=>L=Z-M+T
13
43 U=F+A◢ 44 V=U 45 Goto 2⊿ 46 Lbl 2
47 P=E+LcosV◢ 48 Q=F+LsinV◢ 49 Goto 3⊿ 50 Lbl 3 51 {B}
52 P“XL”=P+Bcos(U+90) ◢ 53 Q“YL”=Q+Bsin(U+90) ◢ 54 Goto 6 55 Lbl 4
56 P=E+TcosW+ VcosQ◢ 57 Q=F+ TsinP+ VsinQ◢ 58 Goto 3 59 Lbl 5 60 P=J+R cosL 61 Q=O+RsinL 62 Goto 3 63 Lbl 6
64 L“L=”Pol（（P-X），（Q-Y））
65 W ≥0=>W“W=”W:≠=> W“W” =W+360⊿◢ 66 L=L◢ 67 Goto 1
[摘要] 经过一定范围的市场调查，从用户的角度对目前投入使用的数字化测图软件在性能和适用方面进行分析和实际工作的坐标法测图的应用，不断对坐标法测图应用于土方测量计算进行探讨优化。细述计算机以清华山维公司研制的EPSW电子平板测图系统为工作平台进行土方量三角网法计算的数据处理及操作应用。 [关键词] 三角网 底高程 顶高程
土石方量的计算是工程费用概算及方案选优的重要因素，故工程施工前的设计阶段必须对土石方量进行预算，工程施工后的竣工阶段必须对土石方量进行结算。土石方量预算是设计高程作为底高程，以施工前现状地形高程作为顶高程。土石方量结算是施工后现状地形高程作为底高程，以施工前现状地形高程作为顶高程。土石方量的计算方法有方格网法，等高线法，断面法和三角网法 。在实际测量时，无论是方格网还是作断面，其工作量都较大，而且一旦选定方格网的起始边或断面便很难反映地形的实际特征，使得测出的土方量与实际的土方量有较大的出入。等高线法其存储数据量大，数据结构复杂和难以建立。而三角网法构成的系统性能最优，并能克服地形起伏不大的地区产生冗余数据的问题，精度高于方格网法，等高线法，断面法精度。 一．计算方法：
在原图（施工前现状地形图）中建立数字模型（DTM）。在收方图（施工后的竣工地形图或设计图）中建立数字模型（DTM）。计算在原图中建立数字模型（DTM）的（顶高程）与在收方图中建立数字模型（DTM）的（底高程）相减的数字模型（DTM）。通过底高程与顶高程相减的数字模型（DTM）来计算土方量。即计算在原图中的各坐标点对应在收方图中各坐标点的高程，求得各坐标点相对应的坐标点高差。同样计算在收方图中的各坐标点对应在原图中各坐标点的高程，求得各坐标点相对应的坐标点高差。在土方计算图中，用两个图的坐标点高差数据进行数字模型（DTM）建立，绘制等高线。为了避免三角网跨过填挖方分界线，加密填挖方分界线（即0等高线）的坐标高程点数据，并加入两个图的坐标点高差数据，重新建立数字模型（DTM），生成三角网。通过三角网来计算每个三棱锥的填挖方量，
14
最后累计得到指定范围内填方和挖方的土方量。
二．以清华山维公司研制的EPSW电子平板测图系统为工作平台的操作过程： １．建立文件夹：原图图 ，施测图 ，土方图 COR ，TF
２．在施测收方图中，进行编图，建模。提取各坐标点的高程，即（800）区域输出（*1.cor），（*1.cor）用写字板替换（编码为（808），不建摸）在文件夹COR（in1.cor）。
３．在原图中，进行编图，建模。提取各坐标点的高程，即（800）区域输出（*2.cor），（*2.cor）用写字板替换（编码为（808），不建摸）在文件夹COR（in2.cor）。
４．求得各坐标点相对应的底高程与顶高程相减高差：调人in1.cor，选线（808）处理（lxyz.in），提取断面数据（lxyz.out），运行out-inzcx.exe程序，获取坐标点高差数据,( lxyz.in)（lxyz.out） ---> (out-in.cor)，得到（out.cor），（out-in.cor）。再改（（lxyz.in）（lxyz.out）（out.cor），（out-in.cor））为文件夹COR（（lxyz1.in）（lxyz1.out）（out1.cor），（out-in1.cor））。再做施测图，调人in2.cor，选线（808）处理（lxyz.in），提取断面数据（lxyz.out），且做( lxyz.in)与（lxyz.out 相对换。运行out-inzcx.exe程序，获取坐标点高差数据，( lxyz.in)（lxyz.out） ---> (out-in.cor)，得到（out.cor），（out-in.cor）。再改（（lxyz.in）（lxyz.out）（out.cor），（out-in.cor））为文件夹COR（（lxyz2.in）（lxyz2.out）（out2.cor），（out-in2.cor））。
５．建立相减的数字模型（DTM）：在土方图中，进行编图，建模。调人（out-in1.cor），（out-in2.cor），加密填挖方分界线坐标高程点，建摸（DTM），得到（GET.pn1），（GET.tr1），转入文件夹TF（GET.pn1），（GET.tr1）
６．土方量计算:（引用坐标法测图应用于土方测量计算）：运行TFcenxi.exe程序，（HH0=0）得到在文件夹TF（TF.cor），（TF.not）。再做土方图，调人（TF.cor）整理出图。 三，算例
１．坐标点相对应的底高程与顶高程高差坐标文件格式及说明： 0 : 606198.991 700305.820 -0.326 1 : 606279.739 700646.650 0.047 2 : 606299.719 700661.439 -0.097 ：：：： ：：：：
77 : 606363.537 700558.572 -0.182 78 : 606350.202 700638.175 -0.013 点号：坐标X 坐标Y 高差dH ：：：： ：：：：
点号：坐标X 坐标Y 高差dH ２．土方计算的数据文件格式及说明： 1：2718253 505340.4 224.015 2：2718242 505334.5 162.938 3：2718224 505336.1 177.211 ：：：： ：：：：
104：2718211 505356.3 -75.002 105：2718220 505374.8 -88.791
三角形序号：中心坐标X 中心坐标Y 土方量Ｕ ：：：： ：：：：
三角形序号：中心坐标X 中心坐标Y 土方量Ｕ
15