主权项 |
1.一种保持边缘的自适应图像滤波方法,其特征在于:a)使用极值法检测噪声图像中被椒盐噪声所污染的像素;b)对被污染的像素使用自适应滤波法修正其灰度值,而未受噪声污染的像素则保持其灰度值不变,得到去除椒盐噪声后的图像;c)对步骤b得到的图像进行平稳小波分解,得到低频分量和不同频带、不同方向的高频分量,低频分量及高频分量的大小与原始噪声图像的大小相同;d)由于小波分解后的低频分量已经比较平滑,故保持其系数值不变;对于高频分量中的噪声和边缘,利用边缘在不同尺度的对应位置上具有较强的相关性,而噪声的相关性则很弱的特性,将高频分量中的像素标记为边缘或噪声;e)如果高频分量的某一像素被标记为边缘,则保持其系数值不变;如果被标记为噪声,则采用以下步骤对其系数值进行收缩:(1)对某一像素(x,y),求以其为中心的某一大小邻域内像素小波系数Dji的平方和<maths id="math0001" num="0001" ><math><![CDATA[ <mrow> <msubsup> <mi>S</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mn>2</mn> </msubsup> <mo>=</mo> <munder> <mi>Σ</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>∈</mo> <mi>ϵ</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> </munder> <msup> <mrow> <mo>[</mo> <msubsup> <mi>D</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>]</mo> </mrow> <mn>2</mn> </msup> </mrow>]]></math></maths>其中(m,n)∈ε(x,y)表示(m,n)是以(x,y)为中心的邻域内的点;(2)计算如下三式的值,得到收缩后的小波系数值<maths id="math0002" num="0002" ><math><![CDATA[ <mrow> <mi>λ</mi> <mo>=</mo> <msqrt> <mn>2</mn> <msubsup> <mi>σ</mi> <mi>n</mi> <mn>2</mn> </msubsup> <mi>log</mi> <mi>L</mi> </msqrt> </mrow>]]></math></maths><maths id="math0003" num="0003" ><math><![CDATA[ <mrow> <mi>β</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <mrow> <mo>(</mo> <msup> <mi>λ</mi> <mn>2</mn> </msup> <mo>/</mo> <msubsup> <mi>S</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> </mtd> <mtd> <msup> <mi>λ</mi> <mn>2</mn> </msup> <mo>/</mo> <msubsup> <mi>S</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mn>2</mn> </msubsup> <mo>≤</mo> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msup> <mi>λ</mi> <mn>2</mn> </msup> <mo>/</mo> <msubsup> <mi>S</mi> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mn>2</mn> </msubsup> <mo>></mo> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>]]></math></maths><maths id="math0004" num="0004" ><math><![CDATA[ <mrow> <mi>new</mi> <msubsup> <mi>D</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>β</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <msubsup> <mi>D</mi> <mi>j</mi> <mi>i</mi> </msubsup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow>]]></math></maths>式中log是对数函数,σn2是噪声方差,L表示图像的大小;f)当噪声强度较大时,用上述步骤e收缩后最小尺度的高频分量中会出现一些孤立的亮点和暗点,但是次小尺度即最小尺度的上一层高频分量中的噪声已经被去除,则采用以下步骤将这些孤立点去除:(1)分别计算次小尺度三个方向的高频分量收缩后系数绝对值的均值mean;(2)设mask为一个标记矩阵,其大小等于次小尺度高频分量图像的大小,如果次小尺度高频分量中某一像素(i,j)小波系数的绝对值小于mean,则置mask(i,j)的值为1;(3)由于平稳小波分解后各尺度分量的大小等于原始图像的大小,因此相邻两层之间很容易找到对应的像素。如果mask(i,j)等于1,则将最小尺度高频分量各小波系数D1(i,j)的值置为0;g)对经过上述处理的高频分量和低频分量进行平稳小波重构得到滤波图像。 |