在我干测量的这么多年中。明白一般的程序大家都有，就是不等长缓和曲线的程序却很少。呵呵。在没有这套程序以前我就有过手算的经历。痛苦呀！！！
仅以此程序送给曾经用手算过不等长缓和曲线坐标的同行们。
一、 文件名（LXQX）主程序：
Lbl 0：｛DZG｝：A“XA”：B“YA”：C“FW”：D“1┛RA”：Z“1┛RB”：F“DKA”：G“DKB”
Lbl 1：｛H｝：H“DKI”：H＞G=＞Goto2 ⊿ P=（Z－D）÷Abs（G－F）：K= Abs（H－F）：S=P×K：Q=C＋（S＋2D）K×90÷π◢ V=C＋（S÷4＋2D）K×45÷（2π）：N=C＋（3S÷4＋2D）K×135÷（2π）：J=C＋（S÷2＋2D）K×45÷π
X=A＋K÷12×（cosC＋4（cosV＋cosN）＋2cosJ＋cosQ）◢ Y=B＋K÷12×（sinC＋4（sinV＋sinN）＋2sinJ＋sinQ）◢
｛T｝：T“LR”：T=2=＞｛LW｝：L“LD=”：W“RA=”：X“XW”=X＋Lcos（Q＋W）◢ Y“YW”=Y＋Lsin（Q＋W）◢⊿｛P｝：P“QH=”：P=3=＞｛L｝：L“LQ=”：X“XQ”=X＋LcosQ◢ Y“YQ”=Y＋LsinQ◢⊿ M=X：U=Y：Prog“XY”：Goto1
Lbl 2：H＞G＋0．0001=＞Goto1⊿ A=X：B=Y：F=G：C=Q：Goto 0
二、 坐标反算（XY）
E“X0”：O“Y0”：M“X1”：U“Y1”：Pol（M－E, U－O）：J＜0 =＞J=J＋360⊿J◢I◢
说明：
1、 该程序对任意曲线都实用；
2、 对该程序要求输入数据的说明：
输入：XA？YA？：要求输入曲线计算起点X，Y坐标；
输入：FW：要求输入曲线计算起点的切线方位角（沿计算方向）；
输入：1┛RA？：要求输入曲线起点半径的倒数（左偏为负，右偏为正）；
输入：1┛RB？：要求输入曲线终点半径的倒数（左偏为负，右偏为正）；
输入：DKA：要求输入曲线计算段起点里程（以米为单位，反算时为负值）；
输入：DKB：要求输入曲线计算段终点里程（以米为单位，反算时为负值）；
输入：DKI：要求输入曲线计算段内待求点里程（以米为单位，反算时为负值）；
输入：LR：是否测设左右点，输入0时只计算线路中心，2时计算左、右；
输入：LD：左、右距离（左时为负，右时为正）；
RA：结构物与线路前进方向夹角（左时为负，右时为正）；
QH：是否计算相对左右点的前后点坐标，输入3时则计算；
LQ：相对左右点的前后距离（前为正，后为负，此式为平行切线方向）；
X0、Y0：为置镜点坐标；
输出说明：
Q：输出求得的曲线（某一里程的）切线方位角；
X、Y：输出曲线上某一里程的中点坐标；
XW、YW：输出曲线上某一里程左、右点坐标；
XQ、YQ：输出曲线上某一里程相对左、右点的前后点坐标；
J：为待测点方位角；
I：为置镜点到待测点距离。
3、 利用此程序计算需分段进行，线形复杂时经几条高速公路匝道实例计算与设计值有误差，最大不到十公分，可采用正算一半、反算一半的方法，使成果与设计符合。(这条还没碰见过)
- 本文固定链接: https://www.chtfs.com/1712/
- 转载请注明: dawang 于 测绘途夫 发表
