(matlab)數(shù)字圖像處理實(shí)驗(yàn)報(bào)告

1、重 慶 交 通 大 學(xué)學(xué) 生 實(shí) 驗(yàn) 報(bào) 告實(shí)驗(yàn)課程名稱數(shù)字圖像處理課程上機(jī)實(shí)驗(yàn)開課實(shí)驗(yàn)室河海學(xué)院仿真實(shí)驗(yàn)室學(xué)院河海學(xué)院年 級(jí) 專 業(yè)08 級(jí)地理信息系統(tǒng)學(xué) 生 姓 名學(xué)號(hào)08260129 時(shí)間2011 至2012 學(xué)年第1 學(xué)期實(shí)驗(yàn)一圖像顯示【實(shí)驗(yàn)內(nèi)容】1）使用 matlab圖像讀取函數(shù)imread 讀取圖像。2）使用 matlab圖像顯示函數(shù)imshow 顯示圖像。3）使用 matlab添加色帶函數(shù)colorbar為圖像添加色帶。【實(shí)驗(yàn)?zāi)康摹?）掌握 matlab 圖像讀取和顯示函數(shù)的應(yīng)用方法。2）了解如何為圖像添加色帶?！緦?shí)驗(yàn)結(jié)果】（放置處理前圖像）（放置處理后圖像）2-1 2-5-3

2、2-10 【程序說明】a=imread(yq.jpg); a=double(a); %a=uint8(a); imshow(a); %save saturn.dat a -ascii; save yu.text a -ascii; % 結(jié)果圖 2-1 e=imread(yq.jpg); imshow(e); iptsetpref(imshowtruesize,manual); figure,imshow(e); iptsetpref(imshowtruesize,auto); bw1=zeros(20,20); bw1(2:2:18,2:2:18)=1; figure,imshow(bw1,n

3、 otruesize); bw1 whos % 結(jié)果圖 2-5-3 % 使用一個(gè)調(diào)色板來顯示一副二進(jìn)制圖像figure,imshow(bw,1 0 0;0 0 1); % 結(jié)果圖 2-10 本次實(shí)驗(yàn)得分實(shí)驗(yàn)二圖像運(yùn)算【實(shí)驗(yàn)內(nèi)容】1）使用 matlab滑動(dòng)鄰域操作函數(shù)nlfilter對(duì)圖像進(jìn)行處理。2）使用 matlab分離鄰域操作函數(shù)blkproc對(duì)圖像進(jìn)行處理?！緦?shí)驗(yàn)?zāi)康摹?）掌握滑動(dòng)鄰域操作函數(shù)的應(yīng)用方法。2）掌握分離鄰域操作函數(shù)的應(yīng)用方法?！緦?shí)驗(yàn)結(jié)果】（放置處理前圖像）（放置處理后圖像） 2-14jpg 【程序說明】（放置程序段）%f3.23(滑動(dòng)鄰域操作）i=imread(2-14.j

4、pg); f=inline(max(x(:); i2=nlfilter(i,3 3,f); subplot(1,2,1),imshow(i); subplot(1,2,2),imshow(i2); %結(jié)果圖： 3-23.jpg %f3.25(分離鄰域操作 1）i=imread(2-14.jpg); f=inline(uint8(round(mean2(x)*ones(size(x); i2=blkproc(i,8 8,f); i3=blkproc(i,1 8,f); subplot(1,3,1),imshow(i); subplot(1,3,2),imshow(i2); subplot(1,3

5、,3),imshow(i3); %結(jié)果圖： 3-25.jpg 本次實(shí)驗(yàn)得分實(shí)驗(yàn)三圖像變換【實(shí)驗(yàn)內(nèi)容】1）構(gòu)造一矩形或正方形或菱形或其他形狀的二進(jìn)制圖像。2）使用 matlab二維傅里葉變換函數(shù)fft2對(duì)圖像進(jìn)行處理?！緦?shí)驗(yàn)?zāi)康摹?）掌握二進(jìn)制圖像的構(gòu)造方法。2）掌握二維傅里葉變換函數(shù)的應(yīng)用方法?！緦?shí)驗(yàn)結(jié)果】（放置處理前圖像）無（放置處理后圖像）4-2.jpg 4-4.jpg 4-5.jpg 4-6.jpg【二維傅里葉變換數(shù)學(xué)計(jì)算】（放置數(shù)學(xué)公式）c=ifft2(fft2(a).*fft2(b); 【程序說明】（放置程序段）%fig_4.6 f=zeros(40,40); f(6:24,13:1

6、7)=1; imshow(f,notruesize); % 結(jié)果圖： 4-2.jpg %fig_4.7 f=fft2(f); f2=log(abs(f); imshow(f2,-1 5,notruesize); colormap(jet); % 結(jié)果圖： 4-4.jpg %fig_4.8 f=fft2(f,256,256); imshow(log(abs(f),-1 5); colormap(jet); % 結(jié)果圖 :4-5.jpg %fig_4.3 f=fft2(f,256,256); f2=fftshift(f); imshow(log(abs(f2),-1 5); colormap(je

7、t); % 結(jié)果圖 :4-6.jpg 本次實(shí)驗(yàn)得分實(shí)驗(yàn)四濾波和濾波器設(shè)計(jì)【實(shí)驗(yàn)內(nèi)容】1）應(yīng)用 matlab二維窗口設(shè)計(jì)函數(shù)fwind2創(chuàng)建一個(gè)近似圓對(duì)稱的帶通濾波器，通頻帶為 0.1 0.5 ，輸出濾波器的頻率響應(yīng)曲線?！緦?shí)驗(yàn)?zāi)康摹?）了解 matlab 濾波器設(shè)計(jì)方法。2）掌握用二維窗口設(shè)計(jì)函數(shù)fwind2 創(chuàng)建帶通濾波器的方法?！緦?shí)驗(yàn)結(jié)果】（放置處理后圖像）【程序說明】（放置程序段） f1,f2=freqspace(34,meshgrid); hd1=ones(34); r=sqrt(f1.2+f2.2); hd1(r0.5)=0; win=fspecial(gaussian,34,2);

8、 win=win./max(win(:); h3=fwind2(hd1,win); figure,freqz2(h3); % 結(jié)果圖： 5-9.jpg 本次實(shí)驗(yàn)得分實(shí)驗(yàn)五二值形態(tài)學(xué)操作【實(shí)驗(yàn)內(nèi)容】1）應(yīng)用 matlab圖像膨脹函數(shù)imdilate對(duì)圖像進(jìn)行膨脹操作。2）應(yīng)用 matlab圖像腐蝕函數(shù)imreode 對(duì)圖像進(jìn)行腐蝕操作?！緦?shí)驗(yàn)?zāi)康摹?）了解 matlab 圖像膨脹和圖像腐蝕基本算法。2）掌握 matlab 圖像膨脹imdilate和圖像腐蝕函數(shù)imreode 的使用方法。【實(shí)驗(yàn)結(jié)果】（放置處理后圖像）6-1.jpg 6-3.jpg 【圖像膨脹和圖像腐蝕數(shù)學(xué)計(jì)算】（放置數(shù)學(xué)公式）【

9、程序說明 】（放置程序段）bw=zeros(9,10); bw(5:6,5:7)=1; se=strel(square,3); bw2=imdilate(bw,se); subplot(1,2,1),imshow(bw); subplot(1,2,2),imshow(bw2); % 結(jié)果圖： 6-1.jpg se=strel(rectangle,40 30); bw1=imread(6-2.bmp); bw2=imerode(bw1,se); imshow(bw2); % 結(jié)果圖： 6-3.jpg 本次實(shí)驗(yàn)得分實(shí)驗(yàn)六圖像的空間變換【實(shí)驗(yàn)內(nèi)容】1）讀取兩幅圖像。2）對(duì)兩幅圖像進(jìn)行匹配?！緦?shí)驗(yàn)?zāi)康?/p>

10、】1）掌握 matlab 圖像匹配方法。【實(shí)驗(yàn)結(jié)果】（放置處理前圖像）（放置處理后圖像）7-9.jpg 7-10.jpg 【程序說明】（放置程序段）orthophoto=imread(23_s.jpg); subplot(1,2,1),imshow(orthophoto) unregistered=imread(23_s.jpg); subplot(1,2,2),imshow(unregistered) cpselect(unregistered(:,:,1),orthophoto) input_points=120.7086 93.9772 319.2222 78.9202 127.9838

11、 291.6312 352.0729 281.1445; input_points_corr=cpcorr(input_points,base_points,. unregistered(:,:,1),orthophoto) input_points_corr=120.7086 93.9772 本次實(shí)驗(yàn)得分319.2222 78.9202 127.1064 289.8935 352.0729 281.1445; mytform=cp2tform(input_points,base_points,projective); registered=imtransform(unregistered,m

12、ytform) %結(jié)果圖： 7-9.jpg %ch7_8_11 lily=imread(7-8.jpg); flowers=imread(7-8.jpg); subplot(1,2,1),imshow(lily); subplot(1,2,2),imshow(flowers); % rect_lily=93 13 81 69; rect_flowers=190 68 235 210; sub_lily=imcrop(lily,rect_lily); sub_flowers=imcrop(flowers,rect_flowers); sub_lily1,rect_lily=imcrop(lily

13、); sub_flowers1,rect_flowers=imcrop(flowers); subplot(1,2,1),imshow(sub_lily1) subplot(1,2,2),imshow(sub_flowers1) % c=normxcorr2(sub_lily(:,:,1),sub_flowers(:,:,1); figure,surf(c),shading flat % sub_lily,rect_lily=imcrop(lily); sub_flowers,rect_flowers=imcrop(flowers); % max_c,imax=max(abs(c(:); yp

14、eak,xpeak=ind2sub(size(c),imax(1); corr_offset=(xpeak-size(sub_lily,2) (ypeak-size(sub_lily,1); rect_offset=(rect_flowers(1)-rect_lily(1) (rect_flowers(2)-rect_lily(2); offset=corr_offset+rect_offset; xoffset=offset(1); yoffset=offset(2); % xbegin=xoffset + 1; xend=xoffset + size(lily,2); ybegin=yof

15、fset+1; yend=yoffset+size(lily,1); extracted_lily=flowers(ybegin:yend,xbegin:xend,:); if isequal(lily,extracted_lily) disp(lily.tif was extracted from flowers.tif) end % recovered_lily=uint8(zeros(size(flowers); recovered_lily(ybegin:yend,xbegin:xend,:)=lily; m,n,p=size(flowers) mask=ones(m,n); i=fi

16、nd(recovered_lily(:,:,1)=0); mask(i)=.2; figure,imshow(flowers(:,:,1) hold on h=imshow(recovered_lily); set(h,alphadata,mask) %結(jié)果圖： 7-10.jpg 實(shí)驗(yàn)七圖像增強(qiáng)【實(shí)驗(yàn)內(nèi)容】1）讀取一副圖像。2）應(yīng)用 matlab直方圖均衡化函數(shù)histeq對(duì)圖像進(jìn)行均衡化操作。【實(shí)驗(yàn)?zāi)康摹?）掌握 matlab 直方圖均衡化方法?！緦?shí)驗(yàn)結(jié)果】（放置處理后圖像）8-4-2.jpg 8-11.jpg 【程序說明 】（放置程序段）（放置處理前后圖像的數(shù)據(jù)矩陣）i=imread(84.jpg); imshow(i) figure,