温州大学计算机网络,数据通信与计算机网络（温州大学）数字图像处理及Matlab实现.doc...

温州大学

WENZHOU UNIVERSITY

实 验 报 告

学 期2015-2016学年第一学期课 程数字图像处理读取图象，显示这幅图象，对图象作傅立叶变换，显示频域振幅图象。作傅立叶逆变换，显示频域振幅图象看是否与原图象相同Ia=imread('D:\girl.jpg');%?-í???

I=rgb2gray(Ia);%oú°×í???

J=fft2(I);%?μà?ò?±??ˉ

F=abs(J);%·ù?è?×

F1=angle(J);%?à???×

subplot(2,4,1);

imshow(Ia);title('?-í???');

subplot(2,4,2);

imshow(I);title('?ò?èí???');

subplot(2,4,3);

imshow(J);title('?μà?ò?í?');

subplot(2,4,4);

imshow(F);title('·ù?èí?');

subplot(2,4,5);

imshow(F1);title('?à??í?');

K=ifft2(I);%?μà?ò?·′±???

IJ=abs(ifft2(I));%·ù?è?×

RJ=angle(ifft2(I));%?à???×

subplot(2,4,6);

imshow(K);title('?μà?ò?·′±???í?');

subplot(2,4,7);

imshow(IJ);title('·′±???·ù?è?×');

subplot(2,4,8);

imshow(RJ);title('·′±????à??í?');

设计一个低通滤波器，截止频率自选，对图象作低通滤波，再作反变换，观察不同的截止频率下反变换后的图象与原图象的区别I=imread('D:\girl.jpg');

I=rgb2gray(I);

I1=fftshift(fft2(I));

[M,N]=size(I1);

n=1;

d0=50;

n1=floor(M/2);

n2=floor(N/2);

for i=1:M

for j=1:N

d=sqrt((i-n1)^2+(j-n2)^2);

H=1/(1+(d/d0)^(2*n));

I2(i,j)=H*I1(i,j);

end

end

I2=ifftshift(I2);

I3=real(ifft2(I2));

subplot(2,2,1);

imshow(I);title('?-í????ò?è′|àí');

subplot(2,2,2);

imshow(I3,[]);title('??2¨oóμ?í???');

subplot(2,2,3),hist(double(I),16),title('?-í?μ??±·?í?');

subplot(2,2,4),hist(double(I3),16),title('??2¨oóμ??±·?í?');

变化后的图像灰度直方图灰度值更加集中。

图像边沿的模糊程度减小。

3.请在excel表中绘制上述各灰度级以及均衡化后的概率密度直方图与累加密度，截图到word里面。

灰度等级原密度累加密度密度等级规划等级规划密度00.1950.1950S3010.160.3550.066667S5020.1470.5020.133333S8030.1060.6080.2s90.19540.0730.6810.266667s10050.0560.7370.333333s110.1660.0490.7860.4s12070.0410.8270.466667s12080.0370.8640.533333s130.14790.0310.8950.6s130.106100.0270.9220.666667s140.073110.0230.9450.733333s140.056120.0190.9640.8s140.09130.0170.9810.866667s150.068140.0120.9930.933333s150.069150.00711s150.036

三、数据记录与分析

四、思考与讨论

指导师(签名) 实验报告成绩 批改日期

作答记录