opencv3.2.0图像离散傅里叶变换

源码：

##名称：离散傅里叶变换
##平台：QT5.7.1+opencv3.2.0
##日期：2017年12月13、

/**** 新建QT控制台程序****/

#include <QCoreApplication>
#include <opencv2/core/core.hpp>
#include <opencv2/highgui/highgui.hpp>
#include <opencv2/imgproc/imgproc.hpp>
#include <iostream>
using namespace std;
using namespace cv;

int main()
{
    //读入原始图像
   Mat srcImage = imread("/home/ttwang/Valley_logo.jpg",);
   imshow("srcImage",srcImage);

    /****将图像扩大到合适尺寸********/
   int m = getOptimalDFTSize(srcImage.rows);
   int n = getOptimalDFTSize(srcImage.cols);
   Mat padded;
   copyMakeBorder(srcImage,padded,,m-srcImage.rows,,n-srcImage.cols,BORDER_CONSTANT,
                  Scalar::all());

   /******为傅里叶变换的结果（实部和虚部)分配存储空间******/
   Mat planes[] = {Mat_<float>(padded),Mat::zeros(padded.size(),CV_32F)};
   Mat complexI;
   merge(planes,,complexI);

   /********进行傅里叶变换*********/
   dft(complexI,complexI);

   /*********将复数转化为幅值***********/
   split(complexI,planes);
   magnitude(planes[],planes[],planes[]);
   Mat magnitudeImage = planes[];

   /**********进行对数尺度**********/
   magnitudeImage += Scalar::all();
   log(magnitudeImage,magnitudeImage);

   /*******剪切和重分布幅度图象限 ******/
   magnitudeImage = magnitudeImage(Rect(,,magnitudeImage.cols & -,magnitudeImage.rows & -));
   int cx = magnitudeImage.cols / ;
   int cy = magnitudeImage.rows /;
   Mat q0(magnitudeImage,Rect(,,cx,cy));
   Mat q1(magnitudeImage,Rect(cx,,cx,cy));
   Mat q2(magnitudeImage,Rect(,cy,cx,cy));
   Mat q3(magnitudeImage,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(magnitudeImage,magnitudeImage,,,NORM_MINMAX);

   /******显示效果图******/
   imshow("magnitudeImage",magnitudeImage);

   waitKey();
   return ;
}

运行结果:

                        

效果图                                                                                              原图