图像处理实验三报告.docxVIP

Lab 3 【实验目的】 1、 了解傅立叶变换在图像处理中的应用； 2、 掌握Mat lab对图像进行滤波处理，并对不同类型的结果进行分析; 【实验内容】 2 Exercise 1 一 2?D Fourier Transforms Using the image files hardware?￡if、 shuttle.tif, and xray.tif, read in each of the images and coiiiputc the 2-D DFT lnagnitucle for cacli of the images?()】i a single plot (using a grid of 2x2 images), plot the following: ? the original gray scale image in the upper left cell the image represent at ion of tlie 2-D DFT magnitude of tlic image being studied in tlic iqjper right cell a cli])pcd and scaled version of the 2-D DFT magnitude that is clipped and scale at a level wlicrc you can sec the nature of the 2-D DFT magnitude for each image; this plot should be ])lcicecl in the lower left cell ? the log transformed 2-D DFT magnitude plotted again on a scale tliat enables you to see the structure of tlic tiansfonn magnitude; this plot should be placed in the lower right cell You will need to use the Matlab routine FS=fftshift (F) to shift the DC magnitude point from the up])er left corner of the image to tho center of the magnitude range. In your Lal)3 lcpoT仁 inchido tlie MATLAB code along with image plots of the three iinagcs for wliicli you rcj^cat tlie analysis? Can you say liow the image properties (in a gross sense) arc seen in tlie 2-D DFT magnitude plots? f=imread(fhardware.tiff); f=double(f); subplot(221),imshow(f,[J) F=fft2(f); S=abs(F); subplot(222),imshow(S,[]) Fc=fftshift(angle(S)); subplot(223),imshow(Fc,[]) F=fftshift(S); S2=log(l+abs(F)); subplot(224),imshow(S2,[]) 4 Exercise 3 — Linear Filtering of Images Using the MATLAB routine H=lpf ilter(type, M, N, DO, n) for designing lowjjass filters, where the input arguments are: type=,gaussian;, ideal or btw (for Butterworth filters) M,N = dimensionality of filter frequency response DO = filter cutoff frequency (normalized to range [0, M] or [0 N]) n = filter order for Butterworth filters and the output H is the frequency response of the filter, write a MATLAB ni-file that accepts as input the f