1、实用测量程序解析中交二航局青荣项目部 何承国摘要曲线坐标计算是测量工作中比较重要的一个环节，本文用浅显易懂的语言配合附图，简明扼要地讲述程序编写及计算式由来，尽可能使掌握该计算方法更容易。关键词测量 程序 解析一、计算思路先设定起始直线方位角为 0，以直缓点为坐标计算起点, 缓和曲线长 L0,圆半径 R,曲线总偏角为 A,E 为曲线偏移值(曲线右偏为正，E 偏向曲线外侧为正)，曲线长为 L。将曲线分为 5 段，即第一直线，第一缓和曲线，圆曲线，第二缓和曲线，第二直线。推得所求点坐标，再进行起始直线方位角换算，得到统一点坐标。二、程序解析（以 CASIO-4500PA 型计算器为例）1、坐标计算
2、程序：ZBTU.CO （用于计算曲线点位坐标）LB1 1:A:M=L”L0”/2-Lxy3/(240RR) （输入曲线总偏角 A，缓和曲线长 L0,圆半径 R，计算切垂距 M，将 L0 值赋予 L）P=LL/(24R)-LxY4/(2688RxY3) （计算圆曲线内移量度p）T”T”=(R+P)tan(A/2)+M （计算切线长 T, 表示显示该值）D”L”=RA/180+L （计算曲线总长 L，赋值予D）B”ZH”:G”Z(-1)Y(1)” （输入 ZH 点里程 B，判别曲线左偏 G=-1,右偏G=1）Lb1 A:Z,E （子程序 A，输入计算点里程 Z，曲线偏移值E）C=Z-B （计算点到
3、 ZH 点曲线长C）C0=X=C （计算点位于 ZH 点前时，如附图 1 所示）Y=-EGoto B （运算结束，进入子程序 B，表示循环结束）CQ=90CC/(RL) （计算点位于第一缓和曲线，如附图 2 所示.,Q 为计算点切线角）X=C-Cxy 5/(40RRLL) +EsinQ （用缓和曲线公式计算点坐标）Y= Cxy 3/(6RL) -Cxy7 /(336(RL) xy3) -ECosQ Goto B （运算结束，进入子程序B）CD-L=Q=90(2C-L)/(R) （计算点位于圆曲线上，如附图 3 所示）X=M+(R+E) SinQ （用坐标投影公式计算点坐标）Y=R+P-(R+E
4、)CosQ Goto B （运算结束，进入子程序 B）CD=S=D-C:Q=90SS/(RL) （计算点位于第二缓和曲线，如附图 4 所示.,Q为计算点切线角,计算先以 HZ 点为坐标原点，S 为自 HZ 点至计算点曲线长）O=S-Sxy 5/(40RRLL) +EsinQ （用缓和曲线公式计算点坐标，相对 HZ 点，O为 X 坐标值，U 为 Y 坐标值） U= Sxy 3/(6RL) -S xy7 /(336(RL) xy3) -ECosQX=-OcosA-UsinA+T(1+CosA) （坐标换算至独立坐标系）Y=-OsinA+UcosA+TsinA Goto B （运算结束，进入子程序
5、B）X=T+(C-D+T)CosA+EsinA （计算点位于第二直线（CD） ，如附图 5 所示.,Q为计算点切线角）Y=(C-D+T)SinA-EcosA （运算结束，自然进入下一程序 B） LB1 B：H”XD”=X （计算独立坐标 XD,YD 值）I”YD”=GYJ”XT”=XcosF-YsinF （计算统一坐标 XT,YT 值，如附图 6 所示 F 为第一直线方位角，将 XT 值赋予 I，将 YT 值赋予J） K“YT”=XsinF+YCosF Goto 1 2、坐标放样程序：ZBFY.COLB1 1：A“XH”:B”YH” （输入后视及置镜点坐标）C”XZ”:D”YZ”Pol(A-C
6、),(B-D) （极坐标转换）E”SH”=V （输出后视距离）F=W:X,Y （输入前视点坐标）K=X:G=YPol(K,G) （极坐标转换）S=V （输出前视距离）P=W-F：PP=P+360=P=P （输出顺拔角）Goto 13、测量坐标程序：ZBTS.COLB1 1：A“XH”:B”YH” （输入后视及置镜点坐标）C”XZ”:D”YZ”Pol(A-C),(B-D) （极坐标转换）E”SH”=V （输出后视距离）F=W:S,P （输入前视点距离及顺拔角）K=S:G=F+PRec(K,G) （极坐标转换）X=V+C （输出 X 坐标）Y=W+D （输出 Y 坐标）Goto 1三、程序说明：1、各小程序可合并为一个程序，但须扩展存储区。2、 程序的使用步骤见 CASIO4500PA 说明书。
