专利名称：基于修正扩展形态学算子的高光谱遥感图像端元提取方法
技术领域：
本发明提出了一种基于修正扩展形态学算子的高光谱遥感图像端元提取方法，结合高光谱遥感图像的空间和光谱信息用于端元提取，属于高光谱遥感数字图像处理和模式识别领域。本发明将形态学概念用于高光谱遥感图像端元自动提取，通过修正排序规则和替换准则，解决了传统扩展形态学算子在端元提取时造成的错误替换现象。适用于高光谱遥感图像的端元提取、混合像元分析、地物分布制图等工作。
背景技术：
高光谱遥感成像技术的出现和发展是对地观测方面的一次革命，高光谱遥感卫星可以捕获某ー场景下地物所反射、散射、吸收和发射的电磁波谱，这些波谱信息反映了地物空间特征、辐射特征和光谱特征，信息量及其丰富。针对高光谱遥感图像的特性，在将其用于具体遥感应用前衍生了ー些特有的图像预处理技术，包括端元提取、混合像元分解等。一幅高光谱遥感图像中由单一地物形成的像元称为端元，而通常情况是，大部分像元所表示的瞬时现场光谱值由存在于该现场的多种地物生成，这样的像元称为混合像元。高光谱成像仪有限的空间分辨率和复杂的地貌情况是混合像元大量存在于高光谱遥感图像的主要原因。在目前高光谱遥感图像相关的研究領域，端元提取技术占有重要的地位，已经发展出很多经典有效的端元提取方法，各种方法因侧重点不同而具有各自的优势和局限性。其中包括基于凸面几何学概念的N-FINDR、最小体积转换(Minimum Volume Transforms,MVT)方法；以平方误差为理论基础结合部分解混的迭代误差分析(Iterative ErrorAnalysis, IEA)方法；独立成份分析(Independent Component Analysis, I CA)方法；向特征空间中向量投影，以极值点数目表示像元纯度的纯像元指数(Pixel Purity Index7PPI)方法等。上述端元提取方法均是从高光谱数据的光谱特征信息和像元光谱间的相关性入手，并没有使用像素之间的空间相关性。将形态学算子引入高光谱遥感图像端元提取中，结合光谱信息和空间信息的自动形态学端元提取(Automated Morphological EndmemberExtraction, AMEE)方法,是将形态学概念扩展应用的ー个实例。形态学算子是ー种非线性图像处理方法，其技术基于空间代数学、拓扑学、随机几何学等基本概念的研究，最初应用于ニ值图像处理，而后扩展到灰度图像。形态学概念应用于高光谱遥感图像涉及如何处理多维光谱向量的问题，在ー些分类应用中，采用经典的形态学算子分别处理各个波段的数据(也可以是各个主成分)然后将处理结果融合成向量格式，这种方式会造成某些像元光谱向量的变化；将光谱向量看作ー个整体是形态学概念在高光谱遥感图像应用中的另ー个研究研究方向，派生出在结构元素(窗ロ)中采用何种排序规则和替换准则等问题，由此产生了扩展形态学算子(Extended-MorphologicalOperator, EM0),此种应用方式多用于端元提取和混合像元分析。AMEE方法是扩展形态学算子应用于高光谱遥感图像端元提取的ー个范本，其初衷是经过形态学滤波进而分离纯像元和混合像元，然而通过研究发现，AMEE在应用形态学概念时，所设定的排序规则和替换准则存在一定局限性，使得其对纯像元和混合像元分离的能力有限，而且容易造成图像中两类地物边界处像元的交叉替换现象，进而影响最终端元提取的精度以及后续混合像元解混的效果。
发明内容
要解决的技术问题为了避免现有技术的不足之处，针对扩展形态学算子设定的排序规则和替换准则存在的局限性，本发明从其公式出发，分析了局限性产生的原因，并通过引入基准向量，提出了修正方案，给出基于修正扩展形态学膨胀和腐蚀运算的开运算、闭运算、开-闭运算、闭-开运算公式。将上述开-闭运算和闭-开运算作用于原始高光谱遥感图像，计算用于端元提取的判定向量，作为本发明提出的AMEMEE方法的主要依据。本发明要解决的技术问 题包括对扩展形态学算子排序规则和替换准则的修正、基准向量的计算以及AMEMEE方法的步骤。技术方案一种基于修正扩展形态学算子的高光谱遥感图像端元提取方法，其特征在于步骤如下步骤I :对于高光谱遥感图像I = f(x,y,N), (x,y) e Z2为空间坐标，N为波段数，f(x，y)为N维光谱向量。计算均值向量作为基准向量
Γ ^Σル，ァ)/
_1] '了/
VRef _/M其中M为像元数目；步骤2 :根据步骤2中计算的Vltef,对原始图像分别做η次本方法中提出的基于修正扩展形态学算子的开-闭运算(/』)广レ，ァΗ(/°句腦パLwo (ぶ，ァ)η介于8 12之间，姆做一次运算,需要重新返回步骤I计算VKef。其中，(foB)MEMQ和び· B)memo为基于修正扩展形态学算子的开运算和闭运算(/oU,ァ) = ((/ UU ァ){f*B)MEMO{x,y) = {{f@B)MEMO B)MEMo{x,y)由修正扩展形态学算子中的膨胀和腐蚀运算定义
{f B)MEM0 {x,y) = Dist\{f B){x,y),VRef] vDist[f{x,y),VRef]{f B)MEM0{x,y) = Dist[{f B){x,y),yKet]ADist]lf{x,y),yKet]其中B代表空间元素SE，SE从3X3方形结构开始，毎次运算时空间尺寸按I阶逐次递增，(/ 句レび)和(/ 句レび)是传统扩展形态学算子定义下的膨胀和腐蚀运算，DistG为高光谱遥感图像中两个光谱向量的距离测度，用以表征光谱间的相似性，使用光谱角
DistSAD (aN,bN ) = = cos—1!卜"
IIaN Il IIdN Il其中aN、bN为N维光谱向量；步骤3 :步骤2完成后，得出每个像元的判定向量DV-y用于端元提取DVopen^close (X，y) = { α 0，L a n_J=Dist[{f,B)iMTe [x,パ人m°=se {x,yr\步骤4 :根据步骤3得出的判定向量DVt5pra^lt5'将含有最大数值元素DV°pm_d°se所 对应的像元作为ー个端元输出，记为ep步骤5 :在原图像中消除步骤4中已提取端元影响，将原图像投影到ち的正交子空间/(x,j) = Pぐ/(x,j)，其中！^=为ち的正交子空间；将受到压制像元所对应
的DV-mi·置为零向量。步骤6 :重复步骤5，直到所有的都变为零向量，得到最終端元集合。有益效果本发明提出的一种适用于高光谱遥感图像端元提取的修正扩展形态学算子，和以此为依据的AMEMEE方法，是ー种结合光谱信息和空间信息的端元提取方法。修正扩展形态学算子通过引入基准向量，修正了排序规则和替换准则，避免了在用扩展形态学算子对高光谱遥感图像进行滤波时产生的错误替换现象，増大了纯像元和混合像元的可分性。AMEMEE方法可以自动提取端元，不须先验知识，运算量适中(主要由形态学滤波次数决定)，端元提取效果优于仅使用光谱信息的方法和基于扩展形态学算子的AMEE方法。本发明在AMEE方法所使用的扩展形态学算子的基础上，引入基准向量的概念，设定了更为严谨的排序规则和替换准则，提出了适用于高光谱遥感图像端元提取的修正扩展形态学算子(Modified Extended-Morphological Operator, MEMO),并给出了基于修正扩展形态学算子的高光谱端元提取方法——自动修正扩展形态学端元提取(AutomatedModified Extended-Morphological Endmember Extraction, AMEMEE)方法。本发明提出的AMEMEE方法由于采用了修正后的扩展形态学算子，避免了上面提到的交叉替换现象，可以更好的分离纯像元和混合像元，从而增强了端元提取的效果。应用于实际高光谱遥感图像的实验数据从光谱曲线相似性、相应地物分布制图等方面证明了本发明提出方法的有效性。
图I :本发明方法AMEMEE的基本流程图。图2 :本发明方法AMEMEE在AVIRIS于美国内华达州Cuprite地区获取的高光谱遥感图像上端元提取的结果与參考光谱的比较图。图2(a)给出了 Cuprite地区5种常见矿物USGS光谱库中的实测光谱曲线，以及AMEMEE提取的端元与实测光谱曲线以及ENVI公布的该地物平均端元光谱曲线的比较图(b) (f)，为便于观察，光谱曲线采用堆叠方式排列，在(b) (f)中，本文提取端元光谱曲线用实线表示，USGS光谱库中參考光谱用点线表示，ENVI给出的平均光谱用虚线表示。从图中可以看出，本文提取的端元光谱曲线与參考光谱有较好的吻合，正确反映了不同矿物在不同波段的吸收特征。
图3 :本发明方法AMEMEE提取的端元集线性解混得到的Cuprite地区三种常见矿物的分布。通过比较可以看出，图3显示的地物分布与USGS公布的Cuprite地区矿物分布能较好的吻合，从侧面证明了本发明端元提取方法的有效性。
具体实施例方式现结合实施例、附图对本发明作进ー步描述用于实施的硬件环境是Pentium-43. OG计算机、IG内存、128M显卡；运行的软件环境是=Window XP操作系统，用IDL7. O程序设计语言结合ENVI实现了本发明提出的方法。I.对于高光谱遥感图像中某一 N维光谱向量f(x，y)，N为波段数，扩展形态学算子给出在结构元素SE (Structure Element)中的排序规则
权利要求
1.一种基于修正扩展形态学算子的高光谱遥感图像端元提取方法，其特征在于步骤如下 步骤I :对于高光谱遥感图像I = f (x,y,N), (x,y) G Z2为空间坐标,N为波段数，f (x,y)为N维光谱向量；计算高光谱遥感图像的均值向量作为基准向量
全文摘要
本发明提出了一种基于修正扩展形态学算子的高光谱遥感图像端元提取方法，其特征是，通过修正扩展形态学算子定义下的开-闭运算和闭-开运算对原始图像进行滤波，以达到结合光谱信息和空间信息进行端元提取的目的。本方法主要着重于解决存在于扩展形态学算子中的局限性问题，通过引入基准向量对其进行修正，提出的修正扩展形态学算子去除了交叉替换现象，保证了正确的替换方向，增强了纯像元和混合像元的分离效果。实验结果从光谱曲线相似性和相应地物分布图等方面证明了本发明提出的方法可以优化端元提取的效果，且计算量适中，是后续高光谱遥感图像混合像元分析以及分类等应用的基础。
文档编号G06K9/46GK102663402SQ20121011886
公开日2012年9月12日 申请日期2012年4月21日 优先权日2012年4月21日
发明者梁楠, 王瀛, 郭雷 申请人:西北工业大学