求5800隧道测量程序！！！

程序浏览:

FileName:RESET初始化程序

Norm 1:50→C:12345→J

"RESET PW"?I:I=J=>500→DimZ 为数据库增加额外变量500个，在SET、SETPFDYS程序根据实际再增减变量

50→Z[C+22]

"PASSWORDS"?I:I→Z[C+39] 重设要素保护密码

Cls:Stop

FileName:DATLOCK 要素保护密码确定认程序

Cls:Norm 1:50→C

"PASSWORDS"?J:Cls:J≠Z[C+39]=>Stop2008-08-07修改

FileName:SHELL(外壳程序)

50→C 在扩充变量预留前50个给别的程序用。如不够就适量加大。RESET，SHELL，SET，SETPFDYS，这几个程序中C值必需一致

Z[C+35]→I:"STATION-N"?I:I→Z[C+35] 设置测站N坐标

Z[C+36]→I:"STATION-E"?I:I→Z[C+36] 设置测站E坐标

Z[C+44]→I:"STATION-Z"?I:I→Z[C+44]　设置测站Z坐标

Z[C+45]→I:"STATION-HI"?I:I→Z[C+45] 设置仪高

Z[C+41]→I:"GC-DH"?I:I→Z[C+41] 放样点高差常数

Norm 1

Z[C+23]→N　当前分段要素N坐标

Z[C+24]→E　当前分段要素E坐标

Z[C+25]→M　当前分段要素起点桩号

Z[C+26]→H　当前分段要素起点方位角（正北）单位：弧度

Z[C+27]→A　当前分段要素起点曲率有左偏负右偏正（注意不半径）

Z[C+28]→R　当前分段要素终点曲率有左偏负右偏正（注意不半径）

Z[C+29]→L　当前分段要素长度

Z[C+3]→W

Z[C+32]→I:"DAT1 2 3"?I:I→Z[C+32] 平曲线要素数库类型选择1为内置式，2文件式，3实时输入(查看当前要素值)

0→I

Z[C+4]→I:"1 2 3"?I:I→Z[C+4] 选择本程序模式默认为坐标正算，1坐标反算，2横向边仰坡放样，3，纵向边仰坡放样(隧道进出口用到)

Lbi 1

Z[C+4]→I

I=1=>Prog"TURNZH"

I=2=>Prog"HXBYP"

I=3=>Prog"ZXBYP"

Lbi A

Deg:Norm 1:Cls

Z[C+1]→G:"ZH"?G:G→Z[C+1] G桩求桩号

Z[C+2]→B:"JL"?B:B→Z[C+2]　 　B横向距离，左正右负

Lbi B

Prog "ZBJS" 坐标正算

Prog "GCJS" 高程计算

Fix 3:Cls 设置三位小数

"ZH＝":Locate 4,1,G　 第一行显示桩号

"X＝":Locate 3,2,X 第二行显示X(N)坐标

"Y＝":Locate 3,3,Y；第三行显示Y(E)坐标

Prog "GCJS" 高程计算

"Z":Locate 2,4,Z+Z[C+41] ；第四行显示Z坐标

Locate 10,4,Z[C+2] ；第四行显示横向距离

0→I

Do:I+1→I:I=1000=>Goto 1:LpWhile Not(Getkey=57 Or Getkey=27) 锁定键盘，并在几十秒后自动返回

Getkey=57=>Goto 1

Z+Z[C+41]+Z[C+49]-Z[C+44]-Z[C+45]→Z[C+48]

Pol(Z[C+35]-X,Z[C+36]-Y)

Cls:"DL":Locate 3,1,I:Locate 12,1,Z[C+45]

"DH":Locate 3,2,Z[C+48]:Locate 12,2,Z[C+49]

"FWJ"

J<0=>J+360→J:J◤DMS◢ 2008-08-07 修改

J→Z[C+46]

I→Z[C+47]

Goto 1

FileName：ZBJS 坐标计算程序

Prog"READDAT"

Rad

G-M→Q

IF AR=0 :Then If A=R :Then 1→J:Else 3→J :IfEnd:Else If A=R :Then 2→J:Else 3→J:IfEnd:IfEnd

If J=1 :Then H→F:N+QCos(H)→X:E+QSin(H)→Y:IfEnd 直线段直接计算

If J=2 :Then H+QR→F:Rec(Abs(2Sin(Abs(0.5QR))÷R),H+QR÷2):N+I→X:E+J→Y:IfEnd 圆弧段直接计算

If J=3 :Then Goto 5 :IfEnd

X+BCos(F+W)→X

Y+BSin(F+W)→Y

Deg:Return

Lbi 5 用五点通用坐标计算计算缓和段

0.5(R-A)÷L→K

AQ→I

KQ2→J

0.0469100770→P:H+IP+JP2→U

0.2307653449→P:H+IP+JP2→V

0.5→P :H+IP+JP2→D

0.7692346551→P:H+IP+JP2→F

0.9530899230→P:H+IP+JP2→T

0.1184634425→I

0.2393143352→J

0.2844444444→O

N+Q(ICos(U)+JCos(V)+OCos(D)+JCos(F)+ICos(T))→X

E+Q(ISin(U)+JSin(V)+OSin(D)+JSin(F)+ISin(T))→Y

H+AQ+KQ2→F

X+BCos(F+W)→X

Y+BSin(F+W)→Y

Deg:Return

FileName:TURNZH(坐标反算)

Norm 1

Z[C+37]→I:"DQD-N"?I:I→Z[C+37] 输入待求点N坐标

I=-1=>Prog"INFWJDLDH"

Z[C+38]→I:"DQD-E"?I:I→Z[C+38] 输入待求点E坐标

0→B：M+L÷2→G：Prog"ZBJS"

Lbi S:Rad

Z[C+37]-X→I:Z[C+38]-Y→J

If I=0 And J=0 :Then Goto A:IfEnd

Pol(I,J)

Lbi A:Rec(I,J-F):J→Z[C+2]

G+I→G:IF Abs(I)>0.0001 :Then Prog"ZBJS":Goto S↙

G→Z[C+1]

FileName：INFWJDLDH　实测坐标（用方位角，距离，高差）输入程序

Deg:Norm 1:Cls

Z[C+46]→J:"DQD-FWJ"?J:J→Z[C+46]

Z[C+47]→I:"DQD-DL"?I:I→Z[C+47]

Z[C+48]→K:"DQD-DZ"?K:K→Z[C+48]

Z[C+49]→P:"RHT"?P:P→Z[C+49]

Rec(I,J)

Z[C+35]+I→Z[C+37]

Z[C+36]+J→Z[C+38]

Z[C+44]+Z[C+45]+K-P→Z[C+43]

FileName:HXBYP 横向边仰坡放样程序

Cls:Norm 1

Z[C+40]→I:"QPD-DL"?I:I→Z[C+40]输入起坡点与中桩距离常数

Z[C+41]→I:"QPD-DZ"?I:I→Z[C+41] 输入起坡点与中桩高差常数

Z[C+42]→I:"i"?I:I→Z[C+42]输入边仰坡坡度，左仰坡(路堑)为正，右仰坡(路堑)为负，左边坡为负，右边坡为正

Cls:Z[C+43]→I:"DQD-Z"?I:I→Z[C+43] 输入实测高程

Prog"TURNZH"

Prog "GCJS"

(Z[C+43]-Z-Z[C+41])Z[C+42]+Z[C+40]→B

Fix 3:B-Z[C+2]◢ 显示与上一次测点的偏差

B→Z[C+2]

FileName:ZXBYP 纵向仰坡放样程序(隧道进口使用)

Cls:Norm 1:Z[C+33]→I:"QPD-ZH"?I:I→Z[C+33]输入起坡点桩号

Z[C+34]→I:"QPD-Z"?I:I→Z[C+34] 输入起坡点高程

Z[C+42]→I:"i"?I:I→Z[C+42]仰坡时(隧道进出口仰坡)，进口为正，出口为负)

Cls:Z[C+43]→I:"DQD-Z"?I:I→Z[C+43] 输入实测高程

Prog"TURNZH"

Z[C+42](Z[C+43]-Z[C+34])+Z[C+33]→G

Fix 3:G-Z[C+1]◢ 显示与上一次测点的偏差

G→Z[C+1]

FileName:SET 设置程序

1\Lbi S

2\Norm 1

3\50→C

4\50→Z[C+22]

5\0→K

6\Z[C+5]→I:"1 Or 2"?I:I→Z[C+5] 选择1线或者2线

7\180Z[C+3]÷π→J:"XZJJ"?J:Jπ÷180→Z[C+3]斜桩夹角

8\If I=1:Then Z[C+8]→Z[C+6]:Z[C+9]→Z[C+7]

9\Z[C+12]→Z[C+16]:Z[C+13]→Z[C+17]:Z[C+18]→Z[C+20]

10\Else Z[C+10]→Z[C+6]:Z[C+11]→Z[C+7]

11\Z[C+14]→Z[C+16]:Z[C+15]→Z[C+17]:Z[C+19]→Z[C+20]

12\IfEnd

13\Z[C+7]+3→Z[C+31]

14\0→I:"1PQX,2SQX,3YSDS"?I 1输入平曲线要素，2输入竖曲线要素，3输入1、2线夹平竖曲线交点(变坡点)数量以让程序规划内存建立要素数据为，其它值退出本程序

15\I≠0=>Prog"DATLOCK" 2008-07-28 18：23修改

16\I=1 =>Goto 1:I=2 =>Goto 2:I=3 =>Goto 3 2008-07-28 18：23修改

17\Prog"SETPFDYS"

18\Stop

19\Lbi 1　 平曲线要素输入部分

20\Z[C+6]→D

21\-1→Z[D+3]

22\For 1→J To Z[C+16]

23\Cls:"PQX":Locate 9,1,J

24\Z[D+1]→I:"N"?I:I→Z[D+1] 请输入N坐标

25\Z[D+2]→I:"E"?I:I→Z[D+2] 请输入E坐标

26\Z[D+3]→I:"R"?I:I→Z[D+3] 请输入半径R

27\Z[D+4]→I:"LS1"?I:I→Z[D+4] 请输入缓和曲线1长度

28\Z[D+5]→I:"LS2"?I:I→Z[D+5]请输入缓和曲线2长度

29\D+5→D

30\Next

31\1→Z[D-2]

32\Goto S

33\Lbi 2 竖曲线要素输入部分

34\Z[C+7]→D

35\-1→Z[D+3]

36\For 1→J To Z[C+17]

37\Cls:"SQX":Locate 9,1,J

38\Z[D+1]→I:"BPDZH"?I:I→Z[D+1] 请输入变坡点桩号2008-08-07修改

39\Z[D+2]→I:"BPDGC"?I:I→Z[D+2] 请输入变坡点高程 2008-08-07修改

40\Z[D+3]→I:"R"?I:I→Z[D+3] 请输入半径R 2008-08-07修改

41\D+3→D

42\Next

43\1→Z[D]

44\Goto S

45\Lbi 3 要素点数输入部分

46\Z[C+12]→I:"1PQXDS"?I:I→Z[C+12] 输入1线平曲线点数,最小值为3

47\Z[C+13]→I:"1SQXDS"?I:I→Z[C+13] 输入1线竖曲线点数,最小值为3

48\Z[C+18]→I:"1JDZH"?I:I→Z[C+18] 输入1线平曲线第一个交点的交点桩号

49\Z[C+14]→I:"2PQXDS"?I:I→Z[C+14] 输入2线平曲线点数,最小值为3

50\Z[C+15]→I:"2SQXDS"?I:I→Z[C+15] 输入2线竖曲线点数,最小值为3

51\Z[C+19]→I:"2JDZH"?I:I→Z[C+19] 输入2线平曲线第一个交点的交点桩号

52\C+Z[C+22]→Z[C+8] 1线路平曲线要素指针偏移基数

53\Z[C+8]+5Z[C+12]→Z[C+9] 1线路竖曲线要素指针偏移基数

54\Z[C+9]+3Z[C+13]→Z[C+10] 2线路平曲线要素指针偏移基数

55\Z[C+10]+5Z[C+14]→Z[C+11] 2线路竖曲线要素指针偏移基数

56\Z[C+11]+3Z[C+15]→Z[C+21] 要素数据库结束指针

57\Z[C+21]+1→Z[C+30]

58\Z[C+21]→DimZ 增加额外变量

59\-1→K

60\Goto S

正算主程序(ZS)FIX4：?S：?Z：Prog “PM-SJ”： Abs(S-O) → W：Prog "SUB1"： "XS="：X◢

"YS="：Y◢F-90→F(需要时可以让他显示)：Prog“SQX-SJ”: Prog“SQX”：“H=”：H◢

反算主程序(FS)FIX4: ?S：?X：？Y：Prog“PM-SJ”： X→ I： Y→ J：Prog "SUB2"： "S="：O+W→S◢"Z="：Z◢

Prog“SQX-SJ”：Prog“SQX” ：“H=”：H◢边坡放样主程序（BP-FY）Lb1 0：Prog“ZS”：“H-BG”（中桩与坡脚起算点高差值，比中桩高正，反之负）?A：H+A →B：?P（实测点高程）：?L（坡脚起算点到中桩的距离）：0.75（挖方时一级坡度） →C：1 （挖方时二级坡度）→D：1.5 （挖方时三级坡度）→E：8（挖方时一级坡高） →G：10（挖方时二级坡高） →M：15（挖方时三级坡高） →N：2 （平台宽度）→K：1 （填方时一级坡度）→I：1.5（填方时二级坡度） →J：2 （填方时三级坡度）→O：2 （填方时一级坡高）→Q：8 （填方时二级坡高）→R：10 （填方时三级坡高）→T：ifP＞B：thenGoto1：ElseGoto2Lb1 1：ifP＞B：thenL+C（P-B）→U：P-B →F：IfEnd←┘ifP＞B+G：thenL+CG+K+D（P-B-G）→U：P-B-G →F：IfEnd←┘ifP＞B+G+M：thenL+GC+2k+MD+E（P-B-G-M）→U：P-B –G-M→F：IfEnd：Goto3←┘Lb1 2：ifP≤B：thenL+I（B -P）→U： B -P →F：IfEnd←┘ifP≤B-Q：thenL+IQ+K+J（B-Q-P）→U：B-Q-P→F：IfEnd←┘ifP≤B-Q-R：thenL+IQ+2k+JR+O（B-Q-R-P）→U：B-Q-R-P→F：IfEnd：Goto3←┘Lb1 3：U-AbsZ→V：ifZ＜0：thenZ-V→Z：Else z+v→z：IfEnd←┘ “Z=”：Z◢计算得出正确的宽度，路线左为负，右为正。“V=”V◢ 宽度调整值，正为远离中线，负为向中线靠近。“PT=”：F◢距离即将到达的填挖方平台的高度。注：由于路基的宽度和超高情况较为复杂，因此该程序的缺点是不能自动计算路基宽度和路基边缘坡脚起算点的设计高程，需要手算坡脚起算点与中桩的宽度及与中桩的高差。算例；S? 700（放样点桩号） Z？-15（试算宽度） X=19839.5817 Y=28515.5010 H=107.144 H-BG？ -0.15（边缘起算点的高程比中桩低0.15） p?150（实测高程）L? 12（起算边缘点距中桩的距离） Z=-69.509（计算得出正确的宽度）V=54.509（远离中线方向距测点54.509米），PT=25.006（距离平台高度）。 然后再次进入下一轮试算，直到V趋近于零。

隧道3心圆放样主程序（SD-3XY）Lb1 0：Prog“FS”：1.5→A：1.65→B：0.89→G：2.27→C：3.3→D：2.41→E：0.51→J：“H1”：F：ifF≤H+A：then E-AbsZ→W：IfEnd：ifF＞H+A And F≤H+A+B：then√（D�0�5-（F-H-A）�0�5）-G-AbsZ→W：IfEnd：ifF＞H+A+B：then√（C�0�5-（F-H-A-J）�0�5）-AbsZ→W：IfEnd：“W=”：W◢ （水平方向的超欠挖：正为欠,负为超）Goto0正算子程序(SUB1)

1÷P→ C： (P-R)÷(2HPR) → D： 180÷π→ E： 0.1739274226→ A： 0.3260725774→ B： 0.0694318442→ K： 0.3300094782→ L： 1-L→ F：1-K→ M： U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+MWD))) → X： V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))) → Y： G+QEW(C+WD)+90→ F： X+ZcosF→ X： Y+ZsinF→ Y反算子程序(SUB2)

G-90→T： Abs((Y-V)cosT-(X-U)sin（T）) → W： 0→ Z：Lbl 0：Prog "SUB1"： T+QEW(C+WD) →L： (J-Y)cosL-(I-X)sinL → Z：ifAbsZ<1E-6：thenGoto1：ElssW+Z → W：Goto 0：IfEnd←┘

Lbl 1： 0→ Z：Prog "SUB1"： (J-Y)÷sinF→ Z子程序（平面数据库）PM-SJifS ≥500（线元起点里程）AndS＜769.256（线元止点里程）：Then19942.837 →U（线元起点X坐标）：28343.561 →V（线元起点Y坐标）：500 →O（线元起点里程）：125°16°31°→G（线元起点方位角）：269.265 →H（线元长度）：10^（45）→P（线元起点曲率半径）：10^（45）→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左负右正）：IfEnd←┘ifS ≥769.256（线元起点里程）AndS＜806.748（线元止点里程）：Then19787.34 →U（线元起点X坐标）：28563.378→V（线元起点Y坐标）：769.256→O（线元起点里程）：125°16°31°→G（线元起点方位角）：37.492→H（线元长度）：10^（45）→P（线元起点曲率半径）：221.75→R（线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘ifS ≥806.748（线元起点里程）AndS＜919.527（线元止点里程）：Then19766.566→U（线元起点X坐标）：28584.574→V（线元起点Y坐标）：806.748→O（线元起点里程）：120°25°54.07°→G（线元起点方位角）：112.779→H（线元长度）：221.75→P（线元起点曲率半径）：221.75→R（线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘ifS ≥919.527（线元起点里程）AndS＜999.812（线元止点里程）：Then19736.072 →U（线元起点X坐标）：28701.893 →V（线元起点Y坐标）：919.527→O（线元起点里程）：91°17°30.63°→G（线元起点方位角）：80.285→H（线元长度）：221.75→P（线元起点曲率半径）：9579.228→R（线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘ifS ≥999.812（线元起点里程）：Then19744.038 →U（线元起点X坐标）：28781.659→V（线元起点Y坐标）：999.812→O（线元起点里程）：80°40°50°→G（线元起点方位角）：100→H（线元长度）：10^（45）→P（线元起点曲率半径）：10^（45）→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左负右正）：IfEnd←┘子程序（竖曲线数据库）SQX-SJifS≤999.812（竖曲线终点里程）：then0.0357（前坡度：上坡正下坡负）→E：0.0227（后坡度：上坡正下坡负）：600→R（竖曲线半径）：780→G（变坡点里程）：110→C（变坡点高程）：IfEnd←┘子程序（竖曲线计算公式）SQXE-F→J：ifJ＞0：then-R→R：ElssR→R：Abs（R*J÷2）→T：ifS≤G-T：thenG-S→L：C-LE→H：IfEnd←┘ifS＞G-T And S≤G： thenG-S→L：（S+T-G）�0�5÷2÷R→N：C-LE+N → H：IfEnd←┘ifS＞G And S≤G+T： then S-G→L：（G +T- S）�0�5÷2÷R→N：C-LF+N → H：IfEnd←┘ifS＞G+T ：thenS-G→L： C+LF → H：IfEnd←┘ 三、使用说明

1、规定

(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，

Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。

(2) 当所求点位于中线时，Z=0；当位于中线左铡时，Z取负值；当位于中线中线右

侧时，Z取正值。

(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次代替。

(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆

弧的半径。

(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45

次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半

径为无穷大，以10的45次代替；与圆曲线相接时，曲率半径等于圆曲线的半径。

(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的

值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等

于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。2、输入与显示说明

S ？ 正算时所求点的里程，反算时估算点里程，反算时尽量不要超出真实测点所在的线元里程终点，即尽量向小里程估计，如超出线元终点计算器将死机。

Z ？正算时所求点距中线的边距(左侧取负，值右侧取正值，在中线上取零)

X ？反算时所求点的X坐标

Y ？反算时所求点的Y坐标

显示部分：

XS=×××正算时，计算得出的所求点的X坐标

YS=×××正算时，计算得出的所求点的Y坐标

S=×××反算时，计算得出的所求点的里程

Z=×××反算时，计算得出的所求点的边距H=××× 正反算时，计算得出所求点的中桩设计高程（如弯道沿边线超高时要另外考虑超高段中桩高程）四、算例

某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组

成，各段线元的要素（起点里程S0、起点坐标X0 Y0、起点切线方位角F0、线元长度

LS、起点曲率半径R0、止点曲率半径RN、线 元左右偏标志Q）如下：

S0 X0 Y0 F0LS R0 RN Q

500.00019942.83728343.561 125 16 31.00269.256 1E45 1E45 0

769.25619787.34028563.378 125 16 31.00 37.492 1E45 221.75 -1

806.74819766.56628594.574 120 25 54.07 112.779221.75 221.75 -1

919.52719736.07228701.893 91 17 30.63 80.285221.75 9579.228 -1

999.81219744.03828781.659 80 40 50.00 100.0001E45 1E45 0

1、正算

（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元

要素）

S=700 Z=-5计算得 XS=19831.41785 YS=28509.72590 H=107.144

S=700 Z=0 计算得 XS=19827.33592 YS=28506.83837 H=107.144

S=700 Z= 5计算得 XS=19823.25398 YS=28503.95084 H=107.144 S=780 Z=-5计算得 XS=19785.25749 YS=28575.02270 H=112.759

S=780 Z=0 计算得 XS=19781.15561 YS=28572.16358 H=112.759

S=780 Z= 5计算得 XS=19777.05373 YS=28569.30446 H=112.759

S=870 Z=-5计算得 XS=19747.53609 YS=28654.13091 H=112.759

S=870 Z=0 计算得 XS=19742.68648 YS=28652.91379 H=112.759

S=870 Z= 5计算得 XS=19737.83688 YS=28651.69668 H=112.759

S=940 Z=-5.123计算得 XS=19741. 59118 YS=28722.05802 H=130.475

S=940 Z=0计算得 XS=19736.47687 YS=28722.35642H=130.475

S=940 Z= 3.009计算得 XS=19733.47298 YS=28722.53168 H=130.475

2、 反算

X=19831.418 Y=28509.726 计算得 S=699.9999974 Z= -5 .00018164

X=19827.336 Y=28506.838 计算得 S=699.9996493 Z= 0.000145136H=107.144

X=19823.25398 Y=28503.95084计算得 S=699.9999985 Z= 5.000003137

X=19785.25749 Y=28575.02270 计算得 S=780.0000035 Z= -5 .000001663

X=19781.15561 Y=28572.16358 计算得 S=780.0000025 Z=- 0.000002979H=112.759

X=19777.05373 Y=28569.30446 计算得 S=780.0000016 Z= 4.99999578

X=19747.536 Y=28654.131 计算得 S=870.0001137 Z= -4.99941049

X=19742.686 Y=28652.914 计算得 S=870.0003175 Z=- 0.00041814 H=130.475

X=19737.837 Y=28651.697 计算得 S=870.0002748 Z= 4.999808656 X=19741.5912 Y=28722.0580 计算得 S=939.9999786 Z= -5.123024937

X=19736.4769 Y=28722.3564 计算得 S=939.9999862 Z=- 0.000027710 H=130.475

X=19733.4730 Y=28722.5317 计算得 S=940.0000238 Z= 3.00898694

---