CN104318572A - 基于极化方向角补偿的polsar图像海上目标检测方法 - Google Patents

Abstract

本发明公开了一种基于极化方向角补偿的POLSAR图像海上目标检测方法，该方法利用海上人造金属目标的极化散射特性进行检测，针对海浪扰动产生的方向角偏移对海杂波反射对称性影响的问题，通过对POLSAR影像进行方向角补偿，来提高海上目标检测精度。首先获取全极化影像，进行相关预处理操作；然后通过相干矩阵计算出极化方向角，并对相干矩阵进行方向角补偿；最后提取出方向角补偿后的相干矩阵中T′₁₃元素的强度值，进行海上人造金属目标检测。该方法克服了常规海上目标检测算法理论复杂、难于实现及稳健性差等缺陷，同时从目标的极化散射特性出发，较准确的从海杂波、目标旁瓣中分离出海上目标。该检测方法原理直观，算法简洁，便于编程实现和扩展。

Description

基于极化方向角补偿的POLSAR图像海上目标检测方法

技术领域

本发明涉及POLSAR图像领域，特别涉及一种基于极化方向角补偿的POLSAR图像海上目标检测方法。

背景技术

海洋这一蓝色国土，拥有丰富的各类资源，在国民经济发展、环境气候变化以及军事国防中有着极其重要的战略意义。

近年来，有关海洋问题频发不断，预防和解决这些问题就迫切需求对海洋进行监测。海上人造金属目标(包括舰船、钻井平台等，以下简称海上目标)与人类生活生产有着紧密关联，因而成为海洋监测的重点。

目前，常规的海上目标监测方法包括岸基、舰载探测装置，光学遥感等。岸基、舰载探测装置存在着成本高、周期长及监测范围小等缺陷，在海上目标监测应用中受到了极大的限制。对于光学遥感，由于其受天气条件影响十分严重，而海洋中天气又经常恶劣多变，致使光学遥感很难做到对海上目标进行实时动态的监测。因此，如何实时、精确、高效的对海上目标进行监测对准确决策、高效管理有着重大意义。

极化合成孔径雷达(Polarimetric SAR，简称POLSAR)技术，集成了极化散射机理、后向散射介质信息及对目标物理特性敏感为一体的优点，具有全天时、全天候、观测周期短，范围大，精度高的特点，在海上目标检测领域有着巨大潜力。

目前基于极化SAR图像的海上目标检测方法主要有：1)相干目标分解(Coherent targetdecomposition,CTD)，其通过引入对称散射特征方法(Symmetric Scattering CharacterizationMethod，简称SSCM)，较好的从相干散射结构框架中提取对称散射成分信息来进行海上目标检测；2)极化交叉熵(Polarimetric cross entropy,PCE)，它是基于海杂波和海上目标所呈现出不同的极化熵值来进行海上目标检测的；3)反射对称性(XC)，其基于海杂波和海上目标所呈现出不同的反射对称性来进行目标检测。在极化SAR影像中，反射对称性是通过同极化通道与交叉极化通道间的相关性来呈现出的。经过大量实验验证，基于反射对称性方法(XC)已经被证明是比较简单优越的检测方法，它适用于X、C和L等不同波段的极化数据，增强了目标与海面之间的差异，同时不受地物的散射机理、风向风况及大范围的入射角的影响，并且检测结果明显优于PCE、CTD的检测结果。为此，该方法是目前检测精确度较高的方法，其理论清晰明了，流程简洁高效。但是该方法忽略了受海浪扰动导致的方向角变化对海杂波反射对称性的影响问题，从而一定程度上影响了目标检测效果。

发明内容

本发明的目的在于提出一种基于极化方向角补偿的POLSAR图像海上目标检测方法。利用海杂波满足反射对称性，而海上目标满足反射非对称性的特点来检测海上目标。同时通过极化方向角补偿解决由于海浪扰动产生的方向角偏移对海杂波反射对称性影响的问题。

一种基于极化方向角补偿的POLSAR图像海上目标检测方法，包括如下步骤：

步骤1：获取海上待检测场景全极化影像数据即POLSAR数据；

所述全极化影像数据包括HH、HV、VH和VV四种极化通道数据；

步骤2：对POLSAR数据进行预处理操作；

所述预处理操作包括辐射校正、几何校正及陆地掩膜操作；

【采用不同的传感器获取POLSAR数据时，对POLSAR数据的预处理方式会不一样，但是所有的预处理方式均属于本领域普通技术人员熟知的技术；】

步骤3：从预处理后的海上待检测场景对应的POLSAR数据中依次提取每个像素对应的S矩阵；

$S=\left[\begin{array}{cc}{S}_{\mathrm{hh}}& S_{\mathrm{hv}}\\ S_{\mathrm{vh}}& S_{\mathrm{vv}}\end{array}\right]=\left[\begin{array}{cc}{S}_{\mathrm{hh}}& S_{\mathrm{hv}}\\ S_{\mathrm{hv}}& S_{\mathrm{vv}}\end{array}\right]$

其中，所述S矩阵满足互易对称定律：S_hv＝S_vh，S_hh表示水平方式发射水平方式接收的极化通道数据，S_hv表示水平方式发射垂直方式接收的极化通道数据，S_vh表示垂直方式发射水平方式接收的极化通道数据，S_vv表示垂直方式发射垂直方式接收的极化通道数据；

步骤4：通过Pauli基进行共轭相乘，将S矩阵转换为相干矩阵T：

$T=<k\&CenterDot;k^H>=\left[\begin{array}{ccc}{T}_{11}& T_{12}& T_{13}\\ T_{21}& T_{22}& T_{23}\\ T_{31}& T_{32}& T_{33}\end{array}\right]$

其中， $k=\frac{1}{\sqrt{2}}\left[\begin{array}{c}{S}_{\mathrm{hh}}+S_{\mathrm{vv}}\\ S_{\mathrm{hh}}-S_{\mathrm{vv}}\\ {2S}_{\mathrm{hv}}\end{array}\right],$

$\left[\begin{array}{ccc}{T}_{11}& T_{12}& T_{13}\\ T_{21}& T_{22}& T_{23}\\ T_{31}& T_{32}& T_{33}\end{array}\right]=\frac{1}{2}\left[\begin{array}{ccc}<{|S_{\mathrm{hh}}+S_{\mathrm{vv}}|}^2>& <(S_{\mathrm{hh}}+S_{\mathrm{vv}}){(S_{\mathrm{hh}}-S_{\mathrm{vv}})}^{*}>& 2<(S_{\mathrm{hh}}+S_{\mathrm{vv}})S_{\mathrm{hv}}^{*}>\\ <(S_{\mathrm{hh}}-S_{\mathrm{vv}}){(S_{\mathrm{hh}}+S_{\mathrm{vv}})}^{*}>& <{|S_{\mathrm{hh}}-S_{\mathrm{vv}}|}^2>& 2<(S_{\mathrm{hh}}-S_{\mathrm{vv}})S_{\mathrm{hv}}^{*}>\\ 2<S_{\mathrm{hv}}{(S_{\mathrm{hh}}+S_{\mathrm{vv}})}^{*}>& 2\left(S_{\mathrm{hv}}{(S_{\mathrm{hh}}-S_{\mathrm{vv}})}^{*}\right)& 4<{\left|S_{\mathrm{hv}}\right|}^2>\end{array}\right]$

其中，上标H表示共轭转置，*表示复共轭，<·>表示邻域平均；

步骤5：采用精致Lee滤波，对相干矩阵T进行滤波处理，得到滤波后的相干矩阵T⁰；

$T^0=\left[\begin{array}{ccc}{T}_{11}^0& T_{12}^0& T_{13}^0\\ T_{21}^0& T_{22}^0& T_{23}^0\\ T_{31}^0& T_{32}^0& T_{33}^0\end{array}\right]$

【消除斑点噪声对影像的影响】；

步骤6：利用滤波后的相干矩阵T⁰计算极化方向角：

$\mathrm{\&theta;}=\frac{1}{4}\arctan \left(\frac{2\&CenterDot;\mathrm{Re}\left(T_{23}^0\right)}{T_{22}^0-T_{33}^0}\right)$

其中， $\mathrm{\&theta;}=\frac{1}{4}\arctan \left(\frac{2\&CenterDot;\mathrm{Re}\left(T_{23}^0\right)}{T_{22}^0-T_{33}^0}\right)+\frac{1}{4}\mathrm{\&pi;},$ 若 $\mathrm{\&theta;}>\frac{1}{4}\mathrm{\&pi;},$ 则 $\mathrm{\&theta;}=\mathrm{\&theta;}-\frac{1}{2}\mathrm{\&pi;};$ 表示的实部；θ∈(-π/4,π/4]；

【极化方向角θ为将入射平面围绕雷达视线旋转直到与表面法线重合所旋转的角度；】

步骤7：利用步骤6中求取出的极化方向角对滤波后的相干矩阵T矩阵进行方向角补偿，得到补偿方向角后的相干矩阵T₀：

$T_0={\mathrm{QT}}^{0}Q=\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos 2\mathrm{\&theta;}& \sin 2\mathrm{\&theta;}\\ 0& -\sin 2\mathrm{\&theta;}& \cos 2\mathrm{\&theta;}\end{array}\right]T^0\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos 2\mathrm{\&theta;}& -\sin 2\mathrm{\&theta;}\\ 0& \sin 2\mathrm{\&theta;}& \cos 2\mathrm{\&theta;}\end{array}\right]=\left[\begin{array}{ccc}{T}_{11}^{\&prime;}& T_{12}^{\&prime;}& T_{13}^{\&prime;}\\ T_{21}^{\&prime;}& T_{22}^{\&prime;}& T_{23}^{\&prime;}\\ T_{31}^{\&prime;}& T_{32}^{\&prime;}& T_{33}^{\&prime;}\end{array}\right]$

其中，T₀表示补偿方向角后的相干矩阵，Q为旋转矩阵；

【极化方向角补偿即为将入射平面围绕雷达视线旋转直到与表面法线重合，从而使得目标同极化响应达到最大；】

步骤8：提取补偿方向角后的相干矩阵T₀中的元素T′₁₃，计算T′₁₃的强度值I；

【 $T_{13}^{\&prime;}=(S_{\mathrm{hh}1}+S_{\mathrm{vv}1})S_{\mathrm{hv}1}^{*}=S_{\mathrm{hh}1}{S}_{\mathrm{hv}1}^{*}+S_{\mathrm{vv}1}{S}_{\mathrm{hv}1}^{*},S_{\mathrm{hh}1}{S}_{\mathrm{hv}1}^{*}$ 和均表示滤波及方向角补偿后的同极化通道与交叉极化通道之间的相关性；】

【由于目标的反射非对称性反映到极化SAR影像中表现为同极化通道与交叉极化通道之间的相关性。海杂波通常满足反射对称性特点，而舰船及钻井平台等人造金属目标通常满足反射非对称性特点，因此可以利用T′₁₃的强度值I来进行海上目标检测。】

步骤9：采用恒虚警率检测器求取每个像素的检测阈值N；

步骤10：遍历海上待检测场景中的每一个像素，计算每个像素对应的T′₁₃的强度值I，比较I与阈值N，若I大于阈值N，则当前像素为海上目标，否则为杂波。

【由于海浪扰动产生的方向角偏移影响着海杂波的反射对称性，因此利用海杂波满足反射对称性的特点检测海上目标的效果同样受其影响。通过极化方向角补偿来消除该现象的影响。】

【人造金属目标的反射非对称性在极化SAR影像中是通过同极化通道与交叉极化通道之间的相关性非零来呈现的，即

元素包含了极化SAR影像中同极化通道与交叉极化通道之间的相关性的所有存在情况，并进行了等权组合。相比单独使用或T′₁₃能够更好的描述目标的反射非对称性，适合于海上目标的精确检测。】

在所述步骤9中采用双参数恒虚警率检测器进行局部CFAR检测获得每个像素的比较阈值N，局部CFAR检测中所使用的滑动窗口以待检测像素为中心的环形窗口，包括目标窗口、保护窗口、背景窗口，其中，所述目标窗口大小为能包含最小检测目标所占像素个数大小，所述保护窗口为能包含最大检测目标所占像素个数大小，所述背景窗口大小为保护窗口大小的1.5-2.5倍；

【最大最小目标是根据实际现实情况来判断的，如最大舰船目标的长轴不会超过500米；最大目标的像素不会超出保护窗口，背景窗口是用于海杂波统计，越大统计越详细，但是太大就有可能其他目标落入背景窗口中，影响统计效果，所以一般约为保护窗口的两倍；】

N＝u_b+σ_b·n

其中，u_b表示背景窗口的灰度均值，σ_b表示背景窗口的灰度标准差，其中，n表示检测器设计参数，由公式计算获得；Pfa为给定的虚警率，小于10^-8的正数。

【双参数恒虚警率检测器的主旨是寻找出比周围区域异常明亮的像素，基于杂波服从高斯统计模型及给定的虚警概率来求出阈值，进行目标检测。】

有益效果

本发明提供了一种基于极化方向角补偿的POLSAR图像海上目标的检测方法，将极化中的方向角补偿与目标反射非对称特性相结合，利用补偿方向角后的相干矩阵中的T′₁₃元素来描述反射非对称性，消除海浪扰动导致的方向角变化影响海杂波反射对称性的问题，更准确地描述了目标的反射非对称性，提高了极化SAR图像海上目标检测的效果。其原理直观，算法简洁，稳定性强，易于实现，存在着实用价值。该方法充分利用考虑海杂波及目标的极化散射特性，有效的消除目标方位模糊和旁瓣的影响，高效地检测出海上目标，并同时减少检测虚警率。

附图说明

图1为本发明的处理流程图；

图2为极化方向角的示意图；

图3为环形窗口的示意图；

图4为原始全极化影像灰度图和检测结果图，其中，(a)为原始全极化影像灰度图，(b)为检测结果图，图中矩形框为检测目标，椭圆框为虚假检测，三角形为漏检目标；

图5为极化方向角补偿前后比较图，其中，(a)为从海上场景全极化影像中提取的局部区域灰度图，(b)为对(a)未进行方向角补偿时得到的检测结果图，(c)为对(a)进行方向角补偿的检测结果图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明了，以下结合附图及实施事例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施事例仅仅用以解释本发明，并不用于限定本发明。

对于获取的全极化SAR影像，用极化散射矩阵(S矩阵)来表达四个通道影像数据，如下式所示。该式中S_pq表示q极化方式发射、p极化方式接收时目标的后向散射，p、q为水平h或垂直v极化方式。一般而言，对互易介质，根据单站互易原理，S矩阵将满足S_hv＝S_vh的条件，此时散射矩阵为对称矩阵。

$S=\left[\begin{array}{cc}{S}_{\mathrm{hh}}& S_{\mathrm{hv}}\\ S_{\mathrm{vh}}& S_{\mathrm{vv}}\end{array}\right]=\left[\begin{array}{cc}{S}_{\mathrm{hh}}& S_{\mathrm{hv}}\\ S_{\mathrm{hv}}& S_{\mathrm{vv}}\end{array}\right]$

如图1所示，为本发明所述方法的流程示意图，一种基于极化方向角补偿的POLSAR图像海上目标检测方法，包括如下步骤：

步骤1：获取海上待检测场景全极化影像数据即POLSAR数据；

所述全极化影像数据包括HH、HV、VH和VV四种极化通道数据；

步骤2：对POLSAR数据进行预处理操作；

所述预处理操作主要包括辐射校正、几何校正及陆地掩膜等操作；

【采用不同的传感器获取POLSAR数据时，对POLSAR数据的预处理方式会不一样，但是所有的预处理方式均属于本领域普通技术人员熟知的技术；】

步骤3：从预处理后的海上待检测场景对应的POLSAR数据中依次提取每个像素对应的S矩阵；

$S=\left[\begin{array}{cc}{S}_{\mathrm{hh}}& S_{\mathrm{hv}}\\ S_{\mathrm{vh}}& S_{\mathrm{vv}}\end{array}\right]=\left[\begin{array}{cc}{S}_{\mathrm{hh}}& S_{\mathrm{hv}}\\ S_{\mathrm{hv}}& S_{\mathrm{vv}}\end{array}\right]$

其中，所述S矩阵满足互易对称定律：S_hv＝S_vh，S_hh表示水平方式发射水平方式接收的极化通道数据，S_hv表示水平方式发射垂直方式接收的极化通道数据，S_vh表示垂直方式发射水平方式接收的极化通道数据，S_vv表示垂直方式发射垂直方式接收的极化通道数据；

步骤4：通过Pauli基进行共轭相乘，将S矩阵转换为相干矩阵T：

$T=<k\&CenterDot;k^H>=\left[\begin{array}{ccc}{T}_{11}& T_{12}& T_{13}\\ T_{21}& T_{22}& T_{23}\\ T_{31}& T_{32}& T_{33}\end{array}\right]$

其中， $k=\frac{1}{\sqrt{2}}\left[\begin{array}{c}{S}_{\mathrm{hh}}+S_{\mathrm{vv}}\\ S_{\mathrm{hh}}-S_{\mathrm{vv}}\\ {2S}_{\mathrm{hv}}\end{array}\right],$

$\left[\begin{array}{ccc}{T}_{11}& T_{12}& T_{13}\\ T_{21}& T_{22}& T_{23}\\ T_{31}& T_{32}& T_{33}\end{array}\right]=\frac{1}{2}\left[\begin{array}{ccc}<{|S_{\mathrm{hh}}+S_{\mathrm{vv}}|}^2>& <(S_{\mathrm{hh}}+S_{\mathrm{vv}}){(S_{\mathrm{hh}}-S_{\mathrm{vv}})}^{*}>& 2<(S_{\mathrm{hh}}+S_{\mathrm{vv}})S_{\mathrm{hv}}^{*}>\\ <(S_{\mathrm{hh}}-S_{\mathrm{vv}}){(S_{\mathrm{hh}}+S_{\mathrm{vv}})}^{*}>& <{|S_{\mathrm{hh}}-S_{\mathrm{vv}}|}^2>& 2<(S_{\mathrm{hh}}-S_{\mathrm{vv}})S_{\mathrm{hv}}^{*}>\\ 2<S_{\mathrm{hv}}{(S_{\mathrm{hh}}+S_{\mathrm{vv}})}^{*}>& 2\left(S_{\mathrm{hv}}{(S_{\mathrm{hh}}-S_{\mathrm{vv}})}^{*}\right)& 4<{\left|S_{\mathrm{hv}}\right|}^2>\end{array}\right]$

其中，上标H表示共轭转置，*表示复共轭，<·>表示邻域平均；

步骤5：采用精致Lee滤波，对相干矩阵T进行滤波处理，得到滤波后的相干矩阵T⁰；

$T^0=\left[\begin{array}{ccc}{T}_{11}^0& T_{12}^0& T_{13}^0\\ T_{21}^0& T_{22}^0& T_{23}^0\\ T_{31}^0& T_{32}^0& T_{33}^0\end{array}\right]$

【精致Lee滤波采用了边沿对齐窗口和最小均方误差(MMSE,Minimum Mean SquareError)方法,能够避免交叉串扰效应，保持好极化信息和影像的清晰度，消除斑点噪声对影像的影响】；

步骤6：利用滤波后的相干矩阵T⁰计算极化方向角：

$\mathrm{\&theta;}=\frac{1}{4}\arctan \left(\frac{2\&CenterDot;\mathrm{Re}\left(T_{23}^0\right)}{T_{22}^0-T_{33}^0}\right)$

其中， $\mathrm{\&theta;}=\frac{1}{4}\arctan \left(\frac{2\&CenterDot;\mathrm{Re}\left(T_{23}^0\right)}{T_{22}^0-T_{33}^0}\right)+\frac{1}{4}\mathrm{\&pi;},$ 若 $\mathrm{\&theta;}>\frac{1}{4}\mathrm{\&pi;},$ 则 $\mathrm{\&theta;}=\mathrm{\&theta;}-\frac{1}{2}\mathrm{\&pi;};$ 表示的实部；θ∈(-π/4,π/4]；

【如图2所示，极化方向角θ为将入射平面围绕雷达视线旋转直到与表面法线重合所旋转的角度；】

由于海浪扰动产生的方向角偏移对海杂波的反射对称性影响问题，这将影响采用反射非对称性进行目标检测的检测效果。通过方向角补偿消除该影响，从而改善目标检测效果。

步骤7：利用步骤6中求取出的极化方向角对滤波后的相干矩阵T矩阵进行方向角补偿，得到补偿方向角后的相干矩阵T₀：

$T_0={\mathrm{QT}}^{0}Q=\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos 2\mathrm{\&theta;}& \sin 2\mathrm{\&theta;}\\ 0& -\sin 2\mathrm{\&theta;}& \cos 2\mathrm{\&theta;}\end{array}\right]T^0\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos 2\mathrm{\&theta;}& -\sin 2\mathrm{\&theta;}\\ 0& \sin 2\mathrm{\&theta;}& \cos 2\mathrm{\&theta;}\end{array}\right]=\left[\begin{array}{ccc}{T}_{11}^{\&prime;}& T_{12}^{\&prime;}& T_{13}^{\&prime;}\\ T_{21}^{\&prime;}& T_{22}^{\&prime;}& T_{23}^{\&prime;}\\ T_{31}^{\&prime;}& T_{32}^{\&prime;}& T_{33}^{\&prime;}\end{array}\right]$

其中，T₀表示补偿方向角后的相干矩阵，Q为旋转矩阵；

【极化方向角补偿即为将入射平面围绕雷达视线旋转直到与表面法线重合，使得海面小波面法线与雷达波入射平面在同一平面内，从而使得目标同极化响应达到最大；】

步骤8：提取补偿方向角后的相干矩阵T₀中的元素T′₁₃，计算T′₁₃的强度值I；

【 $T_{13}^{\&prime;}=(S_{\mathrm{hh}1}+S_{\mathrm{vv}1})S_{\mathrm{hv}1}^{*}=S_{\mathrm{hh}1}{S}_{\mathrm{hv}1}^{*}+S_{\mathrm{vv}1}{S}_{\mathrm{hv}1}^{*},S_{\mathrm{hh}1}{S}_{\mathrm{hv}1}^{*}$ 和均表示滤波及方向角补偿后的同极化通道与交叉极化通道之间的相关性；】

【由于目标的反射非对称性反映到极化SAR影像中表现为同极化通道与交叉极化通道之间的相关性。海杂波通常满足反射对称性特点，而舰船及钻井平台等人造金属目标通常满足反射非对称性特点，因此可以利用T′₁₃的强度值I来进行海上目标检测。】

步骤9：采用双参数恒虚警率检测器进行局部CFAR检测获得每个像素的比较阈值N，局部CFAR检测中所使用的滑动窗口以待检测像素为中心的环形窗口，如图3所示，包括目标窗口、保护窗口、背景窗口，其中，所述目标窗口大小为能包含最小检测目标所占像素个数大小，所述保护窗口为能包含最大检测目标所占像素个数大小，所述背景窗口大小约为保护窗口大小的2倍；

【最大最小目标是根据实际现实情况来判断的，如最大舰船目标的长轴不会超过500米；最大目标的像素不会超出保护窗口，背景窗口是用于海杂波统计，越大统计越详细，但是太大就有可能其他目标落入背景窗口中，影响统计效果，所以一般为保护窗口的两倍多一点；】

N＝u_b+σ_b·n

其中，u_b表示背景窗口的灰度均值，σ_b表示背景窗口的灰度标准差，其中，n表示检测器设计参数，由公式计算获得；Pfa为给定的虚警率，小于10^-8的正数。

【双参数恒虚警率检测器的主旨是寻找出比周围区域异常明亮的像素，基于杂波服从高斯统计模型及给定的虚警概率来求出阈值，进行目标检测。】

步骤10：遍历海上待检测场景中的每一个像素，计算每个像素对应的T′₁₃的强度值I，比较I与阈值N，若I大于阈值N，则当前像素为海上目标，否则为杂波。

【由于海浪扰动产生的方向角偏移影响着海杂波的反射对称性，因此利用海杂波满足反射对称性的特点检测海上目标的效果同样受其影响。通过极化方向角补偿来消除该现象的影响。】

【人造金属目标的反射非对称性在极化SAR影像中是通过同极化通道与交叉极化通道之间的相关性非零来呈现的，即

元素包含了极化SAR影像中同极化通道与交叉极化通道之间的相关性的所有存在情况，并进行了等权组合。相比单独使用或T′₁₃能够更好的描述目标的反射非对称性，适合于海上目标的精确检测。】

应用本发明所述的方法进行海上目标检测，其检测结果如图4所示。为具体呈现本发明的优越性，开辟局部区域，进行分析，如图5所示。从图中可以看出采用方向角补偿前，目标与海杂波难以分离，检测效果差；采用方向角补偿后，目标与海杂波分离明显，该区域的目标全部检测出。

Claims (2)

1.一种基于极化方向角补偿的POLSAR图像海上目标检测方法，其特征在于，包括如下步骤：

步骤1：获取海上待检测场景全极化影像数据即POLSAR数据；

所述全极化影像数据包括HH、HV、VH和VV四种极化通道数据；

步骤2：对POLSAR数据进行预处理操作；

所述预处理操作包括辐射校正、几何校正及陆地掩膜操作；

步骤3：从预处理后的海上待检测场景对应的POLSAR数据中依次提取每个像素对应的S矩阵；

$S=\left[\begin{array}{cc}{S}_{\mathrm{hh}}& S_{\mathrm{hv}}\\ S_{\mathrm{vh}}& S_{\mathrm{vv}}\end{array}\right]=\left[\begin{array}{cc}{S}_{\mathrm{hh}}& S_{\mathrm{hv}}\\ S_{\mathrm{hv}}& S_{\mathrm{vv}}\end{array}\right]$

其中，所述S矩阵满足互易对称定律：S_hv＝S_vh，S_hh表示水平方式发射水平方式接收的极化通道数据，S_hv表示水平方式发射垂直方式接收的极化通道数据，S_vh表示垂直方式发射水平方式接收的极化通道数据，S_vv表示垂直方式发射垂直方式接收的极化通道数据；

步骤4：通过Pauli基进行共轭相乘，将S矩阵转换为相干矩阵T：

$T=<k\&CenterDot;k^H>=\left[\begin{array}{ccc}{T}_{11}& T_{12}& T_{13}\\ T_{21}& T_{22}& T_{23}\\ T_{31}& T_{32}& T_{33}\end{array}\right]$

其中， $K=\frac{1}{\sqrt{2}}\left[\begin{array}{c}{S}_{\mathrm{hh}}+S_{\mathrm{vv}}\\ S_{\mathrm{hh}}-S_{\mathrm{vv}}\\ {2S}_{\mathrm{hv}}\end{array}\right],$

$\left[\begin{array}{ccc}{T}_{11}& T_{12}& T_{13}\\ T_{21}& T_{22}& t_{23}\\ T_{31}& T_{32}& T_{33}\end{array}\right]=\frac{1}{2}\left[\begin{array}{ccc}<{|S_{\mathrm{hh}}+S_{\mathrm{vv}}|}^2>& <(S_{\mathrm{hh}}+S_{\mathrm{vv}}){(S_{\mathrm{hh}}-S_{\mathrm{vv}})}^{*}>& 2<(S_{\mathrm{hh}}+S_{\mathrm{vv}})S_{\mathrm{hv}}^{*}>\\ <(S_{\mathrm{hh}}-S_{\mathrm{vv}}){(S_{\mathrm{hh}}+S_{\mathrm{vv}})}^{*}>& <{|S_{\mathrm{hh}}-S_{\mathrm{vv}}|}^2>& 2<(S_{\mathrm{hh}}-S_{\mathrm{vv}})S_{\mathrm{hv}}^{*}>\\ 2<S_{\mathrm{hv}}{(S_{\mathrm{hh}}+S_{\mathrm{vv}})}^{*}>& 2<S_{\mathrm{hv}}{(S_{\mathrm{hh}}-S_{\mathrm{vv}})}^{*}>& 4<{\left|S_{\mathrm{hv}}\right|}^2>\end{array}\right]$

其中，上标H表示共轭转置，*表示复共轭，<·>表示邻域平均；

步骤5：采用精致Lee滤波，对相干矩阵T进行滤波处理，得到滤波后的相干矩阵T⁰；

$T^0=\left[\begin{array}{ccc}{T}_{11}^0& T_{12}^0& T_{13}^0\\ T_{21}^0& T_{22}^0& T_{23}^0\\ T_{31}^0& T_{32}^0& T_{33}^0\end{array}\right]$

步骤6：利用滤波后的相干矩阵T⁰计算极化方向角：

$\mathrm{\&theta;}=\frac{1}{4}\arctan \left(\frac{2\&CenterDot;\mathrm{Re}\left(T_{23}^0\right)}{T_{22}^0-T_{33}^0}\right)$

其中， $\mathrm{\&theta;}=\frac{1}{4}\arctan \left(\frac{2\&CenterDot;\mathrm{Re}\left(T_{23}^0\right)}{T_{22}^0-T_{33}^0}\right)+\frac{1}{4}\mathrm{\&pi;},$ 若 $\mathrm{\&theta;}>\frac{1}{4}\mathrm{\&pi;},$ 则 $\mathrm{\&theta;}=\mathrm{\&theta;}-\frac{1}{2}\mathrm{\&pi;};$ 表示的实部；θ∈(-π/4,π/4]；

步骤7：利用步骤6中求取出的极化方向角对滤波后的相干矩阵T矩阵进行方向角补偿，得到补偿方向角后的相干矩阵T₀：

$T_0={\mathrm{QT}}^{0}Q=\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos 2\mathrm{\&theta;}& \sin 2\mathrm{\&theta;}\\ 0& -\sin 2\mathrm{\&theta;}& \cos 2\mathrm{\&theta;}\end{array}\right]T^0\left[\begin{array}{ccc}1& 0& 0\\ 0& \cos 2\mathrm{\&theta;}& -\sin 2\mathrm{\&theta;}\\ 0& \sin 2\mathrm{\&theta;}& \cos 2\mathrm{\&theta;}\end{array}\right]=\left[\begin{array}{ccc}{T}_{11}^{\&prime;}& T_{12}^{\&prime;}& T_{13}^{\&prime;}\\ T_{21}^{\&prime;}& T_{22}^{\&prime;}& T_{23}^{\&prime;}\\ T_{31}^{\&prime;}& T_{32}^{\&prime;}& T_{33}^{\&prime;}\end{array}\right]$

其中，T₀表示补偿方向角后的相干矩阵，Q为旋转矩阵；

步骤8：提取补偿方向角后的相干矩阵T₀中的元素T′₁₃，计算T′₁₃的强度值I；

步骤9：采用恒虚警率检测器求取每个像素的检测阈值N；

步骤10：遍历海上待检测场景中的每一个像素，计算每个像素对应的T′₁₃的强度值I，比较I与阈值N，若I大于阈值N，则当前像素为海上目标，否则为杂波。

2.根据权利要求1所述的基于极化方向角补偿的POLSAR图像海上目标检测方法，其特征在于，在所述步骤9中采用双参数恒虚警率检测器进行局部CFAR检测获得每个像素的比较阈值N，局部CFAR检测中所使用的滑动窗口以待检测像素为中心的环形窗口，包括目标窗口、保护窗口、背景窗口，其中，所述目标窗口大小为能包含最小检测目标所占像素个数大小，所述保护窗口为能包含最大检测目标所占像素个数大小，所述背景窗口大小为保护窗口大小的1.5-2.5倍；

N＝u_b+σ_b·n

其中，u_b表示背景窗口的灰度均值，σ_b表示背景窗口的灰度标准差，其中，n表示检测器设计参数，由公式计算获得；Pfa为给定的虚警率，小于10^-8的正数。