C++でMatの型が合わない問題

前提・実現したいこと

画像を読み込んで離散コサイン変換した後、逆離散コサイン変換して出力するプログラムを作成しています。
警告が沢山出て、プログラムが止まってしまいます。
多分、Matの型が合っていないと思うのですが、理解が足りずどうすればいいのか分かりません。
初歩的な質問で申し訳ありませんが、よろしくお願いします。
以下のコードを参考にしました。
https://github.com/yoyoyo-yo/Gasyori100knock/blob/master/Question_31_40/answers_cpp/answer_36.cpp

発生している問題・エラーメッセージ

重大度レベル	コード	説明	プロジェクト	ファイル	行	抑制状態
警告	C26451	Arithmetic overflow: Using operator '*' on a 4 byte value and then casting the result to a 8 byte value. Cast the value to the wider type before calling operator '*' to avoid overflow (io.2).	qim_1	c:\users\ak\downloads\opencv\build\install\include\opencv2\core\cuda.inl.hpp	301	
警告	C26451	Arithmetic overflow: Using operator '+' on a 4 byte value and then casting the result to a 8 byte value. Cast the value to the wider type before calling operator '+' to avoid overflow (io.2).	qim_1	c:\users\ak\downloads\opencv\build\install\include\opencv2\core\cuda.inl.hpp	301	
警告	C26451	Arithmetic overflow: Using operator '*' on a 4 byte value and then casting the result to a 8 byte value. Cast the value to the wider type before calling operator '*' to avoid overflow (io.2).	qim_1	c:\users\ak\downloads\opencv\build\install\include\opencv2\core\cuda.inl.hpp	482	
警告	C26451	Arithmetic overflow: Using operator '+' on a 4 byte value and then casting the result to a 8 byte value. Cast the value to the wider type before calling operator '+' to avoid overflow (io.2).	qim_1	c:\users\ak\downloads\opencv\build\install\include\opencv2\core\cuda.inl.hpp	482	
…

該当のソースコード

#include <opencv2/core.hpp>
#include <opencv2\opencv.hpp>
#include <opencv2/highgui.hpp>
#define _USE_MATH_DEFINES
#include <iostream>
#include <math.h>
#include <complex>
using namespace cv;
using namespace std;

const int height = 512, width = 512;

// DCT hyper-parameter
int T = 8;
int K = 8;

// Discrete Cosine transformation
void dct(Mat img, Mat dct_coef) {
	double I;
	double F;
	double Cu, Cv;

	for (int ys = 0; ys < height; ys += T) {
		for (int xs = 0; xs < width; xs += T) {
			for (int v = 0; v < T; v++) {
				for (int u = 0; u < T; u++) {
					F = 0;

					if (u == 0) {
						Cu = 1. / sqrt(2);
					}
					else {
						Cu = 1;
					}

					if (v == 0) {
						Cv = 1. / sqrt(2);
					}
					else {
						Cv = 1;
					}

					for (int y = 0; y < T; y++) {
						for (int x = 0; x < T; x++) {
							I = (double)img.at<uchar>(ys + y, xs + x);
							F += 2. / T * Cu * Cv * I * cos((2. * x + 1) * u * M_PI / 2. / T) * cos((2. * y + 1) * v * M_PI / 2. / T);
						}
					}

					dct_coef.at<uchar>(ys + v, xs + u) = F;
				}
			}
		}
	}
}

// Inverse Discrete Cosine transformation
void idct(Mat out, Mat dct_coef) {
	double f;
	double Cu, Cv;

	for (int ys = 0; ys < height; ys += T) {
		for (int xs = 0; xs < width; xs += T) {
				for (int y = 0; y < T; y++) {
					for (int x = 0; x < T; x++) {
						f = 0;

						for (int v = 0; v < K; v++) {
							for (int u = 0; u < K; u++) {
								if (u == 0) {
									Cu = 1. / sqrt(2);
								}
								else {
									Cu = 1;
								}

								if (v == 0) {
									Cv = 1. / sqrt(2);
								}
								else {
									Cv = 1;
								}

								f += 2. / T * Cu * Cv * dct_coef.at<uchar>(ys + v, xs + u) * cos((2. * x + 1) * u * M_PI / 2. / T) * cos((2. * y + 1) * v * M_PI / 2. / T);
							}
						}

						f = fmin(fmax(f, 0), 255);
						out.at<uchar>(ys + y, xs + x) = (uchar)f;
					}
				}
			}
		}

}



// Main
int main(int argc, const char* argv[]) {
	// read original image
	Mat img = cv::imread("lenna.jpg", cv::IMREAD_COLOR);
	// グレースケール変換
	Mat gray1;
	cvtColor(img, gray1, COLOR_RGB2GRAY);
	//std::cout << "image.channels(): " << gray1.channels() << std::endl;

	// output image
	Mat out = cv::Mat::zeros(img.rows, img.cols, CV_8U);
	imshow("gray", out);
	// DCT
	Mat dct_coef;
	dct(gray1, dct_coef);

	// IDCT
	idct(out, dct_coef);

	imwrite("out.jpg", out);

	waitKey();

}