直线--缓和曲线-圆曲线坐标计算EXCEL程序1
曲线道路坐标计算(Excel)
曲线道路坐标计算§1 曲线要素计算缓和曲线是在不改变直线段方向和保持圆曲线半径不变的条件下，插入到直线段和圆曲线之间的。
其曲率半径ρ从直线的曲率半径∞（无穷大）逐渐变化到圆曲线的半径R ，在缓和曲线上任意一点的曲率半径ρ与缓和曲线的长度l 成反比，以公式表示为：l1∝ρ 或 C l =⋅ρ（C 为常数，称曲线半径变更率）。
当o l l =时，R =ρ，应有o l R l C ⋅=⋅=ρ以上几式是缓和曲线必要的前提条件。
在实际应用中，可采取符合这一前提条件的曲线作为缓和曲线。
常用的有辐射螺旋线及三次抛物线，我国采用辐射螺旋线。
为了在圆曲线与直线之间加入一段缓和曲线o l ，原来的圆曲线需要在垂直于其切线的方向移动一段距离p ，因而圆心就由'O 移到O ，而原来的半径R 保持不变，如图。
由图中可看出，缓和曲线约有一半的长度是靠近原来的直线部分，而另一半是靠近原来的圆曲线部分，原来圆曲线的两端其圆心角o β相对应的那部分圆弧，现在由缓和曲线所代替，因而圆曲线只剩下缓圆点(HY )到圆缓点(YH )这段长度即y l 。
o β为缓和曲线的切线角，即缓圆点或圆缓点切线与直缓点或缓直点切线的交角，亦即圆曲线HY→YH 两端各延长2ol 部分所对应的圆心角。
γ为缓和曲线总偏角，即从直缓点(ZH )测设缓圆点(HY )或从缓直点(HZ )测设圆缓点(YH )的偏角。
q 为切线增量(切垂距)，即ZH (或HZ )到从圆心O 向ZH (或HZ )的切线作垂线垂足的距离。
p 为圆曲线内移值，即垂线(从圆心O 向ZH (或HZ )的切线作垂线)长与圆曲线半径R 之差。
§1.1 不等长缓和曲线要素计算：在铁路曲线测设中，线路曲线一般是由相等的两条缓和曲线中间加一个圆曲线构成，有时还会出现由两个不等长的缓和曲线中间加一个圆曲线构成的特殊情况，如图：缓和曲线长分别为1o l 、2o l ， 切线长分别为1T 、2T ，曲线偏角(线路转角)为α，圆曲线半径为R ，圆曲线长为y l ，曲线长为L ，外矢距为E ，切曲差为J ，(缓和曲线后)圆曲线内移值分别为1p 、2p ，(缓和曲线)切线增量分别为1q 、2q ，缓和曲线偏角分别为1o β、2o β ， 回旋线参数分别为121o Rl A =、222o Rl A =各曲线要素计算公式如下：231112402R ll q o o -=232222402Rll q o o -=341211268824R lR l p o o -=342222268824RlR l p o o -=ααsin )(2)(12111p p tgp R q T -+++=ααsin )(2)(21222p p tgp R q T -+++=Rl R l o o o ππβ111901802=⋅=Rl R l o o o ππβ222901802=⋅=︒--++=180)(2121Rl l L o o o o πββα从以上公式可以看出，当21o o l l =时，就是等长(对称)缓和曲线的情况。
公路逐桩坐标计算程序
2 )sinαAB +(
P 点为顺时针方向时，其方位角为 αAB+900 P 点为逆时针方向时，其方位角为 αAB−900
第一段缓和曲线时：以直缓点（ZH）为起点计算，αAB 为 ZH 点的坐标方位角，L 为 P 点距 ZH 点的距离。 第二段缓和曲线时：以缓直点（HZ）为起点计算，αAB 为 HZ 点坐标方位角的反 方向即 HZ 点方位角加 180 度，L 为 P 点距 HZ 点的距离，加减 90 度刚好与第一 段缓和曲线相反。图如下：
缓和曲线转角公式：β
=
L2
2RLs
（2） 边桩坐标计算公式
左侧
XA=XP+T1cos(αAB± β-900) YA=YP+T1sin(αAB± β-900)
右侧
XB=XP+T2cos(αAB± β+900) YB=YP+T2sin(αAB± β+900)
*第二段缓和曲线计算边桩坐标时， 注意加减 90 度时， 与第一段缓和曲线相反。
2、 缓和曲线段 （1） 中桩坐标计算公式
Xp=X1+(L − Yp=Y1+(L −
L5 40R 2 Ls L5 40R 2 Ls
2 )cosαAB+(
L3 6RLs L3 6RLs
− −
L7 336R 3 ������s 3 L7 336R 3 ������s 3
)sin(αAB±900) )sin(αAB±900)
右侧
XB=XP+T2cos(αAB± β+900) YB=YP+T2sin(αAB± β+900)
第二章
公路导线测量计算
1
公路导线测量为附合导线测量，按路线前进方向测量右角。β
EXCEL曲线坐标计算公式
公式解析一．坐标转换X =A +N COSα-E SINαY =B +N SINα+E COSαN=(X-A) COSα±(Y-B)SINαE=(Y-B)COSα±(X-A)SINαA,B为施工坐标系坐标原点α为施工坐标系与北京坐标系X轴的夹角（旋转角）即大地坐标系方位角X,Y为北京坐标值 N，E为施工坐标值二．方位角计算1.直线段方位角: α=tanˉ¹ [(Yb-Ya)/(Xb-Xa)]2.交点转角角度: α=2 tanˉ¹ (T/R)计算结果①为﹢且＜360，则用原数;②为﹢且＞360，则减去360;③为﹣，则加上180.3.缓和曲线上切线角: α=ƟZH±90°*Lo²/(π*R* Ls）α= Lo/(2ρ)=Lo²/(2 A²)=Lo²/(2R*Ls)ρ—该点的曲率半径4.圆曲线上切线角: α=ƟHY±180°*Lo/(π*R）ƟZH—直缓点方位角, ƟHY—缓圆点方位角,注：以计算方向为准，左偏，取"﹣"；右偏，取"﹢"。
左偏，则第一段缓和曲线和圆曲线上取"﹣"，第二段缓和曲线上取"﹢" ；右偏，则第一段缓和曲线和圆曲线上取"﹢"，第二段缓和曲线上取"﹣" .。
符号说明:A—回旋线参数（A²=R* Ls） Ls—缓和曲线长度R—曲线半径Lo—曲线长度：计算点位到特殊点(ZH、HY、YH、HZ)的长度三．坐标值计算1．直线段坐标计算公式：直线两端点A.B间距离为S；A点坐标为A(Xa, Ya)；方位角为αXb= Xa+S*cosαYb= Ya+S*sinα2．缓和曲线及圆曲线坐标计算公式：①缓和曲线坐标计算公式：X=XZH+(Lo-Lo^5/(40*R^2*Ls^2)+Lo^9/(3456*R^4*Ls ^4)-Lo^13/(599040*R^6*Ls^6)+Lo^17/(175472640*R ^8*Ls^8))*cosα-(Lo^3/(6*R*Ls)-Lo^7/(336*R^3*L s^3)+Lo^11/(42240*R^5*Ls^5)-Lo^15/(9676800*R^7 *Ls^7)+Lo^19/(3530096640*R^9*Ls^9))*sinαY=YZH+(Lo-^5/(40*R^2*Ls^2)+Lo^9/(3456*R^4*Ls^4 )-Lo^13/(599040*R^6*Ls^6)+Lo^17/(175472640*R^8 *Ls^8))*sinα+(Lo^3/(6*R*Ls)-Lo^7/(336*R^3*Ls^ 3)+Lo^11/(42240*R^5*Ls^5)-Lo^15/(9676800*R^7*L s^7)+Lo^19/(3530096640*R^9*Ls^9))* cosα符号说明:XZH—直缓点X坐标值 YZH—直缓点Y坐标值 A—回旋线参数（A²=R* Ls）Lo—计算点位到特殊点的长度 Ls—缓和曲线长度R—曲线半径α—方位角注:式中，紫色部分为缓和曲线任意点的坐标增量(支距坐标)。
圆曲线、缓和曲线计算方法
● 圆曲线方法一：sin (1cos )180i i i i i i x R y R l R ϕϕϕπ⎧⎪=⎪=-⎨⎪︒⎪=⋅⎩——i l 为待定点i P 至起点间的弧长i ϕ为i l 所对的圆心角R 为曲线半径方法二：11802l A R π︒=⋅⋅ 2sin l R A =⋅00cos(/)sin(/)x x l A y y l A =+⋅+-⎧⎨=+⋅+-⎩起点方位角左减右加起点方位角左减右加——00(,)x y 为圆曲线起点坐标方法三：180l A R π︒=⋅ 00cos(/)sin(/)x x R B A y y R B A =+⋅+-⎧⎨=+⋅+-⎩——l 为圆曲线上任意一点距起点距离00(,)x y 为圆曲线圆心坐标B 为圆心到圆曲线起点的方位角，A 为任意点对应的圆心角● 缓和曲线522030406l x l R l ly Rl ⎧=-⎪⎪⎨⎪=⎪⎩——l 为曲线上任一点至起点的曲线长R 曲线半径0l 为缓和曲线全长圆曲线、缓和曲线计算方法1、直线段：先由JD1以及JD2的坐标算出JD1到JD2的方位角，即直线段方位角A ，故可算出HZ 、ZH 坐标及其直线段各点坐标。
2、缓和曲线：以HZ 、ZH 为起点，缓和曲线上任意一点离HZ 、ZH 距离为l ，利用公式522003040()6l x l R l l R ly Rl ⎧=-⎪⎪⎨⎪=⎪⎩为缓和曲线全长,为圆曲线半径算出该点的相对起点坐标，利用arctan y x算出该点相对起点的方位角B ，再根据线路走向及直线段方位角可算出该点的方位角C （顺时针加，逆时针减），用可求出该点相对起点的距离D ，最后用00cos sin x x D C y y D C =+⎧⎨=+⎩可求出该点的坐标。
（00(,)x y 为缓和曲线起点的坐标）3、圆曲线：用上述方法求出圆曲线两端点HY 、YH 坐标，算出HY 到YH 的方位角F ，以及两点间的距离E ，用12arccos ER可算出两端点连线与起点到圆心连线的夹角G ，根据线路走向求出起点到圆心的方位角H （H=F+/-G ），00(,)x y 圆曲线为起点坐标，根据00cos sin x x R H y y R H=+⎧⎨=+⎩，求出圆心坐标。
坐标计算程序
测量实用技术程序一：线路坐标计算程序采用的公式1：曲线要素公式内移值P=Ls xy2/(24R)切线增量q=Ls/2-Ls xy3/(240R xy2)切线长T=(R+p)tg(A/2)+q切线角B0=Ls/(2R)曲线外距E=(R+P)SEC(A/2)-R2:缓和曲线任意一点与ZH点之间的弦长公式C=L-L xy5/(90R xy2Ls xy2)3:缓和曲线任意一点的偏角公式S=(L/LS)xy2S0S0=1/3B04：其他采用坐标增量公式△X=LCOSW△Y=Lsinw然后与JD的坐标相加X=XJD+△XY=YJD+△Y其中L为计算点至交点JD的长度W为JD到计算点的方位角XJD为交点的X坐标YJD为交点的Y坐标圆曲线增设缓和曲线示意图程序功能可计算直线，圆曲线，缓和曲线上任意一点的中桩坐标及其法线边桩坐标A:曲线偏角 T切线长 E外失距 R 园曲半茎ZH 直缓点 HY缓园点 HY 缓园点 QZ曲中 YH 缓园 HZ 缓直线路直曲线坐标计算程序1：A:C”LS=”:D”JD=”:R”R=”2: P=C xy2/24/R-C xy4/2688/R xy3 3:Q=C/2-C xy3/240/R xy24: B=90C/R/∏5：T=(R+P)Tg(absA/2)+Q6: E=(R+P)/COS(A/2)-R7；L=(absA-2B) ∏R/180+2C 8: G”ZH”=D-T▲ H” HY”=G+C 9: I”QZ”=G+L/2 10:K”YH”=G+L-C 11:M”HZ”=G+L 12:N”JDX=”:W”JDY=”:F”FANGWEIJIAO”:J 13:A＜0→S=-1: S=1⊿14:U=F+A/2+90S 15:V=E+R 16:B=N+VCOSU 17:O=W+VsinU 18:Lb1 19:{Z} 20:Z≤ G →L=D-Z 21:V=F+18022:U=F23:Goto 2 ⊿24:Z≤H→L=Z-G25:V=L-L xy5/90/(RC) xy226:L=30SL xy2/R/∏/C27:P=F+18028:Q=F+L29:U=F+3L30:Goto5⊿31:Z≤K→L=F+A/2+180+90S+180S(Z-I)/R/∏32:U=L+90S33:Goto5⊿34:Z≤M→L=M-Z35:V=L-L xy5/90(RC) xy236:L=30Sl xy2/(R∏C)37:P=F+A38:Q=P+180-L39:U=Q-2L-18040:Goto4⊿41:Z＞M→L=Z-M+T42:V=F+A43:U=V44:Goto2⊿45:Lb1 246:X =N+LcosV▲47:Y =W+LsinV▲48: Goto6⊿49: Lb1 350:{E}51:P”XL”=X+E cos(U-90) ▲52:Q”YL”=Y+E sin(U-90) ▲53:P”XR”=X+E cos(U+90) ▲54:Q”YR”=Y+E sin(U+90) ▲55: Goto 156: Lb1 457:X=N+T cosP+V cosQ▲58:Y=W+T sinP+V sinQ▲59: Goto 660: Lb1 561: X=B+R cosL▲62: Y=O+R sinL▲63; Goto 664: Lb1 665:J=1→ Goto 3⊿66:Goto1 说明：A曲线偏角（左偏为负数）LS缓和曲线长D JD交点里程R园曲线半径N JD的X坐标W JD的Y坐标F是起算方位角ZH→JD Z待算里程 J 输入1时计算中桩及两侧坐标输入0仅计算中桩坐标E输入中桩到两侧的距离XL YL 是中桩左侧的坐标XR LR是右侧坐标已知两点坐标计算方位角程序1：X”XA=”:Y”YH=”:S”XB=”:T”YB=”2：A=(S-X)3：B=(T-Y)4：POL(A,B) ▲5:W≥0→W=W: W=W+360⊿6：F”FWJ”=W▲。
纯Excel公式开发的施工测量数据计算表V1.4
纯Excel公式开发的施工测量数据计算表版本号：1.4功能：可以根据输入的起点桩号、各交点坐标、各交点的圆曲线半径、缓和曲线长（支持不对称缓和曲线）直接生成标准的《直线曲线及转角表》和积木法平曲线参数表。
可以计算非对称基本型（直-缓1-圆-缓2）以及下列由其变化而成共6种线型的任意交角中边桩坐标。
1、单一完整缓和曲线；分缓1和缓2两种情况，此曲线为前直线、圆曲线和另一缓和曲线长度均为0的特殊曲线。
2、单一圆曲线，即前直线和缓和曲线长度均为零的曲线。
3、卵型曲线，即圆曲线长为0的曲线（前直线有无不限，此类曲线在程序中要加入防止R被作除数时出错）4、缓1+圆曲线型曲线，即缓2长度为零的曲线（前直线长度不限）5、圆曲线+缓2型曲线，即前直线和缓1长度为零的曲线。
6、单一直线，即所有缓和曲线和圆曲线的长度均为零的曲线（常用于路线终点直线计算，此类曲线需防止已知量为零并被当作被除数时出错）可以根据输入的桩号自动生成用于AutoCAD脚本绘图的代码，拷贝后保存为*.scr文件，在任意版本AutoCAD中加载展绘中线和桩号。
可以用近似公式计算竖曲线高程。
对于以变量赋值形式作为数据库的CASIOfx-5800P计算器程序，可自动生成交点法、积木法、竖曲线数据库。
使用说明：此表用纯Excel公式编成，理论上可计算多达30000个交点，但考虑到实际情况，只拉出了200个交点的表格，所有的已知数据在此表中输入，结果在《直曲表》中查看。
此表格经多条100公里以上的线路计算，误差均小于1mm，可用于设计数据复核。
因为很多线路的桩号，缓和曲线参数，主点桩号等数值均只取2位小数，而交点坐标一般会给出小数点后三位，半径和缓和曲线长一般均为整数，用此表计算可以帮助有程序的朋友在建立数据库时提高计算精度。
又有朋友手上只有积木法的程序，对于交点法不想另编，本表又提供了一种从交点法表格直接转为积木法的功能，已知数据还是没有变，根据交点个数，拉了800行。
EXCEL曲线坐标计算公式
公式解析一．坐标转换X =A +N COSα-E SINαY =B +N SINα+E COSαN=(X-A) COSα±(Y-B)SINαE=(Y-B)COSα±(X-A)SINαA,B为施工坐标系坐标原点α为施工坐标系与北京坐标系X轴的夹角（旋转角）即大地坐标系方位角X,Y为北京坐标值N，E为施工坐标值二．方位角计算1.直线段方位角: α=tanˉ¹ [(Y b-Y a)/(X b-X a)]2.交点转角角度: α=2 tanˉ¹ (T/R)计算结果①为﹢且＜360，则用原数;②为﹢且＞360，则减去360;③为﹣，则加上180.3.缓和曲线上切线角: α=ƟZH±90°*Lo²/(π*R* Ls）α= Lo/(2ρ)=Lo²/(2 A²)=Lo²/(2R*Ls)ρ—该点的曲率半径4.圆曲线上切线角: α=ƟHY±180°*Lo/(π*R）ƟZH—直缓点方位角, ƟHY—缓圆点方位角,注：以计算方向为准，左偏，取"﹣"；右偏，取"﹢"。
左偏，则第一段缓和曲线和圆曲线上取"﹣"，第二段缓和曲线上取"﹢" ；右偏，则第一段缓和曲线和圆曲线上取"﹢"，第二段缓和曲线上取"﹣" .。
符号说明:A—回旋线参数（A²=R* Ls）Ls—缓和曲线长度R—曲线半径Lo—曲线长度：计算点位到特殊点(ZH、HY、YH、HZ)的长度三．坐标值计算1．直线段坐标计算公式：直线两端点A.B间距离为S；A点坐标为A(X a, Y a)；方位角为αX b= X a+S*cosαY b= Y a+S*sinα2．缓和曲线及圆曲线坐标计算公式：①缓和曲线坐标计算公式：X=X ZH+(Lo-Lo^5/(40*R^2*Ls^2)+Lo^9/(3456*R^4*Ls^4)-Lo^13/( 599040*R^6*Ls^6)+Lo^17/(175472640*R^8*Ls^8))*cosα-(Lo^3/(6 *R*Ls)-Lo^7/(336*R^3*Ls^3)+Lo^11/(42240*R^5*Ls^5)-Lo^15/(9 676800*R^7*Ls^7)+Lo^19/(3530096640*R^9*Ls^9))*sinαY=Y ZH+(Lo-^5/(40*R^2*Ls^2)+Lo^9/(3456*R^4*Ls^4)-Lo^13/(59 9040*R^6*Ls^6)+Lo^17/(175472640*R^8*Ls^8))*sinα+(Lo^3/(6* R*Ls)-Lo^7/(336*R^3*Ls^3)+Lo^11/(42240*R^5*Ls^5)-Lo^15/(96 76800*R^7*Ls^7)+Lo^19/(3530096640*R^9*Ls^9))* cosα符号说明:X ZH—直缓点X坐标值Y ZH—直缓点Y坐标值A—回旋线参数（A²=R* Ls）Lo—计算点位到特殊点的长度Ls—缓和曲线长度R—曲线半径α—方位角注:式中，紫色部分为缓和曲线任意点的坐标增量(支距坐标)。
完整及不完整缓和曲线 (1)
关于不同类型缓和曲线的起点、终点曲率半径判断方法目前在匝道或线路施工坐标计算中经常遇到缓和曲线，实际中相信有很多测友选择用积木法或叫线元法正反算程序进行线路坐标计算，这就牵涉到线元的起点终点曲率半径判断的问题，一般的直线元，圆曲线元的起点终点半径判断，比较容易，可能令大家感觉麻烦的就是缓和曲线起点终点半径判断问题，缓和曲线有时候判断算对了，有时候却坐标算不对，究其原因，其实问题出于该缓和曲线是否是完整缓和曲线引起的。
关于这点，相关的课本教材上没有明确的讲述，网上对此问题的解释也是散见于不同的论文著作中，对于测量新手来说，线元法程序是非常适用上手的，但却往往因为遇到不完整缓和曲线的起点或终点的半径判断计算不出来导致坐标计算错误，的确是件令人恼火的事情，在此我就把自己的判断经验做一论述，给用线元法程序的测友们一同分享，当然高手们请一笑而过，也可留下你的经验与大家一起分享交流学习。
第一：先说说完整缓和曲线和不完整缓和曲线以及不对称缓和曲线与对称缓和曲线的概念问题，以免混为一谈 .1.当对于单独一段缓和曲线从其完整与否来讲是分为完整与不完整两类；当对于一个单交点内的两段缓和曲线(即常说的第一缓和曲线和第二缓和曲线而言)又有对称缓和曲线与不对称缓和曲线之分。
由此看来，完整与对称与否是针对缓和曲线两个方面来看待区分的。
2.缓和曲线我们的测量教材上讲述的其实就是完整缓和曲线，也可以知道缓和曲线上：各个点的半径是不同的，起点到终点的半径值过度是从正无穷大到所接圆曲线半径之过度如从 ZH 向 HY 方向；或者是从所接圆曲线半径值向正无穷大过度的，如从 YH 向 HZ 方向。
那么由此可以不难判断出来，完整缓和曲线就是符合上述特征的，那么不完整的缓和曲线就是不符合上述特征的，但是线路上的平曲线设计时候一般缓和曲线不单独存在的，整体上缓和曲线前或后一般都是要连接一个圆曲线的，那么不完整缓和曲线其实就是在完整缓和曲线上截取的一段，一般就是去掉了半径无穷大的那端而是从某个点开始的半径值向所接圆曲线半径值过度的。
3.1不对称缓和曲线加圆曲线逐桩坐标计算讲解
β2 = l2 180 2R π
主点里程
kZH = K JD T1
K HY K ZH l1
K YH K ZH L l 2
K HZ K ZH L
二、不等长的缓和曲线逐桩坐标计算 1、已知条件： (1)设计文件所给已知条件：
交点坐标及里程、曲线半径R、第一缓和曲线长、第二缓和曲线 长、转向角、ZH点切线方位角。
CF DF ( R P2 ) tan
T 1 AJ AG
2
m1
2
m2
P1 P2 sin
P1 P2 sin
T 2 DJ DF
则切线长
l1 l m1 = 2 240R 2
3 1
l1 p1 = 24R
2 2
2
l2 l m2 = 2 240R 2
子情境3：特殊曲线测量
3.1 不等长的缓和曲线+圆曲线逐桩坐标计算
3.1不等长缓和曲线逐桩坐标计算
JD
T1
缓和曲线
T2
圆曲线
直线
ZH
1 l0
HY
YH
缓和曲线
l02
HZ
直线
不等长的缓和曲线+圆曲线通常敷设在地形复杂地
区以及铁路旧线改造中（如图）。 其组成为：直线-缓和曲线（A1)-圆曲线-缓和曲 线（A2)-直线
X P X HZ x P cos ky P sin P sin ky P cos YP YHZ x
k 1曲线右偏时 K 1曲线左偏时
360 HZ切
(2)地面上的已知条件：
导线点坐标及高程、导线控制点的实际位置。
缓和曲线、 圆曲线计算方法
缓和曲线计算方法(ZH～HY)中线首先计算直线段坐标方位角（即ZH～JD坐标方位角），及ZH点坐标。
备用偏角公式：{30*L/（π*RLS）缓和曲线}1、计算待求点偏角=（（L/10）2*（57296/（RLS））/60。
其中L=待求点至ZH距离、R=圆曲线半径、LS=缓和曲线长。
2、待求点方位角=直线方位角±待求点偏角。
（曲线左转-偏角，曲线右转+偏角）3、待求点至ZH点弦长=L—L5/（90*R2*LS 2），其中L=待求点至ZH距离（里程）、R=圆曲线半径。
4、待求点坐标：X=ZH点X坐标+COS（待求点方位角）*弦长Y= ZH点Y坐标+SIN（待求点方位角）*弦长缓和曲线计算左右边线坐标(ZH～HY)1、左侧方位角=（待求点方位角±2倍偏角=直线方位角±3倍偏角）—边线与中线夹角。
2、右侧方位角=（待求点方位角±2倍偏角=直线方位角±3倍偏角）+边线与中线夹角。
3、左侧边线坐标：X=该点中线X坐标+COS（左侧方位角）*边线至中线距离Y=该点中线Y坐标+SIN（左侧方位角）*边线至中线距离4、右侧边线坐标：X=该点中线X坐标+COS（右侧方位角）*边线至中线距离Y=该点中线Y坐标+SIN（右侧方位角）*边线至中线距离圆曲线计算方法（HY～YH）中线注：（ZY-YZ）同理，方位角=用直线方位角-待求点偏角首先计算直线段坐标方位角（即ZH～JD坐标方位角），及HY点坐标。
1、求出缓圆点（HY）偏角=（LS*90）/（π* R）。
2、求待求点偏角=（L*90）/（π* R）。
其中：L=待求点至HY距离（里程）、R=圆曲线半径、LS=缓和曲线长。
3、待求点至HY点弦长=2* R*SIN（待求点偏角）。
4、待求点方位角=直线方位角±HY点偏角±待求点偏角，（曲线左转-偏角，曲线右转+偏角）。
5、待求点坐标：X=HY点X坐标+COS（待求点方位角）*弦长Y=HY点Y坐标+SIN（待求点方位角）*弦长圆曲线计算左右边线坐标1、左侧方位角=（待求点方位角±偏角—边线与中线夹角）。
EXCEL快速计算任何难度的道路中线、边线坐标
EXCEL快速计算任何难度的道路中线、边线坐标主程序Public i As IntegerPublic pi As DoubleSub TP()Dim ii As IntegerDim k(1000) As DoubleDim xzq, yzq, kq, xzh, yzh, kzh, xjd, yjd, kjd, khy, kyh As Double'直线区域pi = 3.14159265358979xzq = 71862.642yzq = 63474.651kq = 0 '因为直线连接终点为ZH点，与圆曲线起点为同一点，所以在直线区域不定义ZH点参数'直线区域'曲线区域xzh = 71858.3267yzh = 63375.2684kzh = 99.4763xhz = 71909.3687yhz = 63283.8076 '曲线区域定义内容有：ZH（坐标、里程）、HZ（坐标、里程）、JD（坐标、里程）khz = 212.3392 'R（半径）、LS（缓和曲线长度）、HY（里程）、YH（里程）xjd = 71855.658yjd = 63313.806kjd = 160.9966ls = 30r = 75khy = 129.4763kyh = 182.3385'曲线区域i = 2 '从第二格开始读取数据所以定义I=2ii = 1 '桩号从第一个开始启用，所以定义II=2Dok(ii) = Workbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 1) '定义桩号等于读取数据If Workbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 1) = "" Then Exit Do '当没有数据读取时退出循环If k(ii) < kq Then '若计算点超过计算起点给予提示并退出程序MsgBox ("猪啊！！你的输入的桩号居然超过计算起点桩号")Exit SubElseIf k(ii) <= kzh Then '若计算点在ZH点前，则进入直线程序Call zx(xzq, yzq, kq, xzh, yzh, kzh, k(ii))ElseIf kzh < k(ii) And k(ii) <= khy Then '若计算点在ZH和HY之间则调入前段缓和曲线程序Call qhhqx(xzh, yzh, kzh, xhz, yhz, khz, xjd, yjd, kjd, ls, r, k(ii))ElseIf khy < k(ii) And k(ii) <= kyh Then '若计算点在HY和YH之间则调入圆曲线程序Call yqx(xzh, yzh, kzh, xhz, yhz, khz, xjd, yjd, kjd, ls, r, k(ii))ElseIf kyh < k(ii) And k(ii) <= khz Then '若计算点在YH和HZ之间则调入后段缓和曲线程序Call hhhqx(xzh, yzh, kzh, xhz, yhz, khz, xjd, yjd, kjd, ls, r, k(ii))ElseMsgBox ("笨啊！！数据已超出计算范围了") '若出现超出范围的桩号则给与提示并退出程序Exit SubEnd Ifi = i + 1ii = ii + 1LoopEnd Sub直线模块Sub zx(ByV al xzq As Double, ByV al yzq As Double, ByV al kq As Double, ByV al xzh As Double, ByV al yzh As Double, ByV al kzh As Double, ParamArray k())fw = fwj(xzh, xzq, yzh, yzq) '首先调入方位角程序计算直线方位角x = xzq + (k(ii) - kq) * Cos(fw) '然后根据桩号和长度计算出坐标值y = yzq + (k(ii) - kq) * Sin(fw)zdfm = dfm(fw) '将弧度形式的前进方位角转换度分秒形式'输出坐标值以弧度和度分秒形式的前进方位角Workbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 2) = xWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 3) = yWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 4) = fwWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 5) = zdfmEnd Sub圆曲线模块Sub yqx(ByV al xzh As Double, ByV al yzh As Double, ByV al kzh As Double, ByV al xhz As Double, ByV al yhz As Double, ByV al khz As Double, ByV al xjd As Double, ByV al yjd As Double, ByV al kjd As Double, ByV al ls As Double, ByV al r As Double, ParamArray k())l = Abs(k(ii) - kzh) '计算ZH点（因为以直缓点起算）到待求桩号的弧度长度ly = l - ls / 2 '计算圆弧长度p = ls ^ 2 / 24 / r - ls ^ 4 / 2688 / r ^ 3 '曲线内移值m = ls / 2 - ls ^ 3 / 240 / r ^ 2 '曲线切线长增量u = r * Sin(ly / r) + m '偏量坐标计算v = r * (1 - Cos(ly / r)) + p'调入方位角fwq = fwj(xjd, xzh, yjd, yzh) '计算ZH点方位角fwh = fwj(xhz, xjd, yhz, yjd) '计算HZ点方位角（此角作用是用来推算曲线是左偏还是右偏）'调入偏角判定nq = n(fwq, fwh) '计算偏角方向，左偏为-1右偏为1'计算坐标x = u * Cos(fwq) - nq * v * Sin(fwq) + xzhy = u * Sin(fwq) + nq * v * Cos(fwq) + yzhd = (90 * (2 * l - ls) / pi / r) * pi / 180 '计算圆曲线上的偏角（此句要点为角度必须转换为弧度即：pi/180）fw = fwq + d * nq '计算前进方位角zdfm = dfm(fw) '将弧度形式的前进方位角转换度分秒形式'输出坐标值以弧度和度分秒形式的前进方位角Workbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 2) = xWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 3) = yWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 4) = fwWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 5) = zdfmEnd Sub前缓和段Sub qhhqx(ByV al xzh As Double, ByV al yzh As Double, ByV al kzh As Double, ByV al xhz As Double, ByV al yhz As Double, ByV al khz As Double, ByV al xjd As Double, ByV al yjd As Double, ByV al kjd As Double, ByV al ls As Double, ByV al r As Double, ParamArray k())l = Abs(k(ii) - kzh) '计算测点到ZH点的距离u = l - l ^ 5 / 40 / r ^ 2 / ls ^ 2 + l ^ 9 / r ^ 4 / ls ^ 4 / 3456 '计算偏量v = l ^ 3 / 6 / r / ls - l ^ 7 / 336 / r ^ 3 / ls ^ 3'用公式二则以下两计算可省略(圆曲线和后缓和曲线也同理)Rem t = Atn(v / u)Rem s = Sqr(u ^ 2 + v ^ 2)'调入方位角计算fwq = fwj(xjd, xzh, yjd, yzh) '计算ZH点方位角fwh = fwj(xhz, xjd, yhz, yjd) '计算HZ点方位角（此角作用是用来推算曲线是左偏还是右偏）'调入偏角判定nq = n(fwq, fwh) '计算偏角方向，左偏为-1右偏为1'结果计算Rem x = xzh + s * Cos(fwq + nq * t)Rem y = yzh + s * Sin(fwq + nq * t)x = u * Cos(fwq) - nq * v * Sin(fwq) + xzh '经过测试，计算结果中的两种公式计算结果是一样的y = u * Sin(fwq) + nq * v * Cos(fwq) + yzhd = (90 * l * l / pi / r / ls) * pi / 180 '计算缓和曲线上的偏角（此句要点为角度必须转换为弧度即：pi/180）fw = fwq + d * nq '计算前进方位角zdfm = dfm(fw) '将弧度形式的前进方位角转换度分秒形式'输出坐标值以弧度和度分秒形式的前进方位角Workbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 2) = xWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 3) = yWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 4) = fwWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 5) = zdfmEnd Sub后缓和段模块Sub hhhqx(ByV al xzh As Double, ByV al yzh As Double, ByV al kzh As Double, ByV al xhz As Double, ByV al yhz As Double, ByV al khz As Double, ByV al xjd As Double, ByV al yjd As Double, ByV al kjd As Double, ByV al ls As Double, ByV al r As Double, ParamArray k())l = Abs(k(ii) - khz) '计算测点到HZ点的距离（后缓和曲线是以HZ点为起点）u = l - l ^ 5 / 40 / r ^ 2 / ls ^ 2 + l ^ 9 / r ^ 4 / ls ^ 4 / 3456 '计算偏量v = l ^ 3 / 6 / r / ls - l ^ 7 / 336 / r ^ 3 / ls ^ 3Rem t = Atn(v / u)Rem s = Sqr(u ^ 2 + v ^ 2)'调入方位角计算fwq = fwj(xjd, xzh, yjd, yzh) '计算ZH点方位角fwh = fwj(xhz, xjd, yhz, yjd) '计算HZ点方位角（此角作用是用来推算曲线是左偏还是右偏）'调入偏角判定nh = n(fwh, fwq) '计算偏角方向，左偏为-1右偏为1（注意：因为是从后HZ点起算，所以必须将HZ点方位角放在前ZH放在后）'结果计算Rem x = xzh + s * Cos(fwq + nq * t)Rem y = yzh + s * Sin(fwq + nq * t)x = xhz - (u * Cos(fwh) - nh * v * Sin(fwh)) '经过测试，计算结果中的两种公式计算结果是一样的y = yhz - (u * Sin(fwh) + nh * v * Cos(fwh))d = (90 * l * l / pi / r / ls) * pi / 180 '计算缓和曲线上的偏角（此句要点为角度必须转换为弧度即：pi/180）fw = fwh + d * nh '计算前进方位角zdfm = dfm(fw) '将弧度形式的前进方位角转换度分秒形式'输出坐标值以弧度和度分秒形式的前进方位角Workbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 2) = xWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 3) = yWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 4) = fwWorkbooks("单交点平曲线.xls").Worksheets("sheet1").Cells(i, 5) = zdfmEnd Sub偏角模块Function n(ByV al fw1 As Double, ByV al fw2 As Double) As Doublepj = fw1 + pi - fw2 '前进的右角pjIf pj - pi > 0 Then '当右角pj-pi 〉0时为左偏否则为右偏n = -1Else: n = 1End IfEnd Function方位角模块Function fwj(ByV al x1 As Double, ByV al x2 As Double, ByV al y1 As Double, ByV al y2 As Double) As Double'计算增量x0 = x1 - x2y0 = y1 - y2'由增量判断方位角所在象限，不同象限取不同的值If x0 = 0 And y0 > 0 Thenfwj = pi / 2 '当在大地坐标中偏量在X轴上的值时ElseIf x0 = 0 And y0 < 0 Thenfwj = 3 * pi / 2 '当在大地坐标中偏量在负X轴上的值时ElseIf x0 < 0 Thenfwj = Atn(y0 / x0) + pi '当在大地坐标中偏量在第二第三象限上的值时ElseIf x0 > 0 And y0 < 0 Thenfwj = Atn(y0 / x0) + 2 * pi '当在大地坐标中偏量在第四象限上的值时Elsefwj = Atn(y0 / x0) '当在大地坐标中偏量在第一象限上的值时End IfEnd Function度分秒模块Function dfm(ByV al ao As Double) As V ariantao = ao * 180 / pi '将弧度转化为度jd = Int(ao) '计算度jf = Int(ao * 60 - jd * 60) '计算分jmx = (ao - jd - jf / 60) * 3600 '计算秒jm = Left(jmx, 8) '因为拆分出来的秒数经常占到十多位，所以只取秒数的前八位dfm = jd & "°" & jf & "′" & jm & "″" '连接度分秒End Function边桩公式(此公式在电子表格中直接输入)：=B2+J2*COS(D2+RADIANS(L2)+PI()) =C2+J2*SIN(D2+RADIANS(L2)+PI()) =B2+K2*COS(D2+RADIANS(M2)) =C2+K2*SIN(D2+RADIANS(M2))=B3+J3*COS(D3+RADIANS(L3)+PI()) =C3+J3*SIN(D3+RADIANS(L3)+PI()) =B3+K3*COS(D3+RADIANS(M3)) =C3+K3*SIN(D3+RADIANS(M3))=B4+J4*COS(D4+RADIANS(L4)+PI()) =C4+J4*SIN(D4+RADIANS(L4)+PI()) =B4+K4*COS(D4+RADIANS(M4)) =C4+K4*SIN(D4+RADIANS(M4))=B5+J5*COS(D5+RADIANS(L5)+PI()) =C5+J5*SIN(D5+RADIANS(L5)+PI()) =B5+K5*COS(D5+RADIANS(M5)) =C5+K5*SIN(D5+RADIANS(M5))=B6+J6*COS(D6+RADIANS(L6)+PI()) =C6+J6*SIN(D6+RADIANS(L6)+PI()) =B6+K6*COS(D6+RADIANS(M6)) =C6+K6*SIN(D6+RADIANS(M6))=B7+J7*COS(D7+RADIANS(L7)+PI()) =C7+J7*SIN(D7+RADIANS(L7)+PI()) =B7+K7*COS(D7+RADIANS(M7)) =C7+K7*SIN(D7+RADIANS(M7))=B8+J8*COS(D8+RADIANS(L8)+PI()) =C8+J8*SIN(D8+RADIANS(L8)+PI()) =B8+K8*COS(D8+RADIANS(M8)) =C8+K8*SIN(D8+RADIANS(M8))=B9+J9*COS(D9+RADIANS(L9)+PI()) =C9+J9*SIN(D9+RADIANS(L9)+PI()) =B9+K9*COS(D9+RADIANS(M9)) =C9+K9*SIN(D9+RADIANS(M9))=B10+J10*COS(D10+RADIANS(L10)+PI()) =C10+J10*SIN(D10+RADIANS(L10)+PI()) =B10+K10*COS(D10+RADIANS(M10)) =C10+K10*SIN(D10+RADIANS(M10))=B11+J11*COS(D11+RADIANS(L11)+PI()) =C11+J11*SIN(D11+RADIANS(L11)+PI()) =B11+K11*COS(D11+RADIANS(M11)) =C11+K11*SIN(D11+RADIANS(M11))=B12+J12*COS(D12+RADIANS(L12)+PI()) =C12+J12*SIN(D12+RADIANS(L12)+PI()) =B12+K12*COS(D12+RADIANS(M12)) =C12+K12*SIN(D12+RADIANS(M12))=B13+J13*COS(D13+RADIANS(L13)+PI()) =C13+J13*SIN(D13+RADIANS(L13)+PI()) =B13+K13*COS(D13+RADIANS(M13)) =C13+K13*SIN(D13+RADIANS(M13))=B14+J14*COS(D14+RADIANS(L14)+PI()) =C14+J14*SIN(D14+RADIANS(L14)+PI()) =B14+K14*COS(D14+RADIANS(M14)) =C14+K14*SIN(D14+RADIANS(M14))=B15+J15*COS(D15+RADIANS(L15)+PI()) =C15+J15*SIN(D15+RADIANS(L15)+PI()) =B15+K15*COS(D15+RADIANS(M15)) =C15+K15*SIN(D15+RADIANS(M15))=B16+J16*COS(D16+RADIANS(L16)+PI()) =C16+J16*SIN(D16+RADIANS(L16)+PI()) =B16+K16*COS(D16+RADIANS(M16)) =C16+K16*SIN(D16+RADIANS(M16))=B17+J17*COS(D17+RADIANS(L17)+PI()) =C17+J17*SIN(D17+RADIANS(L17)+PI()) =B17+K17*COS(D17+RADIANS(M17)) =C17+K17*SIN(D17+RADIANS(M17))=B18+J18*COS(D18+RADIANS(L18)+PI()) =C18+J18*SIN(D18+RADIANS(L18)+PI()) =B18+K18*COS(D18+RADIANS(M18)) =C18+K18*SIN(D18+RADIANS(M18))=B19+J19*COS(D19+RADIANS(L19)+PI()) =C19+J19*SIN(D19+RADIANS(L19)+PI()) =B19+K19*COS(D19+RADIANS(M19)) =C19+K19*SIN(D19+RADIANS(M19))=B20+J20*COS(D20+RADIANS(L20)+PI()) =C20+J20*SIN(D20+RADIANS(L20)+PI()) =B20+K20*COS(D20+RADIANS(M20)) =C20+K20*SIN(D20+RADIANS(M20))=B21+J21*COS(D21+RADIANS(L21)+PI()) =C21+J21*SIN(D21+RADIANS(L21)+PI()) =B21+K21*COS(D21+RADIANS(M21)) =C21+K21*SIN(D21+RADIANS(M21))=B22+J22*COS(D22+RADIANS(L22)+PI()) =C22+J22*SIN(D22+RADIANS(L22)+PI()) =B22+K22*COS(D22+RADIANS(M22)) =C22+K22*SIN(D22+RADIANS(M22))=B23+J23*COS(D23+RADIANS(L23)+PI()) =C23+J23*SIN(D23+RADIANS(L23)+PI()) =B23+K23*COS(D23+RADIANS(M23)) =C23+K23*SIN(D23+RADIANS(M23))呵呵，这个就是传说已久的EXCEL自动计算曲线坐标．给大家研究研究．欢迎在此基础上开发增加新模块。
复曲线(缓和曲线、圆曲线)坐标自动计算
X
曲线起点ZH 坐标及桩号
Y
桩号（m) 桩号（m)
X Y
与路线方向 右夹角 与中桩 垂直距离 计算 方位角
V( ) S(m)
。
1.570796327 10.000 -0.2671 9.645 -2.640 -9.645 2.640
248.614 239.527 229.324 244.806
1057.505 610.551
交点JD转角α(
。
)
-120.000 1000.000 30.000
α1 ( )
第二圆曲线半径R2(m) 第二缓和曲线长ls2(m)
。
#
70.000
第一圆曲线半径R1(m) 第一缓和曲线长ls1(m)
300.000 100.000 0
0.0548 98.508 49.269 49.954 1.389 0.16667 162.486
(0.0309) ###### ######
待测点桩号(m) 待测点 中桩计 算坐标 边桩 已知JD桩号(m) 内插线 段长 lL(m) --计算 缓曲线
80.000 98.627
0.000 15.000 0.038 0.0150 1206.730
内插线段 调整角θ或偏角△ 内插线段弦长(m) 内插缓曲时切线 调整长度 第二 缓和曲线增值q 2 第二圆曲线内移值p 2 第二缓和曲线角β02 第二圆曲线长Ly2
内插直线时切线 调整长度 第一 缓和曲线增值q 1 第一圆曲线内移值p 1 第一缓和曲线角β01 第一圆曲线长Ly1
X
曲线终点HZ 坐标及桩号
56.166 32.428 1995.320 1900.000 -243.718 1361.813 397.477 1708.231
方位角及坐标计算
方位角及坐标计算公路工程各点方位角及坐标计算公式（一）各点方位角计算：1、第一直线段（k0～zh）：f=arctgδy/δx备注：直线方位角必须考量象限角就可以厘定恰当线路迈向2、第一缓解曲线段（kzh～khy）：δ1=（k0－kzh）2/（2rlh）×180/π3、圆曲线段（khy～kyh）：δ2=[2（k0－kzh）－lh]/2r×180/πδ2=（khy－kzh）/2r×180/π＋（k0－khy）/r×180/π无缓和曲线时：δ2=（k0－khy）/r×180/π（即圆曲线圆心角）4、第二缓和曲线段（kyh～khz）：δ3=（khz－k0）2/（2rlh）×180/π5、第二直线段（khz～kzh）：f±α（左偏时f－α，右偏时f+α）备注：k0――排序点的程α――曲线交点偏角lh――缓和曲线长（注意有时第一和第二缓和曲线长不一样）（二）各点坐标计算xzh=xjd－t?cosfxhz=xjd＋t?cos（f±α）yzh=yjd－t?sinfyhz=yjd＋t?sin（f±α）1、第一直线段：x=xzh＋（k0－kzh）?cosf中桩y=yzh＋（k0－kzh）?sinfx边=x中±b?cos（f－δ）边桩y边=y中±b?sin（f－δ）备注：b――中桩至所求点的距离（左幅时为＋b，右幅时为－b，当设计轴线与线路不横向时b取斜短，即b/sinδ）设计轴线线路方向。
bδ图s-12、第一缓和曲线段：xx=xzh－y′?sinθ＋x′?cosθxx′x′中桩′y=yzh＋y′?cosθ＋x′?sinθyzhyθhzx边=x中±b?cos（f＋μδ1－δ）hyyh边桩y边=y中±b?sin（f＋μδ1－δ）jdy′注：（本公式只适用与图s-2线形）图s-2μ――曲线左转为－1，右转为＋1θ――线路方位角与y轴所缠的锐角，见到图s-2y′=l－l5／（40r2lh2）；x′=l3／（6rlh）－l7／（336r3lh3）；（r―圆曲线半径，l―缓解曲线就任一点至曲线起点长度）3、圆曲线段：x=xhy+2r?sinφ?cos（f+μ（ξ+φ））中桩y=yhy+2r?sinφ?s in（f+μ（ξ+φ））x边=x中±b?cos（f＋μδ2－δ）边桩y边=y中±b?sin（f＋μδ2－δ）备注：φ=（k0－khy）/2r×180/π；ξ=（khy－kzh）/2r×180/π4、第二缓解曲线段：x=xhz－y′?sinθ＋x′?cosθ中桩y=yhz－y′?cosθ－x′?sinθx边=x中±b?cos（f＋μδ1－δ）边桩y边=y中±b?sin（f＋μδ1－δ）注：1、本公式只适用于与图s-2线形，其他线形可以根据本线形公式转换2、式中符号与第一缓解曲线意义相同3、注意有时第一缓和曲线长和第二缓和曲线长不一样4、第二直线段：x=xhz+（k0－khz）?cos（f±α）中桩y=yhz＋（k0－khz）?sin（f±α）x边=x中±b?cos（f±α－δ）边桩y边=y中±b?sin（f±α－δ）备注：f――第一直线段的方位角（三）用casiofx-4500p计算已知坐标点在线路上的里程和距中线距离1、直线段（已知坐标x、y）pol(x-xhz，y-yhz)：k=v?cos（f－w）＋khzb=v?sin（f－w）备注：1、在fx-4500p中计算结果取走变量储存区v和w，必须表明储存区内容时按rclv、w键。
excel绶和曲线计算公式
excel绶和曲线计算公式文字1.应用程序的目标计算最为复杂的诗路中线逐点坐标和切线方位角，既中桩的x值、y值。
excel的主要目标就是计算公路中线的中桩参数。
2．源程序l1 f：n：g：r：u：q：j：k起算要素l2 defm 4 扩大内存l3 m=i/2-i^3/240/r^2 m为加设缓和曲线后使切线增长的距离l4 p=i^2/24/r-i^4/2688/r^3 p为加设缓和曲线，圆曲线相对于切线的内移量l5 e=（r+p）/cos(n/2)-r e外矢距在excel中角度默认为弧度。
l6 l=πrn/180+i l为曲线长含缓和段l7 t=m+（r+p）tan（n/2）切线长l8 a=q-t：b=a+i：d=a+l：c=d-i zh：hy：hz：yh里程桩号l9 rec（t，f+180）zh点坐标增量计算l10 z[1]=v+j：z[2]=w+k zh点坐标（v，w坐标增量，计算机内部运算）l11 rec（t，f+gn）hz点坐标增量计算l12 z[3]=v+j：z[4]=w+k hz点坐标lbi θl14 {h，s，e} 待点要素l15 h<a＝＞goto 1 第一段直线上任一点坐标计算l16 ≠＞h<b＝＞goto 2 第一段缓和曲线上任一点坐标计算l17 ≠＞h<c＝＞goto 3 圆曲线上任一点坐标计算l18 ≠＞h<d＝＞goto 4 第二段缓和曲线上任一点坐标计算l19 ≠＞goto 5第二段直线上任一点坐标计算l20 lbi i 第一段直线坐标计算开始l21 rec（q-h，f+180）中桩坐标增量计算l22 x=v+j：y=w+k 中桩坐标l23 x=x+v：“x”=◢l24 y=y+w：“y”=◢第一段直线上任一点坐标计算结果l25 goto θl26 lbi 2第一段缓和曲线坐标计算开始l27 z=h-a所求点到zh距离l28 o=90z^2/r/i/π 所求点的方位角l29 x=z-z^5/40/r^2/i^2+z^9/3456/r^4/i^4l30 z=z^6/6/r/i-z^7/336/r^3/i^3+z^11/42440/r^5/i^5 第一段缓和曲线上任一点切线支距法坐标计算l31 lbi 6 第一段缓和曲线、圆曲线坐标换算计算l32 rec（x，f）l33 x=z[1]+v：y=z[2]+wl34 rec（x，f+90g）l35 x=x+v：y=y+wl36 lbi 3 圆曲线上任一点切线支距法计算开始l37 z=h-a-i/2l38 o=180z/r/π z所对中心角l39 x=z-z^3/6/r^2+z^5/120/r^4+ml40 z=z^2/2/r-z^4/r^3/24+z^6/720/r^5+pl41 goto 6l42 lbi 4 第二段缓和曲线计算开始l43 z=d-h到hz点长度l44 o=90z^2/r/i/π z所对中心角l45 x=z-z^5/40/r^2/i^2+z^9/3456/r^4/i4l46 z=z^3/6/r/i-z^7/336/r^3/i^3+z^11/42240/r^5/i^5l47 rec（x，f+gn+180）l48 x=z[3]+v：y=z[4]+wl49 rec（x，f+gn+180-90g）l50 x=x+v：y=y+wl51 x=x+v：“x”=◢l52 y=y+v：“y”=◢l53 goto θ程序中符号说明f-第一直线段正方位角（即后视切线上任一点（如zh）至交点的方位角；n—交点转向角（即偏角）；右偏角为正，左偏角为负，输入转向角时不得输入“+、-”g—条件（左偏输入-1，右偏输入+1）；r—半径；i—缓和曲线长度（m）（圆曲线i输入0）；q—交点里程桩号；j—交点x坐标值；k—交点y坐标值；h—所求点里程桩号；x、y—所求点x、y坐标值。
excel求曲线坐标
1．新建一个工作薄，在其中输入如表1所示的内容；2. 选中工作表A列，打开格式菜单，选中“单元格”，在单元格菜单中选中“数字”栏，自定义单元格格式为“K000＋000.000”。
按此方法分别将其他列设置为如表1所示单元格格式；3. 将“4”行做为路线起点数据行，在“5”行中进行公式编辑；4．在“J5”单元格中输入“=IF(C5=4,RADIANS(IF((G5+H5/60+I5/60/60)<180,(G5+H5/60+I5/60/60)+180,(G 5+H5/60+I5/60/60)-180)),IF(C5=5,RADIANS(IF(B5=0,G5+H5/60+I5/60/60-E5/2/D 5*180/PI(),G5+H5/60+I5/60/60+E5/2/D5*180/PI())),RADIANS(G5+H5/60+I5/60/60 6．在“L5”单元格中输入“=$L$4+IF(C5=1,K5*COS(J5),IF(C5=2,2*D5*SIN(K5/2/D5)*IF(B5=0,COS(J5-K5/ 2/D5),COS(J5+K5/2/D5)),IF(C5=3,(K5-K5^5/40/D5^2/E5^2)*COS(J5)+(K5^3/6/D5/ E5-K5^7/336/D5^3/E5^3)*IF(B5=0,COS(J5-90*PI()/180),COS(J5+90*PI()/180)),IF( C5=4,(K5-K5^5/40/D5^2/E5^2)*COS(J5)+(K5^3/6/D5/E5-K5^7/336/D5^3/E5^3)*IF (B5=0,COS(J5-90*PI()/180),COS(J5+90*PI()/180)),2*D5*SIN(K5/2/D5)*IF(B5=0,C 7．在“M5”单元格中输入“=$M$4+IF(C5=1,K5*SIN(J5),IF(C5=2,2*D5*SIN(K5/2/D5)*IF(B5=0,SIN(J5-K5/2 /D5),SIN(J5+K5/2/D5)),IF(C5=3,(K5-K5^5/40/D5^2/E5^2)*SIN(J5)+(K5^3/6/D5/E5 -K5^7/336/D5^3/E5^3)*IF(B5=0,SIN(J5-90*PI()/180),SIN(J5+90*PI()/180)),IF(C5= 4,(K5-K5^5/40/D5^2/E5^2)*SIN(J5)+(K5^3/6/D5/E5-K5^7/336/D5^3/E5^3)*IF(B5= 0,SIN(J5-90*PI()/180),SIN(J5+90*PI()/180)),2*D5*SIN(K5/2/D5)*IF(B5=0,SIN(J5- 8．在“N5”单元格中输入“=IF(B5=0,J5-RADIANS(IF(C5=2,K5/D5*180/PI(),IF(C5=3,K5^2/2/D5/E5*180/PI( ),IF(C5=4,K5^2/2/D5/E5*180/PI(),IF(C5=5,K5/D5*180/PI()))))),J5+RADIANS(IF(C 5=2,K5/D5*180/PI(),IF(C5=3,K5^2/2/D5/E5*180/PI(),IF(C5=4,K5^2/2/D5/E5*180/ PI(),IF(C5=5,K5/D5*180/PI()))))))”；9．在“P5”单元格中输入“=IF(C5=4,L5+O5*COS(N5+90*PI()/180),L5+O5*COS(N5-90*PI()/180))”；10．在“Q5”单元格中输入“=IF(C5=4,M5+O5*SIN(N5+90*PI()/180),M5+O5*SIN(N5-90*PI()/180))”；11．在“S5”单元格中输入“=IF(C5=4,L5+O5*COS(N5-90*PI()/180),L5+O5*COS(N5+90*PI()/180))”；12．在“T5”单元格中输入“=IF(C5=4,M5+O5*SIN(N5-90*PI()/180),M5+O5*SIN(N5+90*PI()/180))”；13．选定“5行”向下复制到表格最后一行；14．点击“文件”菜单中的“另存为”选项，输入文件名为“坐标计算”在文件类型中选择为“模版”并点击“保存”。
利用ExcelVBA程序计算曲线坐标.kdh
利用Excel VBA 程序计算曲线坐标张卫庭高强（中铁济南勘察设计咨询院有限公司山东济南250022）【摘要】Excel 是大家很熟悉的办公软件，相信大家在工作中经常使用。
本文利用Excel VBA 编制程序计算铁路缓和曲线的坐标，并将计算坐标转化为任意平面坐标系下的坐标，最后转化全站仪可识别的文件格式，从而实现了曲线坐标计算和现场放样的一体化。
【关键词】Excel VBA ；曲线要素；坐标转化；全站仪1.曲线坐标计算1.1角度转化为弧度计算曲线坐标，难免要把角度转化为EXCEL所识别的弧度来计算，为此首先要做的是角度转化为弧度。
打开Excel ，按Alt+F11即进入VBAIDE ，在菜单上依次点击[插入]->[模块]，然后输入如下代码Public Const pi =3.14159265359Public Function DEG(n As Double)Dim A As Double,B As Double,C As Double,D As Double,E As Double,F As Double,G As Double,KA As DoubleD =Abs(n)+0.000000000000001F =Sgn(n)A =Int(D)B =Int((D -A)*100)C =D -A -B /100DEG =F *(A +B /60+C /0.36)*pi /180End Function1.2曲线坐标计算1.2.1计算曲线坐标的公式加入缓和曲线后的铁路曲线示意图（见右图）坐标系的建立主要取决于控制点的情况。
如果控制点是为测设曲线而布设的，则坐标系一般采用ZH —XY 坐标系统①ZH~HY 段曲线点的坐标根据缓和曲线方程计算：x A =l A -l A 540R 2l 02y A =±l A 36Rl 0-l A7336R 3l 03!"#%%%%%$%%%%%&当曲线右偏时y A 坐标为正，左偏时y A 坐标为负。
直线 曲线及转角标准计算表-道路工程Excel计算系列
点
交点桩号
号
交点坐标
x
y
转角 度° 分′ 秒″
弧度
半径
缓和曲线 缓和曲线 长度L1 长度L2
△R
曲线要素值
q
切线长度T!
切线长度 T2
曲线长度
外距
校正值
第一缓和曲 线起点
曲线
1
2
3
4
5
67
8
起点 0
3893278.707 380915.695
9
10
11
12
13
14
15
16
17
18
19
0.00
1
175.75 3893447.717 380867.495 21 30 16
直线长度 (m)
交点间距 (m)
计算方位角
20
21
22
23
24
25
26
27
夹角(度)
28
转角(度)
29
30
153.6666 174.9651 196.2636 266.2636 83.66661 175.7487 -0.27782 344.0824 21.50468
右
21.50468
#DIV/0!
#DIV/0!
曲线位置第一缓和曲线终点第二缓和曲线起点第二缓和直线长度直线长度及方向交点间距曲线中点计算方位角夹角度转角度或圆曲线或圆曲线曲线终点mm起点终点202122232425262728291749651196263626626368366661175748702778234408242150468右div0
3
#DIV/0!
###### ##### ####