基于光谱空间映射与纯化的对象光谱特征提取方法技术

本发明专利技术是一种基于光谱空间映射与纯化的对象光谱特征提取方法。首先根据全色与多光谱图像像元之间的空间关系映射得到全色图像初始光谱特征空间S；其次基于全色图像分割结果O定义多光谱像元对象纯度，得到多光谱像元对象纯度图；最后设定对象纯度阈值和与全色图像同大小的标记图像B，对象纯度高于阈值对应B的位置记为1，低于阈值对应B的位置记为0，结合S、B和O统计对象内标记为1的像元初始映射光谱特征的均值得到对象的光谱特征。本发明专利技术适用于高分辨率遥感图像单一数据源全色与多光谱图像之间或者多种分辨率数据源间面向对象分析过程中对象的光谱特征提取，可广泛应用于面向对象的专题目标识别、土地覆盖分类与变化检测等领域。

【技术实现步骤摘要】
基于光谱空间映射与纯化的对象光谱特征提取方法
本专利技术是一种基于光谱空间映射与纯化的对象光谱特征提取方法，适用于高分辨率遥感图像单一数据源全色与多光谱图像之间或者多种不同分辨率数据源间面向对象分析过程中对象的光谱特征提取。本专利技术可广泛应用于面向对象的专题目标识别、土地覆盖分类与变化检测等领域。
技术介绍
随着对地观测技术的快速发展，一系列高分辨率遥感卫星的陆续升空，使得高分辨率对地观测数据的获取能力大幅提升，人类进入到多源高分辨率对地观测数据获取的新时代。传统的遥感图像信息提取都是基于最原始的像元级光谱特征进行的。面向对象的图像处理与分析技术的出现，使得图像处理的基本单元由像元过渡到对象，使得特征空间由单一像元光谱特征扩展到对象的光谱特征。面向对象的对象光谱特征提取是对遥感图像分割对象的光谱特征进行统计，得到各个对象的光谱特征向量，进而构成面向对象专题目标识别和分类等应用的光谱特征空间。面对多源的高分辨率遥感数据，如何实现高分辨率遥感图像空间特征与低分辨率遥感图像光谱特征的有机融合，准确提取分割对象的光谱特征，是高分辨率遥感图像综合应用的热点与难点。不同遥感卫星传感器的波段数目和对应的波谱覆盖范围也不一样。对于中低分辨率遥感数据(例如Landsat-5有7个波段、Landsat-8有11个波段等)具有中等数量的光谱波段数；对于高分辨率光学遥感图像一般包括一个高分辨率全色波段和四个相对低分辨率多光谱波段，光谱测量仅在蓝、绿、红、近红外等范围(例如GF-1、GF-2、SPOT-5等)。高分辨率光学遥感图像提供了更多的空间结构和纹理等细节信息，而较低分辨率的遥感图像则较好地体现了地物的光谱特征，如何实现高分辨率遥感图像空间信息与相对较低分辨率遥感图像光谱信息的有机融合是多源遥感图像面向对象分析与信息提取的关键。常用的处理方法是首先将相对较低分辨率的多光谱图像与较高分辨率的全色图像进行图像融合，然后再基于融合后的高分辨率多光谱图像进行图像分割对象的光谱特征提取。然而，在图像融合的过程中往往会发生光谱信息的失真和丢失，没有最大程度的保留图像的原始光谱特征，图像融合结果的好坏直接决定了特征提取的效果和最终特征应用的精度。同时，基于图像融合的特征提取方法适用于两个数据源图像之间的融合与对象光谱特征提取，对三个及以上不同分辨率数据源的对象光谱特征提取存在一定的局限性。对象的光谱特征是面向对象遥感图像处理与分析的特征基础，研究高效的多源遥感图像分割对象光谱特征提取方法，对于提高多源遥感图像的光谱信息利用效率，促进多源遥感图像综合应用具有十分重要的应用价值。因此，本专利技术从高分辨率与低分辨率遥感图像像元间的空间对应关系出发，提出了基于光谱空间映射与纯化的对象光谱特征提取方法，将高分辨率遥感全色图像与低分辨率遥感多光谱图像像元之间的空间对应关系，扩展到高分辨率全色图像分割结果中对象与低分辨率多光谱图像中像元之间的空间对应关系，基于光谱空间映射和纯化策略，实现高分辨率全色图像空间信息与低分辨率多光谱图像光谱信息的有机融合，减少了图像融合的过程中的光谱丢失和失真，最大程度的保留了地物对象的原始光谱特征，为后续面向对象分析与应用提供更优的光谱特征空间。
技术实现思路
本专利技术是一种实用的面向对象的光谱特征提取方法，将高分辨率全色遥感图像与低分辨率多光谱遥感图像像元间的空间对应关系扩展到高分辨率全色图像分割结果中对象与低分辨率多光谱遥感图像中像元之间的空间对应关系，采用光谱空间映射与纯化的策略，实现多源高分辨率遥感图像空间特征与低分辨率遥感图像光谱特征的有机融合，获得面向对象分析的光谱特征空间。方法具体的步骤为：第一步:生成高分辨率遥感全色图像的初始光谱特征空间S根据高分辨率遥感全色图像与低分辨率遥感多光谱图像之间的分辨率比，建立图像像元之间的空间映射关系。以分辨率比为1:4的遥感数据为例，多光谱图像上的一个像元对应全色图像上一个4×4的图像块，是一种一对多的空间映射关系，因此，给定全色图像上的一个像元坐标就能求得对应的多光谱图像上的像元坐标，这样全色与多光谱图像对的像元之间就建立了空间映射关系，对应的坐标位置用下述公式计算：xi'＝Int((xi-1)/4)+1yi'＝Int((yi-1)/4)+1上述两个公式中，Int()函数表示向下取整，(xi’,yi’)表示多光谱图像中像元的行列位置，(xi,yi)表示全色图像中像元的行列位置。同理，给定一个多光谱像元的坐标位置，能够求得对应的全色图像中像元组的位置。根据高分辨率全色遥感图像与低分辨率多光谱遥感图像像元之间的空间对应关系，可以将低分辨率多光谱图像像元的光谱信息直接映射赋值到全色图像的每个像元上，形成全色图像的初始光谱特征空间S。第二步:基于高分辨率全色图像的分割结果O，定义低分辨率图像多光谱像元的对象纯度低分辨率多光谱图像上的1个像元对应高分辨率全色图像上的16个像元(以全色图像与多光谱图像的空间分辨率比1:4为例)，全色图像上的16个像元可能会隶属于不同的分割对象。如果一个多光谱像素对应的全色图像子块完整的落在同一个分割对象内，那么相应的位置上全色图像像元应该被赋予相同的光谱向量数值，但是，如果一个多光谱像元落在两个甚至多个对象上，那么该多光谱像元认为可能是由多种地物混合而成，直接映射赋值方式会产生光谱的不确定性，因此，为了消除这种地物混合多光谱像元的影响，定义多光谱图像中像元的对象纯度来度量多光谱像元是属于一个对象还是多个对象，计算公式如下：公式中P表示多光谱图像中像元的对象纯度；NMax表示一个低分辨率多光谱像元对应的高分辨率全色图像块的像素个数，由高分辨率全色与低分辨率多光谱图像分辨率的比值确定，以分辨率比1：4为例，则NMax为16；NMaxObject表示多光谱像元对应的全色图像块中，面积最大对象对应的像素数。如果多光谱像元的对象纯度P为1，说明多光谱像元对应全色图像块中所有的像元均在同一个区域内且具有相同的光谱信息；如果多光谱像元的对象纯度P小于1，则说明该多光谱像元会落在两个或者两个以上的分割对象中，可能是多种地物混合的光谱信息。通过对所有的多光谱像元计算，得到多光谱图像的对象纯度图Mp。第三步:基于全色图像的初始光谱特征空间S和多光谱像元的对象纯度图Mp，计算对象的光谱特征在通过光谱空间映射机制得到全色图像的像元光谱特征的过程中，人为设定一个多光谱像元对象纯度阈值Tp和一个与全色图像相同大小的标记图像B，对于Mp中对象纯度大于等于阈值的多光谱像元对应的标记图像像元位置标记为1，对于Mp中小于对象纯度阈值的多光谱像元对应的标记图像像元位置标记为0。最后结合初始光谱特征空间S、标记图像B和全色图像分割结果O，统计对象内部标记值为1的像元初始映射光谱特征的均值来表示该对象的光谱特征。通过这种光谱空间映射与纯化的方式，可以剔除处于两个及以上对象边界处混合多光谱像元的影响，从而更加准确的提取对象的光谱信息。依次对分割结果O上的每个对象，按照第一至三步的方法进行对象的光谱特征提取就得到整个分割图像面向对象的光谱特征空间。对于多种不同分辨率数据源的情况，只需将最高分辨率全色图像依次与多源的多光谱图像分别组成图像对，采用上述方法步骤分别进行对象特征提取，最终组成了对象本文档来自技高网...

【技术保护点】
基于光谱空间映射与纯化的对象光谱特征提取方法，其步骤为：第一步：生成高分辨率遥感全色图像的初始光谱特征空间S根据高分辨率遥感全色图像与低分辨率遥感多光谱图像之间的分辨率比，建立图像像元之间的空间映射关系，以分辨率比为1:4的遥感数据为例，多光谱图像上的一个像元对应全色图像上一个4×4的图像块，是一种一对多的空间映射关系，给定全色图像上的一个像元坐标就能求得对应的多光谱图像上的像元坐标，对应的坐标位置计算公式如下：x'i＝Int((xi‑1)/4)+1y'i＝Int((yi‑1)/4)+1上述两个公式中，Int()函数表示向下取整，(xi’,yi’)表示多光谱像元的行列位置，(xi,yi)表示全色图像中像元的行列位置，根据高分辨率全色遥感图像与低分辨率多光谱遥感图像像元之间的空间对应关系，可以将低分辨率多光谱图像像元的光谱信息直接映射赋值到全色图像的每个像元上，形成全色图像的初始光谱特征空间S；第二步：基于高分辨率全色图像分割结果O，定义低分辨率多光谱图像像元的对象纯度以分辨率比为1:4的遥感数据为例，低分辨率多光谱图像上的1个像元对应高分辨率全色图像上的16个像元，全色图像上的16个像元可能会隶属于不同的分割对象，如果一个多光谱像素对应的全色图像子块完整的落在同一个分割对象内，那么相应的位置上全色图像像元应该被赋予相同的光谱向量数值，但是，如果一个多光谱像元落在两个甚至多个对象上，那么该多光谱像元认为可能是由多种地物混合而成，直接映射赋值方式会产生光谱的不确定性，因此，为了消除这种地物混合多光谱像元的影响，定义多光谱图像中像元的对象纯度来度量多光谱像元是属于一个对象还是多个对象，计算公式如下：...

【技术特征摘要】
1.基于光谱空间映射与纯化的对象光谱特征提取方法，其步骤为：第一步：生成高分辨率遥感全色图像的初始光谱特征空间S根据高分辨率遥感全色图像与低分辨率遥感多光谱图像之间的分辨率比，建立图像像元之间的空间映射关系，以分辨率比为1:4的遥感数据为例，多光谱图像上的一个像元对应全色图像上一个4×4的图像块，是一种一对多的空间映射关系，给定全色图像上的一个像元坐标就能求得对应的多光谱图像上的像元坐标，对应的坐标位置计算公式如下：x'i＝Int((xi-1)/4)+1y'i＝Int((yi-1)/4)+1上述两个公式中，Int()函数表示向下取整，(xi’,yi’)表示多光谱像元的行列位置，(xi,yi)表示全色图像中像元的行列位置，根据高分辨率全色遥感图像与低分辨率多光谱遥感图像像元之间的空间对应关系，可以将低分辨率多光谱图像像元的光谱信息直接映射赋值到全色图像的每个像元上，形成全色图像的初始光谱特征空间S；第二步：基于高分辨率全色图像分割结果O，定义低分辨率多光谱图像像元的对象纯度以分辨率比为1:4的遥感数据为例，低分辨率多光谱图像上的1个像元对应高分辨率全色图像上的16个像元，全色图像上的16个像元可能会隶属于不同的分割对象，如果一个多光谱像素对应的全色图像子块完整的落在同一个分割对象内，那么相应的位置上全色图像像元应该被赋予相同的光谱向量数值，但是，如果一个多光谱像元落在两个甚至多个对象上，那么该多光谱像元认为可能是由多种地物混合而成，直接映射赋值方式会产生光谱的不确定性，...

【专利技术属性】
技术研发人员：何国金，王桂周，袁继颖，张兆明，王猛猛，彭燕，程博，白雅卿，张晓美，刘慧婵，
申请(专利权)人：三亚中科遥感研究所，
类型：发明
国别省市：海南,46