一种极化雷达建筑损毁信息提取方法技术

本发明专利技术提供一种极化雷达建筑损毁信息提取方法，包括：对原始PolSAR影像进行预处理包含选择研究区和提取相干矩阵；对相干矩阵进行特征值‑特征矢量分解，将研究区中小于特征阈值的像素识别为非建筑区，剔除非建筑区；利用统计模型对建筑区的相干矩阵进行建模，基于统计模型纹理参数提取纹理特征，将纹理特征大于预设的纹理阈值的像素归为倒塌建筑物，小于或等于预设的纹理阈值的像素归为完好建筑物；对建筑区进行分块，然后结合上面所得的建筑物损毁信息提取结果，评价建筑物损毁信息提取结果的精度。本发明专利技术仅用统计模型纹理参数提取倒塌建筑物，充分利用了PolSAR影像的信息，同时也提高了建筑物损毁信息提取的精度和效率。

【技术实现步骤摘要】
一种极化雷达建筑损毁信息提取方法
本专利技术涉及极化雷达(PolarimetricSyntheticApertureRadar，PolSAR)图像处理
，具体涉及一种极化雷达建筑损毁信息提取方法。
技术介绍
地震、台风等自然灾害给人们的生命财产带来极大损失。建筑物是人们生活的主要场所，灾害造成的损失大多是有建筑物倒塌造成的。因此准确的提取灾后建筑物的损毁信息能为灾害救援和灾后重建提供决策支持。灾害发生后，交通和通信设备被毁坏，天气条件恶劣，同时还有发生二次灾害的危险，使得及时获取灾区受灾情况非常困难。SAR由于其覆盖范围广、穿透性强、全天时全天候的工作能力成为获取灾害信息的主要途径。针对建筑物损毁信息提取使用的数据不同，现有的提取方法有：基于单通道或多通道SAR影像的建筑物损毁信息提取。该类方法主要利用灾害前后SAR影像的变化如相关性、相干性等，或灾后SAR影像的纹理信息提取倒塌建筑物。但单通道或多通道SAR获取的地物信息不够全面，对倒塌建筑物检测效果较差。基于PolSAR影像的建筑物损毁信息提取。PolSAR可以获得地物在4个极化通道下的回波，对地物更加敏感，为建筑物损毁信息提取提供了丰富的信息。基于多时相PolSAR影像提取倒塌建筑物主要是利用灾害前后建筑物散射机制的变化。但在某些情况下，匹配的灾害前后数据难以获取，且SAR影像的配准较困难，因此利用多时相PolSAR影像的倒塌建筑物提取受到一定的限制。利用灾后单幅PolSAR影像提取倒塌建筑物的方法主要有基于极化特征、基于纹理特征、基于极化特征和纹理特征。目前基于灾后单幅PolSAR影像提取倒塌建筑物的方法使用参数较多，过程比较复杂。且用于PolSAR影像提取倒塌建筑物的纹理信息通常都是基于总功率影像提取的，只利用了相干矩阵中的对角线元素，对PolSAR影像的信息利用不全面。
技术实现思路
有鉴于此，本专利技术的目的在于提供了一种极化雷达建筑损毁信息提取方法，以解决现有提取信息不全、参数多、过程复杂等缺陷。为实现上述目的，本专利技术提供了一种极化雷达建筑损毁信息提取方法,基于PolSAR影像处理技术,包括以下步骤:步骤101：对原始PolSAR影像进行预处理，包括选择研究区和提取相干矩阵；步骤102：对相干矩阵进行特征值-特征矢量分解，基于分解的特征值设定特征阈值，将研究区中小于特征阈值的像素识别为非建筑区，大于特征阈值的像素识别为建筑区，剔除非建筑区；步骤103：利用统计模型对建筑区的相干矩阵进行建模，基于统计模型纹理参数提取纹理特征，将纹理特征大于预设的纹理阈值的像素归为倒塌建筑物，小于或等于预设的纹理阈值的像素归为完好建筑物，得到建筑物损毁信息提取结果；步骤104：对建筑区进行分块，然后结合步骤103所得的建筑物损毁信息提取结果，评价建筑物损毁信息提取结果的精度。进一步地，步骤102中，对相干矩阵进行特征值-特征矢量分解具体为：其中，T为相干矩阵，λn和μn分别表示特征值和特征向量，λ1≥λ2≥λ3，H表示共轭的转置，特征矢量μn表示为：其中，αn表示散射矢量对应的目标散射机理，βn为目标方位角，φn、δn、γn为3个相位角，T代表转置操作，将λ2+λ3小于特征阈值的像素识别为非建筑区，λ1＞0。进一步地，步骤103中纹理特征的提取过程包括对统计模型纹理参数的估算和将估算的统计模型纹理参数取对数作为纹理特征。进一步地，统计模型为G0统计模型，其纹理参数为：M＝tr(Σ-1T)其中，T为相干矩阵，Σ＝E[T]，λ是统计模型的纹理参数，L代表视数，Var{·}表示方差，d代表相干矩阵的维度。进一步地，步骤104的具体过程如下：步骤401：设所得分块中的第j块的建筑物损毁指数为DLIj，则：其中，Cj表示第j个区块内倒塌建筑物的像素个数，Ij表示第j个区块内完好建筑物的像素个数，建筑物损毁指数DLIj代表了建筑物的损毁程度；步骤402：根据DLIj和预设分级阈值将相应区块分为不同等级，得到建筑物损毁程度的分类结果；步骤403：利用混淆矩阵，用检出率、虚警率或总精度评估建筑物信息提取的精度。为实现上述目的，本专利技术还提供了另一种极化雷达建筑损毁信息提取方法，基于PolSAR影像处理技术，包括以下步骤：步骤101：对原始PolSAR影像进行预处理，包括选择研究区和提取相干矩阵；步骤102：基于相干矩阵识别研究区中的建筑区和非建筑区，剔除非建筑区；步骤103：利用统计模型对建筑区的相干矩阵进行建模，估算统计模型的纹理参数，再将估算的纹理参数取对数作为纹理特征，然后将纹理特征大于预设的纹理阈值的像素归为倒塌建筑物，小于或等于预设的纹理阈值的像素归为完好建筑物，得到建筑物损毁信息提取结果；步骤104：对建筑区进行分块，然后结合步骤103所得建筑物损毁信息提取结果，评价建筑物损毁信息提取结果的精度。进一步地，步骤103中，统计模型为G0统计模型，其纹理参数的估算基于相干矩阵：M＝tr(Σ-1T)其中，T为相干矩阵，Σ＝E[T]，λ是统计模型的纹理参数，L代表视数，Var{·}表示方差，d代表相干矩阵的维度。本专利技术与现有技术相比，其有益效果是：用统计模型纹理参数反映建筑区的均匀程度，根据倒塌建筑区的均匀程度较完好建筑区高，利用基于统计模型纹理参数提取的纹理特征实现倒塌建筑和完好建筑的提取。针对倒塌建筑区和道路、水体等非建筑区都具有较高的均匀程度，利用特征值-特征矢量分解的特征值剔除研究区内的非建筑，减少非建筑物与倒塌建筑物的混淆，从而提高建筑物损毁信息提取的精度。这种基于PolSAR影像统计模型纹理参数的极化雷达建筑损毁信息提取方法充分利用PolSAR影像的信息，且减少了对数据的依赖，效率和评估精度都有所提高。附图说明图1是本专利技术极化雷达建筑损毁信息提取方法实施例一的步骤示意图；图2是本专利技术极化雷达建筑损毁信息提取方法的原理示意图。具体实施方式为使本专利技术的目的、技术方案和优点更加清楚，下面将结合附图对本专利技术实施方式作进一步地描述。本专利技术技术方案可采用计算机技术实现自动运行，请参阅图1和图2所示。本专利技术提供了一种极化雷达建筑损毁信息提取方法,基于PolSAR影像处理技术,包括以下步骤：步骤101：对原始PolSAR影像进行预处理，包括选择研究区和提取相干矩阵。在提取建筑物损毁信息之前，需要对原始的PolSAR影像进行预处理，包括选取研究区域和提取相干矩阵。由于PolSAR影像存在大量的相干斑，为了降低噪声的影响，一般需要对原始PolSAR影像进行滤波处理，如精细Lee滤波、SigmaLee滤波等。当PolSAR分辨率较低时，为了保留影像信息也可以不做滤波处理。步骤102：对相干矩阵进行特征值-特征矢量分解，基于分解的特征值设定特征阈值，将研究区中小于特征阈值的像素识别为非建筑区，大于特征阈值的像素识别为建筑区，剔除非建筑区。相干矩阵进行特征值-特征矢量分解具体为：其中，T为相干矩阵，λn和μn分别表示特征值和特征向量，λ1≥λ2≥λ3，H表示共轭的转置，特征矢量μn表示为：其中，αn表示散射矢量对应的目标散射机理，βn为目标方位角，φn、δn、γn为3个相位角，T代表转置操作，将λ2+λ3小于特征阈值的像素识别为非建筑区，将非建筑区剔除。由于本专利技术是根据倒塌建本文档来自技高网...

【技术保护点】
一种极化雷达建筑损毁信息提取方法，基于极化PolSAR影像处理技术,其特征在于，包括以下步骤：步骤101：对原始PolSAR影像进行预处理，包括选择研究区和提取相干矩阵；步骤102：对相干矩阵进行特征值‑特征矢量分解，基于分解的特征值设定特征阈值，将研究区中小于特征阈值的像素识别为非建筑区，大于特征阈值的像素识别为建筑区，剔除非建筑区；步骤103：利用统计模型对建筑区的相干矩阵进行建模，基于统计模型纹理参数提取纹理特征，将纹理特征大于预设的纹理阈值的像素归为倒塌建筑物，小于或等于预设的纹理阈值的像素归为完好建筑物，得到建筑物损毁信息提取结果；步骤104：对建筑区进行分块，然后结合步骤103所得的建筑物损毁信息提取结果，评价建筑物损毁信息提取结果的精度。

【技术特征摘要】
1.一种极化雷达建筑损毁信息提取方法，基于极化PolSAR影像处理技术,其特征在于，包括以下步骤：步骤101：对原始PolSAR影像进行预处理，包括选择研究区和提取相干矩阵；步骤102：对相干矩阵进行特征值-特征矢量分解，基于分解的特征值设定特征阈值，将研究区中小于特征阈值的像素识别为非建筑区，大于特征阈值的像素识别为建筑区，剔除非建筑区；步骤103：利用统计模型对建筑区的相干矩阵进行建模，基于统计模型纹理参数提取纹理特征，将纹理特征大于预设的纹理阈值的像素归为倒塌建筑物，小于或等于预设的纹理阈值的像素归为完好建筑物，得到建筑物损毁信息提取结果；步骤104：对建筑区进行分块，然后结合步骤103所得的建筑物损毁信息提取结果，评价建筑物损毁信息提取结果的精度。2.如权利要求1所述的极化雷达建筑损毁信息提取方法，其特征在于，步骤102中，对相干矩阵进行特征值-特征矢量分解具体为：其中，T为相干矩阵，λn和μn分别表示特征值和特征向量，λ1≥λ2≥λ3，H表示共轭的转置，λ1＞0,特征矢量μn表示为：其中，αn表示散射矢量对应的目标散射机理，βn为目标方位角，φn、δn、γn为3个相位角，T代表转置操作，将λ2+λ3小于特征阈值的像素识别为非建筑区。3.根据权利要求1所述的极化雷达建筑损毁信息提取方法，其特征在于，步骤103中纹理特征的提取过程包括对统计模型纹理参数的估算和将估算的统计模型纹理参数取对数作为纹理特征。4.根据权利要求3所述的极化雷达建筑损毁信息提取方法，其特征在于：统计模型为G0统...

【专利技术属性】
技术研发人员：陈启浩，刘修国，李林林，
申请(专利权)人：中国地质大学武汉，
类型：发明
国别省市：湖北,42