图像的正交变换matlab
《数字图像处理》课程实验报告
实验名 :图像的正交变换实验
1
院系: 自动化测试与控制系
班级: 1201132
姓名: 李丹阳
学号: 1120110113
哈尔滨工业大学
电气工程及自动化学院
光电信息工程
2015 年12月13 日
一、实验原理
二、实验内容
三、实验结果与分析
1、傅立叶变换
A) 绘制一个二值图像矩阵,并将其傅立叶函数可视化。
(傅里叶变换A)的实验结果
B)利用傅立叶变换分析两幅图像的相关性,定位图像特征。读入图像‘cameraman.tif’,
抽取其中的字母‘a’。
1
(傅里叶变换B)的实验结果
离散余弦变换(DCT)
使用dct2对图像 ‘linyichen.jpg’进行DCT变换。
5
0
-5
(离散余弦变换A)的实验结果
将上述DCT变换结果中绝对值小于 10的系数舍弃,使用idct2重构图像并与原图像
比较。
离散余弦变换
2
附主要程序代码:
f=zeros(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(jet);
B)利用傅立叶变换分析两幅图像的相关性,定位图像特征。读入图像‘cameraman.tif’,
抽取其中的字母‘a’。
bw=imread('cameraman.tif');
a=bw(59:71,81:91);
imshow(bw);
figure,imshow(a);
C=real(ifft2(fft2(bw).*fft2(rot90(a,2),256,256)));%求相关性
figure,imshow(C,[]);
thresh=max(C(:));
figure,imshow(C>thresh-10)
figure,imshow(C>thresh-15)
1. 离散余弦变换(DCT)
A) 使用dct2对图像 ‘linyichen.jpg’进行DCT变换。
RGB=imread('linyichen.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('linyichen.jpg');
I=rgb2gray(RGB); %转换为灰度图像
J=dct2(I);
figure,imshow(I)