发明名称 一种基于抗噪的混合像元分解模型的构建方法
摘要 本发明公开了一种基于抗噪的混合像元分解模型的构建方法,包括给出抗噪模型:<img file="DDA0000477240370000011.GIF" wi="469" he="98" />,采用常用的欧几里德距离及交替最小二乘算法推导我们提出的抗噪模型的解混结果,并由均方差准则MSE定量评估抗噪模型混合像元分解性能,采用IS距离对抗噪模型进行优化。本发明通过构建基于抗噪的混合象元分解模型来克服加性噪声、乘性噪声及混合噪声对高光谱图像混合像元分解的影响,本模型根据IS距离的统计特性,有效地克服乘性噪声的影响,从而避免产生混合噪声,获得较为精确地信号重构,提高混合像元分解结果的精度,并对抗噪模型进行理论推导,证明了基于抗噪模型的混合像元分解的可行性及优越性。
申请公布号 CN103886639A 申请公布日期 2014.06.25
申请号 CN201410096568.8 申请日期 2014.03.14
申请人 湖州师范学院 发明人 蒋云良;李春芝;陈晓华
分类号 G06T17/00(2006.01)I 主分类号 G06T17/00(2006.01)I
代理机构 北京天奇智新知识产权代理有限公司 11340 代理人 韩洪
主权项 1.一种基于抗噪的混合像元分解模型的构建方法,依次包括以下步骤:a)在图像处理过程中考虑乘性噪声的影响因素,所述乘性噪声在信号独立性方面满足非高斯分布,且乘性噪声随着图像空间的变化而变化,于是给出抗噪模型:<img file="FDA0000477240350000011.GIF" wi="478" he="101" />其中<img file="FDA0000477240350000012.GIF" wi="379" he="92" />是指包括了多种噪声的混合像元矩阵,<img file="FDA0000477240350000013.GIF" wi="341" he="93" />表示从含有多种噪声的混合像元矩阵X中提取的端元矩阵,<img file="FDA0000477240350000014.GIF" wi="366" he="88" />代表相应的端元丰度矩阵,<img file="FDA0000477240350000015.GIF" wi="358" he="100" />表示随机测量误差矩阵,<img file="FDA0000477240350000016.GIF" wi="356" he="110" />表示乘性噪声;b)采用常用的欧几里德距离及交替最小二乘算法推导我们提出的抗噪模型的解混结果,并由均方差准则MSE定量评估抗噪模型混合像元分解性能:b1)由公式<img file="FDA0000477240350000017.GIF" wi="465" he="89" />可以得到X≈WH<sub>ρ</sub>,从公式X≈WH<sub>ρ</sub>中可推导H:<img file="FDA00004772403500000120.GIF" wi="943" he="122" /><img file="FDA0000477240350000019.GIF" wi="85" he="105" />表示乘性噪声矩阵ρ的伪逆矩阵,令U<sub>W</sub>=(W<sup>T</sup>W)<sup>-1</sup>W<sup>T</sup>为端元矩阵W的伪逆矩阵,那么<img file="FDA00004772403500000121.GIF" wi="544" he="125" />计算抗噪模型<img file="FDA00004772403500000111.GIF" wi="460" he="89" />的丰度矩阵的均方差误差为:<img file="FDA00004772403500000112.GIF" wi="1593" he="189" />b2)若固定丰度矩阵H,求解端元矩阵W,可推导抗噪模型的提取端元的精确性,即首先从公式X≈WH<sub>ρ</sub>推导W,令V<sub>H</sub>=(H<sup>T</sup>H)<sup>-1</sup>H<sup>T</sup>为矩阵H的伪逆矩阵,<img file="FDA00004772403500000115.GIF" wi="585" he="123" />那么<img file="FDA00004772403500000116.GIF" wi="660" he="124" />则计算抗噪模型<img file="FDA00004772403500000117.GIF" wi="468" he="90" />的端元矩阵的均方差误差为:<img file="FDA00004772403500000118.GIF" wi="1468" he="193" />c)基于IS距离将抗噪模型:<img file="FDA00004772403500000119.GIF" wi="463" he="99" />进行优化处理,并转化为新的优化模型:<maths num="0001"><![CDATA[<math><mrow><mi>min</mi><mo>:</mo><msub><mi>D</mi><mi>IS</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>|</mo><mo>|</mo><mi>WH</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>d</mi><mrow><mo>(</mo><msub><mrow><mo>[</mo><mi>X</mi><mo>]</mo></mrow><mi>mn</mi></msub><mo>|</mo><mo>|</mo><msub><mrow><mo>[</mo><mi>WH</mi><mo>]</mo></mrow><mi>mn</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><msub><mi>w</mi><mi>ij</mi></msub><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><msub><mi>h</mi><mi>jt</mi></msub><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo><munder><mi>&Sigma;</mi><mi>j</mi></munder><msub><mi>h</mi><mi>jt</mi></msub><mo>=</mo><mn>1</mn><mo>,</mo></mrow></math>]]></maths>1≤i≤M,1≤j≤K,1≤t≤N;其中D<sub>IS</sub>(X||WH)表示IS距离,d(w||v)表示尺度代价函数,<maths num="0003"><![CDATA[<math><mrow><msub><mi>d</mi><mi>IS</mi></msub><mrow><mo>(</mo><mi>w</mi><mo>|</mo><mo>|</mo><mi>v</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>w</mi><mi>v</mi></mfrac><mo>-log</mo><mfrac><mi>w</mi><mi>v</mi></mfrac><mo>-</mo><mn>1</mn><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mover><mi>L</mi><mo>^</mo></mover><mo>,</mo><mover><mi>F</mi><mo>^</mo></mover></mrow></math>]]></maths>分别表示实测的丰度矩阵及实测的端元矩阵,s.t.表示约束条件,w<sub>ij</sub>、h<sub>jt</sub>分别表示端元矩阵及丰度矩阵中的元素,ζ<sub>H</sub>表示计算的丰度矩阵与实测的丰度矩阵之间的差值,ζ<sub>W</sub>表示计算的端元矩阵与实测的端元矩阵之间的差值。
地址 313000 浙江省湖州市吴兴区学士路1号