c++实现矩阵直接做dct变换

#include <iostream>  
#include <string.h>
#include <opencv2/core/core.hpp>  
#include <opencv2/highgui/highgui.hpp>  
#include <opencv2\core\mat.hpp>
#include<opencv2/opencv.hpp>//如果要包含所有库
#define PI 3.1415926
using namespace cv;
using namespace std;
//将mat类型转换为二维数组的类型
//进行dct变换
double Coun(int x)
{
	if (x != 0) return 1.0;
	return 1.0 / (sqrt(2.0));
}
void mydct(double (*matrix)[8], int M, int N)
{
	double matrix1[4][8];

	int p, q;
	int m, n;
	double ap, aq;
	double num;
	double max;
	//double *point = (double *)A.data;
	int i, j;
	//for (i = 0; i < M; i++)
	//{
	//	for (j = 0; j < N; j++)
	//	{
	//		//cout << point[i*N + j] << "\n";
	//		cout << A.at<float>(i, j) << "\n";
	//	}
	//}


	//for (p = 0; p < M; p++)
	//{
	//	cout << p << "\n";
	//	for (q = 0; q < N; q++)
	//	{
	//		double tmp = 0.0;

	//		for (m = 0; m < M; m++)
	//		{
	//			for (n = 0; n < N; n++)

	//			{
	//				tmp += A.data[m*N + n] * cos((2 * m + 1)*p*PI / (2.0 * M)) * cos((2 * n + 1)*q*PI / (2.0 * N));//每个像素对应的tmp值即累加式计算完毕
	//			}
	//		}
	//		//将像素值赋给B.data
	//		B.at<float>(p,q)= 2.0*Coun(p)*Coun(q)*tmp/M;
	//		cout << "p:" << p << " q:&#