一种基础资料缺乏地区测量控制点成果获取方法
【专利摘要】本发明公开了一种基础资料缺乏地区测量控制点成果获取方法，包括如下步骤：步骤一、在测绘区域选定测量控制点桩位位置，埋设控制点标志，采用GPS设备按照静态模式测量卫星数据，并根据测量数据计算出控制点的WGS84坐标；步骤二、获取并选定测绘区域周边的转换区域中兼具WGS84坐标与国家坐标的标注点作为重合点；步骤三、求解重合点WGS84坐标与国家坐标的换算关系；步骤四、根据控制点的WGS84坐标和前述的换算关系计算出控制点的国家坐标。本发明具有能够方便快捷地获取控制点成果、成果精度高、满足工程项目前期测绘要求的特点，可以广泛应用于测绘【技术领域】。
【专利说明】一种基础资料缺乏地区测量控制点成果获取方法
【技术领域】
[0001] 本发明涉及测绘【技术领域】，特别是涉及一种基础资料缺乏地区测量控制点成果获 取方法。
【背景技术】
[0002] 目前，工程项目测绘地形图需要用到国家测绘控制点坐标成果。而在基础资料缺 乏地区进行项目前期测绘，如何获取满足项目前期测绘要求的控制点成果，是困扰电力测 绘人员的一大技术难题。在大山区或偏远地区进行电力工程测量时，控制点相对稀少甚至 没有测量控制点成果，即使有控制点也需要远距离联测，既费时又费力。同时，某些地区即 使有控制点成果，由于测量控制点坐标属于涉密成果，需要测绘行政主管部门审批，审批周 期较长，难以满足项目前期阶段勘测设计市场化快速运作要求。
【发明内容】
[0003] 本发明的目的是为了克服上述【背景技术】的不足，提供一种基础资料缺乏地区测量 控制点成果获取方法，具有能够方便快捷地获取控制点成果、成果精度高、满足工程项目前 期测绘要求的特点。
[0004] 本发明提供的一种基础资料缺乏地区测量控制点成果获取方法，包括如下步 骤：步骤一、在测绘区域选定测量控制点桩位位置，埋设控制点标志，采用GPS(Global Positioning System，全球定位系统）设备按照静态模式测量卫星数据，并根据测量数据计 算出控制点的WGS84(World Geodetic Systeml984,是为GPS全球定位系统使用而建立的坐 标系统）坐标；步骤二、获取并选定测绘区域周边的转换区域中兼具WGS84坐标与国家坐标 的标注点作为重合点；步骤三、求解重合点WGS84坐标与国家坐标的换算关系；步骤四、根 据控制点的WGS84坐标和前述的换算关系计算出控制点的国家坐标。
[0005] 在上述技术方案中，所述步骤二中，获取并选定重合点的过程如下：A、依据外业技 术总结与点之记相结合的方法选取均匀分布于整个转换区域的标注点；B、通过单点定位计 算方法计算出所述标注点的WGS84坐标，并利用EGM2008 (是近年来由美国国家地理空间情 报局释放的全球超高阶地球重力场模型）重力场模型通过点位内插计算方法计算标注点 的高程异常，将GPS测量的大地高程减去高程异常即得到该点的WGS84高程；C、并对算好的 标注点进行筛选和试算，剔除粗差点，获得重合点。
[0006] 在上述技术方案中，所述步骤三中，重合点WGS84坐标与国家坐标的换算采用二 维坐标转换模式。
[0007] 在上述技术方案中，所述步骤三中，所述二维坐标转换模式为计算重合点的二维 坐标换算参数，所述二维坐标换算参数包括平移参数、旋转参数和尺度参数。
[0008] 在上述技术方案中，所述步骤四中，控制点国家坐标的计算步骤如下：A、根据二维 坐标换算参数计算出控制点的国家平面坐标；B、分别计算各重合点WGS84高程与国家坐标 高程之间的高差，再根据各重合点WGS84平面坐标、前述高差以及控制点的WGS84平面坐标 计算控制点处WGS84高程与国家坐标高程之间的高差，再将控制点WGS84高程与控制点处 的1?差相减得到控制点的国家坐标1?程。
[0009] 在上述技术方案中，所述步骤三中，重合点WGS84坐标与国家坐标的换算采用三 维坐标转换模式。
[0010] 在上述技术方案中，所述步骤三中，所述三维坐标转换模式为计算重合点的三维 坐标换算参数，所述三维坐标换算参数包括平移参数、旋转参数和尺度参数。
[0011] 在上述技术方案中，所述步骤四中，控制点国家坐标的计算步骤如下：A、根据三维 坐标换算参数计算出控制点的三维国家坐标，再将控制点的三维国家坐标换算成控制点的 国家平面坐标；B、分别计算各重合点WGS84高程与国家坐标高程之间的高差，再根据各重 合点WGS84平面坐标、前述高差以及控制点的WGS84平面坐标计算控制点处WGS84高程与 国家坐标高程之间的高差，再将控制点WGS84高程与控制点处的高差相减得到控制点的国 家坐标高程。
[0012] 在上述技术方案中，所述步骤一中，利用测量数据计算出控制点的WGS84坐标的 过程为：将测量数据转换成Rinex(Receiver Independent Exchange Format/与接收机 无关的交换格式）格式文件，将测量数据结合下载的精密星历和精密种差文件进行单点计 算，算出WGS84坐标系下控制点的经纬度、大地高、正常高和三维坐标。
[0013] 在上述技术方案中，所述国家坐标包括1954北京坐标系和1956黄海高程。
[0014] 本发明基础资料缺乏地区测量控制点成果获取方法，具有以下有益效果：能够方 便快捷地获取控制点成果，成果精度高，提高了工程测绘企业的市场反应能力，完全满足工 程项目前期测绘要求。
【专利附图】
【附图说明】
[0015] 图1为本发明基础资料缺乏地区测量控制点成果获取方法的流程示意图；
[0016] 图2为本发明基础资料缺乏地区测量控制点成果获取方法步骤二中重合点数据 库的功能结构；
[0017] 图3为本发明基础资料缺乏地区测量控制点成果获取方法步骤四中GPS大地高与 正常1?的关系图；
[0018] 图4为本发明基础资料缺乏地区测量控制点成果获取方法所选实施例中重合点 与控制点的空间几何关系示意图。
【具体实施方式】
[0019] 下面结合附图及实施例对本发明作进一步的详细描述，但该实施例不应理解为对 本发明的限制。
[0020] 参见图1，本发明基础资料缺乏地区测量控制点成果获取方法，包括如下步骤：
[0021] 步骤一、在测绘区域选定测量控制点桩位位置，埋设控制点标志，采用GPS设备按 照静态模式测量卫星数据，并根据测量数据计算出控制点的WGS84坐标：
[0022] 在测绘区域埋设控制点标志是控制点成果获取的基础，GPS设备按照静态模式测 量控制点的卫星数据后，利用GPS设备随机软件将测量数据转换成Rinex格式文件（后缀 为 〇 文件和 η 文件），如 algol840. 08η 和 algol840. 08〇。
[0023] 利用单点定位软件下载精密星历和精密种差（时间与测量卫星数据记录时间相 同），解压精密星历和精密种差文件（如codl4862. elk和codl4862.印h)。将测量数据结 合下载的精密星历和精密种差文件进行单点计算，算出WGS84坐标系下控制点的经纬度、 大地高、正常高和三维坐标。
[0024] 步骤二、获取并选定测绘区域周边的转换区域中兼具WGS84坐标与国家坐标的标 注点作为重合点：
[0025] A、依据外业技术总结与点之记相结合的方法选取均匀分布于整个转换区域的标 注点：
[0026] 进行两种不同坐标系的转换计算，需要具有两种坐标重合的点。在本发明中，两种 坐标系指WGS84坐标系与国家坐标系，其中国家坐标系包括平面坐标和高程，我国的平面 坐标为1954北京坐标系（简称北京54)或1980西安坐标系（简称西安80)，高程为1956 黄海高程（简称黄海高程）。获取重合点一方面是通过实测获取，另一方面是通过收集获 取。
[0027] 两种不同坐标系坐标转换的精度除取决于坐标转换的数学模型与求解转换参数 的重合点坐标精度外，还与重合点的多少和几何形状结构有关。
[0028] 重合点选取原则是：依据外业技术总结与点之记等方法选取足够的高等级、高精 度且分布均匀的点作为坐标转换的重合点。采用二维坐标转换模式至少选取2个以上的重 合点，采用三维坐标转换模式至少选取3个以上的重合点，重合点的分布要覆盖整个转换 区域且尽量分布均匀。
[0029] B、通过单点定位计算方法计算出所述标注点的WGS84坐标，并利用EGM2008重力 场模型通过点位内插计算方法计算标注点的高程异常，将GPS测量的大地高程减去高程异 常即得到该点的WGS84高程：
[0030] 其中,EGM2008重力场模型先通过点位内插计算方法计算GPS标注点重力异常，再 由重力异常计算出标注点高程异常，然后求出该点的WGS84坐标高程，具体参见步骤四中B 条b款（一）项。
[0031] C、并对算好的标注点进行筛选和试算，剔除粗差点，获得重合点：
[0032] 应对参与求解转换参数的重合点进行认真分析、筛选、试算，剔除粗差点（即误差 较大的点），采用不含粗差、分布均匀且能包围转换区域的一定密度的重合点求解转换参 数。
[0033] 转换参数的精度与重合点分布的密度、均匀性有关，对于普通工程勘测设计单位 来说，不可能在全国范围内布设一定密度的重合点，只能根据原有的工程项目资料情况，选 择一定数量的重合点。
[0034] 可以利用数据库软件SQL ServerfOOS建立本单位测量控制点数据库，该数据库实 现以下4个功能模块：系统管理模块、基准参数管理模块、电子地图数据管理模块、控制点 数据管理模块，参见图2。
[0035] 步骤三、求解重合点WGS84坐标与国家坐标的换算关系：
[0036] 实现不同坐标系的转换包括不同空间直角坐标系的转换和不同大地坐标系的转 换。不同空间直角坐标系的转换既包括不同参心空间直角坐标系的转换，也包括参心空间 直角坐标系和地心空间直角坐标系的转换。不同大地坐标系的转换既包括不同参心大地坐 标系的转换，也包括参心大地坐标系和地心大地坐标系的转换。
[0037] 坐标转换通常采用分区转换法，即将整个转换区域划分成若干个分区，分别对各 分区计算坐标换算参数。
[0038] 不同坐标系之间的坐标转换通常有两类坐标转换模式：一类是二维坐标转换模 式；一类是三维坐标转换模式。对于不同坐标系之间的坐标换算，目前常用的是平面四参数 转换模型与Bursa七参数转换模型两种。其中，平面四参数转换模型适于二维坐标转换模 式，Bursa七参数转换模型适于三维坐标转换模式。
[0039] A、重合点WGS84坐标与国家坐标的换算采用二维坐标转换模式，即计算重合点的 二维坐标换算参数，所述二维坐标换算参数包括平移参数（dx、dy)、旋转参数Θ和尺度参 数(1λ共计四个参数：
[0040] 二维坐标转换模式只适合于小区域转换，且只需要两坐标系的二维坐标成果（高 斯平面直角坐标X、y或大地经度L、大地纬度Β)。
[0041] 相对而言，平面四参数转换模型原理简单，数值稳定可靠。对较小区域，它转换的 精度较高。但当范围较大时，由于受投影变形误差的影响，其转换精度就较差，因而它只适 合于较小区域的坐标转换。
[0042] B、重合点WGS84坐标与国家坐标的换算采用三维坐标转换模式，即计算重合点 的三维坐标换算参数，所述三维坐标换算参数包括平移参数（dx、dy、dz)、旋转参数（Θ X、 Θ y、θ z)和尺度参数d λ共计七个参数：
[0043] 三维坐标转换模式适合任何区域坐标转换，且需要两坐标系的三维坐标成果（空 间直角坐标X、Υ、Ζ或大地纬度Β、大地经度L、大地高Η)。
[0044] Bursa七参数转换模型为三维模型，在空间直角坐标系中，两坐标系之间存在严密 的转换模型。由于理论比较严密，不存在模型误差和投影变形误差，因而它适合于任何区域 的坐标转换。
[0045] 步骤四、根据控制点的WGS84坐标和前述的换算关系计算出控制点的国家坐标：
[0046] A、二维坐标转换模式下控制点国家坐标的计算步骤如下：
[0047] a、根据二维坐标换算参数计算出控制点的国家平面坐标；
[0048] b、分别计算各重合点WGS84高程与国家坐标高程之间的高差，再根据各重合点 WGS84平面坐标、前述高差以及控制点的WGS84平面坐标计算控制点处WGS84高程与国家坐 标高程之间的高差，再将控制点WGS84高程与控制点处的高差相减得到控制点的国家坐标 高程。
[0049] B、三维坐标转换模式下控制点国家坐标的计算步骤如下（下面结合具体实施例 对三维坐标转换模式的计算步骤加以详述）：
[0050] a、根据三维坐标换算参数计算出控制点的三维国家坐标，再将控制点的三维国家 坐标换算成控制点的国家平面坐标：
[0051] 本实施例以襄樊电厂（103)、岳阳电厂（U087)、浠水核电（D3)、广水接地极（D01)4 个项目控制点计算换算参数，利用武汉特高压控制点进行校核比较。
[0052] 原则上，在计算各分区坐标换算参数时，为了保持各分区在接边处坐标换算参数 的连续性，需要各分区之间相互重叠一部分重合点并重复使用以求取坐标换算参数。当然， 为了计算简便，本实施例中每块区域只选择了一个重合点。
[0053] 区域范围约为200千米X200千米。计算襄樊电厂、岳阳电厂、浠水核电、广水接 地极4个项目点三维坐标换算参数如下：
[0054] 平移参数：
[0055] dxO:***· 97138
[0056] dyO:***· 49047
[0057] dzO:**· 74987
[0058] 旋转参数：
[0059] θ χ：〇. 34878190100
[0060] Θ y：〇. 08062259900
[0061] θ z:2. 07731135100
[0062] 尺度参数
[0063] d λ :-〇· 00001262600
[0064] 武汉特高压变电站WGS84三维坐标如下表1 :
[0065] 表1武汉特高压变电站WGS84三维坐标
[0066]
【权利要求】
1. 一种基础资料缺乏地区测量控制点成果获取方法，其特征在于：包括如下步骤： 步骤一、在测绘区域选定测量控制点桩位位置，埋设控制点标志，采用GPS设备按照静 态模式测量卫星数据，并根据测量数据计算出控制点的WGS84坐标； 步骤二、获取并选定测绘区域周边的转换区域中兼具WGS84坐标与国家坐标的标注点 作为重合点； 步骤三、求解重合点WGS84坐标与国家坐标的换算关系； 步骤四、根据控制点的WGS84坐标和前述的换算关系计算出控制点的国家坐标。
2. 根据权利要求1所述的基础资料缺乏地区测量控制点成果获取方法，其特征在于： 所述步骤二中，获取并选定重合点的过程如下： A、 依据外业技术总结与点之记相结合的方法选取均匀分布于整个转换区域的标注 占 . B、 通过单点定位计算方法计算出所述标注点的WGS84坐标，并利用EGM2008重力场模 型通过点位内插计算方法计算标注点的高程异常，将GPS测量的大地高程减去高程异常即 得到该点的WGS84高程； C、 并对算好的标注点进行筛选和试算，剔除粗差点，获得重合点。
3. 根据权利要求2所述的基础资料缺乏地区测量控制点成果获取方法，其特征在于： 所述步骤三中，重合点WGS84坐标与国家坐标的换算采用二维坐标转换模式。
4. 根据权利要求3所述的基础资料缺乏地区测量控制点成果获取方法，其特征在于： 所述步骤三中，所述二维坐标转换模式为计算重合点的二维坐标换算参数，所述二维坐标 换算参数包括平移参数、旋转参数和尺度参数。
5. 根据权利要求4所述的基础资料缺乏地区测量控制点成果获取方法，其特征在于： 所述步骤四中，控制点国家坐标的计算步骤如下： A、 根据二维坐标换算参数计算出控制点的国家平面坐标； B、 分别计算各重合点WGS84高程与国家坐标高程之间的高差，再根据各重合点WGS84 平面坐标、前述高差以及控制点的WGS84平面坐标计算控制点处WGS84高程与国家坐标高 程之间的高差，再将控制点WGS84高程与控制点处的高差相减得到控制点的国家坐标高 程。
6. 根据权利要求2所述的基础资料缺乏地区测量控制点成果获取方法，其特征在于： 所述步骤三中，重合点WGS84坐标与国家坐标的换算采用三维坐标转换模式。
7. 根据权利要求6所述的基础资料缺乏地区测量控制点成果获取方法，其特征在于： 所述步骤三中，所述三维坐标转换模式为计算重合点的三维坐标换算参数，所述三维坐标 换算参数包括平移参数、旋转参数和尺度参数。
8. 根据权利要求7所述的基础资料缺乏地区测量控制点成果获取方法，其特征在于： 所述步骤四中，控制点国家坐标的计算步骤如下： A、 根据三维坐标换算参数计算出控制点的三维国家坐标，再将控制点的三维国家坐标 换算成控制点的国家平面坐标； B、 分别计算各重合点WGS84高程与国家坐标高程之间的高差，再根据各重合点WGS84 平面坐标、前述高差以及控制点的WGS84平面坐标计算控制点处WGS84高程与国家坐标高 程之间的高差，再将控制点WGS84高程与控制点处的高差相减得到控制点的国家坐标高 程。
9. 根据权利要求1至8中任一项所述的基础资料缺乏地区测量控制点成果获取方法， 其特征在于：所述步骤一中，利用测量数据计算出控制点的WGS84坐标的过程为：将测量 数据转换成Rinex格式文件，将测量数据结合下载的精密星历和精密种差文件进行单点计 算，算出WGS84坐标系下控制点的经纬度、大地高、正常高和三维坐标。
10. 根据权利要求1至8中任一项所述的基础资料缺乏地区测量控制点成果获取方法， 其特征在于：所述国家坐标包括1954北京坐标系和1956黄海高程。
【文档编号】G01C7/02GK104111061SQ201410311051
【公开日】2014年10月22日 申请日期:2014年6月30日 优先权日:2014年6月30日
【发明者】胡吉伦, 程正逢, 曾渠丰, 周毅, 王圣祖, 徐辉, 张建, 刘佳莹, 陈功 申请人:中国电力工程顾问集团中南电力设计院