全站仪的坐标放样,那个坐标怎么计算出来的啊？.doc

坐标要已知才能放样呀，如果要计算坐标，可以用CAsio4800编程计算，只要有公式就可以自己编入计算器运用，当然你可直接上网下载

如果是公路的我整理的你可以参考
CASIO4800程序组
1、极坐标法放样
Prog：FY
Lb1 0：A“X0”：B“Y0”：I=0：J=0：Pol（（C“XA”-A），（D“YA”-B）：J<0=>G“FW- OA”=J+360▲L“L0”=I▲Goto 1：≠> G“FW O-A”=J▲L“L0”=I▲
Lb1 1：{EQ}：E“Xi”：Q“Yi”：Pol（（E-A），（Q-B））：J<0=>J=J+360：Goto 2：≠> Goto 2
Lb1 2： F“FW-OB”=J▲L=I▲0=F-G：O<0=>O“BJ”=O+360▲Goto 3：≠> O “BJ” ▲
Lb1 3：P=O-180▲Goto 1
注：a、输入：（X0、Y0）、（XA、YA）——测站点坐标、后视点坐标
Xi、Yi ——放样点坐标
b、输出：FW-OA——测站至后视边方位角、L0——后视边长
FW-OB——测站至放样点方位角、L——放样边长
BJ——后视边置零，放样点顺时针拨角
P——偏角（+为右偏、-为左偏）{本值用于计算路线偏角}
2、公路竖曲线高程计算程序
Prog:SQX
Lbl A:A“+（-）i1”:B“+（-）i2” W=（B-A）÷100：R：T=Abs（RW）÷2：L=T*2：E=T2÷（2R）：K“JD K+”：G“JD H”：C=K-T：D=K+T：
Lbl 0：J“Ki+”：J<0=>Goto 1：≠> Goto 2△△
Lb1 1：“Out QX1”：H=G-(K-J)A÷100▲Goto 5
Lb1 2：J>D=>Goto 4 △W<0=>F=-1△W>0=>F=1△J>K=>Goto 3△H=G-（K-J）A÷100+F(J-C)2÷（2R）▲Goto 5△
Lb1 3：H=G+（J-K）B÷100+F（D-J）2÷（2R）▲Goto 5△
Lb1 4：“OUT QX2”：H=G+（J-K）B÷100▲Goto 5△
Lb1 5：M“DHi”：H=H+M▲
注：a、公式：L=|R(i2-i1)| 、T=L÷2、E=T2÷（2R）、h=l2÷(2R)
b、功能：已知前后坡度%、竖曲线半径，计算各桩高程。
c、输入：i1、i2——前、后坡度%（上坡+、下坡-）
JD K+、 JD H、R——变坡点桩号、高程、竖曲线半径
Ki+——计算桩号
DHi——计算点与设计线高差（高+、低-）
d、输出：H——待计算点高程

3、直线坐标计算
ProgZX
Lb1 0：M“K0”：A“X1”：B“Y1”：C“X2”：D“Y2”：I=0：J=0：
Pol（（C-A），（D-B））：L=I▲J<0=>F”FW”=360+J▲Goto 1：≠>F”FW”=J▲
Lb1 1：{KZQ}:K“Ki”：Z“Lb”：Q=90：Q“ANG”：L=K-M：Rec（L，F）：X=A+I：Y=B+J：Rec（Z，F+Q）：O=I：P=J：X=X+O▲Y=Y+P▲ Goto 1
注：
a、公式：X=X0+L*cos(Fw)、Y=Y0+L*sin(Fw)
b、功能：已知直线起终点坐标，计算中间桩中边桩坐标
c、输入： K0、（X1、Y1）、（X2、Y2）——起点桩号、起点坐标、终点坐标
Ki、Lb、Ang —— 待求桩号、边桩距（左-、右+）、边桩与中线右夹角（90）
d、输出：L、FW——直线长度、起点至终点方位角
X、Y——计算点坐标
a多边形面积计算（一）
Prog：Area
Lb1 0：S=0：L=0：E“X0”：F“Y0”：A=E：B=F
Lb1 1：{CD}：C“Xi： Xi＝X0，Yi=Y0 ＝＞END”：D“Yi”：M=√（（C-E）2+（D-F）2）：S=S+（AD-CB）：L=L+√（（C-A）2+（D-B）2）：A=C：B=D：S=0＝＞Goto 1：
Lb1 3：M＞0=＞Goto 1：≠＞S“Area”=Abs S÷2▲L▲
注：a、公式：Area =（∑（XiYi+1-Xi+1Yi））÷2
b、功能：已知多边形各转点坐标，计算多边形面积。
c、输入：X0、Y0——起点坐标
Xi、Yi——各点坐标（按顺序输入、最后要闭合即最后一点要输入起点坐标）
d、输出：Area——面积
L——周长

4、对称标准曲线中边桩坐标计算
Prog：QX
Lb1 A: M“JDK+”:A“PJ=>R（+），L（-）”:R：W“LS”: A≥0=>G=1:Goto E: ≠>G=-1:
Lb1 E:Defm 10:A=Abs A:B=90W÷R÷π：P=W2÷（6R）-W4÷（336R3）-R（1-cosB）:Q=W-W3÷(40R2)-RsinB：T=（R+P）tan（A÷2）+Q：L=RAπ÷180+W:E=(R+P)÷cos(A÷2)-R:N=M-T:Z[5]=N+W:Z=N+L÷2:C=N+L-W:D=N+L:A=AG:
{U}:U“ZD-2,XY-1”=2=>T▲E▲L▲N“ZH”▲Z[5]“HY”▲Z“QZ”▲C“YH”▲D“HZ”▲Goto A：
≠> Goto B
Lbl B：F“FW”：V“JD-X”：H“JD-Y”：I=0：J=0：Rec（T，F+180）：Z[1]“ZH X”=I+V：Z[2]“ZH Y”=J+H：Rec（T，F+A）：Z[3]“HZ X”=I+V：Z[4]“HZ Y”=J+H：L=0
Lbl 0：E=90：{KSE}：K“JS K+”：S“LB”：E“ANG”：S=0=>U=1：Goto C：≠>U=S÷（Abs S）
Lbl C：S=Abs S：W=0=> Goto D：≠>K Goto 1：≠>K Goto 2：≠>K Goto 4：≠> Goto 5Δ
Lbl D：KGoto 1：≠>KGoto 3：≠> Goto 5
Lbl 1：Rec（M-K，F+180）：X=V+I：Y=J+H：Rec（S，F+180-（180-E）U）：X“Xi”=X+I▲Y“Yi”=Y+J▲S=SU：Goto 7△
Lbl 2：Z=K-N：O=90Z2÷R÷W÷π：X=Z-Z5÷40÷R2÷W2+Z9÷3456÷R4÷W4：Z=Z3÷6÷R÷W-Z7÷336÷R3÷W3+Z11÷42240÷R5÷W5：
Lbl 6：Rec（X，F）：X=Z[1]+I：Y=Z[2]+J：Rec（Z，F+90G）：X=X+I：Y=Y+J：Rec（S，F+OG+EU）：X“Xi”=X+I▲Y“Yi”=Y+J▲S=SU：Goto 7
Lbl 3：Z=K-Z[5]：O=（K-N-W÷2）*180÷π÷R：X=Q+Rsin（B+180Z÷π÷R）：Goto 6
Lbl 4：Z=D-K：O=90Z2÷R÷W÷π：X=Z-Z5÷40÷R2÷W2+Z9÷3456÷R4÷W4：Z=Z3÷6÷R÷W-Z7÷336÷R3÷W3+Z11÷42240÷R5÷W5：
Rec（X，F+A+180）：X=Z[3]+I：Y=Z[4]+J：Rec（Z，F+A+180-90G）：X=X+I：Y=Y+J：Rec（S，F+A+180-OG-（180-E）U）：X“Xi”=X+I▲Y“Yi”=Y+J▲S=SU：Goto 7
Lbl 5：Rec（（K-D+T），（F+A））：X=V+I：Y=H+J：Rec（S，F+A+EU）：X“Xi”=X+I▲Y“Yi“=Y+J▲ S=SU： Goto 7
Lbl 7：Goto 0
注：输入：JDK+、PJ、R、LS、——交点桩号、路线偏角（右+、左-）、半径、缓和曲线长
“ZD-2,XY-1”——要显示主点要素输入2，否则进入坐标计算
FW、（JD X、JD Y）——起始边方位角（即ZH—JD方位角）、交点坐标
JS K+、LB、Ang——待计算桩号、边桩距（右+、左-）、边桩与路中线右夹角（默认90度）
输出：T、E、L、ZH、HY、YH、HZ、——曲线要素、 Xi,Yi——待求点坐标
5、公路任何曲线通用中边桩坐标正反算（复化辛卜生公式——也称积木法
Prog：HQX
Defm 4
V“V=1=>XY： 2=>Ki”
Lbl 0：{KLW}
Lbl 4：Q“one=0：1、2……” =0=> Goto 1△Q=1=>Prog“ZD”（线元要素数据文件名）△ Q=3=>Prog“要素表文件名” △……（注：本句是为了全线及多条线路计算用，如仅单个可省略本句）
Lbl 1：A"K1"：B"K2"：C"FW0"：D"R1"：E"R2"：F"X0"：G"Y0"
D≠0=>I=1÷D：≠=>I=D△E≠0=>J=1÷E：≠=>J=E△
AbsD+AbsE=0=>M=2：≠=>M=16△
V=2=>L=0：W=90△K“Ki”：LL-BZ”：W"ANG"：
N=0：Z[1]=0：Z[2]=0：Z[3]=0：Z[4]=0 进入坐标迭代计算
Lbl 2：N=N+1：H=2(K-A)÷M：R=NH÷2+A：R=C+180÷π×(I+(J-I)÷2÷(B-A)×(R-A))×(R-A)：Int(N÷2)=N÷2=>Z[1]=Z[1]+cosR：Z[2]=Z[2]+sinR：≠=>Z[3]=Z[3]+cosR：Z[4]=Z[4]+ sinR△N=M=>Goto 3：≠>Goto 2
Lbl 3：
X=F+H÷6×(cosC+4Z[3]+2Z[1]-cosR)+Lcos(R+W)：
Y=G+H÷6×(sinC+4Z[4]+2Z[2]-sinR)+Lsin(R+W)：
V=2=>Goto 6△
X"Xi="▲
Y"Yi="▲
R"FW"=R-360Intg(R÷360) ▲
Q=0=>{KLW}：Goto 1 ：≠> Goto 0
Lbl 6： Pol(T"XF"-X,U"YF"-Y) 进入桩号求算
O=Icos(J-R)：AbsO≤0.0001=>O"DL"▲K=K+O▲
O"LP"=Isin(J-R)▲
{TU}：Goto 6：≠=>K=K+O：L=0：Goto 4
线元要素数据文件格式：K>起点桩号=> A"K1"=？：B"K2"=？：C"FW0"=？：D"R1"=？：E"R2"=？：F"X0"=？：G"Y0"=？：K> B=>A=B：B=？……（下曲线要素）

注：a、功能：利用公路通用复化辛卜生积分公式进行中边桩坐标计算及已知坐标反算桩号。
b、输入：
V=1=>XY： 2=>Ki分别进入坐标计算、桩号反算
Q“one=0：0、2……”要素直接辒入按0，否则戉嵌套对应文仲要素表计算，嵌套几个文件自己需要定。
K1、K2 、FW0、R1、 R2-------曲线起点、终点里程、曲线起点方位角、曲线起点、终点半径 (左-、右+，0为直线)
X0、Y0-------曲线起点坐标
Ki------------曲线待求点里程
BL---求炩左右偏距(左-右+)
ANG-------- -求点的右斜交角
XF、YF-------求桩号的点坐标（反算程序）
输出：
X、Y-------,-曲线求得坐怇
FW-----------待求点的切线方位角
DL、K+O、LP分别为桩号误差、求得桩号、左右偏距(左-右+)（反算程序）
当曲线的设计半径较小时，为保证点位计算精度，M(即程序中n的2倍)的取值可适当的大些。M为偶数，直线时M=2即可，经计算M=16即可满足半径为60的小半径曲线精度。

6、 a多边形面积计算（一）
Prog：Area
Lb1 0：S=0：L0：E“X0”：F“Y0”：A=E：B=F
Lb1 1：{CD}：C“Xi： Xi＝X0，Yi=Y0 ＝＞END”：D“Yi”：M=√（（C-E）2+（D-F）2）：S=S+（AD-CB）：L=L+√（（C-A）2+（D-B）2）：A=C：B=D：S=0＝＞Goto 1：
Lb1 3：M＞0=＞Goto 1：≠＞S“Area”=Abs S÷2▲L▲
注：a、公式：Area =（∑（XiYi+1-Xi+1Yi））÷2
b、功能：已知多边形各转点坐标，计算多边形面积。
c、输入：X0、Y0——起点坐标
Xi、Yi——各点坐标（按顺序输入、最后要闭合即最后一点要输入起点坐标）
d、输出：Area——面积
L——周长
b、多边形面积计算（二）
Prog：Area
Lb1 0：S=0：L=0： A“X0”：B“Y0”：：C“X1”：D“Y1”: L=√（（C-A）2+（D-B）2）
Lb1 1：{EF}：E“Xi：Xi＝X0÷Yi=Y0＝＞END”: F“Yi”：L=L+√（（E-C）2+（F-D）2）:M=√（（E-A）2+（F-B）2）:S=S+(C(F-B)+E(B-D)+A(D-F))÷2：C=E:D=F : S=0＝＞Goto 1
Lb1 2：M＞0=＞Goto 1：≠＞S“Area”=Abs S▲ L▲
注：a、公式：Area =∑{（（Xi（Yi+1-Y0）+ Xi+1（（Y0- Yi）+ X0（Yi- Yi+1））÷2}
b、功能：已知多边形各转点坐标，计算多边形面积。
c、输入：（X0、Y0）、（X1、Y1）——起始边坐标
Xi、Yi——各点坐标（按顺序输入、最后要闭合即最后一点要输入起点坐标）
d、输出：Area——面积
L——周长
7、道路土方计算
a、土方（棱台法）
Prog：TFa
Lb1 0：T=0：W=0：F“K0+”：A“AT0”：B“AW0”
Lb1 1：{CDE}：E“Ki=-1＝＞end”：E=-1＝＞Goto 6：≠＞C：“ATi”：D“AWi”：L=E-F
Lb1 2：C>A： ＝＞Goto 3： ≠＞M=C÷A： Goto 4
Lb1 3：M=A÷C： Goto 4
Lb1 4：O=(A+C)L(1+√M÷(1+M)) ÷3：D>B＝＞M=B÷D： Goto 5： ≠＞M=D÷B
Lb1 5：P=(B+D)L(1+√M÷(1+M)) ÷3：T=T+O：W=W+P：A=C：B=D：F=E： Goto 1
Lb1 6：T“TF” ▲W“WF” ▲“end”
注：：a、公式：V =（A1+A2）L（1+√m÷(1+m)）÷3
m=A1÷A2(其中A1>A2)
b、功能：已知路线各桩填表挖断面积，计算路线填挖量。
c、输入：K0+、AT0： AW0——起点桩号、填方面积、挖方面积
Ki、ATi、AWi——各桩号、填方面积、挖方面积（结束运行计算结果Ki输入—1）
d、输出：TF、WF——填方、挖方体积
b、土方（平均断面法）
Prog：TFb
Lb1 0：T=0：W=0：F“K0+”：A“AT0”：B“AW0”
Lb1 1：{CDE}：E“Ki=-1＝＞end”：E=-1＝＞Goto 3：≠＞C：“ATi”：D“AWi”：L=E-F
Lb1 2：T=T+（C+A）L÷2：W=W+（D+B）L÷2：A=C：B=D：F=E： Goto 1
Lb1 3：T“TF” ▲W“WF” ▲“end”
注：：a、公式：V =（A1+A2）L÷2
b、功能：已知路线各桩填挖断面积，计算路线填挖量。
c、输入：K0+、AT0： AW0——起点桩号、填方面积、挖方面积
Ki、ATi、AWi——各桩号、填方面积、挖方面积（结束运行计算结果Ki输入—1）
d、输出：TF、WF——填方、挖方体积