河北工業(yè)大學(xué)數(shù)字圖像處理實(shí)驗(yàn)報(bào)告

1、數(shù)字圖像處理實(shí)驗(yàn)報(bào)告計(jì)算機(jī)科學(xué)與軟件學(xué)院軟件 C103XXXXXXXXX實(shí)驗(yàn)一熟悉MATLAB軟件的開(kāi)發(fā)環(huán)境、基本操作以及圖像處理工具箱，為編寫(xiě)圖像處理程序 奠定基礎(chǔ)。1矩陣：A=1 2 3;4 5 6;2讀取圖像，并進(jìn)行顯示匸imread( xy.jpg ;Imshow(l);二.說(shuō)明以下程序的顯示結(jié)果為何是一幅幾乎全黑的圖像。clear; close all;myi=zeros(20,20);myi(2:2:18,2:2:18)=1;myi=ui nt8(myi);figure, imshow(myi, no truesize);首先1中的clear是清除當(dāng)前內(nèi)存中的變量，close al

2、l是關(guān)閉所有窗口 (程序運(yùn)行產(chǎn)生的， 不包括命令窗，editor窗和幫助窗)2中初始化了一個(gè) 20X20的矩陣，每個(gè)元素的值都是0并且這個(gè)矩陣的名字是 myi。3從矩陣的第2行開(kāi)始，每次走2行，到第18行，從矩陣的第2列開(kāi)始，每次走2列， 到第18列，讓這些地方的元素值為 1.4.uint8這個(gè)函數(shù)是把數(shù)值矩陣轉(zhuǎn)換為圖像矩陣，就是說(shuō)最大值255表示白色，最小值0表示黑色。但是轉(zhuǎn)換后矩陣的值和原來(lái)的一樣。也就是說(shuō)在3中矩陣的值變到 4后沒(méi)變，所以因?yàn)檫@些值要么仍是 0要么仍是1,這樣顯示以后會(huì)看起來(lái)都是黑色。5.imshow里的noturesize參數(shù)就是設(shè)定圖片去適應(yīng)窗口，窗口大小不動(dòng)，自動(dòng)調(diào)

3、節(jié)圖片大小。所以顯示出來(lái)的圖片跟實(shí)際尺寸可能不一致。三.編寫(xiě)程序?qū)D一中orangutan_1.tif圖片拉伸一倍后形成orangutan_2.tif所示的圖片A= imread(C:Usersjaso nDesktopora nguta n_2.tif);imshow(A);figuretform = maketform(affi ne,1.414 0 0; 0 1 0; 0 0 1);B = imtra nsform(A,tform);imshow(B);四、應(yīng)用MATLAB (或C)語(yǔ)言編寫(xiě)一幅灰度圖像直方圖統(tǒng)計(jì)程序，并選擇一幅圖像顯示其直方圖，將結(jié)果與MATLAB圖像處理工具箱中提供的

4、灰度直方圖函數(shù)imhist的處理結(jié)果進(jìn)行比較。clear;i=imread(2.jpg);x=rgb2gray(i);figure;subplot(2,3,1),imshow(i),title(原圖);subplot(2,3,3),imhist(x),title(灰度直方圖)五、利用以上編寫(xiě)的程序，估算圖一所示圖像iris.tif中的瞳孔半徑(以像素為單位)。程序：clear all;A=imread(ya n.j pg);B=rgb2gray(A); imhist(B);m n=size(B); cout=0;for i=1:mfor j=1: nif(B(i,j)20) cout=cout

5、+1;endendenddisp(cout=, nu m2str(cout);r=sqrt(cout/ 3.14); disp(r=, nu m2str(r);結(jié)果圖：面積的結(jié)果是cout=4429，半徑是r=37.5568實(shí)驗(yàn)一.二維傅里葉變換：程序：clear all;img=imread(C:Usersjas onDesktople na512.jpg); img仁img;subplot(2,2,1),imshow(img1);title(灰度圖像);F=fft2(img1);subplot(2,2,2),imshow(F,);title(傅立葉變換后，高頻部分在四角);F=abs(F)

6、;F仁 fftshift(F); subplot(2,2,3),imshow(F1,);title(把高頻變到中間)；F2=1+log(F1);subplot(2,2,4),imshow(F2,);title(指數(shù)優(yōu)化后顯示);處理結(jié)果：二.直方圖均衡化：利用Matlab提供的可進(jìn)行圖像直方圖修正的函數(shù)，自己選擇幾幅直方圖不均勻的圖像（如 圖二pout.tif），對(duì)這些圖像進(jìn)行直方圖均衡處理，顯示處理前后的圖像以及它們的灰度直方 圖，體會(huì)直方圖均衡化算法的特點(diǎn)clear all;A=imread(C:Usersjaso nDesktople na512.jpg);B=A;subplot(2,2

7、,1),imshow(B);title(均衡化之前);subplot(2,2,2),imhist(B);title(均衡化之前的直方圖);C=histeq(B);subplot(2,2,3),imshow(C);title(均衡化之后);subplot(2,2,4),imhist(C);title(均衡化之后的直方圖);三灰度修正技術(shù)clear all;A = imread(C:Usersjas onDesktoppout.tif);RGB仁 A;C=fspecial(gaussia n); subplot(2,4,1),imshow(RGB1); title( 原始圖像); subplot(

8、2,4,2),imhist(RGB1);title(原始圖像的灰度直方圖);RGB2 = imadjust(RGB1,.2 .3 0; .6 .7 1,0.3); subplot(2,4,3),imshow(RGB2);title( 參數(shù)為0.3,灰度修正后); subplot(2,4,4),imhist(RGB2);title( 參數(shù)為0.3灰度修正后的灰度直方圖);RGB2 = imadjust(RGB1,.2 .3 0; .6 .7 1,1); subplot(2,4,5),imshow(RGB2);title( 參數(shù)為1,灰度修正后); subplot(2,4,6),imhist(RG

9、B2);title( 參數(shù)為1灰度修正后的灰度直方圖);RGB2 = imadjust(RGB1,.2 .3 0; .6 .7 1,3); subplot(2,4,7),imshow(RGB2);title( 參數(shù)為3,灰度修正后); subplot(2,4,8),imhist(RGB2);title（ 參數(shù)為3,灰度修正后的灰度直方圖）;原始圖唆的灰塩白方圖曲數(shù)為0一3灰痕憐正啟的灰度白.方圖參數(shù)為1,撫度條正后痙數(shù)為灰度幡正后的灰度直方圖參數(shù)為亂灰:應(yīng)條正后發(fā)第為3灰度懷正后的灰曳玄方圖四圖像平滑(去噪)i=imread(C:Usersjas onDesktople na512.jpg);

10、p=i mno ise(i,gaussia n,0.1); f1=0 1 0;1 0 1;0 1 0/4; % f2=0 1 0;1 0 1;0 1 0/8; % j=imfilter(p,f1,corr,replicate); k=imfilter(p,f2,corr,replicate); subplot(2,2,1),imshow(i),title( subplot(2,2,2),imshow(p),title( subplot(2,2,3),imshow(j),title(4- subplot(2,2,4),imshow(k),title(8-所用模板（4-鄰域）所用模板（8-鄰域）原

11、圖）;噪聲污染后圖）;鄰域處理后圖）;鄰域處理后圖）;File Edit View Insert Tools Desktop Vindaw Help aaibid 噪聲污舉后團(tuán)兒鄰域處理后團(tuán)8-鄰域處理后圖五.邊緣檢測(cè)技術(shù):clear all;A=imread(C:Usersjaso nDesktopora nguta n_ 2.tif); subplot(2,3,1),imshow(A);title(原始圖像);B=im2bw(A,0.95); subplot(2,3,2),imshow(B); title(原始二值圖像);C=double(B);D=edge(C,sobel); subplot(2,3,3),imshow(D); title(sobel邊緣提取方法);E=edge(C,prewitt); subplot(2,3,4),imshow(E); title(prewitt邊緣提取方法);F=edge(C,roberts); subplot(2,3,5),imshow(F); title(roberts邊緣提取方法);G=edge