數(shù)字圖像處理實驗報告(編輯修改稿)
2025-02-14 22:32 本頁面
【文章內(nèi)容簡介】 數(shù)對其進行FFT：① 同屏顯示原圖和的幅度譜圖；② 若令，重復(fù)以上過程，比較二者幅度譜的異同，簡述理由； 實驗圖象f1(x, y)③ 若將順時針旋轉(zhuǎn)45度得到，試顯示的幅度譜，并與的幅度譜進行比較。1.%生成圖形f1f1=zeros(128,128)。f1((6430):(63+30),(6410):(63+10))=1。%FFT變換fft_f1=log(1+abs(fftshift(fft2(f1))))。figure。subplot(121)。imshow(f1)。title(39。Image f139。)。subplot(122)。imshow(fft_f1,[])。title(39。FFT f139。)。2.%計算f2f2=zeros(128,128)。for i=1:128。 for j=1:128。 f2(i,j)=((1)^(i+j))*f1(i,j)。 endendfft_f2A=log(1+abs(fft2(f2)))。fft_f2B=log(1+abs(fftshift(fft2(f2))))。figure。subplot(131)。imshow(f2)。title(39。Image f239。)。subplot(132)。imshow(fft_f2B,[])。title(39。FFT f239。)。subplot(133)。imshow(fft_f2A,[])。title(39。FFT f2 Without FFTShift39。)。分析：根據(jù)傅里葉變換對的平移性質(zhì)：；當且時，有：因此可得到：所以，就是頻譜中心化后的結(jié)果。3.%計算f3f3=imrotate(f2,45,39。nearest39。)。fft_f3=log(1+abs(fftshift(fft2(f3))))。figure。subplot(121)。imshow(f3)。title(39。Image f339。)。subplot(122)。imshow(fft_f3,[])。title(39。FFT f339。)。（2）（256256大小、256級灰度）進行頻域的理想低通、高通濾波，同屏顯示原圖、幅度譜圖和低通、高通濾波的結(jié)果圖。 低通濾波：fid=fopen(39。D:\matlab7\image\39。,39。r39。)。data=(fread(fid,[256,256],39。uint839。))39。subplot(1,2,1)imagesc(data)。colormap(gray)。title(39。LENA39。,39。Color39。,39。r39。)。fft_lena=fft2(data)。f=fftshift(fft_lena)。for i=1:256for j=1:256 if sqrt((i128)^2+(j128)^2)30 f(i,j)=0。 end endendsubplot(1,2,2)。 [x,y]=meshgrid(1:1:256)。surf(x,y,f)高通濾波：fid=fopen(39。D:\matlab7\image\39。,39。r39。)。data=(fread(fid,[256,256],39。uint839。))39。subplot(1,2,1)imagesc(data)。colormap(gray)
公司管理相關(guān)推薦
鄂ICP備17016276號-1