大地坐标转换成施工坐标公式.docx
- 文档编号:3100012
- 上传时间:2022-11-17
- 格式:DOCX
- 页数:13
- 大小:26.40KB
大地坐标转换成施工坐标公式.docx
《大地坐标转换成施工坐标公式.docx》由会员分享,可在线阅读,更多相关《大地坐标转换成施工坐标公式.docx(13页珍藏版)》请在冰豆网上搜索。
大地坐标转换成施工坐标公式
大地(高斯平面)坐标系工程坐标系转换
大地坐标系—-—〉工程坐标系
========================
待转换点为P,大地坐标为:
Xp、Yp
工程坐标系原点o:
大地坐标:
Xo、Yo
工程坐标:
xo、yo
工程坐标系x轴之大地方位角:
a
dX=Xp—Xo
dY=Yp-Yo
P点转换后之工程坐标为xp、yp:
xp=dX*COS(a)+dY*SIN(a)+xo
yp=—dX*SIN(a)+dY*COS(a)+yo
工程坐标系-——>大地坐标系
========================
待转换点为P,工程坐标为:
xp、yp
工程坐标系原点o:
大地坐标:
Xo、Yo
工程坐标:
xo、yo
工程坐标系x轴之大地方位角:
a
dx=xp—xo
dy=yp-yo
P点转换后之工程坐标为xp、yp:
xp=Xo+dx*COS(a)-dy*SIN(a)
yp=Yo+dx*SIN(a)+dy*COS(a)
坐标方位角计算程序
置镜点坐标:
ZX ZY
后视点坐标:
HX HY
方 位 角:
W
两点间距离:
S
Lb1 0←
{A, B, C, D}←
A〝ZX=〞:
B〝ZY=〞:
C〝HX=〞:
D〝HY=〞:
W=tg1((D-B)÷(C-A)):
(D—B)>0=>(C—A)〉0=>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+W∟∟W=W◢
S=√((D-B)2+(C—A)2) ◢
Goto 0←
CASIO fx-4500p坐标计算程序
根据坐标计算方位角
W=W+360△W:
“ALF(1~2)="⇒L1 A“X1=”:
B“Y1=”:
Pol(C“X2"-A,D“Y2”-B:
“S=”▲W<0
直线段坐标计算
L1 X“X(0)”:
Y“Y(0)”:
S“S(0)”:
A“ALF"
L2 Lb1 2
L3 {L}:
L“LX”
L4 M“X(Z)”=X+(L-S)cosA▲
L5 N“Y(Z)"=Y+(L-S)sinA▲
L6 {B}:
B“B(L)”:
Q“Q”
L7 O“X(L)”=M+Bcos(A+Q+180)▲
L8 P“Y(L)"=N+Bsin(A+Q+180)▲
L9 {C}:
C“B(R)”
L10 U“X(R)”=M+Ccos(A+Q)▲
L11 V“Y(R)”=N+Csin(A+Q)▲
L12 Goto 2
园曲线段坐标计算
L1 S“S(0)-Km”:
X“X(0)”:
Y“Y(0)":
A“ALF”:
R“R”:
K“K(L=1,R=2)”
L2 Lb1 2
L3 {L}:
L“L(X)”
L4 V=180/π×(L-S)/R:
W=V/2
L5 C=A+(—1)K×W:
D=2RsinW:
F=A+(-1)K×V
L6 M“X(Z)”=X+DcosC▲
L7 N“Y(Z)”=Y+DsinC▲
L8 {E}:
E“B(L)”:
Q“Q”
L9 O“X(L)"=M+Ecos(F+Q+180)▲
L10 P“Y(L)”=N+Esin(F+Q+180)▲
L11 {G}:
G“B(R)"
L12 T“X(R)"=M+Gcos(F+Q)▲
L13 U“Y(R)”=N+Gsin(F+Q)▲
L14 Goto 2
正向缓和曲线段坐标计算
L1 S“ZH—Km”:
X“X(ZH)”:
Y“Y(ZH)”:
A“ALF”:
R“R":
H“LS”:
K“K(L=1,R=2)”
L2 Lb1 2
L3 {L}:
L“L(X)”
L4 D=30(L-S)2/π/R/H:
C=L-S-(L-S)5/90/(R×H)2:
B=A+D(-1)K:
E=A+3D(-1)K
L5 U“X(Z)”=X+CcosB▲
L6 V“Y(Z)”=Y+CsinB▲
L7 {G}:
G“B(L)":
Q“Q”
L8 F“X(L)”=U+Gcos(E+Q+180)▲
L9 I“Y(L)”=V+Gsin(E+Q+180)▲
L10 {J}:
J“B(R)"
L11 M“X(R)"=U+Jcos(E+Q)▲
L12 N“Y(R)"=V+jsin(E+Q)▲
L13 Goto 2
卵形曲线坐标计算
X=1,D=2)"⇒L1 S“Km—YH”:
E“X(YH)”:
F“Y(YH)”:
G“ALF”:
B“R1”:
D“A”:
K“K(L=1,R=2)”:
Q“R1—R2
L2 Lb1 2
L3 {Z}:
Z“L(X)”
L4 J“L1”=D2/B:
R“RP"=D2B/(D2+(-1)Q(Z-S)B):
L“LP"=D2/R
L5 M=(L-J)-(L5-J5)/40/D4+(L9-J9)/3456/D8
L6 N=(L3-J3)/6/D2-(L7-J7)/336/D6+(L11-J11)/42240/D10
L7 T=G-(—1)Q(-1)K×J2×90/D2/π
L8 X“X(Z)”=E+(-1)QMcosT-(—1)KNsinT▲
L9 Y“Y(Z)”=F+(-1)QMsinT+(—1)KNcosT▲
L10 A“ALF(P)"=G+(-1)K(Z-S)×90×(1/B+1/R)/π
L11 {H}:
H“B(L)”:
U“Q”
L12 W“X(L)”=X+Hcos(A+U+180)▲
L13 V“Y(L)"=Y+Hsin(A+U+180)▲
L14 {C}:
C“B(R)”
L15 I“X(R)"=X+Ccos(A+U)▲
L16 P“Y(R)”=Y+Csin(A+U)▲
L17 Goto 2
公路逐桩坐标计算4800程序
公路逐桩坐标计算程序
(可以计算对称、不对称缓和曲线)
Lb1 0
Z=?
V=?
W=V+2:
Fixm
{K}
Lb1 1
K>Z[W+5Z+4]=〉W=W+1:
Goto 1⊿ (判断桩号在哪个交点范围,就是该交点曲线起点至下一交点曲线起点)
S=K—Z[W+5Z+3] (计算该桩号与曲线起点的距离)
R=Z[W+2Z+2]:
L=Z[W+3Z+2]:
E=Z[W+4Z+2] (读取该交点曲线要素R、Ls1 、Ls2)
Pol(Z[W]—Z[W-1],Z[W+Z+2]—Z[W+Z+1]) (计算该交点与下一交点直线方位角)
J<0=〉J=J+360⊿
A=J
Pol(Z[W—1]-Z[W-2],Z[W+Z+1]-Z[W+Z]) (计算该交点与上一交点直线方位角)
J<0=>J=J+360⊿
C=A-J:
A=J (计算偏角)
W=V+2=〉Goto2⊿ (如果桩号在起点与第一交点曲线起点之间,则转 Lb1 2 )
I=Abs(tan(c÷2))
M=L÷2—L^3÷240R^2:
N=E÷2-E^3÷240R^2
P=L^2÷6R—L^4÷336R^3—R(1-cos(90L÷πR))
Q=E^2÷6R—E^4÷336R^3-R(1—cos(90E÷πR))
D=(P-Q)I÷2 :
F=(P+Q+2R)I÷2
M=F+M-D:
Q=F+N+D
N=πRAbsC÷180+(L+E)÷2
X=Z[W—1]—McosA
Y=Z[W+Z+1]—MsinA
M=Z[W-1]+Qcos(A+C)
V=Z[W+Z+1]+Qsin(A+C)
Q=AbsC÷C
S≤L=〉P=0:
Goto3⊿ (如果桩号在第一缓和曲线内,则转 Lb1 3)
S≤N-E=>S=S—L:
Goto4⊿(如果桩号在圆曲线内,则转 Lb1 4)
S≤N=>S=N—S
Q=-Q:
A=A+C—180:
X=M:
Y=V:
L=E:
P=180:
Goto3 ⊿ (如果桩号在第二缓和曲线内,则转 Lb1 3)
P=A+C:
S=S-N:
D=M+ScosP:
F=V+SsinP
Goto6 (如果桩号在直线内,则转 Lb1 6)
Lb1 2
P=A+C
D=Z[W—1]+ScosP
F=Z[W+Z+1]+SsinP:
Goto6
Lb1 3
I=S-S^5÷40R^2÷L^2+S^9÷3456R^4÷L^4
J=Q(S^3÷6RL—S^7÷336R^3÷L^3)
P=P+A+90QS^2÷πRL:
Goto5
Lb1 4
M=90(2S+L)÷πR
I=RsinM+L÷2-L^3÷240R^2
J=Q(L^2÷24R+R(1-cosM))
P=A+QM
Lb1 5
D=X+IcosA—jsinA:
F=Y+JcosA+IsinA
Lb1 6
D″X=″◢ (结果显示X坐标)
F″Y=″◢ (结果显示Y坐标)
P″AT=″◢ (结果显示该桩号方位角)
{BO}:
B″S″O″⊿″ (输入边桩距离,交角)
P=P+O
L″XB″=D+BcosP◢ (结果显示边桩X坐标)
M″YB″=F+BsinP◢ (结果显示边桩Y坐标)
以上是坐标计算程序,括号内是程序计算的大致原理及说明,中间部分为直线、圆曲线、缓和曲线计算的各种公式,大家也知道,书上也有.
该程序是最基本的,如再加几条语句就可以处理断链。
该程序还需变量输入,把整条线路曲线要素一次输入到4800中,就可以计算线路任意中、边桩坐标。
如有需要,需RMB ¥200,附有变量输入方法及详细使用说明。
联系:
13858423737
曲线任意里程中边桩坐标正反算(CASIOfx—4800P计算器)程序
曲线任意里程中边桩坐标正反算(CASIO fx—4800P计算器)程序
一、程序功能
本程序由一个主程序(TYQXjs)和两个子程——正算子程序(SUB1)、反算子程序(
SUB2)序构成,可以根据曲线段——直线、圆曲线、缓和曲线(完整或非完整型)的线
元要素(起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲
率半径)及里程边距或坐标,对该曲线段范围内任意里程中边桩坐标进行正反算。
另
外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中,用于对曲
线任意里程中边桩坐标进行正反算。
本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运行.
二、源程序
1.主程序(TYQXjs)
"1。
SZ => 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 2Δ←┘
Lbl 1:
{SZ}:
SZ:
W=Abs(S—O):
Prog "SUB1":
X”XS"=X◢
Y"YS"=Y◢
Goto 1←┘
Lbl 2:
{XY}:
XY:
I=X:
J=Y:
Prog "SUB2":
S"S”=O+W◢
Z”Z”=Z◢
Goto 2
2。
正算子程序(SUB1)
A=0.1739274226:
B=0。
326072577
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 大地 坐标 转换 施工 公式