专利名称：基于cqg2000的区域似大地水准面精化方法
技术领域：
本发明涉及地球重力场研究技术领域，尤其是涉及一种基于CQG2000 的区域似大地水准面精化方法。
背景技术：
随着GPS定位技术的广泛应用，人们已经能够在10—7至10—9的精度量 级上，简捷而经济地获得所测点位的平面位置，但是却一直未能以相应的 精度求解点的高程(海拔高)。原因是尽管GPS能给出高精度的大地高， 却由于没有一个具有相应精度和高分辨率的(似)大地水准面模型，至使 在GPS大地高至GPS海拔高的转换中精度严重丢失。为此，目前包括我国
在内的国际大地测量学界都在致力于研究区域性高分辨率高精度似大地 水准面或大地水准面的建立。
(似)大地水准面是获取地理空间信息的高程基准面；高精度(似) 大地水准面是研究海洋动力环境和海洋地球物理、地球动力学等有关地球 科学问题的基础；GPS技术结合高精度高分辨率大地水准面模型，可以取
代传统的水准测量方法测定正高或正常高，真正实现GPS技术在几何和物
理意义上的三维定位功能；在现今GPS定位时代，精化区域性大地水准面 和建立新一代传统的国家或区域性高程控制网同等重要，也是一个国家或 地区建立现代高程基准的主要任务。同时，怎样将已有的(似)大地水准 面成果应用到现实社会生产中，更好地更快地指导成产也是当前测绘学者 在努力研究的主要方面。因此，(似)大地水准面的精化是大地测量学的 一项重要科学目标，也是一个极具实用价值的工程任务。
同时，在确立高精度的似大地水准面精化模型时，需要高精度高分辨
率的重力数据，高分辨率的地形数据以及高精度的GPS水准数据，而实测
4的重力数据不仅测量过程复杂，而且测量费用巨大。实际工作中，有些区 域重力测量工作难以进行，精度和分辨率相对较低，达不到设计要求。地 形数据虽然随着高分辨率的地球监测卫星技术的发展，其分辨率有了相应 的提高，但是往往因为保密等原因，致使其很难应用到工程建设等。因此 我国现有的一些区域似大地水准面是在地面重力测量数据的分辨率偏低 且分布不均匀的情况下建立的，东部地区平均为5' x5',中部地区为5'
x 5'至15' x 15',西部地区为5' x 5'至1° x 1° ,还有一些重力测量 空白区，这就给精化大地水准面工作带了很大的困难，也就使得大地水准 面的成果在很多地区难于应用于实际生产工作中。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一 种基于CQG2000的区域似大地水准面精化方法，其实施简便且测试数据准 确可靠、测试速度快，能有效解决似大地水准面精化时因无法获取重力数 据等所产生的水准面拟合缺陷问题，且能克服利用重力数据精化时所存在 的测量过程复杂、测量费用大的实际问题。
为解决上述技术问题，本发明采用的技术方案是 一种基于CQG2000的 区域似大地水准面精化方法，其特征在于该方法包括以下步骤
步骤一、构建待测区域的区域似大地水准面模型，其主要包括以下步
骤
(a)移去，并计算所需构建区域似大地水准面模型的起算数据值A《 首先，布设待测区域的区域GPS水准网，釆用GPS接收机和水准测量 仪器，对所述待测区域同时进行GPS测量与水准测量，并获取GPS观测数 据和水准测量数据，且相应计算得出GPS水准点的高程异常观测值《；
其次，釆用插值方法，利用CQG2000格网数据插值所述GPS水准点的 高程异常观测值《，得出相应高程异常值&;起伏不大的山区或丘陵地区时， A《=《-^-《T,其中&为利用地形改正公式计算得出的所述GPS水准点对
应的地形改正值；
(b) 拟合以A,作为起算数据，构建数学模型并对所构建的数学模 型进行拟合，即得到A《-AX+S+A,其中，A为由所布设的区域GPS水准 网数据点坐标相应求得的系数矩阵，X为待定参数，S为随机信号，A为 残差；
当A#0, S-O时，所构建的数学模型为函数模型； 当A=0, S邦时，所构建的数学模型为随机模型；
当A#0, S#0时，所构建的数学模型为最小二乘配置模型，釆用最小 二乘模型，估算出X的最或然值f和S的最或然值S;
(c) 恢复
对于待测区域内任一未知点来说，首先釆用常规测量方法对所述未知 点进行测量并获取对应的空间坐标；同时釆用插值方法，利用CQG2000格 网数据插值计算得出所述未知点的高程异常观测值^;再利用步骤(b) 中所构建得出的数学模型A《-AX+S+A，计算得出所述未知点对应的模型值 Ld;之后，再计算得出所述未知点的高程异常最终值么，当所述待测区域 为地形起伏不大的平原地区时，6=Ld+4;当所述待测区域为地形起伏不 大的山区或丘陵地区时，《=Ld+" + ^，其中g为利用地形改正公式计算
得出所述未知点对应的地形改正值；
同时，将所推算得出的高程异常最终值《的计算公式输入至运算控制
器；
步骤二、采用GPS系统获取待测区域内任一未知点的实时大地坐标，
并将所测得数据同步传送至所述运算控制器进行处理运算，所述运算控制 器利用步骤一中推算得出的高程异常最终值^的计算公式，即可对应计算
得出所测未知点的实时正常高即水准高。
步骤(a)中所述的插值方法为Shepard插值法。
6步骤二中所述的GPS系统为GPS接收机，步骤一和步骤二中所述的运 算控制器为与GPS接收机相接的PDA。
步骤(a)中所述的地形改正公式为利用FFT/FTH方法计算地形改正 球面坐标形式的严格卷积公式。
本发明与现有技术相比具有以下优点，1、本发明利用我国最先进的 似大地水准面模型CQG2000成果，避开了重力实际测量的复杂性，同时也 减少工程预算费用，通过将GPS/水准数据与CQG2000数据进行结合，从而 使我国最新的最先进的似大地水准面模型CQG2000的成果数据直接应用到 实践生产中去，直接指导现实的社会生产活动；我国新一代似大地水准面 CQG200G采用了国内外最新的重力场资料、地形资料和卫星测高资料，其 精度总体达到了分米级水平(±0. 3-0. 6m)，在中东部地区其精度优于 0. 3m。其陆地区域的分辨率为5' X5'，海域分辨率为15' Xl5' 。 2、 利用本发明所述的方精化法，建立区域(似)大地水准面模型，并且在此 模型的平台上，结合高精度GPS大地高就可以快速地高精度地获取地面点 的正常高，将极大改善传统高程测量作业模式，满足长距离、长跨度工程 水准测量的精度要求，从而使费用高、难度大、周期长的传统水准作业减 少到最低限度；结合GPS测量可以满足大比例尺测图对精度的要求，为加 快各种基础测绘、国土资源、水利、交通、生态环境保护等工程的建设， 具有重要的科学意义、社会效益和巨大的经济效益。3、以本发明所建立 (似)大地水准面模型，代表着一个高精度的测量基准，其可以满足各种 不同行业用户对精密定位，快速和实时定位、导航的要求，及时地满足城 巿规划、国土测绘、地籍管理、城乡建设、环境监测、防灾减灾、交通监 控等多种现代信息化管理的社会需求，这正是现有的同等精度级别的似大 地水准面精化软件所不具备的。综上，本发明以CQG2000为平台，按照移 去拟合恢复法来确立待测区域的似大地水准面模型，具体而言目前计算 区域重力似大地水准面广泛采用的移去恢复法，主要是利用重力场的可叠 加性原理，分别处理不同波长成分的贡献，再经过简单的叠加恢复所逼近的局部重力场；本发明将现有移去恢复法原理用在基于CQG2000的大地水 准面精化上，以体现了局部(似)大地水准面物理性质的CQG2000作为平 台，求解GPS水准点上高程异常观测值与CQG2000在该点的计算值之残差，
并构建数学模型进行拟合，然后利用该模型计算待定点的同类残差，再与 CQG2000在该点的计算值叠加，最终获得待定点高程异常最终值。综上所 述，本发明利用所构建的被测区域的区域似大地水准面模型，能实时将所 测未知点的GPS大地高转换为正常高即水准高，其实施简便且测试数据准 确可靠、测试速度快，能有效解决似大地水准面精化时因无法获取重力数 据等所产生的水准面拟合缺陷问题，以及利用重力数据精化时所存在的测 量过程复杂、测量费用大的问题，大大了提高似大地水准面的精度和分辨 率，使用价值非常高。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
图1为本发明第一实施方式的流程框图。 图2为本发明第二实施方式的流程框图。
具体实施例方式
实施例1
如图1所示，本发明所述的基于CQG2000的区域似大地水准面精化方
法，包括以下步骤
步骤一、构建待测区域的区域似大地水准面模型，其主要包括以下步
骤
(a)移去，并计算所需构建区域似大地水准面模型的起算数据值A《 首先，布设待测区域的区域GPS水准网，釆用GPS接收机和水准测量 仪器，对所述待测区域同时进行GPS测量与水准测量，并获取GPS观测数 据和水准测量数据，且相应计算得出GPS水准点的高程异常观测值《。本
8实施例中，所述待测区域为地形起伏不大的平原地区。
也就是说，在精化前期，按照《国家一、二等水准测量规范》(国家
技术监督局GB12897-1991 )、《全球定位系统(GPS )测量规范》(GB/T 18314-2001 )及《测绘技术设计规定》(ZBA 75001-89)等相关技术规范 布设区域GPS水准网，同时采用所述GPS系统和水准测量仪器进行GPS测 量与高等级水准测量，以便获取GPS观测数据和高等级水准测量数据，其
网形及强度、密度和最弱点精度都需要按照一定的设计标准和规范执行， 并用高精度的GPS解算软件GAMIT/GLOBK 10. 21软件，选取2005国家GPS 大地控制网的框架和历元，进行数据处理，并充分利用国内外的GPS连续 运行站作为参考，在处理数据时以提高整网的精度。GPS网基线解算时以 同步环为单位，同步观测数据同步处理。网平差釆用整体平差，在WGS-84 椭球上采用长安大学空间定位技术和信息研究所自主研制HPGPSADJ高精 度GPS网平差软件进行整体平差处理，进而得出GPS水准点的相关信 息一一即大地坐标和正常高。
其次，釆用插值方法，利用CQG2000格网数据插值所述GPS水准点的 高程异常观测值《，得出相应高程异常值"即所述GPS水准点的高程异常值&。
本实施例中，所述插值方法为spline插值法，所述spline插值法是 用三次多项式，对相邻数据点之间的各段建模，每个三次多项式的头两个 导数在该相邻数据点上相一致。spline插值法是以计算点为中心，取拟合 半径R以内已知函数值的权中数，数据点上的权按距离计算点的不同范围 釆用不同的权函数式，使靠近中心点的权增大，远离中心的权迅速减少。 在Shepard局部内插模型中，选R-0.25。，并规定
i 内插的函数模型为
<formula>formula see original document page 10</formula>当^ # 0时
当~ =0时
(2)
式中
。=
利用(2)式便可计算得出待定点的高 程异常值，其中(x， y)为待定点坐标，(I, yJ为已知点坐标。
接着，计算起算数据A《由于所述待测区域为地形起伏不大的平原地 区，则
(b) 拟合以^作为起算数据，构建数学模型并对所构建的数学模 型进行拟合，即得到A《=AX+S+A,其中，A为由所布设的区域GPS水准 网数据点坐标相应求得的系数矩阵，X为待定参数，S为随机信号，A为 残差；
当A#0, S-O时，所构建的数学模型为函数模型； 当A^， S邦时，所构建的数学模型为随机模型；
当A#0， S#0时，所构建的数学模型为最小二乘配置模型，釆用最小 二乘模型，估算出X的最或然值J和S的最或然值S;
(c) 恢复
对于待测区域内任一未知点来说，首先釆用常规测量方法对所述未知 点进行测量并获取并获取对应的空间坐标；同时采用插值方法，利用 CQG2000格网数据插值计算得出所述未知点的高程异常观测值^;再利用
步骤(b)中所构建得出的数学模型A《-AX+S+A,计算得出所述未知点对 应的模型值Ld;之后，再计算得出所述未知点的高程异常最终值么，由于 所述待测区域为地形起伏不大的平原地区，则^ = L d +《w 。
步骤二、釆用GPS系统获取待测区域内任一未知点的实时大地坐标，
并将所测得数据同步传送至所述运算控制器进行处理运算，所述运算控制 器利用步骤一中推算得出的高程异常最终值^的计算公式，即可对应计算得出所测未知点的实时正常高即水准高。
本实施例中，步骤二中所述的GPS系统为GPS接收机，步骤一和步骤
二中所述的运算控制器为与GPS接收机相接的PDA (personal digital assistant )。 实施例2
如图2所示，所述待测区域为地形起伏不大的山区或丘陵地区，本发 明所述的基于CQG2000的区域似大地水准面精化方法，其步骤与实施例1
大致相同，主要区别在于需要计算地形改正模型，其主要包括以下步骤 U)移去，并计算所需构建区域似大地水准面模型的起算数据值A^
首先，布设待测区域的区域GPS水准网，釆用GPS接收机和水准测量 仪器，对所述待测区域同时进行GPS测量与水准测量，并获取GPS观测数 据和水准测量数据，且相应计算得出GPS水准点的高程异常观测值《。
本实施例中，设有m个GPS水准重合点且其对应的高程异常观测值分 别为^ (k=l...m)。
其次，釆用插值方法，利用CQG2000格网数据插值所述GPS水准点的 高程异常观测值《，得出相应高程异常值^即所述GPS水准点的高程异常值"。
本实施例中，在上述m个GPS水准重合点上，釆用Shepard插值法， 利用CQG2000格网数据插值m个GPS水准重合点的高程异常观测值^，得 出相应高程异常值《。
接着，计算起算数据A^由于所述待测区域为地形起伏不大的山区或 丘陵地区，则^-《-^-^，其中^为利用地形改正公式计算得出的所述
GPS水准点对应的地形改正值。
本实施例中，利用地形改正公式计算上述m个水准重合点上的地形改 正值"，并计算厶^ = ^-
同时，所述地形改正公式为利用FFT/FTH方法计算地形改正球面坐标 形式的严格卷积公式。在计算地形改正时，用现有高分辨率DEM (Digital Elevation Model 数字高程模型)数据，同时采用先进的FFT/FHT技术计算地形改正。由于 DEM数据是以球坐标形式的格网高程数据，因而为了减少和消除过去计算 地形改正釆用平面坐标形式的二维卷积公式引起的近似误差，选取利用 FFT/FTH方法计算地形改正的球面坐标形式的严格卷积公式。
地形改正可用下面的公式计算
z厂 ln(少+ +少ln(x + " + ——}x2少2z2
X乂5
其中，取2,=0且22=^ (即高差)时，公式(3)所求得的函数值即
为地形改正值电.(.；取z,--(仏+r), (//。+r+4.45//),且用补偿密度A^ 代替P,公式(3)所求得的函数值即为均衡改正值糸/5,其中T为均衡深
度，仏为重力点高程。
另外，地形改正的另 一 种方法就是谱方法确定网格地形改正。
<%『=^7Tl+^rC2 (4)
=} [//2^ ] — 2 Vf间^ ] +《巧-1 (5)
电(2=-糊K-化A]
式中，巧—'表示二维Fourier变换逆算子，仏=《[^], k=0， 1， 2， 3， 4; 《为二维Fourier变换逆算子，k=l， 2, /为两点间平面距离，
h为高程。
(b)拟合以^作为起算数据，构建数学模型并对所构建的数学模 型进行拟合，即得到A《=AX+S+A，其中，A为由所布设的区域GPS水准 网数据点坐标相应求得的系数矩阵，X为待定参数，S为随机信号，A为 残差；
当A#0， S-0时，所构建的数学模型为函数模型； 当A=0， S^)时，所构建的数学模型为随机模型； 当A#0, S邦时，所构建的数学模型为最小二乘配置模型，釆用最小
12的最或然值f和S的最或然值f。 本实施例中，所述起算数据A《为A^。 (C)恢复
对于待测区域内任 一 未知点来说，首先采用常规测量方法对所述未知
点进行测量并获取对应的空间坐标，同时釆用插值方法，利用CQG2000格 网数据插值计算得出所述未知点的高程异常观测值"；再利用步骤(b) 中所构建得出的数学模型A《-AX+S+A，计算得出所述未知点对应的模型值 Ld;之后，再计算得出所述未知点的高程异常最终值《，由于所述待测区 域为地形起伏不大的平原地区，则《,Ld+" + g,其中^为利用地形改正
公式计算得出所述未知点对应的地形改正值。同时，将所推算得出的高程 异常最终值^的计算公式输入至运算控制器。
在本步骤中，所用的地形改正公式为与步骤(a)中相对应的地形改 正公式。
步骤二、采用GPS系统获取待测区域内任一来知点的实时大地坐标，
并将所测得数据同步传送至所述运算控制器进行处理运算，所述运算控制 器利用步骤一中推算得出的高程异常最终值二的计算公式，即可对应计算
得出所测未知点的实时正常高即水准高。
以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是 根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构 变化，均仍属于本发明技术方案的保护范围内。
1权利要求
1.一种基于CQG2000的区域似大地水准面精化方法，其特征在于该方法包括以下步骤步骤一、构建待测区域的区域似大地水准面模型，其主要包括以下步骤(a)移去，并计算所需构建区域似大地水准面模型的起算数据值Δξ首先，布设待测区域的区域GPS水准网，采用GP8接收机和水准测量仪器，对所述待测区域同时进行GPS测量与水准测量，并获取GPS观测数据和水准测量数据，且相应计算得出GPS水准点的高程异常观测值ξ；其次，采用插值方法，利用CQG2000格网数据插值所述GPS水准点的高程异常观测值ξ，得出相应高程异常值ξg；接着，计算起算数据Δξ当所述待测区域为地形起伏不大的平原地区时，Δξ＝ξ-ξg；当所述待测区域为地形起伏不大的山区或丘陵地区时，Δξ＝ξ-ξg-ξT，其中ξT为利用地形改正公式计算得出的所述GPS水准点对应的地形改正值；(b)拟合以Δξ作为起算数据，构建数学模型并对所构建的数学模型进行拟合，即得到Δξ＝AX+S+Δ，其中，A为由所布设的区域GPS水准网数据点坐标相应求得的系数矩阵，X为待定参数，S为随机信号，Δ为残差；当A≠0，S＝0时，所构建的数学模型为函数模型；当A＝0，S≠0时，所构建的数学模型为随机模型；当A≠0，S≠0时，所构建的数学模型为最小二乘配置模型，采用最小二乘模型，估算出X的最或然值<overscore>X</overscore>和S的最或然值<overscore>S</overscore>；(c)恢复对于待测区域内任一未知点来说，首先采用常规测量方法对所述未知点进行测量并获取对应的空间坐标；同时采用插值方法，利用CQG2000格网数据插值计算得出所述未知点的高程异常观测值ξw；再利用步骤(b)中所构建得出的数学模型Δξ＝AX+S+Δ，计算得出所述未知点对应的模型值Ld；之后，再计算得出所述未知点的高程异常最终值ξe，当所述待测区域为地形起伏不大的平原地区时，ξe＝Ld+ξw；当所述待测区域为地形起伏不大的山区或丘陵地区时，<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><msub> <mi>ξ</mi> <mi>e</mi></msub><mo>=</mo><mi>Ld</mi><mo>+</mo><msub> <mi>ξ</mi> <mi>w</mi></msub><mo>+</mo><msubsup> <mi>ξ</mi> <mi>T</mi> <mi>W</mi></msubsup><mo>,</mo> </mrow>]]></math> id="icf0001" file="A2009100214690003C1.tif" wi="29" he="4" top= "58" left = "79" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>其中ξTW为利用地形改正公式计算得出所述未知点对应的地形改正值；同时，将所推算得出的高程异常最终值ξe的计算公式输入至运算控制器；步骤二、采用GPS系统获取待测区域内任一未知点的实时大地坐标，并将所测得数据同步传送至所述运算控制器进行处理运算，所述运算控制器利用步骤一中推算得出的高程异常最终值ξe的计算公式，即可对应计算得出所测未知点的实时正常高即水准高。
2. 按照权利要求1所述的基于CQG2000的区域似大地水准面精化方 法，其特征在于步骤(a)中所述的插值方法为Shepard插值法。
3. 按照权利要求1或2所述的基于CQG2000的区域似大地水准面精化 方法，其特征在于步骤二中所述的GPS系统为GPS接收机，步骤一和步 骤二中所述的运算控制器为与GPS接收机相接的PDA。
4. 按照权利要求1或2所述的基于CQG2000的区域似大地水准面精化 方法，其特征在于步骤(a)中所述的地形改正公式为利用FFT/FTH方 法计算地形改正球面坐标形式的严格巻积公式。
全文摘要
本发明公开了一种基于CQG2000的区域似大地水准面精化方法，包括步骤一、以CQG2000为平台，利用GPS水准数据按照移去拟合恢复法确立待测区域的似大地水准面模型；步骤二、用GPS系统获取待测区域内任一未知点的实时大地坐标，并将测得数据传送至运算控制器，运算控制器利用所构建的区域似大地水准面模型，实时将所测未知点的GPS大地高转换为正常高即水准高。本发明实施简便且测试数据准确可靠、测试速度快，能有效解决似大地水准面精化时因无法获取重力数据等所产生的水准面拟合缺陷问题，以及利用重力数据精化时所存在的测量过程复杂、测量费用大的问题，大大了提高似大地水准面的精度和分辨率，使用价值非常高。
文档编号G01C5/00GK101493324SQ200910021469
公开日2009年7月29日 申请日期2009年3月10日 优先权日2009年3月10日
发明者刘站科, 勤 张, 利 王, 程传录, 黄观文 申请人:长安大学