數(shù)字圖像處理MATLAB上機(jī)作業(yè)

1、數(shù)字圖像處理實(shí)驗(yàn)報(bào)告指導(dǎo)老師:學(xué)號姓名班級1.產(chǎn)生右圖所示圖像fl(m,n),其中圖像大小為256X256,中間亮條為128X32,暗處=0,亮處=100對其進(jìn)行FFT(f2)和 FFT(f3)的關(guān)系,FFT：同屏顯示原圖fl(m,n)和FFT(fl)的幅度譜圖； 若令f2(m,n)=(T)(m+n)fl(m,n),重復(fù)以上過程,比較二者幅度譜的異同，簡述理由； 若將f2(m,n)順時(shí)針旋轉(zhuǎn)90度得到f3(m,n),試顯示FFT(f3)的幅度譜，并與FFT(f2)的幅度譜進(jìn)行比較； 若將fl(m,n)順時(shí)針旋轉(zhuǎn)90度得到f4(m,n),令f5(m,n)=f1(m,n)+f4(m,n),試顯示F

2、FT(f5)的幅度譜，并指出其與FFT(fl)和FFT(f4)的關(guān)系；若令f6(m,n)=f2(m,n)+f3(m,n),試顯示FFT(f6)的幅度譜,并指出其與比較FFT(f6)和FFT(f5)的幅度譜。fl=zeros(256,256);for1=64:1:192forj=122:l:144fl(i,j)=100;endfft_fl=fft2(fl);fft_fl=abs(fft_fl);tmax=fft_fl(1,1);tmin=fft_fl(1,1);fori=l:256forj=l:256iftmax<fft_fl(i,j)tmax=fft_fl(i,j);endiftmin&

3、gt;fft_fl(i,j)tmin=fft_fl(i,j);enddelta=tmax-tmin;fori=l:256forj=l:256j)=255*(fft_fl(i,n)/delta;endendsubplot(1,2,1);imshow(fl);title('原圖')；subplot(1,2,2);imshow(fft_fl);titleC原圖的幅度譜，)：原因I的幅度譜fori=l:256forj=l:256f2(i,j)=(-l)F+j)*fl(i.j);endendfft_f2=fft2(f2);fft_f2=abs(fft_f2);tmax=fft_f2(l,

4、1);tmin=fft_f2(l,1);fori=l:256forj=l:256iftmax<fft_f2(i,j)tmax=fft_f2(i,j);endiftmin>fft_f2(i,j)tmin=fft_f2(i,j);endendenddelta=tmax-tmin;fori=l:256forj=l:256fft_f2(i,j)=255*(fft_f2(i,j)-tmin)/delta;endendsubplot(2,2,1);imshow(fl);title('原圖);subplot(2,2,2);title(原圖的幅度譜，)：subplot(2,2,3);ims

5、how(f2);title。原圖中心化，)：subplot(2,2,4);imshow(fft_f2);title(原圖中心化的幅度譜，):原IS原IS的幅度譜原直I中心化原圖中心化的幅度譜bili near');f3=imrotate(f2,-90,fft_f3=fft2(f3);fft_f3=abs(fft_f3);tmax=fft_f3(l,1);tmin=fft_f3(l,1);fori=l:256forj=l:256iftmax<fft_f3(i,j)tmax=fft_f3(i,j);endiftmin>fft_f3(i,j)tmin=fft_f3(i,j);en

6、dendenddelta=tmax-tmin;fori=l:256forj=l:256fft_f3(i,j)=255*(fft_f3(i,n)/delta;endendsubplot(2,2,1);imshow(f2);title。原圖中心化，)：subplot(2,2,2);imshow(fft_f2);title。原圖中心化的幅度譜，):subplot(2,2,3);imshow(f3);titleC旋轉(zhuǎn)后的圖像，)；subplot(2,2,4);imshow(fft_f3);title(，旋轉(zhuǎn)后的幅度譜，);原圖中心化的幅度譜旋轉(zhuǎn)后的幅度譜圖像旋轉(zhuǎn)90度后，幅度譜也旋轉(zhuǎn)了90度。bili

7、 near* );f4=imrotate(fl,-90,f5=fl+f4;fft_f5=fft2(f5);fft.f5=abs(fft.f5);tmax=fft_f5(l,1);tmin=fft_f5(l,1);fori=l:256forj=l:256iftmax<fft_f5(i,j)tmax=fft_f5(itj);endiftmin>fft_f5(i,j)tmin=fft_f5(i,j);endendenddelta=tmax-tmin;fori=l:256forj=l:256fft.f5(i,j)=255*(fft_f5(i,j)-tmin)/delta;endendsub

8、plot(2,2,1);imshow(fl);，原圖)：subplot(2,2,2);imshow(fft_f1);title。原圖的幅度譜，)：subplot(2,2,3);imshow(f5);titleC相加后的圖像，)：subplot(2,2,4);imshow(fft_f5);title。相加后的幅度譜，)；原國的幅度譜相加后的畫像相加后的幅度譜原圖與旋轉(zhuǎn)90度后的圖像進(jìn)行疊加后，相應(yīng)的幅度譜也進(jìn)行疊加。f6=f2+f3;fft_f6=fft2(f6);fft_f6=abs(fft_f6);tmax=fft_f6(l,1);tmin=fft_f6(L1)；fori=l:256forj

9、=l:256iftmax<fft_f6(i,j)tmax=fft_f6(i,j);endiftmin>fft_f6(i,j)tmin=fft_f6(i,j);endendenddelta=tmax-tmin;fori=l:256forj=l:256fft_f6(i,j)=255*(fft_f6(i,n)/delta;endendsubplot(2,2,1);imshow(fft_f2);title(?F2*);subplot(2,2,2);imshow(fft_f3);title(3,);subplot(2,2,3);imshow(fft_f5);title('25'

10、;);subplot(2,2,4);imshow(fft_f6);title('F6');F3F2FFT和FFT(3)相差90度。FFT(6)是FFT(5)的中心化。2.產(chǎn)生教材104頁題圖4.18(右圖)所示的二值圖像(白為1,黑為0),編程實(shí)現(xiàn)習(xí)題4.18所要求的處理(3*3的平均濾波和中值濾波)功能(圖像四周邊界不考慮，處理結(jié)果按四舍五入仍取?；?),顯示處理前后的圖像，比較其異同。a二zeros(64,64);f=zeros(256,256);fori=l:l:32forj=l:1:32a(i,j)=l：endendfori=33:1:64forj=33:1:64a(i

11、,j)=l；endendfori=l:64:256forj=l:64:256f(i:i+63,j:j+63)=a;endendimshow(f)fori=2:1:255forj=2:1:255j-l),f(i-l,j),f(i-l,j+l),f(i,j),f(i,j+l),f(i+l,j-1),f(i-l,j),f(i+Lj+l):b=f(i-l,j),f(i,j-l),f(i,j),f(i,j+1),f(i+1,j)；c=sort(a);d=sort(b);zlfl(ifj)=c(5);zlf2(i,j)=d(5);endendsubplot(1,3,1);imshow(f);title(原

12、圖像，);subplot(1,3,2);imshow(zlfl);titleC方形中值濾波后圖像，)subplot(1,3,3);imshow(zlf2);titleC十字形中值濾波后圖像)原IS像方形中值濾波后圖像十字形中值濾波后圖像fori=2:1:255forj=2:1:255jlf4(i,j)=(f(i-1,j)+f(i+1,j)+f(i,(i,j+l)/4;jlf8(i,j)=(f(i-l,j-l)+f(i-l,j)+f(i-l,j+l)+f(i,j-l)+f(i,j+l)+f(i+1,j-l)+f(i+l,j)+f(i+l,j+l)/8；jjlf4(i,j)=(f(i-l,j)+f

13、(i+l,j)+f(i,j-l)+f(itj+D+f(i,j)/5;jjlf8(i,j)=(f(i-l,j)+f(i-l,j+l)+f(i,j-l)+f(i,j+l)+f(i+ltj-l)+f(i+l,j)+f(i+Lj+l)+f(i,j)/9;endendjlf4=round(jlf4);jlf8=round(jlf8);jjlf4=round(jjlf4);jjlf8=round(jjlf8);subplot(2,2,1);imshow(jlf4);title(，4鄰域平均濾波后圖像，)：subplot(2,2,2);imshow(jlf8);鄰域平均濾波后圖像，)subplot(2,2,

14、3);imshow(jjlf4);title。4鄰域加權(quán)平均濾波后圖像，)subplot(2,2,4);imshow(jjlf8);title。8鄰域加權(quán)平均濾波后圖像')4鄰域加權(quán)平均濾波后11像塔鄰域加權(quán)平均濾波后圖像放大后可看出區(qū)別4鄰域平均濾波后圖像8鄰域平均濾液后S像4鄰域加權(quán)平均濾波后BI像8鄰域加權(quán)平均濾波后BI像3產(chǎn)生教材104頁題圖4.16所示的灰度圖像(白為255,黑為0),分別加入高斯白噪聲和椒鹽噪聲，再分別進(jìn)行3*3的平均濾波和中值濾波，顯示原圖像、加噪圖像和濾波結(jié)果圖像，并比較四種濾波結(jié)果。f=zeros(256,256);a=255*ones(210,7);

15、fori=0:l:8f(24:233,(29+24*i):(29+24*1+6)=a;endfl=imnoise(f,'gaussian,0,0.1);先加入均fil為零方差為0.1的反J斯噪聲f2=imnoise(f,'salt&pepper',0.1);先加入噪聲密度為0.1的椒鹽噪聲fori=2:1:255forj=2:1:255j+1);b=f2(i-l,j),f2(i-l,j+1),f2(i,j-1),f2(i,j),f2(i,j+1),f2(i+l,j-1),f2(i+l,j)ff2(i+l,j+1);c=sort(a);d=sort(b);zlfl

16、(i,j)=c(5);氣中值濾波j)=cU5);%中值濾波jlf8(i,j)=(fl(i-l,j)+fl(i-l,j+l)+fl(i,j+l)+fl(i+l,j-l)+fl(i+l,j)+fl(i+l,j+l)+fl(i,j)/9;%平均濾波jjlf8(i,j)=(f2(i-l,j-l)+f2(i-l,j)+f2(i-l,j+1)+f2(i,j-l)+f2(i,j+l)+f2(i+l,j-l)+f2(i+l,j)+f2(i+l,j+l)+f2(i,j)/9;%平均濾波endendsubplot(2,4,1);imshow(f);title(原圖像，);subplot(2,4,2);imshow

17、(fl);title。高斯噪聲污染后的圖像，)：subplot(2,4,3);imshow(zlfl);titleC中值濾波后的圖像，)；subplot(2,4,4);imshow(jlf8);titleC平均濾波后的圖像，)；subplot(2,4,5);imshow(f);title(，原圖像，)；subplot(2,4,6);imshow(f2);title(，椒鹽噪聲污染后的圖像')：subplot(2,4,7);imshow(zlf2);title（中值濾波后的圖像，）:subplot（2, 4, 8）;imshow（jjlf8）;title（平均濾波后的圖像，）;高斯嗓再污

18、染后的18像中值速盈后的困怫平均遽戒后的Iffl擦原圖像根鹽唾聲污染石的困像中值遙渡后的Iffl像平均遮波石的圖像5.（第4、5題選做一題）編程實(shí)現(xiàn)教材214頁所給圖像門限化分割的迭代閾值算法，實(shí)現(xiàn)對某一灰度圖像的二值化。Pic=imread(C:UsersAdministratorDesktop242dd42a2834349bfb0c018ccbeal5ce36d3beb5.j九先讀取RGB格式的圖像a=rgb2gray(Pic);與進(jìn)行RGB到灰度圖像的轉(zhuǎn)換T=0.5*(double(min(a(:)+double(max(a(:);do=false;whiled。g=(a>T)；T

19、next=0.5*(mean(a(g)+mean(a(、g);do=abs(T-TnextXO.2;T=Tnext;endfigure;subplot(1,2,1);imshow(a);subplot(1,2,2);imshow(g);6.%ReadInputRetinaImageinlmg=imread('Input.bmp');if(dim=3)dim二ndims(inlmg);%lnputisacolorimageinImg二rgb2gray(inlmg);end%ExtractBloodVesselsThreshold=10;bloodvessels=VesselExtract(inlmg.Threshold);%0utputBloodVesselsimagefigure;Bloodsubplot(121);imshow(inlmg);title('InputImage');subplot(122);imshow(bloodVessels);title('ExtractedVessels");VesselExtract.m:functionbloodvessels二VesselExtract(inlmg,threshold)/Kirsch'sTemplateshl=5