一种基于张量表示的偏振高光谱图像目标检测方法
【专利摘要】一种基于张量表示的偏振高光谱图像目标检测方法，本发明涉及偏振高光谱图像目标检测方法。本发明的目的是为了解决现有偏振高光谱图像目标检测准确率低以及信息利用率低的问题。本发明组成包括：步骤一、建立张量表示下的偏振高光谱图像的数据模型：将偏振高光谱图像表示成四阶张量，则得到偏振高光谱图像的数据模型；步骤二、在张量表示下的偏振高光谱图像的数据模型的基础上，得到基于Tucker分解的四阶张量匹配检测算法的模型，根据基于Tucker分解的四阶张量匹配检测算法的模型判断检测结果。本发明用于偏振高光谱图像目标检测。
【专利说明】
一种基于张量表示的偏振高光谱图像目标检测方法
技术领域
[0001] 本发明涉及偏振高光谱图像目标检测方法。
【背景技术】
[0002] Stokes矢量较为完整的表示了光的偏振特性，比较精确的描述了物体的组成成分 及其表面特性。但是对于偏振高光谱图像而言，其具有的偏振、强度、光谱以及空间的多维 信息，Stokes矢量很难完整的描述出来，它忽略了偏振高光谱图像随偏振角变化的偏振维 信息;并且传统的偏振高光谱图像检测算法也忽略了图像的空间信息，很难实现偏振高光 谱图像多维信息的共同利用。此外，当偏振高光谱图像的数据量较大时，对于Stokes矢量的 计算量也将增加，计算误差大，综合以上，导致偏振高光谱图像目标检测准确率低。
【发明内容】
[0003] 本发明的目的是为了解决现有偏振高光谱图像目标检测准确率低以及信息利用 率低的问题，而提出基于张量表示的偏振高光谱图像目标检测方法。
[0004] -种基于张量表示的偏振高光谱图像目标检测方法具体是按照以下步骤进行的：
[0005] 步骤一、建立张量表示下的偏振高光谱图像的数据模型；
[0006]将偏振高光谱图像表示成四阶张量￥ ，则偏振高光谱图像的数据模型 表示为
[0007] y:=S'+'M (2-3)
[0008] 1表不信号张量，f表不干扰信号，In为图像第η维像素大小，η取值为1、2、3或4， R为实数域空间；
[0009] 步骤二、在张量表示下的偏振高光谱图像的数据模型的基础上，得到基于Tucker 分解的四阶张量匹配检测算法的模型，根据基于Tucker分解的四阶张量匹配检测算法的模 型判断检测结果。
[0010] -种基于张量表示的偏振高光谱图像目标检测方法具体是按照以下步骤进行的：
[0011] 步骤一、建立张量表示下的偏振高光谱图像的数据模型；
[0012] 将偏振高光谱图像表示成张量FeRWW，〗"为图像第n维像素大小，n取值为K 2、3或4，R为实数域空间，则偏振高光谱图像的数据模型表示为：
[0013] y=S + ：N' (2>3)
[0014] JT表不信号张量，表不干扰信号；
[0015] 步骤二、在张量表示下的偏振高光谱图像的数据模型的基础上，构造四阶张量，对 构造的四阶张量进行CP分解重建，然后采用常用的统计异常检测（RX)和匹配滤波检测 (MF)，得到检测结果。
[0016] 本发明的有益效果为：
[0017] 为了解决传统的Stokes矢量对于偏振高光谱图像目标检测信息利用不完整的问 题，本发明研究利用张量表示来完整描述偏振高光谱图像的多维信息，并实现偏振与光谱 在偏振高光谱图像目标检测的联合应用。张量表示完整保留了图像的多维信息，并且在检 测过程中综合利用图像的这些信息，增加了目标与背景的差异，使得目标检测更加准确。本 发明提出的FTMF算法不仅完整地利用了图像光谱维和偏振维的信息，还大大减化了计算过 程，减少了由于计算而带来的误差。由于CP分解具有保留原始信号特征，去除干扰信号的作 用，在相同虚警率下，未经处理的偏振高光谱图像检测效果要低于经过CP分解重建后的检 测效果，实验结果证明了本发明提出的算法的有效性。考虑了偏振高光谱图像随偏振角变 化的偏振维信息以及图像的空间信息，实现偏振高光谱图像多维信息的共同利用，提高了 偏振高光谱图像目标检测准确率。
[0018] 对于提出的四阶张量匹配检测算法，本发明将其检测结果与传统的基于Stokes矢 量提取的强度I、偏振度DoLP以及偏振角orient的匹配检测结果作对比，如图6a、6b、6c、6d、 6e所示，根据其ROC曲线，发现，在相同的虚警率下，本发明提出的四阶张量匹配检测算法的 检测效率高于数据集I、D 〇LP和orient的匹配检测效率。这是由于常规的基于Stokes矢量的 偏振高光谱图像目标检测算法破坏了偏振维的信息，而本发明提出的FTMF算法不仅完整地 利用了图像光谱维和偏振维的信息，还大大减化了计算过程，减少了由于计算而带来的误 差。当虚警概率为0.1时，偏振度Dolp检测准确率为92.5%，强度I检测准确率为99.2%，偏 振角orient检测准确率为98.8%，FTMF检测准确率为99.9% ;当虚警概率为0.15时，偏振度 Dolp检测准确率为94.7%，强度I检测准确率为99.8%，偏振角orient检测准确率为99%， FTMF检测准确率基本已达到100% ;
[0019] 为了进一步直观看到算法的检测结果，这里给出检测结果的三维视图，如图7a、 7b、7c、7d所示，同样可以发现，本发明提出的算法能够有效的抑制背景、突出目标，使得目 标和背景能够更准确地被分离出来。
[0020] 对于本发明提出的基于CP分解重建的偏振高光谱图像目标检测，为了验证其有效 性，分别对CP分解重建后的数据采用RX和MF检测算法，后续将其命名为CPRX和CPMF，检测结 果如图8 &、813、8〇、8(1、86和图9&、％所示，分别展示了两种检测算法对应的检测结果及其1?0〇 的对比曲线，可以发现，由于CP分解具有保留原始信号特征，去除干扰信号的作用，在相同 虚警率下，未经处理的偏振高光谱图像检测效果要低于经过CP分解重建后的检测效果，实 验结果证明了本发明提出的算法的有效性。
[0021] 检测结果的三维视图如图10a、10b、10c、10d、10e所示，从圈内目标的峰值可以看 出，经过CP分解重建的偏振高光谱数据的检测结果的目标与背景基本被分离，能够更精确 地从背景中找出目标，原因在于当光谱维加入偏振维信息后，目标与背景的差异被进一步 放大，使得数据更加契合基于匹配滤波的检测算法和基于统计规则的RX检测算法。
[0022] 如图11所示，当虚警概率为0.01时，CPRX检测准确率为80 %，CPMF检测准确率为 100%，FTMF检测准确率为100% ;当虚警概率为0.02时，CPRX检测准确率为87%，CPMF检测 准确率为100 %，FTMF检测准确率为99.8 % ;当虚警概率为0.05时，CPRX检测准确率为94 %， CPMF检测准确率为100 %，FTMF检测准确率为100 %。
【附图说明】
[0023]图la为一个三阶张量和它的各模纤示意图；
[0024]图lb为三阶张量和它的模-1纤示意图；
[0025]图lc为三阶张量和它的模-2纤示意图；
[0026]图Id为三阶张量和它的模-3纤示意图；
[0027]图2为张量的模-1乘积过程图，A为张量，B为张量；
[0028]图3为三阶张量的Tucker分解图，X为张量，w为一维数据，v为二维数据，u为三维数 据，s为核心张量；
[0029]图4为三阶张量的CP分解图，X为张量，W1为第一个一维数据，V1为第一个二维数据， m为第一个三维数据，w2为第二个一维数据，ν2为第二个二维数据，u2为第二个三维数据，職 为第R个一维数据， VR为第R个二维数据，uR为第R个三维数据，U、V、W为因子矩阵，λ是一个常 数；
[0030]图5为CP分解重建的过程图，y为张量，y'为CP分解重建后的张量，ai为第一个一维 数据，bi为第一个二维数据，&为第一个三维数据，cU为第一个四维数据，&2为第二个一维数 据，b 2为第二个二维数据，C2为第二个三维数据，d2为第二个四维数据，&"为第三个一维数 据，b n为第三个二维数据，Cn为第三个三维数据，dn为第三个四维数据，row为行，column为 列，band为波段，polarization为偏振；
[0031]图6a为偏振高光谱图像的数据集I的传统MF检测结果图；
[0032]图6b为偏振高光谱图像的数据集DoLP的传统MF检测结果图；
[0033]图6c为偏振高光谱图像的数据集orient的传统MF检测结果图；
[0034]图6d为偏振高光谱图像的数据集FTMF的传统MF检测结果图；
[0035]图6e为偏振高光谱图像的传统MF检测结果和FTMF检测结果对比图，false alarm rate为虚警概率，detection rate为检测准确率；
[0036]图7a为偏振高光谱图像的数据集I的MF结果三维视图；
[0037]图7b为偏振高光谱图像的数据集DoLP的MF结果三维视图；
[0038]图7c为偏振高光谱图像的数据集orient的MF结果三维视图；
[0039] 图7d为偏振高光谱图像的原始数据的FTMF结果三维视图；
[0040] 图8a为偏振高光谱图像数据集I的RX结果图；
[0041]图8b为偏振高光谱图像数据集DoLP的RX结果图；
[0042]图8c为偏振高光谱图像数据集orient的RX结果图；
[0043]图8d为CP分解重建的偏振高光谱图像RX结果图；
[0044]图8e为RX检测结果的R0C曲线对比图；
[0045] 图9a为基于CP分解重建的偏振高光谱图像MF检测结果图；
[0046] 图9b为基于CP分解重建的偏振高光谱图像MF检测结果的R0C曲线对比图；
[0047]图10a为基于CP分解重建的偏振高光谱图像数据集I的RX结果三维视图；
[0048]图10b为基于CP分解重建的偏振高光谱图像数据集DoLP的RX结果三维视图；
[0049]图10c为基于CP分解重建的偏振高光谱图像数据集orient的RX结果三维视图； [0050]图l〇d为CP分解重建的偏振高光谱图像RX结果三维视图；
[00511图l0e为CP分解重建的偏振高光谱图像MF结果三维视图；
[0052] 图11为偏振高光谱图像CPRX、CPMF和FTMF检测结果的R0C曲线对比图，FTMF为本发 明基于张量表示的偏振高光谱图像目标检测方法。
【具体实施方式】
[0053]【具体实施方式】一:本实施方式的一种基于张量表不的偏振高光谱图像目标检测方 法具体是按照以下步骤进行的：
[0054] 张量的概念最早是18世纪90年代末由Woldemar Voigt提出的，后来，Ricci和他的 学生在研究中逐渐完善了张量的概念和相关的运算。从代数的角度看，张量可以看成是矩 阵的进一步扩展，矩阵的行是模-1纤，矩阵的列是模-2纤，三阶张量则多出了管纤，如图la、 lb、lc、ld所示，展示了一个三阶张量的格式。对于向量空间U (1)，U(2)'. . .，U(M)，定义他们的 外积空间Vn) ?L/(2) Θ...??/(Μ)为张量。
[0055] 对于后续涉及到的关于张量的基本运算，本节做一个简单介绍。首先是张量的η-模展开，对于张量安照η模进行展开指的是将张量X的η模作为列展开成矩阵形 式，例如张量X为： rl 4 Γιο 13 16λ
[0056] ΧΓ= 2 5 8 Χ2= 11 14 17 ν3 6 9J [?2 15 l.Sy
[0057] 则张量X的η-模展开分别为： ? 4 7 10 13 16、
[0058] 张量的 1-模展开Χ(ι)= 2 5 8 11 14 17 ν3 6 9 12 15 18y· ? 2 3 10 11 12'
[0059] 张量的2-模展开X(2:厂 4 5 6 13: 14 1:5 ,7 8 9 16 17 18y _ (\ 2 3 4 5 6 7 8
[0060] 张量的3_模展开....... m里口j I矢肷力? j〇 11 12 13 14 15 16 17 18;
[0061 ]张量与矩阵的乘法被称为η-模乘积，对于张量X e 和矩阵U e 5/?的n- 模乘积可以表示为
[0062] I = XxB U=UX(I1) (〇.Γ)
[0063] 即每个模-η向量与矩阵U做乘积的结果。图2给出了张量与矩阵1-模乘积过程的图 不。
[0064]基于Tucker分解的四阶张量匹配检测算法
[0065] Tucker分解是指，将一个张量表示成一个核心张量沿着每一个模式乘上一个因子 矩阵的形式，以一个三阶张量X e 的Tucker分解为例，其分解形式为：
[0066] X - 5 X, U x： V X, W = q^l r^\ (2-1)
[0067] 其中，SeRPXQXR是核心张量，U、V、W被称之为因子矩阵，三阶张量Tucker分解如图3 所示
[0068]可以发现，Tucker分解实质上是一种高阶的主成分分析，因子矩阵U、V、W为正交矩 阵，可以视为沿对应模的主成分。因此，在Tucker分解的基础上，本发明提出了一种适合于 偏振高光谱图像目标检测的四阶张量匹配检测算法。
[0069]传统的匹配检测算法(MF)是一种基于线性最优准则假设推导而来的检测算法，在 高光谱图像的目标检测中，主要利用了图像的光谱信息，可以得到比较理想的检测结果，算 法的基本形式为：
丨/乏目.标 背景 （2_2)
[0071] 其中，ε为输出结果的分割阈值，Rx是待测数据X的协方差矩阵，s则表示目标的光 谱信息。
[0072] 偏振高光谱图像不仅包含了高光谱图像的信息，还多出了偏振维度的信息，因此 根据张量乘法与匹配检测算法的原理，本发明采用了一种四阶张量匹配检测算法实现了偏 振高光谱图像的目标检测。
[0073] 步骤一、建立张量表示下的偏振高光谱图像的数据模型；
[0074] 将偏振高光谱图像表示成四阶张量偏振高光谱图像的数据模型 表示为
[0075] y 二S + ：N. (2-3)
[0076] X表示信号张量，f表示干扰信号，1"为图像第11维像素大小，11取值为1、2、3或4， R为实数域空间；
[0077] 步骤二、在张量表示下的偏振高光谱图像的数据模型的基础上，得到基于Tucker 分解的四阶张量匹配检测算法(FTMF)的模型，根据基于Tucker分解的四阶张量匹配检测算 法(FTMF)的模型判断检测结果。
【具体实施方式】 [0078] 二:本实施方式与一不同的是:所述步骤二中在张量 表示下的偏振高光谱图像的数据模型的基础上，得到基于Tucker分解的四阶张量匹配检测 算法(FTMF)的模型，根据基于Tucker分解的四阶张量匹配检测算法(FTMF)的模型判断检测 结果;具体过称为：
[0079] 对训练样本采用Tucker分解，提取目标光谱维和偏振维的信息，选择目标光谱维 和偏振维的因子矩阵的第一主成分(即第一列)来替代目标的光谱和偏振谱，然后根据光谱 维和偏振维的信息构造光谱匹配滤波器和偏振匹配滤波器，得到基于Tucker分解的四阶张 量匹配检测算法(FTMF)的模型，根据四阶张量匹配检测算法(FTMF)的模型判断检测结果。
[0080] 其它步骤及参数与【具体实施方式】一相同。
[0081 ]【具体实施方式】三:本实施方式与【具体实施方式】一或二不同的是:所述根据光谱维 和偏振维的信息构造光谱匹配滤波器和偏振匹配滤波器，得到基于Tucker分解的四阶张量 匹配检测算法(FTMF)的模型，根据四阶张量匹配检测算法(FTMF)的模型判断检测结果;具 体过称为：
[0082] 光谱匹配滤波器H3构造如下：
α-4)
[0084] 其中，S3表示目标光谱，Rs=E(Y⑶Υωτ)为待检测数据Y⑶的协方差矩阵，x':, 表示四阶张量J/的3-模展开，In为图像第η维大小，η取值为1、2、3或4，T为转置；
[0085] 将对光谱匹配滤波器Η3和四阶张量J/.做3-模乘积有：
[0086] ~ -^3 -^3 (.2-5)
[0087]式中，X3为张量与矩阵的3模乘积；
[0088]偏振匹配滤波器Η4的构造如下：
(2-6)
[0090]其中，S4是目标偏振谱，为待测数据Ts⑷的协方差矩阵，Ts⑷表示 四阶张量y的4-模展开矩阵；
[0091 ] 根据式(0-2)~(0-3)，基于Tucker分解的四阶张量匹配检测算法(FI1F)的模型如下：
[0092] - ^ s X4 = J^ x：a Η ； χ4 Η4 (2-.7)..
[0093] 式中，Χ4为张量与矩阵的4模乘积；
[0094] 当μ时，待检测像素被判定为目标；
[0095] 当μ时，待检测像素被判定为背景;μ为阈值。
[0096] 其它步骤及参数与【具体实施方式】一或二相同。
【具体实施方式】 [0097] 四：本实施方式的一种基于张量表示的偏振高光谱图像目标检测方 法具体是按照以下步骤进行的：
[0098] 张量的概念最早是18世纪90年代末由Woldemar Voigt提出的，后来，Ricci和他的 学生在研究中逐渐完善了张量的概念和相关的运算。从代数的角度看，张量可以看成是矩 阵的进一步扩展，矩阵的行是模-1纤，矩阵的列是模-2纤，三阶张量则多出了管纤，如图la、 lb、lc、ld所示，展示了一个三阶张量的格式。对于向量空间U (1)，U(2)'. . .，U(M)，定义他们的 外积空间f/(1) ?"(2) ?为张量。
[0099] 对于后续涉及到的关于张量的基本运算，本节做一个简单介绍。首先是张量的n-模展开，对于张量Xei^x〃···~按照η模进行展开指的是将张量X的η模作为列展开成矩阵 形式，例如张量X为： "1 4 7> (:1.0 Β 1:6、
[0100] XI：= 2 5 ：8 X, = 11 14 17 3 6 9」 [12 15 l：Sy
[0101] 则张量X的η-模展开分别为： ? 4 7 10 13 16'
[0102] 张量的 1-模展开Χ(1)= 2 51 11 14 17 ,3 (S 9 12 15 18y ? 23 10 11 12'
[0103] 张量的2-模展开Χρ)= 4 5 6 13 Η 15 b 8 9 16 17 lSy. (\ 2 3 4 56789)
[0104] 张里的3-模展开X(,「ju) n 12 13 14 15 16 17 %
[0105] 张量与矩阵的乘法被称为n-模乘积，对于张量Xe沪#···^和矩阵U e 的n- 模乘积可以表示为
[0106] !/ =义'xnU=UX⑷ f ,_j )
[0107] 即每个模-η向量与矩阵U做乘积的结果。图2给出了张量与矩阵1-模乘积过程的图 不。
[0108] 基于Tucker分解的四阶张量匹配检测算法
[0109] Tucker分解是指，将一个张量表示成一个核心张量沿着每一个模式乘上一个因子 矩阵的形式，以一个三阶张量X e 的Tucker分解为例，其分解形式为：
[0110] X。s 'U X2V x3W mvjvwh.sPil,. (2-1) p=l r=\
[0111] 其中，SeRPXQXR是核心张量，U、V、W被称之为因子矩阵，三阶张量Tucker分解如图3 所示
[0112] 可以发现，Tucker分解实质上是一种高阶的主成分分析，因子矩阵U、V、W为正交矩 阵，可以视为沿对应模的主成分。因此，在Tucker分解的基础上，本发明提出了一种适合于 偏振高光谱图像目标检测的四阶张量匹配检测算法。
[0113] 传统的匹配检测算法(MF)是一种基于线性最优准则假设推导而来的检测算法，在 高光谱图像的目标检测中，主要利用了图像的光谱信息，可以得到比较理想的检测结果，算 法的基本形式为：
目标 背景 （2-2)
[0115]其中，ε为输出结果的分割阈值，Rx是待测数据X的协方差矩阵，s则表示目标的光 谱信息。
[0116]偏振高光谱图像不仅包含了高光谱图像的信息，还多出了偏振维度的信息，因此 根据张量乘法与匹配检测算法的原理，本发明采用了一种四阶张量匹配检测算法实现了偏 振高光谱图像的目标检测。
[0117]张量的CP分解指的是，将张量分解成有限个秩一张量之和的形式，可以被看作是 Tucker分解的一种特殊形式。同样以三阶张量为例，CP分解可以表示为
(3-1)
[0119]其中Ar是一个常数，Ur、Vr、Wr表示每个维度上的向量。图4展示了一个三阶张量的 CP分解示意图。
[0120] CP分解和Tucker分解一样，每个维度上的主要信息都被保留在其相应的向量上， CP分解后能够在保留原始信号特征的同时舍去信号的冗余信息，降低了噪声的影响，因此 基于CP分解重建的偏振高光谱图像目标检测的核心是利用CP分解重建偏振高光谱数据，使 其成为常规的三维数据，在保留图像多维信息的同时使用常规的检测算法对其进行目标检 测。
[0121] 步骤一、建立张量表示下的偏振高光谱图像的数据模型；
[0122] 将偏振高光谱图像表示成张量y e RV&y4，In为图像第n维像素大小，n取值为i、 2、3或4，则偏振高光谱图像的数据模型表示为：
[0123] y=S + -%r (2-3)
[0124] 表示信号张量，f表示干扰信号，R为实数域空间；
[0125] 步骤二、在张量表示下的偏振高光谱图像的数据模型的基础上，提出基于CP分解 重建的偏振高光谱图像目标检测算法，构造四阶张量，对构造的四阶张量进行CP分解重建， 然后采用常用的统计异常检测(RX)和匹配滤波检测(MF)，得到检测结果。
[0126] 对构造的四阶张量进行CP分解，得到检测结果。
【具体实施方式】 [0127] 五:本实施方式与一不同的是:所述步骤二中在张量 表示下的偏振高光谱图像的数据模型的基础上，提出基于CP分解重建的偏振高光谱图像目 标检测算法，构造四阶张量，对构造的四阶张量进行CP分解重建，然后采用常用的统计异常 检测(RX)和匹配滤波检测(MF)，得到检测结果;具体过称为：
[0128] 步骤二一、对原始偏振高光谱图像张量中的每一个像素选择其周围kXk邻域，构 造四阶张量;:k为分块图像大小；
[0129] 步骤二二、对分块图像进行CP分解；
[0130]步骤二三、将CP分解的第三维和第四维投影向量重组为一维向量替代原始偏振高 光谱图像中心的光谱和偏振，形成一个新的三维数据；
[0131] 步骤二四、采用常规的统计异常检测（RX)算法，匹配滤波(MF)算法对新的三维数 据进行目标检测。
[0132] 其它步骤及参数与【具体实施方式】一相同。
[0133] 采用以下实施例验证本发明的有益效果：
[0134] 实施例一：
[0135] 本实施例基于张量表示的偏振高光谱图像目标检测方法，具体是按照以下步骤制 备的：
[0136] 发明采用实测的偏振高光谱图像来进行实验，并且与传统的基于Stokes矢量的偏 振高光谱图像检测算法进行对比，为了验证算法的有效性和实用性，引入R0C曲线来进行算 法的评价。
[0137] 对于提出的四阶张量匹配检测算法，本发明将其检测结果与传统的基于Stokes矢 量提取的强度I、偏振度DoLP以及偏振角orient的匹配检测结果作对比，如图6a、6b、6c、6d、 6e所示，根据其ROC曲线，可以发现，在相同的虚警率下，本发明提出的四阶张量匹配检测算 法的检测效率高于数据集I、D 〇LP和orient的匹配检测效率。这是由于常规的基于Stokes矢 量的偏振高光谱图像目标检测算法破坏了偏振维的信息，而本发明提出的FTMF算法不仅完 整地利用了图像光谱维和偏振维的信息，还大大减化了计算过程，减少了由于计算而带来 的误差。当虚警概率为0.1时，偏振度Dolp检测准确率为92.5%，强度I检测准确率为 99.2%，偏振角orient检测准确率为98.8%，FTMF检测准确率为99.9% ;当虚警概率为0.15 时，偏振度Dolp检测准确率为94.7%，强度I检测准确率为99.8%，偏振角orient检测准确 率为99%，FTMF检测准确率为100% ;
[0138] 为了进一步直观看到算法的检测结果，这里给出检测结果的三维视图，如图7a、 7b、7c、7d所示，同样可以发现，本发明提出的算法能够有效的抑制背景、突出目标，使得目 标和背景能够更准确地被分离出来。
[0139] 对于本发明提出的基于CP分解重建的偏振高光谱图像目标检测，为了验证其有效 性，分别对CP分解重建后的数据采用RX和MF检测算法，检测结果如图8a、8b、8c、8d、8e和图 9a、9b所示，分别展示了两种检测算法对应的检测检测结果及其ROC的对比曲线，可以发现， 由于CP分解具有保留原始信号特征，去除干扰信号的作用，在相同虚警率下，未经处理的偏 振高光谱图像检测效果要低于经过CP分解重建后的检测效果，实验结果证明了本发明提出 的算法的有效性。
[0140] 检测结果的三维视图如图10a、10b、10c、10d、10e所示，从红圈内目标的峰值可以 看出，经过CP分解重建的偏振高光谱数据的检测结果的目标与背景基本被分离，能够更精 确地从背景中找出目标，原因在于当光谱维加入偏振维信息后，目标与背景的差异被进一 步放大，使得数据更加契合基于匹配滤波的检测算法和基于统计规则的RX检测算法。
[0141] 如图11所示，当虚警概率为0.01时，CPRX检测准确率为80%，CPMF检测准确率为 100%，FTMF检测准确率为100% ;当虚警概率为0.02时，CPRX检测准确率为87%，CPMF检测 准确率为100 %，FTMF检测准确率为99.8 % ;当虚警概率为0.05时，CPRX检测准确率为94 %， CPMF检测准确率为100%，FTMF检测准确率为100%。
[0142] 本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域 技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于 本发明所附的权利要求的保护范围。
【主权项】
1. 一种基于张量表示的偏振高光谱图像目标检测方法，其特征在于：一种基于张量表 示的偏振高光谱图像目标检测方法具体是按照W下步骤进行的： 步骤一、建立张量表示下的偏振高光谱图像的数据模型； 将偏振高光谱图像表示成四阶张量则偏振高光谱图像的数据模型表示 为 (2-3) X表示信号张量，表示干扰信号，In为图像第η维像素大小，η取值为1、2、3或4，R为实 数域空间； 步骤二、在张量表示下的偏振高光谱图像的数据模型的基础上，得到基于化cker分解 的四阶张量匹配检测算法的模型，根据基于化cker分解的四阶张量匹配检测算法的模型判 断检测结果。2. 根据权利要求1所述一种基于张量表示的偏振高光谱图像目标检测方法，其特征在 于:所述步骤二中在张量表示下的偏振高光谱图像的数据模型的基础上，得到基于化cker 分解的四阶张量匹配检测算法的模型，根据基于化cker分解的四阶张量匹配检测算法的模 型判断检测结果;具体过称为： 对训练样本采用Tucker分解，提取目标光谱维和偏振维的信息，选择目标光谱维和偏 振维的因子矩阵的第一主成分来替代目标的光谱和偏振谱，然后根据光谱维和偏振维的信 息构造光谱匹配滤波器和偏振匹配滤波器，得到基于化cker分解的四阶张量匹配检测算法 的模型，根据四阶张量匹配检测算法的模型判断检测结果。3. 根据权利要求2所述一种基于张量表示的偏振高光谱图像目标检测方法，其特征在 于：所述根据光谱维和偏振维的信息构造光谱匹配滤波器和偏振匹配滤波器，得到基于 化cker分解的四阶张量匹配检测算法的模型，根据四阶张量匹配检测算法的模型判断检测 结果;具体过称为： 光谱匹配滤波器出构造如下：其中，S3表示目标光谱，1?3 = 6(￥(3)￥(3)了）为待检测数据￥(3)的协方差矩阵， e 'x"i1表示四阶张量y的3-模展开矩阵，T为转置. 对光谱匹配滤波器H3和四阶张量J/做3-模乘积有： (2-5) 式中，X3为张量与矩阵的3模乘积； 偏振匹配滤波器此的构造如下：其中，S4是目标偏振谱，R4 = E(Ts(4)Ts(4)T)为待测数据Tsw的协方差矩阵，Ts(4康示四阶 张量的4-模展开矩阵； 根据式(2-5)~(2-6)，基于化cker分解的四阶张量匹配检测算法的模型如下： 《習=了.、. X,马=.//X，.申 X4 馬(2-7) 式中，X4为张量与矩阵的4模乘积； 当> μ时，待检测像素被判定为目标； 当< ^1时，待检测像素被判定为背景;μ为阔值。4. 一种基于张量表示的偏振高光谱图像目标检测方法，其特征在于：一种基于张量表 示的偏振高光谱图像目标检测方法具体是按照W下步骤进行的： 步骤一、建立张量表示下的偏振高光谱图像的数据模型； 将偏振高光谱图像表示成张量Fe，?η为图像第η维像素大小，η取值为1、2、3或 4，则偏振高光谱图像的数据模型表示为： (2-3) J'表示信号张量，表示干扰信号，R为实数域空间； 步骤二、在张量表示下的偏振高光谱图像的数据模型的基础上，构造四阶张量，对构造 的四阶张量进行CP分解重建，然后采用统计异常检测和匹配滤波检测，得到检测结果。5. 根据权利要求4所述一种基于张量表示的偏振高光谱图像目标检测方法，其特征在 于:所述步骤二中在张量表示下的偏振高光谱图像的数据模型的基础上，构造四阶张量，对 构造的四阶张量进行CP分解重建，然后采用统计异常检测和匹配滤波检测，得到检测结果； 具体过称为： 步骤二一、对原始偏振高光谱图像张量中的每一个像素选择其周围kXk邻域，构造四 阶张量y € Rtxw'xSk为分块图像大小，In为图像第η维像素大小，η取值为1、2、3或4,； 步骤二二、对分块图像进行CP分解； 步骤二Ξ、将CP分解的第Ξ维和第四维投影向量重组为一维向量替代原始偏振高光谱 图像中屯、的光谱和偏振，形成一个新的Ξ维数据； 步骤二四、采用常规的统计异常检测算法，匹配滤波算法对新的Ξ维数据进行目标检 测。
【文档编号】G06T7/00GK106097321SQ201610393847
【公开日】2016年11月9日
【申请日】2016年6月6日 公开号201610393847.X, CN 106097321 A, CN 106097321A, CN 201610393847, CN-A-106097321, CN106097321 A, CN106097321A, CN201610393847, CN201610393847.X
【发明人】张钧萍, 谭建
【申请人】哈尔滨工业大学