matlab图像处理综合实验实验报告

实验二

实验名称：图像变换 实验目的：

（1）进一步对matlab的了解和使用； （2）学习如何在matlab中对数字图像的处理； 实验原理：

图像和其他信号一样，既能在空间域处理，也能在频率域处理。把图像信息从空域变换到频域，可以更好的分析加工处理。因为图像信息的频域处理具有如下特点： （1）能量守恒，但能量重新分配； （2）有利于提取图像的某些特征；

（3）正交变换具有能量集中作用，可以实现图像的高效压缩编码； （4）频域用快速算法； 实验器材：电脑 matlab软件 实验内容：

（1）了解正交变换基本概念

（2）掌握图像的离散傅里叶和离散余弦变换 （3）熟悉图像的沃尔什及哈达玛变换 实验过程如下： 傅里叶变换

I = imread('E:\\dog.jpg') I = rgb2gray(I)

J = fft2(I)

subplot(2,2,1),imshow(I),title('灰度变换') J = fftshift(J)

subplot(2,2,2),imshow(log(abs(J)),[]),title('傅里叶变换') J(abs(J)<5000)=0

subplot(2,2,3),imshow(log(abs(J)+eps),[]),title('滤波') J = ifftshift(J) K = ifft2(J)

subplot(2,2,4),imshow(K,[0 255]),title('傅里叶逆变换')

灰度变换傅里叶变换滤波傅里叶逆变换余弦变换： clear all;

RGB=imread('E:\\dog.jpg'); I=rgb2gray(RGB);

figure;imshow(I);title('灰度图像');

J=dct2(I);

figure;imshow(log(abs(J)),[]);

colormap(jet(64)),colorbar;title('二维离散余弦变换')

灰度图像二维离散余弦变换 1086420-2-4-6-8

-10沃尔什—哈达玛变换： I=zeros(2.^8);

I(2.^7-2.^4+1:2.^7+2.^4,2.^7-2.^4+1:2.^7+2.^4)=ones(2*2.^4);

subplot(1,2,1);colormap(gray(128)),imagesc(I); [m,n]=size(I) for k=1:n

wht(:,k)=hadamard(m)*I(:,k)/m; end for j=1:m

wh(:,j)=hadamard(n)*wht(j,:)'/n; end wh=wh';

subplot(1,2,2);

colormap(gray(128)),imagesc(wh);

50501001001501502002002505010015020025025050100150200250

Radon变换：

R=radon（I，theta）：I表示待处理的图像，theta表示Radon变换的方向角度，可以是一个标量或向量值。R的每一列对应图像I在theta某祎角度的Radon变换值。 I=zeros(120,120); I(25:75,25:75)=1; figure;imshow(I);

theta=0:180;[R,xp]=radon(I,theta);

figure;imshow(R,[],'Xdata',theta,'Ydata',xp,... 'InitialMagnification','fit')

xlabel('\\theta (degrees)')ylabel('x''') colormap(hot),colorbar

iptsetpref('ImshowAxesVisible','off')