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

大地(高斯平面)坐标系工程坐标系转换大地坐标系-工程坐标系=待转换点为P,大地坐标为：Xp、Yp工程坐标系原点o：大地坐标：Xo、Yo工程坐标：xo、yo工程坐标系x轴之大地方位角：adX=Xp-XodY=YpYoP点转换后之工程坐标为xp、yp:xp=dXCOS(a）+dYSIN（a）+xoyp=dX*SIN(a)+dYCOS(a）+yo工程坐标系大地坐标系=待转换点为P,工程坐标为：xp、yp工程坐标系原点o:大地坐标：Xo、Yo工程坐标:xo、yo工程坐标系x轴之大地方位角：adx=xpxody=ypyoP点转换后之工程坐标为xp、yp:xp=Xo+dxCOS（a)-dySIN（a）yp=Yo+dxSIN（a)+dyCOS（a)坐标方位角计算程序置镜点坐标：ZXZY后视点坐标：HXHY方位角：W两点间距离：SLb10A，B，C，DAZX=:BZY=：CHX=：DHY=：W=tg1（D-B）(C-A)）:(DB)0=（CA）0=W=W：(D-B）0=（C-A)W=W+180：(DB)0=（CA)W=360+WW=WS=（(DB)2+（CA)2)Goto0CASIOfx4500p坐标计算程序根据坐标计算方位角WW360W：“ALF(12）L1A“X1”：B“Y1”:Pol（C“X2”A，D“Y2B：“S”W0直线段坐标计算L1X“X（0）”：Y“Y（0)”：S“S（0)”:A“ALF”L2Lb12L3L：L“LXL4M“X（Z)”X(LS）cosAL5N“Y（Z)Y(LS)sinAL6B：B“B（L）:Q“QL7O“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：WV2L5CA（1)KW：D2RsinW：FA（1）KVL6M“X（Z)”XDcosCL7N“Y（Z)”YDsinCL8E：E“B（L):Q“QL9O“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“ZHKm”：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（-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“KmYH: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)(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（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（ZWZW1，ZW+Z+2ZW+Z+1)（计算该交点与下一交点直线方位角）JGoto2（如果桩号在起点与第一交点曲线起点之间，则转Lb12）I=Abs（tan（c2)）M=L2-L3240R2：N=E2E3240R2P=L26R-L4336R3-R（1-cos(90LR）)Q=E26RE4336R3-R(1-cos(90ER)）D=（P-Q)I2：F=(P+Q+2R）I2M=F+MD：Q=F+N+DN=RAbsC180+（L+E）2X=ZW-1-McosAY=ZW+Z+1MsinAM=ZW1+Qcos（A+C）V=ZW+Z+1+Qsin（A+C）Q=AbsCCSL=P=0:Goto3（如果桩号在第一缓和曲线内,则转Lb13）SNE=S=SL：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=SN:D=M+ScosP:F=V+SsinPGoto6(如果桩号在直线内,则转Lb16）Lb12P=A+CD=ZW1+ScosPF=ZW+Z+1+SsinP：Goto6Lb13I=SS540R2L2+S93456R4L4J=Q（S36RL-S7336R3L3）P=P+A+90QS2RL：Goto5Lb14M=90（2S+L）RI=RsinM+L2L3240R2J=Q(L224R+R（1-cosM)P=A+QMLb15D=X+IcosAjsinA：F=Y+JcosA+IsinALb16DX=(结果显示X坐标）FY=（结果显示Y坐标）PAT=（结果显示该桩号方位角）BO：BSO(输入边桩距离，交角）P=P+OLXB=D+BcosP(结果显示边桩X坐标)MYB=F+BsinP（结果显示边桩Y坐标）以上是坐标计算程序,括号内是程序计算的大致原理及说明，中间部分为直线、圆曲线、缓和曲线计算的各种公式,大家也知道,书上也有。该程序是最基本的,如再加几条语句就可以处理断链。该程序还需变量输入，把整条线路曲线要素一次输入到4800中,就可以计算线路任意中、边桩坐标。如有需要,需RMB￥200，附有变量输入方法及详细使用说明。联系：13858423737曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器）程序曲线任意里程中边桩坐标正反算（CASIOfx4800P计算器)程序一、程序功能本程序由一个主程序（TYQXjs）和两个子程-正算子程序（SUB1)、反算子程序（SUB2)序构成，可以根据曲线段-直线、圆曲线、缓和曲线（完整或非完整型)的线元要素(起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIOfx4500P计算器及CASIOfx-4850P计算器上运行。二、源程序1.主程序（TYQXjs）”1。SZ=XY”：”2。XY=SZ”：N：U”X0”：VY0”：O”S0：GF0”：H”LS：P”R0：R”RN”：Q：C=1P:D=（PR）（2HPR)：E=180：N=1=Goto1：Goto2Lbl1：SZ：SZ：W=Abs（SO）：ProgSUB1：X”XS”=XYYS=YGoto1Lbl2：XY:XY:I=X:J=Y：Prog”SUB2”:SS”=O+WZ”Z=ZGoto22。正算子程序（SUB1）A=0。1739274226：B=0.3260725774：K=0.0694318442：L=0。3300094782:F=1-L：M=1K：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=G90:W=Abs(（YV）cosT-（X-U)sinT)：Z=0:Lbl0:Prog”SUB1：L=T+QEW（C+WD）：Z=（JY)cosL-（IX）sinL:AbsZ1E-6=Goto1：W=W+Z:Goto0Lbl1：Z=0:ProgSUB1：Z=(JY）sinF三、使用说明1、规定(1）以道路中线的前进方向（即里程增大的方向）区分左右;当线元往左偏时，Q=1；当线元往右偏时,Q=1；当线元为直线时，Q=0。(2）当所求点位于中线时，Z=0；当位于中线左铡时,Z取负值；当位于中线中线右侧时，Z取正值.（3）当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。（4)当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。（5)当线元为完整缓和曲线时,起点与直线相接时,曲率半径为无穷大，以10的45次代替；与圆曲线相接时,曲率半径等于圆曲线的半径.止点与直线相接时，曲率半径为无穷大，以10的45次代替;与圆曲线相接时,曲率半径等于圆曲线的半径。(6)当线元为非完整缓和曲线时,起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时,曲率半径等于圆曲线的半径.止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明输入部分：1.SZ=XY2。XY=SZN?选择计算方式,输入1表示进行由里程、边距计算坐标;输入2表示由坐标反算里程和边距.X0?线元起点的X坐标Y0？线元起点的Y坐标S0?线元起点里程F0？线元起点切线方位角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、起点曲率半径R0、止点曲率半径RN、线元左右偏标志Q)如下：S0X0Y0F0LSR0RNQ500。00019942.83728343.5611251631。00269。2561E451E450769.25619787.34028563.3781251631。0037.4921E45221。751806.74819766。56628594。5741202554.07112。779221.75221。75-1919.52719736。07228701。893911730。6380。285221.759579.2281999。81219744。03828781。659804050。00100。0001E451E4501、正算（注意:略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素)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=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=28506.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。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自己写的fx4800小程序一、 ZBFS(坐标反算）NEXY，Fixm.F,Pol（N-X,EY).J0J=J。J0J=J+360.说明:“，”表示回车符；“。表示终止符。二、ZBFS2（坐标反算）NEXY，Pol(NX，EY)。J0J=J。J0J=J+360.S=I。R=J.说明：“,表示回车符；“。表示终止符。三、ZBZS(坐标正算)SRXY，V=X+Rec（S，R），W=Y+J,N=V.E=W用全站仪测设公路中桩和边桩程序程序说明:本程序用于用全站仪测设中桩和边桩,采用坐标变换的方法编制而成,在使用fx4500计算器时，可直接输入以下程序便可使用;而在使用fx4800计算器时，见fx4800程序清单.注意：因程序的0(零)和O易混淆，所以字母O用较大的全角符号，零仍用标准字。程序输入时应连续输入，无空格。而且输入时文件名需和给出的一样，否则会出现调用文件错误.变量说明：Goto4FracZ22=0=Z6=Z：Z6=HA=90T/AbsTL3BN=1B=AbsBL4=Abs(KZ6：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=Z3Z5：Z4=R+L2/24R：ProgDL11Z8=U+A(2(L)/R1+L/R）：B0ProgG:P=RProgE：Goto0L12Lbl4：M=C:=D:C=E：D=F：E=M：F=:Z2=Z2+1：ProgB:Goto6文件：F3文件名：0L1Pol(M,）：W0W=W+360文件:F4文件名：BL1L2M=EI：=FJ：Prog0:T=W:M=IC：=J-D:Prog0：U=W：T=T-WTT=T-360文件：F5文件名：CL1A=90（1）Xy(N+Z2）文件：F6文件名：DL1Z7=IQcosU+Z4cos(U+A):Z1=JQsinU+Z4sin（U+A）文件：F7文件名：EL1L2M=Z7+PcosZ8X：=Z1+PsinZ8Y：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）22V：SNTD=O+（SN）GU(TS+N）22V:D=O+(SN）GD“H=L4：RE=:M=Goto1：RE=：M=Goto1：L5：Lbl1:MLC=L:Goto2：C=ML6:Lbl2：PX0Goto3：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：Z1=0。04（EPX+0。02）C：SZ+Z1I=0。02：Goto5：SZ+CI=0。02+（EPX0。02）(CZ1）(SZZ1）：Goto5：SHCI=EPX:Goto5:SHZ1I=EPX（EPX0.02)（CZ1)(SH+C)：Goto5：I=0。02:Goto5L9：Lbl5：B:W“HB”=DIBL10：Goto6说明:R:平曲线半径L“LS”:缓和曲线长P：偏转，Z=1，Y=1X“XL”：线路(Z,Y）左线1，右线1Z“ZH”：直缓点里程H“HZ”：缓直点里程V“RS：竖曲线半径F“I1”:前坡堵（带符号）G“I2”:后坡度值（带符号）O“BPH”：变坡点高程N“BP：变坡点里程T：切线长E：最大超高值C：渐变段长度(LC)I：横坡值M:最大渐变长度（Lmax）B：距中心距离W“HB”：边部高程公路中线坐标计算程序 CASIO 4500公路中线坐标计算程序casio4500ZXCLCXDefm16Lbl1X=1：X:X“ZJF1ZBF2TT-3SJD4P2L5B2A6”X=1Prog”ZJFX=2Prog“ZBF”X=3Prog“TT”X=4Prog“SJDX=5Prog“P2L”X=6Prog“B2A”Goto1ZBF“-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（PT)Y0=Y=Y+360Y“I”C“L”=Abs（PT)X=0=V=P：Goto2V“JDB（X+Yi)=Arg（VP)D0=D=D+360D“IZ“L”=Abs（V-P）ProgFP”:A=BLbl2W“QJD(X+Yi）”=Arg（W-V）DD=D+360D“I”C“L”=Abs(WV）Lbl4Prog“FPX=1=B=BAX=0=A=B：B=0：Z=0:E=RepP:F=ImpP:Goto3M=Rep（PT）:N=Imp（P-T）:I=Rep(WV）:J=Imp（W-V）：F“Y”=（MImpTN-IImpVJ+RepV-RepT)（MNIJ）E“X”=M（F-ImpT）N+RepTC“L1”=（(ERepP)2+(FImpP）2)U=U+C：C“L2”=（E-RepV）2+（FImpV）2)Lbl3X=0:X：X“ZBF（0）ZJF(1）X=1=O=ZrogZJF”A=A+BrogZBjs”ZBjsI=0：J=0:B=0:RGH:R：G“LS1”：H“LS2”rog“TWE”Z=0=E“JD(X)”:F“JD（Y）”:U“JDLCRec（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=CZ6Z0=Rec（Z，Y）：Z11“X=I+Z15Z12“Y”=J+Z16Prog“BZ”：Goto1ZG=S=G=Z1:Q=Z3rog“H”ol(Z11，Z12)：Rec（I，YJO）:Z11“X”=I+Z15Z12“Y=J+Z16Prog“BZ：Goto1ZW-H=S=G=Z1：Q=Z3rog“Y”ol(Z11，Z12):Rec（I,YJO）：Z11“X”=I+Z15Z12“Y=J+Z16Prog”BZ”:Goto1ZW=Z=WZ=H=Z2:Q=Z4rogHol（Z11，Z12）：Rec(I,D+180+JO）：Z11“X”=I+Z13Z12“Y”=J+Z14Z=WZrog“BZ:Goto1ZW=Z=ZW：Rec(Z,D）:Z11“X”=I+Z13Z12“Y”=J+Z14Z=W+Zrog“BZ”：Goto1FPB=YDB180=B=B360BB=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(tan1(（R+Z2）(R+Z1)cos（A+B）)sin(A+B)）)RZ9=Z7-OsinBsin(A+B）：Z10=Z8OsinAsin（A+B)B0=Z9“TAZ10“TB