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

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-x

2、ody=yp-yoP点转换后之工程坐标为xp、yp:xp=Xo+dx*COS(a)-dy*SIN(a)yp=Yo+dx*SIN(a)+dy*COS(a) 坐标方位角计算程序置镜点坐标：ZX  ZY后视点坐标：HX  HY方  位  角：W两点间距离: SLb1 0A,  B,  C,  DAZX=:BZY=:CHX=:DHY=:W=tg1(D-B)÷(C-A):(D-B)>0=>(C-A)>0=&g

3、t;W=W:(D-B)>0=>(C-A)<0=>W=W+180:(D-B)<0=>(C-A)<0=>W=W+180:(D-B)<0=>(C-A)>0=>W=360+WW=WS=(D-B)2+(C-A)2) Goto 0CASIO  fx4500p坐标计算程序根据坐标计算方位角WW360W：“ALF(12)”ÞL1  A“X1”：B“Y1”：Pol(C“X2”A，D“Y2”B：“S”W0直线段坐标计算L1   X“X(0)”：

4、Y“Y(0)”：S“S(0)”：A“ALF”L2   Lb1 2L3   L：L“LX”L4   M“X(Z)”X(LS)cosAL5   N“Y(Z)”Y(LS)sinAL6   B：B“B(L)”：Q“Q”L7   O“X(L)”MBcos(AQ180)L8   P“Y(L)”NBsin(AQ180)L9   C：C“B(R)”L10 

5、; U“X(R)”MCcos(AQ)L11  V“Y(R)”NCsin(AQ)L12  Goto  2园曲线段坐标计算L1   S“S(0)-Km”：X“X(0)”：Y“Y(0)”：A“ALF”：R“R”：K“K(L1,R2)”L2   Lb1 2L3   L：L“L(X)”L4   V180×（LS）R：WV2L5   CA(-1)K×W：D2R

6、sinW：FA(-1)K×VL6   M“X(Z)”XDcosCL7   N“Y(Z)”YDsinCL8   E：E“B(L)”：Q“Q”L9   O“X(L)”MEcos(FQ180)L10  P“Y(L)”NEsin(FQ180)L11  G：G“B(R)”L12  T“X(R)”MGcos(FQ)L13  U“Y(R)”NGsin(FQ)L14  Goto

7、60; 2正向缓和曲线段坐标计算L1   S“ZH-Km”：X“X(ZH)”：Y“Y(ZH)”：A“ALF”：R“R”：H“LS”：K“K(L1,R2)”L2   Lb1 2L3   L：L“L(X)”L4   D30（LS）2RH：CLS（LS）590（R×H）2：BAD(-1)K：EA3D(-1)KL5   U“X(Z)”XCcosBL6   V“Y(Z)”YCsinBL7 

8、;  G：G“B(L)”：Q“Q”L8   F“X(L)”UGcos(EQ180)L9   I“Y(L)”VGsin(EQ180)L10  J：J“B(R)”L11  M“X(R)”UJcos(EQ)L12  N“Y(R)”Vjsin(EQ)L13  Goto  2卵形曲线坐标计算X1,D2)”ÞL1   S“Km-YH”：E“X(YH)”：F“Y(YH)”：G“ALF”：B“

9、R1”：D“A”：K“K(L1,R2)”：Q“R1-R2L2   Lb1 2L3   Z：Z“L(X)”L4   J“L1”D2B： R“RP”D2B(D2(-1)Q(ZS)B)：L“LP”D2RL5   M(LJ)(L5J5)40D4(L9J9)3456D8L6   N(L3J3)6D2(L7J7)336D6(L11J11)42240D10L7   TG(-1)Q(-1)K×J2

10、15;90D2L8   X“X(Z)”E(-1)QMcosT(-1)KNsinTL9   Y“Y(Z)”F(-1)QMsinT(-1)KNcosTL10  A“ALF(P)”G(-1)K(ZS)×90×(1B1R)L11  H：H“B(L)”：U“Q”L12  W“X(L)”XHcos(AU180)L13  V“Y(L)”YHsin(AU180)L14  C：C“B(R)”L15  I“X(R)”

11、XCcos(AU)L16  P“Y(R)”YCsin(AU)L17  Goto  2公路逐桩坐标计算4800程序公路逐桩坐标计算程序(可以计算对称、不对称缓和曲线)Lb1 0Z=.  V=.W=V+2：FixmKLb1 1K>ZW+5Z+4=>W=W+1:Goto 1    (判断桩号在哪个交点范围，就是该交点曲线起点至下一交点曲线起点)S=K-ZW+5Z+3      &#

12、160;                       （计算该桩号与曲线起点的距离） 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）     &#

13、160; （计算该交点与下一交点直线方位角）J<0=>J=J+360A=JPol（ZW-1-ZW-2,ZW+Z+1-ZW+Z）      （计算该交点与上一交点直线方位角）J<0=>J=J+360C=A-J:A=J                        &#

14、160;   （计算偏角）W=V+2=>Goto2（如果桩号在起点与第一交点曲线起点之间，则转 Lb1 2 ）I=Abs（tan（c÷2）M=L÷2-L3÷240R2:N=E÷2-E3÷240R2P=L2÷6R-L4÷336R3-R（1-cos（90L÷R）Q=E2÷6R-E4÷336R3-R（1-cos（90E÷R）D=（P-Q）I÷2 : F=（P+Q+2R）I÷2M=F+M-

15、D:Q=F+N+DN=RAbsC÷180+（L+E）÷2X=ZW-1-McosAY=ZW+Z+1-MsinAM=ZW-1+Qcos（A+C）V=ZW+Z+1+Qsin（A+C）       Q=AbsC÷CSL=>P=0:Goto3（如果桩号在第一缓和曲线内，则转 Lb1 3）SN-E=>S=S-L:Goto4（如果桩号在圆曲线内，则转 Lb1 4）SN=>S=N-SQ=-Q:A=A+C-180:X=M:Y=V:L=E：P=180:Got

16、o3 （如果桩号在第二缓和曲线内，则转 Lb1 3）P=A+C:S=S-N:D=M+ScosP:F=V+SsinPGoto6                                     &

17、#160;     （如果桩号在直线内，则转 Lb1 6）Lb1 2P=A+CD=ZW-1+ScosPF=ZW+Z+1+SsinP:Goto6Lb1 3I=S-S5÷40R2÷L2+S9÷3456R4÷L4J=Q（S3÷6RL-S7÷336R3÷L3）P=P+A+90QS2÷RL:Goto5Lb1 4M=90（2S+L）÷RI=RsinM+L÷2-L3÷240R2J=Q（L2÷

18、24R+R（1-cosM）P=A+QMLb1 5D=X+IcosA-jsinA:F=Y+JcosA+IsinA Lb1 6DX=（结果显示X坐标）FY=（结果显示Y坐标）PAT=（结果显示该桩号方位角）BO：BSO      （输入边桩距离，交角）P=P+OLXB=D+BcosP（结果显示边桩X坐标）MYB=F+BsinP（结果显示边桩Y坐标）以上是坐标计算程序，括号内是程序计算的大致原理及说明，中间部分为直线、圆曲线、缓和曲线计算的各种公式，大家也知道，书上也有。   &#

19、160; 该程序是最基本的，如再加几条语句就可以处理断链。    该程序还需变量输入，把整条线路曲线要素一次输入到4800中，就可以计算线路任意中、边桩坐标。如有需要，需RMB ￥200，附有变量输入方法及详细使用说明。联系：   曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序一、程序功能        本程序由一个主程序(

20、TYQXjs)和两个子程正算子程序(SUB1)、反算子程序(SUB2)序构成，可以根据曲线段直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。另外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运行。二、源程序1.主程序(TYQXjs)"1.SZ =>

21、0;XY"："2.XY => SZ"：N：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"RN"：Q：C=1÷P：D=(P-R)÷(2HPR)：E=180÷：N=1=>Goto 1：>Goto 2Lbl 1：SZ：SZ：W=Abs(S-O)：Prog "SUB1"：X&qu

22、ot;XS"=XY"YS"=YGoto 1Lbl 2：XY：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+WZ"Z"=ZGoto 22. 正算子程序(SUB1)A=0.1739274226：B=0.3260725774：K=0.0694318442：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+Q

23、EMW(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：Lbl 0：Prog "SUB1"：L=T+QEW(C+WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：>W=W+Z：Goto

24、0;0Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF三、使用说明1、规定        (1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。        (2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右侧时，Z取正值。 &#

25、160;      (3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。        (4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。         (5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的

26、半径。止点与直线相接时，曲率半径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。         (6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明       输入部分：     

27、 1. SZ => XY      2. XY = > SZ       N "  选择计算方式，输入1表示进行由里程、边距计算坐标 ；输入2表示由坐标反算              里程

28、和边距。     X0 .线元起点的X坐标      Y0 .线元起点的Y坐标      S0 .线元起点里程     F0 .线元起点切线方位角     LS .线元长度     R0 .线元起点曲率半径 &

29、#160;   RN .线元止点曲率半径     Q .  线 元左右偏标志(左偏Q=-1，右偏Q=1，直线段Q=0）      S . 正算时所求点的里程      Z .正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)      X

30、0;.反算时所求点的X坐标      Y .反算时所求点的Y坐标       显示部分：       XS=×××    正算时，计算得出的所求点的X坐标       YS=×××   &#

31、160;正算时，计算得出的所求点的Y坐标       S=×××    反算时，计算得出的所求点的里程       Z=×××    反算时，计算得出的所求点的边距四、算例        某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完

32、整缓和曲线+直线）组成，各段线元的要素（起点里程S0、起点坐标X0 Y0、起点切线方位角F0、线元长度LS、起点曲率半径R0、止点曲率半径RN、线 元左右偏标志Q）如下：S0   X0       Y0                   F0    

33、;        LS R0 RN     Q500.000    19942.837    28343.561     125 16 31.00    269.256     1E45 1E45&#

34、160;     0769.256    19787.340    28563.378     125 16 31.00      37.492     1E45 221.75     -1806.748  &

35、#160; 19766.566    28594.574     120 25 54.07     112.779    221.75 221.75     -1919.527    19736.072    28701.893  

36、;     91 17 30.63       80.285    221.75 9579.228   -1999.812    19744.038    28781.659       80 40 50.00

37、     100.000    1E45 1E45       0         1、正算          （注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）    &#

38、160;       S=700  Z=-5    计算得  XS=19831.41785   YS=28509.72590            S=700  Z=0     计算得  XS=19827.33

39、592   YS=28506.83837            S=700  Z= 5    计算得  XS=19823.25398   YS=28503.95084            S=780 

40、; Z=-5    计算得  XS=19785.25749   YS=28575.02270            S=780  Z=0     计算得  XS=19781.15561   YS=28572.16358   

41、60;        S=780  Z= 5    计算得  XS=19777.05373   YS=28569.30446                       

42、S=870  Z=-5    计算得  XS=19747.53609   YS=28654.13091            S=870  Z=0     计算得  XS=19742.68648   YS=28652.91379  

43、;          S=870  Z= 5    计算得  XS=19737.83688   YS=28651.69668                   S=940 

44、60;Z=-5.123    计算得  XS=19741. 59118  YS=28722.05802            S=940  Z=0            计算得  XS=19736.47687  

45、; YS=28722.35642            S=940  Z= 3.009    计算得  XS=19733.47298   YS=28722.53168               &#

46、160;   2、 反算             X=19831.418   Y=28509.726   计算得  S=699.9999974  Z= -5 .00018164           &

47、#160; X=19827.336   Y=28506.838   计算得  S=699.9996493  Z= 0.000145136             X=19823.25398   Y=28503.95084    计算得  S=699.99999

48、85  Z= 5.000003137           X=19785.25749   Y=28575.02270  计算得  S=780.0000035  Z= -5 .000001663           X=19781.15561

49、   Y=28572.16358   计算得  S=780.0000025  Z=- 0.000002979           X=19777.05373   Y=28569.30446   计算得  S=780.0000016  Z= 4.99999578 

50、;                    X=19747.536   Y=28654.131   计算得  S=870.0001137  Z= -4.99941049          &#

51、160;X=19742.686   Y=28652.914   计算得  S=870.0003175  Z=- 0.00041814           X=19737.837   Y=28651.697   计算得  S=870.0002748  Z= 4.99980

52、8656          X=19741.5912   Y=28722.0580   计算得  S=939.9999786  Z= -5.123024937          X=19736.4769   Y=28722.3564  &#

53、160;计算得  S=939.9999862  Z=- 0.000027710          X=19733.4730   Y=28722.5317   计算得  S=940.0000238  Z= 3.00898694 point=10/point自己写的fx-4800小程序一、 ZBFS(坐标反算)NEXY ,

54、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用全站仪测设公路中桩和边桩程序程序说明： 本程序用于用全

55、站仪测设中桩和边桩，采用坐标变换的方法编制而成，在使用fx-4500计算器时，可直接输入以下程序便可使用；而在使用fx-4800计算器时，见fx-4800程序清单。 注意：因程序的0(零)和O易混淆,所以字母O用较大的全角符号，零仍用标准字。程序输入时应连续输入，无空格。而且输入时文件名需和给出的一样，否则会出现调用文件错误。变量说明： <程序运算符定义>X 输入摆站点的x坐标 C 输入第一个曲线角桩x坐标 I 输入中间曲线角桩的x坐标Y 输入摆站点的y坐标 D 输入第一个曲线角桩y

56、坐标 J 输入中间曲线角桩的y坐标 E 最后一个曲线角桩x坐标 R 中间曲线半径 L 中间曲线缓和曲线长 F 最后一个曲线角桩y坐标 K 第一个K为中间角桩桩号 S 中间角桩曲线总长，检验用 K 第二个K为要测设的桩号 B 测中桩时输入0，测边桩时输入宽度（左输负值，右输正值） 程 序 清 单注 释文件:F1 文件名：1 L1 L2

57、60;L3Defm 8 XYCDIJEFRLK:Fixm:Z2=0:Prog B Z5=0.5L-L Xy 3/240R2:Z3=Z5+(R+L2/24R)Abs tan T2:Z=K-Ans:S=R Abs T/180+L:SH=Z+Ans:G=RL:Prog 2文件:F2 文件名：2 L1 Lbl 0:Fixm:K L2 Lbl 6:B:Norm:K(-1) Xy Z2>(H-L)(-1

58、) Xy Z2>Goto 4Frac Z22=0=>Z6=Z:>Z6=HA=90T/Abs T L3 B<0>N=0:>N=1B=Abs B L4 =Abs(K-Z6:K(-1) Xy Z2>Z6(-1) Xy Z2>Goto 2Q=Z3+:Z4=B L5B0>Prog CProg D:P=0 L6 Prog E:Goto 

59、0 L7 Lbl 2:L>Goto 3 L8 Q=Z3-+ Xy 5/40G2:Z4= Xy 3/6G- Xy 7/336G Xy 3:Prog D:P=B L9 B0>Prog FProg E:Goto 0 L10 Lbl 3:Q=Z3-Z5:Z4=R+L2/24R:Prog D L11 Z8=U+A(2(-L)/R-1+L

60、/R):B0>Prog G:>P=RProg E:Goto 0 L12Lbl 4:M=C:=D:C=E:D=F:E=M:F=:Z2=Z2+1:Prog B:Goto 6文件:F3 文件名：0 L1Pol(M,):W<0>W=W+360文件:F4 文件名：B L1 L2M=E-I:=F-J:Prog 0:T=W:M=I-C:=J-D:Prog 0:U=W:T=T-W T<-180>T=T+360T>180&

61、gt;T=T-360文件:F5 文件名：C L1 A=-90(-1) Xy (N+Z2)文件:F6 文件名：D L1 Z7=I-QcosU+Z4cos(U+A):Z1=J-QsinU+Z4sin(U+A)文件:F7 文件名：E L1 L2M=Z7+PcosZ8-X:=Z1+PsinZ8-Y:Prog 0 WFix 3:VNorm文件:F8 文件名：F L1Z8=U+A+A2/G+90+A(-1) Xy (N+Z2)文

62、件:F9 文件名：G L1P=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: Lbl 6 : SL3: S+TN<0 D=O+(SN)F: SN<0

63、60;D=O+(SN)FU(SN+T)2÷2÷V: SNT D=O+(SN)GU(TS+N)2÷2÷V: D=O+(SN)G D“ H=”  L4: R< E= : M= Goto 1: R< E= : M= Goto 1:L5: Lbl 1 : M  L C=L: Goto&

64、#160;2: C=ML6: Lbl 2 : P X >0 Goto 3 : Goto 4L7: Lbl 4 : S Z I=0.02 : Goto 5: S<Z+C I=0.02(0.02EPX)÷C×(SZ) : Goto 5: S<HC I=EPX : G

65、oto 5: S<H I=EPX+(0.02EPX)÷C×(SH+C) : Goto 5: I=0.02 : Goto 5L8: Lbl 3 : Z1=0.04÷(EPX+0.02)×C : S Z + Z1 I=0.02 : Goto 5: S<Z+C I=0.02+(EPX0.02)÷

66、;(CZ1)×(SZZ1) : Goto 5: S<HC I=EPX : Goto 5: S<HZ1 I=EPX(EPX0.02)÷(CZ1)×(SH+C) : Goto 5: I=0.02 : Goto 5L9: Lbl 5 : B : W“HB”=DIBL10: Goto 6说明： R:平曲线

67、半径 L“LS”:缓和曲线长 P：偏转，Z=1,Y=1X“XL”：线路（Z,Y）左线1,右线1 Z“ZH”: 直缓点里程H“HZ” ：缓直点里程 V“RS”：竖曲线半径 F“I 1” ：前坡堵（带符号）G“I2”： 后坡度值（带符号） O“BPH”：变坡点高程 N“BP”： 变坡点里程T： 切线长 E：最大超高值 C ：渐变段长度（LC） I： 横坡值 M：最大渐变长度（Lmax） B ：距中心距离 W“HB”：边部高程   公路中线坐标计算程序 CASIO 4500公路中线坐标计算程序 casio 4500ZXCLCX Defm16 Lbl1 X