主权项 |
1.一种确定图像中斑点噪声强度的方法,其特征在于它包括下述步骤:(1)选取灰度相对均匀区域在含噪图像中选取一块灰度相对均匀的区域;(2)计算灰度相对均匀区域象素的灰度均值和该区域内象素的不同阶高斯厄米特矩;上述的灰度均值是灰度相对均匀区域内各象素点的灰度平均值;灰度相对均匀区域象素的不同阶高斯厄米特矩按下式计算:<maths num="0001"><![CDATA[<math><mrow><msub><mi>M</mi><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>t</mi><mo>=</mo><mo>-</mo><mn>1</mn></mrow><mn>1</mn></munderover><munderover><mi>Σ</mi><mrow><mi>v</mi><mo>=</mo><mo>-</mo><mn>1</mn></mrow><mn>1</mn></munderover><mi>G</mi><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>v</mi><mo>,</mo><mi>σ</mi><mo>)</mo></mrow><msub><mi>H</mi><mrow><mi>P</mi><mo>,</mo><mi>q</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>/</mo><mi>σ</mi><mo>,</mo><mi>v</mi><mo>/</mo><mi>σ</mi><mo>)</mo></mrow><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mi>t</mi><mo>,</mo><mi>y</mi><mo>+</mo><mi>v</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(1)中,G(t,v,σ)为二维高斯函数,且有<maths num="0002"><![CDATA[<math><mrow><mi>G</mi><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>v</mi><mo>,</mo><mi>σ</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msup><mi>πσ</mi><mn>2</mn></msup></mrow></mfrac><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mrow><mo>(</mo><msup><mi>t</mi><mn>2</mn></msup><mo>+</mo><msup><mi>v</mi><mn>2</mn></msup><mo>)</mo></mrow><msup><mrow><mn>2</mn><mi>σ</mi></mrow><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>H<sub>P,q</sub>(t/σ,v/σ)为二维(p,q)阶厄米特多项式,且有H<sub>P,q</sub>(t/σ,v/σ)=H<sub>p</sub>(t/σ)H<sub>q</sub>(v/σ) (3)H<sub>n</sub>(t)=(-1)<sup>n</sup>exp(t<sup>2</sup>)(d<sup>n</sup>/dtn)exp(-t<sup>2</sup>) (4)式(1)~(4)中t为-1或0或1,v为-1或0或1,σ为0.7;(3)构造特征向量根据式(1)计算的不同阶高斯厄米特矩,通过下式构造特征向量:<maths num="0003"><![CDATA[<math><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>M</mi><mi>u</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>M</mi><mi>v</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>λM</mi><mn>1,0</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>λ</mi><mo>)</mo></mrow><msub><mi>M</mi><mn>3.0</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>λM</mi><mn>0.1</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>λ</mi><mo>)</mo></mrow><msub><mi>M</mi><mn>0.3</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(5)中,λ是不同阶的高斯厄米特矩的联合权重系数,0<λ<1,对图像中所选均匀区域的每一个点(x,y),通过式(5)可得到一个特征向量[M<sub>u</sub>,M<sub>v</sub>]<sup>T</sup>;(4)计算检测图像的噪声强度特征值按式(6)计算噪声强度特征值:<maths num="0004"><![CDATA[<math><mrow><msub><mi>M</mi><mi>uv</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>N</mi></mrow></mfrac><mi>Σ</mi><mrow><mo>(</mo><mo>|</mo><msub><mi>M</mi><mi>u</mi></msub><mo>|</mo><mo>+</mo><mo>|</mo><msub><mi>M</mi><mi>v</mi></msub><mo>|</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中M<sub>uv</sub>为噪声强度特征值,N是特征向量的个数,N为正整数;(5)转换噪声强度特征值将步骤2的灰度相对均匀区域内各象素点的灰度平均值和步骤4的检测图像噪声强度特征值代入式(7):<img file="FSA00000006267100021.GIF" wi="1444" he="118" />转换噪声强度特征值;(6)确定图像中斑点噪声的强度值按下式确定图像中斑点噪声的强度值:y=99.1468x<sup>4</sup>-32.0896x<sup>3</sup>+51.8903x<sup>2</sup>-0.1269x+0.0013 (8)式中x为步骤(5)转换后的噪声强度特征值。 |