sobel算子锐化C++

Sobel算子是一种常用的边缘检测算子，可以用于图像锐化。下面是C语言实现Sobel算子锐化的示例代码：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#define WIDTH 512
#define HEIGHT 512

int main()
{
    int i, j, k, l;
    int sobel_x, sobel_y, sobel;
    int image[HEIGHT][WIDTH];
    int sharp_image[HEIGHT][WIDTH];
    int gx[3][3] = {{-1, 0, 1}, {-2, 0, 2}, {-1, 0, 1}};
    int gy[3][3] = {{-1, -2, -1}, {0, 0, 0}, {1, 2, 1}};

    // 读取原始图像
    FILE *fp = fopen("lena512.raw", "rb");
    fread(image, sizeof(unsigned char), WIDTH * HEIGHT, fp);
    fclose(fp);

    // Sobel算子锐化
    for (i = 1; i < HEIGHT - 1; i++) {
        for (j = 1; j < WIDTH - 1; j++) {
            sobel_x = sobel_y = 0;
            for (k = -1; k <= 1; k++) {
                for (l = -1; l <= 1; l++) {
                    sobel_x += (image[i + k][j + l] * gx[k + 1][l + 1]);
                    sobel_y += (image[i + k][j + l] * gy[k + 1][l + 1]);
                }
            }
            sobel = (int) sqrt(sobel_x * sobel_x + sobel_y * sobel_y);
            if (sobel > 255) sobel = 255;
            sharp_image[i][j] = sobel;
        }
    }

    // 将结果写入文件
    fp = fopen("lena512_sharp.raw", "wb");
    fwrite(sharp_image, sizeof(unsigned char), WIDTH * HEIGHT, fp);
    fclose(fp);

    return 0;
}

该示例代码使用Sobel算子对lena512.raw图像进行锐化，并将结果写入lena512_sharp.raw文件中。在该示例代码中，使用了两个3x3的矩阵gx和gy来表示Sobel算子的两个卷积核。在Sobel算子锐化过程中，对于每个像素点，先将其周围的像素点与卷积核进行卷积，得到横向和纵向的梯度值，然后计算其模长，即得到锐化后的像素值。