sobel算子进行边缘跟踪的c++代码（不利于OpenCV与MATLAB）

以下是使用Sobel算子进行边缘检测的C++代码示例：

#include <iostream>
#include <cmath>
using namespace std;

// Sobel算子
int sobel_operator(int Gx[3][3], int Gy[3][3], int pixel[3][3]) {
    int sumx = 0, sumy = 0;
    for (int i = 0; i < 3; ++i) {
        for (int j = 0; j < 3; ++j) {
            sumx += Gx[i][j] * pixel[i][j];
            sumy += Gy[i][j] * pixel[i][j];
        }
    }
    return sqrt(sumx * sumx + sumy * sumy);
}

// 边缘检测
void edge_detection(int height, int width, unsigned char *input, unsigned char *output) {
    // Sobel算子
    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}};

    for (int y = 1; y < height - 1; ++y) {
        for (int x = 1; x < width - 1; ++x) {
            // 获取3x3像素矩阵
            int pixel[3][3];
            for (int i = -1; i <= 1; ++i) {
                for (int j = -1; j <= 1; ++j) {
                    pixel[i+1][j+1] = input[(y+i) * width + (x+j)];
                }
            }
            // Sobel算子边缘检测
            int edge = sobel_operator(Gx, Gy, pixel);
            // 截断边缘灰度值
            output[y * width + x] = edge > 255 ? 255 : edge < 0 ? 0 : edge;
        }
    }
}

int main() {
    // 读入图像
    int width = 512, height = 512;
    unsigned char *input = new unsigned char[width * height];
    unsigned char *output = new unsigned char[width * height];
    FILE *fp = fopen("lena.raw", "rb");
    fread(input, sizeof(unsigned char), width * height, fp);
    fclose(fp);

    // 边缘检测
    edge_detection(height, width, input, output);

    // 输出图像
    FILE *fout = fopen("lena_edge_detection.raw", "wb");
    fwrite(output, sizeof(unsigned char), width * height, fout);
    fclose(fout);

    delete[] input;
    delete[] output;

    return 0;
}

该代码使用Sobel算子进行边缘检测，并将结果输出到文件中。其中，`edge_detection`函数对每个像素进行边缘检测，`sobel_operator`函数计算Sobel算子的响应值。