测量程序解析.doc

实用测量程序解析中交二航局青荣项目部 何承国摘要曲线坐标计算是测量工作中比较重要的一个环节，本文用浅显易懂的语言配合附图，简明扼要地讲述程序编写及计算式由来，尽可能使掌握该计算方法更容易。关键词测量 程序 解析一、计算思路先设定起始直线方位角为0，以直缓点为坐标计算起点, 缓和曲线长L0,圆半径R,曲线总偏角为A,E为曲线偏移值(曲线右偏为正，E偏向曲线外侧为正)，曲线长为L。将曲线分为5段，即第一直线，第一缓和曲线，圆曲线，第二缓和曲线，第二直线。推得所求点坐标，再进行起始直线方位角换算，得到统一点坐标。二、程序解析（以CASIO-4500PA型计算器为例）1、坐标计算程序：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 （计算点到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)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 （运算结束，进入子程序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),(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) （极坐标转换）