基于张量匹配子空间的高光谱图像目标检测方法
技术领域
本发明涉及高光谱图像目标检测方法。
背景技术
高光谱传感器主要通过在不同的光学或近红外波段进行地物反射信息提取,获取地面的图-谱特征,高光谱图像包含了地物的光谱特征,地物的空间特征,这种几乎连续的光谱采样信息,可以反映地物在光谱上很小的差异变化,因此这种特征常被称为诊断特征,作为对地物进行分类和检测的依据。研究高光谱图像目标检测新技术,具有重要的理论意义和应用价值。在军事方面,可以揭露敌方目标的伪装、隐藏和欺骗。在民用方面,公共安全、食品安全、质量监控、森林着火点检测、失事点搜寻与营救等方面已经有重要应用。
现有典型的目标检测方法有基于有限冲积响应滤波器的约束能量最小化方法(constrained energy minimization,CEM),正交子空间目标检测方法(orthogonalsubspace projection,OSP),匹配子空间检测方法(matched subspace detector,MSD),自适应子空间检测方法(adaptive subspace detector,ASD)等经典方法以及近几年提出的稀疏表示检测方法(sparse representation,SR)等。
当前目标检测主要利用光谱匹配特性,通过光谱匹配程度判断单点光谱的属性,并未充分考虑在高分条件下空间约束增强的特性,即局部相关性增大的特性,空间-光谱联合检测方法可以提高光谱和空间的利用能力,但一般的空间-光谱联合方法仅仅从光谱或空间的简单组合操作进行分析,而不能从三维数据整体进行信息挖掘,导致高光谱图像目标检测精度低以及空间信息利用不足。
发明内容
本发明的目的是为了解决现有高光谱图像目标检测精度低以及空间信息利用率低的问题,而提出基于张量匹配子空间的高光谱图像目标检测方法。
具体过程为:
步骤一、建立张量表示下的目标样本H1和背景样本H0的信号表示模型;
步骤二、根据高光谱样本数据和窗口大小,分别建立目标样本H1和背景样本H0的空间X、空间Y、光谱和原子的四阶张量矩阵;
步骤三、依据张量匹配子空间投影算法,求取目标样本H1和背景样本H0的空间X、空间Y、光谱和原子的四阶张量矩阵中空间X、空间Y和光谱三个背景方向和三个目标方向的正交投影矩阵;
待检测信号在三个背景方向正交投影矩阵映射下的数据空间为背景样本投影子空间,待检测信号在三个目标方向正交投影矩阵映射下的数据空间为目标样本投影子空间;
步骤四、提取待检测信号的空间X、空间Y和光谱的三阶张量矩阵,分别利用步骤三得到的正交投影矩阵将待检测信号映射到步骤三得到的目标样本投影子空间和背景样本投影子空间中;
步骤五、分别计算待检测信号与步骤四得到的目标样本和背景样本投影子空间投影下的信号的误差,依据张量表达下的广义似然比检测模型和给定的阈值,判定待检测信号是否为检测目标,如果张量表达下的广义似然比检测模型值大于等于η时,为检测目标;否则为背景目标;其中η为阈值。
本发明的有益效果为:
由于本发明方法将高光谱数据转换成张量形式进行处理,能够对多维数据整体信息进行挖掘,相比于当前目标检测单独利用光谱匹配特性,以及一般的仅通过光谱和空间的简单组合操作的空间-光谱联合方法,克服了不能从三维数据整体进行信息挖掘的缺点。提高了空间信息利用率以及高光谱图像目标检测的精度。
为了验证本发明所提出的方法的性能,针对一组机载可见光/红外成像光谱仪(Airborne Visible Infrared Imaging Spectrometer,AVIRIS)采集的一组机场数据进行了实验,实验结果验证了本发明提出的基于张量匹配子空间的高光谱图像目标检测方法的有效性。当虚警概率为10-2时,本发明方法的检测精度达到89%,比传统方法检测精度提高了10%以上。
附图说明
图1是本发明的实现流程示意图;
图2a是信号张量表示条件下背景H0的表示模型图,为待检测信号,和为从背景库中选取的一些背景样本的张量原子,a1、a2和an分别为和相对应的权重系数;
图2b是信号张量表示条件下目标H1的表示模型图,和为从目标库中选取的一些目标样本的张量原子,b1、b2和bm为和相对应的权重系数;
图3a是飞机场真实高光谱数据的一个波段的灰度图;
图3b是飞机场数据基于张量匹配子空间的高光谱图像目标检测结果图;
图3c是飞机场数据高光谱图像匹配子空间检测结果图;
图4是飞机场数据上述两种方法检测结果的接收机工作ROC曲线图,Tensor MSD为本发明基于张量匹配子空间的高光谱图像目标检测方法,MSD为传统的高光谱图像目标检测方法,Detection power为检测精度,False alarm rate为虚警概率。
具体实施方式
具体实施方式一:结合图1说明本实施方式,本实施方式的基于张量匹配子空间的高光谱图像目标检测方法,具体过称为:
步骤一、建立张量表示下的目标样本H1和背景样本H0的信号表示模型;
步骤二、根据高光谱样本数据和窗口大小,分别建立目标样本H1和背景样本H0的空间X、空间Y、光谱和原子的四阶张量矩阵;
步骤三、依据张量匹配子空间投影算法,求取目标样本H1和背景样本H0的空间X、空间Y、光谱和原子的四阶张量矩阵中空间X、空间Y和光谱三个背景方向和三个目标方向的正交投影矩阵;
待检测信号在三个背景方向正交投影矩阵映射下的数据空间为背景样本投影子空间,待检测信号在三个目标方向正交投影矩阵映射下的数据空间为目标样本投影子空间;
步骤四、提取待检测信号的空间X、空间Y和光谱的三阶张量矩阵,分别利用步骤三得到的正交投影矩阵将待检测信号映射到步骤三得到的目标样本投影子空间和背景样本投影子空间中;
步骤五、分别计算待检测信号与步骤四得到的目标样本和背景样本投影子空间投影下的信号的误差,依据张量表达下的广义似然比检测模型和给定的阈值,判定待检测信号是否为检测目标,如果张量表达下的广义似然比检测模型值大于等于η时,为检测目标;否则为背景目标;其中η为阈值。
具体实施方式二:结合图2a、图2b说明本实施方式,本实施方式与具体实施方式一不同的是:所述步骤一中建立张量表示下的目标样本H1和背景样本H0的信号表示模型;具体过程为:
其中,表示待检测信号的三阶张量表示;表示背景样本形成的四阶张量子空间;表示目标样本形成的四阶张量子空间;x4表示在张量的第四维进行加权求和运算;α、β表示相应的丰度系数,即相应的权重;为一个高斯随机噪声的三阶张量表示;
其它步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述步骤二中根据高光谱样本数据和窗口大小,分别建立目标样本H1和背景样本H0的空间X、空间Y、光谱和原子的四阶张量矩阵;具体过称为:
步骤二一、从高光谱样本数据库中随机选取高光谱样本数据及设定高光谱样本数据的窗口尺寸;
步骤二二、将每个选取的高光谱样本数据的背景样本转换成三阶张量的形式,然后将所有的三阶张量形式的背景样本组成四阶张量;
步骤二三、将每个选取的高光谱样本数据的目标样本转换成三阶张量的形式,然后将所有的三阶张量形式的目标样本组成四阶张量;
所述三阶张量为空间X、空间Y和光谱;四阶张量为空间X、空间Y、光谱和原子。
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:所述步骤三中依据张量匹配子空间投影算法,求取目标样本H1和背景样本H0的空间X、空间Y、光谱和原子的四阶张量矩阵中空间X、空间Y和光谱三个背景方向和三个目标方向的正交投影矩阵;
待检测信号在三个背景方向正交投影矩阵映射下的数据空间就是背景样本投影子空间,待检测信号在三个目标方向正交投影矩阵映射下的数据空间就是目标样本投影子空间;投影矩阵的具体过程为:
其中,Ui表示四阶背景字典相关的正交投影矩阵,i=1,2,3,4;Vi表示四阶目标背景字典相关的正交投影矩阵,i=1,2,3,4;表示两个张量在i,j,k维的内积运算,i,j,k为信号维度,i,j,k取值为1,2,3,4,为正交投影矩阵大小,I1、I2为空间尺寸,I3为光谱维数;表示目标-背景样本形成的四阶张量子空间,表示背景样本形成的四阶张量子空间。
其它步骤及参数与具体实施方式一或二相同。
具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:所述步骤四中提取待检测信号的空间X、空间Y和光谱的三阶张量矩阵,分别利用步骤三得到的正交投影矩阵将待检测信号映射到步骤三得到的目标样本投影子空间和背景样本投影子空间中;具体过程为:
其中,为待检测信号映射到步骤三得到的背景样本投影子空间中;
为待检测信号映射到步骤三得到的目标-背景样本投影子空间中。
其它步骤及参数与具体实施方式一至三之一相同。
具体实施方式六:本实施方式与具体实施方式一至五之一不同的是:所述步骤五中分别计算待检测信号与步骤四得到的目标样本和背景样本投影子空间投影下的信号的误差,依据张量表达下的广义似然比检测模型和给定的阈值,判定待检测信号是否为检测目标,如果张量表达下的广义似然比检测模型值大于等于η时,为检测目标;否则为背景目标;其中η为阈值;具体过程为:
其中,TD(·)表示待检测信号的检测结果,E(·)表示信号能量,表示Frobenius(弗罗贝尼乌斯范数)范数的平方,η表示设定的阈值。
采用以下实施例验证本发明的有益效果:
实施例一:
实验所用数据是AVIRIS传感器获得的美国圣地亚哥高光谱图像,数据包含224个光谱波段,波长范围为0.4~1.8μm,地面分辨率3.5m,灰度范围为0~10000,图像大小150×150。数据已经经过了大气、几何校正等预处理,并去除了低信噪比和水汽吸收波段,保留了126个波段。图3a为原始数据第4波段图像,图3b为基于张量匹配子空间的高光谱图像目标检测结果图,图3c为高光谱图像匹配子空间检测结果图,用来作为对比实验,图4为两种方法的接收机工作特性曲线。从图3b和图3c的检测结果图和图4所示的与之对应的接收机工作特性曲线可以看出:与对比试验相比,本发明方法能够在低虚警率的情况下达到更好的检测结果,证明了本发明方法的有效性;图4当虚警概率为10-2时,本发明方法的检测概率达到89%,比传统方法提高程度大于10%。
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。