canny算子四部曲之二（求梯度）

梯度求法和sobel之类的算子雷同，甚至更简单，就是一个离散差分，不清楚的童鞋可以百度，一大堆资料呢，从源码也可清晰的看出原理。

// 方向导数,求梯度
/*
 *  @parameter sz:     图像大小
 *  @parameter pGray:    图像灰度值
 *  @parameter pGradx:  图像x轴梯度
 *  @parameter pGrady:  图像y轴梯度
 *  @parameter pMag   :  图像梯度幅值
 */

void Grad(SIZE sz, LPBYTE pGray,int *pGradX, int *pGradY, int *pMag)
{
    LONG y,x;

    //中间变量
    double dSqt1;
    double dSqt2;

    //x方向的方向导数
    for(y=1;y<sz.cy-1;y++)
    {
        for(x=1;x<sz.cx-1;x++)
        {
            pGradX[y*sz.cx +x] = (int)( pGray[y*sz.cx+x+1]-pGray[y*sz.cx+ x-1]  );
        }
    }

    //y方向方向导数
    for(x=1;x<sz.cx-1;x++)
    {
        for(y=1;y<sz.cy-1;y++)
        {
            pGradY[y*sz.cx +x] = (int)(pGray[(y+1)*sz.cx +x] - pGray[(y-1)*sz.cx +x]);
        }
    }

    //求梯度
    for(y=0; y<sz.cy; y++)
    {
        for(x=0; x<sz.cx; x++)
        {
            //二阶范数求梯度
            dSqt1 = pGradX[y*sz.cx + x]*pGradX[y*sz.cx + x];
            dSqt2 = pGradY[y*sz.cx + x]*pGradY[y*sz.cx + x];
            pMag[y*sz.cx+x] = (int)(sqrt(dSqt1+dSqt2)+0.5);//四舍五入
        }
    }
}

用到的类型：

typedef struct {
    int cy;
    int cx;
} SIZE;

typedef unsigned char *LPBYTE;

typedef long LONG;