[发明专利]一种空间引力描述下的高光谱图像亚像元定位方法

摘要

本发明提供的是一种空间引力描述下的高光谱图像亚像元定位方法。包括高光谱图像的预处理；SPSAM方法的初始化；完全贯彻空间相关性的空间引力描述。亚像元/像元空间引力方法初始化是利用亚像元/像元空间引力方法获得初步的亚像元定位结果，以作为下一步的输入；完全贯彻空间相关性的空间引力描述方法是指既考虑像元间的空间相关性，又考虑像元内的亚像元间的相关性，同时像元间的空间相关性的描述是通过计算亚像元与邻域像元内的亚像元之间的空间引力实现的。该方法具有完全贯彻空间相关性，亚像元定位精度高的优点。

主权项

1.一种空间引力描述下的高光谱图像亚像元定位方法，其特征是包括以下三个步骤：步骤一，高光谱图像的预处理：输入高光谱图像，通过光谱混合分析技术，得到高光谱图像混合像元中各个类别的混合比例值，根据设定的亚像元定位的放大比例，对混合比例值进行量化，得到各类别的分量图；步骤二，SPSAM方法的初始化：根据各类别的分量图，用SPSAM方法求得混合像元内各个类别对应的亚像元的初始空间分布，即得初始的亚像元定位结果Initial_result；步骤三，完全贯彻空间相关性的空间引力描述：1、选取待分析混合像元P_ab，操作步骤如下：1)根据各类分量图中的混合比例值，产生R种像元P_ab可能的空间分布，其中σ_copyR种分布源自Initial_result中对应像元的直接复制；2)记为P_ab内的亚像元p_ij受到邻域像元中属于类c的亚像元的空间引力，其中：i，j＝1，2，…，S，S为放大比例，c＝1，2，…，C，C为类别总数；为p_ij受到邻域亚像元中属于类c的亚像元的空间引力，对R种分布均计算空间引力：JPab=Σc=1CΣi=1SΣi=1S(k1Ac,pij+k2Bc,pij)]]>其中，k₁和k₂为两个加权系数；3)根据式JPab=Σc=1CΣi=1SΣi=1S(k1Ac,pij+k2Bc,pij)]]>找出对应引力最大的一种空间分布，并选出对应最小的ρ₁R种分布，被引力最大的分布取代；4)随机选择ρ₂R种分布，对每种分布均交换其中任意两个亚像元的类别；5)重复步骤2)～4)，进行多次迭代；在所有迭代中，找出对应最大的一种空间分布，用于更新P_ab内各个类别对应的亚像元的空间分布；2、对每个待分析混合像元，进行步骤1；3、重复步骤1～2，多次更新。