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

《大地坐标转换成施工坐标公式.doc》由会员分享，可在线阅读，更多相关《大地坐标转换成施工坐标公式.doc（12页珍藏版）》请在装配图网上搜索。

1、大地(高斯平面)坐标系工程坐标系转换大地坐标系-工程坐标系=待转换点为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工程坐标：xo、yo工程坐标系x轴之大地方位角：adx=xp-xody=yp-yoP点转换后之工程坐标为xp、yp:xp=Xo+dx*COS(a)-d

2、y*SIN(a)yp=Yo+dx*SIN(a)+dy*COS(a)坐标方位角计算程序置镜点坐标：ZXZY后视点坐标：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)Goto0CASIOfx4500p坐标计算程序根据坐标计算方位角WW360W：“ALF(12)”L1A“X1”：B“Y1”：Pol(C“X2”A，D“Y2”B：“S”W0

3、直线段坐标计算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)-Km”：X“X(0)”：Y“Y(0)”：A“ALF”：R“R”：K“K(L1,R2)”L2Lb12L3L：L“L(X)”L4V180（LS）R：WV2L5C

4、A(-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”：H“LS”：K“K(L1,R2)”L2Lb12L3L：L“L(X)”L4D30（LS）2RH：CLS（LS）590（RH）2：BAD(-1)K：EA3D(

5、-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)”：Q“R1-R2L2Lb12L3Z：Z“L(X)”L4J“L1”D2B：R“RP”D2B(D2(-1)Q(ZS)B)：L“LP”D2RL5M(LJ)(L5

6、J5)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(AU180)L14C：C“B(R)”L15I“X(R)”XCcos(AU)L16P“Y(R)”YCsin(AU)L17Goto2公路逐桩坐标计算4800程

7、序公路逐桩坐标计算程序(可以计算对称、不对称缓和曲线)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=J+360A=JPol（ZW-1-ZW-2,ZW+Z+1-ZW+Z）（计算该交点与上一交点直线方位角）JJ=J+360C=A-J:A=J（计算偏角）W

8、=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+C）V=ZW+Z+1+Qsin（A+C）Q=AbsCCSL=P=0:Goto3（如果桩号在第一缓和曲线内，则转Lb13）SN-E=S=S-L:Goto

9、4（如果桩号在圆曲线内，则转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-S7336R3L3）P=P+A+90QS2RL:Goto5Lb14M=90（2S+L）RI=RsinM+L2-L3240R2J=Q（L224R+R（1-cosM）P=

10、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坐标）以上是坐标计算程序，括号内是程序计算的大致原理及说明，中间部分为直线、圆曲线、缓和曲线计算的各种公式，大家也知道，书上也有。该程序是最基本的，如再加几条语句就可以处理断链。该程序还需变量输入，把整条线路曲线要素一次输入到4800中，就可以计算线路任意中、边桩坐标。如有需要，需RMB￥200，附有变量输

11、入方法及详细使用说明。联系：13858423737曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序曲线任意里程中边桩坐标正反算(CASIOfx-4800P计算器)程序一、程序功能本程序由一个主程序(TYQXjs)和两个子程正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进

12、行正反算。本程序也可以在CASIOfx-4500P计算器及CASIOfx-4850P计算器上运行。二、源程序1.主程序(TYQXjs)1.SZ=XY：2.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)A=0.1739274226：B=0.3260725774：K=0.069431844

13、2：L=0.3300094782：F=1-L：M=1-K：X=U+W(Acos(G+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

14、(C+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZGoto1：W=W+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)当线元为完整缓和曲线时，起点与直线

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

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

17、率半径R0、止点曲率半径RN、线元左右偏标志Q）如下：S0X0Y0F0LSR0RNQ500.00019942.83728343.5611251631.00269.2561E451E450769.25619787.34028563.3781251631.0037.4921E45221.75-1806.74819766.56628594.5741202554.07112.779221.75221.75-1919.52719736.07228701.893911730.6380.285221.759579.228-1999.81219744.03828781.659804050.00100.0001

18、E451E4501、正算（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）S=700Z=-5计算得XS=19831.41785YS=28509.72590S=700Z=0计算得XS=19827.33592YS=28506.83837S=700Z=5计算得XS=19823.25398YS=28503.95084S=780Z=-5计算得XS=19785.25749YS=28575.02270S=780Z=0计算得XS=19781.15561YS=28572.16358S=780Z=5计算得XS=19777.05373YS=28569.30446S=870Z=-5计算得XS

19、=19747.53609YS=28654.13091S=870Z=0计算得XS=19742.68648YS=28652.91379S=870Z=5计算得XS=19737.83688YS=28651.69668S=940Z=-5.123计算得XS=19741.59118YS=28722.05802S=940Z=0计算得XS=19736.47687YS=28722.35642S=940Z=3.009计算得XS=19733.47298YS=28722.531682、反算X=19831.418Y=28509.726计算得S=699.9999974Z=-5.00018164X=19827.336Y=28

20、506.838计算得S=699.9996493Z=0.000145136X=19823.25398Y=28503.95084计算得S=699.9999985Z=5.000003137X=19785.25749Y=28575.02270计算得S=780.0000035Z=-5.000001663X=19781.15561Y=28572.16358计算得S=780.0000025Z=-0.000002979X=19777.05373Y=28569.30446计算得S=780.0000016Z=4.99999578X=19747.536Y=28654.131计算得S=870.0001137Z=-4.

21、99941049X=19742.686Y=28652.914计算得S=870.0003175Z=-0.00041814X=19737.837Y=28651.697计算得S=870.0002748Z=4.999808656X=19741.5912Y=28722.0580计算得S=939.9999786Z=-5.123024937X=19736.4769Y=28722.3564计算得S=939.9999862Z=-0.000027710X=19733.4730Y=28722.5317计算得S=940.0000238Z=3.00898694point=10/point自己写的fx-4800小程序一、

22、 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。R=J。说明：“，”表示回车符；“。”表示终止符。三、ZBZS(坐标正算）SRXY,V=X+Rec(S,R),W=Y+J,N=V。E=W用全站仪测设公路中桩和边桩程序程序说明：本程序用于用全站仪测设中桩和边桩，采用坐标变换的方法编制而成，在使用fx-4500计算器时，可直接输入以下程序便可使用；而在使用fx-4800计算器时，见fx-4800程序清

23、单。注意：因程序的0(零)和O易混淆,所以字母O用较大的全角符号，零仍用标准字。程序输入时应连续输入，无空格。而且输入时文件名需和给出的一样，否则会出现调用文件错误。变量说明：X输入摆站点的x坐标C输入第一个曲线角桩x坐标I输入中间曲线角桩的x坐标Y输入摆站点的y坐标D输入第一个曲线角桩y坐标J输入中间曲线角桩的y坐标E最后一个曲线角桩x坐标R中间曲线半径L中间曲线缓和曲线长F最后一个曲线角桩y坐标K第一个K为中间角桩桩号S中间角桩曲线总长，检验用K第二个K为要测设的桩号B测中桩时输入0，测边桩时输入宽度（左输负值，右输正值）程序清单注释文件:F1文件名：1L1L2L3Defm8XYCDIJE

24、FRLK: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:Fixm: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:LG

25、oto3L8Q=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+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+360T1

26、80T=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“R

27、S”: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)22V:SNTD=O+(SN)GU(TS+N)22V:D=O+(SN)GD“H=”L4:RE=:M=Goto1:R0Goto3:Goto4L7:Lbl4:SZI=0.02:Goto5:SZ+CI=0.02(0.02EPX)C(SZ):Goto5:SHCI=EPX:Goto5:SHI=EPX+(0.02EPX)C(SH+C):Goto5:I=0.02:Goto5L8:Lbl3

28、:Z1=0.04(EPX+0.02)C:SZ+Z1I=0.02:Goto5:SZ+CI=0.02+(EPX0.02)(CZ1)(SZZ1):Goto5:SHCI=EPX:Goto5:SProgZJFX=2Prog“ZBF”X=3Prog“TTX=4Prog“SJDX=5Prog“P2LX=6Prog“B2AGoto1ZBF“-ZXZBjs-”:YD:Y“I1”“I2”YD=Z=0rog“FP”:A=B:B=0:Goto3X=0:TPXVWU:T“HJD(X+Yi)”Lbl1X“XJA=1”“JD(X+Yi)”:U“JDLC”:Y=Arg(P-T)YY=Y+360Y“I”C“L”=Abs(P-T

29、)X=0=V=P:Goto2V“JDB(X+Yi)”=Arg(V-P)DD=D+360D“I”Z“L”=Abs(V-P)ProgFP:A=BLbl2W“QJD(X+Yi)”=Arg(W-V)DD=D+360D“I”C“L”=Abs(W-V)Lbl4Prog“FP”X=1=B=B-AX=0=A=B:B=0:Z=0:E=RepP:F=ImpP:Goto3M=Rep(P-T):N=Imp(P-T):I=Rep(W-V):J=Imp(W-V):F“Y”=(MImpTN-IImpVJ+RepV-RepT)(MN-IJ)E“X”=M(F-ImpT)N+RepTC“L1”=(E-RepP)2+(F-ImpP

30、)2)U=U+C:C“L2”=(E-RepV)2+(F-ImpV)2)Lbl3X=0:X:X“ZBF(0)ZJF(1)”X=1=O=ZrogZJFA=A+BrogZBjsZBjsI=0:J=0:B=0:RGH:R:G“LS1”:H“LS2”rog“TWE”Z=0=E“JD(X)”:F“JD(Y)”:U“JDLC”Rec(-Z7,Y):Z15=I+E:Z16=J+F:Rec(Z8,D):Z13=I+E:Z14=J+F:C=U-Z7:Z6“ZH”=CProg“ZDZ”C=Z6Lbl1C:C“JIAZHUANG”:Z=C-Z6Z0=Rec(Z,Y):Z11“X”=I+Z15Z12“Y”=J+Z16P

31、rog“BZ”:Goto1ZG=S=G=Z1:Q=Z3rog“H”ol(Z11,Z12):Rec(I,Y-JO):Z11“X”=I+Z15Z12“Y”=J+Z16Prog“BZ”:Goto1ZW-H=S=G=Z1:Q=Z3rog“Y”ol(Z11,Z12):Rec(I,Y-JO):Z11“X”=I+Z15Z12“Y”=J+Z16ProgBZ:Goto1ZW=Z=W-Z=H=Z2:Q=Z4rogHol(Z11,Z12):Rec(I,D+180+JO):Z11“X=I+Z13Z12“Y=J+Z14Z=W-Zrog“BZ”:Goto1ZW=Z=Z-W:Rec(Z,D):Z11“X=I+Z13Z12“

32、Y=J+Z14Z=W+Zrog“BZ”:Goto1FPB=Y-DB180=B=B-360BB=B+360O=BAbsB:B=AbsBTWES=Grog“PQ”:Z1=P:Z3=Q=Hrog“PQ”:Z2=P:Z4=QZ7“T1”=(R+Z2-(R+Z1)cos(A+B)sin(A+B)+Z3Z8“T2”=(R+Z1-(R+Z2)cos(A+B)sin(A+B)+Z4W“L”=R(A+B)180+.5G+.5HX“E”=(R+Z1)cos(tan-1(R+Z2)(R+Z1)-cos(A+B)sin(A+B)-RZ9=Z7-OsinBsin(A+B):Z10=Z8-OsinAsin(A+B)B0=Z9“TA”Z10“TB”