图像增强-图像锐化

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

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

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

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

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

1.微分法

2.高通滤波法

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

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

 

1.梯度锐化

基本理论

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

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

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

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

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

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

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

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

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

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

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

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

另外拉普拉斯算子也可用于图像锐化，这里不再赘述。

算法实现：

 

 1 /*************************************************************************
 2 * 函数名称：
 3 *   GradSharp()
 4 * 参数:
 5 *   BYTE  bThre        - 阈值
 6 * 返回值:
 7 *   BOOL            - 成功返回TRUE，否则返回FALSE。
 8 * 说明:
 9 *   该函数用来对图像进行梯度锐化,设定梯度锐化的阈值为30
10 /************************************************************************/
11 void CImgEnhance::GradSharp(unsigned char Thre)
12 {
13        unsigned char*    pSrc;       // 指向源图像的指针
14     unsigned char*    pDst;     
15     unsigned char*    pSrc1;
16     unsigned char*    pSrc2;    
17     LONG    i,j;                // 循环变量
18     int    bTemp;
19     if(m_pImgDataOut != NULL)
20     {
21         delete []m_pImgDataOut;
22         m_pImgDataOut = NULL;
23     }
24     
25     if(m_nBitCount != 8)
26     {
27         AfxMessageBox("只能处理8位灰度图像!");
28         return ;
29     }
30        int lineByte = (m_imgWidth * m_nBitCount / 8 + 3) / 4 * 4;
31 
32     //创建要复制的图像区域
33     m_nBitCountOut = m_nBitCount;
34     int lineByteOut = lineByte;
35     if (!m_pImgDataOut)
36     {
37         m_pImgDataOut = new unsigned char[lineByteOut * m_imgHeight];
38     }
39     
40     int pixelByte = m_nBitCountOut / 8; // 此处实际就是1，因为只处理8位图像
41     for(i = 0; i < m_imgHeight; i++)
42     {
43         for(j = 0; j < m_imgWidth * pixelByte; j++)
44             *(m_pImgDataOut + i * lineByteOut + j) = *(m_pImgData + i * lineByteOut + j);
45     }
46     
47     for(i = 0; i < m_imgHeight; i++)        // 每行
48     {        
49         for(j = 0; j < m_imgWidth; j++)        // 每列
50         {
51                 //指向新DIB第i行第j列的像素的指针
52             pDst = m_pImgDataOut + lineByte * i + j;
53             
54             // 进行梯度运算
55             // 指向DIB第i行，第j个象素的指针
56             pSrc  = (unsigned char*)m_pImgData + lineByte * i + j;            
57             // 指向DIB第i+1行，第j个象素的指针
58             pSrc1 = (unsigned char*)m_pImgData + lineByte * (i+1) + j;            
59             // 指向DIB第i行，第j+1个象素的指针
60             pSrc2 = (unsigned char*)m_pImgData + lineByte * i + j + 1; 
61             bTemp = abs((*pSrc)-(*pSrc1)) + abs((*pSrc)-(*pSrc2));
62             
63             // 判断是否小于阈值
64             if ((bTemp+120) < 255)
65             {
66                 // 判断是否大于阈值，对于小于情况，灰度值不变。
67                 if (bTemp >= Thre)
68                 {
69                     *pSrc = bTemp + 120;
70                 }
71             }
72             else
73             {
74                 *pSrc = 255;
75             }
76             //生成新的DIB像素值
77             *pDst = *pSrc;
78         }
79     }
80     
81 }