CN107239768A - 一种基于张量主成分分析降维的高光谱图像目标检测方法 - Google Patents
一种基于张量主成分分析降维的高光谱图像目标检测方法 Download PDFInfo
- Publication number
- CN107239768A CN107239768A CN201710433705.6A CN201710433705A CN107239768A CN 107239768 A CN107239768 A CN 107239768A CN 201710433705 A CN201710433705 A CN 201710433705A CN 107239768 A CN107239768 A CN 107239768A
- Authority
- CN
- China
- Prior art keywords
- tensor
- mrow
- block
- detected
- test sample
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/194—Terrestrial scenes using hyperspectral data, i.e. more or other wavelengths than RGB
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
Abstract
一种基于张量主成分分析降维的高光谱图像目标检测方法,本发明涉及高光谱图像目标检测方法。本发明的目的是为了解决现有高光谱图像目标检测方法中未充分考虑在高分条件下空间约束增强的特性,不能从三维数据整体进行信息挖掘,检测精度低的问题。过程为:一:获得三阶目标、三阶背景和待检测的三阶测试样本张量块;二:获得目标、背景和待检测的测试样本三个维度上的投影矩阵;三:将目标、背景和待检测的测试样本投影到预设的张量子空间中;四:计算每一个待检测的测试样本到背景和目标模板的总距离;五:将距离的比值作为灰度值,如果灰度值大于阈值,则确定该中心点的像元为目标,否则认为该中心点的像元为背景。本发明用于图像处理领域。
Description
技术领域
本发明涉及高光谱图像目标检测方法。
背景技术
高光谱图像传感器能够通过上百个光谱通道获取地物的反射辐射信息,其波段范围覆盖了从可见光到近红外乃至长波红外区域。高光谱图像同时包含了地物的空间信息、反射或辐射信息及光谱信息,其特性通常被称为“图谱合一”。而且,高光谱图像数据提供了近乎连续的光谱采样信息,可以记录地物在光谱上很小的反射差异,作为地物分类和检测的依据。研究高光谱图像目标检测新技术,具有重要的理论意义和应用价值。在军事方面,可以揭露敌方目标的伪装、隐藏和欺骗;在民用方面,在公共安全、质量监控、失事点搜寻与营救等方面已经有重要应用。
当前的目标检测手段主要有3种:异常检测、已知样本下的光谱匹配检测和一类分类方式检测。而常用的高光谱图像光谱匹配检测模型一般有三种,欧式距离模型,概率统计模型和子空间模型。典型的方法有正交子空间投影检测方法(Orthogonal SubspaceProjection,OSP),特征子空间投影检测方法(Signature Subspace Projection,SSP),匹配子空间检测方法(Matched Subspace Detector,MSD),光谱匹配滤波检测方法(SpectralMatched filter,SMF),自适应余弦估计检测方法(Adaptive Cosine Estimator,ACE)等。
目前针对高光谱图像的目标检测主要是利用光谱匹配特性,通过光谱匹配程度来判断单点的光谱属性,并未充分考虑在高分条件下空间约束增强的特性,即局部相关性增强的特性。空谱联合检测方法可以提高光谱和空间的利用能力,但一般的空谱联合方法仅仅从空间或光谱的简单组合操作进行分析,不能从三维数据的整体进行信息挖掘。
张量(tensor),即多维数组,更正式地说,一个N阶张量就是N个向量的张量积的结果。一个向量是一阶张量,一个矩阵是二阶张量,三阶或更高阶张量称为高阶张量。利用张量模型对图像进行表示,可以有效地利用图像的内在几何结构信息,大幅提升了遥感图像分类等的性能。目前基于张量表示的遥感图像处理方法已经成功应用到高光谱图像的分类和目标检测当中,并取得了更好的分类和检测结果。
发明内容
本发明的目的是为了解决现有高光谱图像目标检测方法中未充分考虑在高分条件下空间约束增强的特性,不能从三维数据整体进行信息挖掘,检测精度低的问题,而提出一种基于张量主成分分析降维的高光谱图像目标检测方法。
一种基于张量主成分分析降维的高光谱图像目标检测方法具体过程为:
步骤一:对待检测的高光谱图像进行张量块的选择和划分,获得空X-空Y-光谱三阶目标模板张量块、空X-空Y-光谱三阶背景模板张量块和待检测的空X-空Y-光谱三阶测试样本张量块;
步骤二:设定目标模板张量块、背景模板张量块和待检测的测试样本张量块每一维投影后维度的大小,利用所有待检测的测试样本张量块获得目标模板张量块、背景模板张量块和待检测的测试样本张量块三个维度上的投影矩阵;
步骤三:根据步骤二获得的目标模板张量块、背景模板张量块和待检测的测试样本张量块三个维度上的投影矩阵,将步骤一得到的目标模板张量块、背景模板张量块和待检测的测试样本张量块投影到预设的张量子空间中;
步骤四:采用TSAM距离度量方式,计算在投影后的张量子空间内,每一个待检测的测试样本张量块到背景模板张量块的总距离angle_b(m),m=1,2,...M,以及每一个待检测的测试样本张量块到目标模板张量块的总距离angle_t(m),m=1,2,...M;
步骤五:建立张量距离比检测模型,将步骤四中得到的距离angle_b(m)和距离angle_t(m)的比值ratio(m)作为每个待检测的测试样本张量块中心点对应的灰度值,设定阈值η,如果任意一个待检测的测试样本张量块中心点对应的灰度值ratio(m)大于阈值η,则确定该中心点的像元为目标,否则认为该中心点的像元为背景。
本发明的有益效果为:
本发明基于一种张量主成分分析降维的高光谱图像目标检测方法,求取张量块在三个维度上的投影矩阵,通过投影矩阵将所有待检测的张量块以及背景和目标的模板张量块映射到张量子空间中,在投影后的子空间中求取待检测的测试样本张量块到目标模板张量块和背景模板张量块的距离,采用张量距离比模型来刻画其中心点的地物类别属性,最终充分利用高光谱数据的三维整体信息,实现对目标的有效检测。其中,待检测的测试样本张量块中心点地物类别属性指的其中心点的像元是背景还是目标。
张量(tensor)作为对多维数据进行分析的有力数学工具,选用张量模型对高光谱图像进行表示,有效地利用了图像内在的几何结构,能够充分挖掘数据整体的信息。
1、由于本发明方法对张量形式的高光谱图像进行处理,对三阶张量块进行投影、距离度量等操作,相比于传统的基于光谱匹配特性进行目标检测的方法,能够充分考虑在高分条件下空间约束增强的特性,从三维数据整体进行信息挖掘,提高检测精度,实现对高光谱图像中目标的有效检测。
2、本发明方法中,对传统的光谱角映射(Spectral Angle mapping,SAM)方法进行了张量形式的推广,提出用TSAM方法来衡量张量块之间的相似性,用两种距离的比值来刻画测试样本的属性,能够较好地实现目标和背景的划分。
结合图4a、5a和图4b、5b所示的高光谱图像1、2的匹配子空间检测结果和基于张量主成分分析降维的检测结果,以及图4c、5c所示的两种检测方法得到的检测结果对应的ROC曲线,从图4c可得出虚警率为0.1时,现有匹配子空间检测方法发现概率为0.72,本发明张量主成分分析降维检测方法发现概率为0.98;虚警率为0.2时,现有匹配子空间检测方法发现概率为0.82,本发明张量主成分分析降维检测方法发现概率为0.99;从图5c可得出虚警率为0.2时,现有匹配子空间检测方法发现概率为0.56,本发明张量主成分分析降维检测方法发现概率为0.62;虚警率为0.4时,现有匹配子空间检测方法发现概率为0.62,本发明张量主成分分析降维检测方法发现概率为0.95;可得出本发明方法能够在低虚警率的情况下达到好的检测结果,验证了本发明提出的一种基于张量主成分分析降维的高光谱图像目标检测方法的有效性。
附图说明
图1为本发明所述方法流程示意图;
图2a为圣地亚哥局部高光谱图像1第5波段图像示意图;
图2b为圣地亚哥局部高光谱图像1第5波段地物真值图;
图3a圣地亚哥局部高光谱图像2第5波段图像示意图;
图3b圣地亚哥局部高光谱图像2第5波段地物真值图;
图4a圣地亚哥局部高光谱图像1匹配子空间检测结果图;
图4b圣地亚哥局部高光谱图像1基于张量主成分分析降维的检测结果图;
图4c圣地亚哥局部高光谱图像1两种检测方法得到的检测结果对应的ROC曲线,横坐标为虚警率,纵坐标为检测概率;
图5a圣地亚哥局部高光谱图像2匹配子空间检测结果图;
图5b圣地亚哥局部高光谱图像2基于张量主成分分析降维的检测结果图;
图5c圣地亚哥局部高光谱图像2两种检测方法得到的检测结果对应的ROC曲线,横坐标为虚警率,纵坐标为检测概率;
具体实施方式
具体实施方式一:结合图1说明本实施方式,本实施方式的一种基于张量主成分分析降维的高光谱图像目标检测方法,具体过程为:
步骤一:对待检测的高光谱图像进行张量块的选择和划分,获得空X-空Y-光谱三阶目标模板张量块、空X-空Y-光谱三阶背景模板张量块和待检测的空X-空Y-光谱三阶测试样本张量块;
步骤二:设定目标模板张量块、背景模板张量块和待检测的测试样本张量块每一维投影后维度的大小,利用所有待检测的测试样本张量块获得目标模板张量块、背景模板张量块和待检测的测试样本张量块三个维度上的投影矩阵;
步骤三:根据步骤二获得的目标模板张量块、背景模板张量块和待检测的测试样本张量块三个维度上的投影矩阵,将步骤一得到的目标模板张量块、背景模板张量块和待检测的测试样本张量块投影到预设的张量子空间中;
步骤四:采用TSAM距离度量方式,计算在投影后的张量子空间内,每一个待检测的测试样本张量块到背景模板张量块的总距离angle_b(m),m=1,2,...M,以及每一个待检测的测试样本张量块到目标模板张量块的总距离angle_t(m),m=1,2,...M;
步骤五:建立张量距离比检测模型,将步骤四中得到的距离angle_b(m)和距离angle_t(m)的比值ratio(m)作为每个待检测的测试样本张量块中心点对应的灰度值,设定阈值η,如果任意一个待检测的测试样本张量块中心点对应的灰度值ratio(m)大于阈值η,则确定该中心点的像元为目标,否则认为该中心点的像元为背景。
具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤一中对待检测的高光谱图像进行张量块的选择和划分,获得空X-空Y-光谱三阶目标模板张量块、空X-空Y-光谱三阶背景模板张量块和待检测的空X-空Y-光谱三阶测试样本张量块;具体过程为:
给定一个3×3的窗口,将待检测的高光谱图像转换成三阶张量的形式,滑动取样窗口,当取样窗口中心点对应真值图的值为1的时候确定为空X-空Y-光谱三阶目标模板张量块T-tensor(i),i=1,2,...NT;当窗口内所有点对应真值图的值均为0的时候确定为空X-空Y-光谱三阶背景模板张量块B-tensor(j),j=1,2,...NB;滑动过程中一共从张量化的高光谱图像中获得M个待检测的、未知类别属性的空X‐空Y‐光谱三阶测试样本张量块NT、NB和M分别为目标模板张量块、背景模板张量块和待检测的测试样本张量块的个数,取值均为正整数。
其它步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述步骤二中设定目标模板张量块、背景模板张量块和待检测的测试样本张量块每一维投影后维度的大小,利用所有待检测的测试样本张量块获得目标模板张量块、背景模板张量块和待检测的测试样本张量块三个维度上的投影矩阵;具体过程为:
步骤二一、设定目标模板张量块、背景模板张量块和待检测的测试样本张量块每一维度投影后维度的大小Pn,n=1,2,3;
步骤二二、采用张量主成分分析算法利用所有待检测的测试样本张量块获得目标模板张量块、背景模板张量块和待检测的测试样本张量块三个维度上的投影矩阵,其步骤如下:
步骤二二一:预处理,将M个待检测的测试样本张量中心化其中 为M个待检测的测试样本张量的均值;
步骤二二二:初始化,计算取的前Pn个最大的特征值对应的特征向量按列组成
其中,为中心化后的待检测测试样本张量块的n模展开矩阵,m=1,...,M,n=1,...,N,上标T表示转置;是初始化时,用于求取特征向量来得到的矩阵,其上标*表示初始化的过程,用于和后面局部最优化过程中的Φ(n)进行区分;为初始化时得到的第n个方向上的投影矩阵,上标n表示对应所有张量第n维的投影矩阵,下标0表示初始化的过程,即此时迭代次数k=0;
步骤二二三:局部最优化;
【1】:根据步骤二二二得到的将进行投影:
计算的初始离散度:
其中,为中心化后的待检测测试样本张量块经过投影的结果,为初始化时得到的N个方向上的投影矩阵,m=1,...,M,n=1,...,N;×n表示在中心化后的待检测测试样本张量块的第n维进行投影运算,n=1,...,N,表示Frobenius范数的平方,□表示
【2】:
(1)初始化迭代次数k=1,取Φ(n)的前Pn个最大的特征值对应的特征向量按列组成其表达式为:
其中,为第k-1次迭代得到的第n维对应的投影矩阵,上标n表示对应所有张量第n维的投影矩阵,下标k‐1表示第k‐1次迭代;Xm(n)为待检测测试样本张量块的n模展开矩阵,其中m=1,...,M,n=1,...,N,上标T表示转置,为所有Xm(n)的均值,为除了外,其他N-1个投影矩阵的Kronecker积,上标T表示转置,符号表示矩阵的Kronecker积;Φ(n)是迭代过程中,用于求取特征向量来得到的矩阵,Φ(n)和投影后的待检测测试样本张量块总离散度间的关系为:其中trace(□)表示矩阵的迹,□为
(2)计算和
其中,为第k次迭代时得到的第n维对应的投影矩阵,上标n表示对应所有张量第n维的投影矩阵,n=1,...,N,下标k表示第k次迭代;为中心化后的待检测测试样本张量块经过当前第k次迭代得到的矩阵投影的结果;×n表示在待检测测试样本张量的第n维进行投影运算,为当前第k次迭代得到的的总离散度,m=1,...,M, 为的均值,表示Frobenius范数的平方;
(3)重复执行步骤(1)、(2),直至n从1依次取到N,针对本方法中张量块3个方向上投影矩阵的求解,N为张量的维数,N为3;
(4)判断是否成立,如果成立则终止循环,输出当前第k次迭代得到的投影矩阵作为最终的投影矩阵反之,令k=k+1,重复执行步骤(1)~(4),直至迭代次数k=Kmax;
其中,Kmax为预设的最大迭代次数。
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:所述Kmax取100。
其它步骤及参数与具体实施方式一至三之一相同。
具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:所述步骤三中根据步骤二获得的目标模板张量块、背景模板张量块和待检测的测试样本张量块三个维度上的投影矩阵,将步骤一得到的目标模板张量块、背景模板张量块和待检测的测试样本张量块投影到预设的张量子空间中;具体的过程为:
预设的张量子空间的大小对应步骤二一中设定的每一维度投影后维度的大小Pn;
其中,PT-tensor(i)为目标模板张量块T-tensor(i)经三个方向投影矩阵投影后得到的张量子空间里的目标模板张量块;PB-tensor(j)为背景模板张量块B-tensor(j)经三个方向投影矩阵投影后得到的张量子空间里的背景模板张量块;为待检测的测试样本张量块经三个方向投影矩阵投影后得到的张量子空间里的待检测的测试样本张量块; 和分别为空-X、空-Y和光谱维度上实际用于投影的投影矩阵。
其它步骤及参数与具体实施方式一至四之一相同。
具体实施方式六:本实施方式与具体实施方式一至五之一不同的是:所述步骤四中采用TSAM距离度量方式,计算在投影后的张量子空间内,每一个待检测的测试样本张量块到背景模板张量块的总距离angle_b(m),m=1,2,...M,以及每一个待检测的测试样本张量块到目标模板张量块的总距离angle_t(m),m=1,2,...M;具体过程为:
TSAM为光谱角映射(Spectral Angle Mapping)方法的张量形式推广,对于张量和按照TSAM距离度量方式计算和间距离的公式为:
其中,示张量和在(1,2,3)维的内积,示张量在(1,2,3)维同自身的内积,示张量在(1,2,3)维同自身的内积。
angle_b(m)和angle_t(m)的计算公式为:
其中,angle_b(m)为在投影后的张量子空间内,任意一个待检测的测试样本张量块到NB个背景模板张量块PB-tensor(j)的总距离;angle_t(m)为在投影后的张量子空间内,任意一个待检测的测试样本张量块到NT个目标模板张量块PT-tensor(i)的总距离;表示和PB-tensor(j)两个张量在(1,2,3)维的内积,表示同自身两个张量在(1,2,3)维的内积,<PB-tensor(j),PB-tensor(j)>[1,2,3]表示PB-tensor(j)同自身两个张量在(1,2,3)维的内积,表示和PT-tensor(i)两个张量在(1,2,3)维的内积,<PT-tensor(i),PT-tensor(i)>[1,2,3]表示PT-tensor(i)同自身两个张量在(1,2,3)维的内积。
其它步骤及参数与具体实施方式一至五之一相同。
具体实施方式七:本实施方式与具体实施方式一至六之一不同的是:所述步骤五中建立张量距离比检测模型,将步骤四中得到的距离angle_b(m)和距离angle_t(m)的比值ratio(m)作为每个待检测的测试样本张量块中心点对应的灰度值,设定阈值η,如果任意一个待检测的测试样本张量块中心点对应的灰度值ratio(m)大于阈值η,则确定该中心点的像元为目标,否则认为该中心点的像元为背景;具体过程为:
建立张量距离比检测模型,将步骤四中得到的距离angle_b(m)和距离angle_t(m)的比值ratio(m)作为每个待检测的测试样本张量块中心点对应的灰度值,用于衡量每个待检测的测试样本张量块的中心点的类别属性,通过与设定的阈值η相比来得到最终的检测结果,实施的具体过程为:
设定阈值η,如果任意一个待检测的测试样本张量块中心点对应的灰度值ratio(m)大于阈值η,则确定该中心点的像元为目标,否则认为该中心点的像元为背景。
其它步骤及参数与具体实施方式一至六之一相同。
采用以下实施例验证本发明的有益效果:
实施例一:
本例实施一种基于张量主成分分析降维的高光谱图像目标检测方法具体是按照以下步骤制备的,同时采用高光谱图像匹配子空间检测方法来作为对比实验:
本文实验所用的数据是AVIRIS传感器获得的美国圣地亚哥地区的高光谱图像的两个截取部分,其大小为分别为100×100和150×150,数据经过大气、几何校正等预处理,并去除了低信噪比和水汽吸收波段,保留了126个光谱波段,波长范围为0.4~1.8μm,地面分辨率3.5m,灰度范围为0~10000。为了验证本发明方法的有效性,利用高光谱图像匹配子空间检测作为对比实验。图2a为第1幅局部高光谱图像的第5波段图像,图2b为其对应的真值图;图3a为第2幅局部高光谱图像的第5波段图像,图3b为其对应的真值图;图4a为第1幅局部高光谱图像匹配子空间检测结果,图4b为第1幅局部高光谱图像基于张量主成分分析降维的检测结果;图4c为两种检测方法得到的检测结果对应的ROC曲线;从图4c可得出虚警率为0.1时,现有匹配子空间检测方法发现概率为0.72,本发明张量主成分分析降维检测方法发现概率为0.98;虚警率为0.2时,现有匹配子空间检测方法发现概率为0.82,本发明张量主成分分析降维检测方法发现概率为0.99;图5a为第2幅局部高光谱图像匹配子空间检测结果,图5b为第2幅局部高光谱图像基于张量主成分分析降维的检测结果;图5c为两种检测方法得到的检测结果对应的ROC曲线;从图5c可得出虚警率为0.2时,现有匹配子空间检测方法发现概率为0.56,本发明张量主成分分析降维检测方法发现概率为0.62;虚警率为0.4时,现有匹配子空间检测方法发现概率为0.62,本发明张量主成分分析降维检测方法发现概率为0.95;从检测结果可以及ROC曲线可以看出,本发明相比于匹配子空间目标检测,能够在虚警概率相同的情况下有更高的检测概率,证明了本发明方法的有效性。
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (7)
1.一种基于张量主成分分析降维的高光谱图像目标检测方法,其特征在于:所述方法具体过程为:
步骤一:对待检测的高光谱图像进行张量块的选择和划分,获得空X-空Y-光谱三阶目标模板张量块、空X-空Y-光谱三阶背景模板张量块和待检测的空X-空Y-光谱三阶测试样本张量块;
步骤二:设定目标模板张量块、背景模板张量块和待检测的测试样本张量块每一维投影后维度的大小,利用所有待检测的测试样本张量块获得目标模板张量块、背景模板张量块和待检测的测试样本张量块三个维度上的投影矩阵;
步骤三:根据步骤二获得的目标模板张量块、背景模板张量块和待检测的测试样本张量块三个维度上的投影矩阵,将步骤一得到的目标模板张量块、背景模板张量块和待检测的测试样本张量块投影到预设的张量子空间中;
步骤四:采用TSAM距离度量方式,计算在投影后的张量子空间内,每一个待检测的测试样本张量块到背景模板张量块的总距离angle_b(m),m=1,2,...M,以及每一个待检测的测试样本张量块到目标模板张量块的总距离angle_t(m),m=1,2,...M;
步骤五:建立张量距离比检测模型,将步骤四中得到的距离angle_b(m)和距离angle_t(m)的比值ratio(m)作为每个待检测的测试样本张量块中心点对应的灰度值,设定阈值η,如果任意一个待检测的测试样本张量块中心点对应的灰度值ratio(m)大于阈值η,则确定该中心点的像元为目标,否则认为该中心点的像元为背景。
2.根据权利要求1所述一种基于张量主成分分析降维的高光谱图像目标检测方法,其特征在于:所述步骤一中对待检测的高光谱图像进行张量块的选择和划分,获得空X-空Y-光谱三阶目标模板张量块、空X-空Y-光谱三阶背景模板张量块和待检测的空X-空Y-光谱三阶测试样本张量块;具体过程为:
给定一个3×3的窗口,将待检测的高光谱图像转换成三阶张量的形式,滑动取样窗口,当取样窗口中心点对应真值图的值为1的时候确定为空X-空Y-光谱三阶目标模板张量块T-tensor(i),i=1,2,...NT;当窗口内所有点对应真值图的值均为0的时候确定为空X-空Y-光谱三阶背景模板张量块B-tensor(j),j=1,2,...NB;滑动过程中一共从张量化的高光谱图像中获得M个待检测的、未知类别属性的空X‐空Y‐光谱三阶测试样本张量块NT、NB和M分别为目标模板张量块、背景模板张量块和待检测的测试样本张量块的个数,取值均为正整数。
3.根据权利要求2所述一种基于张量主成分分析降维的高光谱图像目标检测方法,其特征在于:所述步骤二中设定目标模板张量块、背景模板张量块和待检测的测试样本张量块每一维投影后维度的大小,利用所有待检测的测试样本张量块获得目标模板张量块、背景模板张量块和待检测的测试样本张量块三个维度上的投影矩阵;具体过程为:
步骤二一、设定目标模板张量块、背景模板张量块和待检测的测试样本张量块每一维度投影后维度的大小Pn,n=1,2,3;
步骤二二、利用所有待检测的测试样本张量块获得目标模板张量块、背景模板张量块和待检测的测试样本张量块三个维度上的投影矩阵,其步骤如下:
步骤二二一:预处理,将M个待检测的测试样本张量中心化其中 为M个待检测的测试样本张量的均值;
步骤二二二:初始化,计算取的前Pn个最大的特征值对应的特征向量按列组成
其中,为中心化后的待检测测试样本张量块的n模展开矩阵,m=1,...,M,n=1,...,N,上标T表示转置;是初始化时,用于求取特征向量来得到的矩阵,为初始化时得到的第n个方向上的投影矩阵,上标n表示对应所有张量第n维的投影矩阵,下标0表示初始化的过程,即此时迭代次数k=0;
步骤二二三:局部最优化;
【1】:根据步骤二二二得到的将进行投影:
计算的初始离散度:
其中,为中心化后的待检测测试样本张量块经过投影的结果,为初始化时得到的N个方向上的投影矩阵,m=1,...,M,n=1,...,N;×n表示在中心化后的待检测测试样本张量块的第n维进行投影运算,n=1,...,N,表示Frobenius范数的平方,□表示
【2】:
(1)初始化迭代次数k=1,取Φ(n)的前Pn个最大的特征值对应的特征向量按列组成其表达式为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>&Phi;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>X</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>&Phi;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>&Phi;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>T</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>X</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>,</mo>
<msub>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>&Phi;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>=</mo>
<mo>(</mo>
<msubsup>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>&CircleTimes;</mo>
<msubsup>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msubsup>
<mo>&CircleTimes;</mo>
<mn>...</mn>
<mo>&CircleTimes;</mo>
<msubsup>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<mo>&CircleTimes;</mo>
<msubsup>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mn>1</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&CircleTimes;</mo>
<msubsup>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>&CircleTimes;</mo>
<mn>...</mn>
<mo>&CircleTimes;</mo>
<msubsup>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,为第k-1次迭代得到的第n维对应的投影矩阵,上标n表示对应所有张量第n维的投影矩阵,下标k‐1表示第k‐1次迭代;Xm(n)为待检测测试样本张量块的n模展开矩阵,其中m=1,...,M,n=1,...,N,上标T表示转置,为所有Xm(n)的均值,为除了外,其他N-1个投影矩阵的Kronecker积,上标T表示转置,符号表示矩阵的Kronecker积;Φ(n)是迭代过程中,用于求取特征向量来得到的矩阵,Φ(n)和投影后的待检测测试样本张量块总离散度间的关系为:其中trace(□)表示矩阵的迹,□为
(2)计算和
其中,为第k次迭代时得到的第n维对应的投影矩阵,上标n表示对应所有张量第n维的投影矩阵,n=1,...,N,下标k表示第k次迭代;为中心化后的待检测测试样本张量块经过当前第k次迭代得到的矩阵投影的结果;×n表示在待检测测试样本张量的第n维进行投影运算,为当前第k次迭代得到的的总离散度,m=1,...,M, 为的均值,表示Frobenius范数的平方;
(3)重复执行步骤(1)、(2),直至n从1依次取到N,N为张量的维数,N为3;
(4)判断是否成立,如果成立则终止循环,输出当前第k次迭代得到的投影矩阵作为最终的投影矩阵反之,令k=k+1,重复执行步骤(1)~(4),直至迭代次数k=Kmax;
其中,Kmax为预设的最大迭代次数。
4.根据权利要求3所述一种基于张量主成分分析降维的高光谱图像目标检测方法,其特征在于:所述Kmax取100。
5.根据权利要求4所述一种基于张量主成分分析降维的高光谱图像目标检测方法,其特征在于:所述步骤三中根据步骤二获得的目标模板张量块、背景模板张量块和待检测的测试样本张量块三个维度上的投影矩阵,将步骤一得到的目标模板张量块、背景模板张量块和待检测的测试样本张量块投影到预设的张量子空间中;具体的过程为:
预设的张量子空间的大小对应步骤二一中设定的每一维度投影后维度的大小Pn;
<mrow>
<mi>P</mi>
<mi>T</mi>
<mo>-</mo>
<mi>t</mi>
<mi>e</mi>
<mi>n</mi>
<mi>s</mi>
<mi>o</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>T</mi>
<mo>-</mo>
<mi>t</mi>
<mi>e</mi>
<mi>n</mi>
<mi>s</mi>
<mi>o</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msub>
<mo>&times;</mo>
<mn>1</mn>
</msub>
<msup>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</msup>
<msub>
<mo>&times;</mo>
<mn>2</mn>
</msub>
<msup>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</msup>
<msub>
<mo>&times;</mo>
<mn>3</mn>
</msub>
<msup>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</msup>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<msub>
<mi>N</mi>
<mi>T</mi>
</msub>
</mrow>
<mrow>
<mi>P</mi>
<mi>B</mi>
<mo>-</mo>
<mi>t</mi>
<mi>e</mi>
<mi>n</mi>
<mi>s</mi>
<mi>o</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>B</mi>
<mo>-</mo>
<mi>t</mi>
<mi>e</mi>
<mi>n</mi>
<mi>s</mi>
<mi>o</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<msub>
<mo>&times;</mo>
<mn>1</mn>
</msub>
<msup>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</msup>
<msub>
<mo>&times;</mo>
<mn>2</mn>
</msub>
<msup>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</msup>
<msub>
<mo>&times;</mo>
<mn>3</mn>
</msub>
<msup>
<mover>
<mi>U</mi>
<mo>~</mo>
</mover>
<msup>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</msup>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<msub>
<mi>N</mi>
<mi>B</mi>
</msub>
</mrow>
其中,PT-tensor(i)为目标模板张量块T-tensor(i)经三个方向投影矩阵投影后得到的张量子空间里的目标模板张量块;PB-tensor(j)为背景模板张量块B-tensor(j)经三个方向投影矩阵投影后得到的张量子空间里的背景模板张量块;为待检测的测试样本张量块经三个方向投影矩阵投影后得到的张量子空间里的待检测的测试样本张量块; 和分别为空-X、空-Y和光谱维度上实际用于投影的投影矩阵。
6.根据权利要求5所述一种基于张量主成分分析降维的高光谱图像目标检测方法,其特征在于:所述步骤四中采用TSAM距离度量方式,计算在投影后的张量子空间内,每一个待检测的测试样本张量块到背景模板张量块的总距离angle_b(m),m=1,2,...M,以及每一个待检测的测试样本张量块到目标模板张量块的总距离angle_t(m),m=1,2,...M;具体过程为:
TSAM为光谱角映射方法的张量形式推广;
angle_b(m)和angle_t(m)的计算公式为:
其中,angle_b(m)为在投影后的张量子空间内,任意一个待检测的测试样本张量块到NB个背景模板张量块PB-tensor(j)的总距离;angle_t(m)为在投影后的张量子空间内,任意一个待检测的测试样本张量块到NT个目标模板张量块PT-tensor(i)的总距离;表示和PB-tensor(j)两个张量在(1,2,3)维的内积,表示同自身两个张量在(1,2,3)维的内积,<PB-tensor(j),PB-tensor(j)>[1,2,3]表示PB-tensor(j)同自身两个张量在(1,2,3)维的内积,表示和PT-tensor(i)两个张量在(1,2,3)维的内积,<PT-tensor(i),PT-tensor(i)>[1,2,3]表示PT-tensor(i)同自身两个张量在(1,2,3)维的内积。
7.根据权利要求6所述一种基于张量主成分分析降维的高光谱图像目标检测方法,其特征在于:所述步骤五中建立张量距离比检测模型,将步骤四中得到的距离angle_b(m)和距离angle_t(m)的比值ratio(m)作为每个待检测的测试样本张量块中心点对应的灰度值,设定阈值η,如果任意一个待检测的测试样本张量块中心点对应的灰度值ratio(m)大于阈值η,则确定该中心点的像元为目标,否则认为该中心点的像元为背景;具体过程为:
建立张量距离比检测模型,将步骤四中得到的距离angle_b(m)和距离angle_t(m)的比值ratio(m)作为每个待检测的测试样本张量块中心点对应的灰度值,用于衡量每个待检测的测试样本张量块的中心点的类别属性,通过与设定的阈值η相比来得到最终的检测结果,实施的具体过程为:
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>t</mi>
<mi>i</mi>
<mi>o</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>a</mi>
<mi>n</mi>
<mi>g</mi>
<mi>l</mi>
<mi>e</mi>
<mo>_</mo>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>a</mi>
<mi>n</mi>
<mi>g</mi>
<mi>l</mi>
<mi>e</mi>
<mo>_</mo>
<mi>t</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mi>M</mi>
</mrow>
设定阈值η,如果任意一个待检测的测试样本张量块中心点对应的灰度值ratio(m)大于阈值η,则确定该中心点的像元为目标,否则认为该中心点的像元为背景。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710433705.6A CN107239768B (zh) | 2017-06-09 | 2017-06-09 | 一种基于张量主成分分析降维的高光谱图像目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710433705.6A CN107239768B (zh) | 2017-06-09 | 2017-06-09 | 一种基于张量主成分分析降维的高光谱图像目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107239768A true CN107239768A (zh) | 2017-10-10 |
CN107239768B CN107239768B (zh) | 2020-10-20 |
Family
ID=59987411
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710433705.6A Active CN107239768B (zh) | 2017-06-09 | 2017-06-09 | 一种基于张量主成分分析降维的高光谱图像目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107239768B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108510013A (zh) * | 2018-07-02 | 2018-09-07 | 电子科技大学 | 基于低秩核心矩阵的改进稳健张量主成分分析方法 |
CN109190506A (zh) * | 2018-08-13 | 2019-01-11 | 北京市遥感信息研究所 | 一种基于核稀疏和空间约束的高光谱目标检测方法 |
CN109447073A (zh) * | 2018-11-08 | 2019-03-08 | 电子科技大学 | 一种基于张量鲁棒主成分分析的红外弱小目标检测方法 |
CN110175581A (zh) * | 2019-05-29 | 2019-08-27 | 哈尔滨工业大学 | 一种基于累积张量分解的高光谱视频图像气体检测方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101807301A (zh) * | 2010-03-17 | 2010-08-18 | 北京航空航天大学 | 一种基于高阶统计量的高光谱图像目标检测方法 |
CN103065160A (zh) * | 2013-01-23 | 2013-04-24 | 西安电子科技大学 | 基于局部协同表示和邻域信息约束的高光谱图像分类方法 |
CN103336968A (zh) * | 2013-05-28 | 2013-10-02 | 中国矿业大学 | 基于张量距离补丁校准的高光谱数据降维方法 |
CN103729854A (zh) * | 2014-01-16 | 2014-04-16 | 重庆邮电大学 | 一种基于张量模型的红外弱小目标检测方法 |
US20160094827A1 (en) * | 2014-09-26 | 2016-03-31 | Xerox Corporation | Hyperspectral imaging devices using hybrid vector and tensor processing |
CN105825200A (zh) * | 2016-03-31 | 2016-08-03 | 西北工业大学 | 基于背景字典学习和结构稀疏表示的高光谱异常目标检测方法 |
CN105913448A (zh) * | 2016-05-25 | 2016-08-31 | 哈尔滨工业大学 | 基于张量匹配子空间的高光谱图像目标检测方法 |
CN106097321A (zh) * | 2016-06-06 | 2016-11-09 | 哈尔滨工业大学 | 一种基于张量表示的偏振高光谱图像目标检测方法 |
-
2017
- 2017-06-09 CN CN201710433705.6A patent/CN107239768B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101807301A (zh) * | 2010-03-17 | 2010-08-18 | 北京航空航天大学 | 一种基于高阶统计量的高光谱图像目标检测方法 |
CN103065160A (zh) * | 2013-01-23 | 2013-04-24 | 西安电子科技大学 | 基于局部协同表示和邻域信息约束的高光谱图像分类方法 |
CN103336968A (zh) * | 2013-05-28 | 2013-10-02 | 中国矿业大学 | 基于张量距离补丁校准的高光谱数据降维方法 |
CN103729854A (zh) * | 2014-01-16 | 2014-04-16 | 重庆邮电大学 | 一种基于张量模型的红外弱小目标检测方法 |
US20160094827A1 (en) * | 2014-09-26 | 2016-03-31 | Xerox Corporation | Hyperspectral imaging devices using hybrid vector and tensor processing |
CN105825200A (zh) * | 2016-03-31 | 2016-08-03 | 西北工业大学 | 基于背景字典学习和结构稀疏表示的高光谱异常目标检测方法 |
CN105913448A (zh) * | 2016-05-25 | 2016-08-31 | 哈尔滨工业大学 | 基于张量匹配子空间的高光谱图像目标检测方法 |
CN106097321A (zh) * | 2016-06-06 | 2016-11-09 | 哈尔滨工业大学 | 一种基于张量表示的偏振高光谱图像目标检测方法 |
Non-Patent Citations (3)
Title |
---|
LEFEI ZHANG等: "《Sparse Transfer Manifold Embedding for Hyperspectral Target Detection》", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
TAO LIN等: "《Small Target Detection Improvement in Hyperspectral Image》", 《INTERNATIONAL CONFERENCE ON ADVANCED CONCEPTS FOR INTELLIGENT VISION SYSTEM》 * |
谷延锋等: "《高分辨率航空遥感高光谱图像稀疏张量目标检测》", 《测绘通报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108510013A (zh) * | 2018-07-02 | 2018-09-07 | 电子科技大学 | 基于低秩核心矩阵的改进稳健张量主成分分析方法 |
CN108510013B (zh) * | 2018-07-02 | 2020-05-12 | 电子科技大学 | 基于低秩核心矩阵的改进稳健张量主成分分析的背景建模方法 |
CN109190506A (zh) * | 2018-08-13 | 2019-01-11 | 北京市遥感信息研究所 | 一种基于核稀疏和空间约束的高光谱目标检测方法 |
CN109447073A (zh) * | 2018-11-08 | 2019-03-08 | 电子科技大学 | 一种基于张量鲁棒主成分分析的红外弱小目标检测方法 |
CN109447073B (zh) * | 2018-11-08 | 2021-11-02 | 电子科技大学 | 一种基于张量鲁棒主成分分析的红外弱小目标检测方法 |
CN110175581A (zh) * | 2019-05-29 | 2019-08-27 | 哈尔滨工业大学 | 一种基于累积张量分解的高光谱视频图像气体检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107239768B (zh) | 2020-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109636742B (zh) | 基于对抗生成网络的sar图像和可见光图像的模式转换方法 | |
Xu et al. | Joint reconstruction and anomaly detection from compressive hyperspectral images using Mahalanobis distance-regularized tensor RPCA | |
Zhao et al. | Hyperspectral anomaly detection via a sparsity score estimation framework | |
Xie et al. | Background learning based on target suppression constraint for hyperspectral target detection | |
CN105427300B (zh) | 一种基于低秩表示和学习字典的高光谱图像异常探测方法 | |
Gu et al. | A selective KPCA algorithm based on high-order statistics for anomaly detection in hyperspectral imagery | |
CN107239768A (zh) | 一种基于张量主成分分析降维的高光谱图像目标检测方法 | |
CN103955926B (zh) | 基于Semi-NMF的遥感图像变化检测方法 | |
CN103456020B (zh) | 基于treelet特征融合的遥感图像变化检测方法 | |
CN107038436A (zh) | 一种基于张量光谱匹配滤波的高光谱图像目标检测方法 | |
CN107194936A (zh) | 基于超像素联合稀疏表示的高光谱图像目标检测方法 | |
CN105631480A (zh) | 一种基于多层卷积网络与数据重组折叠的高光谱数据分类方法 | |
CN106886760B (zh) | 一种基于空谱信息结合的高光谱舰船检测方法 | |
CN101807301A (zh) | 一种基于高阶统计量的高光谱图像目标检测方法 | |
CN107992891A (zh) | 基于光谱矢量分析多光谱遥感图像变化检测方法 | |
CN105913448A (zh) | 基于张量匹配子空间的高光谱图像目标检测方法 | |
CN107316009A (zh) | 基于张量线性判别分析降维的高光谱图像目标检测方法 | |
Xiang et al. | Visual attention and background subtraction with adaptive weight for hyperspectral anomaly detection | |
CN107292258A (zh) | 基于双边加权调制与滤波的高光谱图像低秩表示聚类方法 | |
CN105989597A (zh) | 基于像素选择过程的高光谱图像异常目标检测方法 | |
CN111986162A (zh) | 基于粗定位和协同表示的高光谱异常点快速检测方法 | |
CN110276746A (zh) | 一种鲁棒性遥感图像变化检测方法 | |
CN106548031A (zh) | 一种结构模态参数识别方法 | |
CN115661069A (zh) | 一种高光谱异常检测方法及计算机装置 | |
CN105760857B (zh) | 一种高光谱遥感影像目标探测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |