数字图像处理--图像的傅里叶变换.doc

_实验名称图像的傅立叶变换课程名称数字图像姓 名专业、班级学 号班内序号实验时间实验地点实验内容1.傅立叶变换A)傅立叶变换基本操作B)利用MATLAB软件实现数字图像傅立叶变换的程序C)绘制一个二值图像矩阵,并将其傅立叶函数可视化2，离散余弦变换A) 使用dct2对图像autumn.tif进行DCT变换。B) 将上述DCT变换结果中绝对值小于10的系数舍弃，使用idct2重构图像并与原图像比较。C) 利用DCT变换进行图像压缩。实验过程及结果分析实验过程：1.傅立叶变换A) 傅立叶变换基本操作A=imread(N.JPG); I=rgb2gray(A);imshow(I);title(原图像);J=fft2(I);figure,imshow(J);title(傅里叶变换);JSH=fftshift(J);figure,imshow(JSH);title(傅里叶变换频移);Ji=ifft2(J);figure,imshow(Ji/256);title(直接傅里叶变换);JA=abs(J);iJA=ifft2(JA);figure,imshow(iJA/256);title(幅度傅里叶反变换);JP=angle(J);iJP=ifft2(JP);figure,imshow(abs(iJP)*100);title(相应的傅里叶反变换); B) 利用MATLAB软件实现数字图像傅立叶变换的程序B=imread(M.JPG); I=rgb2gray(B);imshow(I);fftI=fft2(I);sfftI=fftshift(fftI);RR=real(sfftI);II=imag(sfftI);A=sqrt(RR.2+II.2);A=(A-min(min(A)/(max(max(A)-min(min(A)*225;figure;imshow(A);C) 绘制一个二值图像矩阵,并将其傅立叶函数可视化f=zeros(30,30);f(5:24:13:17)=1;imshow(f,notruesize)F=fft2(f);F2=log(abs(F);figure,imshow(F2,-1 5,notruesize);colormap(jet);F=fft2(f,256,256);figure,imshow(log(abs(F),-1 5,notruesize);colormap(jet);F2=fftshift(F);figure,imshow(log(abs(F2),-1 5,notruesize);colormap(jet);1. 离散余弦变换A) 使用dct2对图像N.jpg进行DCT变换。RGB=imread(N.jpg);imshow(RGB)I=rgb2gray(RGB);figure,imshow(I)J=dct2(I);figure,imshow(log(abs(J),),colormap(jet(64);colorbar;B) 将上述DCT变换结果中绝对值小于10的系数舍弃，使用idct2重构图像并与原图像比较。RGB=imread(N.jpg);%BI=rgb2gray(RGB);J=dct2(I);figure,imshow(I);J(abs(J)10)=0;K=idct2(J);figure,imshow(K,0 255)J(abs(J)10)=0;K2=idct2(J);figure,imshow(K2,0 255)C) 利用DCT变换进行图像压缩。I=imread(N.jpg);I=im2double(I);T=dctmtx(8);B=blkproc(I,8,8,P1*x*P2,T,T);mask=1 1 1 1 0 0 0 0 1 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0B2=blkproc(B,8 8,P1.*x,mask);I2=blkp