图像锐化——高斯—拉普拉斯算子

/**************************************************
 * 功能: 设定指定位置的像素灰度
 * 参数: imageBuf为目标图像 x,y为要设定像素的坐标
 **************************************************/
 void SetPixelXY(uchar** imageBuf1, int x, int y,int nc, int a)
 {
	 if (nc ==1 )
	 {
		 imageBuf1[y][x *nc] = a;
	 }
	 else if (nc ==3)
	 {
		 imageBuf1[y][x *nc] = a;
		 imageBuf1[y][x *nc + 1] = a;
		 imageBuf1[y][x *nc + 2] = a;
	 }
	 else if (nc == 4)
	 {
		 imageBuf1[y][x * nc] = a;
		 imageBuf1[y][x * nc + 1] = a;
		 imageBuf1[y][x * nc + 2] = a;
		 imageBuf1[y][x * nc + 3] = 255;
	 }
 }

 int GetAsh(uchar** imageBuf, int x, int y, int nc)
 {
	 int clr = 0;
	 if (nc == 1)
	 {
		 clr = imageBuf[y][x * nc] ;
	 }
	 else
	 {
		 clr = (imageBuf[y][x * nc] + imageBuf[y][x * nc + 1] + imageBuf[y][x * nc + 2]) / 3;
	 }

	 return clr;
 }

 /********************************************************
 * 把线形存储的像素转化为二维数组形式
 * 参数: image 线形存储的像素， width,height 图象的长宽
 ********************************************************/
 uchar** CreatImage(uchar * image, unsigned int width, unsigned int height, int nc)
 {
	 int imgWidthStep = ((width*nc + 3) / 4) * 4;

	 uchar** imageBuf = (uchar**)malloc(sizeof(uchar*)*(height));
	 for (int y = 0; y < height; y++)
	 {
		 //使imageBuf中每个指针分别指向其下标表示的行的行首地址
		 imageBuf[y] = image + y*imgWidthStep;
	 }
	 return imageBuf;
 }

 /**************************************************
 * 功能: 使用模板对彩色图邻域进行运算
 * 参数: imageBuf为目标图像 w、h为图像大小
 *       templt为模板 tw为邻域大小
 *		x,y为要取得像素的坐标
 *       cn为颜色分量编号 0为蓝色 1为绿色 2为红色
 **************************************************/
 int TempltExcuteCl(uchar** imageBuf0, int w, int h, int nc,int* templt, int tw, int x, int y, int cn)
 {
 	int i, j;                        //循环变量
 	int m = 0;                      //用来存放加权和
 	int px, py;
 	//依次对邻域中每个像素进行运算
 	for (i = 0; i < tw; i++)
 	{
 		for (j = 0; j < tw; j++)
 		{
 			//计算对应模板上位置的像素在原图像中的位置
 			py = y - tw / 2 + i;
 			px = x - tw / 2 + j;
 			//加权求和
 			m += imageBuf0[py][px * nc + cn] * templt[i*tw + j];
 		}
 	}
 	return m;                     //返回结果
 }

 /******************************************************************
 * 功能: 彩色图像的拉普拉斯锐化处理（scale = 3）
 * 参数: image0为原图形，image1锐化结果，
 *		w、h为图象的宽和高
 ******************************************************************/
 void SharpLaplacianCl(uchar* image0, uchar* image1, unsigned int w, unsigned int h, int nc, float factor = 1.0)
 {
 	//将图像转化为矩阵形式
 	uchar** imageBuf0 = CreatImage(image0, w, h, nc);
 	uchar** imageBuf1 = CreatImage(image1, w, h, nc);
 	//设定模板
 	//int templt[9] = { -1, -1, -1, -1, 8, -1, -1, -1, -1 };

	int templt[25] = { 0, 0, -1, 0, 0,
		0, -1, -2, -1, 0,
		-1, -2, 16,-2, -1,
		0, -1, -2, -1, 0,
		0, 0, -1, 0, 0};

 	int x, y, c;
 	int a;
 	//设定衰减因子
 
 	//依次对原图像的每个像素进行处理
 	for (y = 2; y < h - 2; y++)
 	{
 		for (x = 2; x < w/2 - 2; x++)
 		{
 			for (c = 0; c < nc; c++)
 			{
 				//利用拉普拉斯模板对邻域进行处理
 				a = TempltExcuteCl(imageBuf0, w, h,nc, templt, 5, x, y, c);
 				a = (int)((float)a / factor);
 				//对中心像素进行增强
 				a = imageBuf0[y][x * nc + c] + a;
 				//过限处理
 				a = a > 255 ? 255 : a;
 				a = a < 0 ? 0 : a;
 				imageBuf1[y][x * nc + c] = a;
 			}
 		}
 	}
 	//清理内存
 	free(imageBuf0);
 	free(imageBuf1);
 }