Ошибка Opencv

#include <opencv2/highgui/highgui.hpp>
#include <opencv2/highgui/highgui_c.h>
#include <math.h>
#include <iostream>
using namespace cv;
 
void NewHouph(const Mat& in, Mat& out)
{
    Mat Gx = Mat::zeros(in.size(), CV_64FC1);
    Mat Gy = Mat::zeros(in.size(), CV_64FC1);
    Mat Kb = Mat::zeros(in.size(), CV_64FC1);
    Mat Xy = Mat::zeros(in.size(), CV_64FC1);
    out = Mat::zeros(in.size(), CV_64FC1);
 
 
    double dx, dy;
    for (int i = 1; i < in.rows - 1; ++i)
    {
        for (int j = 1; j < in.cols - 1; ++j)
        {
            dx = ((double)in.at<unsigned char>(i, j + 1) - in.at<unsigned char>(i, j - 1)) / 2;
            dy = ((double)in.at<unsigned char>(i + 1, j) - in.at<unsigned char>(i - 1, j)) / 2;
            Gx.at<double>(i, j) = dx;
            Gy.at<double>(i, j) = dy;
            Xy.at<double>(i, j) = sqrt(dx*dx + dy*dy);
        }
    }   
 
    int count = 0;
    int sum;
    for (int k = 1; k < in.rows - 1; ++k)
    {
        for (int b = 1; b < in.cols - 1; ++b)
        {
            count = 0;
            sum = 0;
            for (int x = 1; x < in.rows - 1; ++x)
            {
                for (int y = 1; y < in.cols - 1; ++y)
                {
                    if (y == k*x + b)
                    {
                        sum+= Xy.at<double>(x, y);
                        count++;
                    }                   
                }
            }
            if (count!=0) Kb.at<double>(k, b) = sum / count;        
        }
    }
 
    for (int k = 1; k < in.rows - 1; ++k)
    {
        for (int b = 1; b < in.cols - 1; ++b)
        {
            if ((Kb.at<double>(k, b) > Kb.at<double>(k, b-1))
                &(Kb.at<double>(k, b) > Kb.at<double>(k, b+1))
                &(Kb.at<double>(k, b) > Kb.at<double>(k-1, b-1))
                &(Kb.at<double>(k, b) > Kb.at<double>(k-1, b))
                &(Kb.at<double>(k, b) > Kb.at<double>(k-1, b+1))
                &(Kb.at<double>(k, b) > Kb.at<double>(k+1, b-1))
                &(Kb.at<double>(k, b) > Kb.at<double>(k+1, b))
                &(Kb.at<double>(k, b) > Kb.at<double>(k+1, b+1)))
            {
                out.at<unsigned char>(k, b) = 255;
            }
            else
            {
                out.at<unsigned char>(k, b) = 0;
            }
        }
    }
 
}
 
 
int main(int argc, char* argv[])
{
    Mat orig, src;
 
    orig = imread("gray1.png");
 
    
    NewHouph(orig, src);
    imwrite("gray1_result.png", src);
 
    waitKey(0);
    return 0;
}