图像锐化主要影响图像中的低频分量,不影响图像中的高频分量。
图像锐化的主要目的有两个:
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 }