发明名称 一种基于鲁棒低秩张量的高光谱图像去噪方法
摘要 本发明提出了一种基于鲁棒低秩张量的高光谱图像去噪方法,包括建立高光谱图像噪声的数学模型,构造高光谱图像鲁棒低秩张量(RLRTR)去噪模型,求解RLRTR去噪优化模型。本发明充分利用高光谱图像(HSI)的先验知识,高光谱图像被不同的噪声污染,如高斯噪声、脉冲噪声、死像素和条带噪声等。利用干净的高光谱图像数据具有潜在的低秩张量特性以及异常和非高斯噪声具有稀疏性的特性,同时分别采用核范数和<i>l</i><sub>2,1</sub>范数来表征低秩和稀疏特性;本发明的技术方案充分利用高光谱图像的先验信息和内在结构特征,可以同时去除高斯噪声、异常和非高斯噪声。
申请公布号 CN105069758A 申请公布日期 2015.11.18
申请号 CN201510521057.0 申请日期 2015.08.21
申请人 武汉大学 发明人 马佳义;马泳;黄珺;梅晓光;樊凡
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 张火春
主权项 一种基于鲁棒低秩张量的高光谱图像去噪方法,其特征在于,包括以下步骤:步骤1,建立高光谱图像噪声的数学模型所述的高光谱图像的噪声模型如下:Y=X+S+N其中,Y表示带噪的高光谱图像,X表示干净的高光谱图像,S表示异常和非高斯噪声N表示高斯噪声;<img file="FDA00007865638000000110.GIF" wi="432" he="71" />其中I<sub>1</sub>和I<sub>2</sub>分别表示高光谱图像在空间维度上的宽度和高度,I<sub>3</sub>表示高光谱图像在光谱维度上的波段数;R为实数;步骤2,构造高光谱图像鲁棒低秩张量RLRTR去噪优化模型,去噪优化模型的数学表达式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><mi>X</mi><mo>,</mo><mi>S</mi></mrow></munder><mo>{</mo><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><mo>|</mo><mo>|</mo><msub><mrow><mo>(</mo><msub><mi>X</mi><mi>n</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msub><mo>|</mo><msub><mo>|</mo><mo>*</mo></msub><mo>+</mo><mi>&lambda;</mi><mo>|</mo><mo>|</mo><msub><mrow><mo>(</mo><msub><mi>S</mi><mi>n</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msub><mo>|</mo><msub><mo>|</mo><mrow><mn>2</mn><mo>,</mo><mn>1</mn></mrow></msub><mo>}</mo><mo>,</mo><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mo>|</mo><mo>|</mo><msub><mi>Y</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msub><mo>-</mo><msub><mrow><mo>(</mo><msub><mi>X</mi><mi>n</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msub><mo>-</mo><msub><mrow><mo>(</mo><msub><mi>S</mi><mi>n</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msub><mo>|</mo><msub><mo>|</mo><mi>F</mi></msub><mo>&le;</mo><mi>&delta;</mi><mo>,</mo><mi>n</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>3</mn><mo>,</mo></mrow>]]></math><img file="FDA0000786563800000011.GIF" wi="2055" he="173" /></maths>其中,<img file="FDA0000786563800000012.GIF" wi="97" he="92" />是最小化算子,<img file="FDA0000786563800000013.GIF" wi="73" he="149" />表示n从1到3求和,s.t.表示使得条件满足;其中<img file="FDA0000786563800000014.GIF" wi="266" he="81" />和<img file="FDA0000786563800000015.GIF" wi="276" he="81" />n=1,2和3;X<sub>n</sub>,S<sub>n</sub>分别是X和S引入的辅助变量;X中的每个元素表示成<img file="FDA0000786563800000016.GIF" wi="124" he="68" />其中索引i<sub>1</sub>=1,…,I<sub>1</sub>;i<sub>2</sub>=1,…,I<sub>2</sub>;i<sub>3</sub>=1,…,I<sub>3</sub>,张量的第n维称为n‑模,Z<sub>(n)</sub>表示各种张量Z的n‑模展开矩阵,Z代表任何张量;其中<img file="FDA0000786563800000017.GIF" wi="293" he="87" />n=1,2,3,M<sub>1</sub>=I<sub>2</sub>×I<sub>3</sub>,M<sub>2</sub>=I<sub>1</sub>×I<sub>3</sub>,M<sub>3</sub>=I<sub>1</sub>×I<sub>2,</sub><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>Z</mi><mo>&Element;</mo><msup><mi>R</mi><mrow><msub><mi>I</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>I</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>I</mi><mn>3</mn></msub></mrow></msup><mo>,</mo></mrow>]]></math><img file="FDA0000786563800000018.GIF" wi="262" he="69" /></maths>Y<sub>(n)</sub>表示张量Y的n‑模展开矩阵,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mo>|</mo><mo>|</mo><msub><mrow><mo>(</mo><msub><mi>X</mi><mi>n</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msub><mo>|</mo><msub><mo>|</mo><mo>*</mo></msub><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>I</mi></msubsup><msub><mi>&sigma;</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mrow><mo>(</mo><msub><mi>X</mi><mi>n</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000786563800000019.GIF" wi="595" he="99" /></maths>表示矩阵(X<sub>n</sub>)<sub>(n)</sub>的奇异值的和,||(S<sub>n</sub>)<sub>(n)</sub>||<sub>2,1</sub>表示矩阵(S<sub>n</sub>)<sub>(n)</sub>的每行的l<sub>2</sub>范数的和,||.||<sub>F</sub>表示矩阵的弗罗贝尼乌斯Frobenius范数,δ是一个常数,表示高斯噪声的标准偏差,λ是正则化参数;步骤3,求解高光谱图像鲁棒低秩张量RLRTR去噪优化模型,获得降噪后的高光谱图像。
地址 430072 湖北省武汉市武昌区珞珈山武汉大学