济业2014-46 既有铁路平面曲线拟合软件

1、既有铁路平面曲线拟合软件2.3介绍及使用说明既有铁路平面曲线拟合软件说明报告编制： 陈礼坤审核： 贾克永审定： 赵振刚共 25 页中铁工程设计咨询集团有限公司济南设计院2015年01月 济南1. 软件开发背景既有铁路经过长期运营后会发生不规则变形，影响行车速度和行车安全，为进行曲线整正，需要根据以各测点的拨距值的平方和最小或绝对值之和最小作为设计条件，拟合出最佳的平面曲线，计算出各个点的拨距量。现阶段采用的有偏角法和坐标法，其中偏角法源于渐伸线法计算原理，是以弦长代替弧长，具有一定的近似性误差且没有理论严密的计算公式，使得误差具有隐蔽性。并且偏角法在实施过程中，必须将仪器架设在既有线上，随着火

2、车运行速度和运行密度的提高，现场的作业条件已经很难满足这样的测量要求，这样就造成外业测量工作的难以实施和效率的低下。相对于传统的渐伸线法，坐标法具有理论严密、测量计算成果精度高等优点，既有的相关算法对于数据点的间距以及数据采集格式有较高的要求，部分软件没有曲线自动拟合功能。并且在对缓和曲线以及半径进行自动最佳匹配时，存在对半径和缓和曲线长初始值依赖性较高、收敛速度慢、适用范围窄的特点，现实中，也没有合适的免费软件可供使用，通过了解济南铁路局既有线测量队以及其他设计院在曲线拨距整正方面使用的软件和测量作业模式，可以考虑编写具有自主知识产权的平面曲线拟合软件，以适应测量形势的发展。这里有两张拨距施

3、工现场的图片，可以看出拨距施工的难度，那么怎样才能更好的拟合出最佳的平面曲线，获得各个点的最小拨距值，减少实际施工工作量，就显得尤为重要了。2. 软件主要功能2.1、软件可以根据既有铁路曲线上测量的点坐标（包含缓和曲线、圆曲线和曲线两端直线点坐标），通过自动模式或者手动模式，拟合出最佳平面曲线，得到每个点对应曲线上点的坐标、里程、拨距量，并给出拨距量是否合格的评判结果。既有线测量的位置可以自由选择（内、外轨以及线路中心）。2.2、可以生成CSV格式的既有曲线拨距成果表，画出DXF格式的曲线线位文件、以及曲线信息数据（线路HRCAD设计软件支持的dat格式）。2.3、支持缓和曲线长为0的曲线以及

4、缓和曲线不等长的曲线计算。2.4、如果导入的数据格式不正确或者没有正确选取直线边以及起始里程点，程序会给出错误提示。在软件自动模式参数设置里面，可以固定缓和曲线长，程序自动匹配最佳半径，也可以固定半径，程序自动匹配最佳缓长。2.5、支持设定单个点以及多个点的拨距限差值，程序在计算过程中，将按照设置的限差值进行曲线自动拟合以及质量结果评判，拟合出满足要求的曲线数据。3. 软件界面及操作按钮说明3.1 软件界面如下：左侧主要是参数设置区域，右侧表格区域是显示计算出的相关点的坐标、里程数据以及曲线的各个曲线要素，允许显示120行数据。现在先从左侧参数设置区域开始介绍。3.2 在左上角的参数设置程序框

5、内，有四个弹出按钮，点击，弹出如下对话框：项目相关信息将在曲线拨距成果表中显示。点击，弹出如下对话框：可以选择是50轨还是60轨，其中50轨的单轨顶面宽70mm ，60轨的单轨顶面宽73mm，43轨和50轨的单轨顶面宽式一样的。点击，弹出如下对话框： 可以选择你所测量的位置是内、外轨中心还是线路中心。点击按钮，点开后如下弹出如下对话框：可以在对话框里面根据点名指定两端固定的直线边和确定起始里程，也可以在导入的数据文件中，采用标识符来进行直线边固定和确定起始里程。如果同时进行了两种直线边设置方法，程序会优先以对话框的设置为准。在待导入的测量数据文件中，第一条直线边两个点在第四列以ZG1和ZG2表

6、示，第二条直线边两个点在第四列以ZG3和ZG4表示，里程在对应里程点的第四列以ZL加里程表示。具体格式见下表，如果起始里程点和ZG1是同一个点，那么将第一行的ZL179780改为“ZG1ZL179780”就可以了。120803.0819567.712ZL179780220788.86119581.809179800320774.64619595.823179820420760.49719609.868179840520746.19219623.932ZG1620732.02519637.946ZG2720717.7119652.04179900820703.52219666.099-100Z0

7、50920689.25919680.1541799405819802.07820006.245ZG35919782.12720005.5611809406319702.17820002.6341810206419682.24820001.9411810406519662.28620001.239ZG4 在上面的文件中的中，如果需要设定单个或者多个点的拨距限差值，例如：如果想让自动拟合出来的最佳曲线，在某个点的右拨距量不超过50mm。那么将第四列的点值设置成-100Z050；“Z”前面的-100为默认的左侧拨距限差值为100mm，Z后面的“050” 意思是自设的右侧拨距限差值为50mm。3.3在

8、模式选择程序框里，可单选计算模式，有手动模式和自动模式两种方式，其中默认是自动模式。若要切换成手动模式，请点击前面的。3.4 在自动模式程序框里面，第一个按钮，计算曲线的大偏角以及半径初始值，在没有台账，或者不知道此条曲线初始半径的情况下，可以点此按钮，计算曲线的初始半径以及曲线大偏角，根据初始半径值来设置“半径搜索范围”。需要按照中铁咨询新规则中的第二篇（测绘）中的P81页的取值要求来设置“半径搜索步长”。这样既可以满足半径取整要求，拟合出更合适的平面曲线，同时减少计算时间。在计算大半径曲线（半径>=10000m）以及小半径曲线（半径<=300m)曲线最好先进行此项计算。第二个按

9、钮，点开后弹出如下参数设置对话框：一般设置里面是的半径搜索步长是每次半径的更新值，缓长搜索范围最大支持500m。可以在高级设置中，固定半径，进行缓长拟合，也可以固定缓长进行半径拟合,对话框中的编辑框的数据均是默认的设置，一般情况下，默认的设置就可以满足我们的计算需要了。3.5 在表格的右下方是操作按钮，主要用于数据的导入、导出和计算处理的。其中：点此按钮，软件读入数据文件，文件格式如下：11580.7051999.609zl45036021600.6681999.027zg131620.6411998.41345040041640.6641997.82245042051660.6021997.

10、213zg261670.6371996.9445045071680.5761996.61645046081700.581996.03345048091720.5731995.423450500101740.571994.846 450520402350.4612083.034zg3412369.6742088.654451160422408.0152099.881451200432455.9722113.933zg4616219722.1620003.361810006319702.17820002.6341810206419682.24820001.9411810406519662.2862

11、0001.239ZG4 数据文件共有四列，分别对应着“点名，北坐标，东坐标，备注”。第一列和第四列支持任意字符，曲线上的坐标数据宜采用全站仪导线法测量获取。数据总数应不少于13行。其中圆曲线上至少3个点，如果曲线较短，应该多测几个点。 第一行和第二行以及最后两行的直线点是曲线两切线上的A、B、C、D点：示意图如下： A、B、C、D四个点必须分别选在曲线两端的直线上。AB、CD的长度以100200米为宜，靠近曲线的B、C两点必须设在曲线起、终点外方约4080米的直线上，如果夹直线很短（技规中规定不小于20米，但实际上有的曲线不能达到这个规定）直线不明显示，应在现场仔细判定，必要时可多测量几个点。

12、 文件后缀名为“XY”。文件支持“，”以及空格、制表符隔开的数据，每行之间允许有空行。：点此按钮，软件根据导入的数据，计算出曲线要素以及每个点对应的曲线上点的坐标、里程、拨距量等拨距计算成果。并在右侧显示出来，具体结果见下表，其中：第一列“点名”是导入的原始数据的点名；第二列“里程”是计算后测量点投影到拟合曲线上对应点的里程；第三列“北坐标(X)”和第四列“东坐标(Y)”是计算后测量点投影到拟合曲线上在曲线上的坐标值。第五列“内拨(mm)”和第六列 “外拨(mm)”表示的测量点到拟合曲线的距离，点在曲线外侧，需要往线路内侧移动，就是内拨，反之，就是外拨，第七列“是否超限(mm)”是拨距量是否合

13、格的评判标准，如果不合格，就会给出提示。第八列“曲线要素”是给出了拟合曲线的曲线成果表：其中第二行“Aa”表示曲线大偏角，左偏为负，右偏为正， 第三行“R”表示半径，第四行 “L1”、“ L2”分别代表第一缓长和第二缓长，第五行 “T1”和“T2”分别代表第一和第二切线长。第六行“l”表示曲线全长，第七到第十行“ZHLC”、“HYLC”、“YHLC”、“HZLC”分别代表直缓里程、缓圆里程、圆缓里程、缓直里程。第十一行“BoJu” 表示所有点的拨距值总和，第十三和第十四行数据时项目信息，其中“CurNO” “Pro” “Actor” “Check”分别表示曲线编号、项目名称、计算者、复核者，并

14、且一次设定，不论计算多少个曲线，只要不关闭软件，就一直是你设定的信息，保持项目设置信息的一致性。：点此按钮，将拨距成果保存成CSV数据格式文件，表格的数据内容和前面表格显示的结果一致，做到所见即所得，下面是保存计算结果后，直接打开的文件。：点此按钮，将拟合出来的最佳曲线保存成DXF数据格式文件，这种Autodesk公司的绘图交换文件可在AutoCAD里面直接打开，打开后的数据结果如下所示：点此按钮，保存成HRCAD线路设计软件所支持的.dat数据格式文件，提供此数据接口，方便与线路设计部门的数据交换。数据格式如下:：退出当前程序，和右上角的的功能一致，执行退出操作。4. 软件核心算法4.1在程

15、序设计过程中，关键是如何自动进行最佳曲线的拟合，为此采用改进的遗传算法，在自动匹配缓和曲线和半径的过程中，根据不同的计算阶段采用不同的适应度函数，兼顾计算效率和计算精度。主要分为下面俩个步骤：步骤一：采用最小二乘法以曲线中部一定范围的数据来初始化既有线路的半径R，对缓长L不进行初始化，下一步中将进行全面匹配（0m-500m），对观测精度较低或者实际线路变形较大的时候，现有的一些算法将在初始化缓长时会与实际缓长产生较大误差，在我的算法中对缓长不进行初始化，虽然增大了初始群体的数量，但避免对缓长的优化陷入局部最优的情况。提高了算法的适应性。同时根据半径R和搜索步长RB确定初始群体P(0)的个数M。

16、步骤二：利用GA算法来更新半径和缓和曲线长，其中半径的搜索步长既可以采取默认，也可以根据实际工作要求，通过提供的外部数据交换接口（编辑框）进行自定义，支持（0m-1000m）任意数据，在进行最佳的（R, L）匹配时，通过两代进化计算，以各测点的拨距值i 的平方和最小或绝对值之和最小作为目标函数，自动拟合出最佳的平面曲线，其中在每代的迭代计算中，采取不同的目标函数和计算方向，通过选择和变异算子，快速筛选出部分符合条件的设计曲线，在得到的下一代数据中，再根据精确算法，计算出每一点的拨距量，筛选出最佳半径和缓和曲线长的设计曲线。既兼顾效率和精度，同时由于在计算过程中，采用了逐点计算。还能满足单个控制

17、点对拨距量的特殊要求。4.2 主要算法框图如下所示：确定半径的搜步长和拨距限差并编码实现以半径R和缓长L组合的个体来构造初始化群体以5000-为适应度函数计算种群的适应度并寻找出群体中的最优值,即寻找出种群中满足适应度函数的最佳R和 L。是否获得满足适应度函数值最大的（R L）重新选择部分半径和缓长组合成的群体得到新种群并且采取精确的算法来计算每点拨距值NOYES结束5、编程实现5.1在拨距成果计算过程中，主要新增FReadCouSave类和LineCal类，其中FReadCouSave类负责数据的读入，LineCal类主要实现数据的计算，包含了多个成员函数，以实现不同的功能，其中部分代码如下

18、：if(pDlg->m_PanDclwz<0)/外轨m_x1=pDlg->m_x1-pDlg->m_PanDclwz*cos(m_JD1+zypj*3.141592654/2); m_x2=pDlg->m_x2-pDlg->m_PanDclwz*cos(m_JD1+zypj*3.141592654/2);m_y1=pDlg->m_y1-pDlg->m_PanDclwz*sin(m_JD1+zypj*3.141592654/2); m_y2=pDlg->m_y2-pDlg->m_PanDclwz*sin(m_JD1+zypj*3.141

19、592654/2); m_x3=pDlg->m_x3-pDlg->m_PanDclwz*cos(m_JD2+zypj*3.141592654/2); m_x4=pDlg->m_x4-pDlg->m_PanDclwz*cos(m_JD2+zypj*3.141592654/2); m_y3=pDlg->m_y3-pDlg->m_PanDclwz*sin(m_JD2+zypj*3.141592654/2); m_y4=pDlg->m_y4-pDlg->m_PanDclwz*sin(m_JD2+zypj*3.141592654/2);m_x1qs=ato

20、f(pDlg->Res4.GetAt(pDlg->m_lcdh)*4+1)-pDlg->m_PanDclwz*cos(m_JD1+zypj*3.141592654/2); m_y1qs=atof(pDlg->Res4.GetAt(pDlg->m_lcdh)*4+2)-pDlg->m_PanDclwz*sin(m_JD1+zypj*3.141592654/2);else if(pDlg->m_PanDclwz>0)/内轨m_x1=pDlg->m_x1+pDlg->m_PanDclwz*cos(m_JD1-zypj*3.141592654

21、/2); m_x2=pDlg->m_x2+pDlg->m_PanDclwz*cos(m_JD1-zypj*3.141592654/2);m_y1=pDlg->m_y1+pDlg->m_PanDclwz*sin(m_JD1-zypj*3.141592654/2); m_y2=pDlg->m_y2+pDlg->m_PanDclwz*sin(m_JD1-zypj*3.141592654/2); m_x3=pDlg->m_x3+pDlg->m_PanDclwz*cos(m_JD2-zypj*3.141592654/2); m_x4=pDlg->m_

22、x4+pDlg->m_PanDclwz*cos(m_JD2-zypj*3.141592654/2); m_y3=pDlg->m_y3+pDlg->m_PanDclwz*sin(m_JD2-zypj*3.141592654/2); m_y4=pDlg->m_y4+pDlg->m_PanDclwz*sin(m_JD2-zypj*3.141592654/2);m_x1qs=atof(pDlg->Res4.GetAt(pDlg->m_lcdh)*4+1)+pDlg->m_PanDclwz*cos(m_JD1-zypj*3.141592654/2); m_

23、y1qs=atof(pDlg->Res4.GetAt(pDlg->m_lcdh)*4+2)+pDlg->m_PanDclwz*sin(m_JD1-zypj*3.141592654/2); JDX=(m_y2-m_y1)*m_x1/(m_x2-m_x1)-(m_y4-m_y3)*m_x3/(m_x4-m_x3)+m_y3-m_y1) /(m_y2-m_y1)/(m_x2-m_x1)-(m_y4-m_y3)/(m_x4-m_x3); JDY=(m_y2-m_y1)*(JDX-m_x1)/(m_x2-m_x1)+m_y1;/计算交点X、Y坐标pDlg->m_sx=m_x1;p

24、Dlg->m_mx=m_x4;pDlg->m_sy=m_y1;pDlg->m_my=m_y4;double Fa1=0.0,Fa2=0.0,Fa3=0.0,Fa4=0.0,Fa5=0.0,Faxx1=0.0,Fayy1=0.0,Faxx2=0.0,Fayy2=0.0;/Fa:方位角的意思CStringArray FaDH;/存点号的CStringArray数组int FaHH300=0;/存行号的数组，赋初值0 intFai=0;CString Biaosf("CLK");FaDH.SetSize(300,-1);/给CStringArray数组分配内存大

25、小，需要足够大的空间 for(i=0;i<300;i+) FaDH.SetAtGrow(i,Biaosf);/将每一个数组值赋值为Biaosf("CLK") for(j=1;j<hsz-1;j+)/j此时是数据个数循环,注意hszFaxx1=atof(pDlg->Res4.GetAt(j*4+1)-atof(pDlg->Res4.GetAt(j-1)*4+1); Fayy1=atof(pDlg->Res4.GetAt(j*4+2)-atof(pDlg->Res4.GetAt(j-1)*4+2);Fa1=QJD(Fayy1,Faxx1);F

26、axx2=atof(pDlg->Res4.GetAt(j+1)*4+1)-atof(pDlg->Res4.GetAt(j*4+1); Fayy2=atof(pDlg->Res4.GetAt(j+1)*4+2)-atof(pDlg->Res4.GetAt(j*4+2);Fa2=QJD(Fayy2,Faxx2);if(fabs(Fa1-Fa2)<3.141592654)/20140502新加,解决第四象限和第一象限夹角的情况 Fa3=fabs(Fa1-Fa2);/m_JD为大曲线偏角,以弧度为单位,取正,为角度拼音的简拼 else if(fabs(Fa1-Fa2)&g

27、t;3.141592654&&Fa1>Fa2) Fa3=fabs(2*3.141592654-Fa1+Fa2);else if(fabs(Fa1-Fa2)>3.141592654&&Fa1<Fa2)Fa3=fabs(-2*3.141592654-Fa1+Fa2);5.2对于曲线图的绘制，采用CadIO的动态链接库。通过OnDxfsave() 函数调用其中的一部分类，结合在LineCal类中已经计算好的曲线数据，进行画图数据的计算，最终自动绘出DXF格式的曲线图。具体代码示例如下：if(m_BZPD2=2)UpdateData(TRUE);/对话

28、框中修改了输入的值，刷新程序中变量的值。LineCal ALC2;ALC2.CalData(); static char szFilter = "DXF file format(*.dxf) |*.dxf|All Files (*.*)|*.*|"CFileDialog SaveDlg( FALSE, NULL, NULL/*LastFilePath*/, OFN_HIDEREADONLY | OFN_EXPLORER, szFilter, NULL );SaveDlg.m_ofn.lpstrTitle = "Save DXF File As"if(Sa

29、veDlg.DoModal()=IDOK)/ Saving sample dxf file dataCString DxfFileName( SaveDlg.GetPathName() );CDrawing drw;if(SaveDlg.GetFileExt().IsEmpty()DxfFileName += ".dxf"/ Tables Section -/ LTYPE table type -LTYPEltype;OBJHANDLEobjhandle1;/ Layer1ZeroMemory(&layer, sizeof(layer);strcpy(layer.N

30、ame, "Layer1");layer.Color = 2;/1红色、2黄、3绿、4青，5蓝，6紫色，7白色layer.LineTypeObjhandle = objhandle1; / Continuousdrw.AddLayer(&layer);/ Layer2ZeroMemory(&layer, sizeof(layer);strcpy(layer.Name, "Layer2");layer.Color = 1;layer.LineTypeObjhandle = objhandle1; / Continuousdrw.AddLay

31、er(&layer);/ STYLE table type -STYLEstyle;ZeroMemory(&style, sizeof(style);strcpy(style.Name, "Style1");strcpy(style.PrimaryFontFilename, "TIMES.TTF");style.LastHeightUsed = 0.3;style.WidthFactor = 1;objhandle1 = drw.AddTextStyle(&style);result &= drw.SetLayer(&qu

32、ot;Layer3");drw.SetTextStyle("Style1");result &= drw.Text(ALC2.ResHuaTu.GetAt(4), atof(ALC2.ResHuaTu.GetAt(3),atof(ALC2.ResHuaTu.GetAt(2), 7, atof(ALC2.ResHuaTu.GetAt(5);/ZHLC备注内容，上面5个参数分别表示：文字内容、起点X坐标（CAD与测量坐标系统相反）、起点Y坐标、文字大小、文字方向（CAD右手定则）result &= drw.Text(ALC2.ResHuaTu.GetA

33、t(10), atof(ALC2.ResHuaTu.GetAt(9), atof(ALC2.ResHuaTu.GetAt(8), 7, atof(ALC2.ResHuaTu.GetAt(57);/HYLC备注内容result &= drw.Text(ALC2.ResHuaTu.GetAt(15), atof(ALC2.ResHuaTu.GetAt(14), atof(ALC2.ResHuaTu.GetAt(13), 7, atof(ALC2.ResHuaTu.GetAt(58);/YHLC备注内容result &= drw.Text(ALC2.ResHuaTu.GetAt(21

34、), atof(ALC2.ResHuaTu.GetAt(20),atof(ALC2.ResHuaTu.GetAt(19), 7, atof(ALC2.ResHuaTu.GetAt(16);/HZLC备注内容result &= drw.Text(m_quxbianh, atof(ALC2.ResHuaTu.GetAt(45), atof(ALC2.ResHuaTu.GetAt(44), 7, atof(ALC2.ResHuaTu.GetAt(24);/曲线编号result &= drw.Text(ALC2.ResHuaTu.GetAt(27), atof(ALC2.ResHuaT

35、u.GetAt(26), atof(ALC2.ResHuaTu.GetAt(25), 7, atof(ALC2.ResHuaTu.GetAt(24);/大偏角result &= drw.Text(ALC2.ResHuaTu.GetAt(30), atof(ALC2.ResHuaTu.GetAt(29), atof(ALC2.ResHuaTu.GetAt(28), 7, atof(ALC2.ResHuaTu.GetAt(24);/半径result &= drw.Text(ALC2.ResHuaTu.GetAt(33), atof(ALC2.ResHuaTu.GetAt(32),

36、atof(ALC2.ResHuaTu.GetAt(31), 7, atof(ALC2.ResHuaTu.GetAt(24);/曲线长result &= drw.Text(ALC2.ResHuaTu.GetAt(36), atof(ALC2.ResHuaTu.GetAt(35), atof(ALC2.ResHuaTu.GetAt(34), 7, atof(ALC2.ResHuaTu.GetAt(24);/切线长result &= drw.Text(ALC2.ResHuaTu.GetAt(39), atof(ALC2.ResHuaTu.GetAt(38), atof(ALC2.Re

37、sHuaTu.GetAt(37), 7, atof(ALC2.ResHuaTu.GetAt(24);/缓和曲线长 atof(ALC2.ResHuaTu.GetAt(22)>0; ENTVERTEX vertex5;ZeroMemory(vertex, sizeof(ENTVERTEX)*5);vertex0.Point.x = 1.2125; vertex0.Point.y = 5.1811;vertex1.Point.x = 2.5888; vertex1.Point.y = 4.9520;vertex2.Point.x = 2.7910; vertex2.Point.y = 7.11

38、88;vertex3.Point.x = 4.7204; vertex3.Point.y = 6.6414;vertex4.Point.x = 5.0015; vertex4.Point.y = 8.0484;result &= drw.PolyLine(vertex, 5);*/result &= drw.SetLayer("Layer1");/设置图层 intk=0;ENTVERTEX vertex100;/画上缓和曲线,容器数据个数的上限if(result>0) MessageBox("保存 DXF Flie 失败.", &q

39、uot;Error", MB_ICONERROR | MB_OK);elseMessageBox("保存 DXF Flie 成功.", "LineCal", MB_ICONINFORMATION | MB_OK);/if纠错elseAfxMessageBox(" 请先计算数据!");6、软件使用流程介绍及现场测试6.1、进行一个平面曲线的具体拟合计算之前，需要先在参数设置里面进行项目信息、轨道型号、测量位置、直线选择相关信息的查看和设置，可以不做修改，采用默认值，也可以按照实际情况修改，其中直线选择这个对话框采用非模态对话框

40、，可以和主窗口并列显示，就方便了直线边的选择和调整。这也是在给济南铁路局既有线测量队的试用时，他们提出的细节改进要求。然后可以选择自动模式还是手动模式。如果选择自动模式计算，想得到拟合的最佳曲线的半径需要精确到米的话，可以在计算完初始半径和大偏角后，在“自动模式参数设置里面”进行进一步的设置，这样可以更有针对性，也会节省一些计算时间（现场演示搜索步长为1m的曲线计算）。6.2、软件在导入数据之后，不论采取手动模式还是自动模式，均只需导入一次数据，中间可以进行不同的参数调整，然后可以多次计算，分别存储，这样可以减少一些操作步骤，以实现程序操作的流程简洁、清楚。（现场演示，自动模式计算时，在直线边

41、选择上的更改）。 具体测试数据见后面附表一7、外业测量流程及注意事项7.1 在每测量一个曲线之前，均需要按照规范要求直线每50m，曲线每20m进行线路丈量，然后按照丈量的里程点进行曲线测量，其中点名可以是任意数值，应该采用全站仪法进行导线测量、里程推算和曲线计算，其他具体测量要求可以参照中铁咨询铁路勘察细则的第二篇测绘P79至P82的相关规定进行。7.2本软件每次导入的数据文件不超过300行，采用自动模式计算大半径曲线数据时（半径>=10000m）或者长缓和曲线（缓和曲线>=300m）数据时，计算时间可能稍微长点。如果曲线上有少量粗差大的点（少于总数的10%），程序在计算的过程中会

42、对其进行自动判断和剔除，不会影响最终计算结果。但有可能会需要将“拨距允许值”和“半径搜索范围”设置大些。在极少部分情况下（当测量的数据质量较差时），采用自动模式无法计算或者计算不出满意数据，可以根据既有线的台账数据，在选取直线边后，采用手动模式，计算获得每一个点的拨距量。8、总结非常感谢李总、王总等院领导对本项目的高度重视和具体指导，同时也感谢贾克勇、高强、刘成伟、李厚霖、李耐振在软件编制过程中的测试和帮助。在给济南局既有线测量队试用过程中，发现软件还存在一些待改进的地方（比如对复心曲线的拟合计算等），需要进一步的完善，如在使用的过程中还有其他建议和要求，请和陈礼坤联系（QQ：28524289

43、5，TEL附表一：以2014年12月4号进行外业测量的蓝翔路下穿立交桥为例下面的数据是测量的线位资料（外轨中心）534063886.386495008.4129.284524063885.416495044.929.259514063884.032495097.42229.2504063882.925495139.06729.188494063882.329495162.23229.193484063881.814495181.98229.185474063881.381495198.21929.167464063881.086495209.24329.1474540

44、63880.631495225.59829.118444063880.121495242.5629.109424063879.179495268.35429.101414063878.25495288.29629.103404063877.038495308.8229.089394063875.721495327.08629.093254063875.605495328.62429.428274063874.529495341.4429.433284063871.949495368.04429.439294063871.762495369.8829.44304063869.794495387.

45、33929.445314063867.471495405.85629.459324063864.708495425.78829.468724063860.891495450.37428.361334063859.316495460.00129.466734063858.091495467.02528.367344063856.169495478.08429.454744063854.707495486.01528.37754063850.127495510.55928.37764063846.243495530.83228.366774063843.664495544.21128.355804

46、063842.839495548.50228.354834063838.52495570.78228.36874063833.522495596.57828.363下表是拟合的曲线数据以及每个点对应的拨距距离。点名里程北坐标(X)东坐标(Y)内拨(mm)外拨(mm)是否超限(100mm)曲线要素15-89.65244063898.855494992.835441.08Aa=9.312901404063896.517495082.45730R:1500.01345.55354063895.329495127.99530L1=70.0L2=70.01266.70634063894.77749514

47、9.140916.65T1=159.9771T2=159.97711686.43544063894.262495168.86531.71l=319.357017106.63334063893.705495189.05554.86ZHLC=79.667318146.4274063892.226495228.821314.45HYLC=149.667319166.49774063891.144495248.860612.13QZLC=239.345820186.52724063889.797495268.844610.01YHLC=329.024321206.72284063888.1684952

48、88.974310.27HZLC=399.024322226.71884063886.289495308.88165.39BoJu=0.317624246.90814063884.123495328.95427.92CurNO: JD1 Pro: CREC26260.50894063882.511495342.45921.34Actor: CLK Check: LDJ61266.69064063881.738495348.59238.3862286.43444063879.1495368.15918.1863306.58414063876.142495388.090312.1564326.70

49、844063872.922495407.9558.4765346.73184063869.462495427.67520.9766366.71494063865.813495447.322126.2667386.99894063861.986495467.241828.2168406.59314063858.243495486.477235.2369426.64954063854.409495506.163921.4970446.71714063850.574495525.86157.8971466.72064063846.75495545.49624.2378470.41354063846.044495549.121082493.0574063841.7