大地坐标转换成施工坐标公式.docx

1、大地坐标转换成施工坐标公式 Company number：【0089WT-8898YT-W8CCB-BUUT-202108】大地坐标转换成施工坐标公式大地(高斯平面)坐标系工程坐标系转换大地坐标系-工程坐标系=待转换点为P，大地坐标为：Xp、Yp工程坐标系原点o:大地坐标：Xo、Yo工程坐标：xo、yo工程坐标系x轴之大地方位角：adX=Xp-XodY=Yp-YoP点转换后之工程坐标为xp、yp:xp=dX*COS(a)+dY*SIN(a)+xoyp=-dX*SIN(a)+dY*COS(a)+yo工程坐标系-大地坐标系=待转换点为P，工程坐标为：xp、yp工程坐标系原点o:大地坐标：Xo、Yo

2、工程坐标：xo、yo工程坐标系x轴之大地方位角：adx=xp-xody=yp-yoP点转换后之工程坐标为xp、yp:xp=Xo+dx*COS(a)-dy*SIN(a)yp=Yo+dx*SIN(a)+dy*COS(a)坐标方位角计算程序置镜点坐标：ZX?ZY后视点坐标：HXHY方位角：W两点间距离:SLb10A,B,C,DAZX=:BZY=:CHX=:DHY=:W=tg1(D-B)(C-A):(D-B)0=(C-A)0=W=W:(D-B)0=(C-A)W=W+180:(D-B)(C-A)W=W+180:(D-B)(C-A)0=W=360+WW=WS=(D-B)2+(C-A)2)Goto?0CAS

3、IO?fx4500p坐标计算程序根据坐标计算方位角WW360W：“ALF(12)”L1A“X1”：B“Y1”：Pol(C“X2”A，D“Y2”B：“S”W0直线段坐标计算L1X“X(0)”：Y“Y(0)”：S“S(0)”：A“ALF”L2Lb12L3L：L“LX”L4M“X(Z)”X(LS)cosAL5N“Y(Z)”Y(LS)sinAL6B：B“B(L)”：Q“Q”L7O“X(L)”MBcos(AQ180)L8P“Y(L)”NBsin(AQ180)L9C：C“B(R)”L10U“X(R)”MCcos(AQ)L11V“Y(R)”NCsin(AQ)L12Goto2园曲线段坐标计算L1S“S(0)-

4、Km”：X“X(0)”：Y“Y(0)”：A“ALF”：R“R”：K“K(L1,R2)”L2Lb12L3L：L“L(X)”L4V180（LS）R：WV2L5CA(-1)KW：D2RsinW：FA(-1)KVL6M“X(Z)”XDcosCL7N“Y(Z)”YDsinCL8E：E“B(L)”：Q“Q”L9O“X(L)”MEcos(FQ180)L10P“Y(L)”NEsin(FQ180)L11G：G“B(R)”L12T“X(R)”MGcos(FQ)L13U“Y(R)”NGsin(FQ)L14Goto2正向缓和曲线段坐标计算L1S“ZH-Km”：X“X(ZH)”：Y“Y(ZH)”：A“ALF”：R“R”

5、：H“LS”：K“K(L1,R2)”L2Lb12L3L：L“L(X)”L4D30（LS）2RH：CLS（LS）590（RH）2：BAD(-1)K：EA3D(-1)KL5U“X(Z)”XCcosBL6V“Y(Z)”YCsinBL7G：G“B(L)”：Q“Q”L8F“X(L)”UGcos(EQ180)L9I“Y(L)”VGsin(EQ180)L10J：J“B(R)”L11M“X(R)”UJcos(EQ)L12N“Y(R)”Vjsin(EQ)L13Goto2卵形曲线坐标计算X1,D2)”L1S“Km-YH”：E“X(YH)”：F“Y(YH)”：G“ALF”：B“R1”：D“A”：K“K(L1,R2)

6、”：Q“R1-R2L2Lb12L3Z：Z“L(X)”L4J“L1”D2B：R“RP”D2B(D2(-1)Q(ZS)B)：L“LP”D2RL5M(LJ)(L5J5)40D4(L9J9)3456D8L6N(L3J3)6D2(L7J7)336D6(L11J11)42240D10L7TG(-1)Q(-1)KJ290D2L8X“X(Z)”E(-1)QMcosT(-1)KNsinTL9Y“Y(Z)”F(-1)QMsinT(-1)KNcosTL10A“ALF(P)”G(-1)K(ZS)90(1B1R)L11H：H“B(L)”：U“Q”L12W“X(L)”XHcos(AU180)L13V“Y(L)”YHsin

7、(AU180)L14C：C“B(R)”L15I“X(R)”XCcos(AU)L16P“Y(R)”YCsin(AU)L17Goto2公路逐桩坐标计算4800程序公路逐桩坐标计算程序(可以计算对称、不对称缓和曲线)Lb10Z=V=W=V+2：FixmKLb11KZW+5Z+4=W=W+1:Goto1(判断桩号在哪个交点范围，就是该交点曲线起点至下一交点曲线起点)S=K-ZW+5Z+3（计算该桩号与曲线起点的距离）R=ZW+2Z+2：L=ZW+3Z+2:E=ZW+4Z+2（读取该交点曲线要素R、Ls1、Ls2）Pol（ZW-ZW-1,ZW+Z+2-ZW+Z+1）（计算该交点与下一交点直线方位角）JJ

8、=J+360A=JPol（ZW-1-ZW-2,ZW+Z+1-ZW+Z）（计算该交点与上一交点直线方位角）JJ=J+360C=A-J:A=J（计算偏角）W=V+2=Goto2（如果桩号在起点与第一交点曲线起点之间，则转Lb12）I=Abs（tan（c2）M=L2-L3240R2:N=E2-E3240R2P=L26R-L4336R3-R（1-cos（90LR）Q=E26R-E4336R3-R（1-cos（90ER）D=（P-Q）I2:F=（P+Q+2R）I2M=F+M-D:Q=F+N+DN=RAbsC180+（L+E）2X=ZW-1-McosAY=ZW+Z+1-MsinAM=ZW-1+Qcos（A

9、+C）V=ZW+Z+1+Qsin（A+C）Q=AbsCCSL=P=0:Goto3（如果桩号在第一缓和曲线内，则转Lb13）SN-E=S=S-L:Goto4（如果桩号在圆曲线内，则转Lb14）SN=S=N-SQ=-Q:A=A+C-180:X=M:Y=V:L=E：P=180:Goto3（如果桩号在第二缓和曲线内，则转Lb13）P=A+C:S=S-N:D=M+ScosP:F=V+SsinPGoto6（如果桩号在直线内，则转Lb16）Lb12P=A+CD=ZW-1+ScosPF=ZW+Z+1+SsinP:Goto6Lb13I=S-S540R2L2+S93456R4L4J=Q（S36RL-S7336R3

10、L3）P=P+A+90QS2RL:Goto5Lb14M=90（2S+L）RI=RsinM+L2-L3240R2J=Q（L224R+R（1-cosM）P=A+QMLb15D=X+IcosA-jsinA:F=Y+JcosA+IsinALb16DX=（结果显示X坐标）FY=（结果显示Y坐标）PAT=（结果显示该桩号方位角）BO：BSO（输入边桩距离，交角）P=P+OLXB=D+BcosP（结果显示边桩X坐标）MYB=F+BsinP（结果显示边桩Y坐标）以上是坐标计算程序，括号内是程序计算的大致原理及说明，中间部分为直线、圆曲线、缓和曲线计算的各种公式，大家也知道，书上也有。该程序是最基本的，如再加几

11、条语句就可以处理断链。该程序还需变量输入，把整条线路曲线要素一次输入到4800中，就可以计算线路任意中、边桩坐标。曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序曲线任意里程中边桩坐标正反算(CASIOfx-4800P计算器)程序一、程序功能本程序由一个主程序(TYQXjs)和两个子程正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的

12、两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIOfx-4500P计算器及CASIOfx-4850P计算器上运行。二、源程序1.主程序(TYQXjs)=XY：=SZ：N：UX0：VY0：OS0：GF0：HLS：PR0：RRN：Q：C=1P：D=(P-R)(2HPR)：E=180：N=1=Goto1：Goto2Lbl1：SZ：SZ：W=Abs(S-O)：ProgSUB1：XXS=XYYS=YGoto1Lbl2：XY：XY：I=X：J=Y：ProgSUB2：SS=O+WZZ=ZGoto22.正算子程序(SUB1)M=1-K：X=U+W(Acos(G+

13、QEKW(C+KWD)+Bcos(G+QELW(C+LWD)+Bcos(G+QEFW(C+FWD)+Acos(G+QEMW(C+MWD)：Y=V+W(Asin(G+QEKW(C+KWD)+Bsin(G+QELW(C+LWD)+Bsin(G+QEFW(C+FWD)+Asin(G+QEMW(C+MWD)：F=G+QEW(C+WD)+90：X=X+ZcosF：Y=Y+ZsinF3.反算子程序(SUB2)T=G-90：W=Abs(Y-V)cosT-(X-U)sinT)：Z=0：Lbl0：ProgSUB1：L=T+QEW(C+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZGoto1：W=W

14、+Z：Goto0Lbl1：Z=0：ProgSUB1：Z=(J-Y)sinF三、使用说明1、规定(1)以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。(2)当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。(3)当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。(4)当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。(5)当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半

15、径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。(6)当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明输入部分：1.SZ=XY2.XY=SZN选择计算方式，输入1表示进行由里程、边距计算坐标；输入2表示由坐标反算里程和边距。X0线元起点的X坐标Y0线元起点的Y坐标S0线元起点里程F0线元起点切线方位角LS线元长度R0线元起点曲率半径RN线元止点曲率半径Q线元左右偏标志(左偏Q

16、=-1，右偏Q=1，直线段Q=0）S正算时所求点的里程Z正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)X反算时所求点的X坐标Y反算时所求点的Y坐标显示部分：XS=正算时，计算得出的所求点的X坐标YS=正算时，计算得出的所求点的Y坐标S=反算时，计算得出的所求点的里程Z=反算时，计算得出的所求点的边距四、算例某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组成，各段线元的要素（起点里程S0、起点坐标X0Y0、起点切线方位角F0、线元长度LS、起点曲率半径R0、止点曲率半径RN、线元左右偏标志Q）如下：S0X0Y0F0LSR0RNQ125161E451E45

17、0125161E45-112025-19117-180401E451E4501、正算（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）S=700Z=-5计算得XS=YS=S=700Z=0计算得XS=YS=S=700Z=5计算得XS=YS=S=780Z=-5计算得XS=YS=S=780Z=0计算得XS=YS=S=780Z=5计算得XS=YS=S=870Z=-5计算得XS=YS=S=870Z=0计算得XS=YS=S=870Z=5计算得XS=YS=S=940Z=计算得XS=19741.59118YS=S=940Z=0计算得XS=YS=S=940Z=计算得XS=YS=2、反算

18、X=Y=计算得S=Z=-5.00018164X=Y=计算得S=Z=X=Y=计算得S=Z=X=Y=计算得S=Z=-5.000001663X=Y=计算得S=Z=-X=Y=计算得S=Z=-X=Y=计算得S=Z=-X=Y=计算得S=Z=point=10/point自己写的fx-4800小程序一、ZBFS(坐标反算)二、NEXY,三、Fixm。四、F,五、Pol(N-X,E-Y)。六、J0J=J。七、J0J=J+360。八、九、说明：“，”表示回车符；“。”表示终止符。十、十一、二、ZBFS-2(坐标反算)十二、NEXY,十三、Pol(N-X,E-Y)。十四、J0J=J。十五、J0J=J+360。S=I

19、。R=J。说明：“，”表示回车符；“。”表示终止符。三、ZBZS(坐标正算）SRXY,V=X+Rec(S,R),W=Y+J,N=V。E=W用全站仪测设公路中桩和边桩程序程序说明：本程序用于用全站仪测设中桩和边桩，采用坐标变换的方法编制而成，在使用fx-4500计算器时，可直接输入以下程序便可使用；而在使用fx-4800计算器时，见fx-4800程序清单。注意：因程序的0(零)和O易混淆,所以字母O用较大的全角符号，零仍用标准字。程序输入时应连续输入，无空格。而且输入时文件名需和给出的一样，否则会出现调用文件错误。变量说明：X输入摆站点的x坐标C输入第一个曲线角桩x坐标I输入中间曲线角桩的x坐标

20、Y输入摆站点的y坐标D输入第一个曲线角桩y坐标J输入中间曲线角桩的y坐标E最后一个曲线角桩x坐标R中间曲线半径L中间曲线缓和曲线长F最后一个曲线角桩y坐标K第一个K为中间角桩桩号S中间角桩曲线总长，检验用K第二个K为要测设的桩号B测中桩时输入0，测边桩时输入宽度（左输负值，右输正值）程序清单注释文件:F1文件名：1L1L2L3Defm8XYCDIJEFRLK:Fixm:Z2=0:ProgBZ5=0.5L-LXy3/240R2:Z3=Z5+(R+L2/24R)AbstanT2:Z=K-Ans:S=RAbsT/180+L:SH=Z+Ans:G=RL:Prog2文件:F2文件名：2L1Lbl0:Fi

21、xm:KL2Lbl6:B:Norm:K(-1)XyZ2(H-L)(-1)XyZ2Goto4FracZ22=0=Z6=Z:Z6=HA=90T/AbsTL3BN=0:N=1B=AbsBL4=Abs(K-Z6:K(-1)XyZ2Z6(-1)XyZ2Goto2Q=Z3+:Z4=BL5B0ProgCProgD:P=0L6ProgE:Goto0L7Lbl2:LGoto3L8Q=Z3-+Xy5/40G2:Z4=Xy3/6G-Xy7/336GXy3:ProgD:P=BL9B0ProgFProgE:Goto0L10Lbl3:Q=Z3-Z5:Z4=R+L2/24R:ProgDL11Z8=U+A(2(-L)/R-1

22、+L/R):B0ProgG:P=RProgE:Goto0L12Lbl4:M=C:=D:C=E:D=F:E=M:F=:Z2=Z2+1:ProgB:Goto6文件:F3文件名：0L1Pol(M,):WW=W+360文件:F4文件名：BL1L2M=E-I:=F-J:Prog0:T=W:M=I-C:=J-D:Prog0:U=W:T=T-WTT=T+360T180T=T-360文件:F5文件名：CL1A=-90(-1)Xy(N+Z2)文件:F6文件名：DL1Z7=I-QcosU+Z4cos(U+A):Z1=J-QsinU+Z4sin(U+A)文件:F7文件名：EL1L2M=Z7+PcosZ8-X:=Z1+PsinZ8-Y:Prog0WFix3:VNorm文件:F8文件名：FL1Z8=U+A+A2/G+90+A(-1)Xy(N+Z2)文件:F9文件名：GL1P=R+AB(-1)Xy(N+Z2)/90此处L2句加入：M+X+Y可显示测点坐标值CASIO 4800超高及高程（源程序）超高及高程（源程序）L1:R:L“LS”:Z“ZH”:H“HZ”:V“RS”:F“I1”:G“I2”:O“BPH”:N“BP”:T=Vabs(GF)2:X“X(R+1,L-1)”:P“P(R+,L-1)”L2:Lbl6:SL3:S+TN0D=O+(SN)F:SN0D=O+(SN)FU(SN+T)2