基于MATLAB語(yǔ)言的圖像幾何變換

1、實(shí)驗(yàn)總成績(jī):報(bào)告份數(shù):! 西安郵電學(xué)院i通信與信息工程學(xué)院1科研訓(xùn)練報(bào)告專業(yè)班級(jí)：通工0907班學(xué)生姓名：學(xué)號(hào)（班內(nèi)序號(hào)）：2012 年 4 月 11 h一. 實(shí)驗(yàn)?zāi)康?1) 查閱相關(guān)數(shù)字圖像處理的文獻(xiàn)資料，熟悉matlab編程環(huán)境；(2) 學(xué)習(xí)兒種常見的圖像兒何變換，并通過(guò)實(shí)驗(yàn)體會(huì)兒何變換的 效果，對(duì)其作出分析；(3) 編寫圖像的位置變換(平移、旋轉(zhuǎn)、鏡像)，圖像的形狀變換(縮小、放大、錯(cuò)切)；等幾何變換的算法原理及編程實(shí)現(xiàn)；(4) 掌握matlab編程環(huán)境中基本的圖像處理函數(shù)；(5) 整理實(shí)驗(yàn)結(jié)果，撰寫科研報(bào)告。二. 實(shí)驗(yàn)原理(1) 圖像的幾何變換，是指原始圖像按照需要產(chǎn)生大小、形狀和位

2、置的變化。從圖像類型來(lái)分，圖像的幾何變換可以有二維平面 圖像的幾何變換和三維圖像的幾何變換以及由三維向二維平面 投影的變換等。從變換的性質(zhì)分，圖像的兒何變換有平移、比 例縮放、旋轉(zhuǎn)和錯(cuò)切等基本變換。(2) j=imreadc12.pngf),讀取圖片信息;subplot(a,b,c),分配一個(gè) a 行b列的圖像窗口，將現(xiàn)在的圖片放在第c個(gè)子模塊中； imshow(j),顯示圖像j；(3) 圖像縮放 imresize 函數(shù)，語(yǔ)法格式為：b = imresize(a,m,method) 這里參數(shù)method用于指定插值的方法，口j選用的值為nearest1(最鄰近法'bilinear1 (

3、雙線性插值),'bicubic1 (雙三次插值), 默認(rèn)為nearest1 o 即 b 二 imresize(a,m,method)返回原圖 a 的 m 倍 放人的圖像(當(dāng)m小丁 1時(shí)效果是縮小原圖的m倍)。三. 實(shí)驗(yàn)代碼及結(jié)果1平移：將圖像向左下方和右下方水平移動(dòng)，顯示結(jié)果在一個(gè)圖像窗口;j=imread('c:usersthinkdesktop12.png,);figure;subplot( 1,3,1);imshow(j);titlec 原圖 j;s=translate(strel(l),5o,-5o);jl=imdilate(j,s);subplot(l,3,2);im

4、show(jl);title(左下方平移j;s=translate(strel(l),50,50);j2=imdilate(j,s);subplot(l,3,3);imshow(j2);title。右下方平移j;(1)2縮小放大：將圖像分別放大2倍和縮小0.5倍，插值方法使用雙線性插值法，分別顯示圖像;i=imread('blobs.png,);figure;imshow(i);title(源圖 j;j l=imresize(1,0.5,'nearest1);figure;imshow(jl);title(縮小圖j;j2=imresize(1,2,'nearest

5、9;); figure;imshow(j2);title c放大圖)2鏡像：將圖像水平鏡像，垂直鏡像，顯示結(jié)果在一個(gè)圖像窗口;imreadcblobs.png1);figure;subplot (1,3 j)；imshow(i);title。原圖 j;tform 1 =maketform(,affine,-1 0 0;0 1 0;l 0 1);j1 =imtransform(i,tform 1 /nearest1);subplot(l,3,2);imshow(jl); title。水平鏡像j;tform2=maketform(,affine, 1 0 0;0 -1 0;0 1 1);j2=im

6、transform(i,tform2,'nearest1);subplot( 1,3,3);imshow(j2);title(唾直鏡像)；原像水平覆像垂直鏡像4錯(cuò)切：將圖像分別進(jìn)行垂直45度錯(cuò)切，顯示在同一個(gè)圖像窗口; f=imread(,c:usersthinkdesktop12.png,);subplot( 1,2,1);imshow(f)；title(原圖 j;h=size(f);fl=zeros(h(l)+round(h(2)*tan(pi/6),h(2),h(3);for m=l:h(l)for n=l:h(2)fl (m4-round(n*tan(pi/6),n, 1:h(

7、3)=f(m,n,l:h(3);endendsubplot(l,2,2);imshow(uint8(fl);title(唾直 45 度);原圖5旋轉(zhuǎn)：將圖像分別進(jìn)行順時(shí)針30度旋轉(zhuǎn)，圖像顯示在一個(gè)圖像窗口;imreadcblobs.png1);figure;subplot( 1,2,1);imshow(i);title。原圖 j;m,n=size(i);ang=30;j=i；for i=l:mfor j=l:nx=floor(i-m/2)*cos(ang*pi/180)-(j-n/2)*(-sin(ang*pi/180)4-m/2);y=floor(j-n/2)*cos(ang*pi/180)

8、+(i-m/2)*(-sin(ang*pi/180)+n/2);if(xvm)&(y<n)&(x>0)&(y>0)j(i,j)=i(x,y);elsej(ij)=o;endendendsubplot(l,2,2);title。圖像旋轉(zhuǎn)j;四.作業(yè)1. 將一副50*50的數(shù)字圖像逆吋針旋轉(zhuǎn)30度，顯示旋轉(zhuǎn)后圖像 并指出旋轉(zhuǎn)后圖像的大小。將該圖像平移10習(xí)，顯示平移后 圖像，要求將畫布擴(kuò)大。將該平移后圖像再逆時(shí)針旋轉(zhuǎn)15度。i=imread(,blobs.png,);figure;imshow(i);title(原圖 j;figure;j=imcrop(i

9、,l00,100,50,50);imshow(j);title('錯(cuò)切后)；figure;h=imrotate(j,30,'bilinear1);title(30度旋轉(zhuǎn)后);figure;s=translate(strel( 1), 10,5);w=imdilate(i,s);imshow(w);figure;r=imrotate(w, 15/bilinear1);imshow(r);titlefl5度旋轉(zhuǎn)后*);30度旋轉(zhuǎn)后2. 將一副圖像先旋轉(zhuǎn)一定角度然后再平移和將該圖像先平移再 選轉(zhuǎn)是否一樣？給出實(shí)驗(yàn)結(jié)果。(1).imreadc'blobs.png1);figur

10、e;subplot( 1,3,1);imshow(i);title。原圖 j;s=translate(strel(l),5o,5o);w=imdilate(i,s);subplot(l,3,2);imshow(w);titled平移后j;r=imrotate( w, 30, 'bilinear');subplot(l,3,3);titlec30度旋轉(zhuǎn)后')；(2).i=imread(,blobs.png,);figure;subplot( 1,3,1);imshow(i);title(隙圖)；r=imrotate(1,30,bilinear');subplot( 1,3,2);imshow(r);title(*30度旋轉(zhuǎn)后j;s=translate(strel( 1),50,50);w=imdilate(r,s);subplot(l,3,3);imsh