CN100514085C - 高光谱图像的空间-光谱信息协同提高分辨率的方法 - Google Patents
高光谱图像的空间-光谱信息协同提高分辨率的方法 Download PDFInfo
- Publication number
- CN100514085C CN100514085C CNB2007101444495A CN200710144449A CN100514085C CN 100514085 C CN100514085 C CN 100514085C CN B2007101444495 A CNB2007101444495 A CN B2007101444495A CN 200710144449 A CN200710144449 A CN 200710144449A CN 100514085 C CN100514085 C CN 100514085C
- Authority
- CN
- China
- Prior art keywords
- information
- spectral
- resolution
- pixel
- spatial
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
高光谱图像的空间-光谱信息协同提高分辨率的方法,它涉及利用高光谱图像信息提高空间分辨率的方法。它解决了现有高光谱图像处理中存在的不能充分利用空间信息和光谱信息以改善图像分辨率的问题。本发明步骤为:一、输入高光谱图像数据;A、空间信息提取;A一、特征波段选择;A二、空间局部分析及判断;B、光谱信息提取;B一、光谱端元提取;B二、混合像元分解;C、空-谱协同超分辨;D、得到分辨率提高的高光谱图像。本发明突破图像获取时的空间分辨率极限;利用支持向量机解混和局域空间相关性进行空间-光谱信息协同技术提高高光谱图像的空间分辨率,可极大提高目标的探测和定位能力,突破图像获取手段的限制,弥补硬件的不足。
Description
技术领域
本发明涉及一种利用高光谱图像信息提高空间分辨率的方法。
背景技术
空间分辨率简单来说就是成像系统对图像细节分辨能力的一种度量,也是图像中目标细微程度的指标,它表示景物信息详细程度,是评价传感器性能和遥感信息的重要指标之一,也是识别地物目标形状大小的重要依据。遥感图像的空间分辨率的高低与成像光学系统有着直接的关系,如果其分辨率较低,将使得遥感图像中存在较多的混合像元,严重影响图像的分析和理解,这对于目标分类、检测和识别来说,是非常不利的。
光谱分辨率是指传感器在一定波长范围内对地物光谱进行离散采样的精细程度。光谱分辨率是表征传感器获取地物光谱信息性能的主要指标。相对于空间图像信息,作为刻划地物特征的另一种方式,通过远程探测得到的光谱信息同样可实现对地物的辨识,并且光谱信息直接与目标的物质组成有关,特别是对于目标识别、植被的精细分类、海洋水色定量监测以及军事上对伪装的辨识等从光谱的角度比空间的图像更适合。
高光谱遥感(Hyperspectral Remote Sensing)是在20世纪80年代发展起来的一种全新的遥感技术。该技术利用星载或机载的成像光谱仪设备对地面进行成像,成像光谱仪在对目标的空间特征成像的同时,对每个空间像元经过色散形成几十个乃至几百个窄波段以进行连续的光谱覆盖,从而形成谱分辨率达到纳米数量级的遥感数据。这种数据由于光谱分辨率高,通常称为高光谱数据或高光谱图像。高光谱数据的光谱分辨率在10纳米左右,比多光谱(Multispectral)图像高出几十甚至上百倍。伴随着成像光谱技术的不断发展,高光谱数据已经被应用到了众多领域中。从民用领域广泛应用的环境监测、城市规划、农作物估产、洪涝灾害调查、国土资源调查,到军事领域的卫星侦察、目标检测识别等等。
高光谱图像的突出特点是在获得目标图像二维空间景像信息的同时,还可以获得高分辨率的一维表征其物理属性的光谱信息,即图谱合一。通过处理高光谱图像中目标图像的空间特征和光谱特征,可以以较高的可信度辨别和区分地物目标。这对遥感图像军事侦察、真/假目标识别、农林的精细化分类等都具有重要应用意义和巨大的潜力。遥感技术长期以来的两个主要发展趋势就是向高空间分辨率和高光谱分辨率方向发展,但二者的发展往往是相互矛盾、相互制约的,这主要是由于成像光学系统在设计和实现上的限制。高光谱图像的光谱分辨率一般较高,但其空间分辨率却偏低,这对于目标识别算法较为不利。在遥感技术快速发展的今天,我们对遥感图像的分辨率有着越来越高的要求,但对于现有的成像设备,由于其制作工艺和现有技术的约束,还远远不能满足各方面的要求。因此,运用软件方法来有效提高高光谱成像系统的空间分辨率具有重要的实用价值和现实意义。
现有方法是单纯从高光谱图像的空间几何关系出发,通过空间域或者频率域的处理技术来进行空间超分辨。然而在此过程中,很少利用到光谱信息,使得这类算法还局限于原有的一般数字图像超分辨,没有充分发挥高光谱的优势。或者只考虑光谱信息,获得混合像元中各组成成分的比例,而很少利用空间信息进行进一步细化处理。总之,以往的算法没有将光谱和空间信息协同利用,所以高光谱的特点没有充分得以体现,存在信息没有充分利用的不足。
发明内容
本发明的目的是为了解决现有高光谱图像处理技术中存在的不能充分利用空间信息和光谱信息以改善图像分辨率的问题,进而提出了一种高光谱图像的空间-光谱信息协同提高分辨率的方法。
本发明的步骤如下:
步骤一:输入高光谱图像数据:根据数据格式读入高光谱图像数据,将所述高光谱图像数据同时输入到计算机的步骤A中的空间信息提取模块和步骤B中的光谱信息提取模块中进行信息提取;
步骤二:本步骤包括步骤A和步骤B,
步骤A:空间信息提取:通过空间信息模块进行特征波段选择和空间局部分析及判断;
步骤A一:特征波段选择:特征波段选择是对高光谱图像数据信息量大的特征谱段进行提取或选择;
步骤A二:空间局部分析及判断:利用微分算子边缘检测算法提取步骤A一中所提取或选择的特征谱段的边缘信息,通过边缘信息进行边缘像元判断;是非边缘像元,利用复制技术进行插值处理,得到空间类别信息和相互关系信息,输入到步骤三中;是边缘像元,即为混合像元所在,获取混合像元,将混合像元数据输入到步骤B二中;
步骤B:光谱信息提取:通过光谱信息提取模块进行光谱端元提取和混合像元分解;
步骤B一:光谱端元提取:利用PPI或者N-FINDR算法提取高光谱图像中的光谱端元,记录端元的光谱信息和位置信息;
步骤B二:混合像元分解:用得到的端元的光谱信息和位置信息通过支持向量机方法,从混合像元中获取各组成成分的比例;
步骤三:空-谱协同超分辨:对步骤B二中得到的混合象元解混后的各组成成分的比例分布和步骤A二中得到的非边缘像元所取得的空间类别信息和相互关系信息,按照相关性准则进行分布;
步骤四:得到分辨率提高的高光谱图像。
本发明与现有技术相比具有如下优点:
(1)克服通常方法单纯从高光谱图像的空间几何关系出发进行空间超分辨,或者只考虑光谱信息进行混合像元的分解来获得子像元构成,有效解决了现有算法信息利用率低的问题,可突破图像获取时的空间分辨率极限;
(2)利用支持向量机解混和局域空间相关性进行空间-光谱信息协同技术提高高光谱图像的空间分辨率,可极大提高目标的探测和定位能力,突破图像获取手段的限制,弥补硬件的不足。
附图说明
图1是本发明的流程图;图2是混合像元与其他相邻像元的示意图;图3是解混分量图;图4是解混分量图所对应的第一种空间分布关系图;图5是解混分量图所对应的第二种空间分布关系图。
具体实施方式
具体实施方式一:结合图1说明本实施方式,本实施方式的步骤如下:
步骤一:输入高光谱图像数据:根据数据格式读入高光谱图像数据,将所述高光谱图像数据同时输入到计算机的步骤A中的空间信息提取模块和步骤B中的光谱信息提取模块中进行信息提取;
步骤二:本步骤包括步骤A和步骤B,
步骤A:空间信息提取:通过空间信息模块进行特征波段选择和空间局部分析及判断;
步骤A一:特征波段选择:特征波段选择是对高光谱图像数据信息量大的特征谱段进行提取或选择;
步骤A二:空间局部分析及判断:利用微分算子边缘检测算法提取步骤A一中所提取或选择的特征谱段的边缘信息,通过边缘信息进行边缘像元判断;是非边缘像元,利用复制技术进行插值处理,得到空间类别信息和相互关系信息,输入到步骤三中;是边缘像元,即为混合像元所在,获取混合像元,将混合像元数据输入到步骤B二中;
步骤B:光谱信息提取:通过光谱信息提取模块进行光谱端元提取和混合像元分解;
步骤B一:光谱端元提取:利用PPI或者N-FINDR算法提取高光谱图像中的光谱端元,记录端元的光谱信息和位置信息;
步骤B二:混合像元分解:用得到的端元的光谱信息和位置信息通过支持向量机方法,从混合像元中获取各组成成分的比例;
步骤三:空-谱协同超分辨:对步骤B二中得到的混合象元解混后的各组成成分的比例分布和步骤A二中得到的非边缘像元所取得的空间类别信息和相互关系信息,按照相关性准则进行分布;
步骤四:得到分辨率提高的高光谱图像。
具体实施方式二:本实施方式与具体实施方式一不同点在于步骤A一中的特征波段选择采用的方法有主成分变换提取特征波段或者基于信息熵最大原则选择特征波段;其它组成和连接方式与具体实施方式一相同。
具体实施方式三:结合图2说明本实施方式,本实施方式与具体实施方式一不同点在于步骤A二中的空间局部分析及判断是分析每一个特征谱段的高光谱图像,其中xi,j为单波段高光谱图像中待处理的像元,W为以该像元为中心的局域窗,以3×3的窗为例加以说明:
在局域窗W内应用边缘检测算法进行边缘像元判断;xi,j不是边缘像元,则将其直接进行像元复制处理;xi,j是边缘像元,则输入数据到步骤B二进行混合像元分解。
式中,x为窗口内的像元灰度均值。具体地说,对于每一个待插值像元点xi,j,判断它是否为边缘点的方法为:计算像元xi,j所在的邻域窗W内像元的标准方差当且仅当的值超过预先设定的阈值时,才认为xi,j为边缘像元点,从而得到准确的边缘像元点。其它组成和连接方式与具体实施方式一相同。
具体实施方式四:本实施方式与具体实施方式一不同点在于步骤B二中的混合像元分解采用支持向量机方法,线性光谱混合模型可以表示为
其中x为观测光谱向量,s1,s2,...,sM为M个线性无关的光谱端元,a1,a2,...,aM为响应的混合比例,称为“丰度”,w代表加性噪声。通常,丰度a1,a2,...,aM满足归一和非负的约束,即
在两类混合的情况下,同时考虑到归一和非负的约束,混合模型可以改写为
x=a1s1+a2s2+w=as1+(1-a)s2+w (5)
通常我们假设w为高斯白噪声,此时x的最小均方估计可以写作
支持向量机的输出可以写为
其中xi为支持向量,它位于分类边界附近;yi为相应的支持向量xi所属类别,值取+1或-1;αi为拉格朗日乘子;K(x,xi)表示输入向量x与支持向量xi的核函数输出;b为偏置。
公式(7)表明,对于输入向量x,其支持向量机输出为K(x,xi)的加权和,且这些权值是已知的。在利用支持向量机进行光谱解译的过程中,需要寻找支持向量机输出g(x)与所求丰度之间的关系。进一步的,我们用x的最小均方估计代替x,同时考虑两类情况,那么目前的主要问题即为寻找与的关系。
根据公式(6),我们得到
如果两个端元与的核函数与支持向量xi的核函数由下式表示,分别对应于 和
那么公式(9)中的待定系数就可以表示为
其它组成和连接方式与具体实施方式一相同。
具体实施方式五:结合图3、图4和图5说明本实施方式,本实施方式与具体实施方式一不同点在于步骤三,在与图3光谱解混分量图相对应的图4和图5中的高分辨率图像窗中的两种空间分布关系中,图5体现出更强的空间相关性,因此成为更加合理的超分辨结果。BP神经网络可以实现输入空间到输出空间的非线性映射,基于此,利用BP网络训练输入输出之间的函数关系,即用该模型表达解混分量窗中各像元值与其中心像元内部的子像元空间分布之间的函数关系。设窗口W包含3×3个像元,xij为中心像元,放大因子为z=2,该函数关系可以用下式表达:
其中
该模型是建立于空间相关性假设基础上的,即邻近的像元值比远离的像元值更接近。子像元值可由对应的像元值和其邻域像元值确定。邻近的像元值越大,子像元被置为目标的可能性就越大,这使得空间相关性达到最大值。BP网络通过训练可以得到高分辨图像的空间分布规律并将之存储于连接权中。
表达式(12)给出了一对学习模式。解混分量图窗口作为输入模式,以向量X=(xi-1,j-1,xi-1,j,xi-1,j+1,xi,j-1,xi,j,xi,j+1,xi+1,j-1,xi+1,j,xi+1,j+1)T的形式输入到网络。高分辨图像中的一个子像元作为输出模式,以 的形式输出。因此,表达式(12)需要建立4个BP网络模型。可以看出BP网络输入节点的个数仅由解混分量图窗口大小决定。
学习模式被确定以后,网络开始训练。需要大量样本训练网络,并且这些样本均满足空间相关性原则。设一个子像元的邻域内包含N个像元,空间相关性可由下面的公式来度量:
其中权值wk可定义为子像元到第k个像元中心的距离的倒数或该距离平方的倒数,αk为由步骤B二获得的解混分量,通过这个度量标准可以分析出不合理的样本,并将之淘汰,由此逐像元操作,即可获得分辨率提高的超分辨图像。其它组成和连接方式与具体实施方式一相同。
Claims (1)
1、高光谱图像的空间-光谱信息协同提高分辨率的方法,其特征在于它的步骤如下:
步骤一:输入高光谱图像数据:根据数据格式读入高光谱图像数据,将所述高光谱图像数据同时输入到计算机的步骤A中的空间信息提取模块和步骤B中的光谱信息提取模块中进行信息提取;
步骤二:本步骤包括步骤A和步骤B,
步骤A:空间信息提取:通过空间信息模块进行特征波段选择和空间局部分析及判断;
步骤A一:特征波段选择:特征波段选择是对高光谱图像数据信息量大的特征谱段进行提取或选择;
步骤A二:空间局部分析及判断:利用微分算子边缘检测算法提取步骤A一中所提取或选择的特征谱段的边缘信息,通过边缘信息进行边缘像元判断;是非边缘像元,利用复制技术进行插值处理,得到空间类别信息和相互关系信息,输入到步骤三中;是边缘像元,即为混合像元所在,获取混合像元,将混合像元数据输入到步骤B二中;
步骤B:光谱信息提取:通过光谱信息提取模块进行光谱端元提取和混合像元分解;
步骤B一:光谱端元提取:利用PPI或者N-FINDR算法提取高光谱图像中的光谱端元,记录端元的光谱信息和位置信息;
步骤B二:混合像元分解:用得到的端元的光谱信息和位置信息通过支持向量机方法,从混合像元中获取各组成成分的比例;
步骤三:空-谱协同超分辨:对步骤B二中得到的混合象元解混后的各组成成分的比例分布和步骤A二中得到的非边缘像元所取得的空间类别信息和相互关系信息,按照相关性准则进行分布;
步骤四:得到分辨率提高的高光谱图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007101444495A CN100514085C (zh) | 2007-10-16 | 2007-10-16 | 高光谱图像的空间-光谱信息协同提高分辨率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007101444495A CN100514085C (zh) | 2007-10-16 | 2007-10-16 | 高光谱图像的空间-光谱信息协同提高分辨率的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101140325A CN101140325A (zh) | 2008-03-12 |
CN100514085C true CN100514085C (zh) | 2009-07-15 |
Family
ID=39192351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007101444495A Expired - Fee Related CN100514085C (zh) | 2007-10-16 | 2007-10-16 | 高光谱图像的空间-光谱信息协同提高分辨率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100514085C (zh) |
Families Citing this family (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101477197B (zh) * | 2009-01-24 | 2011-12-21 | 哈尔滨工业大学 | 用于林地复杂场景高光谱遥感数据的仿真方法 |
CN101540049B (zh) * | 2009-04-29 | 2012-02-08 | 中国人民解放军海军装备研究院信息工程技术研究所 | 一种高光谱图像的端元提取方法 |
CN101697008B (zh) * | 2009-10-20 | 2012-07-04 | 北京航空航天大学 | 一种自动估计正则化参数的高光谱解混方法 |
CN101702021B (zh) * | 2009-11-06 | 2011-10-26 | 华中科技大学 | 一种图谱一体化的时变对象光谱信息获取方法与装置 |
CN102117483B (zh) * | 2009-12-31 | 2014-11-05 | 核工业北京地质研究院 | 不同空间分辨率的多光谱遥感图像融合方法 |
CN102074009B (zh) * | 2011-01-06 | 2012-07-11 | 哈尔滨工程大学 | 高光谱图像多端元模式的光谱混合分析方法 |
US8989501B2 (en) * | 2012-08-17 | 2015-03-24 | Ge Aviation Systems Llc | Method of selecting an algorithm for use in processing hyperspectral data |
CN103268502B (zh) * | 2013-06-05 | 2016-04-20 | 中国科学院对地观测与数字地球科学中心 | 高光谱图像波段选择方法及装置 |
CN103471715B (zh) * | 2013-09-02 | 2015-09-09 | 北京航空航天大学 | 一种共光路组合式光场光谱成像方法及装置 |
CN104036298A (zh) * | 2013-09-23 | 2014-09-10 | 苏州工业职业技术学院 | 基于核Fisher自适应学习的高光谱遥感图像端元分类方法 |
CN104574283B (zh) * | 2013-10-23 | 2018-07-27 | 核工业北京地质研究院 | 一种用于叶蜡石信息提取的高光谱影像处理方法 |
CN104732488B (zh) * | 2013-12-20 | 2018-07-06 | 核工业北京地质研究院 | 一种用于阳起石信息提取的高光谱影像处理方法 |
CN104236710B (zh) * | 2014-09-29 | 2015-12-02 | 杭州彩谱科技有限公司 | 一种手持式光源颜色照度测量仪的光谱超分辨方法 |
CN105260993A (zh) * | 2015-10-23 | 2016-01-20 | 哈尔滨工业大学 | 基于边缘约束和自适应形态学滤波的高光谱图像条带缺失修复方法 |
CN106097252B (zh) * | 2016-06-23 | 2019-03-12 | 哈尔滨工业大学 | 基于图Graph模型的高光谱图像超像素分割方法 |
CN106595873B (zh) * | 2017-01-03 | 2019-03-29 | 哈尔滨工业大学 | 基于长波红外大气底层辐射和可见光波段线性混合模型的亚像元温度反演方法 |
CN108593104B (zh) * | 2018-03-30 | 2020-06-26 | 北京化工大学 | 一种小型高信噪比手持式光谱检测系统 |
CN109447951B (zh) * | 2018-09-30 | 2021-06-22 | 中国科学院西安光学精密机械研究所 | 基于吉文斯旋转的高光谱图像端元提取方法 |
CN109708754A (zh) * | 2018-11-11 | 2019-05-03 | 上海卫星工程研究所 | 针对点目标的高光谱探测系统及方法 |
CN109580495B (zh) * | 2018-11-28 | 2021-08-24 | 重庆工商大学 | 一种基于高光谱图像的解混装置及方法 |
CN110111276B (zh) * | 2019-04-29 | 2021-09-10 | 西安理工大学 | 基于空-谱信息深度利用的高光谱遥感影像超分辨率方法 |
CN110119780B (zh) * | 2019-05-10 | 2020-11-27 | 西北工业大学 | 基于生成对抗网络的高光谱图像超分辨重建方法 |
CN111638185B (zh) * | 2020-05-09 | 2022-05-17 | 哈尔滨工业大学 | 基于无人机平台的遥感探测方法 |
CN113280919A (zh) * | 2021-05-08 | 2021-08-20 | 中国科学院长春光学精密机械与物理研究所 | 目标物体的光谱探测方法、光谱探测装置及存储介质 |
CN113446998B (zh) * | 2021-06-29 | 2022-09-30 | 哈尔滨工业大学 | 一种基于高光谱目标探测数据的动态解混方法 |
-
2007
- 2007-10-16 CN CNB2007101444495A patent/CN100514085C/zh not_active Expired - Fee Related
Non-Patent Citations (1)
Title |
---|
基于多特征多分辨率融合的高光谱图像分类. 张钧萍,张晔.红外与毫米波学报,第23卷第5期. 2004 * |
Also Published As
Publication number | Publication date |
---|---|
CN101140325A (zh) | 2008-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100514085C (zh) | 高光谱图像的空间-光谱信息协同提高分辨率的方法 | |
Im et al. | A change detection model based on neighborhood correlation image analysis and decision tree classification | |
Gómez-Chova et al. | Multimodal classification of remote sensing images: A review and future directions | |
Wan et al. | An object-based hierarchical compound classification method for change detection in heterogeneous optical and SAR images | |
CN103971364B (zh) | 基于加权Gabor小波特征和两级聚类的遥感图像变化检测方法 | |
CN104331698A (zh) | 一种遥感图像城区提取方法 | |
Lefebvre et al. | Object-oriented approach and texture analysis for change detection in very high resolution images | |
Kumar et al. | Improving image classification in a complex wetland ecosystem through image fusion techniques | |
CN111639587A (zh) | 基于多尺度谱空卷积神经网络的高光谱图像分类方法 | |
Rangzan et al. | Supervised cross-fusion method: a new triplet approach to fuse thermal, radar, and optical satellite data for land use classification | |
CN104680184A (zh) | 基于深度rpca的极化sar地物分类方法 | |
Hasanlou et al. | A sub-pixel multiple change detection approach for hyperspectral imagery | |
Rayegani et al. | An effective approach to selecting the appropriate pan-sharpening method in digital change detection of natural ecosystems | |
US7283664B2 (en) | Interactive computer aided analysis of images from sensor combinations | |
Mansourmoghaddam et al. | Evaluating the capability of spatial and spectral fusion in land-cover mapping enhancement | |
Boldt et al. | SAR image segmentation using morphological attribute profiles | |
Padrón-Hidalgo et al. | Unsupervised anomaly and change detection with multivariate gaussianization | |
Zhang et al. | Forest roads extraction through a convolution neural network aided method | |
Hosseiny et al. | A convolutional neural network for flood mapping using sentinel-1 and Srtm Dem data: case study in Poldokhtar-Iran | |
Ebert-Uphoff et al. | Evaluation, tuning and interpretation of neural networks for meteorological applications | |
Samadzadegan et al. | Evaluating the sensitivity of image fusion quality metrics to image degradation in satellite imagery | |
Gupta et al. | A comparative investigation of image fusion in the context of classification | |
Karan et al. | Dual-tree complex wavelet transform-based image enhancement for accurate long-term change assessment in coal mining areas | |
Sengar et al. | Soft computing approach for liquefaction identification using Landsat-7 temporal indices data | |
Iervolino et al. | Land Classification using a novel Multispectral and SAR data Fusion in Doha area |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090715 Termination date: 20111016 |