Operadores de Sabel para processamento digital de imagens (Filtro 3x3)

operadorSabel.cpp

Mat operadoresSabel(Mat img){
    Mat imgResultante = img.clone();

    int gx,gy;

    for(int i=0; i< img.rows; i++){
        for(int j=0; j < img.cols; j++){
            gx =    (img.at<uchar>(i+1,j-1) + 2*(img.at<uchar>(i+1,j)) + img.at<uchar>(i+1,j+1)) -
                    (img.at<uchar>(i-1, j-1) + 2*(img.at<uchar>(i-1, j)) + img.at<uchar>(i-1, j+1));

            gy =    (img.at<uchar>(i-1,j+1) + 2*(img.at<uchar>(i,j+1)) + img.at<uchar>(i+1,j+1)) -
                    (img.at<uchar>(i-1, j-1) + 2*(img.at<uchar>(i, j-1)) + img.at<uchar>(i+1, j-1));

            imgResultante.at<uchar>(i,j) = sqrt(pow(gx, 2) + pow(gy, 2));
        }
    }

    return imgResultante;
}