隧道线元法坐标计算excel
测量坐标计算EXCEL技巧（适用范围交通、建筑、水利、电力等土木工程专业）
测量坐标计算EXCEL技巧（适用范围交通、建筑、水利、电力等土木工程专业）测量坐标计算EXCEL技巧(适用范围交通、建筑、水利、电力等土木工程专业）内容摘要：[ft=#0000ff,,楷体_GB2312]结合公路工程实践，充分利用Excel电子表格中的公式与函数功能，使用简单易懂的嵌套公式的方法，制作公路工程中线逐桩坐标计算的文档模版，应用于公路中线一般线型的坐标计算。
关键词：Excel 公路中线坐标计算现阶段我国公路工程中已普遍使用大地坐标进行线型的控制及测设，在施工中经常要对中线坐标进行复核、加密，才能满足公路工程施工的需要。
本文是结合公路工程的实际需要，利用Exce电子表格制作的用于由直线、圆曲线、缓和曲线组成的一般公路线型中桩、边桩等计算的通用模版，用于减轻计算工作的劳动强度和提高计算结果的准确度。
一、采用公式1 直线段1.1 中桩坐标计算公式1.2 边桩坐标计算公式2 缓和曲线段2.1 中桩坐标计算公式当P点位于顺时针方向时,其方位角为αE→p=αA→B+90°；当P 点位于逆时针方向时,其方位角为αE→p=αA→B-90°。
2.2 边桩坐标计算公式3 圆曲线段3.1 中桩坐标计算公式当E点位于顺时针方向时取“+”,当E点位于逆时针方向时取“-”。
3.2 边桩坐标计算公式XP、YP——未知点P的坐标X1、Y1——各线型起点的坐标（第二曲线段为终点）XA、YA、XB、YB——P点边桩A点、B点的坐标（A为左侧、B 为右侧）α1→2——直线段起点的方位角αA→B——各线形起点的切线方位角（第二曲线段为终点）L——P点距各线形起点的长度LS——缓和曲线段缓和曲线长R——各曲线段的半径β——P点的切线角（曲线左转时取“-”、曲线右转时取“+”）T1、T2——P点至边桩A、B的距离（A为T1、B为T2）二、计算模板的建立步骤表1－1表1－21．新建一个工作薄，在其中输入如表1所示的内容；2. 选中工作表A列，打开格式菜单，选中“单元格”，在单元格菜单中选中“数字”栏，自定义单元格格式为“K000＋000.000”。
利用EXCEL表格计算线路中边桩坐标
利用EXCEL表格计算线路中边桩坐标高速公路施工中，桥梁、隧道施工以及路沿石施工时对测量要求相当严格，要求总体宽度误差在1cm以内。
这就要求测量工程师必须把线路边桩都计算出来，每隔10米（曲线段）或者20米（直线段）放样线路边桩，以便指导施工。
如果采用常规的计算器计算，不仅繁琐、费力而且容易出现差错。
本人在计算中边桩坐标时，试着利用EXCEL表格功能，编辑函数大批量计算线路中边桩坐标，达到了高效、准确的目的。
由于直线部分相对简单，复曲线又特别繁琐。
现就圆曲线举例计算如下：例：在浙江省龙丽一级公路施工中，左线圆曲线起点里程为K86+966.6，曲线半径为1500m，曲线左偏，起始方位角为147°43 ′58.2″，线路中心和隧道中心偏差15cm。
计算隧道中心线和衬砌边线（半径5.1m）。
1.在B4方格中输入：=$B$3+3000*SIN((A4-$A$3)/3000)*COS(2.5784244-(A4-$A$3)/3000)+ 0.15*COS(4.149220679-2*(A4-$A$3)/3000)按回车键确认即可以计算出该里程隧道中心线X轴坐标。
2.在C4方格中输入：=$C$3+3000*SIN((A4-$A$3)/3000)*SIN(2.5784244-(A4-$A$3)/3000)+0.15*SIN(4.149220679-2*(A4-$A$3)/3000)按回车键确认即可以计算出该里程隧道中心线Y轴坐标。
同理在3.D4方格中输入：=$B$3+3000*SIN((A4-$A$3)/3000)*COS(2.5784244-(A4-$A$3)/3000)+ COS(4.149220679-2*(A4-$A$3)/3000) *（0.15+5.1）4.在E4方格中输入：=$C$3+3000*SIN((A4-$A$3)/3000)*SIN(2.5784244-(A4-$A$3)/3000)+SIN(4.149220679-2*(A4-$A$3)/3000)*（0.15+5.1）其中 2.5784244为方位角147°43 ′58.2″的弧度表达方式，4.149220679为起始方位角147°43 ′58.2″加90°后的弧度值。
交点法、线元法坐标计算
3、交点法、线元法坐标计算坐标计算是根据图纸中“直线及曲线转角一览表”提供的数据计算道路中桩坐标，然后和图纸提供的“逐桩坐标表”比对，如果一样则说明输入平曲线参数输入正确，可以计算边桩坐标和其他结构物坐标了；如果中桩坐标不一样，一般是平曲线参数输入有误，需要重新检查输入，另一种结果是图纸有错，这种情况少见，但不代表没有。
“直线及曲线转角一览表”和“逐桩坐标表”见附件1、附件2。
线元法是以路线的起点坐标、方位角、起终点桩号等节点元素来计算出要求的坐标；交点法是以路线的交点要素和路线的主要要素来求得坐标。
①交点法交点：路线的转折点，路线改变方向是相邻两直线的延长线相交的点。
用JD表示，有些图纸上用IP 表示。
看下图：交点是针对曲线的（包含圆曲线和缓和曲线），一段曲线就有一个交点。
交点参数有：坐标（X,Y）、交点桩号、转角值、圆曲线半径R、缓和曲线长度。
教学提供软件（轻松测量、双心软件、测量工具）交点法曲线要素输入说明：1、QD起点坐标：起点坐标必须在直线段上，或填写前一交点的坐标。
2、JD交点曲线要素：（1）交点桩号（2）交点坐标（X,Y）（3）曲线半径R（4）第一缓和曲线长度LS1,若为0,输入0,不能为空。
（5）第二缓和曲线长度LS2,若为0,输入0,不能为空。
3、ZD终点坐标：终点坐标也必须在直线段上，或填写后一交点的坐标。
检核数据是否输入正确的方法：软件生成的圆曲线要素中切线长、外距、交点里程：注意校正起点里程、等与设计图纸是否一致。
如果上述数据和图纸不一样，请认真检查有错误的交点处的数据输入是否正确，如果输入没有错误，请考虑是否包含不完整缓和曲线，使用公式A2=R*Ls检查是否包含不完整缓和曲线。
如果包含不完整缓和曲线，那就需要用线元法也叫积木法计算了。
有的设计院给出的直曲表是整条设计线路的直曲表的一部分，以其中某个交点作为起始点的话，起始里程有时候需要校正，当然，并不是每个图纸给出的起点里程都需要校正，大多数图纸的起点里程已经被设计院校正过，我们输入平曲线的时候需要验证一下。
EXCLE函数和隧道工程应用实例
2019/9/18
5
3、逻辑函数
• 本类我们重点看一下条件选择函数的用法。
• IF(logical_test,value_if_true,value_if_false) • 其中Logical_test 表示计算结果为 TRUE 或 FALSE 的任意值或
表达式。此参数可使用任何比较运算符。 • Value_if_true 是 logical_test 为 TRUE 时返回的值或表达式 • Value_if_false 是 logical_test 为 FALSE 时返回的值或表达式。 • 此函数可以嵌套使用，最多可以嵌套64层。
• 时间值为日期值的一部分，并用十进制数来表示（例如 12:00 PM 可表示为 0.5，因为此时是一天的一半）。
• DATE(year,month,day) 返回特定日期的序列号 • 如果 year 位于 0（零）到 1899（包含）之间，则 Excel 会
将该值加上 1900，再计算年份。例如：DATE(108,1,2) 将返 回 2008 年 1 月 2 日 (1900+108)。 • 如果 year 位于 1900 到 9999（包含）之间，则 Excel 将使 用该数值作为年份。例如：DATE(2008,1,2) 将返回 2008 年 1 月 2 日。 • 如果 year 小于 0 或大于等于 10000，则 Excel 将返回错误 值 #NUM!。
表中两个竖曲线把线路设计标高分成了5个不同的区段，每个区段 的曲线要素和算法都有不同，这给我们现场计算带来很多麻烦， 我们拿到一个里程，先好看它位于哪个区段，在去查这个区段 的曲线参数，然后定算法，这样就比较麻烦。
五个区段的具体参数和算法如下表所示：
2019/9/18
坐标计算表格(1)
单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段549054955500单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段5505551055155520552555305535554055455550555555605565557055755580单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段5585559055955600560556105615562056255630563556405645565056555660单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段5665567056755680568556905695570057055710571557205725573057355740单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段5745575057555760576557705775578057855790579558005805581058155820单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段5825583058355840584558505855586058655870587558805885589058955900单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段5905591059155920592559305935594059455950595559605965597059755980单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段59855987.74459905995600060056010601560206025603060356040604560506055单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段6060606560706075608060856090609561006105611061156120612561306135单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段61406142.74461456150615561606165617061756180618561906195620062056210单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段6215622062256230623562406245625062556260626562706275628062856290单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段629563006305631063156320632563306335634063456346.1636350635563606365单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段6370637563806385639063956400640564106415642064256430643564406445单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段6446.163645064556460646564706475648064856490649565006505651065156520单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段6525653065356540654565506555656065656570657565806585659065956600单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段6605661066156620662566306635664066456650665566606665667066756680单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段66856690669567006705671067156720672567306732.58967356740674567506755单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段6760676567706775678067856790679568006805681068156820682568306832.589单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段6835684068456850685568606865687068756880688568906895690069056910单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段6915692069256930693569406945695069556960696569706975698069856990单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段6995700070057010701570207025703070357040704570507055706070657070单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段7075708070857090709571007105711071157120712571307135714071457150单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段7155716071657170717571807185719071957200720572107215722072257230单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段7235724072457250725572607265727072757280728572907295730073057310单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段7315732073257330733573407345735073557360736573707375738073857390单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段7395740074057410741574207425743074357440744574507455746074657470单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段7475748074857490749575007505751075157520752575307535754075457550单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段7555756075657570757575807585759075957600760576107615762076257630单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段7635764076457650765576607665767076757680768576907695770077057710单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段7715772077257730773577407745775077557760776577707775778077857790单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段7795780078057810781578207825783078357840784578507855786078657870单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段7875788078857890789579007905791079157920792579307935794079457950单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段7955796079657970797579807985799079958000800580108015802080258030单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段8035804080458050805580608065807080758080808580908095810081058110单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段8115812081258130813581408145815081558160816581708175818081858190单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段819582008205821082158220822582308235824082458250825582608264.9818265单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段8270827582808285829082958300830583108315832083258330833583408345单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段8350835583608365837083758380838583908395840084058410841584208425单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段8430843584408445845084558460846584708475848084858490849585008505单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段8510851585208525853085358540854585508555856085658570857585808585单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段8590859586008605861086158620862586308635864086458650865586608665单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段8670867586808685869086958700870587108715872087258730873587408745单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段8750875587608765877087758780878587908795880088058810881588208825单位名称：工程名称：坐 标 计 算 表线元法**集团有限公司****路线**标段8830883588408845885088558860886588708875888088858890889589008905。
公路路线不同曲线线元坐标计算方法
公路路线不同曲线线元坐标计算方法李自康【摘要】由构成公路工程（或铁路等土木工程）路线的不同线元——直线、圆曲线、缓和曲线（包括卵形曲线）的基本特征出发，推导出该线元上任意点的坐标及切线方位角计算方法，并针对公路工程常见基本型曲线各要素计算方法及各主点坐标及切线方位角计算方法进行了推导，以期为全路线坐标计算提供计算依据。
【期刊名称】《交通世界（运输车辆）》【年(卷),期】2016(000)005【总页数】3页(P19-21)【关键词】直线;圆曲线;缓和曲线;坐标;方位角【作者】李自康【作者单位】贵阳市城市发展投资集团股份有限公司，贵州贵阳 550001【正文语种】中文【中图分类】U412.3针对公路工程（或铁路等土木工程）路线设计来说，无论曲线设计多复杂，无外乎是由直线、圆曲线、缓和曲线按照一定的组合方式连接而成。
对于测量人员来说，就是要将设计人员已经确定的路线在施工现场还原，将各种构造物、路基各控制点在施工现场通过各种放样方法确定下来，然后再指导施工人员按要求施工形成路基、桥涵、隧道、路面等等。
因此，如何通过简单、快捷的方法对路线进行放样至关重要。
在教科书中常见的有切线支距法、弦线支距法、偏角法等方法，计算相对繁琐，现场放样也不太准确。
在计算机技术日益普及、全站仪应用广泛的今天，坐标法大量应用于公路曲线坐标计算及放样当中，大大提高了计算和放样的速度，同时也大大提高了放样精度。
笔者就构成公路路线的不同线元（直线、圆曲线、缓和曲线）的坐标计算方法进行推导并给出相应的坐标及方位角计算公式，在实际的放样过程中，对于拟放样坐标点，根据该点桩号判断出该点所处的线元类型，然后根据该线元公式进行计算即可得出待求点坐标及切线方位角，即可快速地通过全站仪或GPS测量等方法放样到施工现场。
已知直线起点坐标X0、Y0及方位角α0，直线上任一点至起点的距离为D，则该点的坐标X、Y以及切线方位角α的计算公式为（参见图1）。
线元法任意线型坐标计算
边桩பைடு நூலகம்
N (m) #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? E (m) #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME?
坐标计算表
中桩
方位角 #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? 偏距 (m) -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 -12.25 夹角 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90
备注
公司: 项目:
计算里程 (m)
K209+580.000 K209+600.000 K209+620.000 K209+640.000 K209+660.000 K209+680.000 K209+700.000 K209+720.000 K209+740.000 K209+760.000 K209+780.000 K209+800.000 K209+820.000 K209+840.000 K209+860.000 K209+880.000 K209+900.000 K209+920.000 K209+940.000 K209+960.000 K209+980.000 K210+000.000 N (m) #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? E (m) #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME? #NAME?
Excel在公路线路坐标计算中的应用
l 9/
(
3
456R
4l
4 0
)
-
l
13/
(
599
040R
6l
6 0
)
( 1)
y1 = l3/ ( 6Rl 0) - l 7/ ( 336R 3l 03)
+ l 11/ ( 42 240R5l 05)
( 2)
m=
l0-
l
3 0
/
(
40R
2)
+
l 05/ ( 3 456R 4) -
l
7 0
/
(
599
040R 6) - R sin( l 0/ 2R )
角
左侧边桩 坐标增量
xy
KL
左侧边 桩坐标
XY
2 A2
D2 E2 F2
3
4
5
第二行输入 HY 点对应的里程A 2、到 HY 点的线 路中线方位角 D2 及 H Y 点坐标( E2, F2) , C 列 C3 = B3/ 2R , 其他计算同前面中桩、边桩计算。
应用 Excel 计算时, 由于 Excel 表格能够满足中桩 及左右边桩一次计算需要, 计算时可一次完成。在计 算熟练后, 可找出一些技巧, 由于本文篇幅关系, 不在 多述。
2. 2 曲线段计算
( 1) 中桩坐标计算 根据曲线段坐标计算思路, 应用 Excel 生成表 3。
表 3 曲线段中桩坐标计算
AB CD E F
G
H I JK
单曲线坐
中线坐标
1
里程
里程 差
标增量
弦长
反正 切值
线路中 线方位角
增量
xy
EXCEL电子表格在隧道轴线计算中的简单应用(修正)
EXCEL电子表格在隧道轴线计算中的简单应用张子敏一、前言在盾构施工中，隧道轴线数据一般采用可编程计算器进行计算，操作繁琐且效率不高，对于后期的数据整理造成困难，工作量较多。
本文通过将设计轴线计算和平面偏差计算的FX-5800P计算器运算程序转化成EXCEL电子表格形式，将计算器的运算语言转化成公式，其中结合相关函数的运用，可有效降低计算工作量，提高工作效率，便于后期数据整理。
二、设计轴线计算1、设计平面坐标计算（1）直线段Fx-5800P计算器程序：“N”?→N （N：计算里程）L 0-N→L或N-L→L沿计算方向里程变小或变大（L：直线段起始里程）ω→C （ω：直线段方位角）L·COS(C)+“X0”◢（X:已知起始里程坐标X值）L·SIN(C)+“Y0”（Y:已知起始里程坐标Y值）EXCEL电子表格：N——输入计算里程值；C——输入已知直线段方位角的弧度值；L——输入该直线段起始里程；L——输入公式：N-L0或L-N。
得到计算里程至起始里程的距离；X——输入公式：L·COS(C)+“起始里程坐标X值”，求得计算里程坐标X值；Y——输入公式：L·SIN(C)+“起始里程坐标Y值”，求得计算里程坐标Y值。
（2）缓和曲线段Fx-5800P计算器程序：“N”?→NL 0-N→L或N-L→L沿计算方向里程变小或变大f 1 (L)→E（f1(L):已知相对坐标X值方程）f 2 (E)→F（f2(E):已知相对坐标Y值方程）Pol(E,F)（生成I值：距离； J值：角度）ω±J→C 根据计算距离起始方向定“±”（ω:直线段方位角）I·COS(C)+ “X0”◢（X:已知起始里程坐标X值）I·SIN(C)+ “Y0”（Y:已知起始里程坐标Y值）EXCEL电子表格：N——输入计算里程值；L——输入该缓和曲线段起始里程；L——输入公式：N-L0或L-N。
交点法和线元法曲线要素输入简介
测量坐标计算程序V5输入简介本程序运用Office Excel 软件VBE标准模块编写，其功能基本全面集成了以往所更新的Excel程序，程序适用于公路、铁路等线路坐标计算，程序主要包括（交点法、线元法、直线坐标正反算，竖曲线计算，平面控制网“导线、高程”平差，隧道超欠挖，超高加宽，测量工具箱等，还可以全自动生成卡西欧5800、9750程序数据库，其中包括：隧道超欠挖、交点法、线元法、竖曲线一系列数据库），已知数据输入明确，操作简单易懂，是工程测量人员的好帮手！交点法曲线要素输入简介一、适用平曲线类型交点法计算坐标适用的平曲线为对称或不对称缓和曲线、圆曲线。
注意：对于非普通的三单元曲线，本程序交点法不适用。
非普通的三单元曲线体现在本程序中的《直线、曲线及转角表》内，点击“生成要素”之后，计算值与设计图纸《直线、曲线及转角表》上的切线长和曲线主点位置等不一致，此时只能采用线元法进行坐标计算。
例如：下表的JD18及JD19处的平曲线，经本程序交点法计算之后发现，为非普通的三单元曲线，交点法不适用该类曲线的坐标计算，故只能采用线元法进行坐标计算。
二、交点法曲线要素输入说明本程序交点法输入的要素有7个（程序不限制输入行数）：1、QD起点坐标：起点坐标必须在直线段上，或填写前一交点的坐标。
2、JD交点曲线要素：（1）交点桩号K，注意：当起始平曲线上的ZH点（缓和曲线）或ZY点（圆曲线）的桩号为负数时，交点桩号K统一加上100000（即增加100Km），以避免坐标正算时出现桩号计算范围错误（但是，线元法计算坐标时可以输入负坐标，坐标正算与反算都不会出现错误）。
（2）交点桩号（X,Y）（3）曲线半径R（4）第一缓和曲线长度LS1,若为0,输入0,不能为空。
（5）第二缓和曲线长度LS2,若为0,输入0,不能为空。
3、ZD终点坐标：终点坐标也必须在直线段上，或填写后一交点的坐标。
三、操作流程：1、根据设计图纸《直线、曲线及转角表》输入第一个交点坐标，作为QD起点坐标。
利用EXCEL进行道路中边线逐桩坐标计算
利用EXCE L 进行道路中边线逐桩坐标计算张宝利1,薛志勤2(1.西北农林科技大学资源与环境学院,陕西杨凌 712100;2.华南农业大学信息与软件学院,广州 510642)摘要:利用坐标增量的基本思想,将道路曲线坐标的计算归结为3种最基本的形式———直线段、圆曲线段、缓和曲线段,直接计算出中、边线的测量坐标,省去了切线支距法坐标向测量坐标的转换过程,简化了计算思路,节约了计算量。
并利用办公软件Excel 的强大计算功能进行实例计算。
关键词:Excel;道路边线;道路中线;逐桩坐标;圆曲线;缓和曲线中图分类号:U452.1文献标志码:B文章编号:1672-741X (2008)03-0302-03Calcul ati on of Coordi n ates of St akes on Centerli n e and Si deli n e Usi n g ExcelZHANG Baoli 1,XUE Zhiqin2(1.College of R esou rces and Environm en t,N orthw est A &F U niversity,Yangling 712100,Shanxi,Ch ina;2.College of Infor m ation and S oft w are,S outh China A gricultural U n iversity,Guangzhou 510642,China )Abstract:The coordinate calculati on of the high way curves is divided int o three basic types,i .e .,straight secti on,cir 2cular curve and transiti on curve,on basis of the coordinate incre ment p rinci p le .W ithout coordinate transf or mati on,the survey coordinates of the stakes on the alignment of the centerline and sideline are directly worked out,which si m p lifies the calculati on p rinci p le and reduces the calculati on quantity .Further more,coordinate calculati on is made f or p r oject cases by using the str ong calculati on functi ons of Excel .Key words:Excel;high way sideline;highway centerline;coordinates of stakes;circular curve;transiti on curve不论是在铁路还是公路的设计和建设中都不可避免地遇到大量的曲线计算。
9860线元法隧道测量(坐标正反计算、高程)计算程序
CASIO 9860G SD线元法隧道三维（坐标正反计算、高程）计算程序1. A（此为主程序）Lbl 0：〝1.LC=>XY〝：〝2.XY=>LC〝：〝3.ZHZL=>GC〝：〝PB=>V=1,2,3〞？→V：If V=1：Then GOTO 1 ：IfEnd ：If V=2：Then GOTO 2 ：IfEnd ：If V=3：Then GOTO 3 ：Else GOTO 0 ：IfEnd：Lbl 3：〝ZH=H〝？→H ：〝SDZF=Z〝？→Z：Prog〝ZGCZCX〝：GOTO 0：Lbl 1 ：〝ZH=L〝？→L：If L>173000 And L<174661.96：Then GOTO 4 ：Else GOTO 0 ：IfEnd ：Lbl 4：L→L：〝SDZF=Q〝？→Q：〝XLZJ,-Z+Y=Q〝：Q+0.125→Q：Prog 〝ZBQXYS〝：〝JSJD=J〝：90→J▲Prog 〝ZSZB〝：〝X=〝：X ▲ 〝Y=〝：Y▲〝FWJ=O〝：O▼DMS▲L→H：Q-0.125→Z：Prog 〝ZGCZCX〝：GOTO 0：Lbl 2：〝XO=M〝?→M：〝YO=R〝?→R：173300→L：If M>3845505.273 And M<3846506.099And R>499371.832 And R<500352.224 ：Then GOTO 5：Else GOTO 2：IfEnd ：Lbl 5：0→Q：0→J：Prog 〝ZBFS〝：〝LC=L〝：L ▲〝JL=Q〝：Q▲ 〝SDZJ，-Z+Y=Q〝：Q-0.125→Q ▲L→H：Q→Z：Prog 〝ZGCZCX〝：GOTO 02.正算坐标ZBZS( L-S ) / 4→H：90/π→F：HHF（1/T-1/I）/（K-S）→U：2HF/ I→D：C+4D+16 U→Ｏ：O+J→P ：C+ D+ U→E：C+2D+4U→W：C+3D+9U→G：A+AbsH/3*(cosC+4(cosG+cosE)+2cosW+cosO)+Qcos P→X ：B+AbsH/3*(sinC+4(sinG+sinE)+2sinW+sinO)+Qsin P→Y3.反算坐标：ZBFSLbl 0：Prog 〝ZBQXYS〝：Prog 〝ZBZS〝：O-90→Z：(R-Y)cosZ-(M-X) sinZ→P ：L+P→L：If Abs P≥0.001：Then GOTO 0 ：Else GOTO 1：IfEnd ：Lbl 1： (R-Y)cosO-(M-X) sinO→Q4. 曲线元要素数据库：ZBQXYSIf L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：If L≥S And L＜K ：Th en **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ：……………………………If L≥S And L＜K ：Then **→ S：**→ A：**→ B：**→ C：**→ I：**→ K：* *→ T IfEnd ┘(注：如有多个曲线元要素继续添加入数据库ZBQXYS中)5 高程计算主程序ZGCZCX (后有修改说明)Lbi1 ：〝SCGC=X〝？→X：〝R=M〝：5.98→M：〝CS=N〝：1.603→N：Prog＂GC SJK＂：C-D→F：Abs(RF÷2) →T：R AbsF÷F→R：If H≤B-T ：Then 0→K：GOTO 2：IfEnd ：If H>B-T And H<B+T ：Then H-B+T→K ：GOTO 2：IfEnd ：If H≥B+T ：Then 0→K：D→C：GOTO 2：IfEnd ：Lbi 2 ：〝XLZG=G〝：A-（B-H）C-K^2÷2R→G▲If Z≥0：Then If X>100 Then 〝YO1XGC=J〝：G+N→J▲〝YGCFSKD=J〝：Abs√(M^2-(X-（G+N）)^2) →J▲〝YKDCQ,+C,-Q=J〝：J- Abs (Z+0.000) →J▲〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J▲〝YGCCQ,C+,Q-=J〝：X-J→J▲Else〝YO1XGC=J〝：G+N→J▲〝YSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J▲ IfEnd ：Else 〝ZO1XGC=J〝：G+N→J▲If X>100 ：Then〝ZGCFSKD=J〝：Abs√(M^2-(X-（G+N）)^2) →J▲〝ZKDCQ,+C,-Q=J〝：J- Abs (Z+0.000) →J▲〝ZSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J▲〝ZGCCQ,C+,Q-=J〝：X-J→J▲Else〝ZSBSJGC=J〝：G+N+√(M^2-(Z+0.000)^2) →J ▲IfEnd ：IfEnd6高程计算主程序子程序：GCSJKIf H>起点桩号 And H≤第一个竖曲线起点桩号Then 第一竖曲线交点高程→A：第一竖曲线交点桩号→B：第一竖曲线前坡度→C：第一竖曲线后坡度→D：第一竖曲线半径→R：IfEnd：If H>第一竖曲线止点桩号 And H≤第二竖曲线起点桩号Then 第一竖曲线交点高程→A：第一竖曲线交点桩号→B：第一竖曲线前坡度→C：第一竖曲线后坡度→D：第一竖曲线半径→R：IfEnd：If H>第一竖曲线止点桩号 And H≤第二竖曲线起点桩号Then 第一竖曲线交点高程→A：第一竖曲线交点桩号→B：第一竖曲线前坡度→C：第一竖曲线后坡度→D：第一竖曲线半径→R：IfEnd………………………(继续添加要素)说明：第一部分坐标部分（1、2、3、4）V=1进入坐标正算 V=2进入坐标反算V=3进入单独的高程计算当V不等于1、2、3时，则返回程序，要求再次输入V值。
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自动计算曲线坐标．给大家研究研究．欢迎在此基础上开发增加新模块。
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—曲线半径α—方位角注:式中，紫色部分为缓和曲线任意点的坐标增量(支距坐标)。
Excel VBA自定义函数批量计算隧道线路坐标正算方法
•I l l-建 筑技术第51卷第7期2020年7月Architecture Technology_______________________Vol.51 No.7Jul.2020Excel V B A自定义函数批量计算隧道线路坐标正算方法谢江伟，贾云朋，孙青龙(北京建工土木工程有限公司，100015,北京）摘要：线路测量通常包括公路、随道、输电线等线状建筑物的测量放线工作，结合地铁隧道施工经验，运用Excel VBA二次幵发的基础上批量处理计算数据，编写自定义函数，将道路全线或部分线路所有数 据一次性输入到一张表格中，通过调用自定义函数调用参数自动识别相应里程点对应线元，自动计算坐标。
Excel函数方便填充序列，可实现批量计算线路坐标，也为后续隧道线路反算函数打下基础。
关键词：Excel VBA编程；线路正算；自动计算坐标；自定义函数；批量计算中图分类号：TU 745 文献标志码： A 文章编号：1000-4726(2020)07-0772-04EXCEL VBA USER DEFINED FUNCTION BATCH CALCULATION OF TUNNEL LINECOORDINATE FORWARD CALCULATION METHODXIE Jiang-wei,JIA Yun-peng,SUN Qing-long(Beijing Construction Engineering Civil Engineering Co.,Ltd., 100015, Beijing,China)Abstract: The line survey usually includes the survey and setting out of linear buildings such as highway, tunnel and transmission line. In combination with several years of subway tunnel construction experience, this paper use Excel on the basis of the secondary development of VBA, the interior data is processed in batches, and the user-defined functions are compiled. All the data of the whole road or a section are input into a table at one time. By calling the parameters of the user-defined functions, the corresponding line elements of the corresponding mileage points are automatically identified and the coordinates are automatically calculated. Because excel function is convenient to fill sequence, it can realize the batch calculation of line coordinates, and lay a foundation for the later part of tunnel line back calculation function.Keywords: Excel VBA programming; line direct calculation; automatic coordinate user; defined function; batch calculation对大多数从事隧道、道路工程测量的人员而言，放线仪器主要依靠全站仪或RTK,而带有线路转折 的曲线几乎是每条隧道、道路设计必不可少的一部分，对计算里程点坐标数据目前仍大量依赖专业测量软 件，而内业资料数据处理用的最多的还是Excel软件，若能将两者轻松结合则更加完美。