一种高光谱图像非线性光谱混合模型及地质填图应用研究制造技术

本发明专利技术属于高光谱遥感影像处理技术领域，具体为一种高光谱遥感影像光谱解混方法。针对线性光谱混合模型解混精度不高的问题，提出一种基于二次散射的高光谱遥感影像非线性光谱混合模型，并定义了一种模型中二次散射系数的确定方法。研究了基于该模型解混结果的遥感地质填图，提出了3种填图规则，通过与线性模型填图结果对比，利用二次散射模型进行解混能得到更精确的填图结果。

【技术实现步骤摘要】
一种高光谱图像非线性光谱混合模型及地质填图应用研究
本专利技术属于高光谱遥感影像处理
，具体为一种高光谱遥感影像光谱解混方法。
技术介绍
20世纪80年代兴起了新型的对地观测技术-高光谱遥感技术。它把人们研究地表信息的能力由陆地推向太空，拓宽了人们的研究视野，极大地提高了人们宏观、准确、及时、综合地对地观测与监测能力。成像光谱技术把遥感波段从几个、几十个推向数百个、上千个，其覆盖的光谱范围也从可见近红外、短波红外，逐步向紫外和中红外拓展。高光谱遥感数据每个像元可以提供几乎连续的地物光谱曲线，使人们利用高光谱反演陆地细节成为可能。其应用领域已涵盖地球科学的各个方面，在地质找矿和制图、大气和环境监测、农业和森林调查、海洋生物和物理研究等领域发挥着越来越重要的作用。高光谱遥感技术已经成为当前遥感领域的前沿技术。然而，在遥感成像系统、高光谱遥感成像光谱仪的设计中，光谱分辨率和空间分辨率是不可兼得的(参考对比文件1)。高光谱遥感光谱分辨率高，一般小于10nm，这使得其空间分辨率较低。较低的空间分辨率导致混合像元大量存在，纯净像元相对较少，给数据处理如精准匹配、光谱解混和目标检测与识别等技术带来了巨大的困难。针对高光谱图像中混合像元大量存在的问题，混合像元分解技术应运而生。其目的是求解各种不同物质(即端元)在混合像元中所占的比例(丰度)。高光谱图像混合像元分解一般称为光谱解混(SpectralUnmixing)，因为高光谱图像混合像元的光谱为各个纯像元光谱以某种方式的混合。在进行混合像元分解之前需要进行端元提取，以确定各种纯像元的光谱信息，再利用各端元的光谱信息确定其在混合像元内的比例。从理论上讲，光谱解混首先需要确定光谱混合模型，当前研究最广泛最深入的光谱混合模型是线性光谱混合模型(参考对比文件2)，然而它不能正确反映光谱的混合方式，使得解混精度不高，非线性光谱混合模型能更好的体现光谱混合方式，能得到更精确的解混结果(参考对比文件3)。对比文件1：马艳华，高空间分辨率和高光谱分辨率遥感图像的融合[J]，红外，2003，10：11-16对比文件2：GautamGhosh，SureshKumarandS.K.Saha，“HyperspectralSatelliteDatainMappingSalt-AffectedSoilsUsingLinearSpectralUnmixingAnalysis”，JOURNALOFTHEINDIANSOCIETYOFREMOTESENSING，Volume40，Number1，129-136，2012对比文件3：A.Halimi，Y.Altmann，N.Dobigeon，andJ.-Y.Tourneret，“Nonlinearunmixingofhyperspectralimagesusingageneralizedbilinearmodel，”IEEETrans.Geosci.andRemoteSensing，49(11)：4153-4162，2011
技术实现思路
根据地物反射机理，专利技术一种基于二次散射的高光谱遥感影像非线性光谱混合模型，该模型考虑了光子在地物间发生二次散射作用，模型如下：y＝c(Ma)+d(Ma)2+nc+d＝1，c≥0，d≥0其中，y＝[y1，…yL]T为L波段的高光谱图像，每个波段用矩阵的一行表示；M为L×R的混合矩阵，每一列为一端元(即纯地物)的光谱；R为端元的个数；a＝[a1，…aR]T为丰度向量；n为误差项，一般认为0均值，方差为σ2的高斯噪声；d叫做二次散射系数。考虑到不同混合像元的二次散射系数应该是不同的，提出二次散射模型的改进模型，模型如下：yij＝cij(Ma)+dij(Ma)2+n，0≤cij∈C≤1，0≤dij∈D≤1，cij+dij＝1其中，yij表示遥感图像第i行第j列混合像元的混合光谱向量，dij为该混合像元的二次散射系数。提出一种基于丰度方差的确定二次散射系数的方法，定义的二次散射系数如下：其中，dmax是最大的二次散射系数(设为0.5)，var(i，j)为地物丰度方差，varmax为各混合像元中地物丰度方差的最大值。由此定义式可知，每个混合像元的二次散射系数为0到dmax之间。根据该模型进行光谱解混，得到不同地物的丰度信息。提出了一种基于光谱解混结果的地质填图方法(见具体实施方式)。附图说明图1为两种地物二次散射情况图2为地物分布复杂度：(a)两种地物均等分布，(b)不均等分布图3为地质填图规则1流程图图4为地质填图规则2流程图图5为地质填图规则3流程图图6为cuprite地区第100波段遥感影像图7为USGS1995年绘制的cuprite地区地质填图图8为二次散射模型解混结果图9为利用规则1分别基于线性解混和非线性解混得到的地质填图。图10为利用规则2分别基于线性解混和非线性解混得到的地质填图。图11为利用规则3分别基于线性解混和非线性解混得到的地质填图。具体实施方式●模型推导：假设某一混合像元内只存在A、B两种地物(图1)，两种地物的空间位置分布及比例都未知，设两种地物的光谱向量分别为m1和m2，其丰度向量分别为a1和a2，设M＝[m1，m2]，a＝[a1，a2]T。混合方式分线性和二次散射两大类：线性情况：(a)光子入射到地物A后直接反射到传感器中；(b)光子在入射到地物B后直接反射到传感器中。二次散射情况：光子在到达传感器前的二次撞击包含4种情况，(c)入射光子射到地物A，散射到地物A后反射到传感器；(d)入射光子射到地物A，散射到地物B后反射到传感器；(e)入射光子射到地物B，散射到地物A后反射到传感器；(f)入射光子射到地物B，散射到地物B后反射到传感器(图2)。光子在入射到某种地物上后，将以概率c直接反射到传感器。以概率d(d＝1-c)散射到其他地物，散射的这部分光子又分别以概率a1和概率a2碰撞到地物A和地物B，之后反射到传感器中。该混合像元得到的光谱为：y＝ca1m1+da1m1a1m1+da1m1a2m2+ca2m2+da2m2a1m1+da2m2a2m2+n上式中的前6项分别对应图1中的(a)、(c)、(d)、(b)、(e)、(f)六种情况，第7项为误差项。上式可进一步化简为：y＝ca1m1+ca2m2+2da1m1a2m2+d(a1m1)2+d(a2m2)2+n＝c(a1m1+a2m2)+d(a1m1+a2m2)2+n＝c(Ma)+d(Ma)2+n●二次散射系数确定：二次散射模型是半确定性模型，需要首先确定二次散射系数d(c＝1-d)。二次散射系数直观反映的是地物分布复杂程度，地物分布越简单二次散射系数越低，极限情况d＝0(线性混合)，相反地物分布越复杂，非线性程度越强，则二次散射系数越高。对于遥感影像，完全靠人力实地调查确定地物复杂程度是不现实的。而只根据一幅遥感图像，缺乏必要的先验知识，很难确定混合像元地物复杂程度。二次散射情况主要发生在不同地物相交界的地方，交界地带本文档来自技高网...

【技术保护点】
在线性光谱混合模型的基础上，提出一种基于二次散射的非线性光谱混合模型，该模型考虑了光子在地物间发生二次散射作用，模型如下：y＝c(Ma)+d(Ma)2+nc+d＝1，c≥0，d≥00≤ai≤1,Σi=1Rai=1其中，y＝[y1，…yL]T为L波段的高光谱图像，每个波段用矩阵的一行表示；M为L×R的混合矩阵，每一列为一端元(即纯地物)的光谱；R为端元的个数；a＝[a1，…aR]T为丰度向量；n为误差项，一般认为0均值，方差为σ2的高斯噪声；d叫做二次散射系数。

【技术特征摘要】
1.一种基于光谱解混的地质填图方法，其特征在于，所述光谱解混是采用基于二次散射的非线性光谱混合模型；填图规则1：将混合像元内比例最大的地物作为该像元的地物类型进行填图；填图规则2：首先对混合像元内所有地物的丰度从大到小排序，判断最大比例是否大于80％，若是，则将该像元标为对应该比例的地物类型，若否，则将最大比例的两种地物作为该像元类型进行填图；填图规则3：首先对混合像元内所有地物的丰度从大到小排序，首先判断最大比例是否大于80％，若是，则将该像元标为对应该比例的地物类型，若否，则将第一和第二比例相加，若大于80％，则将该像元标为对应该两种比例的地物混合类型，若否，则继续加入第三比例，将该像元标为对应该三种比例的地物混合类型；所述非线性光谱混合模型考虑了光子在地物间发生二次散射作用，所述非线性光谱混合模型如下：y＝c(Ma)+d(Ma)2+nc+d＝1，c≥0，d≥0其中，y＝[y1，...yL]T为L波段的高光谱图像，每...

【专利技术属性】
技术研发人员：余先川，李建广，徐金东，
申请(专利权)人：北京师范大学，
类型：发明
国别省市：