資源描述：
《圖像處理分析_FFT》由會員上傳分享，免費在線閱讀，更多相關內(nèi)容在工程資料-天天文庫。
1、實驗名稱實驗5：頻域濾波實驗目的掌握圖像進行頻域濾波的方法和步驟。1、掌握圖像頻域DFT變換和反變換的方法。2、掌握圖像頻域濾波的步驟實驗內(nèi)容1、灰度圖像的DFT和IDFTo具體內(nèi)容：利用OpenCV提供的cvDFT函數(shù)對圖像進行DFT和IDFT變換2、利用理想高通和低通濾波器對灰度圖像進行頻域濾波具體內(nèi)容：利用cvDFT函數(shù)實現(xiàn)DFT,在頻域上利用理想高通和低通濾波器進行濾波，并把濾波過后的圖像顯示在屏幕上（觀察振鈴現(xiàn)象），要求截止頻率可輸入。3、利用布特沃斯高通和低通濾波器對灰度圖像進行頻域濾波。具體內(nèi)容：利用cvDFT函數(shù)實現(xiàn)DFT,
2、在頻域上進行利用布特沃斯高通和低通濾波器進行濾波，并把濾波過后的圖像顯示在屏幕上（觀察振鈴現(xiàn)象），要求截止頻率和n可輸入。實驗完成情況1、實驗步驟：利用OpenCV提供的cvDFT函數(shù)對圖像進行DFT和IDFT變換核心代碼如下：//DFT變換Ipllmage*DFT(IplImage*src){Ipllmage*fourier=cvCreateImage(cvGetSize(src),IPL_DEPTH_64F,2);intdft_H,dft_W;dft_H=src->height;dft_W=src->width;CvMat*src_Re
3、=cvCreateMat(dft_H,dft_W,CV_64FC1);//doubleRe,Im;CvMat*src」m=cvCreateMat(dft_H,dft_W,CV_64FC1);//ImaginarypartCvMat*sum_src=cvCreateMat(dft_H,dft_W,CV_64FC2);//2channels(src_Re,src_Im)CvMat*sum_dst=cvCreateMat(dft_H,dft_W,CV_64FC2);//2channels(dst_Re,dst_Im)cvConvert(src,s
4、rc_Re);cvZero(src_Im);cvMerge(src_Re,src_Im,0,0,sum_src);cvDFT(sum_src,sum_dst,CV_DXT_FORWARD,0);cvConvert(sum_dst,fourier);cvReleaseMat(&src_Re);cvReleaseMat(&src_【m);cvReleaseMat(&sum_s『c);cvReleaseMat(&sum_dst);returnfourier;}//DFT反變換Ipllmage*IDFT(IplImage*fourier){Ipll
5、mage*dst=cvCreateImage(cvGetSize(fourier),IPL_DEPTH_8U,1);intdft_H,dft_W;dft_H=fourier->height;dft_W=fourier->width;CvMat*dst_Re二cvCreateMat(dft_H,dft_W,CV_64FC1);//doubleRe,Im;CvMat*dst_Im=cvCreateMat(dft_H,dft_W,CV_64FC1);//ImaginarypartCvMat*sum_dst=cvCreateMat(dft_H,df
6、t_W,CV_64FC2);//2channels(dst_Re,dst_Im)CvMat*sum_src=cvCreateMat(dft_H,dft_W,CV_64FC2);cvConvert(fourier,sum_src);cvDFT(sum_src,sum_dst,CV_DXT_INV_SCALE,O);cvSplit(sum_dst,dst_Re,dst_Im,O,O);cvConvert(dst_Re,dst);cvReleaseMat(&dst_Re);cvReleaseMat(&dst_Im);cvReleaseMat(&s
7、um_src);cvReleaseMat(&sum_dst);returndst;1〃歸一化，將灰度映射到0~255之間，并將能量最高的四角移到中心，生成圖片頻域能量圖voidBuildDFTImage(IplImage*fourier,Ipllmage*dst){Ipllmage*image_Re=0,*image_Im=0;image_Re=cvCreateImage(cvGetSize(fourier),IPL_DEPTH_64F,1);image_Im=cvCreateImage(cvGetSize(fourier),1PL_DEP
8、TH_64F,1);//ImaginarypartcvSplit(fourier,image_Re,image_Im,0,0);//Computethemagnitudeofth
