实验二图像的频域变换处理.docxVIP

实验二图像的频域变换处理 一、 实验目的： 1、 掌握Fourier, DCT和Radon变换、反变换的算法实现，并初步理解FourierRadon 和DCT变换的物理意义。 2、 利用傅里叶变换、离散余弦变换处理图像。 二、 实验要求： 用Fourier、Radon和DCT变换算法对图像作二维Fourier、Radon和DCT变换，并 输出变换结果，分析图像频域和空域Z间的关系。 用Fourier、Radon和DCT反变换算法对图像作二维Fourier. Radon和DCT反变换。 分别选择频域低通、高通滤波器对同一幅图像进行滤波，观察分析滤波后的图像变 化特征。 三、 实验步骤： 用Fourier Radon和DCT变换算法，对图像作二维FourierRadon和DCT变换和 反变换(考虑图像为256*256,可以采用快速Fourier变换方法)，并比较反变换后 的图像与原图像的视觉效果； 用DCT变换实现图像的压缩，计算压缩比。 利用高斯滤波器对图像进行滤波处理后输出。 实验方法： 设计matlab程序，在计算机中处理图像，输出结果。 1?设计程序生成一幅图像图像，利用函数F=fft2(MmcF,n)对图像进行傅里叶变换。 利用函数imshow(log(abs(F),[])输出图像的幅频特性, 利用函数F二ifft2(FameJ完成傅里叶反变换。 利用函数B=dct2(tname\[m,n])对图像进行离散余弦变换。利用函数 imshow(log(abs⑻,[])输出图像的幅频特性， 利用函数B=idct2(name,)完成傅里叶反变换。 利用det变换对图像进行压缩处理，计算压缩前后的图像大小。 实验要求： 1、 进实验室前必须设计好程序。 2、 认真调试，记录结果。 3、 实验报告必须有实验结论分析。 实验设备： MATLAB实验环境，图像输入输出设备，图像处理设备。 实验结果； l=imread ( E: \qiaole. jpg); A=rgb2gray(I); subplot (1, 2, 1), imshow(A); xlabel (原图); B二dct2(A); subplot (1, 2, 2), imshow(log(abs (B)), [0 5]); xlabel (二维离散余弦变换)； Q「32 1 Rlr 2 VWw Teds Drilflop Help ■ Daa* * % % o ? x a a q ■口 I=imread (J E: \qiaole. jpg); subplot (1, 3, 1) ; imshow(T) ; xlabel (原图); fftl=fft2(l); subplot (1, 3, 2) ; imshow(fftl); xlabel ( FFT 变换); sfftl-fftshift(fftl); subplot (1, 3, 3) ; imshow(sfftl) ; xlabel (中心化); ED「w i _ He Teels Desktop HHp □ -Ji 4 \ x ? ? Z ? 2 0 U □口 Uffi FE?掛 中心化 Uffi FE?掛 中心化 I=imread(,E:\qiaole.jpg,); A=rgb2gray(I); subplot( 1,3 J ),imshow( A); xlabel(原 S); B=dct2(A); subplot(l ,3,2),imshow(log(abs(B)),[0 5]); xlabelC二维离散变换)； I=idct2(B )/200; subplot(l ,3,3),imshow(I);xlabel(反变换)； file trntti Teds Desktop MUow Help ■ l=imread(E:\qiaole.jpg); 匸 double ⑴; g=fft2(f); g=fftshift(g); [M,N]=size(g); d()= 1530100; m=fix(M/2);n=fix(N/2); for i=l:M forj=l:M d=sqrt((i-m)A2+(j-n)A2); if(d=dO); h=l； else h=0; end result(i,j)=h*g(i,j); end end result=ifftshift(result); Jl=ifft2(result); J2=uint8(real(J 1)); imshow(J2) title(理想低通滤波) 4x、 ? ? Z ?a 4 x、 ? ? Z ? a ED「w i He Fd? Teels Desktop ^Mow HHp Id Ji ■ 口 ? ? c ?… l=imread(E:\qiaole.jp