发明名称 一种基于同质指数的高光谱遥感图像稀疏解混方法
摘要 一种基于同质指数的高光谱遥感图像稀疏解混方法,包括对高光谱遥感图像的每个像元进行同质分析提取同质指数,根据整幅高光谱遥图像中像元同质指数的值域对每个像元的同质指数进行线性拉伸,将其作为基于全变分的变量分离增广拉格朗日稀疏解混算法中空间正则项的权重,并重新定义算法中邻接像元间丰度的差商算子。本方法提高了稀疏解混的精度,使得解混丰度的空间平滑性更接近图像真实的空间平滑性分布,并且更有效的抑制了噪音对解混结果的影响。本方法在基于高光谱遥感图像的高精度地物分类以及地面目标检测和识别方面具有重要的应用价值。
申请公布号 CN103761742A 申请公布日期 2014.04.30
申请号 CN201410034872.X 申请日期 2014.01.24
申请人 武汉大学 发明人 邵振峰;王毓乾
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 严彦
主权项 1.一种基于同质指数的高光谱遥感图像稀疏解混方法,其特征在于,包括以下步骤:步骤a,对高光谱遥感图像的每个像元进行同质分析提取同质指数,提取方式如下,对于高光谱遥感图像的某像元P<sub>i,j</sub>,设选择以像元P<sub>i,j</sub>为中心、以n为半径的不包含像元P<sub>i,j</sub>的方形邻域窗口为Ω(P<sub>i,j</sub>),计算像元P<sub>i,j</sub>与它的邻域窗口内像元P<sub>r,c</sub>之间的光谱相似度S(P<sub>r,c</sub>,P<sub>i,j</sub>),按下式将加权平均的光谱相似度作为该像元的同质指数,<maths num="0001"><![CDATA[<math><mrow><mi>hi</mi><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><msub><mi>P</mi><mrow><mi>r</mi><mo>,</mo><mi>c</mi><mo>&Element;</mo><mi>&Omega;</mi><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow></mrow></msub></munder><mi>&omega;</mi><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>r</mi><mo>,</mo><mi>c</mi></mrow></msub><mo>)</mo></mrow><mo>&times;</mo><mi>S</mi><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>r</mi><mo>,</mo><mi>c</mi></mrow></msub><mo>,</mo><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>其中,ω(P<sub>r,c</sub>)表示(r,c)处的空间权重参数;步骤b,根据整幅高光谱遥图像中像元同质指数的值域对每个像元的同质指数进行线性拉伸,线性拉伸方式如下,对于高光谱遥感图像的某像元P<sub>i,j</sub>,按下式拉伸后的同质指数记为<img file="FDA0000461587130000012.GIF" wi="178" he="89" /><img file="FDA0000461587130000019.GIF" wi="752" he="146" />其中,参数<img file="FDA0000461587130000014.GIF" wi="357" he="93" />M为预设参数;步骤c,定义邻接差商算子如下,<img file="FDA0000461587130000015.GIF" wi="366" he="165" />其中,<img file="FDA0000461587130000016.GIF" wi="1300" he="94" />像元P<sub>i,j+1</sub>是像元P<sub>i,j</sub>右边的邻接像元,像元P<sub>i+1,j</sub>是像元P<sub>i,j</sub>下边的邻接像元,x<sub>i,j</sub>、x<sub>i,j+1</sub>和x<sub>i+1,j</sub>分别是像元P<sub>i,j</sub>、像元P<sub>i,j+1</sub>、像元P<sub>i+1,j</sub>对应的端元丰度;且若丰度分布图像是由k行l列像元组成,<img file="FDA0000461587130000017.GIF" wi="284" he="91" /><img file="FDA0000461587130000018.GIF" wi="295" he="94" />步骤d,将步骤c中定义的邻接差商算子作为空间平滑项引入到稀疏解混的目标优化函数中,然后用基于交替方向乘子法策略迭代求解,得到高光谱遥感图像中每个端元的丰度。
地址 430072 湖北省武汉市武昌区珞珈山武汉大学