在图像处理中,通过当前位置的邻域像素计算新的像素值是很常见的操作。当邻域包含图像的上几行和下几行时,就需要同时扫描图像的若干行,这就是图像的邻域操作了。至于模板操作是实现空间滤波的基础,通常是使用一个模板(一个的矩形)滑过整幅图像产生新的像素。下面介绍通过使用OpenCV2实现Laplace算子锐化图像,来介绍OpenCV2中对邻域和模板的操作。
锐化处理主要的目的是突出灰度的过渡部分,通常由微分来定义和实现锐化算子的各种方法。Laplace算子是最贱的各向同性微分算子,常用的Laplace模板如下:
使用Laplace算子锐化图像时需要注意模板中心的系数,如果中心系数是负的,就需要将原图像减去经过Laplace算子处理后的图像,得到锐化后的结果。如果中心系数是正的,则相反。
锐化图像时不能以In-Place的方式来完成,需要提供一个输出图像。在对图像遍历时使用两个3个指针:一个指向当前行,一个指向当前行的上一行,一个指向当前行的下一行。而且,由于每个像素值的计算都需要它的上下左右四个相邻像素,所以无法对图像的边界进行计算,需要另作处理。具体实现代码如下:
1: /*
2: 0 -1 0
3: -1 4 -1
4: 0 -1 0
5: */
6: void sharpen(const Mat & image,Mat & result)
7: {
8: CV_Assert(image.depth() == CV_8U);
9:
10: result.create(image.size(),image.type());
11:
12: const int channels = image.channels() ;
13: for(int j = 1 ;j < image.rows - 1 ; j ++){
14: const uchar * previous = image.ptr<const uchar>(j - 1) ; // 当前行的上一行
15: const uchar * current = image.ptr<const uchar>(j) ; //当前行
16: const uchar * next = image.ptr<const uchar>(j + 1) ; //当前行的下一行
17:
18: uchar * output = result.ptr<uchar>(j) ; // 输出行
19: for(int i = channels ; i < channels * (image.cols - 1) ; i ++) {
20: * output ++ = saturate_cast<uchar>( 4 * current[i] - previous[i] - next[i] - current[i - channels] - current[ i + channels]) ;
21: }
22: }
23:
24: //对图像边界进行处理
25: //边界像素设置为0
26: result.row(0).setTo(Scalar(0));
27: result.row(result.rows-1).setTo(Scalar(0)) ;
28: result.col(0).setTo(Scalar(0));
29: result.col(result.cols-1).setTo(Scalar(0));
30: }
这里使用指针遍历整个图像,使用三个指针同时扫描图像的三行,另外使用一个指针指向输出行。在计算输出像素时,使用模板函数saturate_cast<uchar>对计算结果进行调整。这是因为对像素值的计算有可能导致结果超出了像素允许的范围,即小于0或者大于255,当计算结果是浮点数时,该函数会将结果取整至最近的整数。
由于边界的像素没有完整的邻域,无法使用模板计算其值,需要单独处理。这里只是简单的将其值设为0.
测试代码:
1: Mat image = imread("d:\\lenna.jpg") ;
2: Mat result ;
3: sharpen(image,result) ;
4: imwrite("d:\\lenna1.jpg",result) ;
5: imwrite("d:\\lenna2.jpg",image + result) ;
结果如下:
这里需要说明,在OpenCV2中对Mat进行了大量的运算符重载,例如上面,两幅图像相加直接使用image + result即可。另外,如位操作符:&,|,^,~;函数max,min,abs;比较操作符:<,<=,>,>=,比较操作符返回一个8位二进制图像。另外矩阵乘法m1 * m2,矩阵求逆 m1.inv(),矩阵转置m.t(),举证的行列式m.determinate(),向量的模v.norm(),向量叉乘v.corss(v1),向量的点乘v.dot(v1)等。
由于使用模板实现空间滤波在图像处理中非常的常用,在OpenCV中专门定义了一个特殊的函数来完成该处理:
filter2D.
函数原型:
1: void filter2D(InputArray src, OutputArray dst, int ddepth, InputArray kernel, Point anchor=Point(-1,-1), double delta=0, int borderType=BORDER_DEFAULT )
参数:
src,输入图像
dst 输出图像,和输入图像有相同的大小和通道
ddepth
输出图像的depth,如果为负数,则和输入图像的depth相同
kernel 模板
anchor
进行卷积运算的中心位置,默认的是kernel的中心
delta 可选值,加到输出像素上的值
bordertype
对输出图像边界的处理。
使用filter2D实现Laplace算子
1: Mat kern = (Mat_<char> (3,3) << 1,1,1,
2: 1,-8,1,
3: 1,1,1) ;
4: filter2D(image,result,image.depth(),kern) ;
这里只需要定义好kernel调用filter2D即可,而且OpenCV对该函数进行了优化,其效率要比上面使用指针实现的要高。
原文:http://www.cnblogs.com/wangguchangqing/p/3710397.html