第六章 matlab在數(shù)字圖像處理中的應(yīng)用

第五章matlab在數(shù)字圖像處理中的應(yīng)用數(shù)字圖像的讀入A=imread(filename,fmt):將文件名為filename,擴(kuò)展名為fmt(圖像文件格式)表示的圖像文件中的數(shù)據(jù)讀到矩陣A中；若filename所指為灰度圖像，則A為二維矩陣，若filename為真彩色圖像，則A為m×n×3的三維矩陣。A=imread(filename,fmt):字符串fmt的可能值有：“jpg”or“jpeg”“tif”or“tiff”“gif”“bmp”“png”“hdf”“pcx”“xwd”“cur”“ico”圖像文件的讀入圖像文件的寫入：imwrite(A，filename,fmt):將表示圖像文件數(shù)據(jù)的矩陣A寫入到文件名為filename,擴(kuò)展名為fmt(圖像文件格式)表示的圖像文件中；A、filename、fmt意義與inread命令相同。圖像類型轉(zhuǎn)換在matlab系統(tǒng)下，要將一個類型的圖像文件轉(zhuǎn)換成另一個類型的圖像文件也是一件輕而易舉的事情，只要將前一個文件的圖像數(shù)據(jù)讀出，再寫到后一個圖像文件中去即可。bitmap=imread('mybitmap.bmp','bmp');imwrite(bitmap,'mybitmap.png','png');圖像的顯示：在matlab系統(tǒng)下，最常用的顯示圖像的函數(shù)為imshow函數(shù)，可用于顯示各類圖像，它的調(diào)用方法分別如下：imshow(X,map):顯示索引圖像，X為索引圖像矩陣，map為色彩圖示；imshow（I）：顯示灰度圖像，I為二值圖像矩陣；imshow（BW）：顯示二值圖像，BW為黑白二值圖像矩陣；imshow（RGB）：顯示RGB圖像，RGB為RGB圖像矩陣；imshowfilename:顯示圖像文件，如imshow

rice.pngI=imread(‘rice’,’png’);%圖像的讀入imshow(I);％圖像的顯示imshow

rice.png在不同顯示板顯示多個圖像imshow(I(:,:,:,1))figure;imshow(I(:,:,:,2))figure;imshow(I(:,:,:,3))在一個顯示板顯示多個圖像[X1,map1]=imread('forest.tif');[X2,map2]=imread('trees.tif');subplot(1,2,1),imshow(X1,map1)subplot(1,2,2),imshow(X2,map2)用相同的colormap導(dǎo)致圖像不清晰不建議用imshow函數(shù)顯示多個圖像[X1,map1]=imread('forest.tif');[X2,map2]=imread('trees.tif');subplot(1,2,1),subimage(X1,map1)subplot(1,2,2),subimage(X2,map2)在一個顯示板顯示多個圖像圖像的一些特殊顯示手段：對圖像加顏色條：對于了解被顯示圖像的灰度級別特別有用。I=imread(‘rice’,’png’);imshow(I);colorbar;一次顯示多層圖像的所有層mri=uint8(zeros(128,128,1,27));forframe=1:27[mri(:,:,:,frame),map]=imread('mri.tif',frame);endmontage(mri,map);動畫顯示多層圖像mri=uint8(zeros(128,128,1,27));forframe=1:27[mri(:,:,:,frame),map]=imread('mri.tif',frame);endmov=immovie(mri,map);%構(gòu)建動畫文件movie(mov);%播放動畫紋理成圖功能前面提到的圖像顯示手段都只能使圖像在二維平面上顯示，matlab的一個強(qiáng)大功能是能將平面圖像顯示在空間三維曲面上，由warp函數(shù)來實(shí)現(xiàn)。I=imread(‘flowers’,’jpg’);[x,y,z]=cylinder;subplot(1,2,1);warp(x,y,z,I);[x,y,z]=sphere(50);subplot(1,2,2);warp(x,y,z,I)圖像的空間轉(zhuǎn)換改變圖像的尺寸：B=imresize(A,m,method):返回一個m倍于原圖像A的圖像B，使用的差值方法由method指定。B=imresize(A,[mrows,mcols],method):返回一個尺寸大小為mrows行乘以mcols列的圖像B。B=imresize(I,2);Imshow(B);圖像的旋轉(zhuǎn)B=imrotate(A,angle,method):angle為正，則逆時針旋轉(zhuǎn)，為負(fù)，則順時針旋轉(zhuǎn)。圖像B由圖像A旋轉(zhuǎn)angle角度得到，通常大于圖像AB=imrotate(A,angle,method,’crop’):圖像B由圖像A旋轉(zhuǎn)angle角度得到，尺寸與A相等。原圖B=imrotate(I,45);B=imrotate(I,45,’crop’);B=imcrop(I,rect):rect為可選參數(shù)，指定選取的圖像區(qū)域，其數(shù)據(jù)格式為：[xmin,ymin,width,height];也可使用不帶參數(shù)的imcrop函數(shù)，則等待用戶利用鼠標(biāo)在待裁減圖像上進(jìn)行選擇。圖像的裁剪B=crop(I,[40,40,80,80]);Imshow(B);imshowcircuit.tifI=imcrop;imshow(I);圖像分析與提高圖像直方圖及密度均衡；圖像的噪聲去除。圖像的邊緣檢測；圖像直方圖定義：如果將圖像中像素亮度（灰度級別）看成是一個隨機(jī)變量，則其分布情況就反映了圖像的統(tǒng)計特性，可用PDF(ProbabilityDensityFunction)來刻畫和描述，表現(xiàn)為灰度直方圖。灰度直方圖是灰度級的函數(shù)，它表示圖像中具有某種灰度級的像素的個數(shù)，反映了圖像中每種灰度出現(xiàn)的頻率?；叶戎狈綀D的橫坐標(biāo)是灰度級，縱坐標(biāo)是該灰度出現(xiàn)的頻率?；叶戎狈綀D的計算編程實(shí)現(xiàn)灰度直方圖的計算與繪制I=imread(‘rice.png’);[a,b]=size(I);H=zeros(1,256);fori=1:aforj=1:bfork=0:255ifI(i,j)==kH(k+1)=H(k+1)+1;endendendendplot(H);利用matlab函數(shù)繪制直方圖I=imread(‘rice.png’);imhist(I);%繪制直方圖直方圖均衡當(dāng)某些圖像灰度分布集中于某一個區(qū)域時，圖像會出現(xiàn)過亮或過暗的情況，在這種情況下可以利用直方圖均衡的方法使得圖像灰度值均勻分布在全部灰度值范圍內(nèi)，或者使其灰度值分布與一已知的灰度分布一致。J=histeq(I,hgram):將圖像I進(jìn)行灰度轉(zhuǎn)換，使得轉(zhuǎn)換后的圖像J的灰度直方圖具有hgram的分布。J=histeq(I,n):將圖像I進(jìn)行灰度轉(zhuǎn)換，使得轉(zhuǎn)換后的圖像J的灰度分布具有n個灰度級別，且每一個灰度級別內(nèi)的像素數(shù)相等。I=imread(pout.tif);subplot(2,2,1);imshow(I);subplot(2,2,2);imhist