一种基于全极化SAR新特征的海洋溢油检测方法技术

本发明专利技术公开了一种基于全极化SAR新特征的海洋溢油检测方法，包括以下步骤：极化相干矩阵的提取；特征值和特征向量分解；构建基尼系数新特征；阈值分割；本发明专利技术方法基于特征值分解理论提出一种新的极化特征，称为极化基尼系数；该特征不仅能够反映集合中不同目标之间的极化状态，还能够描述不同散射类型在统计意义上的不纯度，当不纯度为0时表明对应集合中的散射机制类别一致；新特征通过目标分解理论得到不同散射机制类别，通过计算该像素的主导散射机制在统计意义上的不纯度提取溢油信息，通过对全极化SAR溢油影像利用该方法提取特征，解决了不能区分生物油膜与矿物膜的问题，且新的特征提取技术具有抑制相干斑噪声的能力。

A new method of oil spill detection based on new features of full polarimetric SAR

【技术实现步骤摘要】
一种基于全极化SAR新特征的海洋溢油检测方法
本专利技术涉及海洋溢油检测
，尤其涉及一种基于全极化SAR新特征的海洋溢油检测方法。
技术介绍
海上溢油是造成海洋环境污染损害的主要因素之一，由管道破裂，油轮碰撞和钻井平台爆炸引起的海洋石油污染是全球关注的问题，自1970年以来，全世界每年发生的石油污染超过10000吨，因此如何有效地检测海面溢油就显得尤为重要。海上溢油发生后，通常要了解发生的位置、溢油量和扩散趋势，在已投入的监测系统中，卫星遥感是最重要和最有效的手段之一，在海上溢油监测和应急响应中发挥着越来越重要的作用。目前，最常用的卫星遥感传感器包括可见光、红外、紫外、激光荧光和微波传感器等，其中又以合成孔径雷达(SAR)能够穿云透雾，且不受昼夜影响的特点而被广泛应用；浮油通常表现为SAR图像中的暗斑，这是由于SAR海面成像的主要机理造成的。当海面上有油膜存在时，海面上的毛细波和短重力波会不断地受到阻尼作用，降低海面粗糙度，从而削弱雷达后向散射回波信号。一般来说，传统的单极化溢油检测可以分为以下三个步骤：(1)暗斑检测：例如，感兴趣区域(ROI)的识别；(2)特征提取：为了进一步区分暗斑的组成，有必要在暗斑检测之后提取图像的各种特征，挑选合适的特征用于随后的分类任务；(3)分类：减少油膜和类油膜的错分现象；随着SAR的极化方式由单极化向双极化、全极化的发展，基于极化SAR数据的新的溢油检测研究已经进行，相比于单极化，全极化SAR图像不仅含有强度信息，而且还包含相位信息，能够有效地获取海面目标的全极化散射特性，更全面地反映海面目标的几何形状和物理特性；在现有的极化特征提取技术中，大都能够很好的区分油膜与海水，但对于全极化SAR溢油影像来说，影像并不是只有油膜与海水，还包含有低风速区、雨团、生物油膜等类油膜现象，特别是生物油膜与矿物油膜的区分在近几年来越来越被重视，而现有的特征提取技术并不能够很好的区分生物油膜与矿物油膜，且容易受到海面相干斑噪声的影响，因此，本专利技术提出一种基于全极化SAR新特征的海洋溢油检测方法，以解决现有技术中的不足之处。
技术实现思路
针对上述问题，本专利技术提出一种基于全极化SAR新特征的海洋溢油检测方法，基于特征值分解理论提出一种新的极化特征，称为极化基尼系数，该特征不仅能够反映集合中不同目标之间的极化状态，还能够描述不同散射类型在统计意义上的不纯度，新特征通过目标分解理论得到不同散射机制类别，通过计算该像素的主导散射机制在统计意义上的不纯度提取溢油信息，解决了不能区分生物油膜与矿物膜的问题，且新的特征提取技术具有抑制相干斑噪声的能力。本专利技术提出一种基于全极化SAR新特征的海洋溢油检测方法，包括以下步骤：步骤一：提取极化相干矩阵，包括原始散射矩阵S与相干矩阵T3的提取，极化散射矩阵S能完整描述SAR目标的电磁散射特性，其定义如公式(1)所示，相干矩阵T3由目标矢量k与自身共轭转置矢量的外积生成，相干矩阵T3定义如公式(2)所示；公式(2)中，*T，*分别代表共轭转置和共轭，<·>表示整体平均值；步骤二：特征值和特征向量分解。基于目标分解理论，将3×3埃尔米特平均相干矩阵T3计算得到的特征值和特征向量生成相干矩阵的一种对角化形式，相干矩阵T3的定义如公式(3)所示：其中，∑是3×3的非负实对角阵，特征值(λ1≥λ2≥λ3≥0)；U3＝[e1e2e3]是3×3特殊酉矩阵SU(3)，其中e1,e2,e3分别是相干矩阵T3的特征矢量，然后相干矩阵T3被重新定义，新的定义如公式(4)所示；步骤三：构建基尼系数新特征。基于步骤二中的特征值与特征向量分解，定义极化基尼系数如公式(5)所示：其中，pi对应于特征值λi获得的伪概率，满足p1+p2+p3＝1，λ1≥λ2≥λ3≥0，且由于特征值为旋转不变的，极化基尼系数也具有旋转不变性；步骤四：阈值分割。通过目标分解理论得到不同散射机制类别以及对应的散射机制，根据基尼系数新特征判断不同散射机制在统计意义上的不纯度与有序度进行溢油信息的提取。进一步改进在于：所述公式(1)中SHV代表水平极化发射，垂直极化接收；SVH代表垂直极化发射，水平极化接收；SHH代表水平极化发射与接收；SVV代表垂直极化发射与接收。进一步改进在于：所述公式(1)中当满足互易性定理时，SHV＝SVH，极化散射矩阵S用Pauli基分解方法得到目标散射矢量k，则极化散射矩阵S的定义如公式(6)所示：其中，因子用于保证目标矢量的范数与散射总功率相等。进一步改进在于：所述步骤一中散射矩阵S的对角元素代表相同极化方式的入射场和散射场关系，用“同极化”项表示；非对角元素代表正交极化方式的入射角和散射场关系，用“交叉极化”项表示。进一步改进在于：所述步骤二中求解特征矢量e1,e2,e3，则获得3个互不相关的目标用于构建统计模型，将相干矩阵T3展开为3个相互独立的目标之和，每个目标对应一种确定的散射机制，然后由一个等价的简单散射矩阵表示，再确定的散射机制成分i在整个散射过程中所占的权重，由特征值λi描述，散射机制的类型与归一化特征矢量ei相关。进一步改进在于：所述步骤三中公式(5)中，当特征值只有一个非零时(λ1≠0,λ2＝λ3＝0)，则相干矩阵T3就退化为由极化散射矩阵S定义的点目标的散射情况。进一步改进在于：所述步骤三中公式(5)中，当所有的特征值均不为零且均相等(λ1＝λ2＝λ3≠0)，则相干矩阵T3就代表一种完全去极化的随机散射状态，介于这两种极端情况之间的目标是部分极化的，此时相干矩阵T3的非零特征值的数量大于1且彼此不完全相等。进一步改进在于：所述步骤四中当不纯度越低，则相应散射机制的有序度就越高，极化基尼系数也越小；当不纯度越高，则相应散射机制的有序度就越低，极化基尼系数就越大；当不纯度达到最小值时，表明对应集合中的散射机制类别一致；当不纯度达到最大值时，表明集合中的目标散射呈现随机噪声状态。进一步改进在于：所述步骤四中从极化状态角度，当极化基尼系数的值较低时，则系统可以被认为是弱去极化的，主导散射机制被视为特定等效点的目标散射机制，根据该机制则进行选择对应于最大特征值的特征向量，而忽略其他特征向量；当极化基尼系数的值较高时，集合内的平均散射体将呈现去极化状态，不再存在单个散射目标，则需要考虑来自整个特征值分布谱的所有可能点目标散射类型的混合比例；当极化基尼系数进一步增大时，从极化测量数据中识别的散射机制的数目也将逐渐减少，当极化基尼系数达到最大值时，极化信息将为零，则此时目标散射完全是随机噪声过程。本专利技术的有益效果为：本专利技术方法基于特征值分解理论提出一种新的极化特征，称为极化基尼系数，该特征不仅能够反映集合中不同目标之间的极化状态，还能够描述不同散射类型在统计意义上的不纯度，当不纯度为0时表明对应集合中的散射机制类别一致，新特征通过目标分解理论得到不同散射机制类别，通过计算该像本文档来自技高网...

【技术保护点】
1.一种基于全极化SAR新特征的海洋溢油检测方法，其特征在于，包括以下步骤：/n步骤一：提取极化相干矩阵，包括原始散射矩阵S与相干矩阵T

【技术特征摘要】
1.一种基于全极化SAR新特征的海洋溢油检测方法，其特征在于，包括以下步骤：
步骤一：提取极化相干矩阵，包括原始散射矩阵S与相干矩阵T3的提取。极化散射矩阵S能完整描述SAR目标的电磁散射特性，其定义如公式(1)所示，相干矩阵T3由目标矢量k与自身共轭转置矢量的外积生成，相干矩阵T3定义如公式(2)所示；






公式(2)中，*T，*分别代表共轭转置和共轭，<·>表示整体平均值；
步骤二：特征值和特征向量分解。基于目标分解理论，将3×3埃尔米特平均相干矩阵T3计算得到的特征值和特征向量生成相干矩阵的一种对角化形式，相干矩阵T3的定义如公式(3)所示：



其中，∑是3×3的非负实对角阵，U3是3×3特殊酉矩阵SU(3)，如果λ1,λ2,λ3分别代表相干矩阵T3中的三个非负特征值，且λ1≥λ2≥λ3≥0；e1,e2,e3分别代表相干矩阵T3的特征矢量，则相干矩阵T3可以进一步展开如公式(4)所示：



步骤三：构建基尼系数新特征。基于步骤二中的特征值与特征向量分解，定义极化基尼系数如公式(5)所示：



其中，pi对应于特征值λi获得的伪概率，满足p1+p2+p3＝1，λ1≥λ2≥λ3≥0，且由于特征值为旋转不变的，极化基尼系数也具有旋转不变性；
步骤四：阈值分割。通过目标分解理论得到不同散射机制类别以及对应的散射机制，根据基尼系数新特征判断不同散射机制在统计意义上的不纯度与有序度进行溢油信息的提取。


2.根据权利要求1所述的一种基于全极化SAR新特征的海洋溢油检测方法，其特征在于：所述公式(1)中SHV代表水平极化发射，垂直极化接收；SVH代表垂直极化发射，水平极化接收；SHH代表水平极化发射与接收；SVV代表垂直极化发射与接收。


3.根据权利要求2所述的一种基于全极化SAR新特征的海洋溢油检测方法，其特征在于：所述公式(1)中当满足互易性定理时，SHV＝SVH，极化散射矩阵S用Pauli基分解方法得到目标散射矢量k，则极化散射矩阵S的定义如公式(6)所示：



其中，因子用于保证目标矢量的范数与散射总功率相等。


4.根据权利要求1所述的一种基于全极化SAR新特征的海洋溢油检测方法，其特征在于：所述步骤一中散射矩阵S的...

【专利技术属性】
技术研发人员：宋冬梅，任慧敏，王斌，甄宗晋，崔建勇，许明明，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：山东;37