高光谱图像_高光谱图像分类 - CSDN
精华内容
参与话题
  • 高光谱图像基础知识(1)

    万次阅读 多人点赞 2019-03-15 12:29:53
    本部分主要介绍什么是高光谱图像,成像原理,以及高光谱图像的特点。 高光谱图像 所谓高光谱图像就是在光谱的维度进行了细致的分割,不仅仅是传统的黑,白,或者R、G、B的区别,而是在光谱维度上也有N个通道,例如:...

    本部分主要介绍什么是高光谱图像,成像原理,以及高光谱图像的特点。

    高光谱图像

    所谓高光谱图像就是在光谱的维度进行了细致的分割,不仅仅是传统的黑,白,或者R、G、B的区别,而是在光谱维度上也有N个通道,例如:我们可以把400nm-1000nm分为300个通道,一次,通过高光谱设备获取的是一个数据立方,不仅有图像的信息,并且在光谱维度上进行展开,结果不仅可以获得图像上每个点的光谱数据,还可以获得任意一个谱段的影像信息。
    高光谱成像技术是基于非常多的窄波段的影像数据技术,它将成像技术与光谱技术相结合,探测目标的二维几何空间和光谱信息,获取高分辨率的连续、窄波段的图像数据。

    高光谱成像原理

    以光栅色分光原理为例:

    在经典物理学中,光波穿过狭缝、小孔或者圆盘之类的障碍物时,不同波长的光会发生不同程度的弯散传播,再通过光栅进行衍射分光,形成一条条谱带。也就是说:空间中的一维信息通过镜头和狭缝后,不同波长的光按照不同程度的弯散传播,这一维图像上的每个点,再通过光栅进行衍射分光,形成一个谱带,照射到探测器上,探测器上的每个像素位置和强度表征光谱和强度。一个点对应一个谱段,一条线就对应一个谱面,因此探测器每次成像是空间一条线上的光谱信息,为了获得空间二维图像再通过机械推扫,完成整个平面的图像和光谱数据采集。

    在这里插入图片描述

    数据表现形式

    三维数据中每一个像素,从光谱维这个角度来看,可以看成一种物质在不同波段光谱信号下的表现,如下图中的土壤、水体、植被等。

    分类原理

    不同物质在不同波段光谱信号下的不同表现,可以绘制成一条关于光谱波段和光谱值之间的曲线,根据曲线的差异,我们可以高光谱图像中不同物质进行分类。

    在这里插入图片描述

    展开全文
  • 高光谱数据集

    万次阅读 多人点赞 2018-06-06 17:32:24
    摘自《基于深度卷积神经网络的高光谱图像分类》徐敏1.Indian Pines 数据集Indian Pines 是最早的用于高光谱图像分类的测试数据,由机载可视红外成像光谱仪(AVIRIS)于 1992 年对美国印第安纳州一块印度松树进行...

    摘自《基于深度卷积神经网络的高光谱图像分类》徐敏

    1.Indian Pines 数据集

    Indian  Pines 是最早的用于高光谱图像分类的测试数据,由机载可视红外成像光谱仪(AVIRIS)于 1992 年对美国印第安纳州一块印度松树进行成像,然后截取尺寸为 145×145 的大小进行标注作为高光谱图像分类测试用途。

    AVIRIS 成像光谱仪成像波长范围为 0.4-2.5μm,是在连续的 220 个波段对地物连续成像的,但是由于第 104-108,第 150-163 和第 220 个波段不能被水反射,因此,我们一般使用的是剔除了这 20 个波段后剩下的 200 个波段作为研究的对象。

    该光谱成像仪所成的图像的空间分辨率约为 20m,因此,其较容易产生混合像元,给分类带来难度。

    图 2.1(a)是抽取其中 3 个波段的二维矩阵叠加作为 RGB 图像三个通道所成的伪彩色图像。该数据总共有 21025个像素,但是其中只有 10249 个像素是地物像素,其余 10776 个像素均为背景像素,在实际分类中,这些像素是需要剔除的,由于截取的这块区域均是庄稼,总共有 16类,因此不同的地物具有较为相似的光谱曲线,而且这 16 类中,样本的分布极不均匀

    145*145=21025=10776(背景,黑)+10249(多种地物,彩)  10249=各类地物占得像素个数。


    2.Pavia University 数据集

    Pavia  University 数据是由德国的机载反射光学光谱成像仪(Reflective  Optics Spectrographic Imaging System,ROSIS-03)在 2003 年对意大利的帕维亚城所成的像的一部分高光谱数据。该光谱成像仪对 0.43-0.86μm 波长范围内的 115 个波段连续成像,所成图像的空间分辨率为 1.3m。其中 12 个波段由于受噪声影响被剔除,因此一般使用的是剩下 103 个光谱波段所成的图像。该数据的尺寸为 610×340,因此共包含2207400 个像素,但是其中包含大量的背景像素,包含地物的像素总共只有 42776 个,这些像素中共包含 9 类地物,包括树、沥青道路(Asphalt)、砖块(Bricks)、牧场(Meadows)等。

    3.Salinas数据集

    同 Indian  Pines 图像一样,Salinas 数据也是由 AVIRIS 成像光谱仪所拍摄的,它是对美国加利福尼亚州的 Salinas 山谷所成的像。同 Indian Pines 不同的是,它的空间分辨率达到了 3.7m。该图像原本也有 224 个波段,同样的,我们一般使用的也是剔除了第 108-112,154-167,和第 224 个不能被水反射的波段后剩下的 204 个波段的图像。该图像的尺寸为 512×217,因此共包含 111104 个像素,其中有 56975 个像素是背景像素,可应用于分类的像素有 54129 个,这些像素总共分为 16 类,包括休耕地(Fallow)、芹菜(Celery)等。



    展开全文
  • 高光谱图像处理之目标检测技术一、高光谱图像处理之目标检测1、高光谱图像目标检测的发展趋势和研究现状: 20世纪80年代末,美国的一些研究机构开始利用高光谱图像数据进行目标检测方面的研究。自上世纪九十年代,...

    高光谱图像处理之目标检测技术

    一、高光谱图像处理之目标检测


    1、高光谱图像目标检测的发展趋势和研究现状:

       20世纪80年代末,美国的一些研究机构开始利用高光谱图像数据进行目标检测方面的研究。自上世纪九十年代,国外出现了进行高光谱图像目标检测算法理论研究的研究组。由Reed和Yu提出了基于广义似然比检验的恒虚警RX 检测器(RXD)。Chang课题组提出了基于正交子空间投影的OSP检测方法, Harsanyi提出了基于约束能量最小化的CEM算法。未来高光谱目标检测的发展将会越来越重视实用性,算法的性能将会进一步提高,同时更也加适合使用FPGA硬件对其进行加速从而具有更高的实际应用价值。

    2、高光谱目标检测技术的应用范围:

       高光谱目标检测具有较强的实用性,可应用于公共安全、环境检测、城市规划、食品卫生、地质岩矿的识别等众多方面。


                             图1.1 高光谱图像数据结构

    3、高光谱图像数据的特点:

       高光谱图像数据“图谱合一”,具有丰富的光谱维信息。高光谱图像数据共有三个维度,如图2.1所示,其中,图像空间维信息x、y用于表示物体的实际空间分布,而光谱波段L用于表示每个像素的光谱属性。

    4、高光谱图像目标检测原理:

       高光谱图像的各波段在成像范围内都是连续成像,因此高光谱的光谱曲线一般是平滑的、连续的曲线。高光谱图像的波段L中涵盖了物质的光谱信息,而每种物质的光谱信息都不一样,我们可以利用图像像素的光谱波段L所包含的特定的光谱信息来判断该像素所代表的特定的物质种类。

    5、高光谱图像目标检测流程:


                        图1.2 高光谱图像目标检测流程图

       如图所示,对于拍摄得到的原始高光谱图像数据,需要先对数据进行预处理,包括数据格式化、无用数据剔除以及亮度到反射率的转化等。同时,对于遥感仪拍摄的高光谱图像还需要进行辐射校正,在目标检测前,应对数据进行调整,包括数据归一化等。最后根据已知的先验信息选择相应的目标检测算法进行检测。

    6、现场可编程门阵列(FPGA)在高光谱图像处理中的应用:

       一个FPGA可以大致定义为一系列互连的逻辑块。这些器件的主要优点之一是,为了实现不同的组合或时序逻辑功能,可以根据需要多次重新配置逻辑块及其互连。这一特性为FPGA提供了软件和硬件系统的优势。因为FPGA具有比专用集成电路(ASIC)更多的灵活性和更短的开发时间,更接近GPU提供的性能,但同时功耗要比GPU低得多。FPGA供应商通过改进FPGA体系结构(包括优化的硬件模块)并利用最新的硅技术不断提升FPGA的功耗和能效。可重构性、低功耗的特点以及FPGA对空间电离辐射耐受性的提高,这些因素已经使得FPGA成为目前板载高光谱遥感的最佳选择之一。


    二、算法分类、比较和选择:


    1、方法分类:

       高光谱目标检测方法按照先验信息是否已知分为监督方法和非监督方法。前者用于目标光谱已知的情况下,利用目标光谱与图像像元光谱进行匹配,从而完成目标检测,比如CEM算法、OSP算法;后者多用于异常目标检测,一般不需要目标和背景的先验信息,根据高光谱图像数据获取目标检测所需要的数据,然后根据数据的大小来判断是否为异常目标,比如RXD算法。

    2、CEM、OSP、RXD算法的区别:

    (1)CEM(Constrained Energy Minimization)算法:

       CEM算法主要思想是设计一个FIR线性滤波器,使得在满足约束条件式的情况下滤波器的输出能量最小。该算法不需要图像的背景信息,只需要知道要检测的先验光谱信息(目标向量)即可,具体方法是通过高光谱图像数据和先验已知的待检测目标确定一个滤波向量,让图像经过该滤波向量即可得到检测结果,其中滤波向量的作用是滤除图像中的非目标像素,让感兴趣的目标能够通过,同时抑制由其他信号带来的滤波器输出能量。

    (2)OSP(Orthogonal Subspace Projection)算法:

       OSP算法与CEM算法相比,最大区别在于不仅需要目标的先验知识,还需要图像中背景的先验知识,但在实际中中这些先验信息很难全部得到。在高光谱检测中我们一般用其来检测异常。该算法需要前提条件:图像信息、目标像元、非目标像元(异常目标)信息。

    (3)RXD(Reed-XiaoliDetector)算法:

       RXD算法是异常目标检测领域中最基础的算法,不需要目标光谱的先验知识,而是基于背景服从多元正态分布的假设,通过检测与背景分布中心相比属于异常像元,并在这些感兴趣区域进一步查找可能存在的目标。该算法主要针对的是小目标检测问题。

    (4)确定所采用的目标检测算法:

       由于我们的应用场景多为有特定目标的目标检测,CEM正是针对未知场景中可能存在的特定目标的检测,只需要知道目标的光谱信息即可,而RXD算法适应于对特定场景的异常(未知目标)检测,而OSP算法除了需要已知目标光谱还需要背景信息。综上,我决定采用CEM算法进行高光谱目标检测的实现。


    三、CEM算法分析:


    1、算法步骤

    (1)对高光谱图像进行预处理,得到二维化和归一化后的数据r(L*N);

    (2)根据图像数据r,求得图像的自相关矩阵:

                   

    (3)确定目标光谱向量dd大小为L*1);

    (4)根据公式:


                                       

       设计FIR线性滤波器:


    (5)将归一化后的数据经过FIR滤波器,得到输出信号y

               

    2、问题分解:

       CEM算法的实现可分为三部分:自相关矩阵、矩阵求逆、线性FIR滤波器。


    图 1.3 CEM算法分解流程图

       如图所示,首先根据高光谱图像r求得自相关矩阵,再利用矩阵求逆模块求得自相关矩阵的逆矩阵,结合从光谱库获取到的目标向量的先验信息求得FIR滤波器的滤波向量,最后将高光谱图像r通过FIR滤波器即可得到最终的检测结果。


    四、CEM算法实现:


    1、算法流程:

       在MATLAB和C语言中实现CEM算法的具体流程如图所示,因为语言特性是串行执行命令,所以在编写程序时与硬件设计比较更加直接明了。

    图1.4 CEM算法流程图

    2、数据预处理:

       对前期得到的高光谱图像在MATLAB平台上进行预处理。这一过程主要对原始的200*200*189大小的高光谱图像进行操作:

    (1)二维化:通过调用MATLAB里面的reshape()函数实现。

    (2)归一化:采用“min-max”方法。 具体步骤是先找到数据的最大值
    (maxA)和最小值(minA), 通过 y = (x-minA) /(maxA-minA) 计算的 y 即为归一化后的数据。高光谱数据量大,也造成了数据的存储比较困难,数据在程序中的存储等处理

    3、CEM算法的MATLAB实现:

       首先,根据 2 对高光谱图像数据进行预处理,而后求其自相关矩阵。高光谱图像数据的大小为200*200*189,其中200*200为像素数,189为波段数。

                  

       首先用reshape()函数和transpose()函数将200*200*189的三维高光谱图像数据转为189*40000的二维矩阵,再用自己写的Normalize()函数对矩阵数据用“min-max”方法进行归一化,然后根据公式式利用MATLAB中的矩阵求转置函数transpose()和矩阵相乘操作得到一个189*189的矩阵,最后对189*189的矩阵除以像素数N即可得到自相关矩阵。

    最后,根据先验已知的目标向量和自相关矩阵R的逆矩阵(在MATLAB中矩阵R求逆即为(1/R))求得FIR滤波器,将高光谱图像模拟数据通过FIR滤波器,即可得到最终的检测结果。(本次测试中选取的目标向量 d 是模拟图像中的 C 物质)

    4、CEM算法的C语言实现:

       在Visual Studio平台上完成CEM算法C程序的编写。相比于MATLAB实现,CEM算法的C语言实现主要难点在于矩阵运算。因为在C语言中,矩阵的转置、相乘、求逆等操作均没有现有的函数,需要编写相应的函数。在实现矩阵基本运算过程中,通过动态分配内存运用二维指针传递参数,完成矩阵加减法运算、矩阵转置运算、矩阵相乘运算以及矩阵求逆运算,这样可以节省存储空间,使用完后释放空间即可。

    CEM算法的C语言实现主要包括如下步骤:

    (1)在MATLAB中,将归一化后的高光谱图像数据转为189*40000的二维形式,保存为CEM.mat;

    (2)编写矩阵初始化、矩阵转置、矩阵相乘、矩阵求逆、内存释放函数;

    (3)编写main函数,读CEM.mat,调用上述函数,求得FIR滤波器,进而获得输出信号y,将其写入CEM.txt;

    (4)在MATLAB中显示CEM.txt中数据所代表的图像。

    5、难点解决:

    (1)矩阵转置

       矩阵转置函数的输入是大小为m*n的矩阵A,输出是大小n*m为的矩阵B。采用嵌套的for循环分别遍历矩阵的行和列,将输入的二维矩阵A按列读出,重新按行写入新矩阵B中,即可实现矩阵的转置。

    (2)矩阵相乘

       矩阵相乘函数的输入是大小m*n为的矩阵A和n*k的矩阵B,输出是大小为m*k的矩阵C。采用嵌套的for循环分别遍历矩阵的行和列,首先将矩阵A按行读出,矩阵B按列读出,然后将读出的矩阵A的第i行和矩阵B的第j列对应位相乘求和,将计算的结果写入矩阵C第i行第j列元素中。

    (3)矩阵求逆——QR分解求逆

       QR分解求逆的原理是:对于可逆矩阵A,首先利用QR分解将A矩阵分解为Q矩阵和R矩阵。即A=QR,其中Q是正交矩阵,R是上三角矩阵。然后将公式左右同时求逆,可以求得A的逆矩阵,其中Q矩阵的逆矩阵和转置矩阵相同,R求逆有特定的公式。QR分解求逆的运算较为简单,且数据稳定度较高,可以得到误差较小的求逆结果。

    6、算法实现结果和对比分析:

    1)MATLAB和C语言实现CEM算法的检测结果如图所示:

           

          (a)               (b)            (c)

    图1.6 原图CEM算法检测结果:(a)原图(每一行是一种物质)(b)MATLAB检测C物质(c)C语言检测C物质

     

    通过对MATLAB和C语言实现结果进行对比,发现二者均成功实现了CEM算法并完成了对目标的检测,观察检测结果基本无差异。

    ######### 转载请注明出处:https://blog.csdn.net/Gentleman_Qin/article/details/79992314 #########


    展开全文
  • 高光谱遥感图像相关知识梳理大全

    万次阅读 多人点赞 2017-05-17 17:09:57
    第一部分介绍高光谱图像的一般性原理和知识,第二部分介绍了高光谱图像的噪声问题;第三部分介绍高光谱图像数据冗余问题以及数据降维解决冗余的方法;第四部分介绍高光谱图像的混合像元问题,对光谱解混做了一定介绍...

    前言

    ​ 本资料整理了高光谱遥感图像概念定义、分析处理与分类识别的基本知识。第一部分介绍高光谱图像的一般性原理和知识,第二部分介绍了高光谱图像的噪声问题;第三部分介绍高光谱图像数据冗余问题以及数据降维解决冗余的方法;第四部分介绍高光谱图像的混合像元问题,对光谱解混做了一定介绍;第五部分和第六部分分别介绍了高光谱图像的监督分类和分监督分类的特点、流程和常用算法。

    1.基本介绍

    高光谱遥感(Hyperspectral remote sensing) 是将成像技术和光谱技术相结合的多维信息获取技术,同时探测目标的二维集合空间与一维光谱信息,获取高光谱分辨率的连续、窄波段图像数据。

    高光谱图像高分辨率图像多光谱图像不同。

    高光谱识别优势:

    • 光谱分辨率高、波段众多,能够获取地物几乎连续的光谱特征曲线,并可以根据需要选择或提取特定波段来突出目标特征;

    • 同一空间分辨率下,光谱覆盖范围更宽,能够探测到地物更多对电磁波的响应特征;

    • 波段多,为波段之间的相互校正提供了便利;

    • 定量化的连续光谱曲线数据为地物光谱机理模型引入图像分类提供了条件;

    • 包含丰富的辐射、空间和光谱信息,是多种信息的综合载体。

    高光谱在识别方面的困难:

    1. 数据量大,图像包含几十个到上百个波段,数据量是单波段遥感图像的几百倍;数据存在大量冗余,处理不当,反而会影响分类精度;
    2. 对高光谱图像的分类一方面要求更高的光谱定标和反射率转换的精度,另一方面又因为成像机理复杂,数据量巨大而导致对图像数据预处理困难,包括大气矫正、几何校正、光谱定标和反射率转换等;
    3. 波段多、波段间的相关性高,因此分类需要的训练样本数目大大增多,往往因训练样本不足导致得到的训练参数不可靠(维数灾难);
    4. 针对常规遥感的处理模型和方法不能满足高光谱图像分类的需要。主要问题之一是统计学分类模型的参数估计问题,其对光谱特征的选择要求很高。

    高光谱图像分类中的Hughes 现象:

    Hughes现象:在机器学习问题中,需要在高维特征空间(每个特征都能够取一系列可能值)的有限数据样本中学习一种“自然状态”(可能是无穷分布),要求有相当数量的训练数据含有一些样本组合。给定固定数量的训练样本,其预测能力随着维度的增加而减小。

    ​ 在高光谱遥感图像中,当训练样本数目有限时,分类精度随着图像波段数目的增加先增加,在到达一定极值后,分类精度随这波段数目的增加而下降。

    ​ 传统遥感图像分析是利用图像空间信息,高光谱图像分析的核心是光谱分析。高光谱的遥感数据是一个光谱图像立方体,其最主要的特点将图像空间维与光谱维信息合为一体,与单波段相比,多出了一维光谱信息。在获取地表空间图像同时,会得到每个像元对应的地物光谱信息。

    单波段遥感影像图

    像元光谱特征图

    上图分别是某波段的遥感影像图和某像元的光谱特征图。

    2.高光谱数据噪声处理

    ​ 高光谱图像精细光谱数据能够充分反映地物光谱的细微特征。依据不同地物在光谱特征上的相对差异就可实现地物类别区分,实现目标探测和精细分类。

    ​ 但是,由于成像光谱仪波段通道很密而造成光成像能量不足,故相对全色图而言,高光谱图像的信噪比(SNR)提高比较困难。在图像数据的获取过程中,地物光谱特征在噪声的影响下容易产生“失真”。另外由于高光谱数据量大,在精细分类过程中,往往需要对数据进行降维处理,而在降维过程中需要最大限度保留信号和压缩噪声,所以精确的噪声评估很有必要。另外噪声对精细分类的结果也有最直接影响。所以需要对高光谱数据进行噪声评估。

    ​ 主要的噪声认为是高斯白噪声,分为加性噪声和乘性噪声。

    ​ 噪声评估主要有三种方法:实验室法、暗电流法和图像法。由于前两种方法在实验中难以实现,广泛应用的是第三种——图像法

    图像法分为以下几种:

    • 均匀区域法(HA)

      主要思想:从图像中选取四个以上均匀区域,通过计算这些均匀区域标准差的平均值获取图像噪声的估计值。

      不足:需要人工进行均匀区域选择,无法自动化;满足条件的均匀区域在大部分遥感图中并不存在;子区域噪声估计并不能代表整幅图像的噪声。

    • 地学统计法(GS)

      主要思想:从图像中选择几条均匀的窄条带,通过对这些窄条带的半方差函数的计算实现对图像噪声的估算。一定程度上利用了成像光谱数据的空间相关性。

      不足:与均匀区域法相似。

    • 局部均值与局部标准差法(LMLSD)

      主要思想:假定图像由大量均匀的小块构成,且噪声以加性噪声为主。使用局部均值和局部标准差的概念,将图像分割为很多小块,然后计算这些子块的标准差作为局部噪声大小,并选择包含子块数最多区间的局部标准差的平均值作为整个图像的最佳噪声估计。

      不足:只对高斯白噪声有效,对高斯随机噪声的图像,信号被噪声干扰。

    • 空间/光谱维去相关法(SSDC)

      主要思想:是一种专门针对高光谱图像的噪声评估方法,利用高光谱图像空间维和光谱维存在高相关性的特点,通过多元线性回归去除具有高相关性的信号,利用得到的残差图像对噪声进行估算。

      评价:该方法受地物覆盖类型影响小,并且可以自动执行,是目前较为稳定的高光谱图像噪声评价方法。

    3.高光谱图像数据降维的常用方法

    ​ 高光谱遥感图像所具有的大量光谱波段为地物信息提取提供了极其丰富的信息,有利于更精细的地物分类,然而波段的增多也会导致信息的冗余和数据处理复杂性的提高。

    数据降维满足下面条件:尽可能保留数据的特征信息;去除数据冗余与相关性。

    主要从特征选择特征提取两方面进行降维处理。

    3.1 特征选择

    特征选择是光谱组合,即从原光谱波段数为N的波段中选择M个波段(N>M),可能的光谱组合数为

    N!M!(NM)!

    在这里可以采用Band Index 方法进行光谱波段选择降维。

    *Band Index 方法:*Hyperspectral 遥感图像根据相关性分为K组(如shortwave light、visible light、near-infrared),设每个组中的波段数为 nl(l=1,2,,k)。用pij 表示波段i 与波段j的相关系数 ,σi 表示波段i的方差,Ra表示波段i与不同组其他波段相关系数的绝对值的和,Rw 表示波段i与同组其他波段的相关系数的绝对值的均值,则有

    Rw=1nl|pij|

    从而波段i的Band Index 可表示为:
    Pi=σiRa+Rw

    通过观察,我们知道随着方差增大,波段包含的信息越多;随着相关系数降低,波段的独立性越高。

    Band Index 是一个重要参数,其反映了波段总体包含特征信息和相关性。

    ​ 通过结合Band Index目标物体有效的光谱范围(effiective spectral scope of object),我们可以进行波段选择,进而下一步识别分类。

    3.2 特征提取

    ​ 高光谱的数据降维技术是以图像特征提取为目的,利用低维数据来有效地表达高维数据的特征,同时也压缩了数据量,更有利于信息的快速提取。数据降维包含的内容非常广泛,高光谱遥感图像主要以降低光谱维度提取光谱维度特征为主。
    ​ 上一节的特征选择就讲到了在原始特征空间进行特征选择形成原始空间的一个子空间的特征选择方法,接下来介绍线性变换方法Y=BX ,从高维数据空间中,产生一个合适的低维子空间(不是简单的特征选择组合),使数据在这个空间的分布可以在某种最优意义上描述原来的数据。

    3.2.1 主成分分析

    主成分分析(PCA)是最基本的高光谱数据降维方法,在高光谱数据压缩、去相关、消噪和特征提取中发挥了巨大的作用。PCA 变换又称为霍特林变换(hotelling transform)K-L (karhunen-loeve)变换。变换后的各主成分分量彼此不相关,且随主成分编号的增加该分量包含的信息量减少。

    在高光谱遥感数据的PCA变换中,一般将每个波段当作一个向量来处理,设高光谱遥感数据有p个波段,图像空间维度为m×n,则具体处理流程:

    1. 图像向量化:输入图像数据可以表示成X=(x1,x2,,xp)T,其中 xi 表示为一个 N×1 列向量,这里有N=m×n。即将图像按行或按列展开有规则连接起来,称为一个向量。
    2. 向量中心化:将向量组中的所有向量减去向量组的均值向量,即 Y=XE(X)
    3. 计算向量组 Y 的协方差矩阵 Σ
    4. 求协方差矩阵 Σ的特征值矩阵 Λ 和特征向量矩阵 A
    5. 进行主成分变换, Z=ATY

    PCA变换是基于信息量的一种正交线性变换,变换后的图像信息主要集中在前几个主成分分量中,在变换域中丢弃信息量小的成分分量,经过反变换后仍能得到复原图像的近似图像。

    在PCA变换的基础上提出了分块主成分分析方法定向主成分分析方法(DPCA)选择主成分方法(selective PCA)

    但PCA 变换存在两个明显的缺陷:一是图像数值变换影响明显;二是变换后的信噪比并不一定随着主成分编号的增加而降低。针对这两个问题分别发展了标准化的PCA(standardized PCA, SPCA)残差调整的PCA(residual-scaled PCA, RPCA)

    3.2.2 最大噪声分数变换

    ​ 当噪声方差或噪声在图像各波段分布不均匀时,基于方差最大化的PCA方法并不能保证图像质量随着主成分的增大而降低。所以这里引入最大噪声分数(maximum noise fraction,MNF) 变换 ,该变换根据图像质量排列成分。MNF 方法主要采用SNR和噪声比例来描述图像质量参数。

    假设X=[x1,x2,,xp]Tp×N 矩阵,行向量组的均值向量 E(X)=0,协方差矩阵 D(X)=Σ ,假设

    X=S+N

    其中SN分别指图像中的信号和噪声,且两者不相关。

    ΣSΣN 分别为SN的协方差矩阵。这里假设噪声为加性噪声,则噪声比例可表示为:

    Var{N}/Var{X}

    MNF 变换是一种线性变换,则有
    Zi=aTiX,i=1,,p

    Zi 的噪声比例在所有正交于Zj(j=1,,i1) 的成分中最大,将 ai标准化,
    aTiΣai=1

    由此,MNF变换表示为:
    Z=ATX

    式中,线性变换系数矩阵 A=[a1,a2,,ap]为矩阵 Σ1ΣN的特征向量矩阵,则有
    Σ1ΣNA=ΛA

    式中,对角线矩阵 Λ 为特征值矩阵,第i个元素为特征值 λi,对应成分的噪声比例为

    Var{aTiN}Var{aTiZ}=aTiΣNaiaTiΣai

    MNF 变换最后变换结果的成分按照信噪比的大小排序。

    3.2.3 其他常见变换

    其他常见的变换包含 最小/最大自相关因子分析(minimum/maximum autocorrelation factor, MAF)噪声调整的主成分分析(NPCA)典型相关分析(CCA)独立成分分析(independent compnent analysis, ICA)投影寻踪(projection pursuit, PP)非负矩阵分解非线性主成分分析(Kernel PCA, KPCA)等。

    这些方法都各有其局限性和适用范围,需要根据数据质量和不同的应用需要选择适宜的方法。

    通过对常用的 CPCA、SPCA、MAF 和MNF 的优缺点进行总结,比较在目标探测的高光谱数据降维中的方法,得到下表性能分析:

    性能分析表

    可以由表总结为:

    • CPCA 的优势在于信息损失小,变换后数据结构变化小,但是该方法受数值和噪声影响大;
    • SPCA的优势在于受数值影响小,且在信息保留和数据结构保留保持两方面也不错,但该方法受噪声影响大;
    • MNF 受数值和噪声影响小,同时信息损失小,但变换后数据结构影响很大;
    • MAF 性能最差,不适用于目标探测中的高光谱数据降维

    4.高光谱图像混合像元

    ​ 遥感器所获取的地面反射或发射光谱信号是以 像元(pixel) 为单位记录的。由于高光谱成像光谱仪在获取大量波段时,会导致其每个波段的辐射信号较弱,为了提高信噪比、保证图像质量,就需要保证一定角度的瞬时视场角(IFOV),因此,相比全色和多光谱图像,高光谱图像的空间分辨率低,使得混合像元的问题尤为突出。

    ​ 所以高光谱图像无法利用传统的图像空间信息进行分析的方法,而是利用其上百个波段丰富的光谱信息,从光谱维来挖掘图像的信息。

    4.1 非线性光谱混合模型和线性光谱混合模型

    ​ 物体的混合和物理分布的空间尺度决定了非线性的程度,大尺度的光谱混合完全可以被看作是一种线性混合,二小尺度的内部无核混合是非线性的。在高光谱应用中,利用非线性模型计算出的结果比线性模型计算出的结果要好,但是其需要输入众多的参数,这个实际的应用带来了困难。

    非线性模型

    1. Hapke 混合光谱模型
    2. K-M(KUBELK-MUNK) 混合光谱理论模型
    3. 基于辐射通道密度理论的植被、土壤光谱混合模型
    4. SAIL 模型
    5. PROSPECT 模型和PROSAIL 模型

    线性模型(LSMM)

    当入射光在地物之间不存在多次散射时,在一定IFOV内所形成的混合像元可以通过线性混合模型(LSMM)进行描述。

    4.2 光谱解混

    ​ 解决混合像元问题的过程称为混合像元分解光谱解混,就是根据遥感图像提供的信息判断每个混合像元是由哪些纯像元以怎样的方式混合的。但是,严格意义上的纯像元实际是不存在的,所以在进行解混时,通常是用图像中包含某种比例很高特征地物的像元代替纯像元。这些代替纯像元的近似纯像元,被称为端元

    线性光谱解混是利用LSMM 将遥感图像X 中每个混合像元分解成其端元和对应丰度,从而得到端元矩阵 E丰度矩阵(即所占比例) A 的过程。即

    XL×n=EL×mAm×n+ϵL×n

    基本工作分两步:端元提取,即确定E,丰度反演,即确定 A

    1. 确定端元数目,在少数波段的多光谱图像中可使用PCA 根据协方差特征值的大小确定;在高光谱中使用Neyman-Pearson 探测理论的特征阈值分析分析方法确定,(确定数目中还需要对噪声进行评估,可采用光谱维去相关(SSDC)法);
    2. 端元提取:MaxV 方法直接用于端元提取,N-FINDER 方法需先降维再提取端元,采用降维时,提取效果:MaxV>MNF>MAF>CPCA>SPCA
    3. 端元光谱变异性:由于非线性因素的影响,端元的光谱存在变异性,即同一类型存在多种端元,这也是线性光谱解混误差的主要因素之一。混合像元以及端元的光谱变异性,都是对客观事物的不确定性描述,可用模糊子集和模糊测度的思想来解决;

    常用的端元提取方法还有:

    • 纯像元指数(pure pixel index, PPI) 算法;
    • 改进的快速迭代PPI(faster iterative PPI, FIPPI) 算法;
    • 内部最大体积(N-FINDER)算法;
    • 顶点成分分析(vertex component analysis,VCA) 算法;
    • 基于最大距离(maximum distance, MaxD) 算法;
    • 基于最大体积(maximum volume, MaxV)算法;
    • 定量化独立成分(ICA)及其衍生算法等。

    当得到端元矩阵 E 后,就要通过丰度反演求解高光谱图像中每个像元里各个端元所占比例,即求 丰度矩阵A 的过程。

    常用的方法有:

    • 最小二乘法(LS):
      • 无约束(UCLS)
      • “和为1”约束(SCLS)
      • “非负”约束(NCLS)
      • 全约束(FCLS)
    • 滤波向量(FV)法
    • 迭代光谱混合分析(ISMA)法等

    5.高光谱遥感数据监督分类

    高光谱图像数据将地物光谱信息和图像信息融为一体,其数据具有几何空间、光谱特征空间两类表达方式。

    几何空间:直观表达每个像元在图像中的空间位置以及它周边像元之间的相互关系,为高光谱图像处理提供空间信息。

    光谱特征空间: 高光谱遥感图像每个像元对应着多个成像波段的反射值。近似连续的光谱曲线表示为一个N维向量,向量在不同波段值的变化反映了其代表的目标的辐射光谱信息,其优势是特征维度的变化和扩展性。我们处理将高光谱像元作为高维特征空间的数据点,根据数据的统计特性来建立分类模型,但是它的弱点是无法表达像元间的几何位置关系。

    监督分类常用于高光谱图像数据的定量分析,其主要流程是:首先,利用分类器对已知类别机器对应的训练样本进行学习,获取各图像上各类别像元的分类特征;然后,选择适当的分类判据,根据分类的决策准则进行分类。

    一般按以下步骤进行:

    1. 高光谱数据选择。
    2. 图像的预处理。即几何配准、校正等,确保获取正确光谱和几何信息。
    3. 确定地物种类。即根据提取的训练数据特征确定分类类别。
    4. 选择训练样本。
    5. 特征提取和特征选择。选择各类地物可区分性最强的特征从而提高分类精度。
    6. 选择合适的分类方法进行分类。
    7. 分类后处理。传统基于像元分类方法分类后可能存在大量噪声及孤立像元,根据地物的连续性,利用主成分滤波等方法减少该因素影响,从蹄提高分类精度。
    8. 分类结果评价。根据已知类别的测试数据类别与分类结果比较,确认分类的精度与可靠性。

    5.1 基于光谱特征空间分类

    ​ 基于光谱特征空间的高光谱图像分类方法主要是建立在对高光谱图像光谱特征提取和变换的基础上,分为两种思路:一种是基于地物物理光学性质的光谱曲线来进行地物识别,代表性方法是光谱特征匹配方法;另一种是基于特征空间的分类方法,主要利用数据的统计特征来建立分类模型,主要方法有传统的遥感分类方法以及神经网络,支持向量机等复杂的图像分类方法。

    5.1.1 光谱特征匹配分类方法

    基于光谱曲线进行分类识别是利用光谱库中已知的光谱数据,采用匹配的算法来识别图像中地面覆盖类型。这种匹配既可以是全谱段范围内比较,也可以是感兴趣波段的光谱比较。

    基于光谱间最小距离的匹配算法是在计算未知光谱与参考光谱距离后,再根据最最小二乘原则进行匹配的分类方法。该方法对噪声敏感,所以匹配前需要去噪预处理。

    光谱角度填图(spectral angle mapping, SAM)方法 把光谱看作是光谱空间的多维矢量,计算两光谱向量的广义夹角,夹角越小,光谱越相似,最终根据相似性阈值对未知像元光谱进行分类。SAM 方法的显著特征是夹角值与光谱向量模无关,只比较光谱在形状上的相似性,这也是SAM方法与前者的区别与优势。

    高光谱图像的光谱波段间隐含这特定的物理含义,光谱曲线的形状特征是地物内在物理化学性质的外在反映。

    5.1.2 遥感图像统计模型分类方法

    基于统计的极大似然分类是传统遥感图像分类中应用最广泛的分类方法,在极大似然分类器中,首先为简化分类过程,假设高光谱图像每类地物的概率密度都服从多维正态分布,然后利用分类训练样本分别对统计参数进行估计获得其概率密度函数,最后利用Bayes 公式将最大似然概率公式变为

    p(x|wi)p(wi)p(x|wj)p(wj)

    式中, p(wi)p(wj) 分别是类i与类j的先验概率分布,实际中由于无法估计,假设其相等。所以在实际分类过程中,将各像元划分到其属于图像上不同类别的概率中较大的那一类中去。

    不足:要想获得好的效果,统计参数的估计变得十分重要;为了获得可靠的参数,每个类别必须有足够的训练样本,这对于上百个波段的高光谱图像是很困难的。

    5.1.3 高光谱图像神经网络的方法

    后面结合深度学习最新知识再详细介绍。

    不足:高光谱图像分类时,经常遇到“同物异谱”想象,这使得神经网络的分类算法难于收敛,严重降低分类精度。

    5.1.4 高光谱图像支持向量机的方法

    支持向量机方法用于高光谱图像分类的优点是能够直接对高维数据处理,不必经过降维处理,而采用全部波段进行分类,保证了光谱信息的充分性。

    原始SVM 算法是二分类器,在高光谱图像的多分类以及精细分类中,采用多个SVM 组合或级联的方式实现多分类。

    5.2 几何空间与光谱特征空间结合的高光谱图像分类

    基于光谱特征空间的高光谱图像分类方法是将高光谱数据作为一种无序的待聚类数据集合,没考虑像元点的几何空间特征,综合高光谱图像的几何空间特征和光谱特征空间进行分类,能够有效提高分类精度。

    5.2.1 基于同质地物提取的高光谱图像分类方法(ECHO)

    前提:图像空间分辨率较高(优于5m),常用实验数据集如AVIRIS Data 大约20m左右

    主要思想:首先将图像划分为不同的图像对象,图像对象为形状与光谱特征具有相似性的同质区域;对于没有划入同质区域的像元利用极大似然分类器对它们进行分类最终获取分类结果。

    5.2.2 纹理信息辅助下的高光谱图像分类

    利用灰度共生矩阵进行纹理统计。

    5.3 面向对象的高光谱遥感图像分类

    ​ 面向对象分类的特点即分类的最基本对象从像元转换到图像对象。分类的核心是高光谱图像的分割,在这个阶段应该结合光谱信息空间信息

    面向对象的高光谱图像分类框架

    上图是面向对象的高光谱图像分类框架。

    5.4 高光谱图像的分类精度评价

    5.4.1 误差矩阵

    误差矩阵的主体是一个 k×k 的方阵,其中列为地面参考验证信息,行为分类结果,矩阵的对角元表示被分到正确类别的像元个数,对角线以外的元素表示错分的误差。

    其中,生产者精度(PA)

    PAi=xi,ix+i

    用户精度(UA)
    UAi=xi,ixi+

    5.4.2 漏分误差和多分误差

    漏分误差(OE)是指类别i在误差矩阵中,有多少被错误分到了其他类别,类别i的漏分误差

    OEi=1xi,ix+i

    多分误差(CE)是将其他类别像元错误划分到i类中,某类别i的多分误差为
    CEi=1xi,ixi+

    5.4.3 Kappa 分析

    Kappa 分析是一种定量评价遥感分类图与参考数据之间一致性或精度的方法,能识别整体图像的分类误差性。

    Kappa=1

    计算得到
    Kappa=Nki=1xi,iki=1(xi+×x+i)N2ki=1(xi+×x+i)

    6.高光谱图像非监督分类

    高光谱图像的非监督分类算法还是采用传统遥感中的非监督方法,如 K-means 算法ISODATA 算法FCM 算法等算法以及它们的优化算法。

    高光谱图像的非监督分类使用过程中需要注意这几个问题:

    1. 高光谱图像的波段众多,在自学习中会产生大量冗余,分类前进行特征提取和选择有利于提高自学习过程效率;
    2. 可以充分利用高光谱图像中端元提取、光谱分析等手段获得较精确的类别初值,辅助非监督分类;
    3. 选择非相似性度量时,需要充分考虑地物光谱特征,适当选择;
    4. 利用搜索算法求解非监督分类目标函数时,尝试利用全局搜索算法代替局部搜索算法提高非监督分类精度。
    展开全文
  • 高光谱/遥感图像常用数据集

    热门讨论 2020-07-30 23:33:23
    内有常用的高光谱图像(HSI)数据集,有常用的Indian,Pavia。每类数据集里包含图像原始信息,以及相对应的地面真实数据的类别标签。因文件大小限制,只能传这两个数据集,还有Salinas等数据集,大家需要的话可留言。
  • 高光谱图像处理

    热门讨论 2020-07-29 14:20:53
    高光谱图像处理 matlab,用matlab处理高光谱图像
  • 2 全色图像、多光谱图像高光谱图像 2.1 全色图像 2.2 多光谱图像 2.3 高光谱图像 参考资料 1 光谱分辨率、空间分辨率、时间分辨率 遥感(Remote Sensing),可以理解为遥远的感知。遥感技术利用搭载在遥感平台...
  • 高光谱图像是什么?

    千次阅读 2017-12-28 20:53:47
    10 的幂,eg.10和1000差了100倍,相差两个数量级)范围内的光谱图像称为高光谱图像(Hyperspectral Image)。遥感技术经过20世纪后半叶的发展,无论在理论上、技术上和应用上均发生了重大的变化。其中,高光谱图像...
  • 高光谱图像分类(一)入门

    万次阅读 多人点赞 2016-10-18 14:04:35
    一、什么是高光谱图像?如何获取?高光谱图像并非普通的二维图像,而是一种三维数据。如下图所示: 1.成像原理:星载光谱成像仪会在沿着轨道飞信时,对地表发出不同波段的光谱信号。而地表不同物质对每一个波段的...
  • 参考: https://www.zhihu.com/question/54688096 http://blog.sina.com.cn/s/blog_15183f5750102vxgm.html http://www.docin.com/p-980446793-f4.html ...
  • 高光谱图像选择波段的研究

    万次阅读 多人点赞 2016-08-25 16:20:39
    高光谱波段选择
  • 高光谱图像分类(三)分类流程

    万次阅读 多人点赞 2017-04-11 19:04:36
    如何利用稀疏表示进行高光谱图像分类呢?前面我们已经了解了高光谱图像分类的一些基本概念,那这篇文章当中将讲解高光谱图像分类具体的流程是怎么样的。以下是高光谱图像分类的具体详细步骤: 1.导入indian_pines...
  • 高光谱遥感图像的读取---使用matlab

    千次阅读 2015-05-25 10:24:34
    函数multibandread  im_hyper = multibandread(filename, size, precision, offset, interleave, byteorder) ; 注意:precision 指数据的存储类型比如:hydice 的urban 为: 'uint16'  offset 一般为 :0 ...
  • 高光谱遥感图像数据下载

    千次阅读 2018-09-14 13:37:19
    包括常用的Indian Pines,Pavia Centre and University等的 .mat 数据:下载地址
  • 是MAT-file格式,不是M,高光谱图像是matlab代码,怎么打开呢?
  • 一、影像剖面线  X轴(水平)、y轴(竖直...主图像窗口,tools-profiles-xprofile  图像波谱图(z)  主图像窗口,tools-profiles-z profile 在spectral profile中 options-collect spectra 绘制多条 二、图像像元
  • 全色和多光谱影像 1.ENVI工具栏Spectral----Gram-Schmidt Spectral Sharpening ...3.然后选择分辨率的正射图像 4.保存图像 选择要保存的融合后的图像右击 save selected fi...
  • 论文地址基于深度卷积神经网络的光谱遥感图像分类西华大学学报同时利用高光谱图像的光谱信息和空间信息的深度卷积神经网路分类模型。基于深度学习到的深度特征,用逻辑回归分类器进行分类训练。高光谱图像面临的...
  • 用于高光谱图像分类的深度卷积神经网络     摘要:近来,卷积神经网络已经在各种视觉任务中表现出优异的性能,包括常见二维图像的分类。 在本文中,采用深卷积神经网络直接在光谱域中对高光谱图像进行分类。更...
  • 高光谱图像基于MATLAB的PCA降维

    千次阅读 2018-09-17 13:45:54
    工具:matlab,本人使用的是2016a 使用数据集用公开的Salinas数据集为例,Salinas数据集为512*217*204,即它有204个波段,我们要把它从204维降至3维。 matlab中内含了进行PCA降维的函数,但这个函数输进去的数据...
1 2 3 4 5 ... 20
收藏数 8,747
精华内容 3,498
关键字:

高光谱图像