资源描述：
《数字图像处理--图像的傅里叶变换》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验名称图像的傅立叶变换课程名称数字图像姓名专业、班级学号班内序号实验时间实验地点实验内容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)
2、;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(ab
3、s(iJP)*100);title('相应的傅里叶反变换');A)利用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);A)绘制一个二值图像矩阵,并将其傅立叶函数可视化f=ze
4、ros(30,30);f(5:24:13:17)=1;imshow(f,'notruesize')F=fft2(f);F2=log(abs(F));figure,imshow(F2,[-15],'notruesize');colormap(jet);F=fft2(f,256,256);figure,imshow(log(abs(F)),[-15],'notruesize');colormap(jet);F2=fftshift(F);figure,imshow(log(abs(F2)),[-15],'notruesize');colormap(je
5、t);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;A)将上述DCT变换结果中绝对值小于10的系数舍弃，使用idct2重构图像并与原图像比较。RGB=imread('N.jpg');%BI=rgb2gray(RGB);J=dct2(I);figure,imshow(I)
6、;J(abs(J)<10)=0;K=idct2(J);figure,imshow(K,[0255])J(abs(J)<10)=0;K2=idct2(J);figure,imshow(K2,[0255])B)利用DCT变换进行图像压缩。I=imread('N.jpg');I=im2double(I);T=dctmtx(8);B=blkproc(I,[8,8],'P1*x*P2',T,T');mask=[1111000011100000110000001000000000000000000000000000000000000000]B2=blkpro
7、c(B,[88],'P1.*x',mask);I2=blkproc(B2,[88],'P1.*x*P2','T',T);imshow(I)figure,imshow(I2)
