图像锐化主要影响图像中的低频分量,不影响图像中的高频分量。

图像锐化的主要目的有两个:

1.增强图像边缘,使模糊的图像变得更加清晰,颜色变得鲜明突出,图像的质量有所改善,产生更适合人眼观察和识别的图像;

2.希望通过锐化处理后,目标物体的边缘鲜明,以便于提取目标的边缘、对图像进行分割、目标区域识别、区域形状提取等,进一步的图像理解与分析奠定基础。

图像锐化一般有两种方法:

1.微分法

2.高通滤波法

这里主要介绍一下两种常用的微分锐化方法:梯度锐化和拉普拉斯锐化

注意:由于锐化使噪声受到比信号还要强的增强,所以要求锐化处理的图像有较高的信噪比;否则,锐化后的图像的信噪比更低。

1.梯度锐化

基本理论

邻域平均法或加权平均法可以平滑图像,反过来利用对应的微分算法可以锐化图像。微分算法是求信号的变化率,有加强高频分量胡作用,从而使图像轮廓清晰。

由于图像模糊胡实质是图像受到平均或积分运算造成的,所以为了把图像中任何方向伸展的边缘肯模糊的轮廓变得清晰,可以对图像进行逆运算如微分运算,从而使图像清晰化。

在图像处理中,一阶微分是通过梯度算法来实现的,对于一幅图像用函数f(x,y)表示,定义在f(x,y)在点(x,y)处的梯度是一个矢量,定义为:

图像基础(1)图像锐化

梯度的方向在函数f(x,y)最大变化率的方向上,梯度的幅度G[f(x,y)]可以由以下公式算出:

图像基础(1)图像锐化

由上式可知:梯度的数值就是f(x,y)在其最大变化率方向上的单位距离所增加的量。

对于数字图像而言,微分可用差分来近似。因此上式可以写成:

图像基础(1)图像锐化

为了便于编程和提高运算速度,在计算精度允许的情况下,可采用绝对差算法近似为:

图像基础(1)图像锐化

这种算法又称为水平垂直差分法,另一种梯度算法是交叉的进行查分计算,称为罗伯特梯度法,表示为:

图像基础(1)图像锐化

同样可采用绝对差算法近似:

图像基础(1)图像锐化

运用上述两种梯度近似算法,在图像的最后一行后最后一列无法计算像素的梯度时,一般用前一行或前一列的梯度值近似代替。

为了不破坏图像背景的前提下更好地增强边缘,也可以对上述直接用梯度值代替灰度值的方法进行改进,可以引入一个阈值来判断是否对某一像素点进行锐化。具体公式如下:

图像基础(1)图像锐化

对于图像而言,物体与物体之间,背景与背景之间的梯度变化很小,灰度变化较大的地方一般集中在图像的边缘上,也就是物体和背景交接的地方。当我们设定一个阈值时,G[f(i,j)]大于阈值就认为该像素点处于图像的边缘,对结果加上常数C,以使边缘变亮;而对于G[f(i,j)]不大于阈值就认为该像素点为同类像素,即同为物体或同为背景,常数C的选取可以根据具体的图像特点。这样既增亮了图像的边界,同时又保留了图像背景原来的状态,比传统的梯度锐化具有更好的增强效果和适用性。

另外拉普拉斯算子也可用于图像锐化,基于拉普拉斯算子的图像锐化

 

  • 对于求一个锐化后的像素点(sharpened_pixel),这个基于拉普拉斯算子的简单算法主要是遍历图像中的像素点,根据领域像素确定其锐化后的值
  • 计算公式:sharpened_pixel = 5 * current – left – right – up – down ;

见下图:

图像基础(1)图像锐化

当一个运算是通过领域像素进行的时候,我们通常用一个矩阵来表示这种运算关系,也就是我们经常所说的 核 (Kernel) 。那么上面的 锐化滤波器 (Sharpening Filter) 就可以用这个矩阵表示为它的核:

图像基础(1)图像锐化

因为 滤波 在图像处理中是一个非常普通且常用的操作,所以OpenCV里面已经定义了一个特殊的函数用来执行这个操作。要使用它的话只需要定义一个 核 ,然后作为参数传递就行了。

 

 


[cpp] view plaincopy

  1. //手动实现拉普拉斯算子图像锐化  
  2. void sharpenImage1(const Mat &image, Mat &result)  
  3. {  
  4.     result.create(image.size(),image.type());//为输出图像分配内容  
  5.     /*拉普拉斯滤波核3*3 
  6.          0  -1   0 
  7.         -1   5  -1 
  8.          0  -1   0  */  
  9.     //处理除最外围一圈外的所有像素值  
  10.     for(int i=1; i<image.rows-1; i++)  
  11.     {  
  12.         const uchar * pre = image.ptr<const uchar>(i-1);//前一行  
  13.         const uchar * cur = image.ptr<const uchar>(i);//当前行,第i行  
  14.         const uchar * next = image.ptr<const uchar>(i+1);//下一行  
  15.         uchar * output = result.ptr<uchar>(i);//输出图像的第i行  
  16.         int ch = image.channels();//通道个数  
  17.         int startCol = ch;//每一行的开始处理点  
  18.         int endCol = (image.cols-1)* ch;//每一行的处理结束点  
  19.         for(int j=startCol; j < endCol; j++)  
  20.         {  
  21.             //输出图像的遍历指针与当前行的指针同步递增, 以每行的每一个像素点的每一个通道值为一个递增量, 因为要  
  22.   
  23. 考虑到图像的通道数  
  24.             //saturate_cast<uchar>保证结果在uchar范围内  
  25.             *output++ = saturate_cast<uchar>(5*cur[j]-pre[j]-next[j]-cur[j-ch]-cur[j+ch]);  
  26.         }  
  27.     }  
  28.     //将最外围一圈的像素值设为0  
  29.     result.row(0).setTo(Scalar(0));  
  30.     result.row(result.rows-1).setTo(Scalar(0));  
  31.     result.col(0).setTo(Scalar(0));  
  32.     result.col(result.cols-1).setTo(Scalar(0));  
  33.     /*/或者也可以尝试将最外围一圈设置为原图的像素值 
  34.     image.row(0).copyTo(result.row(0)); 
  35.     image.row(image.rows-1).copyTo(result.row(result.rows-1)); 
  36.     image.col(0).copyTo(result.col(0)); 
  37.     image.col(image.cols-1).copyTo(result.col(result.cols-1));*/  
  38. }  
  39.   
  40. //调用OpenCV函数实现拉普拉斯算子图像锐化  
  41. void sharpenImage2(const Mat &image, Mat &result)  
  42. {  
  43.     Mat kernel = (Mat_<float>(3,3) << 0,-1,0,-1,5,-1,0,-1,0);  
  44.     filter2D(image,result,image.depth(),kernel);  
  45. }  
  46.   
  47. //main函数  
  48. void main(){      
  49.     Mat mat = imread("images/lena.jpg");  
  50.     Mat result1;  
  51.     Mat result2;  
  52.   
  53.     sharpenImage1(mat,result1);  
  54.     sharpenImage2(mat,result2);  
  55.   
  56.     namedWindow("src");  
  57.     namedWindow("result1");  
  58.     namedWindow("result2");  
  59.     imshow("src",mat);  
  60.     imshow("result1",result1);  
  61.     imshow("result2",result2);  
  62. }  



 

 

结果:

图像基础(1)图像锐化

原图

 

图像基础(1)图像锐化

手动实现

 

图像基础(1)图像锐化

调用OpenCV中函数实现

 

参考文章地址:

http://ggicci.blog.163.com/blog/static/210364096201262123236955/

http://www.cnblogs.com/liu-jun/archive/2012/08/12/2635373.html

分类:

技术点:

相关文章: