主权项 |
1.一种基于小波低频的贝叶斯去噪方法,包括如下步骤:(1)输入一幅待去噪的自然图像;(2)选取待估计像素块在待去噪的自然图像中,逐行扫描选取一个像素点为待估计像素点,以待估计像素点为中心,以固定长度2~5个像素点为块半径,选取一个正方形的待估计像素块;(3)选取中心低频系数块将待去噪的自然图像进行db1小波分解,在其得到的低频图中对应于待去噪像素块的位置上,选取一个与待去噪像素块大小相同的低频系数块,以此低频系数块作为中心低频系数块;(4)确定搜索窗在低频图上,以选取的中心低频系数块的中心系数为中心,以固定长度7~15个像素点为搜索窗半径,选取一个正方形的搜索窗;(5)选取低频系数块在搜索窗中选取一个没被选取过的低频系数,以此低频系数为中心,确定一个与中心低频系数块等大小的低频系数块;(6)判断是否满足约束条件判断低频系数块是否同时满足均值约束条件和方差约束条件,若满足,则进行下一步骤,否则执行步骤(5);(7)计算相似性权值按照下式计算相似性权值:<maths num="0001"><![CDATA[<math><mrow><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>Z</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></mfrac><mi>exp</mi><mo>[</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mrow><mo>(</mo><msqrt><mn>2</mn><mi>dist</mi><mrow><mo>(</mo><msubsup><mi>W</mi><mi>i</mi><mi>k</mi></msubsup><mo>,</mo><msubsup><mi>W</mi><mi>j</mi><mi>k</mi></msubsup><mo>)</mo></mrow></msqrt><mo>-</mo><msqrt><mn>2</mn><mi>N</mi><mo>-</mo><mn>1</mn></msqrt><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo></mrow></math>]]></maths>其中,w(i,j)为相似性权值;i为待估计像素块的中心像素;j为在待去噪的自然图像中,与低频系数块的中心系数所对应的像素点;Z(i)为归一化因子;exp为指数函数;<img file="FSB00001123564200021.GIF" wi="292" he="79" />为<img file="FSB00001123564200022.GIF" wi="72" he="63" />与<img file="FSB00001123564200023.GIF" wi="72" he="69" />之间的相似度距离;<img file="FSB00001123564200024.GIF" wi="71" he="62" />为小波第k层分解的低频图中,以低频系数i为中心的中心低频系数块;<img file="FSB00001123564200025.GIF" wi="71" he="67" />为小波第k层分解的低频图中,以低频系数j为中心的低频系数块;N为待估计像素块中像素点的个数,其值为(2f+1)<sup>2</sup>,f为待估计像素块的半径;(8)判断搜索窗内的点是否搜索完在搜索窗中选取一个没被选取过的低频系数,若选取低频系数的次数大于搜索窗内点的个数,则进行下一步,否则执行步骤(5),直至获得搜索窗中所有低频系数块与中心低频系数块之间的相似性权值;(9)计算待估计像素块的恢复值使用贝叶斯非局部滤波方法中的加权求和方法,将每个低频系数块的相似性权值与待去噪自然图像中每个低频系数块所对应的像素块加权求和,得出待估计像素块的恢复值;(10)判断待去噪自然图像是否搜索完在待去噪自然图像中,选取一个没被选取过的像素点,若选取待估计像素点的次数大于图像中像素点的个数,则进行下一步,否则执行步骤(2),直至获得整幅自然图像中所有待估计像素块的恢复值;(11)整合恢复值按照下式整合所有的待估计像素块的恢复值,获得去噪后的自然图像;<maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>L</mi></mfrac><munderover><mi>Σ</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msubsup><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi><mi>l</mi></msubsup></mrow></math>]]></maths>其中,<img file="FSB00001123564200027.GIF" wi="36" he="57" />为自然图像中第i个像素点的恢复值;L为第i个像素点在步骤(10)中被重复恢复的次数;l为<img file="FSB00001123564200028.GIF" wi="36" he="58" />被恢复的次序;<img file="FSB00001123564200029.GIF" wi="39" he="62" />为步骤(10)中第l次恢复出的包含像素点i的像素块在i点处的恢复值。 |