大区域InSAR形变结果的自动化平差及镶嵌方法及装置的制造方法
【专利摘要】本发明实施例提供大区域InSAR形变结果的自动化平差及镶嵌方法及装置，所述方法包括：对大区域InSAR相邻景SAR影像上同名高相干点进行自动识别；利用自动识别出的同名高相干点，建立所述大区域形变结果的平差模型，生成平差方程；对所述平差方程进行解算，获取所述大区域中每个单景覆盖区的形变速率的平差改正量；利用解算出的所述大区域中多个单景覆盖区的形变速率的平差改正量，对所述大区域形变结果进行自动镶嵌。上述技术方案具有如下有益效果：能自动化地实现，可快速得到大区域协调一致的地表形变信息，从而更好的实现大区域地表形变监测。
【专利说明】
大区域I nSAR形变结果的自动化平差及镶嵌方法及装置
技术领域
[0001 ]本发明涉及InSAR(Synthetic Aperture Radar Interferometry，合成孔径雷达 干涉测量)技术领域，尤其涉及一种大区域InSAR形变结果的自动化平差及镶嵌方法及装 置。
【背景技术】
[0002] 基于卫星重复轨道的差分合成孔径雷达干涉测量（D-InSAR，Dif ferential Interferometric SAR)技术是上世纪80年代出现的地表形变测量领域的新技术，不仅可取 得厘米级的精度，而且能获取面上的信息，相较于分层标、GPS(Global Positioning System，全球定位系统)及水准测量等地面沉降技术有很大的优势。传统的D-InSAR受制于 空间失相干、时间失相干和大气干扰等因素，难以保障其应用的稳定性。在20世纪90年代末 出现的时间序列InSAR技术则完全突破了这些限制因素，是对D-InSAR技术的一次革新。时 间序列InSAR技术的特点是:基于时间序列SAR影像进行处理;处理的对象不是整幅影像的 全部像元，而是其中具有稳定散射特性或在较长时间间隔内保持高相干的像元子集(下称 "高相干点"）。时间序列InSAR技术总体上可以概括为两类：以永久散射体干涉(Permanent Scatterer或者Persistent Scatterer Interferometry或PS-InSAR)为代表的单一主影像 时间序列InSAR技术和以小基线集技术（Small baseline subset interferometry，或SBAS InSAR)为代表的多主影像时间序列InSAR技术。时间序列InSAR技术已经在火山、地震、滑 坡、地面沉降等引起的形变监测中得到大量应用。
[0003] 地面沉降是在自然和人为因素作用下，由于地壳表层土体压缩而导致区域性地面 标高降低的一种地质现象，是一种特殊的地表形变，在我国广泛发生。20世纪90年代初，上 海、天津、北京、江苏、浙江、河北等16省（区、市)地面沉降面积约为48700km 2，到2003年已达 至Ij93855km2,形成了长江三角洲、华北平原、珠江三角洲及汾渭盆地等4个地面沉降灾害严 重区。可见，地面沉降已经成为一种区域性灾害现象。目前国内外已经开展的时间序列 InSAR形变监测研究，绝大多数仅限于一景SAR影像所覆盖的相对较小的区域。针对需要由 多景SAR影像特别是多轨(track) SAR影像覆盖的大区域的InSAR形变监测研究很少。
[0004] 如何在单景时间序列InSAR形变信息提取的基础上，进行多景SAR影像形变结果的 平差及镶嵌，以实现大区域形变监测，这是本领域技术人员亟待解决的关键问题。
【发明内容】
[0005] 本发明实施例提供一种大区域InSAR形变结果的自动化平差及镶嵌方法及装置， 以解决多景SAR影像形变结果的自动镶嵌问题，从而更好的实现大区域地表形变监测。
[0006] -方面，本发明实施例提供了一种大区域InSAR形变结果的自动化平差及镶嵌方 法，所述方法包括：
[0007] 对大区域InSAR相邻景SAR影像上同名高相干点进行自动识别；
[0008] 利用自动识别出的同名高相干点，建立所述大区域形变结果的平差模型，生成平 差方程；
[0009] 对所述平差方程进行解算，获取所述大区域中每个单景覆盖区的形变速率的平差 改正量；
[0010] 利用解算出的所述大区域中多个单景覆盖区的形变速率的平差改正量，对所述大 区域形变结果进行自动镶嵌。
[0011] 另一方面，本发明实施例提供了一种大区域InSAR形变结果的自动化平差及镶嵌 装置，所述装置包括：
[0012] 自动识别单元，用于对大区域InSAR相邻景SAR影像上同名高相干点进行自动识 别；
[0013] 平差模型建立单元，用于利用自动识别出的同名高相干点，建立所述大区域形变 结果的平差模型，生成平差方程；
[0014] 平差方程解算单元，用于对所述平差方程进行解算，获取所述大区域中每个单景 覆盖区的形变速率的平差改正量；
[0015] 自动镶嵌单元，用于利用解算出的所述大区域中多个单景覆盖区的形变速率的平 差改正量，对所述大区域形变结果进行自动镶嵌。
[0016] 上述技术方案具有如下有益效果：能自动化地实现，可快速得到大区域协调一致 的地表形变信息，从而更好的实现大区域地表形变监测。
【附图说明】
[0017] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以 根据这些附图获得其他的附图。
[0018] 图1为本发明实施例一种大区域InSAR形变结果的自动化平差及镶嵌方法流程图；
[0019] 图2是本发明实施例使用的4景SAR影像覆盖大区域示意图（斜线、点、黑色等填充 区分别表示两景影像、三景影像、4四景影像的重叠区域）；
[0020] 图3是本发明实施例使用的Radarsat-2影像覆盖情况及水准点分布；
[0021 ]图4是本发明实施例得到的京津冀大区域地表沉降速率镶嵌图；
[0022]图5是常规处理得到的京津冀大区域地表沉降速率镶嵌图；
[0023] 图6是本发明应用实例142个检核点上InSAR和水准测量获取的沉降速率之差的直 方图；
[0024] 图7为本发明实施例一种大区域InSAR形变结果的自动化平差及镶嵌装置结构示 意图。
【具体实施方式】
[0025]下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0026] 如图1所示，为本发明实施例一种大区域InSAR形变结果的自动化平差及镶嵌方法 流程图，所述方法包括：
[0027] 101、对大区域InSAR相邻景SAR影像上同名高相干点进行自动识别；
[0028] 优选地，所述对大区域InSAR相邻景SAR影像上同名高相干点进行自动识别之前， 所述方法包括:对所述大区域InSAR的SAR影像使用同种类型的数字高程模型数据DEM进行 地理编码。
[0029] 同名高相干点的识别是平差及镶嵌的基础。虽然利用影像匹配技术能通过相邻 SAR振幅影像的匹配完成同名高相干点的自动识别，但耗时太长。事实上，利用地理编码后 同名高相干点有相同大地坐标这一原理即可快速解决这一问题。必须注意的是所有SAR影 像的地理编码应该使用同种类型的DEM，比如SRTM DEM，这样能保证相邻景SAR影像地理编 码后有相同的平面精度。图2中显示了4景SAR影像覆盖大区域的示意图，有更多影像覆盖更 大区域的情形可以依此类推。基于地理编码后同名高相干点的判断条件为：
[0030] |Pl(x，y)-P2(x，y)|<(Lx，L y) (1)
[0031] 其中（x，y)表示地理编码后的大地坐标，？1(^7)，？2(^7)代表相邻5六1?影像上位 于(x，y)的高相干点，（x，y)表示地理编码后的大地坐标，L x，Ly分别表示地理编码后东西向 和南北向的输出像元间隔。满足条件(1)的高相干点被认为是同名点。
[0032] 102、利用自动识别出的同名高相干点，建立所述大区域形变结果的平差模型，生 成平差方程；
[0033] 优选地，所述平差方程由两类方程构成:控制点平差方程和同名高相干点平差方 程;当所述大区域InSAR内没有控制点时，则所述平差方程由同名高相干点平差方程构成。
[0034] 在每个单景覆盖区通过时间序列InSAR技术求得的是相对于该景中某个点的相对 形变速率。在没有指定形变速率的基准点（参考点)情况下，利用时间序列InSAR技术求得的 线性形变速率与真实值之间存在一个系统偏差。这个系统偏差可以在大区域形变结果的平 差过程中求得。对上述4景影像覆盖的情形，假设单景影像的形变速率为Vi(x，y)，平差改正 量为Ai，? = 1，···，4。平差模型中的方程由两类构成，一类是控制点平差方程，另一类是同 名高相干点平差方程。在外部地面控制点实测数据支持下，建立控制点平差方差，即：
[0035] Vi(x，y)+Δ i=GCp_V(x，y) (2)
[0036] 其中GCp_V(x，y)表示在控制点(x，y)处经由水准或者其他手段测量得到的形变速 率值。它们是InSAR获取的形变速率的参考数据。
[0037] 对于相邻景影像的同名高相干点，建立同名高相干点平差方程，即：
[0038] Vi(x，y)+Δ i =Vj(x，y)+Δ j，i，j E {I，2，3，4} (3)
[0039] 这个方程表示平差后不同影像上的同名点应该具有一致的形变速率。图2中分别 以斜线、点、黑色等模式分别标记了两景影像、三景影像、四景影像的重叠区域。当点（X，y) 位于两景、三景、四景影像的重叠区内时，形如(3)式的平差方程分别有1个、2个、3个。
[0040] 103、对所述平差方程进行解算，获取所述大区域中每个单景覆盖区的形变速率的 平差改正量；
[0041]优选地，所述对所述平差方程进行解算，获取所述大区域中每个单景覆盖区的形 变速率的平差改正量，包括:根据最小二乘原理对所述平差方程进行解算，获取所述大区域 中每个单景覆盖区的形变速率的平差改正量。
[0042] 将上述两类方程联合在一起，构成了关于4个未知数（Δ 1，Δ 2Δ 3, Δ4)的线性方 程组。因为同名点数目很多，观测方程的个数远远多于未知数个数，根据最小二乘原理即可 求解未知数，这样同时完成了每景覆盖区的形变速率的平差与绝对定标。此处没有考虑控 制点的权重。在实际应用会出现不同的控制点有不同可信度的情形。此时，可根据控制点的 可信度在式(2)中设置不同的权重，在最小二乘解算中予以考虑，以提高解算精度。
[0043] 当监测区没有控制点时，则只有第二类方程。此时的平差可以任一覆盖区的形变 速率值为基准，未知数中有一个自由度，因此个数减为3个。这样平差完成后的形变速率和 真实值之间依然存在一个系统性偏差。
[0044] 104、利用解算出的所述大区域中多个单景覆盖区的形变速率的平差改正量，对所 述大区域形变结果进行自动镶嵌。
[0045] 优选地，所述利用解算出的所述大区域中多个单景覆盖区的形变速率的平差改正 量，对所述大区域形变结果进行自动镶嵌，包括:利用解算出的所述大区域中多个单景覆盖 区的形变速率的平差改正量，根据地理编码后的大地坐标，及同名点的重叠次数，对所述大 区域形变结果进行自动镶嵌。
[0046] 在完成了多个单景覆盖区的形变结果平差之后，根据地理编码后的大地坐标，可 完成大区域形变结果的自动镶嵌。在空间位置(X，y)处的形变结果为
[0047] (A)
[0048] 式中N表示在位置(x，y)处高相干点的重叠次数。对应于图2中，在斜线、点状、黑色 区域内，N的取值分别是2、3、4。而在其他区域，没有重叠，只有一个测量值，N的取值为1。这 样就获得了大区域内全部高相干点上的形变速率值。
[0049] 为了更好地说明本发明技术方案的有效性和优越性，现对本发明应用实例与常规 处理方法进行如下对比分析:如图3所示，是本发明实施例使用的Radarsat-2影像覆盖情况 及水准点分布图，监测区为京津冀地区，面积约10万平方公里，共8景Radarsat-2宽模式影 像覆盖，单景影像的覆盖区域约为150km*150km。这8个单景覆盖区上的时间序列SAR影像数 目及起始获取时间如表1所示。在利用InSAR进行地表沉降监测的同时，从北京市及天津市 的测绘部门获得了一些水准测量数据。北京和天津的水准测量测量工作每年10月初至11月 底完成，共获得了2012年至2014年的三期水准数据。请注意水准点的有效时段与SAR影像的 有效时段并不完全重合。这些水准点的位置以实心圆点标示于图3中，其中1/3的点表示参 与平差计算的水准点，共70个，2/3的点表示用于精度验证的水准点，共142个。
[0050] 表1京津冀地表沉降监测所用的Radarsat-2影像信息
[0052]利用时间序列InSAR技术得到每景覆盖区上雷达视线向的平均形变速率，然后转 换成竖直方向的地表沉降速率。其中负值表示地表下沉，正值表示地表抬升。然后按所述平 差方法进行处理，得到标定后的大区域镶嵌的平均沉降速率图，如图4所示，不同单景覆盖 区之间的沉降速率具有很好的一致性，这表明本发明方法是有效和可靠的。
[0053]图5为常规处理方法得到的京津冀大区域形变速率镶嵌图，由于仅有北京和天津 市域的水准测量数据，因此仅能对图3中标号为1，2,3,4的三景形变速率进行直接定标，余 下的4景形变速率只能依靠与已经定标的数据之间的同名点进行传递实现间接定标。最后 在生成大区域的镶嵌图时由于不同景的定标精度不一致将会而出现明显的"拼接缝"。可以 看见，图5与图4结果相比存在明显的误差。
[0054]利用142个水准点对图3所示的大区域沉降速率结果进行了精度验证。二者差值的 直方图如图6所示，差值的均方差为4.5毫米每年，差值最大值为12.1毫米每年，最小值为-11.2毫米每年。其中90个点上的沉降速率之差位于[-3mm/year，3mm/year]，这也说明本发 明应用实例得到的京津冀大区域沉降速率结果具有很高的精度。当然，SAR影像的起始时间 与水准点的起始时间并不完全重合，这也会对二者计算得到的沉降速率的差异造成一定影 响。
[0055] 综上可见，本发明应用实例大区域InSAR形变结果的自动化平差及镶嵌方法有两 大优点，其一是能自动化地实现，可快速得到大区域协调一致的地表形变信息，其二是将同 名点和外部控制点联合考虑，同步实现InSAR形变监测结果的平差和绝对定标。该方法在大 范围的地表形变监测中具有重要的应用价值。
[0056] 如图7所示，为本发明实施例一种大区域InSAR形变结果的自动化平差及镶嵌装置 结构示意图，所述装置包括：
[0057] 自动识别单元71，用于对大区域InSAR相邻景SAR影像上同名高相干点进行自动识 别；
[0058]平差模型建立单元72，用于利用自动识别出的同名高相干点，建立所述大区域形 变结果的平差模型，生成平差方程；
[0059] 平差方程解算单元73,用于对所述平差方程进行解算，获取所述大区域中每个单 景覆盖区的形变速率的平差改正量；
[0060] 自动镶嵌单元74,用于利用解算出的所述大区域中多个单景覆盖区的形变速率的 平差改正量，对所述大区域形变结果进行自动镶嵌。
[0061] 优选地，所述装置还包括:地理编码单元70,用于对所述大区域InSAR的SAR影像使 用同种类型的数字高程模型数据DEM进行地理编码。
[0062] 优选地，所述平差方程由两类方程构成:控制点平差方程和同名高相干点平差方 程;当所述大区域InSAR内没有控制点时，则所述平差方程由同名高相干点平差方程构成。 [0063]优选地，所述平差方程解算单元73,进一步具体用于根据最小二乘原理对所述平 差方程进行解算，获取所述大区域中每个单景覆盖区的形变速率的平差改正量。
[0064] 优选地，所述自动镶嵌单元74,进一步具体用于利用解算出的所述大区域中多个 单景覆盖区的形变速率的平差改正量，根据地理编码后的大地坐标，及同名点的重叠次数， 对所述大区域形变结果进行自动镶嵌。
[0065] 上述技术方案具有如下有益效果：能自动化地实现，可快速得到大区域协调一致 的地表形变信息，从而更好的实现大区域地表形变监测。
[0066] 应该明白，公开的过程中的步骤的特定顺序或层次是示例性方法的实例。基于设 计偏好，应该理解，过程中的步骤的特定顺序或层次可以在不脱离本公开的保护范围的情 况下得到重新安排。所附的方法权利要求以示例性的顺序给出了各种步骤的要素，并且不 是要限于所述的特定顺序或层次。
[0067] 在上述的详细描述中，各种特征一起组合在单个的实施方案中，以简化本公开。不 应该将这种公开方法解释为反映了这样的意图，即，所要求保护的主题的实施方案需要比 清楚地在每个权利要求中所陈述的特征更多的特征。相反，如所附的权利要求书所反映的 那样，本发明处于比所公开的单个实施方案的全部特征少的状态。因此，所附的权利要求书 特此清楚地被并入详细描述中，其中每项权利要求独自作为本发明单独的优选实施方案。
[0068] 为使本领域内的任何技术人员能够实现或者使用本发明，上面对所公开实施例进 行了描述。对于本领域技术人员来说;这些实施例的各种修改方式都是显而易见的，并且本 文定义的一般原理也可以在不脱离本公开的精神和保护范围的基础上适用于其它实施例。 因此，本公开并不限于本文给出的实施例，而是与本申请公开的原理和新颖性特征的最广 范围相一致。
[0069] 上文的描述包括一个或多个实施例的举例。当然，为了描述上述实施例而描述部 件或方法的所有可能的结合是不可能的，但是本领域普通技术人员应该认识到，各个实施 例可以做进一步的组合和排列。因此，本文中描述的实施例旨在涵盖落入所附权利要求书 的保护范围内的所有这样的改变、修改和变型。此外，就说明书或权利要求书中使用的术语 "包含"，该词的涵盖方式类似于术语"包括"，就如同"包括，"在权利要求中用作衔接词所解 释的那样。此外，使用在权利要求书的说明书中的任何一个术语"或者"是要表示"非排它性 的或者"。
[0070] 本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块 (illustrative logical block)，单元，和步骤可以通过电子硬件、电脑软件，或两者的结 合进行实现。为清楚展示硬件和软件的可替换性（：!^丨6代1^1^6313；[1；^7)，上述的各种说明 性部件（illustrative components)，单元和步骤已经通用地描述了它们的功能。这样的功 能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员 可以对于每种特定的应用，可以使用各种方法实现所述的功能，但这种实现不应被理解为 超出本发明实施例保护的范围。
[0071] 本发明实施例中所描述的各种说明性的逻辑块，或单元都可以通过通用处理器， 数字信号处理器，专用集成电路(ASIC)，现场可编程门阵列或其它可编程逻辑装置，离散门 或晶体管逻辑，离散硬件部件，或上述任何组合的设计来实现或操作所描述的功能。通用处 理器可以为微处理器，可选地，该通用处理器也可以为任何传统的处理器、控制器、微控制 器或状态机。处理器也可以通过计算装置的组合来实现，例如数字信号处理器和微处理器， 多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或任何其它类似的配置 来实现。
[0072] 本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软 件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储 器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介 中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信息，并 可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储媒介可 以设置于ASIC中，ASIC可以设置于用户终端中。可选地，处理器和存储媒介也可以设置于用 户终端中的不同的部件中。
[0073] 在一个或多个示例性的设计中，本发明实施例所描述的上述功能可以在硬件、软 件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的 媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电 脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以 是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但 不限于RAM、ROM、EEPROM、⑶-ROM或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任 何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理 器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例 如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双 绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的 电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、错射盘、光盘、DVD、软盘 和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述的组合 也可以包含在电脑可读媒介中。
[0074]以上所述的【具体实施方式】，对本发明的目的、技术方案和有益效果进行了进一步 详细说明，所应理解的是，以上所述仅为本发明的【具体实施方式】而已，并不用于限定本发明 的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含 在本发明的保护范围之内。
【主权项】
1. 一种大区域InSAR形变结果的自动化平差及镶嵌方法，其特征在于，所述方法包括： 对大区域InSAR相邻景SAR影像上同名高相干点进行自动识别； 利用自动识别出的同名高相干点，建立所述大区域形变结果的平差模型，生成平差方 程； 对所述平差方程进行解算，获取所述大区域中每个单景覆盖区的形变速率的平差改正 量； 利用解算出的所述大区域中多个单景覆盖区的形变速率的平差改正量，对所述大区域 形变结果进行自动镶嵌。2. 如权利要求1所述大区域InSAR形变结果的自动化平差及镶嵌方法，其特征在于，所 述对大区域InSAR相邻景SAR影像上同名高相干点进行自动识别之前，所述方法包括： 对所述大区域InSAR的SAR影像使用同种类型的数字高程模型数据DEM进行地理编码。3. 如权利要求1所述大区域InSAR形变结果的自动化平差及镶嵌方法，其特征在于，所 述平差方程由两类方程构成：控制点平差方程和同名高相干点平差方程；当所述大区域 InSAR内没有控制点时，则所述平差方程由同名高相干点平差方程构成。4. 如权利要求1所述大区域InSAR形变结果的自动化平差及镶嵌方法，其特征在于，所 述对所述平差方程进行解算，获取所述大区域中每个单景覆盖区的形变速率的平差改正 量，包括： 根据最小二乘原理对所述平差方程进行解算，获取所述大区域中每个单景覆盖区的形 变速率的平差改正量。5. 如权利要求1所述大区域InSAR形变结果的自动化平差及镶嵌方法，其特征在于，所 述利用解算出的所述大区域中多个单景覆盖区的形变速率的平差改正量，对所述大区域形 变结果进行自动镶嵌，包括： 利用解算出的所述大区域中多个单景覆盖区的形变速率的平差改正量，根据地理编码 后的大地坐标，及同名点的重叠次数，对所述大区域形变结果进行自动镶嵌。6. -种大区域InSAR形变结果的自动化平差及镶嵌装置，其特征在于，所述装置包括： 自动识别单元，用于对大区域InSAR相邻景SAR影像上同名高相干点进行自动识别； 平差模型建立单元，用于利用自动识别出的同名高相干点，建立所述大区域形变结果 的平差模型，生成平差方程； 平差方程解算单元，用于对所述平差方程进行解算，获取所述大区域中每个单景覆盖 区的形变速率的平差改正量； 自动镶嵌单元，用于利用解算出的所述大区域中多个单景覆盖区的形变速率的平差改 正量，对所述大区域形变结果进行自动镶嵌。7. 如权利要求6所述大区域InSAR形变结果的自动化平差及镶嵌装置，其特征在于，所 述装置还包括： 地理编码单元，用于对所述大区域InSAR的SAR影像使用同种类型的数字高程模型数据 DEM进行地理编码。8. 如权利要求6所述大区域InSAR形变结果的自动化平差及镶嵌装置，其特征在于，所 述平差方程由两类方程构成：控制点平差方程和同名高相干点平差方程；当所述大区域 InSAR内没有控制点时，则所述平差方程由同名高相干点平差方程构成。9. 如权利要求6所述大区域InSAR形变结果的自动化平差及镶嵌装置，其特征在于， 所述平差方程解算单元，进一步具体用于根据最小二乘原理对所述平差方程进行解 算，获取所述大区域中每个单景覆盖区的形变速率的平差改正量。10. 如权利要求6所述大区域InSAR形变结果的自动化平差及镶嵌装置，其特征在于， 所述自动镶嵌单元，进一步具体用于利用解算出的所述大区域中多个单景覆盖区的形 变速率的平差改正量，根据地理编码后的大地坐标，及同名点的重叠次数，对所述大区域形 变结果进行自动镶嵌。
【文档编号】G01S13/90GK105842691SQ201610154331
【公开日】2016年8月10日
【申请日】2016年3月17日
【发明人】张永红, 吴宏安, 康永辉
【申请人】中国测绘科学研究院