发明名称 一种图像的同步去噪增强方法
摘要 图像获取过程中,由于受到环境、设备等客观因素的制约,常得到含有噪声、对比度差的图像,这些对于后续图像分割与配准的准确性方面都产生较大影响。对于解决去噪与增强的问题,通常的作法是先去除噪声、后增强图像,或者先增强图像、后去除噪声。然而由于去噪与增强矛盾性的存在,导致大多算法都不能很好地解决这一问题,为此提出一种图像的同步去噪增强模型。首先引入RGB与Ycbcr的变换及逆变换,构造了适用于彩色图像的自适应直方图均衡化偏微分方程。其次,通过计算本征矢与x轴的夹角统一梯度变化率、消除法相扩散,给出了彩色图像的PM模型。最后将上述两种模型通过增加系数调节项的方式进行融合,达到同步去噪增强的目的。
申请公布号 CN103279930B 申请公布日期 2015.10.21
申请号 CN201310201177.3 申请日期 2013.05.27
申请人 辽宁工程技术大学 发明人 贾迪;刘影;王伟;董娜;齐浩然
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 代理人
主权项 一种图像的同步去噪增强方法,采用偏微分方程方式,其特征在于通过迭代完成图像的同步去噪与增强,整体处理流程包括:令<img file="dest_path_image001.GIF" wi="20" he="25" />为输入影像,<img file="dest_path_image002.GIF" wi="25" he="25" />为输出结果,<img file="dest_path_image003.GIF" wi="14" he="16" />为迭代次数,步长为<img file="dest_path_image004.GIF" wi="21" he="20" />,切向平滑系数为<img file="dest_path_image005.GIF" wi="22" he="25" />,增强比例系数为<img file="dest_path_image006.GIF" wi="24" he="26" />,算法流程如下:(1) 令<img file="dest_path_image007.GIF" wi="57" he="25" />;(2) 设定滑动窗口大小,计算<img file="dest_path_image008.GIF" wi="93" he="25" />;(3) 将<img file="dest_path_image009.GIF" wi="49" he="25" />与<img file="243160dest_path_image006.GIF" wi="24" he="26" />相乘得到<img file="dest_path_image010.GIF" wi="26" he="25" />;(4) 计算去噪微分项<img file="dest_path_image011.GIF" wi="22" he="21" />,再与<img file="799781dest_path_image005.GIF" wi="22" he="25" />相乘得到<img file="dest_path_image012.GIF" wi="25" he="26" />;(5) 计算<img file="dest_path_image013.GIF" wi="162" he="26" />;(6) 统计迭代次数<img file="dest_path_image014.GIF" wi="18" he="16" />,如果<img file="dest_path_image015.GIF" wi="42" he="16" />转步骤(2);(7) 获得结果<img file="907415dest_path_image002.GIF" wi="25" he="25" />;其中<img file="dest_path_image016.GIF" wi="20" he="25" />的计算过程如下:可将彩色图像看做3维矢量图像<img file="dest_path_image017.GIF" wi="48" he="22" />,根据黎曼几何观点,将其看做Euclidean空间中的一个以<img file="dest_path_image018.GIF" wi="40" he="22" />为参数的超曲面,设<img file="dest_path_image019.GIF" wi="21" he="20" />为曲面<img file="817471dest_path_image017.GIF" wi="48" he="22" />上给定方向的弧长微元,则表达式如下:<img file="dest_path_image020.GIF" wi="379" he="70" /><img file="dest_path_image021.GIF" wi="409" he="45" />其中<img file="dest_path_image022.GIF" wi="182" he="154" /><img file="dest_path_image023.GIF" wi="269" he="53" />则矩阵A的特征值为:<img file="dest_path_image024.GIF" wi="196" he="52" />令<img file="dest_path_image025.GIF" wi="17" he="16" />为归一化梯度矢量,<img file="dest_path_image026.GIF" wi="17" he="22" />为切矢量,给出彩色图像下的PM方程:<img file="dest_path_image027.GIF" wi="343" he="49" /><img file="dest_path_image028.GIF" wi="347" he="29" /><img file="dest_path_image029.GIF" wi="208" he="29" />其中<img file="dest_path_image030.GIF" wi="25" he="26" />为切向二阶偏导差分项,<img file="dest_path_image031.GIF" wi="25" he="28" />为法向二阶偏导差分项,<img file="dest_path_image032.GIF" wi="34" he="22" />为切向边缘停止函数,<img file="dest_path_image033.GIF" wi="36" he="22" />为法向边缘停止函数,<img file="dest_path_image034.GIF" wi="12" he="20" />为通道号,对于RGB彩色图像而言,<img file="487968dest_path_image034.GIF" wi="12" he="20" />取值为0、1、2,由于<img file="dest_path_image035.GIF" wi="53" he="26" />具有与灰度图像梯度模值相似的含义,因此可替代彩色图像分块直方图均衡化过程中使用到的梯度模值,令<img file="dest_path_image036.GIF" wi="49" he="25" />,其中<img file="920217dest_path_image002.GIF" wi="25" he="25" />在步骤(1)中获得,通过如下方程进行<img file="515147dest_path_image016.GIF" wi="20" he="25" />的迭代求解:<img file="dest_path_image037.GIF" wi="321" he="49" />其中<img file="dest_path_image038.GIF" wi="55" he="26" />是第<img file="822369dest_path_image034.GIF" wi="12" he="20" />块的最大灰度值,<img file="dest_path_image039.GIF" wi="52" he="26" />是第<img file="835456dest_path_image034.GIF" wi="12" he="20" />块最小灰度值,<img file="dest_path_image040.GIF" wi="20" he="21" />为当前时刻的图像、<img file="dest_path_image041.GIF" wi="28" he="21" />为下一时刻的图像,<img file="dest_path_image042.GIF" wi="24" he="25" />为迭代步长,<img file="dest_path_image043.GIF" wi="16" he="18" />为分块数,<img file="dest_path_image045.GIF" wi="9" he="18" />为分块编号,<img file="dest_path_image046.GIF" wi="60" he="25" />为取当前时刻图像<img file="756767dest_path_image040.GIF" wi="20" he="21" />第<img file="838993dest_path_image045.GIF" wi="9" he="18" />分块的颜色分量,<img file="dest_path_image047.GIF" wi="36" he="21" />为分块直方图均衡化函数,最后将迭代完成后的<img file="123475dest_path_image041.GIF" wi="28" he="21" />结果作为<img file="974756dest_path_image016.GIF" wi="20" he="25" />;步骤(4)中的<img file="185027dest_path_image011.GIF" wi="22" he="21" />计算过程如下:将彩色PM微分模型与AHE自适应直方图均衡化微分模型进行融合得到:<img file="dest_path_image048.GIF" wi="432" he="51" />其中<img file="dest_path_image049.GIF" wi="45" he="22" />为散度计算函数,<img file="dest_path_image050.GIF" wi="24" he="20" />为图像梯度,<img file="dest_path_image051.GIF" wi="40" he="25" />为分块的中心坐标,<img file="dest_path_image052.GIF" wi="39" he="28" />是以<img file="269395dest_path_image051.GIF" wi="40" he="25" />为分块中心坐标的最大灰度值,<img file="dest_path_image053.GIF" wi="39" he="28" />是以<img file="419885dest_path_image051.GIF" wi="40" he="25" />为分块中心坐标的最小灰度值,使用梯度下降法对其进行求解,得到如下的显式数值计算方程:<img file="dest_path_image054.GIF" wi="396" he="49" />其中,<img file="250306dest_path_image004.GIF" wi="21" he="20" />为迭代步长,<img file="851052dest_path_image011.GIF" wi="22" he="21" />与<img file="dest_path_image055.GIF" wi="22" he="25" />分别为图像的切向与法向的显示求解项,<img file="688296dest_path_image043.GIF" wi="16" he="18" />为分块数,设图像大小为<img file="dest_path_image056.GIF" wi="46" he="20" />,窗口大小为<img file="dest_path_image057.GIF" wi="32" he="20" />,并采用重叠式窗口移动方法,则<img file="642476dest_path_image043.GIF" wi="16" he="18" />取值为<img file="dest_path_image058.GIF" wi="113" he="22" />,为了增强算法的灵活性,引入三个调节参数改进方程,使方程的去噪与增强方面得到不同程度调节,扩展后的方程如下所示:<img file="dest_path_image059.GIF" wi="413" he="50" />其中<img file="576672dest_path_image005.GIF" wi="22" he="25" />为切向平滑系数,<img file="dest_path_image060.GIF" wi="21" he="26" />为法向平滑系数,<img file="426947dest_path_image006.GIF" wi="24" he="26" />为增强系数,通过这三个参数的比例设置,可以较为灵活地调节图像的去噪增强效果,由于法向扩散将使图像边缘沿着法方向扩散,从而导致弱边缘被模糊化,因此去除<img file="dest_path_image061.GIF" wi="40" he="28" />项得到下述方程:<img file="dest_path_image062.GIF" wi="343" he="46" /><img file="dest_path_image063.GIF" wi="349" he="51" />平滑系数<img file="564537dest_path_image005.GIF" wi="22" he="25" />与增强系数<img file="774938dest_path_image006.GIF" wi="24" he="26" />用于调节去噪与增强的比例,实际应用中可根据需求调节两个比例系数,当噪声较强时增加<img file="393001dest_path_image005.GIF" wi="22" he="25" />的取值,需提高对比度时增加<img file="912713dest_path_image006.GIF" wi="24" he="26" />的取值,<img file="dest_path_image064.GIF" wi="21" he="22" />为当前颜色通道下<img file="dest_path_image065.GIF" wi="18" he="17" />轴方向一阶偏导限制差分项,<img file="dest_path_image066.GIF" wi="22" he="25" />为当前颜色通道下<img file="dest_path_image067.GIF" wi="16" he="17" />轴方向一阶偏导限制差分项,<img file="dest_path_image068.GIF" wi="25" he="22" />为当前颜色通道下<img file="662232dest_path_image065.GIF" wi="18" he="17" />轴方向二阶偏导限制差分项,<img file="dest_path_image069.GIF" wi="25" he="25" />为当前颜色通道下<img file="427057dest_path_image067.GIF" wi="16" he="17" />轴方向二阶偏导限制差分项,<img file="dest_path_image070.GIF" wi="14" he="16" />用于防止分母项为0,取值通常为<img file="dest_path_image071.GIF" wi="29" he="21" />,其中<img file="dest_path_image072.GIF" wi="14" he="18" />为正整数。
地址 123000 辽宁省阜新市细河区中华路47号