c语言数学形态边缘检测图像的代码实现

以下是使用C语言实现数学形态边缘检测的代码示例：

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

#define WIDTH 256
#define HEIGHT 256

int erosion(int img[WIDTH][HEIGHT], int se[3][3], int x, int y)
{
    int i, j;
    for (i = -1; i <= 1; i++) {
        for (j = -1; j <= 1; j++) {
            if (se[i + 1][j + 1] && img[x + i][y + j] == 0) {
                return 0;
            }
        }
    }
    return 255;
}

int dilation(int img[WIDTH][HEIGHT], int se[3][3], int x, int y)
{
    int i, j;
    for (i = -1; i <= 1; i++) {
        for (j = -1; j <= 1; j++) {
            if (se[i + 1][j + 1] && img[x + i][y + j] == 255) {
                return 255;
            }
        }
    }
    return 0;
}

void erosion_image(int img[WIDTH][HEIGHT], int se[3][3])
{
    int temp[WIDTH][HEIGHT];
    int i, j;
    for (i = 1; i < HEIGHT - 1; i++) {
        for (j = 1; j < WIDTH - 1; j++) {
            temp[i][j] = erosion(img, se, i, j);
        }
    }
    memcpy(img, temp, WIDTH * HEIGHT * sizeof(int));
}

void dilation_image(int img[WIDTH][HEIGHT], int se[3][3])
{
    int temp[WIDTH][HEIGHT];
    int i, j;
    for (i = 1; i < HEIGHT - 1; i++) {
        for (j = 1; j < WIDTH - 1; j++) {
            temp[i][j] = dilation(img, se, i, j);
        }
    }
    memcpy(img, temp, WIDTH * HEIGHT * sizeof(int));
}

void edge_detection(int img[WIDTH][HEIGHT])
{
    int se[3][3] = {
        {0, 1, 0},
        {1, 1, 1},
        {0, 1, 0}
    };
    dilation_image(img, se);
    erosion_image(img, se);
}

int main()
{
    int img[WIDTH][HEIGHT];
    int i, j;

    // initialize image with random values
    for (i = 0; i < HEIGHT; i++) {
        for (j = 0; j < WIDTH; j++) {
            img[i][j] = rand() % 256;
        }
    }

    // perform edge detection
    edge_detection(img);

    // print result
    for (i = 0; i < HEIGHT; i++) {
        for (j = 0; j < WIDTH; j++) {
            printf("%d ", img[i][j]);
        }
        printf("\n");
    }

    return 0;
}

在以上代码中，我们定义了一个`erosion()`函数和一个`dilation()`函数来执行形态学腐蚀和膨胀操作。这些函数采用一个3x3的结构元素和一个像素位置作为参数，并返回输出像素的值。我们还定义了`erosion_image()`和`dilation_image()`函数来在整个图像上执行形态学腐蚀和膨胀操作。最后，我们定义了一个`edge_detection()`函数来执行边缘检测操作，它首先对图像进行膨胀操作，然后进行腐蚀操作。在`main()`函数中，我们首先初始化一个随机图像，然后执行边缘检测，并将结果打印出来。