數(shù)字圖像處理

1、精選優(yōu)質(zhì)文檔-傾情為你奉上內(nèi)蒙古科技大學(xué)本科生課程設(shè)計(jì)論文 數(shù)字圖像處理課設(shè) 學(xué)生姓名：于鵬飛 學(xué) 號(hào)： 專 業(yè)：電子信息工程 班 級(jí)：13信息二班 指導(dǎo)教師：張明 2016年 5 月 實(shí)驗(yàn)一、計(jì)算圖象統(tǒng)計(jì)參數(shù)及直方圖變換編寫程序：1.讀取圖像（文件名為cameraman.tif）;2.最大值、最小值、均值3.直方圖（要求不能調(diào)用imhist函數(shù)，只能用該函數(shù)來驗(yàn)證自編直方圖程序塊的正確性）。4.計(jì)算并分析圖像直方圖（imhist）5.利用直接灰度變換法對(duì)圖像進(jìn)行灰度變換（imadjust）6.直方圖均衡化處理及顯示（histeq）程序f=imread('C:UsersYSSADesk

2、top數(shù)字圖像課設(shè)于鵬飛圖片cameraman.tif');%求最大值、最小值、均值max(f(:) %最大值min(f(:) %最小值mean(f(:) %均值%顯示其直方圖m,n=size(f); %圖像大小X=zeros(256,1); %灰度級(jí)范圍for i=1:m %對(duì)應(yīng)的統(tǒng)計(jì)位加1 for j=1:n x=f(i,j); X(x)=X(x)+1; endend;%計(jì)算并分析圖像直方圖（imhist）g=imhist(f,256); %利用直接灰度變換法對(duì)圖像進(jìn)行灰度變換（imadjust）g2=imadjust(f,0 1,1 0);%對(duì)原圖像進(jìn)行直方圖均衡化處理g3=hi

3、steq(f); subplot(2,3,1),imshow(f);title('原始圖像');subplot(2,3,2);bar(X); title('原圖像直方圖') subplot(2,3,3);bar(g); title('驗(yàn)證原圖像直方圖') subplot(2,3,4),imshow(g2);title('灰度變換');subplot(2,3,6),imshow(g3);title('直方圖均衡化');照片實(shí)驗(yàn)二傅里葉變換及空間域和頻率域?yàn)V波1、利用 MATLAB 實(shí)現(xiàn)數(shù)字圖像的傅里葉變換讀入并顯示圖

4、“test1.tif”，作該圖的二維 FFT 變換 F，將其直流分量移到頻譜中心 F1（fftshift()函數(shù)），計(jì)算其實(shí)部 RR、虛部 II，用兩種方法計(jì)算幅值 A1=abs(F1)和 A2=sqrt(RR.2+II.2)，分別顯示 A1 和 A2，并加以比較。2、近似沖擊函數(shù)二維傅里葉變換A=zeros(99,99); A(49:51,49:51)=1;作A的二維傅里葉變換 B，將B直流分量移到頻譜中心B1，分別用函數(shù)imshow和mesh顯示A和B1模的對(duì)數(shù)（log(1+abs(B1)）。程序g=imread(' C:UsersYSSADesktop數(shù)字圖像課設(shè)于鵬飛圖片tes

5、t1.tif');%對(duì)圖像進(jìn)行傅立葉變換G=fft2(g);%將其直流分量移到頻譜中心F=fftshift(G);%實(shí)部 RR、虛部 IIRR=real(F); II=imag(F);%顯示模 A1 和 A2，并加以比較A1=abs(F);A2=sqrt(RR.2+II2);aa1=log(1+A1);aa2=log(1+A2);% 生成一個(gè)99*99的零矩陣A=zeros(99,99); A(49:51,49:51)=1;%B=fft2(A);%fftshift就是對(duì)換數(shù)據(jù)的左右兩邊B1=fftshift(B);%顯示A和B1模的對(duì)數(shù)（log(1+abs(B1)）figure;sub

6、plot(2,4,1),imshow(g);title('原圖像'); subplot(2,4,2),imshow(G);title('傅立葉變換');subplot(2,4,3),imshow(aa1,);title('a1 abs');subplot(2,4,4),imshow(aa2,);title('a2 sqrt');subplot(2,4,5),imshow(A);title('沖擊函數(shù)imsow');subplot(2,4,6),mesh(A);title('沖擊函數(shù)mesh');su

7、bplot(2,4,7),imshow(log(1+abs(B1);title('沖擊函數(shù)頻域模mesh');subplot(2,4,8),mesh(log(1+abs(B1);title('沖擊函數(shù)頻域模mesh');實(shí)驗(yàn)2_1-2圖片3、空間濾波與頻域?yàn)V波將圖“test2.tif”分別進(jìn)行空間與頻域?yàn)V波?？臻g濾波：用 fspecial 產(chǎn)生 9×9，標(biāo)準(zhǔn)差為 2 的高斯濾波器 w，用函數(shù) imfilter 對(duì) f 進(jìn)行空間濾波，得到濾波后的圖像 fi1。頻域?yàn)V波:將上述高斯濾波器 w 用函數(shù) freqz2 獲得其256×256 頻域形式

8、W，用 W 在頻域?yàn)V波圖像 f（注意 W 的直流分量在頻譜中心），得到 Ff，求其傅里葉逆變換（ifft2），得到濾波后的圖像 fi2。用函數(shù) imshow 顯示 w，用函數(shù) mesh 顯示 W，用 mshow 顯示 f，fi1，F(xiàn)f 的模的對(duì)數(shù)，fi2。比較 fi1 與 fi2。程序：f=imread(' C:UsersYSSADesktop數(shù)字圖像課設(shè)于鵬飛圖片test2.tif');%傅立葉變化G=fft2(f);%用 fspecial 產(chǎn)生 9*9，標(biāo)準(zhǔn)差為 2 的高斯濾波器 ww=fspecial('gaussian',9,2);%用函數(shù) imfilt

9、er 對(duì) f 進(jìn)行空間濾波fi1=imfilter(f,w);%高斯濾波器 w 用函數(shù) freqz2 獲得其256*256 頻域形式 WW=freqz2(w,256,256);%ifftshift長(zhǎng)度不是偶數(shù)，對(duì)換數(shù)據(jù)的左右兩邊W1=ifftshift(W);Ff=G.*W1;%傅里葉逆變換（ifft2）的實(shí)部fi2=real(ifft2(Ff); % imshow 顯示 w, mesh 顯示 Wsubplot(2,3,1),imshow(f);title('原圖像');subplot(2,3,2),imshow(G);title('傅立葉變換');subplo

10、t(2,3,3),imshow(w);title('w');subplot(2,3,4),mesh(W);title('W');subplot(2,3,5),imshow(fi1,);title('空間濾波');subplot(2,3,6),imshow(fi2,);title('頻域?yàn)V波')實(shí)驗(yàn)2_3圖片實(shí)驗(yàn)三彩色圖像處理3_1.使用 RGB、Indexed和Gray圖像間轉(zhuǎn)換函數(shù)對(duì)圖像“strawberries.tif”使用rgb2ind 分別產(chǎn)生8色抖動(dòng)和非抖動(dòng)圖像；使用rgb2gray實(shí)現(xiàn)圖像轉(zhuǎn)換，并使用函數(shù)dither產(chǎn)

11、生其抖動(dòng)形式的圖像。程序f=imread(' C:UsersYSSADesktop數(shù)字圖像課設(shè)于鵬飛圖片strawberries.tif');%rgb2ind 將真彩色圖像轉(zhuǎn)換為索引圖像，產(chǎn)生8 色非抖動(dòng)圖像（彩色映射矩陣）x1,map1=rgb2ind(f,8,'nodither');%產(chǎn)生8 色抖動(dòng)圖像x2,map1=rgb2ind(f,8,'dither');%使用 rgb2gray 實(shí)現(xiàn)RGB到灰度圖像g=rgb2gray(f);%函數(shù)dither產(chǎn)生其抖動(dòng)形式的圖像g1=dither(g);subplot(1,2,1),imshow(x

12、1,map1);title('ind nodither');subplot(1,2,2),imshow(x2,map1);title('ind dither');figure,subplot(1,2,1),imshow(g);title('gray nodither');subplot(1,2,2),imshow(g1);title('gray dither');實(shí)驗(yàn)3_1圖片3_2.彩色圖像平滑與銳化對(duì)圖像“flower.tif”在RGB空間實(shí)現(xiàn)彩色圖像的平滑（w=ones（25）./(25*25)）與銳化（w= 1 1 1 1

13、 1; 1 1 1 1 1;1 1 -24 1 1; 1 1 1 1 1; 1 1 1 1 1;）（函數(shù) imfilter）;在HSI空間實(shí)現(xiàn)彩色圖像的平滑（w=ones（25）./(25*25)），（函數(shù)imfilter, rgb2hsi, cat），觀察僅平滑亮度分量（intensity）和平滑全部三個(gè)分量結(jié)果的差別。程序f=imread(' C:UsersYSSADesktop數(shù)字圖像課設(shè)于鵬飛圖片flower.tif ');%在RGB空間實(shí)現(xiàn)彩色圖像的平滑w1=ones(25)./(25*25);%實(shí)現(xiàn)彩色圖像的銳化w2=1 1 1 1 1;1 1 1 1 1;1 1

14、-24 1 1;1 1 1 1 1; 1 1 1 1 1;i1=imfilter(f,w1,'replicate');i2=imfilter(f,w2,'replicate'); %在HSI空間實(shí)現(xiàn)彩色圖像的平滑（w=ones（25）./(25*25)） 函數(shù)imfilter, rgb2hsi, cath=rgb2hsi(f);w1=ones(25)./(25*25);H=h(:,:,1);S=h(:,:,2);I=h(:,:,3);h1=imfilter(h,w1,'replicate');%平滑全部三個(gè)分量I2=imfilter(I,w1,&

15、#39;replicate');%僅平滑亮度分量h2=cat(3,H,S,I2);hi1=hsi2rgb(h1);hi2=hsi2rgb(h2);hi1=min(hi1,1);hi2=min(hi2,1);subplot(2,3,1),imshow(f);title('原圖像');subplot(2,3,2),imshow(i1);title('RGB平滑');subplot(2,3,3),imshow(i2);title('RGB銳化');subplot(2,3,4),imshow(f);title('原圖像');sub

16、plot(2,3,5),imshow(hi1);title('平滑全部三個(gè)分量')subplot(2,3,6),imshow(hi2);title('僅平滑亮度分量')實(shí)驗(yàn)3_2圖片實(shí)驗(yàn)四 車牌識(shí)別圖像預(yù)處理技術(shù)要求：1.對(duì)原始車牌圖像做增強(qiáng)處理；2.對(duì)增強(qiáng)后的彩色圖像進(jìn)行灰度變換；3.對(duì)灰度圖像進(jìn)行直方圖均衡處理；4.選取自適應(yīng)的閾值，對(duì)圖像做二值化處理；5.顯示每步處理后的圖像；6.分析此種圖像預(yù)處理的優(yōu)缺點(diǎn)及改進(jìn)措施，簡(jiǎn)要敘述車牌字符識(shí)別。程序i = imread('C:UsersYSSADesktop數(shù)字圖像課設(shè)于鵬飛圖片sd.jpg');figure,imshow(i);title('原始圖像');% 灰度變換g=rgb2gray(i);figure,imshow(g);title('灰度變換');% 均衡化h=histeq(g);figure,imshow(h);title('直方圖均衡化');% 二值化處理e = im2bw(h,0.9);figure,imshow(e);title('二值化圖像(閾值0.9)');實(shí)驗(yàn)4圖片專心-專注-專業(yè)內(nèi)蒙古