专利名称：基于无缝分区技术的超大区域似大地水准面精化方法
技术领域：
本发明涉及测绘技术领域，尤其是一种基于无缝分区技术的超大区域似大地水准面精化方法。
背景技术：
似大地水准面精化对于测绘工作有重要意义(I)似大地水准面是获取地理空间信息的高程基准面；(2) GPS (全球定位系统)技术结合高精度、高分辨率的似大地水准面模型，可以取代传统的水淮测量方法测定正常高，真正实现GPS技术对几何和物理意义上的三维定位功能；(3)在现今GPS定位时代，精化区域性似大地水准面与建立国家高精度GPS 定位网或传统的国家高程控制网同等重要，是一个国家发展测绘事业的一项大地测量基础建设，也是一个国家或地区建立现代高程基准的主要任务，以此满足国家经济建设和测绘科学技术的发展以及相关地学研究的需要。确定似大地水准面的方法可分为几何方法、重力学方法和组合法。几何方法就是根据几何关系测定一点的高程异常或者两点间似大地水准面高程差，例如天文水准、卫星测高和GPS/水准等。重力学方法是以一种或多种重力数据为边值，建立关于扰动位的相应重力边值问题，通过求解边值问题确定扰动位函数，再由Burns公式转换为大地水准面高， 按照莫洛金斯基理论求定高程异常。组合法则是同时利用几何水准数据和重力数据来确定似大地水准面，在计算过程中一般采用移去-恢复技术方法。就目前来说，确定区域高程异常采用最多的方法是GPS/水准法和组合法。(I)利用GPS/水准计算高程异常。如果在一个点上采用GPS观测技术精确地求定该点的大地高， 同时又通过精密水准求得该点正常高，两者之差即为该点的高程异常，具体精度取决于GPS 测定大地高的精度和水准测量正常高的精度。如果在一个区域布设有足够多的GPS水准点，那么通过计算这些点的高程异常值，选择合适的数学模型，便可以建立该区域高程异常计算模型。这是目前确定省市级(区域)高程异常的常用方法。(2)组合法目前重力学方法确定的重力似大地水准面具有较高分辨率但精度较低，而GPS/水准确定的几何似大地水准面精度高但是分辨率较低。组合法以GPS/水准等确定的高精度但分辨率较低的几何大地水准面作为控制，将移去-恢复方法确定的高分辨率但精度较低的重力似大地水准面与之拟合，以达到精化局部大地水准面的目的。对于超大区域的似大地水准面精化工作，为了提高拟合精度，一般采用分区拟合法。我国广泛采用GPS/水准方法来计算高程异常，此方法应用在平坦或者地形不复杂的测区效果良好，并获得了令人满意的结果；但高程异常拟合的精度受到测区大小、地形类别及已知高程点数量、分布等诸多条件的限制。超大区域似大地水准面地形复杂，高程异常的变化程度比较大，单一的拟合模型难以满足精度要求和实际工程的需要。对于超大区域的似大地水准面精化工作，常常采用分区拟合的方法。对于面积很大的超大区域，单一的拟合模型难以满足精度要求和实际工程的需要。因此，采用分区拟合的方法是解决这一问题的有效手段，其中，相邻区间的似大地水准面的平滑连接(区域接
4边问题)是非常关键的，有利于形成区域性高程异常的数字模型。分区拟合法，即将整个测区分成两个或两个以上的小区，根据每个小区的实际情况，选用适合自己的数学模型进行拟合，以满足求取高程异常的精度要求，然后将拟合后相邻小区域的高程异常模型用数学函数平滑连接成一体，形成区域性高程异常的数字模型。(I)分区应遵循一定的分区原则，分区的选择具有不确定性。目前分区主要有两类方法，一是根据测区的实际地形起伏情况划分，二是根据控制点的高程异常值来划分。分区工作需要较丰富的工程经验，且分区结果不同，对整个区域的拟合结果和拟合精度有很大影响；(2)分区拟合之后的“区域接边问题”，很难处理。相邻小区间的似大地水准面的连接(即区域接边问题)，对于形成完整平滑的区域性高程异常的数字模型是非常关键的， 由于小区间的实际情况不同导致使用的拟合模型种类多样，目前仍未出现较为通用的平滑连接方法。区域接边问题是分区拟合法急需解决的技术难题。针对目前分区拟合方法存在的不足提供一种基于无缝分区技术的超大区域似大地水准面精化方法，通过充分研究和利用格网技术的优点，提出基于格网技术，建立一套特殊的格网搜索机制，采用一定的工作流程来进行超大区域的似大地水准面精化工作。一种基于无缝分区技术的超大区域似大地水准面精化方法，包括以下步骤步骤I :数据采集；先确定区域范围，再在区域内收集或测量得到S个测量控制点的信息，包括点位平面坐标(X，Y)，GPS高程He，水准高程Htl,并使S个测量控制点均匀分布于整个区域。计算各测量控制点的高程异常ξ =Hg-H0 ;测量控制点的平均密度要求为每 IOOOkm2的测量控制点数量应大于等于I. O ;整个区域的测量控制点数量S应大于等于40 ；步骤2 :整个区域格网划分；采用等间距的正方形格网划分；格网大小应该根据已知控制点密度情况，以及整个区域高程异常ξ变化情况而定；步骤3:依次计算格网模型中各格网点的高程异常ξ ；步骤4 :建立格网模型；以格网形式建立超大区域的似大地水准面模型；在计算出所有格网节点的高程异常值之后，将所有格网节点的平面位置及其高程异常(Xk，Yk，Ik)保存在一个数据库文件之中；步骤5 :格网内插；格网模型建立起来后，采用双线性内插的方法，可以内插计算出该超大区域内任意点P的高程异常；具体过程为先找到该点所在的网格，再利用该网格 4个格网点的高程异常值，便可以内插得到P点的高程异常；由1、2两点内插出P1的值，由
3、4两点内插出P2的值，再由PpP2两点内插出P点的值；计算公式为
发明内容
权利要求
1.一种基于无缝分区技术的超大区域似大地水准面精化方法，其特征在于，包括以下步骤步骤I:数据采集；先确定区域范围，再在区域内收集或测量得到S个测量控制点的信息，包括点位平面坐标(X，Y)，GPS高程He，水准高程Htl,并使S个测量控制点均匀分布于整个区域。计算各测量控制点的高程异常I = He-Htl;测量控制点的平均密度要求为每 IOOOkm2的测量控制点数量应大于等于I. O ;整个区域的测量控制点数量S应大于等于40 ；步骤2 :整个区域格网划分；采用等间距的正方形格网划分；格网大小应该根据已知控制点密度情况，以及整个区域高程异常I变化情况而定；步骤3:依次计算格网模型中各格网点的高程异常ξ ；步骤4 :建立格网模型；以格网形式建立超大区域的似大地水准面模型；在计算出所有格网节点的高程异常值之后，将所有格网节点的平面位置及其高程异常(Xk，Yk，Ik)保存在一个数据库文件之中；步骤5 :格网内插；格网模型建立起来后，采用双线性内插的方法，可以内插计算出该超大区域内任意点P的高程异常；具体过程为先找到该点所在的网格，再利用该网格4个格网点的高程异常值，便可以内插得到P点的高程异常；由1、2两点内插出P1的值，由3、4 两点内插出P2的值，再由PpP2两点内插出P点的值；计算公式为
2.根据权利要求I所述的方法，其特征在于，所述步骤3，对于每个格网点，其高程异常计算过程分为3步(1)格网搜索；假设现在要计算的格网点为N，以N点为中心，按照格网编号逐步扩展进行搜索，目标是在所搜索的所有格网之内，其包含的测量控制点数量3,大于等于Stl 的取值范围为8-16，应根据已知控制点密度情况而定；具体搜索方法为以N点为中心，首先在 I号格网区域中搜索控制点，判断测量控制点数量是否满足条件，如满足则停止搜索；不满足的话，继续扩展到2号格网区域中搜索；依此类推，直到搜索到的测量控制点数量满足要求，才停止搜索；(2)小区域拟合；搜索结束后，该搜索区域即为格网点N的拟合区域；小区域拟合采用二次多项式拟合，其拟合方程为ξ i = a0+a1Xi+a2Yi+a3Xi2+a4Yi2+a5XiYi (i = 1,2, SN)式中，aj(j = 0,1,…，5)为拟合方程的待定系数；SN为已知测量控制点的个数；(Xi, Yi)为已知测量控制点的坐标值；(ξ D为已知测量控制点的高程异常；利用SNf已知测量控制点信息(Xi, Yi, ξ,)按上式进行回归分析，可以得到拟合方程的待定系数aj(j = 0,1,…，5)，从而可以得到该小区域的拟合模型；(3)计算格网点高程异常；利用上面刚刚建立好的二次多项式模型，计算该格网点N的高程异常。设格网节点N的坐标为(X，Y)，其高程异常值计算公式为ζ W - &Ο+&1Χ+&2Υ+&3Χ +&4Υ +&5ΧΥ式中，aj(j = 0，l，…，5)为上一步通过回归分析得到的拟合方程的系数；(X，Y)为格网节点N的坐标值；(ξ N)为格网节点N的高程异常。
全文摘要
本发明公开了一种基于无缝分区技术的超大区域似大地水准面精化方法，包括以下步骤步骤1数据采集；步骤2整个区域格网划分；步骤3依次计算格网模型中各格网点的高程异常ξ；步骤4建立格网模型；步骤5格网内插；采用本发明方法，无需进行人工分区，也不存在“区域接边问题”。该模型在应用时，直接进行格网内插即可。如图2所示，采用双线性内插的方法可以内插计算出该区域某网格内任意点的高程异常。
文档编号G01C5/00GK102607513SQ20121003197
公开日2012年7月25日 申请日期2012年2月14日 优先权日2012年2月14日
发明者王浩, 胡伍生, 胡杰 申请人:东南大学