OpenCV学习4：图像锐化

前言

开局一张图，内容全靠编。

简介

图像锐化(image sharpening)是补偿图像的轮廓，增强图像的边缘及灰度跳变的部分，使图像变得清晰。

函数声明

void sharpen(const cv::Mat &image, cv::Mat &result);

函数定义

void sharpen(const cv::Mat &image, cv::Mat &result) {
	// 判断是否需要分配图像数据。如果需要，就分配
	result.create(image.size(), image.type());
	int nchannels = image.channels(); // 获得通道数
									  // 处理所有行（除了第一行和最后一行）
	for (int j = 1; j < image.rows - 1; j++) {
		const uchar* previous = image.ptr<const uchar>(j - 1); // 上一行
		const uchar* current = image.ptr<const uchar>(j); // 当前行
		const uchar* next = image.ptr<const uchar>(j + 1); // 下一行
		uchar* output = result.ptr<uchar>(j); // 输出行
		for (int i = nchannels; i < (image.cols - 1)*nchannels; i++) {
			// 应用锐化算子
			*output++ = cv::saturate_cast<uchar>(
				5 * current[i] - current[i - nchannels] -
				current[i + nchannels] - previous[i] - next[i]);
		}
	}
	// 把未处理的像素设为 0 
	result.row(0).setTo(cv::Scalar(0,0,0));
	result.row(result.rows - 1).setTo(cv::Scalar(0,0,0));
	result.col(0).setTo(cv::Scalar(0,0,0));
	result.col(result.cols - 1).setTo(cv::Scalar(0,0,0));
}

图像扫描中使用了三个指针，一个表示当前行、一个表示上面的行、一个表示下面的行。另外，因为在计算每一个像素时都需要访问与它相邻的像素，所以有些像素的值是无法计算的，比如第一行、最后一行和第一列、最后一列的像素。

实现原理

在计算输出像素的值时，我们调用了 cv::saturate_cast 模板函数，并传入运算结果。这是因为计算像素的数学表达式的结果经常超出允许的范围（即小于 0 或大于 255）。使用这个函数可把结果调整到 8 位无符号数的范围内，具体做法是把小于 0 的数值调整为 0，大于 255 的数值调整为 255——这就是 cv::saturate_cast函数的作用。此外，如果输入参数是浮点数，就会得到最接近的整数。可以在调用这个函数时显式地指定其他数据类型，以确保结果在该数据类型定义的范围之内。
由于边框上的像素没有完整的相邻像素，因此不能用前面的方法计算，需要另行处理。这里简单地把它们设置为 0。

测试代码：

#include <opencv2/core.hpp> 
#include <opencv2/highgui.hpp>
#include <opencv2/opencv.hpp>
#include <iostream>

void sharpen(const cv::Mat &image, cv::Mat &result);

int main(int argc, char* argv[])
{
	cv::Mat img;
	cv::Mat result;
	img=cv::imread("test.bmp");
	if(img.empty())
	{
		std::cout << "read error"<<std::endl;
	}
	else
	{
		
		sharpen(img, result);
		//显示结果
		cv::imshow("img", img);
		cv::imshow("result", result);
		cv::waitKey(0);
	}	
	system("pause");
	return 0;
}

void sharpen(const cv::Mat &image, cv::Mat &result) {
	// 判断是否需要分配图像数据。如果需要，就分配
	result.create(image.size(), image.type());
	int nchannels = image.channels(); // 获得通道数
									  // 处理所有行（除了第一行和最后一行）
	for (int j = 1; j < image.rows - 1; j++) {
		const uchar* previous = image.ptr<const uchar>(j - 1); // 上一行
		const uchar* current = image.ptr<const uchar>(j); // 当前行
		const uchar* next = image.ptr<const uchar>(j + 1); // 下一行
		uchar* output = result.ptr<uchar>(j); // 输出行
		for (int i = nchannels; i < (image.cols - 1)*nchannels; i++) {
			// 应用锐化算子
			*output++ = cv::saturate_cast<uchar>(
				5 * current[i] - current[i - nchannels] -
				current[i + nchannels] - previous[i] - next[i]);
		}
	}
	// 把未处理的像素设为 0 
	result.row(0).setTo(cv::Scalar(0,0,0));
	result.row(result.rows - 1).setTo(cv::Scalar(0,0,0));
	result.col(0).setTo(cv::Scalar(0,0,0));
	result.col(result.cols - 1).setTo(cv::Scalar(0,0,0));
}

测试结果：

原图：

处理后图片：