貴州大學(xué)實(shí)驗(yàn)報(bào)告學(xué)院：專(zhuān)業(yè)：班級(jí)：姓名學(xué)號(hào)實(shí)驗(yàn)組實(shí)驗(yàn)時(shí)間指導(dǎo)教師成績(jī)實(shí)驗(yàn)項(xiàng)目名稱(chēng)實(shí)驗(yàn)一：圖像的基本運(yùn)算及常用變換實(shí)驗(yàn)?zāi)康耐ㄟ^(guò)本實(shí)驗(yàn)的學(xué)習(xí)使學(xué)生熟悉MATLAB圖像處理工具箱的編程和使用；掌握數(shù)字圖像的基本算術(shù)運(yùn)算和邏輯運(yùn)算的原理和應(yīng)用；掌握常用的變換：二維離散傅立葉變換、二維離散余弦變換的正反變換，及一些性質(zhì)的驗(yàn)證。實(shí)驗(yàn)要求集中授課的教學(xué)形式實(shí)驗(yàn)原理1、圖像文件的讀取：imread函數(shù)語(yǔ)法：A=imread(‘filename’,fmt)[X,map]=imread(‘filename’,fmt)[...]=imread(‘filename’)[...]=imread(...,idx)(CUR,ICO,andTIFFonly)[...]=imread(...,ref)(HDFonly)[...]=imread(...,'BackgroundColor',BG)(PNGonly)[A,map,alpha]=imread(...)(PNGonly)2、圖像文件的顯示：imshow函數(shù)當(dāng)用戶(hù)調(diào)用imshow函數(shù)顯示一幅圖像時(shí)，該函數(shù)將自動(dòng)設(shè)置圖像窗口、坐標(biāo)軸和圖像屬性。這些自動(dòng)設(shè)置的屬性具體包括圖像對(duì)象的CData屬性和CDataMapping屬性、坐標(biāo)軸對(duì)象的CLim屬性、圖像窗口對(duì)象的Colormap屬性。另外，根據(jù)用戶(hù)使用參數(shù)的不同，imshow函數(shù)在調(diào)用時(shí)除了完成前面提到的屬性設(shè)置外還可以完成以下的操作：設(shè)置其它的圖形窗口對(duì)象的屬性和坐標(biāo)軸對(duì)象的屬性以?xún)?yōu)化顯示效果，如可以設(shè)置隱藏坐標(biāo)軸及其標(biāo)示等。包含和隱藏圖像邊框。調(diào)用truesize函數(shù)以顯示沒(méi)有彩色漸變效果的圖像。語(yǔ)法：imshow(I,n)imshow(I,[lowhigh])imshow(BW)imshow(X,map)imshow(RGB)imshow(...,display_option)imshow(x,y,A,...)imshowfilenameh=imshow(...)3、在圖像處理中一種最簡(jiǎn)單的圖像恢復(fù)算法是利用其鄰域像素的值求和再取平均來(lái)得到改點(diǎn)的一個(gè)恢復(fù)值：4、運(yùn)動(dòng)檢測(cè)原理：在序列圖像中，通過(guò)逐像素比較可直接求取前后兩幀圖像之間的差別。假設(shè)照明條件在多幀圖像間基本不變化，那么差圖像中的不為零處表明該處的像素發(fā)生了移動(dòng)。換句話說(shuō)，對(duì)時(shí)間上相鄰的兩幅圖像求差可以將圖像中目標(biāo)的位置和形狀變化突出出來(lái)。5、Matlab中的fft2函數(shù)語(yǔ)法：Y=fft2(X)Y=fft2(X,m,n)6、Matlab中的ifft2函數(shù)語(yǔ)法：Y=ifft2(X)Y=ifft2(X,m,n)7、二維傅立葉變換的平移性質(zhì)。8、Matlab中二維離散余弦變換函數(shù)dct2語(yǔ)法：B=dct2(A)B=dct2(A,m,n)B=dct2(A,[mn])9、Matlab中的函數(shù)dctmtx此函數(shù)可以用來(lái)求離散余弦變換的變換矩陣D。語(yǔ)法：D=dctmtx(n)10、DCT域的圖像壓縮原理。離散余弦變換DCT（DiscreteCosineTransform）是數(shù)碼率壓縮的一個(gè)常用的變換編碼方法。DCT是先將整體圖像分成N×N像素塊，然后對(duì)N×N像素塊逐一進(jìn)行DCT變換。在變換編碼中一個(gè)很重要的因素是塊的大小。最通用的大小是8×8和16×16，由于計(jì)算方面的原因，兩個(gè)都是2的乘方。若是采用8×8的變換塊，則得到的變換為稱(chēng)為DCT系數(shù)。此矩陣的左上角系數(shù)相當(dāng)于最低頻率的系數(shù)，簡(jiǎn)稱(chēng)為子圖像的直流（DC）系數(shù)或直流成分；除此之外的均為交流系數(shù)。在右下角的系數(shù)是最高頻率的系數(shù)；在右上角的系數(shù)，表示水平方向頻率最高，垂直方向頻率最低的系數(shù)；在左下角的系數(shù)，表示垂直方向頻率最高，水平方向頻率最低的系數(shù)。另外在進(jìn)行DCT變換以前，二維空間圖像亮度數(shù)據(jù)通常較高，應(yīng)該先進(jìn)行向下電平的移位。如果是256灰度級(jí)，則將亮度值剪掉128。由于大多數(shù)圖像的高頻分量較少，相應(yīng)于圖像高頻分量的系數(shù)經(jīng)常為零，加之人眼對(duì)高頻成分的失真不太敏感，所以可用更粗的量化。因此，傳送變換系數(shù)的數(shù)碼率要大大小于傳送圖像像素所用的數(shù)碼率。到達(dá)接收端后通過(guò)反離散余弦變換回到原值。雖然會(huì)有一定的失真，但人眼是可以接受的。例：有一幅灰度圖像為：電平移位后為：DCT系數(shù)從上面的DCT變換系數(shù)矩陣可以看出高頻分量占得很少。所以對(duì)于系數(shù)小的我們可以分配較少的比特，大的系數(shù)可以分配較多的比特，對(duì)其進(jìn)行量化編碼。實(shí)驗(yàn)儀器計(jì)算機(jī)一臺(tái)；Matlab軟件實(shí)驗(yàn)步驟讀取一個(gè)彩色圖像文件并顯示。將此圖像轉(zhuǎn)換為256色的灰度圖像。用subimage函數(shù)將彩色圖像和灰度圖像同時(shí)顯示。讀取灰度圖像。用imnosie函數(shù)加入0.02的椒鹽噪聲；采用8鄰域像素平均的方法對(duì)有噪聲的圖像去噪。同時(shí)顯示原圖像，噪聲圖像，和恢復(fù)后的圖像。加入不同參數(shù)的椒鹽噪聲比較實(shí)驗(yàn)結(jié)果。編寫(xiě)程序。產(chǎn)生一個(gè)50*50的矩陣。其中第（20,10）到（20,15）的值為255。其余矩陣元素為0。對(duì)該矩陣做二維傅立葉變換，并顯示其頻譜圖。將該圖像中的目標(biāo)向y方向平移。即將（20,35）到（20,40）的元素置為255。其余值為0。對(duì)該矩陣做二維傅立葉變換，并顯示其頻譜圖。對(duì)兩頻譜圖進(jìn)行比較，得出結(jié)論。依照以上方法驗(yàn)證傅立葉變換的X方向的平移性質(zhì)。對(duì)給定圖像進(jìn)行二維DCT變換，并將DCT系數(shù)矩陣中值小于給定閾值的元素置為0，然后使用反離散余弦變換重建圖像。同時(shí)顯示原圖像和重建圖像。實(shí)驗(yàn)內(nèi)容1、彩色圖像文件的讀取和顯示。2、利用8鄰域平均法對(duì)圖像去噪。3、編寫(xiě)程序。讀取圖像manface1.jpg和manface2.jpg。用兩圖像相減的方式獲取運(yùn)動(dòng)軌跡4、二維離散傅立葉變換的平移性質(zhì)驗(yàn)證。5、二維離散余弦變換的應(yīng)用。實(shí)驗(yàn)數(shù)據(jù)步驟1代碼部分clearall;closeall;clc;A=imread('Lena.bmp');%讀取彩色圖像figure(1)%產(chǎn)生新的圖形窗口，用于顯示imshow(A)%顯示圖像title('Lena彩色圖像顯示')%給圖像加標(biāo)題B=rgb2gray(A);%彩色圖轉(zhuǎn)灰度圖figure(2)subplot(121),subimage(A)%分為1×2兩個(gè)窗口，當(dāng)前選擇在第一個(gè)窗口subplot(122),subimage(B)gtext('Lena彩色圖像和灰度圖像顯示')步驟2代碼部分clearall;closeall;clc;A=imread('Cameraman.bmp');%給圖片加噪聲密度為0.02的椒鹽噪聲B=imnoise(A,'salt&pepper',0.02);%制定卷積核h=ones(3,3)/5;%ones產(chǎn)生全部元素都為1的矩陣h(1,1)=0;%矩陣特定元素賦值h(1,3)=0;h(3,1)=0;h(3,3)=0;%平滑運(yùn)算C=imfilter(B,h);subplot(131),imshow(A)subplot(132),imshow(B)subplot(133),imshow(C)噪聲密度為0.02噪聲密度為0.06步驟3代碼部分clearall;closeall;clc;A1=zeros(50,50);%zeros產(chǎn)生全部元素都為0的矩陣A1(20,10:15)=255;%二維傅立葉變換F1=fft2(A1);F11=abs(F1);A2=zeros(50,50);A2(20,35:45)=255;F2=fft2(A2);F21=abs(F2);figure(1)subplot(121),imshow(F11),title('平移前頻譜');subplot(122),imshow(F21),title('平移后頻譜');colormap(jet);colorbar;%色圖設(shè)置步驟4代碼部分clearall;closeall;clc;A=imread('Lena.bmp');B=rgb2gray(A);subplot(221),imshow(B),title('原圖像');%對(duì)給定圖像進(jìn)行二維DCT變換Y=dct2(double(B)/255);subplot(222),imshow(mat2gray(Y)),title('DCT圖像');%轉(zhuǎn)換矩陣為灰度圖像fori=1:185forj=1:255if(abs(Y(i,j)<0.05))Y(i,j)=0;%并將DCT系數(shù)矩陣中值小于給定閾值的元素置為0endendendsubplot(22