发明名称 一种基于锐化强度与灰度分布结合的图像增强方法
摘要 一种基于锐化强度与灰度分布结合的图像增强方法,先对彩色图像的亮度信息进行拉普拉斯变换,统计图像中每个像素点的突变程度即锐化强度,其次对锐化强度进行量化,并将量化值融入到图像亮度分量的灰度直方图中对其进行修正,然后进行灰度对比度拉伸,最后完成对彩色图像的增强。本发明在增强图像细节信息的同时,对整幅图像进行对比度拉伸,进而获得更好的图像增强效果。
申请公布号 CN102194220A 申请公布日期 2011.09.21
申请号 CN201110120101.9 申请日期 2011.05.10
申请人 西安理工大学 发明人 王栋;赵文平;朱虹;康凯;刑楠;刘薇
分类号 G06T5/40(2006.01)I 主分类号 G06T5/40(2006.01)I
代理机构 西安弘理专利事务所 61214 代理人 李娜
主权项 1.一种基于锐化强度与灰度分布结合的图像增强方法,其特征在于,具体步骤如下:步骤1、确定彩色图像亮度信息设输入的彩色图像大小为m×n,该彩色图像的红、绿、蓝三色矩阵分量分别为R=[r(i,j)]<sub>m×n</sub>,G=[g(i,j)]<sub>m×n</sub>,B=[b(i,j)]<sub>m×n</sub>,则其亮度矩阵分量L=[L(i,j)]<sub>m×n</sub>由下式得到:l(i,j)=max(r(i,j),g(i,j),b(i,j)),i=1,2,...,m,j=1,2,...,n;(1)其中,r(i,j),g(i,j),b(i,j)以及l(i,j)分别为该彩色图像的红、绿、蓝三色分量,以及亮度分量在点(i,j)处像素的值;步骤2、计算锐化强度设<img file="FDA0000060266370000011.GIF" wi="380" he="56" />为采用Laplacian算子检测出的锐化强度幅值,a.初始化<maths num="0001"><![CDATA[<math><mrow><msup><mo>&dtri;</mo><mn>2</mn></msup><mi>L</mi><mo>=</mo><msub><mrow><mo>[</mo><msup><mo>&dtri;</mo><mn>2</mn></msup><mi>l</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow></msub><mo>:</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msup><mo>&dtri;</mo><mn>2</mn></msup><mi>l</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths>i=1,2,...,m,j=1,2,...,n;(2)b.求解锐化强度矩阵<maths num="0003"><![CDATA[<math><mrow><msup><mo>&dtri;</mo><mn>2</mn></msup><mi>L</mi><mo>=</mo><msub><mrow><mo>[</mo><msup><mo>&dtri;</mo><mn>2</mn></msup><mi>l</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow></msub><mo>:</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msup><mo>&dtri;</mo><mn>2</mn></msup><mi>l</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mn>4</mn><mi>l</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>l</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>l</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mi>l</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>l</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>|</mo><mo>,</mo></mrow></math>]]></maths>其中,i=2,3...,m-1,j=2,3,...,n-1;(3)步骤3、量化锐化强度a.求<img file="FDA0000060266370000016.GIF" wi="380" he="56" />的均值μ:<maths num="0005"><![CDATA[<math><mrow><mi>&mu;</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mo>&dtri;</mo><mn>2</mn></msup><mi>l</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>b.选取k<sub>μ</sub>μ为量化单位,对<img file="FDA0000060266370000018.GIF" wi="82" he="44" />进行量化,得到锐化强度矩阵W=[w(i,j)]<sub>m×n</sub>:<maths num="0006"><![CDATA[<math><mrow><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msup><mo>&dtri;</mo><mn>2</mn></msup><mi>l</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>/</mo><msub><mi>k</mi><mi>&mu;</mi></msub><mi>&mu;</mi><mo>,</mo></mrow></math>]]></maths>i=1,2,...,m,j=1,2,...,n;(5)其中,k<sub>μ</sub>为量化因子;步骤4、计算图像亮度分量的锐化强度融合灰度直方图a.设图像亮度分量L=[l(i,j)]<sub>m×n</sub>的灰度直方图为H=[h(i)]<sub>256×1</sub>,并初始化H=[h(i)]<sub>256×1</sub>:h(k)=0,k=0,1,...,255;(6)b.求解锐化强度融合灰度直方图:k=l(i,j),i=1,2,...,m,j=1,2,...,n;(7)h(k)=h(k)+1+w(i,j);         (8)步骤5、增强亮度分量a.计算步骤4得到的锐化强度融合灰度直方图的概率分布对于步骤3得到的锐化强度矩阵W=[w(i,j)]<sub>m×n</sub>,计算其总强度值N<sub>w</sub>:<maths num="0007"><![CDATA[<math><mrow><msub><mi>N</mi><mi>w</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>w</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>i=1,2,...,m,j=1,2,...,n;(9)则锐化强度融合灰度的概率分布为P<sub>H</sub>=[p<sub>h</sub>(i)]<sub>256×1</sub>为:P<sub>H</sub>(k)=h(k)/(m×n+N<sub>w</sub>);           (10)b.计算锐化强度融合灰度直方图的累积分布设该累积分布为P<sub>HS</sub>=[p<sub>hs</sub>(i)]<sub>256×1</sub>,则<maths num="0008"><![CDATA[<math><mrow><msub><mi>P</mi><mi>HS</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>k</mi></munderover><msub><mi>p</mi><mi>h</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>k=0,1,...,255;(11)c.设对亮度分量进行增强处理后的结果为L*=[l*(i,j)]<sub>m×n</sub>,则k=l(i,j),i=1,2,...,m,j=1,2,...,n;(12)l*(i,j)=255·P<sub>HS</sub>(k);           (13)步骤6、对彩色图像进行增强a.对输入的彩色图像的红、绿、蓝三个颜色分量计算其颜色比例系数矩阵A<sub>R</sub>=[α<sub>R</sub>(i,j)]<sub>m×n</sub>,A<sub>G</sub>=[α<sub>G</sub>(i,j)]<sub>m×n</sub>,A<sub>B</sub>=[α<sub>B</sub>(i,j)]<sub>m×n</sub>:<maths num="0009"><![CDATA[<math><mrow><mover><mi>l</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>l</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><mn>1</mn><mo>,</mo></mrow></math>]]></maths>i=1,2,...,m,j=1,2,...,n;(14)<maths num="0010"><![CDATA[<math><mrow><msub><mi>&alpha;</mi><mi>R</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>r</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>/</mo><mover><mi>l</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>i=1,2,...,m,j=1,2,...,n;(15)<maths num="0011"><![CDATA[<math><mrow><msub><mi>&alpha;</mi><mi>G</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>g</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>/</mo><mover><mi>l</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>i=1,2,...,m,j=1,2,...,n;(16)<maths num="0012"><![CDATA[<math><mrow><msub><mi>&alpha;</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>b</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>/</mo><mover><mi>l</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>i=1,2,...,m,j=1,2,...,n;(17)b.求增强后的彩色图像设对彩色图像进行增强处理后的红、绿、蓝三个颜色分量为R*=[r*(i,j)]<sub>m×n</sub>,G*=[g*(i,j)]<sub>m×n</sub>,B*=[b*(i,j)]<sub>m×n</sub>,则:r*(i,j)=α<sub>R</sub>(i,j)·l*(i,j),i=1,2,...,m,j=1,2,...,n;(18)g*(i,j)=α<sub>G</sub>(i,j)·l*(i,j),i=1,2,...,m,j=1,2,...,n;(19)b*(i,j)=α<sub>B</sub>(i,j)·l*(i,j),i=1,2,...,m,j=1,2,...,n;(20)即得到增强后的彩色图像。
地址 710048 陕西省西安市金花南路5号