OpenCV实现傅立叶变换源代码详解

基于OpenCV的傅里叶变换及逆变换
#include <stdio.h>
#include <cv.h>
#include <cxcore.h>
#include <highgui.h>

/**************************************************************************
//src IPL_DEPTH_8U
//dst IPL_DEPTH_64F
**************************************************************************/

//傅里叶正变换
void fft2(IplImage *src, IplImage *dst)
{   //实部、虚部
 IplImage *image_Re = 0, *image_Im = 0, *Fourier = 0;
 //   int i, j;
 image_Re = cvCreateImage(cvGetSize(src), IPL_DEPTH_64F, 1);  //实部
 //Imaginary part
 image_Im = cvCreateImage(cvGetSize(src), IPL_DEPTH_64F, 1);  //虚部
 //2 channels (image_Re, image_Im)
 Fourier = cvCreateImage(cvGetSize(src), IPL_DEPTH_64F, 2);
 // Real part conversion from u8 to 64f (double)
 cvConvertScale(src, image_Re, 1, 0);
 // Imaginary part (zeros)
 cvZero(image_Im);
 // Join real and imaginary parts and stock them in Fourier image
 cvMerge(image_Re, image_Im, 0, 0, Fourier);
 // Application of the forward Fourier transform
 cvDFT(Fourier, dst, CV_DXT_FORWARD);
 cvReleaseImage(&image_Re);