坐标反算里程公式,excel坐标反算里程公式( 二 )


(注:如有多个曲线元要素,还要继续添加到数据库DAT中)
I= 线元起点的X坐标:S= 线元起点的Y坐标:O= 线元起点里程:
G= 线元起点切线方位角:H= 线元终点里程P= 线元起点曲率半径
R= 线元终点曲率半径Q= 线元左右偏标志(注: 左偏为-1,右偏为+1 )
(注:如有多个曲线元要素,还要继续添加到数据库DAT中,曲率半径直径输入半径值)
5.坐标反算程序名: ZBFS
LBI 0:“X1=”?X:“Y1=”?Y:“X2=”?A:“Y2=”?B◢
POL(A-X,B-Y):J〈0
=〉J+360→J◢
“I=”:I◢
“J=”:J▲DMS ◢
Goto 0 ↙(注:↙表示按EXE键即可)
6. 任意多边形的面积程序名:RYDBXS
0→S :?A :?B :?C :?D↙(注:↙表示按EXE键即可)
Lbl 1:“X”?→X : “Y” ?→Y ↙(注:↙表示按EXE键即可)
0→I:0→J : P0l(C-A , D-B) : I→G : J→H :
P0l(X-A , Y-B) : X→C : Y→D : I→K : J→N :
“M=”:0.5GKsin(Abs(N-H))→M ↙(注:↙表示按EXE键即可)
“S=”:S+M→S ▲平方米
“U=”:S÷666.667→U ▲亩
Goto 1
(注:0表示数字零)
说明:点位必须按顺序输入成封闭形图型!
A B C D 为第一,二两点坐标(常量),X Y……为第三,四,五,六点坐标(变量) 。
例:X1=10.000, Y1=20.000
X2=500.000 , Y2=600.000
X3=700.000 , Y3=400.000三角形面积 107000M2160.50亩
X4=800.000 , Y4=300.000四边形面积 160500M2240.75亩
X5=650.000 , Y5=100.000五边形面积 218500M2327.75亩
X6=550.000 , Y6=50.000六边形面积 230500M2345.75亩
说明:
一、程序功能及原理
1.功能说明:
本程序由两个主程序——正算主程序(GSZS)、反算主程序(GSFS)和两个子程——正算子程序(SUB1)、线元数据库(DAT-M)构成,可以根据曲线段——直线、圆曲线、缓和曲线(完整或非完整型)的线元要素(起点坐标、起点里程、起点切线方位角、终点里程、起点曲率半径、止点曲率半径)及里程边距或坐标,对该曲线段范围内任意里程中边桩坐标进行正反算 。本程序可以在CASIO fx-4800P计算器及 CASIO fx-4850P计算器上运行 。由于加入了数据库(DAT-M),可实现坐标正反算的全线贯通 。
组合程序5可实现M线的正算贯通,组合程序7可实现M线的反算贯通,组合程序6可实现坐标计算到放样一体化 。
2.计算原理:
利用Gauss-Legendre 5点通用公式计算线路中边桩坐标并计算放样数据 。
利用待求点至线元起点切线作垂线,逐次迭代趋近原理反算里程及边距 。
二、使用说明
1、规定
(1) 以道路中线的前进方向(即里程增大的方向)区分左右;当线元往左偏时,Q= -1;当线元往右偏时,Q=1;当线元为直线时,Q=0 。
(2) 当所求点位于中线时,L=0;当位于中线左侧时,L取负值;当位于中线右侧时,L取正值 。
(3) 当线元为直线时,其起点、止点的曲率半径为无穷大,以10的45次代替 。
(4) 当线元为圆曲线时,无论其起点、止点与什么线元相接,其曲率半径均等于圆弧的半径 。
(5) 当线元为完整缓和曲线时,起点与直线相接时,曲率半径为无穷大,以10的45次代替;与圆曲线相接时,曲率半径等于圆曲线的半径 。止点与直线相接时,曲率半径为无穷大,以10的45次代替;(6) 当线元为非完整缓和曲线时,起点与直线相接时,曲率半径等于设计规定的值;与圆曲线相接时,曲率半径等于圆曲线的半径 。止点与直线相接时,曲率半径等于设计规定的值;与圆曲线相接时,曲率半径等于圆曲线的半径 。
(7)曲线元要素数据库(DAT-M)可根据线型不同分为各个线元段输入到DAT-M中,即分为直线段、缓和曲线、圆曲线等 。