基于快速Huynen-Euler分解的受灾建筑物损坏等级估计方法及系统

技术领域

本发明涉及全极化合成孔径雷达数据的极化分解领域，特别涉及基于快速Huynen-Euler分解的受灾建筑物损坏等级估计方法及系统。

背景技术

随着微波技术的发展，能够有效获取目标的全极化散射特性的极化合成孔径雷达作为微波遥感的重要手段之一，被广泛应用于目标检测与识别、地物分类、地表参数反演等方面。近年来，由于自然灾害的频发及其程度的加剧，灾害监测和建筑区域损坏评估作为极化合成孔径雷达的新兴应用引起了广泛关注，因为其能够全天时、全天候、不受云雨气象条件的影响、经由遥感的方式迅速的获取大范围的有关受灾区域的数据。经由比较并分析灾前与灾后该区域的散射特性的变化，可以将其准确识别为损坏区域、水淹区域、未受损区域，等等。模型分解提取的二次散射功率是目前最常用的进行受灾区域状态判断的参数，因为其被认为是由建筑竖直的墙体和地表构成的二面角结构所产生的。故已有算法通过比较二次散射功率的变化来判断建筑区域受损程度，如文献1[S.-W.Chen,X.-S.Wang,andM.Sato,“Urban damage level mapping based on scattering mechanisminvestigation using fully polarimetric SAR data for the 3.11east japanearthquake,”IEEE Trans.Geosci.Remote Sens.,vol.54,no.12,pp.6919–6929,2016.]。然而，采用模型分解得到二次散射功率的过程相对耗时，因为在对这一参数进行求解的过程中需要经过去取向、同极化功率相对大小判断、根据判断结果选择求解方法、消除负功率影响等多个步骤，而灾害应用对参数提取的时效性要求很高——为了保证能够迅速的检测到灾害、对城区受灾程度进行快速评估以制定高效及时的救援措施，挽救人民群众生命安全和财产。

Huynen-Euler分解是Huynen提出的对目标极化散射矩阵进行的一种分解算法(文献2[J.R.Huynen,“Phenomenological theory of radar targets,”Ph.D.dissertation,Dept.Elect.Eng.,Math.Comput.Sci.,Tech.Univ.Delft,Delft,The Netherlands,1970.])，分解得到的五个欧拉参数能够有效反应目标不同方面的散射特性。其中跳跃角参数v与信号的反射次数相关，其本质上反应了与二次散射功率相同的信息。尽管理论上可行，Huynen提出的跳跃角参数在实际应用中存在角度模糊问题。采用快速Huynen-Euler分解可以通过解析的方式提取改进的跳跃角参数v

发明内容

本发明的目的在于解决自然灾害突发时对城区受灾程度进行评估的过程耗时较大的问题。

为了实现上述目的，本发明提供了一种基于快速Huynen-Euler分解的受灾建筑物损坏等级估计方法，所述方法包括：

根据全极化合成孔径雷达的极化数据，分别获得目标灾前的相干矩阵[T

根据目标灾前的相干矩阵[T

根据目标灾前Huynen参数和目标灾后的Huynen参数，利用快速Huynen-Euler分解，分别计算目标灾前的跳跃角参数v

根据目标灾前的跳跃角参数v

根据跳跃角变化率参数Δv

作为上述方法的一种改进，所述目标灾前的相干矩阵[T

其中，T1

目标灾后的相干矩阵[T

其中，T2

作为上述方法的一种改进，所述根据目标灾前的相干矩阵[T

将目标灾前和灾后的Huynen参数分别记为A

C

D

E

F

G

H

将目标灾后的Huynen参数分别记为A

C

D

E

F

G

H

其中，Re{·}表示取实部，Im{·}表示取虚部。

作为上述方法的一种改进，所述根据目标灾前Huynen参数和目标灾后的Huynen参数，利用快速Huynen-Euler分解，分别计算目标灾前的跳跃角参数v

计算目标灾前的跳跃角参数v

计算目标灾后的跳跃角参数v

其中，v

作为上述方法的一种改进，所述根据目标灾前的跳跃角参数v

计算跳跃角变化率参数Δv

其中，Δv

作为上述方法的一种改进，所述根据跳跃角变化率参数Δv

受灾建筑物损坏等级DL为：

其中，损坏等级DL∈[0,1]。

本发明还提出了一种基于快速Huynen-Euler分解的受灾建筑物损坏等级估计系统，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述的方法。

本发明的优势在于：

本发明提出的方法和系统解决了自然灾害突发时对城区受灾程度进行评估的过程耗时较大的问题，既避免了旧跳跃角参数的角度模糊问题，通过快速分解方法实现对建筑区域的快速识别，解决了自然灾害突发时对城区受灾程度进行评估的过程耗时较大的问题，能够得到对受灾地区受灾程度快速准确的估计。

附图说明

图1为本发明的基于快速Huynen-Euler分解的受灾建筑物损坏等级估计方法的流程图；

图2(a)为受灾地区日本石卷市的灾前的光学图像；

图2(b)为受灾地区日本石卷市的灾后的光学图像；

图3为建筑物损坏等级估计结果的示意图；

图4为利用模型分解方法得到的建筑物损坏等级估计结果。

具体实施方式

下面结合附图和具体实施例对本发明进行详细的说明。

Huynen-Euler分解提取的跳跃角参数v与信号的反射次数相关，理论上能够对城区进行识别，但在实际应用中其存在角度模糊的问题。采用快速Huynen-Euler分解可以通过解析的方式提取改进的跳跃角参数v

如图1所示，本发明的实施例1提出了的基于快速Huynen-Euler分解的受灾建筑物损坏等级估计方法，包括以下步骤：

步骤1)根据灾前和灾后的全极化合成孔径雷达的极化数据，分别获得目标灾前和灾后的相干矩阵[T

步骤2)根据步骤1)得到的灾前和灾后的Huynen参数A

步骤3)根据步骤2)得到的目标灾前和灾后的改进的跳跃角v

步骤4)根据步骤3)得到的目标的跳跃角变化率Δv

下面通过实例进一步描述。

在一个实例中，受灾地区日本石卷市的灾前与灾后的光学图像如图2(a)和图2(b)所示，该图像由Google Earth分别在2010年6月25日和2011年4月6日获取。在2011年3月11日，发生了9.0级的日本东北地方太平洋近海地震，并引发最大爬高40.1米的海啸，是日本史上最大的天然灾难。地震造成至少15893人死亡、2553人失踪、伤者6152人，遭受破坏的房屋1292417栋，为日本二战后伤亡最惨重的自然灾害。由图2可以看出，灾前灾后的沿海城区发生了明显变化，远离海边的地区则没有明显的变化。这是因为沿海城区受灾极其严重，多数房屋倒塌且残骸被随之而来的海啸冲走，而内陆地区受这两种灾害的影响则较小。该实施例中采用的全极化合成孔径雷达数据是由JAXA的星载ALOS PALSAR系统分别在2009年4月2日和2011年4月8日获取。

在步骤1)中，根据灾前和灾后的全极化合成孔径雷达的极化数据，分别获得目标灾前和灾后的相干矩阵[T

目标的相干矩阵[T]表示如下：

将目标灾前灾后的相干矩阵分别记为[T

Huynen参数A

其中，Huynen参数A

在步骤2)中，根据步骤1)得到的灾前和灾后的Huynen参数A

快速Huynen-Euler算法是本发明提出的一种能够快速实现Huynen-Euler分解并提取相应的Huynen-Euler参数的算法。Huynen-Euler分解是Huynen基于包含目标的全部极化信息的散射矩阵[S]提出的一种分解方法，表示如下：

其中，

各参数都具备特定的与目标散射特性相关的含义：m

已有的算法通过对[S]进行伪特征分解或对

然后再根据式(2)便可以获得目标的Huynen参数。根据式(2)、(3)和(6)，我们可以将Huynen参数用五个Huynen-Euler参数表示，如下：

利用式(7)可以推导出Huynen-Euler参数的解析解如下：

需要说明的是，原Huynen-Euler分解提取的跳跃角参数v对于建筑物的识别存在角度模糊问题，例如，对于两个散射矩阵为如下所示对角阵的目标[S

原Huynen-Euler分解提取的跳跃角参数v

[S

v

本发明采用v

因此，本发明根据式(8)-(11)和(16)，实现对目标的快速Huynen-Euler分解，并利用式(16)分别计算得到目标灾前和灾后的改进的跳跃角参数v

在步骤3)中，根据步骤2)得到的目标灾前和灾后的改进的跳跃角v

Δv

在步骤4)中，根据步骤3)得到的目标的跳跃角变化率Δv

受灾建筑物损坏等级DL的估计方法如下：

其中，Δv

在本实施例中，利用快速Huynen-Euler分解估计得到的受灾建筑物损坏等级如图3所示，可以看出，其能够准确表示严重受损的沿海城区的与非沿海区域的无损特点，而且不受海洋、山区和植被区域的干扰，即这些非城区得到的受灾程度估计结果都近乎为0。已有的利用模型分解得到的二次散射功率进行城区受损程度估计的算法的估计结果如图4所示，可以看出，其结果与本发明的结果很相近，二者的相干性高达0.7471，但很显然，已有算法不能避免海洋、山区和植被区域的干扰，其在这些地区都有较高的DL值。

除此之外，本发明的快速Huynen-Euler分解另一巨大优势在于求解的高效性。通过对本组实施例数据的计算时间记录(本实施例数据为复数矩阵1248×18432)，可以得到：本发明提出的算法得到计算DL值所需参数v

本发明的实施例2提出了一种基于快速Huynen-Euler分解的受灾建筑物损坏等级估计系统，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现实施例1的方法。

最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。