发明名称 确定图像中斑点噪声强度的方法
摘要 一种确定图像中斑点噪声强度的方法,包括选取灰度相对均匀区域、计算相对均匀区域中象素的灰度均值和各象素不同阶高斯厄米特矩、构造特征向量、计算检测图像的噪声强度特征值、转换噪声强度特征值、确定图像中斑点噪声的强度值6个步骤。本发明直接以乘性噪声模型为基础,在没有任何先验知识的情况下,进行噪声强度的确定,具有精度高、速度快、实用性和通用性强的特点,可用于具有灰度均匀区的可见光图像、合成孔径雷达图像以及医学超声图像等含斑点噪声的图像噪声强度的确定。
申请公布号 CN101807298A 申请公布日期 2010.08.18
申请号 CN201010100890.5 申请日期 2010.01.22
申请人 陕西师范大学 发明人 马苗;丁生荣;张艳宁;郭敏
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 西安永生专利代理有限责任公司 61201 代理人 申忠才
主权项 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>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mo>-</mo><mn>1</mn></mrow><mn>1</mn></munderover><munderover><mi>&Sigma;</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>&sigma;</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>&sigma;</mi><mo>,</mo><mi>v</mi><mo>/</mo><mi>&sigma;</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>&sigma;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msup><mi>&pi;&sigma;</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>&sigma;</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>&lambda;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>&lambda;</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>&lambda;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>&lambda;</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>&Sigma;</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)转换后的噪声强度特征值。
地址 710062 陕西省西安市长安南路199号