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

实验名称 图像的傅立叶变换 课程名称 数字图像 姓 名 专业、班级 学 号 班内序号 实验时间 实验地点 实 验 内 容 1. 傅立叶变换 A) 傅立叶变换基本操作 B) 利用MATLAB软件实现数字图像傅立叶变换的程序 C) 绘制一个二值图像矩阵,并将其傅立叶函数可视化 2，离散余弦变换 使用dct2对图像‘autumn.tif’进行DCT变换。 将上述DCT变换结果中绝对值小于10的系数舍弃，使用idct2重构图像并与原图像比较。 利用DCT变换进行图像压缩。 实 验 过 程 及 结 果 分 析 实验过程： 1.傅立叶变换 傅立叶变换基本操作 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(相应的傅里叶反变换); 利用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); 绘制一个二值图像矩阵,并将其傅立叶函数可视化 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); 离散余弦变换 使用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; 将上述DCT变换结果中绝对值小于10的系数舍弃，使用idct2重构图像并与原图像比较。 RGB=imread(N.jpg);%B I=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]) 利用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 0] B2=blkproc(B,[8 8],P1.*x,mask); I2=blkproc(B2,[8 8],P1.*x*P2,T,T); imshow(I) figure,imshow(I2)