c - 实现 Sobel 过滤器时出现奇怪的结果-6ren

c - 实现 Sobel 过滤器时出现奇怪的结果

转载作者：太空宇宙更新时间：2023-11-04 02:56:23

27

4

<分区>

我一直在学习计算机视觉，想用 C 语言实现一些简单的技术。对于第一种技术，我正在做 Sobel 边缘检测滤波器。我了解它的工作原理，因此我认为编码应该相当容易，但我得到的结果却很奇怪。

我正在使用下面的图片:

Box

并得到这个结果

Filtered box

新结果!

Filtered box...a little more correct

请注意，我使用的是.ppm图片格式(链接是jpgs，因为我找不到支持.ppm的图片主机)

无论如何，这是我实现 Sobel 的代码部分:

/**********************************************************
This program takes in an image file and applies the Sobel
Filter edge detection technique to it.
**********************************************************/

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

void sobelFilter(){


    //Sobel kernels dx (horizontal) and dy (vertical)
    int horizFilter[3][3] = {{ 1,   0,  -1},
                             { 2,   0,  -2},
                             { 1,   0,  -1}};
    int vertFilter[3][3]  = {{ 1,   2,   1},
                             { 0,   0,   0},
                             {-1,  -2,  -1}};
    int pixVal = 0; 
    int horizPixVal = 0;
    int vertPixVal = 0;
    int x, y, i, j;

    //Quick check to make sure dimensions are correct
    printf("Using a Width of: %d\n", width);
    printf("Using a Height of: %d\n\n", height);

    //Start filtering process here
    for(x = 0; x < width; x++){
        for(y = 0; y < height; y++){
            pixVal = 0;
            horizPixVal = 0;
            vertPixVal = 0;
            if(!((x == 0) || (x == width-1) || (y == 0) || (y == height-1))){           //If the current pixel is along the border, ignore it and set to zero
                for(i = -1; i <= 1; i++){                                               //because the kernel does not align to it
                    for(j = -1; j <= 1; j++){
                        horizPixVal += (int)(image[y + j][x + i][0]) * horizFilter[i + 1][j + 1];       //Only need to focus on one of the RGB values since the output is
                        vertPixVal  += (int)(image[y + j][x + i][0]) * vertFilter[i + 1][j + 1];        //greyscale and all three values are the same
                    }
                }
            }
            pixVal = sqrt((horizPixVal * horizPixVal) + (vertPixVal * vertPixVal));     //Calculate magnitude
            pixVal = sqrt(horizPixVal * horizPixVal);
            if(pixVal > 255) pixVal = 255;                                              //Clamp value within 8-bit range
            filteredImage[y][x][0] = (unsigned char)pixVal;                             
        }
    }
}

这是读取 .ppm 文件的代码:

unsigned char image[MAX_IMAGE_HEIGHT][MAX_IMAGE_WIDTH][3];
unsigned char filteredImage[MAX_IMAGE_HEIGHT][MAX_IMAGE_WIDTH][3];

void readPPMImageData(){
    char fileName[MAX_NAME];
    char imageBuff[MAX_BUFF];
    width = 0;
    height = 0;
    maxColor = 0;
    int x;
    int y;

    FILE* file;

    printf("------------------------------------------------------------\n");
    printf("Now attempting to read in the .ppm image file data...\n");
    printf("------------------------------------------------------------\n\n");
    printf("What is the image file name (*.ppm)?  : ");
    scanf("%s", fileName);
    file = fopen(fileName, "rb");               //open the file specified by the user in binary read mode
    if(file == NULL){                           //but if the file was not found, terminate program
        printf("\nThe file %s could not be found! Terminating program...\n", fileName);
        exit(1);
    }

    //The first step is to read in the file type and check it agains P6 (file type of .ppm images)
    fgets(imageBuff, MAX_BUFF, file);
    if(imageBuff[0] != 'P' || imageBuff[1] != '6'){
        printf("\nInvalid image type! Acceptable type is: %s --- Received type is: %c%c\n\n", "P6", imageBuff[0], imageBuff[1]);
    }
    printf("Magic Number is: %c%c\n", imageBuff[0], imageBuff[1]);

    while(width == 0 || height == 0){
        fgets(imageBuff, MAX_BUFF, file);
        if(imageBuff[0] != '#') {
            sscanf(imageBuff, "%d %d", &width, &height);
        }
    }
    printf("Width is: %d\n", width);
    printf("Height is: %d\n", height);
    //if(feof(file)){
    //
    //}

    while(maxColor == 0){
        fgets(imageBuff, MAX_BUFF, file);
        if(imageBuff[0] != '#') {
            sscanf(imageBuff, "%d", &maxColor);
        }
    }
    printf("Maximum color value is: %d\n", maxColor);

    for(x = 0; x < width; x++){
        for(y = 0; y < height; y++){
            image[y][x][0] = (unsigned char)fgetc(file); //Get Red value
            image[y][x][1] = (unsigned char)fgetc(file); //Get Green value
            image[y][x][2] = (unsigned char)fgetc(file); //Get Blue value
        }
    }
    printf("Finished reading image data!\n\n");

    fclose(file);
}

下面是过滤后创建新 .ppm 文件的代码:

void createPPMImage(){
    char fileName[MAX_NAME]; 
    FILE* file;
    int x;
    int y;

    printf("------------------------------------------------------------\n");
    printf("Now attempting to create new .ppm image file...\n");
    printf("------------------------------------------------------------\n\n");
    printf("What is the name of the output image file (*.ppm)?  : ");
    scanf("%s", fileName);

    printf("Width is: %d\n", width);
    printf("Height is: %d\n", height);
    printf("Maximum color value is: %d\n", maxColor);

    file = fopen(fileName, "wb");
    fputs("P6\n", file);
    fprintf(file, "%d %d\n", width, height);
    fprintf(file, "%d\n", maxColor);

    for(x = 0; x < width; x++){
        for(y = 0; y < height; y++){
            fputc(filteredImage[y][x][0], file); //Write Red value
            fputc(filteredImage[y][x][0], file); //Write Green value
            fputc(filteredImage[y][x][0], file); //Write Blue value
        }
    }
    printf("Finished creating new filtered image!\n\n");

    fclose(file);
}

我 100% 确定问题不在于图像的读取或写入，因为我在没有应用过滤器的情况下测试了这些函数，并且只有在我使用上述函数后才会出现问题。

感谢任何帮助，因为据我所知，索引/公式似乎已正确实现，但这显然不是真的。

编辑:正如 Dave 和其他人所指出的，我不再 100% 确定错误在 Sobel 函数内，看起来这只是我在使用.ppm 格式。我继续发布了我的 .ppm 读取器/写入器函数的代码以及应用下面 anatolyg 提出的 [y][x][color] 方案后我得到的新结果。如果我的帖子太长，我很抱歉，如果是的话，请告诉我，因为这是我的第一篇帖子，我还不完全确定什么是合适的。

27

4

0

文章推荐： javascript - 将 python 请求转移到 node.js 请求时遇到问题

文章推荐： jquery - 网站上的图像高度错误

文章推荐： javascript - 等待查询结果

文章推荐： javascript - 根据内容的长度使 div 成为浏览器高度的 100%

java - Sobel 边缘检测会产生噪声
我尝试用java实现Sobel边缘检测。它有点有效，但我收到了很多看似随机的噪音...... 我将图像加载为 BufferedImages 并首先将其转换为灰度图像(通过我在网上找到的算法)。之后我计
java - Sobel 运算符不适用于矩形图像
我尝试在 Java 中实现 Sobel 运算符，但结果只是一些像素的混合。 int i, j; FileInputStream inFile = new FileInputStream
c++ - Sobel 的梯度是什么意思？
我有每个像素的 Sobel 算子的梯度。在我的例子中是 320x480。但是我怎样才能将它们与方向联系起来呢？例如，我打算绘制指纹方向图。那么，我该如何开始呢？是将梯度分成 block (例如 16
java - Sobel 过滤器无法正常工作
我为 Sobel 运算符编写了一个用于边缘检测的类，但是当我使用示例图像时，我的边缘消失了。如果有人可以帮助我，我将不胜感激。 import java.awt.image.BufferedImage;
image-processing - Sobel 过滤器是否意味着要归一化？
x-导数 Sobel 看起来是这样的: -1 0 +1 -2 0 +2 -1 0 +1 假设我的图像有两个样本看起来像这样(0=黑色，1=白色): 0 0 1 1 0 0 0 0
java - Sobel 边缘检测、图像卷积
我正在使用图像上的索贝尔边缘检测进行作业。我目前正在努力进行渐变操作。编译时收到“二元运算符 '*' 的操作数类型错误”错误。我认为这可能是因为我将所有像素定义为字母，并且我不确定下一步应该是什么。任
java - 我的 Sobel 边缘检测运算符输出很奇怪
我的 Sobel 边缘检测算子的输出很奇怪。这是我的代码: BufferedImage temp = img; float kernelx[][] = {{-1, 0, 1},{-2,
android - 如何使用一个像素的 Sobel 矩阵
我想在我的 android 应用程序中使用 Sobel 运算符。但是我不明白如何使用一个像素。 int sobel_x[][] = {{-1, 0, 1}, {-2, 0, 2},
c++ - 使用高斯模糊和 Sobel 算子保留颜色
我想在图像上使用高斯模糊作为使用 sobel 边缘检测过滤器的预处理步骤。我以前在灰度图像上有效地实现了 sobel 和高斯模糊运算符，但是，我从未尝试过在彩色图像上使用它们。之前，我一直采用像素
c++ - 没有缓冲区的 Sobel 边缘检测
对于嵌入式设计，我试图在不使用缓冲器的情况下在板上实现索贝尔边缘检测。即我直接从屏幕上阅读和写作。但是，我可以存储大约一两个图像宽度的数据以供以后引用。这是由于董事会规定的限制。但是我陷入了一些问题。
c - 实现 Sobel 过滤器时出现奇怪的结果
这个问题不太可能帮助任何 future 的访问者；它只与一个小的地理区域、一个特定的时间点或一个非常狭窄的情况有关，这些情况并不普遍适用于互联网的全局受众。为了帮助使这个问题更广泛地适用，visit
python - cv2.Sobel 中的比例是什么意思？
我试图理解 cv2.Sobel 中的 scale 参数.将 scale 设置为 1/8，我沿 x 轴得到如下输出: 但是对于 scale = 10 或 scale = 100，输出非常相似。上面两幅
c# - OpenCV:Sobel 滤波器不是旋转不变的
我正在使用 sobel 运算符通过 EmguCV 3.0(OpenCV 的 .NET 包装器)检测灰度图像中的边缘。我使用的代码是 Image gray = new Image(@"C:\gray.
python - 过滤掉微弱的 sobel 过滤线？
我正在使用 skimage 创建一个与此类似的 sobel 过滤器图像... 我想知道有没有办法锐化这个 sobel 滤镜图像？比如说，把那些比较淡的白线去掉，比如气球后面的淡淡的线条？我用 Ski
opencv - 什么是 Sobel 算子的命令？
OpenCV 文档说，(order == 0) 表示不会在这个方向上应用导数，即不会执行此内核的计算。 (Order == 1) 意味着这个方向的图像和内核只是一个简单的卷积。但是 (order =
opencv - 梯度角的 Sobel 算子
我想在文本中找到笔划的方向。 Sobel 算子如何用于此目的？这张图显示的是dp，也就是梯度方向。我想知道如何应用 Sobel 运算符找到要选择的像素，从 p 到 q，沿着路径 sp，到找到边缘上的
python - 在图像上应用 sobel 过滤器
我正在尝试对墙壁图像使用 sobel 过滤器，但它不起作用。我的代码是: im=scipy.misc.imread('IMG_1479bis.JPG') im = im.astype('int32'
algorithm - 3x3 Sobel 算子和梯度特征
阅读一篇论文，我很难理解所描述的算法: 给定手写样本的黑白数字图像，裁剪出单个字符进行分析。由于这可以是任意大小，因此算法需要考虑到这一点(如果更简单，我们可以假设大小为 2^n x 2^m)。现在
c++ - 实现 Sobel 算子
我正在尝试在水平和垂直方向上实现 sobel 运算符。但不知何故我得到了反向输出。我在下面附上的代码。对于水平蒙版 char mask [3][3]= {{-1,-2,-1},{0,0,0},{1,
c++ - sobel 算子 - 假边自然图像
我在使用 Sobel 算子进行边缘检测时遇到问题:它会产生太多假边缘，效果如下图所示。我正在使用 3x3 sobel 运算符 - 首先提取垂直然后水平，最终输出是每个滤波器输出的幅度。合成图像的边

首页

博学

6Ren·AI

商城

c - 实现 Sobel 过滤器时出现奇怪的结果