【转】Sobel 算子

最新推荐文章于 2023-10-14 09:26:58 发布

volkswageos

最新推荐文章于 2023-10-14 09:26:58 发布

阅读量910

点赞数

文章标签： null float

本文链接：https://blog.csdn.net/volkswageos/article/details/6475359

版权

Sobel 算子

计算机视觉 2009-05-17 18:54 阅读70 评论0

字号：大中小小

由于 Sobel 算子结合了 Gaussian 平滑和微分，因此其结果或多或少对噪声有一定的鲁棒性。采用以下算子分别计算一阶 x 方向和 y 方向的图像差分：

#include <math.h>
// Sobel 算子
// 1. pImageData   图像数据
// 2. nWidth       图像宽度
// 3. nHeight      图像高度
// 4. nWidthStep   图像行大小
BOOL Sobel(unsigned char *pImageData, int nWidth, int nHeight, int nWidthStep)
{
    int            i            = 0;
    int            j            = 0;
    int            nDx          = 0;
    int            nDy          = 0;
    int            nValue       = 0;
    unsigned char *pLine[3     = { NULL, NULL, NULL };
    for (j = 1; j < nHeight - 1; j++)
    {
        pLine[0  = pImageData + nWidthStep * (j - 1);
        pLine[1  = pImageData + nWidthStep * j;
        pLine[2  = pImageData + nWidthStep * (j + 1);
        for (i = 1; i < nWidth - 1; i++)
        {
            nDx =
                pLine[0][i+1     - pLine[0][i-1     +
                pLine[1][i+1 * 2 - pLine[1][i-1 * 2 +
                pLine[2][i+1     - pLine[2][i-1];
            nDy =
                pLine[2][i-1     - pLine[0][i-1     +
                pLine[2][i   * 2 - pLine[0][i   * 2 +
                pLine[2][i+1     - pLine[0][i+1];
            nValue = (int) sqrt((float) (nDx * nDx + nDy * nDy));
            if (nValue > 0xFF)
            {
                nValue = 0xFF;
            }
            pLine[0][i-1 = (unsigned char) nValue;
        }
    }
    return TRUE;
}

Sobel 边缘检测效果：