2D-FFT及IFFT（C语言实现（转载）.pdf

2D-FFT及及IFFT （（C语语⾔⾔实实现现 （（转转载载）） 2D-FFT及IFFT （C语⾔实现） 图像处理 2010-06-01 16:48:52 阅读727 评论0 字号：⼤中中⼩ 订阅 前⾯编过2D-FFT的程序 现在把2D-IFFT的程序整合到⼀起 便于后⾯做图像变换反变换使 ⽤。 FFT与IFFT有如下关系： 相应的2D-FFT与2D-IFFT的关系如下： 所以可以利⽤⼀个FFT核⼼函数实现2D-FFT与2D-IFFT。代码如下： #include stdio.h #include stdli .h #include math.h #define intsize sizeof(int) #define complexsize sizeof(complex) #define PI 3.1415926 int *a,* ; int nLen,init_ nLen,mLen,init_ mLen,N,M; FILE *dataFile; typedef struct{ float real; float image; }complex ; complex *A ,*A_ In,*W; complex Add(complex, complex); complex Su (complex, complex); complex Mul(complex, complex); int calculate_ M(int); void reverse(int,int); void readData(); void fft(int,int); void Ifft(); void printResult_fft(); void printResult_ Ifft(); int main() { int i,j ; readData(); A = (complex *)malloc(complexsize*nLen); reverse(nLen,N); for(i=0; imLen; i++) { for(j =0; j nLen; j ++) { A[j].real = A_ In[i*nLen+ [j]].real; A[j].image = A_ In[i*nLen+ [j]].image; } fft(nLen,N); for(j =0; j nLen; j ++) { A_ In[i*nLen+j].real = A[j].real; A_ In[i*nLen+j].image = A[j].image; } } free(a); free( ); free(A); A = (complex *)malloc(complexsize*mLen); reverse(mLen,M); for(i=0; inLen; i++) { for(j =0; j mLen; j ++) { A[j].real = A_ In[ [j]*nLen+i].real; A[j].image = A_ In[ [j]*nLen+i].image; } fft(mLen,M); for(j =0; j mLen; j ++) { A_ In[j*nLen+i].real = A[j].real; A_ In[j*nLen+i].image = A[j].image; } } free(A);