离散傅里叶变换---准找图片旋转的角度

 
	Mat i=imread("..//images//imageTextR.png",CV_LOAD_IMAGE_GRAYSCALE);
	if(i.empty())return -1;

	//1 转化到 傅里叶 最佳尺寸
	Mat padded;
	int m=getOptimalDFTSize(i.rows);//获取 DFT(离散傅里叶变换) 的最佳尺寸
	int n=getOptimalDFTSize(i.cols);
	copyMakeBorder(i,padded,0,m-i.rows,0,n-i.cols,BORDER_CONSTANT,Scalar::all(0));//(左右上下方向)填充边框


	//2 为傅立叶变换的结果(实部和虚部)分配存储空间
	Mat planes[]={Mat_<float>(padded),Mat::zeros(padded.size(),CV_32F)};//数组：1相位图像（原图--转化为浮点型）2 置零的浮点型图像
	Mat complexI;//复合图像
	merge(planes,2,complexI);//将2个单通道图像 合并成一个 2通道图像

	//3 进行离散傅立叶变换
	dft(complexI,complexI);

	//4 将复数转换为幅度
	split(complexI,planes);
	magnitude(planes[0],planes[1],planes[0]);//将复数转换为幅度
	Mat magI=planes[0];

	//5 对数尺度(logarithmic scale)缩放
	magI+=Scalar::all(1);//转换到对数尺度
	log(magI,magI);//计算每个元素的自然对数

	//6 剪切和重分布幅度图象限
	magI=magI(Rect(0,0,magI.cols-2,magI.rows-2));
	int cx=magI.cols/2;
	int cy=magI.rows/2;

	Mat q0(magI,Rect(0,0,cx,cy));//
	Mat q1(magI,Rect(cx,0,cx,cy));
	Mat q2(magI,Rect(0,cy,cx,cy));
	Mat q3(magI,Rect(cx,cy,cx,cy));

	Mat tmp;
	q0.copyTo(tmp);
	q3.copyTo(q0);
	tmp.copyTo(q3);

	q1.copyTo(tmp);
	q2.copyTo(q1);
	tmp.copyTo(q2);

	normalize(magI,magI,0,1,CV_MINMAX);

	imshow("src",magI);
	waitKey(0);