申请/专利权人：中南林业科技大学

申请日：2019-06-18

公开（公告）日：2023-01-10

公开（公告）号：CN110310370B

主分类号：G06T17/05

分类号：G06T17/05;G06T17/20

优先权：

专利状态码：有效-授权

法律状态：2023.01.10#授权;2019.11.01#实质审查的生效;2019.10.08#公开

摘要：一种基于积分平差实现GPS与SRTM点面融合的方法，包括以下步骤：1通过地面GPS观测，获得目标地区SRTM的DEM格网上任一点p的GPS高程值hp；2用测量平差方法把hp和p所在DEM格网的高程值进行融合；得到p点所在格网的格网值和四个角点a、b、c、d新的高程值；作为第一层融合；3以a、b、c、d四个点为新的观测值，分别对其周围的网格进行融合，得到外围8格网融合后新的格网值及外围12个节点的新的高程；作为第二层融合；依次类推，得到全格网的高程值；当有格网的角点得到多个高程值时，取平均值。本发明的基于积分平差实现GPS与SRTM点面融合的方法的精度高。

主权项：1.一种GPS与SRTM点面融合的方法，其特征在于，包括以下步骤：1通过地面GPS观测，获得目标地区SRTM的DEM格网上任一点p的GPS高程值hp；2用测量平差方法把hp和p所在DEM格网的高程值进行融合；得到p点所在格网的格网值和四个角点a、b、c、d新的高程值；作为第一层融合；3以a、b、c、d四个点为新的观测值，分别对其周围的格网进行融合，得到外围8格网融合后新的格网值及外围12个节点的新的高程；作为第二层融合；计算第二层8个格网平差融合后的新值与原格网值的差值，若各差值的均值阀值，则计算停止；若各差值的均值阀值的格网，则修正原DEM的高程值，直到各差值的均值阀值；依次类推，得到全格网的高程值；当有格网的角点得到多个高程值时，取平均值。

全文数据：一种基于积分平差实现GPS与SRTM点面融合的方法技术领域本发明涉及一种数字高程模型的处理方法，尤其涉及一种基于积分平差实现GPS与SRTM点面融合的方法。背景技术数字高程模型DigitalElevationModel，简称DEM是通过有限的地形高程数据实现对地形曲面的数字化模拟，它是对二维地理空间上具有连续变化特征地理现象的模型化表达和过程模拟。如今，DEM已被广泛应用，如土木工程、军事应用、地质灾害和可视化分析等。建立DEM实现地表重建的前提是对连续地球表面进行有限离散点的精确采样和三维坐标测量。目前，DEM数据获取方法主要包括野外测量、摄影测量和遥感、制图数字化、合成孔径雷达和机载激光扫描。如今，使用卫星图像生成DEM有一个显著的优势因为它相对便宜而且生成DEM所需的时间更少。然而使用光谱范围的缺点是它需要高分辨率，良好的光照条件和无云视线才能得到较好的DEM精度。近年来，干涉合成孔径雷达InSAR因其不依赖于自然光照而作为一个主动系统在提取高程数据方面变成流行技术。2000年美国航天飞机雷达地形测绘任务Shuttleradartopographymission,SRTM利用机载Insar技术，仅11d便获取了全球80％陆地1″分辨率的三维地形信息。2014年9月，1″分辨率的SRTM数据逐步向全球用户免费开放。SRTM在平坦区域的准确度范围为2.9-3.4米，半平坦区域为5.4-6.29米[1]A.KarwelandI.Ewiak,"EstimationoftheaccuracyoftheSRTMterrainmodelontheareaofPoland,"Theinternationalarchivesofthephotogrammetry,remotesensingandspatialinformationsciences,vol.37,no.B7,pp.169-172,2008.。已经有很多方法用来提高SRTM的精度，例如融合Aster和SRTM数据、融合ICESat数据和SRTM数据等；在我国，每个省市都建有CORS站，每天收集全省各地大量点的GPS高程或北斗高程，充分利用这些GPS高程数据具有重要的实际和实用价值。发明内容本发明要解决的技术问题是克服现有技术的不足，提供一种精度高的基于积分平差实现GPS与SRTM点面融合的方法。为解决上述技术问题，本发明提出的技术方案为：如图1所示，某地区SRTM的DEM，假设任意一点p所在DEM格网的高程值为通过地面GPS观测，获得p点的GPS高程值hp，我们的目标是把p点的GPS高程值hp与该地区的SRTMDEM进行融合，以达到提高DEM的精度。融合的基本思想是，首先用测量平差的方法把hp与进行融合，求得p点所在格网的格网值和四个角点a,b,c,d新的高程值。此为以p点为中心的中心融合。a,b,c,d四个点的高程又会分别影响与它们相邻外围的8个相邻网格的高程，因此可以把这4个点的高程看作它们外围相邻格网的新观测值，分别与相邻格网进一步融合，求得外围8个格网融合后新的格网值及外围12个节点的新的高程。此为以p点为中心的第二层融合。外围12个节点的新的高程又会影响围绕第二层的16个格网，逐个格网融合，就可完成以p点为中心的第三层融合。依此类推。p点的高程值hp通过这个方式可以提高其他格网的高程值精度。由于一个格网会有多个相邻格网，有些格网的角点会得到多个新值，可取平均作为最后值。以图1为例，在第二层融合过程中，格网i-1,j与格网i-1,j+1有共同的节点，里面的节点a的高程使用第一层的结果，外围节点的高程值则取第二层融合时，相邻两格网计算结果的平均值。具体的平差模型可按如下方法建立：一般而言，DEM中每个格网的高程值表示这个格网内所有点高程的平均值。用hx,y表示任意一点x,y的高程，表示格网i,j的格网值,那么可以用hx,y在这个格网内的积分的平均值表示，即：其中，S为格网面积。DEM的各格网值一般是通过各种测绘手段获取的的观测值hi,j。考虑观测误差，有：p点的高程可表示为hxp，yp，考虑到观测误差的存在，p点高程的观测值可表示为：hp+v＝hxp，yp3因为实际地形一般都非常复杂，故难以用数学函数准确描述，所以很难求得hx，y的精确表述。我们可以考虑一个非常小的局部，比方说，一个围绕格网i,j的3*3的格网范围，在这个小区域内，找一个数学函数来近似代替hx，y是可行的，例如可用一个二次函数代替hx，y等。利用这个区域9个格网点的观测值与P点的观测值，用测量平差的方法可联合求解这个区域的数学函数一个高程面估计，根据这个高程面，可以求得融合后格网i，j和四个角点的高程值，再由式1可以重新确定融合后格网i，j的高程。按式1可构成这个区域9个点的观测方程，再与式2联合，可形成如下平差模型：hp+vp＝hx，y5如果一个二次函数代替hx，y，就有hx，y＝a1+a2x+a3y+a4xy+a5x2+a6y26hp+v＝a1+a2xp+a3yp+a4xpyp+a5xp2+a6yp28最后可构成平差模型：L+V＝AX9其中L＝hi-1，j-1hi-1，j...hi+1，jhi+1，j+1hpTV＝Vi-1，j-1Vi-1，j...Vi+1，jVi+1，j+1VpTX＝a1a2a3a4a5a6T根据最小二乘平差，有其中P为权值，定权为900。平差融合后，小区域内各点的高程值可表示为格网i,j平差融合后的值分别为：格网i,j四个角点的融合后的高程分别为格网i,j平差融合后的值即为该格网融合后的格网值。格网i,j四个角点融合后的高程值对于相邻的格网就是新的观测值，如对于格网i-1,j，是2个新的观测，新观测值的权定为900。按与上述一样的方法，取格网i-1,j及周边9个格网值与这2个新观测一起进行平差，可以求得该格网融合后的高程值及格网4个角点融合后的高程值。具体可按如下步骤进行：1取p点及所在格网与周边8个格网点一起平差，定权为900，求得所在格网及4个角点融合后的高程值，融合后的格网值作为该格网最后的高程值，此为以p点为中心的第一层融合。2经过第一层融合后，围绕格网i,j的8个网格会受到格网i,j的4个角点高程变化的影响。其中格网i-1,j有a、c两个新观测点，格网i,j+1有a、b两个新观测点，格网i+1,j有b、c两个新观测，格网i,j-1有c、d两个新观测点。格网i-1,j+1有新观测点a，格网i+1,j+1有新观测点b，格网i+1,j有新观测点d，格网i-1,j-1有新观测点c。按上述平差方法分别将新观测点的观测与相应格网融合，可求得第二层8个格网的格网融合后的新值，及外围节点的新的高程值。对于相邻网格共同的节点取两个网格各自计算出的节点新高程值的平均值作为最后的新值。3计算第二层8个网格平差融合后的新值与原格网值的差值，若各差值的均值阀值0.5米的格网，表明原DEM有误差需要修正，则与该格网相关的计算可以进一步往下走进行，直至算法收敛。本发明提出的实现GPS与SRTM点面融合的积分平差算法，该算法的特点是把任意一点的高程表示成一个函数，然后DEM格网的高程可表达为这个函数的积分平均表达式，利用已知GPS实测点与相邻格网建立平差模型，再利用最小二乘解算出格网的高程改正值和四个角点的高程值，接着把四个角点当做下一次迭代的已知高程值进行联合解算，从而实现提高SRTM精度的目的。附图说明图1为DEM示意图。图2为实施例1中平坦区的SRTM影像图。图3为实施例1中半平坦区的SRTM影像图。图4为实施例1中陡峭区的SRTM影像图。图5为实施例1中已知GPS观测点对应的像素位置图。图6为实施例1中平坦区SRTM高程与GPS高程直线拟合图。图7为实施例1中平坦区算法改正高程与GPS高程直线拟合图。图8为实施例1中半平坦区SRTM高程与GPS高程直线拟合图。图9为实施例1中半平坦区算法改正高程与GPS高程直线拟合图。图10实施例1中为陡峭区SRTM高程与GPS高程直线拟合图。图11实施例1中为陡峭区算法改正高程与GPS高程直线拟合图。具体实施方式为了便于理解本发明，下文将结合较佳的实施例对本发明作更全面、细致地描述，但本发明的保护范围并不限于以下具体的实施例。除非另有定义，下文中所使用的所有专业术语与本领域技术人员通常理解的含义相同。本文中所使用的专业术语只是为了描述具体实施例的目的，并不是旨在限制本发明的保护范围。实施例1为了证明该算法的可行性与有效性，选取中国某地区为研究区域，从网站http:earthexplorer.usgs.gov上面下载该地区的SRTM数据，然后从该地区截取三个地理特征不同的区域作为实验数据，它们的大小均为64*64，三个区域的坡度分别在0-5度，5-20度，20度以上，如图2所示，分别对应平坦区、半平坦区、陡峭区，地理特征极具代表性。该区域已知高程点的高程来自于CORS站记录的该区域大量的GPS流动站点的高程。首先对三个研究区域分别任意选择一个已知的GPS观测点，已知GPS观测点与研究区域的位置关系如图3所示，然后利用该观测点获得的GPS高程值采用点面融合积分平差算法改进对应研究区域的SRTM影像。对于研究区域已知的其他GPS高程值则作为精度评定样本，则三个研究区域对应的GPS校验点分别为48、71、37个。为了客观地评估改进的DEM精度，考虑到系统偏差和标准差等精度估计会受到异常值和误差非正态分布的影响，因此这里考虑采用稳健统计方法，采用误差均值MeanError，ME、均方根误差RootMeanSquareError，RMSE、归一化绝对偏差中位数NormalizedMedian，NMAD和误差的标准差Standarddeviationoferror，SDE来进行精度评定。NMAD＝1.4826×medianj|Δhj-mΔh|其中，N是GPS精度评定样本的个数；Hi，Href分别是提出算法改进的GPS高程和GPS的高程；Δhi,mΔh分别是第i个精度评定点的高程差和高程差的中值；这些值都是精度的度量，越小表明改进算法效果越好。为避免偶然性，对每一个研究区域将进行多次独立重复实验。对研究区域GPS精度评定样本点与改进高程值比对，得到精度评定如表1所示，表中的精度为多次实验的平均值。表1三个研究区域精度评定由表1可得，在平坦区、过渡区和陡峭区采用所提算法得到的DEM具有明显优势，与SRTMDEM比较，改进的RMSE提高了1-2米，精度分别提高了19.5％、20.3％和11.6％，表明所提出算法能够有效改善SRTMDEM，特别是在平坦区域和和半平坦区域，但是在陡峭区域，因为受到地势的影响，融合的精度将受到一定程度的影响。为了更好地展示和分析所提出算法的优势，这里对精度评定GPS点的测量高程和采用所提出算法对这些点融合之后的高程之间的测量关系进行详细分析，分别列出GPS测量高程与相应SRTM的高程的直线拟合图，和GPS测量高程与提出算法改正高程的直线拟合图，如图6-图11所示。图6、图8和图10为三个研究区域表明精度评定点的SRTM高程值与GPS观测高程之间的相关性。图7、图9和图11反映的是算法改正高程与GPS观测高程之间的相关性。这六个图显示的均是线性与正斜率的关系，表明两两变量之间朝着同一个方向发展。图7、图9和图11分别与图6、图8和图10比拥有更高的相关性，而且散点的分布更加具有一致性，说明提出的点面融合的平差算法改进DEM效果明显。本实施例中，结果表明所提出点面融合的积分平差算法能够有效改善平坦区和和半平坦区的SRTMDEM，精度大约提高20％，但是在陡峭区域，因为受到地势的影响，融合的精度将受到一定程度的影响，精度大约能提高11％。

权利要求：1.一种基于积分平差实现GPS与SRTM点面融合的方法，其特征在于，包括以下步骤：1通过地面GPS观测，获得目标地区SRTM的DEM格网上任一点p的GPS高程值hp；2用测量平差方法把hp和p所在DEM格网的高程值进行融合；得到p点所在格网的格网值和四个角点a、b、c、d新的高程值；作为第一层融合；3以a、b、c、d四个点为新的观测值，分别对其周围的网格进行融合，得到外围8格网融合后新的格网值及外围12个节点的新的高程；作为第二层融合；计算第二层8个网格平差融合后的新值与原格网值的差值，若各差值的均值阀值的格网，则修正原DEM的高程值，直到各差值的均值阀值；依次类推，得到全格网的高程值；当有格网的角点得到多个高程值时，取平均值。2.根据权利要求1所述的基于积分平差实现GPS与SRTM点面融合的方法，其特征在于：所述DEM中每个格网的高程值表示这个格网内所有点高程的平均值；用hx,y表示任意一点x,y的高程，表示格网i,j的格网值,那么可以用hx,y在这个格网内的积分的平均值表示，即：其中，S为格网面积。3.根据权利要求2所述的基于积分平差实现GPS与SRTM点面融合的方法，其特征在于：所述DEM的各格网值一般是通过各种测绘手段获取的的观测值hi,j；考虑观测误差，有：p点的高程可表示为hxp，yp，考虑到观测误差的存在，p点高程的观测值可表示为：hp+v＝hxp，yp。34.根据权利要求1所述的基于积分平差实现GPS与SRTM点面融合的方法，其特征在于：将目标地区分成N个单元，每个单元独立的进行hp和p所在DEM格网的高程值的融合。

百度查询： 中南林业科技大学 一种GPS与SRTM点面融合的方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD