发明名称 一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器
摘要 本发明所提出的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器是基于一种特殊分数阶偏微分方程去噪算法来实现对纹理图像的分数阶、非线性、多尺度、高精度去噪。该滤波器是采用缓存器、差值平方器一、-v<sub>3</sub>次方幂方器、乘法器一、λ<sup>n</sup>发生器、<img file="DSA0000092596410000011.GIF" wi="146" he="64" />发生器、乘法器二、乘法器三和加法器一以级联方式构成的。该滤波器具有在去除纹理图像噪声的同时,既能尽量保留平滑区域中的低频轮廓,同时又能非线性保留灰度值跃变幅度相对较大的高频边缘,而且还能非线性保留灰度值跃变幅度变化相对不大的纹理细节的显著优点。该滤波器特别适用于对富含复杂纹理细节特征的图像进行去噪的应用场合。本发明属于应用数学、数字图像处理和数字电路交叉学科的技术领域。
申请公布号 CN103337060B 申请公布日期 2015.09.23
申请号 CN201310298382.6 申请日期 2013.07.17
申请人 蒲亦非 发明人 蒲亦非
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 代理人
主权项 一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器,其特征在于:它是由缓存器(2)、差值平方器一(3)、‑v<sub>3</sub>次方幂方器(4)、乘法器一(5)、λ<sup>n</sup>发生器(11)完成的计算是<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mi>&lambda;</mi><mi>n</mi></msup><mo>=</mo><mfrac><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>v</mi><mn>3</mn></msub><mo>)</mo></mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mn>2</mn><mo>-</mo><msub><mi>v</mi><mn>3</mn></msub><mo>)</mo></mrow></mrow><mrow><msup><mi>&sigma;</mi><msup><mi>n</mi><mn>2</mn></msup></msup><msubsup><mi>s</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>n</mi></msubsup></mrow></mfrac><munder><mi>&Sigma;</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></munder><mi>P</mi><mrow><mo>(</mo><msubsup><mi>s</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>n</mi></msubsup><mo>)</mo></mrow><msup><mrow><mo>(</mo><msubsup><mi>s</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mn>0</mn></msubsup><mo>-</mo><msubsup><mi>s</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>n</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow>]]></math><img file="FSB0000141793410000011.GIF" wi="1264" he="200" /></maths><img file="FSB00001417934100000119.GIF" wi="178" he="97" />发生器(10)完成的计算是<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><msubsup><mi>s</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>n</mi></msubsup><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>-</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>v</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>v</mi><mn>1</mn></msub><mo>)</mo></mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>v</mi><mn>3</mn></msub><mo>)</mo></mrow></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mn>1</mn></munderover><mfrac><mrow><munderover><mi>&Pi;</mi><mrow><mi>&tau;</mi><mo>=</mo><mn>1</mn></mrow><mrow><mn>2</mn><mi>k</mi></mrow></munderover><mrow><mo>(</mo><msub><mi>v</mi><mn>2</mn></msub><mo>-</mo><mi>&tau;</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><mn>2</mn><mi>k</mi><mo>)</mo></mrow><mo>!</mo></mrow></mfrac><mo>[</mo><msubsup><mi>D</mi><mi>x</mi><mn>1</mn></msubsup><mrow><mo>(</mo><msup><mrow><mo>|</mo><msup><mi>D</mi><msub><mi>v</mi><mn>1</mn></msub></msup><msubsup><mi>s</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>n</mi></msubsup><mo>|</mo></mrow><mrow><msub><mi>v</mi><mn>2</mn></msub><mo>-</mo><mn>2</mn><mi>k</mi><mo>-</mo><mn>2</mn></mrow></msup><msubsup><mi>D</mi><mi>x</mi><msub><mi>v</mi><mn>1</mn></msub></msubsup><msubsup><mi>s</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>n</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msubsup><mi>D</mi><mi>y</mi><mn>1</mn></msubsup><mrow><mo>(</mo><msup><mrow><mo>|</mo><msup><mi>D</mi><msub><mi>v</mi><mn>1</mn></msub></msup><msubsup><mi>s</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>n</mi></msubsup><mo>|</mo></mrow><mrow><msub><mi>v</mi><mn>2</mn></msub><mo>-</mo><mn>2</mn><mi>k</mi><mo>-</mo><mn>2</mn></mrow></msup><msubsup><mi>D</mi><mi>y</mi><msub><mi>v</mi><mn>1</mn></msub></msubsup><msubsup><mi>s</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>n</mi></msubsup><mo>)</mo></mrow><mo>]</mo><mo>,</mo></mrow>]]></math><img file="FSB0000141793410000012.GIF" wi="1969" he="229" /></maths>乘法器二(8)、乘法器三(9)和加法器一(6)以级联方式构成的,其中,阶次v<sub>1</sub>、v<sub>2</sub>和v<sub>3</sub>取分数或有理小数,且v<sub>3</sub>≠1,2,3,s代表数字图像,x代表x轴上的变量,y代表y轴上的变量,<img file="FSB0000141793410000013.GIF" wi="97" he="116" />代表在数字图像s中对于变量x与y求和,n取任意正整数,λ为正则化参数,整数k取0和1,参数τ取1到2k之间的整数,第0次迭代的原始数字图像为<img file="FSB0000141793410000014.GIF" wi="109" he="81" />第n次迭代的数字图像为<img file="FSB0000141793410000015.GIF" wi="109" he="84" />Gamma函数为Γ,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mi>&sigma;</mi><msup><mi>n</mi><mn>2</mn></msup></msup><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></munder><msup><mrow><mo>(</mo><msubsup><mi>s</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mn>0</mn></msubsup><mo>-</mo><msubsup><mi>s</mi><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mi>n</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow>]]></math><img file="FSB0000141793410000016.GIF" wi="600" he="195" /></maths>方差初始值<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msup><mi>&sigma;</mi><msup><mn>1</mn><mn>2</mn></msup></msup><mo>=</mo><mn>0.01</mn><mo>,</mo></mrow>]]></math><img file="FSB00001417934100000120.GIF" wi="299" he="111" /></maths><img file="FSB0000141793410000017.GIF" wi="77" he="67" />完成的计算是在x轴方向上的v<sub>1</sub>阶分数阶微分,<img file="FSB0000141793410000018.GIF" wi="77" he="57" />完成的计算是在x轴方向上和y轴方向上同时进行v<sub>1</sub>阶分数阶微分,<img file="FSB0000141793410000019.GIF" wi="94" he="74" />完成的计算是在y轴方向上的v<sub>1</sub>阶分数阶微分,<img file="FSB00001417934100000110.GIF" wi="86" he="93" />完成的计算是在x轴方向上的1阶微分,<img file="FSB00001417934100000111.GIF" wi="71" he="92" />完成的计算是在y轴方向上的1阶微分,参数μ=0.005,单位迭代时间间隔Δt在(0,0.1]内取任意的正实数。
地址 610066 四川省成都市锦江区沙河堡菱窠东路122号16栋1单元11号