发明名称 基于对数图像处理模型的改进的反锐化掩模图像增强方法
摘要 一种基于对数图像处理模型的改进的反锐化掩模图像增强方法,属于图像处理领域,方法为:将在传统模型下基于灰度的输入图像转换成对数图像处理模型下基于灰色调的图像;对图像进行改进的反锐化掩模增强:将反锐化掩模增强图像转换为传统模型下基于灰度的图像;本发明在对数图像处理模型下实现反锐化掩模图像增强方法,克服了传统模型下反锐化方法出现灰度“溢出”的不足;本发明利用对数图像处理模型的操作与传统处理模型的操作间的同构关系,通过同构变换和反同构变换,实现更简单、高效;本发明利用图像的梯度信息对图像的边缘区域和平滑区域实施不同程度增强,弥补了传统方法利用常数因子对整个图像进行相同强度增强的不足,所得结果更好。
申请公布号 CN102096909B 申请公布日期 2012.07.04
申请号 CN201010613607.9 申请日期 2010.12.30
申请人 东北大学 发明人 姜慧研;冯锐杰;高熙和
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 沈阳东大专利代理有限公司 21109 代理人 梁焱
主权项 1.一种基于对数图像处理模型的改进的反锐化掩模图像增强方法,其特征在于:包括以下步骤:步骤1:将在传统模型下基于灰度的输入图像转换成对数图像处理模型下基于灰色调的图像,公式为:<maths num="0001"><![CDATA[<math><mrow><mi>I</mi><mo>=</mo><msub><mi>M</mi><mn>1</mn></msub><mo>-</mo><mover><mi>I</mi><mo>&OverBar;</mo></mover></mrow></math>]]></maths>式中,I表示灰色调图像,M<sub>1</sub>表示最大灰度值矩阵,<img file="FDA0000133939760000012.GIF" wi="31" he="42" />表示灰度图像;步骤2:对图像进行改进的反锐化掩模增强,方法为:步骤2-1:计算图像的梯度因子g,首先利用高斯平滑算子对传统模型下基于灰度的图像进行滤波,公式为:I<sub>s</sub>=I*G<sub>σ</sub>式中,I<sub>s</sub>为对传统模型下基于灰度的输入图像进行高斯滤波以后的图像,G<sub>σ</sub>为方差为σ的高斯滤波器;然后,利用差分方法计算图像的梯度信息Gray;最后,计算图像的梯度因子g,公式为:<maths num="0002"><![CDATA[<math><mrow><mi>g</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><mi>Gray</mi><mo>/</mo><mi>m</mi></mrow></mfrac></mrow></math>]]></maths>对图像进行反锐化掩模增强,公式为:I<sub>um</sub>=I+g×(I-I<sub>us</sub>)式中,I<sub>us</sub>表示经低通滤波后的图像,公式为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>I</mi><mi>us</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>M</mi><mo>&times;</mo><mi>N</mi></mrow></mfrac><mo>&times;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>x</mi><mo>-</mo><mrow><mo>(</mo><mi>M</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mrow><mrow><mi>x</mi><mo>+</mo><mrow><mo>(</mo><mi>M</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>y</mi><mo>-</mo><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mrow><mrow><mi>y</mi><mo>+</mo><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mrow></munderover><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></math>]]></maths>式中,N,M为图像的高和宽,I(i,j)表示原始图像,x,y分别表示当前所处理图像像素点的横纵坐标;最后,利用反同构变换对上述图像进行处理,公式如下:<maths num="0004"><![CDATA[<math><mrow><mi>I</mi><mo>=</mo><msup><mi>&phi;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mover><mi>I</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>=</mo><msub><mi>M</mi><mn>1</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mover><mi>I</mi><mo>&OverBar;</mo></mover><mo>/</mo><msub><mi>M</mi><mn>1</mn></msub></mrow></msup><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA0000133939760000016.GIF" wi="127" he="56" />表示对灰色调图像进行反同构变换后的图像;步骤2-2:对图像进行改进的反锐化掩模增强:首先,根据同构变换对灰色调图像进行处理,变换公式为:<maths num="0005"><![CDATA[<math><mrow><mover><mi>I</mi><mo>&OverBar;</mo></mover><mo>=</mo><mi>&phi;</mi><mrow><mo>(</mo><mi>I</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msub><mi>M</mi><mn>1</mn></msub><mi>ln</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mi>I</mi><msub><mi>M</mi><mn>1</mn></msub></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>式中,φ(I)表示对灰色调图像进行同构变换后的图像;然后,利用均值滤波器对传统模型下基于灰度的输入图像进行低通滤波处理,公式为:<maths num="0006"><![CDATA[<math><mrow><msub><mi>I</mi><mi>us</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>M</mi><mo>&times;</mo><mi>N</mi></mrow></mfrac><mo>&times;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>x</mi><mo>-</mo><mrow><mo>(</mo><mi>M</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mrow><mrow><mi>x</mi><mo>+</mo><mrow><mo>(</mo><mi>M</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>y</mi><mo>-</mo><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mrow><mrow><mi>y</mi><mo>+</mo><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mrow></munderover><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></math>]]></maths>式中,N,M为图像的高和宽,I(i,j)表示原始图像;之后,利用边缘信息作为权重因子对原图像进行反锐化掩模增强处理,此时的操作运算为传统的线性运算,其公式为:I<sub>um</sub>=I+g×(I-I<sub>us</sub>)最后,利用反同构变换对上述图像进行处理,公式如下:<maths num="0007"><![CDATA[<math><mrow><mi>I</mi><mo>=</mo><msup><mi>&phi;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mover><mi>I</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>=</mo><msub><mi>M</mi><mn>1</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mover><mi>I</mi><mo>&OverBar;</mo></mover><msub><mrow><mo>/</mo><mi>M</mi></mrow><mn>1</mn></msub></mrow></msup><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA0000133939760000024.GIF" wi="127" he="56" />表示对灰色调图像进行反同构变换后的图像;所述的步骤2-1和步骤2-2同时执行;步骤3:将步骤2得到的反锐化掩模增强图像转换为传统模型下基于灰度的图像,公式为:<maths num="0008"><![CDATA[<math><mrow><msub><mover><mi>I</mi><mo>&OverBar;</mo></mover><mi>um</mi></msub><mo>=</mo><msub><mi>M</mi><mn>1</mn></msub><mo>-</mo><msub><mi>I</mi><mi>um</mi></msub></mrow></math>]]></maths>式中,<img file="FDA0000133939760000026.GIF" wi="57" he="58" />表示增强结果的灰度图像,M<sub>1</sub>表示最大灰度值矩阵,该矩阵与图像矩阵大小相同而且每一位的值都是图像的最大灰度值,I<sub>um</sub>表示增强结果的灰色调图像;
地址 110819 辽宁省沈阳市和平区文化路3号巷11号