matlab圖像幾何變換

1、最新資料推 姑斤希 段工學(xué)卩免學(xué) 光電圖像處理實驗報告 實驗一：圖像幾何變換 作 者肖劍洪 專 業(yè)電弄I學(xué)皺術(shù) 學(xué) 院物理與電子學(xué)院指導(dǎo)老師 王曉明 完成時間 2013.12.2 實驗一圖像幾何變換 一、實驗?zāi)康?圖像兒何變換是圖像處理與分析的基礎(chǔ)內(nèi)容之一。使學(xué)生掌握圖像平 移、 鏡像、縮放和旋轉(zhuǎn)等兒何變換的方法以及程序?qū)崿F(xiàn)。 二、實驗設(shè)備 聯(lián)想圖像處理工作站 三、實驗內(nèi)容及要求 圖像兒何變換不改變圖像的像素值，而是改變像素所在的幾何位置，具有獨 特的方法。本實驗的內(nèi)容如下： 1. 學(xué)會使用Mat lab的圖像處理工具箱仃mage Processing Toolbox

2、) o使學(xué) 生初步具備使用該軟件處理圖像信息能力，并可以利用該軟件完成本課程規(guī)定 的其他實驗和作業(yè)。 式取整。 圖像的比例縮放所產(chǎn)生的圖像中的像素可能在源圖像中找不到相應(yīng)的像素 點，這時就需要進(jìn)行插值處理，通?？捎米罱彶逯岛碗p線性插值算法來實現(xiàn)。 (1) 最近鄰插值 對通過反變換得到的浮點坐標(biāo)進(jìn)行簡單的取整操作，得到一個整數(shù)型坐標(biāo)， 這個整數(shù)型坐標(biāo)對應(yīng)的原圖像的像素值就是新圖像的像素值。 (2) 雙線性插值 對通過反變換得到的浮點坐標(biāo)，求出此浮點坐標(biāo)對應(yīng)原圖像四個鄰域整數(shù)型 坐標(biāo)的灰度值，然后經(jīng)過雙線性插值得到新圖像的像素值。 (3) 雙三次插值 對通過反變換得到的浮點坐標(biāo)，求出此浮點坐標(biāo)對

3、應(yīng)原圖像16個鄰域整數(shù) 型坐標(biāo)的灰度值，然后經(jīng)過雙三次插值得到新圖像的像素值。 4.圖像旋轉(zhuǎn) 一般圖像的旋轉(zhuǎn)是以圖像的中心為原點，將圖像上的所有像素都旋轉(zhuǎn) 個相同的角度。如圖所示： 新圖像坐標(biāo)與原圖像坐標(biāo)對應(yīng)關(guān)系: 兩個常量: 圖像旋轉(zhuǎn)之后會出現(xiàn)許多空洞點，須對這些空洞點進(jìn)行插值處理?？捎媒彶?值，均值插值算法。 五、實驗步驟 1. 根據(jù)實驗內(nèi)容的要求在Matlab軟件中編寫相應(yīng)程序； 2. 調(diào)試運行程序，并記錄結(jié)果： 3. 完成實驗報告。 六、實驗程序清單和實驗結(jié)果及分析 1.圖像平移 I 二 imread(f:/l .jpg)； se = translate(strel(l)z 25 2

4、5); J = imdila te(l,se); imshow(l), title(Originar) figure, imshow(J)z title(Translated); UriqinalTranslated 2.圖像鏡像 A=imread(f:/1 .jpg)； height,width,dim二 size( A); tformmaketformfaffineJ-l 0 0;0 1 0;width 0 1); B=imtransform (Aztform/neoest); tform2=maketform(affine,l 0 0;0 -1 0;0 heigh11); C=imtr

5、on sform (Aztform2/neorest); subplotfl ,3,1 )Jmshow(A); title(原圖像)； subplot(l ,3,2),imshow(B); title fTk 平鏡像)； subplot(l ,3,3),imshow(C); title (*!g 直鏡像)； 原團像 3.圖像比例縮放 水平鏡像 豎直鏡像 T1 二 maketform(affine,2 0 0; 0 1 0; 0 0 1); T2 二 maketformfaffineJl 0 0; 0 2 0; 0 0 1); T3 = maketform(affine,4/5 0 0; 0 3

6、/5 0; 0 0 1); I 二 imread(f:/l .jpg); II = imtr on sform(IJl); 12 = imtr on sform(IJ2); 13 = imtransform(IJ3); imshowfl), figure, imshow(ll)z figure, imshow(l2), figure, imshow(l3); 4. 圖像旋轉(zhuǎn) lmage=imread(,f:/1 .jpg1)； angle=30; pai=3.14; iml=(lmage);subplot(211); imshow(iml); titlefr); 最新資料推 示心3_K3_ r

7、oteang=angle/180*pai; imlXziml Y=size(iml); rotmatrix=cos(rote on g)z-sin(roteang);si n(oteong),cos(oteo ng)； rrotmatrix= cos(-roteang)z-sin(-roteang):sin (-rotea ng) ,cos (-ot eang); im2XY 1 = im 1 Xzim 1Y *rotmatrix; im2XY2= imlXzl*rotmatrix; im2XY3=l z 1 *rotmatrix; im2XY4= Liml Y*rotmatrix; im2p

8、l l=min(im2XYl (1 )zim2XY2(l )Jm2XY3(l)jm2XY4(l )zmin (im2XYl(2)zim2XY2(2)zim2XY3(2)zim2XY4(2)J; im2X=round(max(im2XYl (1 ),im2XY2(l) Jm2XY3(l) Jm2XY4(l) -im2pll(l); im2Y=round(max(im2XYl(2)jm2XY2(2)/im2XY3(2)zim2XY4(2) -im2pll(2); im2=uint8(zeros(im2Xjm2Y)+255; for i=l:im2X forj=l:im2 Y pij=round(i+im2pl 1 (l)-l/j+im2pl 1 (2)-l*rrotmatrix); if (pij(1)二 1 ) end end end; subplot(212); imshow(im2);