第一篇：數(shù)字圖像處理——彩色圖像實驗報告

6.3實驗步驟

(1)對彩色圖像的表達和顯示

* * * * * * * * * * * *顯示彩色立方體* * * * * * * * * * * * *

rgbcube(0,0,10);

%從正面觀察彩色立方體

rgbcube(10,0,10);%從側(cè)面觀察彩色立方 rgbcube(10,10,10);%從對角線觀察彩色立方體

%* * * * * * * * * *索引圖像的顯示和轉(zhuǎn)換* * * * * * * * * f=imread('D:PictureFig0604(a)(iris).tif');figure,imshow(f);

%f是RGB真彩圖像

%rgb圖像轉(zhuǎn)換成8色索引圖像，不采用抖動方式

[X1,map1]=rgb2ind(f,8,'nodither');figure,imshow(X1,map1);

%采用抖動方式轉(zhuǎn)換到8色索引圖像

[X2,map2]=rgb2ind(f,8,'dither');figure,imshow(X2,map2);%顯示效果要好一些 g=rgb2gray(f);

%f轉(zhuǎn)換為灰度圖像

g1=dither(g);

%將灰色圖像經(jīng)過抖動處理，轉(zhuǎn)換打二值圖像figure,imshow(g);

%顯示灰度圖像

figure,imshow(g1);

%顯示抖動處理后的二值圖像

程序運行結(jié)果：

*

彩色立方體原圖

不采用抖動方式轉(zhuǎn)換到8色索引圖像采用抖動方式轉(zhuǎn)換到8色索引圖像

灰度圖像抖動處理后的二值圖像

(2)彩色空間轉(zhuǎn)換

f=imread('D:PictureFig0604(a)(iris).tif');figure,imshow(f);

%轉(zhuǎn)換到NTSC彩色空間

%f是RGB真彩圖像

ntsc_image=rgb2ntsc(f);

figure,imshow(ntsc_image(:,:,1));%顯示亮度信息 figure,imshow(ntsc_image(:,:,2));%顯示色差信息 figure,imshow(ntsc_image(:,:,3));%顯示色差信息

%轉(zhuǎn)換到HIS彩色空間

hsi_image=rgb2hsi(f);

figure,imshow(hsi_image(:,:,1));%顯示色度信息 figure,imshow(hsi_image(:,:,2));%顯示飽和度信息 figure,imshow(hsi_image(:,:,3));%顯示亮度信息

程序運行結(jié)果：

原圖

轉(zhuǎn)換到NTSC彩色空間

顯示亮度信息顯示色差信息顯示色差信息

轉(zhuǎn)換到HIS彩色空間

顯示色差信息顯示飽和度信顯示亮度信息

(3)彩色變換

f=imread('D:PictureFig0614(a)(Chalk Original).tif');

G=ice('image',f);

%打開ice窗口對圖像進行調(diào)整

%在窗口中執(zhí)行以下操作：

%a）得到圖像的補色

%b）拖動映射曲線，對圖像顯示效果進行修改

%c）在顏色通道中選中某一顏色，然后對映射曲線進行修改

程序運行結(jié)果(1)：

全彩色圖片ICE窗口

它的補色ICE窗口

拖動映射曲線，圖像的顯示效果

ICE窗口

f2=imread('D:Picture JLK Magenta.tif');figure,imshow(f2);%在CMYK彩色空間內(nèi)打開圖像 選擇Red

ICE窗口

g2=ice('image',f2,'space','CMYK');

%f2的圖像色彩偏紅，拖動映射曲線，%調(diào)整映射參數(shù)，使圖像的色彩看起來比較正常。

程序運行結(jié)果(2)

原圖

在CMYK彩色空間內(nèi)打開圖像

調(diào)整映射參數(shù)后

ICE窗口

(4)彩色空間濾波

%* * * * * * * * * * * * * * * * * * * * 彩色圖像平滑 * * * * * * * * * * * * * * * * * * * * * *

fc=imread('D:PictureFig0604(a)(iris).tif');figure,imshow(fc);h=rgb2hsi(fc);

H=h(:,:,1);

S=h(:,:,2);

I=h(:,:,3);

%轉(zhuǎn)換到HIS彩色空間

%色度分量

%飽和度分量

%亮度分量

%25×25的方形平滑算子 w=fspecial('average',25);

I_filtered=imfilter(I,w,'replicate');%對亮度分量進行平滑操作 H1=cat(3,H,S,I_filtered);%將H、S、I三個分量重新組合起來f=hsi2rgb(H1);

f=min(f,1);

%轉(zhuǎn)換到RGB彩色空間

%f的值控制在[0,1]之間

figure,imshow(f);

%彩色圖像平滑處理 %* * * **** * * * * * * * * * * * * 彩色圖像銳化 * * * * * * * * * * * * * * * * * * * * * *

fb=imread('D:PictureFig0604(a)(iris).tif');lapmask=[1 1 1;1-8 1;1 1 1];%拉普拉斯算子 g=imfilter(fb,lapmask,'replicate');%拉普拉斯算子處理 figure,imshow(g);

%圖像銳化后與原圖像相減

fen=imsubtract(fb,g);

figure,imshow(fen);%顯示差值圖像

程序運行結(jié)果：

原圖彩色圖像平滑處理

拉普拉斯算子處理

銳化處理

(5)彩色圖像分割

%* * * * * * * * * * * * * * * * * 彩色邊緣檢測* * * * * * * * * * * *

close all;clear;clc;

f=imread('iris.tif');%f如圖6-4所示

%直接計算彩色邊緣VG，并得到對單獨彩色平面的二維梯度求和形成的梯度PPG

[VG,A,PPG]=colorgrad(f);

figure,imshow(VG,[ ]);

%顯示彩色梯度圖像 figure,imshow(A,[ ]);

%顯示彩色梯度對應(yīng)的角度圖像

figure,imshow(PPG,[ ]);

%顯示分量梯度圖像

g=VG-PPG;

%計算兩種梯度方法得到的梯度之差 figure,imshow(g,[ ]);

程序運行結(jié)果：

%顯示差值圖像

彩色梯度圖像 彩色梯度對應(yīng)的角度圖像

分量梯度圖像兩種梯度差值圖像

%* * * * * * * * * * * * * * * * * 彩色圖像分割 * * * * * * * * * * * * * * * * * * * * * * * * *

f=imread('jupitermoon.tif');%f如圖6-4所示 mask=roipoly(f);

%取得選定區(qū)域圖像的紅色分量

%手工選定圖中紅色區(qū)域

R=immultiply(mask,f(:,:,1));

%取得選定區(qū)域圖像的綠色分量

G=immultiply(mask,f(:,:,2));%取得選定區(qū)域圖像的藍色分量

B=immultiply(mask,f(:,:,3));

g=cat(3,R,G,B);

%取得選定區(qū)域的RGB圖像 figure,imshow(g);

%顯示取得選定區(qū)域

[M,N,K]=size(g);

%選定區(qū)域的大小

I=reshape(g,M*N,3);

%對g中的像素進行重排，按列序拉成一行 idx=find(mask);

%找到mask中的像素在重排后的圖像g中位置 I=double(I(idx,1:3));

%轉(zhuǎn)換到double類型 [C,m]=covmatrix(I);

%計算協(xié)方差矩陣C主對角線上的元素，即方差 %>>sd=sqrt(d)%求得這些元素的平方根，即均方差

均方差

%用歐式距離對圖像進行彩色分割，閾值為25

E25=colorseg('euclidean',f,25,m);E251=f;E251(E25==0)=0 figure,imshow(E251,[]);

%用馬氏距離對圖像進行彩色分割，閾值為25 M25=colorseg('mahalanobis',f,25,m,C);figure , imshow(M25 , [ ]);程序運行結(jié)果：

原圖選定區(qū)域

歐氏距離彩色分割馬氏距離彩色分割

6.4 思考題

(1)請將圖6-1所示的彩色圖像轉(zhuǎn)換到灰度圖像，并根據(jù)灰度值對此灰度圖像進行索引化，手工指定索引表中的RGB值，并把此索引圖像顯示出來。

程序代碼如下

%思路：把原圖想轉(zhuǎn)換為128位索引圖

close all;clear;clc;%將f從灰度圖像轉(zhuǎn)換為索引圖像，x是索引圖像； %map是索引矩陣，R，G，B的值都相等

f=imread('D:PictureRGB_iris.tif');figure,imshow(f);

%f是RGB圖像 f1=rgb2gray(f);

%將f轉(zhuǎn)換成灰度圖像 figure,imshow(f1);[x,map]=gray2ind(f1,128);

figure,imshow(x);%將map中的值做下面改變則R，G，B的值不相等，產(chǎn)生彩色

%把灰度圖像轉(zhuǎn)化成128位的索引圖

map(1,:)=0.5;map(: ,1)=0.5;map(2,:)=0.4;map(1,1)=0;figure,imshow(f1,map);

程序運行結(jié)果：

原RGB圖像

灰度圖像

索引圖 修改索引表后的圖像

(2)6.3節(jié)的（4）是將圖像轉(zhuǎn)換到HIS空間后再對亮度分量進行平滑操作，試直接對原圖像的RGB三個分量分別進行平滑操作，并計算兩種操作方法的圖像差值。

程序代碼如下：

f=imread('D:PictureRGB_iris.tif');figure,imshow(f);%直接對原圖像的RGB三個分量分別進行平滑操作 H=f(:,:,1);

S=f(:,:,2);

I=f(:,:,3);

%色度分量

%飽和度分量

%亮度分量

%25×25的方形平滑算子 w=fspecial('average',25);

H_filtered=imfilter(H,w,'replicate');

%對色度分量進行平滑操作 S_filtered=imfilter(S,w,'replicate');

%對飽和度分量進行平滑操作 I_filtered=imfilter(I,w,'replicate');

%對亮度分量進行平滑操作 %將H、S、I三個分量重新組合起來

Ha=cat(3,H_filtered,S_filtered,I_filtered);Ha=im2double(Ha);

%轉(zhuǎn)化成double類型 figure,imshow(Ha);%轉(zhuǎn)換到HIS彩色空間 h=rgb2hsi(f);

%轉(zhuǎn)換到HIS彩色空間

H=h(:,:,1);

S=h(:,:,2);

I=h(:,:,3);

%色度分量

%飽和度分量 %亮度分量 %25×25的方形平滑算子 w=fspecial('average',25);

I_filtered=imfilter(I,w,'replicate');%對亮度分量進行平滑操作 %將H、S、I三個分量重新組合起來 H1=cat(3,H,S,I_filtered);fb=hsi2rgb(H1);

%轉(zhuǎn)換到RGB彩色空間 fb=min(fb,1);

figure,imshow(fb);

%f的值控制在[0,1]之間

%彩色圖像平滑處理

fen=imsubtract(Ha,fb);figure,imshow(fen);程序運行結(jié)果：

%計算兩中方式處理后圖像的差值

原圖

直接對RGB三分量進行平滑操作轉(zhuǎn)換到HIS空間對亮度進行平滑操作

差值圖像

試驗中的問題和心得體會：

本次實驗是彩色圖像處理。

(1)通過這個實驗我了解了RGB圖像以及索引圖像在MATLAB中的存儲和詢問方式

(2)了解RGB空間到其他彩色空間的轉(zhuǎn)換方式和轉(zhuǎn)換方法(3)知道了一些彩色圖像的空間濾波方法(4)了解在RGB向量空間中檢測彩色邊緣的方法(5)了解在RGB向量空間中對圖像進行分割的方法

遇到的問題：在做思考題(2)的時候，計算兩種平滑處理的得到圖像的差值的中程序出現(xiàn)了問題，經(jīng)過看出錯信息，認(rèn)真檢查程序后，發(fā)現(xiàn)是因為直接平滑處理的圖像是unit8類型的，而轉(zhuǎn)換到HIS空間后，得到的圖像是double類型，二者的類型不同，導(dǎo)致使用imsubtract函數(shù)的時候出現(xiàn)了問題。感覺做matlab的實驗和平時編寫C代碼的時候有很多地方是相同的，都會經(jīng)歷：編碼，debug，正確這三個過程，這個過程本身就具有挑戰(zhàn)性和趣味性。

第二篇：數(shù)字圖像處理圖像變換實驗報告

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

實 驗

實驗名稱：圖像處理姓名：劉強

班級：電信

學(xué)號：

報 告

1102

1404110128

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

實驗一 圖像變換實驗——圖像點運算、幾何變換及正交變換

一、實驗條件

PC機 數(shù)字圖像處理實驗教學(xué)軟件

大量樣圖

二、實驗?zāi)康?/p>

1、學(xué)習(xí)使用“數(shù)字圖像處理實驗教學(xué)軟件系統(tǒng)”，能夠進行圖像處理方面的簡單操作；

2、熟悉圖像點運算、幾何變換及正交變換的基本原理，了解編程實現(xiàn)的具體步驟；

3、觀察圖像的灰度直方圖，明確直方圖的作用和意義；

4、觀察圖像點運算和幾何變換的結(jié)果，比較不同參數(shù)條件下的變換效果；

5、觀察圖像正交變換的結(jié)果，明確圖像的空間頻率分布情況。

三、實驗原理

1、圖像灰度直方圖、點運算和幾何變換的基本原理及編程實現(xiàn)步驟

圖像灰度直方圖是數(shù)字圖像處理中一個最簡單、最有用的工具，它描述了一幅圖像的灰度分布情況，為圖像的相關(guān)處理操作提供了基本信息。

圖像點運算是一種簡單而重要的處理技術(shù)，它能讓用戶改變圖像數(shù)據(jù)占據(jù)的灰度范圍。點運算可以看作是“從象素到象素”的復(fù)制操作，而這種復(fù)制操作是通過灰度變換函數(shù)實現(xiàn)的。如果輸入圖像為A(x,y)，輸出圖像為B(x,y)，則點運算可以表示為：

B(x,y)＝f[A(x,y)] 其中f(x)被稱為灰度變換（Gray Scale Transformation，GST）函數(shù)，它描述了輸入灰度值和輸出灰度值之間的轉(zhuǎn)換關(guān)系。一旦灰度變換函數(shù)確定，該點運算就完全確定下來了。另外，點運算處理將改變圖像的灰度直方圖分布。點運算又被稱為對比度增強、對比度拉伸或灰度變換。點運算一般包括灰度的線性變換、閾值變換、窗口變換、灰度拉伸和均衡等。

圖像幾何變換是圖像的一種基本變換，通常包括圖像鏡像變換、圖像轉(zhuǎn)置、圖像平移、圖像縮放和圖像旋轉(zhuǎn)等，其理論基礎(chǔ)主要是一些矩陣運算，詳細原理可以參考有關(guān)書籍。

實驗系統(tǒng)提供了圖像灰度直方圖、點運算和幾何變換相關(guān)內(nèi)容的文字說明，用戶在操作過程中可以參考。下面以圖像點運算中的閾值變換為例給出編程實現(xiàn)的程序流程圖，如下：

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

2、圖像正交變換的基本原理及編程實現(xiàn)步驟 數(shù)字圖像的處理方法主要有空域法和頻域法，點運算和幾何變換屬于空域法。頻域法是將圖像變換到頻域后再進行處理，一般采用的變換方式是線性的正交變換（酉變換），主要包括傅立葉變換、離散余弦變換、沃爾什變換、霍特林變換和小波變換等。正交變換被廣泛應(yīng)用于圖像特征提取、圖像增強、圖像復(fù)原、圖像壓縮和圖像識別等領(lǐng)域。

正交變換實驗的重點是快速傅立葉變換（FFT），其原理過于復(fù)雜，可以參考有關(guān)書籍，這里不再贅述。至于FFT的編程實現(xiàn)，系統(tǒng)采用的方法是：首先編制一個一維FFT程序模塊，然后調(diào)用該模塊對圖像數(shù)據(jù)的列進行一維FFT，再對行進行一維FFT，最后計算并顯示幅度譜。程序流程圖如下：

四、實驗內(nèi)容

圖像灰度直方圖

點運算：圖像反色、灰度線性變換、閾值變換、窗口變換、灰度拉伸和灰度

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

均衡

幾何變換：圖像鏡像變換、圖像轉(zhuǎn)置、圖像平移、圖像縮放和圖像旋轉(zhuǎn) 正交變換：傅立葉變換、離散余弦變換、沃爾什變換、霍特林變換和小波正反變換

注意：

1、所有實驗項目均針對8位BMP灰度圖像進行處理，其它格式（如JPG）的圖像可以利用系統(tǒng)提供的圖像格式轉(zhuǎn)換工具進行轉(zhuǎn)換，再進行處理；

2、本次實驗的重點是圖像的灰度直方圖和點運算，幾何變換和正交變換只作一般性了解。

五、實驗步驟

以圖像灰度閾值變換為例說明實驗的具體步驟，其它實驗項目的步驟與此類似。

1、打開計算機，在系統(tǒng)桌面上雙擊“數(shù)字圖像處理實驗教學(xué)軟件系統(tǒng)”的可執(zhí)行文件“圖象處理”的圖標(biāo)，進入實驗系統(tǒng)；

2、執(zhí)行文件→打開，在OPEN對話框中選擇待處理的圖像，按【OK】后系統(tǒng)顯示出圖像；

3、執(zhí)行查看→圖像基本信息，將顯示圖像基本信息對話框，如圖所示；

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

4、執(zhí)行查看→灰度直方圖，查看圖像的灰度直方圖，如圖所示；

5、執(zhí)行圖像變換→正交變換→傅立葉變換，查看圖像的頻率域分布情況，如圖所示；

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

6、執(zhí)行圖像變換→正交變換→小波變換，查看圖像經(jīng)過小波變換的效果，如圖所示；

7、執(zhí)行圖像變換→點運算→閾值變換，修改閾值變換對話框中的閾值參數(shù)，如圖所示；

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

8、設(shè)置完閾值參數(shù)后按【OK】,系統(tǒng)顯示閾值變換后的圖像，與原圖像進行比較，觀察閾值變換的效果，如圖所示；

9、重復(fù)步驟4，查看閾值變換后圖像的直方圖分布情況；

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

10、重復(fù)步驟5，查看閾值變換后圖像的頻率域分布情況；

11、執(zhí)行文件→保存或另存為，保存處理后的圖像；

12、執(zhí)行文件→重新加載，重新加載原始圖像，但要注意先前對圖像的處理將會丟失； 注意：

13、在執(zhí)行步驟2時可能會出現(xiàn)有些圖像文件不能打開的情況，如圖所示，此時可以先利用圖像格式轉(zhuǎn)換工具將圖像文件轉(zhuǎn)換為8位BMP圖像，再利用系統(tǒng)進行處理。步驟14和15是使用圖像格式轉(zhuǎn)換工具的方法；

14、在桌面上雙擊圖像格式轉(zhuǎn)換工具Jpg2bmp的圖標(biāo)，進入轉(zhuǎn)換工具界面，如圖所示；

15、按照界面提示，把JPG格式的圖像文件轉(zhuǎn)換成8位BMP圖像。

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

步驟13示意圖

步驟14示意圖

六、思考題

1、圖像灰度線性變換、閾值變換、窗口變換、灰度拉伸和灰度均衡之間有何區(qū)別？

灰度線性變換就是將圖像的像素值通過指定的線性函數(shù)進行變換，以此增強或者減弱圖像的灰度。

灰度的閾值變換可以讓一幅圖像變成黑白二值圖。

灰度的窗口變換也是一種常見的點運算。它的操作和閾值變換類似。從實現(xiàn)方法上可以看作是灰度折線變換的特列。窗口灰度變換處理結(jié)合了雙固定閾值法，與其不同之處在于窗口內(nèi)的灰度值保持不變。

灰度拉伸又叫做對比度拉伸，它與線性變換有些類似，不同之處在于灰度拉伸使用的是分段線性變換，所以它最大的優(yōu)勢是變換函數(shù)可以由用戶任意合成。

灰度均衡是增強圖像的有效方法之一。灰度均衡同樣屬于改進圖像的方法，灰度均衡的圖像具有較大的信息量。從變換后圖像的直方圖來看，灰度分布更加均勻。

2、利用圖像鏡像和旋轉(zhuǎn)變換可以實現(xiàn)圖像轉(zhuǎn)置嗎？如果可以，應(yīng)該怎樣實現(xiàn)？

可以。進行一次鏡像變換，順（逆）時針旋轉(zhuǎn)兩次，再以與第一次相反的方向鏡像變換。

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

實驗二 圖像增強及復(fù)原實驗

七、實驗條件

PC機 數(shù)字圖像處理實驗教學(xué)軟件

大量樣圖

八、實驗?zāi)康?/p>

1、熟練使用“數(shù)字圖像處理實驗教學(xué)軟件系統(tǒng)”；

2、熟悉圖像增強及復(fù)原的基本原理，了解編程實現(xiàn)的具體步驟；

3、觀察圖像中值濾波、平滑、銳化和偽彩色編碼的結(jié)果，比較不同參數(shù)條件下的圖像增強效果；

4、觀察圖像退化和復(fù)原的結(jié)果，比較不同復(fù)原方法的復(fù)原效果。

九、實驗原理

1、圖像增強和復(fù)原的基本原理

對降質(zhì)圖像的改善處理通常有兩類方法：圖像增強和圖像復(fù)原。

圖像增強不考慮圖像降質(zhì)的原因，只將圖像中感興趣的特征有選擇地進行突出，并衰減圖像的次要信息，改善后的圖像不一定逼近原始圖像，只是增強了圖像某些方面的可讀性，如突出了目標(biāo)輪廓，衰減了各種噪聲等。圖像增強可以用空域法和頻域法分別實現(xiàn)，空域法主要是在空間域中對圖像象素灰度值直接進行運算處理，一般包括中值濾波、模板平滑和梯度銳化等，空域法可以用下式來描述：

g(x,y)＝f(x,y)*h(x,y)其中f(x,y)是處理前圖像，g(x,y)表示處理后圖像，h(x,y)為空間運算函數(shù)。圖像增強的頻域法是在圖像的頻率域中對圖像的變換值進行某種運算處理，然后變換回空間域，系統(tǒng)涉及的各種濾波器屬于頻域法增強，這是一種間接處理方法，可以用下面的過程模型來描述：

其中：F(u,v)=?[ f(x,y)]，G(u,v)= F(u,v)H(u,v)，g(x,y)=??1[ G(u,v)]，?和??1分別表示頻域正變換和反變換。實驗系統(tǒng)提供了圖像增強相關(guān)內(nèi)容的文字說明，用戶在操作過程中可以參考。

圖像復(fù)原是針對圖像降質(zhì)的原因，設(shè)法去補償降質(zhì)因素，使改善后的圖像盡可能逼近原始圖像，提高了圖像質(zhì)量的逼真度。關(guān)于圖像復(fù)原的詳細原理可以參考相關(guān)書籍，這里不再贅述。本系統(tǒng)提供了圖像的噪聲退化、卷積退化和運動模糊退化操作，并提供了相應(yīng)的逆濾波復(fù)原、維納復(fù)原和運動模糊復(fù)原操作。本次

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

實驗中圖像復(fù)原只作一般性了解。

2、編程實現(xiàn)步驟

下面以圖像增強中的中值濾波操作為例給出編程實現(xiàn)的程序流程圖，如下：

十、實驗內(nèi)容

圖像增強：中值濾波、圖像模板平滑、理想低通濾波器平滑、巴特沃斯低通濾波器平滑、梯度銳化、拉普拉斯銳化、理想高通濾波器銳化、巴特沃斯高通濾波器銳化和偽彩色編碼

圖像復(fù)原：圖像的噪聲退化、卷積退化、卷積加噪聲退化、運動模糊退化、逆濾波復(fù)原、維納復(fù)原和運動模糊復(fù)原

注意：

3、所有實驗項目均針對8位BMP灰度圖像進行處理；

4、本次實驗的重點是圖像增強中的中值濾波和模板平滑，圖像復(fù)原只作一般性了解。

十一、實驗步驟

以圖像中值濾波操作為例說明實驗的具體步驟，其它實驗項目的步驟與此類似。

11、打開計算機，在系統(tǒng)桌面上雙擊“數(shù)字圖像處理實驗教學(xué)軟件系統(tǒng)”的可執(zhí)行文件“圖象處理”的圖標(biāo)，進入實驗系統(tǒng)；

12、執(zhí)行文件→打開，在OPEN對話框中選擇待處理的圖像，按【OK】后系統(tǒng)顯示出圖像；

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

13、執(zhí)行查看→圖像基本信息，將顯示圖像基本信息對話框，如圖所示；

14、執(zhí)行查看→灰度直方圖，查看圖像的灰度直方圖，如圖所示；

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

15、執(zhí)行圖像變換→正交變換→傅立葉變換，查看圖像的頻率域分布情況，如圖所示；

16、執(zhí)行圖像增強→中值濾波，選擇或自定義對話框中的濾波器參數(shù)，如圖所示；

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

17、設(shè)置完濾波器參數(shù)后按【OK】,系統(tǒng)顯示中值濾波后的圖像，與原圖像進行比較，觀察中值濾波的效果，如圖所示；

18、重復(fù)步驟4，查看中值濾波后圖像的直方圖分布情況；

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

19、重復(fù)步驟5，查看中值濾波后圖像的頻率域分布情況；

10、執(zhí)行文件→保存或另存為，保存處理后的圖像；

11、執(zhí)行文件→重新加載，重新加載原始圖像，但要注意先前對圖像的處理將會丟失。

數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

十二、思考題

1、圖像中值濾波和模板平滑之間有何區(qū)別？

圖像平滑處理就是用平滑模板對圖像進行處理，以減少圖像的噪聲。而中值濾波是一種非線性的信號處理方法。

2、圖像增強和圖像復(fù)原之間有何區(qū)別？

圖像增強：利用一定的技術(shù)手段，不用考慮圖像是否失真（即原 始圖像在變換后可能會失真）而且不用分析圖像降質(zhì)的原因。針對給定圖像的應(yīng)用場合，有目的地強調(diào)圖像的整體或局部特性，將原來不清晰的圖像變得清晰或強調(diào)某些感興趣的特征，擴大圖像中不同物體特征之間的差別，抑制不感興趣的特征，使之改善圖像質(zhì)量、豐富信息量，加強圖像判讀和識別效果，滿足某些特殊分析的需要。

圖像復(fù)原：針對質(zhì)量降低或者失真的圖像，恢復(fù)圖像原始的內(nèi)容或者質(zhì)量。圖像復(fù)原的過程包含對圖像退化模型的分析，再對退化的圖像進行復(fù)原。圖像退化是由于成像系統(tǒng)受各種因素的影響，導(dǎo)致了圖像質(zhì)量的降低，稱之為圖像退化。這些因素包括傳感器噪聲、攝像機聚焦不佳、物體與攝像機之間的相對移動、隨機大氣湍流、光學(xué)系統(tǒng)的象差、成像光源和射線的散射等。圖像復(fù)原大致可以分為兩種方法：

一種方法適用于缺乏圖像先驗知識的情況，此時可對退化過程建立模型進行描述，進而尋找一種去除或消弱其影響的過程，是一種估計方法；

另一種方法是針對原始圖像有足夠的先驗知識的情況，對原始圖像建立一個數(shù)學(xué)模型并根據(jù)它對退化圖像進行擬合，能夠獲得更好的復(fù)原效果。

3、圖像維納復(fù)原為什么比逆濾波復(fù)原效果好？

維納濾波復(fù)原的原理可表示為

對于維納濾波，由上式可知，當(dāng)

時，由于存在 項，所以數(shù)字圖象處理實驗指導(dǎo)書

實驗一 圖象變換實驗

不會出現(xiàn)被0除的情形，同時分子中含有項，在處。當(dāng)時，此時維納濾波就變成 了逆濾波；當(dāng)時，表明維納濾波避免了逆濾波中 出現(xiàn)的對噪聲過多的放大作用,也就是說圖像維納復(fù)原比逆濾波復(fù)原效果好。

第三篇：數(shù)字圖像處理實驗報告

實 驗 報 告 書

系部學(xué)生專業(yè)班

實驗

名稱姓名名稱 級時間：

： ： ：

：

實驗一 直方圖均衡

一、實驗?zāi)康?/p>

在學(xué)習(xí)圖像直方圖的概念、計算方法、性質(zhì)和相關(guān)應(yīng)用基礎(chǔ)上，生成、繪制圖像的直方圖，并應(yīng)用MATLAB編程實現(xiàn)圖像直方圖均衡化程序。

二、實驗內(nèi)容

（1）計算并繪制圖像直方圖；

（2）編程實現(xiàn)圖像的直方圖均衡化處理，顯示均衡前后的直方圖和圖像；

三、實驗運行結(jié)果

四、實驗中遇到的問題及解決方法

1、顯示無法找到圖像文件，應(yīng)將圖片與xx.m文件置于同一文件夾；

2、編程過程中應(yīng)注意標(biāo)點的輸入法，應(yīng)該用英文輸入，否則會報錯。

3、編程完成后運行時輸入文件名與保存時文件名相同，區(qū)分大小寫。

五、思考題

（1）、灰度直方圖可以反映一幅圖像的哪些特征？ 答：

1、表征了圖像的一維信息。只反映圖像中像素不同灰度值出現(xiàn)的次數(shù)（或頻數(shù)）而未反映像素所在位置。

2、與圖像之間的關(guān)系是多對一的映射關(guān)系。一幅圖像唯一確定出與之對應(yīng)的直方圖，但不同圖像可能有相同的直方圖。

3、子圖直方圖之和為整圖的直方圖。(2)均衡化后的直方圖有何特點？

答：經(jīng)直方圖均衡化處理后，可以得到一副改善了質(zhì)量的新圖像。這幅圖像的灰度層次將不再是呈黑暗色調(diào)的圖像，而是一副灰度層次較為適中的、比原始圖像清晰、明快得多的圖像。處理的結(jié)果使圖像更適合與人的視覺特征或機器的識別系統(tǒng)。

六、實驗心得體會

本次實驗中，因為初學(xué)這個軟件，我學(xué)習(xí)到了在程序中關(guān)于圖像的運用，以及也復(fù)習(xí)了課本上的許多知識，加深了對直方圖均衡化的理解。

七、程序清單

clear all;I=imread('lena_gray_256.tif');%打開一幅灰度圖像 [m,n]=size(I);p=m*n;J=imhist(I)./p;

%計算圖像的歸一化直方圖 subplot(1,3,1),imshow(I)； subplot(1,3,2),imhist(I,64)； subplot(1,3,3),plot(J)；

（2）直方圖均衡化

clear all;Im=imread('region.jpg');J=histeq(Im);%均衡化 subplot(2,2,1);imshow(Im);title('原圖');%顯示原圖 subplot(2,2,2);imhist(Im);title('原圖直方圖');%顯示原圖的直方圖 subplot(2,2,3);imshow(J);title('均衡化結(jié)果');%顯示均衡化后的圖像 subplot(2,2,4);imhist(J);title('均衡化結(jié)果的直方圖');%顯示均衡化后的直方圖

實驗二 頻域圖像增強

一、實驗?zāi)康?/p>

1、頻域圖像增強

2、掌握基于頻域的圖像增強方法。

二、實驗內(nèi)容

（1）編程實現(xiàn)圖像的理想低通和高通濾波；（2）編程實現(xiàn)圖像的巴特沃斯低通和高通濾波。

三、實驗運行結(jié)果

四、實驗中遇到的問題及解決方法

顯示圖像無法打開，最終查出來時圖像格式弄錯了。

五、思考題

分析為什么圖像通過低通濾波器后變得模糊？為什么通過高通濾波器后得到銳化結(jié)果？

答：圖像的精細結(jié)構(gòu)及突變部分主要由高頻成分起作用，故經(jīng)低通濾波后圖像的精細結(jié)構(gòu)消失，變得模糊；經(jīng)高通濾波后圖像得到銳化。

六、實驗心得體會

本實驗中遇到很多問題及錯誤，例如圖像打不開、處理后圖像模糊等，都是經(jīng)常容易發(fā)生的錯誤，最后實驗幾次，就能夠逐一自己解決了。使自己對數(shù)字圖像處理課程中的許多問題有了更實際和確切的深入了解。

七、程序清單

clc;clear;data4=imread('lena.gif');subplot(3,2,1);imshow(data4);title('原圖');i=fft2(data4);subplot(3,2,2);i=fftshift(i);z=log(abs(i));x=0:1:255;y=0:1:255;[x,y]=meshgrid(x,y);mesh(z);%以三維坐標(biāo)顯示該圖像頻譜圖title('原圖頻譜');[n,m]=size(i);%for k=1:1:n for l=1:1:m if(k^2+l^2)>=190^2 % result(k,l)=0;else result(k,l)=i(k,l);end end end subplot(3,2,4);

對該圖進行低通濾波 選取D=190

z=log(abs(result));%三維方式顯示低通濾波后的頻譜圖 x=0:1:255;y=0:1:255;[x,y]=meshgrid(x,y);mesh(z);title('理想低通濾波后的頻譜');subplot(3,2,3);%新建圖像顯示窗口 result=fftshift(result);%濾波后的數(shù)據(jù)去中心化 b=ifft2(result);%逆傅里葉變換 imshow(uint8(abs(b)));title('理想低通濾波后的圖像');subplot(3,2,6);%新建圖像顯示窗口 % [n,m]=size(c);%對原圖進行高通濾波 for k=1:1:n for l=1:1:m if(k^2+l^2)<=190^2 %選取D=190 result(k,l)=0;else result(k,l)=i(k,l);end end end z=log(abs(result));x=0:1:255;%三維方式顯示高通濾波前的頻譜圖 y=0:1:255;[x,y]=meshgrid(x,y);mesh(z);title('理想高通濾波后的頻譜');subplot(3,2,5);result=fftshift(result);%濾波后的數(shù)據(jù)去中心化 d=ifft2(result);%逆傅里葉變換

imshow(uint8(abs(d)));title('理想高通濾波后的圖像');%頻域增強（巴特沃斯原型）

%二階巴特沃斯（Butterworth）低通濾波器 %clc;%clear;Figure;J1=imread('lena.gif');subplot(3,2,1);imshow(J1);title('原圖');f=double(J1);g=fft2(f);% 傅立葉變換 g=fftshift(g);% 轉(zhuǎn)換數(shù)據(jù)矩陣 subplot(3,2,2);x=0:1:255;y=0:1:255;[x,y]=meshgrid(x,y);z=log(abs(g));%取幅度 mesh(z);%以三維坐標(biāo)顯示該圖像頻譜圖 title('原圖頻譜');[M,N]=size(g);nn=2;% 二階巴特沃斯(Butterworth)低通濾波器 d0=20;m=fix(M/2);n=fix(N/2);for i=1:M for j=1:N d=sqrt((i-m)^2+(j-n)^2);h=1/(1+0.414*(d/d0)^(2*nn));% 計算低通濾波器傳遞函數(shù) result(i,j)=h*g(i,j);

end end subplot(3,2,4);x=0:1:255;y=0:1:255;[x,y]=meshgrid(x,y);z=log(abs(result));%取幅度 mesh(z);%以三維坐標(biāo)顯示該圖像頻譜圖 title('低通濾波后的頻譜');result=ifftshift(result);J2=ifft2(result);J3=uint8(abs(J2));subplot(3,2,3);imshow(J3);title('低通濾波后的圖像');

%利用二階巴特沃斯（Butterworth）高通濾波器

nn=2;% 二階巴特沃斯(Butterworth)高通濾波器 d0=5;m=fix(M/2);n=fix(N/2);for i=1:M for j=1:N d=sqrt((i-m)^2+(j-n)^2);if(d==0)h=0;else h=1/(1+0.414*(d0/d)^(2*nn));% 計算傳遞函數(shù) end result(i,j)=h*g(i,j);

end end subplot(3,2,6);x=0:1:255;y=0:1:255;[x,y]=meshgrid(x,y);z=log(abs(result));%取幅度 mesh(z);%以三維坐標(biāo)顯示該圖像頻譜圖 title('高通濾波后的頻譜');result=ifftshift(result);J2=ifft2(result);J3=uint8(abs(J2));subplot(3,2,5);imshow(J3);title('高通濾波后的圖像');

實驗三 圖像邊緣檢測與連接

一、實驗?zāi)康?/p>

圖像邊緣檢測與連接

二、實驗內(nèi)容

（1）編程實現(xiàn)一階差分邊緣檢測算法，包括Robert梯度算子、Prewitt算子、Sobel算子等；（2）編程實現(xiàn)二階差分拉普拉斯邊緣檢測算法以及LoG檢測法和Canny檢測法；（3）分析與比較各種邊緣檢測算法的性能；

（4）編程實現(xiàn)Hough變換提取直線

（5）分析Hough變換檢測性能；

三、實驗運行結(jié)果

四、實驗中遇到的問題及解決方法

拷貝文件后沒改文件名，直接執(zhí)行時出現(xiàn)錯誤，最后重新修改后重新編譯，使之成功。

五、思考題

（1）邊緣的方向是什么意思？為什么要考慮邊緣的方向？

答：邊緣常常意味著一個區(qū)域的終結(jié)和另一個區(qū)域的開始，圖像的邊緣也包含了物體的形狀的重要信息，他不僅在分析圖像時大幅度的減少了要處理的信息量，而且還保護了目標(biāo)的邊界結(jié)構(gòu)。所以考慮邊緣的方向很重要。

（2）Hough變換原理是什么？

答：Hough變換的基本原理在于利用點與線的對偶性，將原始圖像空間的給定的曲線通過曲線表達形式變?yōu)閰?shù)空間的一個點。這樣就把原始圖像中給定曲線的檢測問題轉(zhuǎn)化為準(zhǔn)找參數(shù)空間的峰值問題。也即把檢測整體特性轉(zhuǎn)化為檢測局部特性。比如直線、橢圓、圓、弧線等。

六、實驗心得體會

對于一些圖像處理的函數(shù)不是很了解，只能夠按課本的參照函數(shù)拷貝做實驗，對于其中的一些函數(shù)問題理解不是很透徹，有些甚至完全不懂。還得繼續(xù)努力。

七、程序清單

1、邊緣檢測

由edge函數(shù)實現(xiàn)各算子對圖像的邊緣檢測

clear all;I = imread('d:office.bmp');I=rgb2gray(I);BW1 = edge(I,'sobel');

%利用Sobel算子進行邊緣檢測 BW2 = edge(I,'roberts');%利用roberts算子進行邊緣檢測 BW3 = edge(I,'prewitt');%利用prewitt算子進行邊緣檢測 BW4 = edge(I,'log');

%利用log算子進行邊緣檢測 BW5 = edge(I,'canny');

%利用canny算子進行邊緣檢測 subplot(2,3,1),imshow(I)subplot(2,3,2),imshow(BW1)subplot(2,3,3),imshow(BW2)subplot(2,3,4),imshow(BW3)

subplot(2,3,5),imshow(BW4)subplot(2,3,6),imshow(BW5)

2、邊緣連接

使用Hough變換作線檢測和連接

clear all;

RGB = imread('d:M_M.bmp');I=RGB;%I = rgb2gray(RGB);BW = edge(I,'canny');

% 利用Canny算子提取圖像邊緣 [H,T,R] = hough(BW,'RhoResolution',0.5,'ThetaResolution',0.5);figure(1), imshow(T,R,H,[],'notruesize'), axis on, axis normal xlabel('T'), ylabel('R')p = houghpeaks(H,5,'threshold',ceil(0.3*max(H(:))));

%找到5個較明顯的Hough變換峰值

hold on plot(T(p(:,2)),R(p(:,1)),'s','color','white');lines = houghlines(BW,T,R,p,'FillGap',10,'MinLength',10);

%查找并鏈接線段

figure, imshow(BW), hold on %在二值圖中疊加顯示這些線段 for k = 1:length(lines)

xy = [lines(k).point1;lines(k).point2];

plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','green');end

第四篇：數(shù)字圖像處理實驗報告

數(shù)字圖像處理

實驗報告

目錄

1.數(shù)字圖像處理簡介

2.實驗?zāi)康?.實驗內(nèi)容

4.實驗結(jié)果及代碼展示

5.算法綜述

6.Matlab優(yōu)勢

7.總結(jié)

8.存在問題

一、數(shù)字圖像處理簡介

圖像處理，是對圖像進行分析、加工、和處理，使其滿足視覺、心理以及其他要求的技術(shù)。圖像處理是信號處理在圖像域上的一個應(yīng)用。目前大多數(shù)的圖像是以數(shù)字形式存儲，因而圖像處理很多情況下指數(shù)字圖像處理。此外，基于光學(xué)理論的處理方法依然占有重要的地位。

圖像處理是信號處理的子類，另外與計算機科學(xué)、人工智能等領(lǐng)域也有密切的關(guān)系。

傳統(tǒng)的一維信號處理的方法和概念很多仍然可以直接應(yīng)用在圖像處理上，比如降噪、量化等。然而，圖像屬于二維信號，和一維信號相比，它有自己特殊的一面，處理的方式和角度也有所不同。

二、實驗?zāi)康?/p>

鞏固所學(xué)知識，提高所學(xué)能力

三、實驗內(nèi)容

利用matlab的GUI程序設(shè)計一個簡單的圖像處理程序，并含有如下基本功能： 1.讀入一幅RGB圖像，變換為灰度圖像和二值圖像，并在同一個窗口內(nèi)分成三個子窗口來分別顯示RGB圖像和灰度圖像，注上文字標(biāo)題 2.對給定圖像進行旋轉(zhuǎn)

3．對給定的圖像添加噪聲（椒鹽噪聲、高斯噪聲）

四、實驗結(jié)果及代碼展示

1.軟件設(shè)計界面

2.各模塊功能展示以及程序代碼

（1）讀入一幅RGB圖像，變換為灰度圖像和二值圖像，并在同一個窗口內(nèi)分成三個子窗口來分別顯示RGB圖像和灰度圖像，注上文字標(biāo)題

效果展示：

代碼：

a = imread('C:Documents and SettingsAdministrator桌面數(shù)字圖像舞美.JPG');

i = rgb2gray(a);I = im2bw(a,0.5);

subplot(3,1,1);imshow(a);title('源圖像')subplot(3,1,2);imshow(i);title('灰度圖像')subplot(3,1,3);imshow(I);title('二值圖像')

（2）圖像旋轉(zhuǎn) 原圖

效果展示：

代碼：

clc;clear all;close all;

Img=imread('D:My DocumentsMy Pictures5.JPG');Img=double(Img);[h w]=size(Img);alpha=pi/4;

wnew=w*cos(alpha)+h*sin(alpha);hnew=w*sin(alpha)+h*cos(alpha);wnew=ceil(wnew);

hnew=ceil(hnew);u0=w*sin(alpha);

T=[cos(alpha),sin(alpha);-sin(alpha),cos(alpha)];Imgnew2=zeros(hnew,wnew);Imgnew1=zeros(hnew,wnew);for u=1:hnew

for v=1:wnew

tem=T*([u;v]-[u0;0]);x=tem(1);y=tem(2);if x>=1&&x<=h&&y>=1&&y<=w x_low=floor(x);x_up=ceil(x);y_low=floor(y);y_up=ceil(y);if(x-x_low)<=(x_up-x)x=x_low;

else

x=x_up;

end

if(y-y_low)<=(y_up-y)y=y_low;

else

y=y_up;

end

p1=Img(x_low,y_low);p2=Img(x_up,y_low);p3=Img(x_low,y_low);p4=Img(x_up,y_up);s=x-x_low;t=y-y_low;Imgnew1(u,v)=Img(x,y);

Imgnew2(u,v)=(1-s)*(1-t)*p1+(1-s)*t*p3+(1-t)*s*p2+s*t*p4;end

end end

figure;imshow(Imgnew2,[]);B=imrotate(Img,alpha/pi*180);figure;imshow(B,[]);

（3）對給定的圖像添加噪聲（斑點噪聲、高斯噪聲）效果展示：

代碼：

I= imread('D:My DocumentsMy Pictures5.JPG');figure,subplot(211);imshow(I);title('原圖');J1=imnoise(I,'gaussian',0,0.02);

subplot(223);imshow(J);title('添加高斯噪聲');J=imnoise(I,'speckle',0.04);

subplot(224);imshow(J);title('添加斑點噪聲');

五、算法綜述 灰度圖像：

一幅完整的圖像，是由紅色、綠色、藍色三個通道組成的。紅色、綠色、藍色三個通道的縮覽圖都是以灰度顯示的。用不同的灰度色階來表示“ 紅，綠，藍”在圖像中的比重。通道中的純白，代表了該色光在此處為最高亮度，亮度級別是255。

通道是整個Photoshop顯示圖像的基礎(chǔ)。色彩的變動，實際上就是間接在對通道灰度圖進行調(diào)整。通道是Photoshop處理圖像的核心部分，所有的色彩調(diào)整工具都是圍繞在這個核心周圍使用的。

在計算機領(lǐng)域中，灰度數(shù)字圖像是每個像素只有一個采樣顏色的圖像。這類圖像通常顯示為從最暗黑色到最亮的白色的灰度，盡管理論上這個采樣可以任何顏色的不同深淺，甚至可以是不同亮度上的不同顏色。灰度圖像與黑白圖像不同，在計算機圖像領(lǐng)域中黑白圖像只有黑色與白色兩種顏色；灰度圖像在黑色與白色之間還有許多級的顏色深度。但是，在數(shù)字圖像領(lǐng)域之外，“黑白圖像”也表示“灰度圖像”，例如灰度的照片通常叫做“黑白照片”。在一些關(guān)于數(shù)字圖像的文章中單色圖像等同于灰度圖像，在另外一些文章中又等同于黑白圖像。灰度圖像經(jīng)常是在單個電磁波頻譜如可見光內(nèi)測量每個像素的亮度得到的。

用于顯示的灰度圖像通常用每個采樣像素 8 位的非線性尺度來保存，這樣可以有 256 級灰度。這種精度剛剛能夠避免可見的條帶失真，并且非常易于編程。在醫(yī)學(xué)圖像與遙感圖像這些技術(shù)應(yīng)用中經(jīng)常采用更多的級數(shù)以充分利用每個采樣 10 或 12 位的傳感器精度，并且避免計算時的近似誤差。在這樣的應(yīng)用領(lǐng)域每個采樣 16 位即 65536 級得到流行。

二值圖像：

是指每個像素不是黑就是白，其灰度值沒有中間過渡的圖像。二值圖像一般用來描述文字或者圖形，其優(yōu)點是占用空間少，缺點是，當(dāng)表示人物，風(fēng)景的圖像時，二值圖像只能描述其輪廓，不能描述細節(jié)。這時候要用更高的灰度級。

二值圖像是每個像素只有兩個可能值的數(shù)字圖像。人們經(jīng)常用黑白、B&W、單色圖像表示二值圖像，但是也可以用來表示每個像素只有一個采樣值的任何圖像，例如灰度圖像等。

二值圖像中所有的像素只能從0和1這兩個值中取，因此在MATLAB中，二值圖像用一個由0和1組成的二維矩陣表示。這兩個可取的值分別對應(yīng)于關(guān)閉和打開，關(guān)閉表征該像素處于背景，而打開表征該像素處于前景。以這種方式來操作圖像可以更容易識別出圖像的結(jié)構(gòu)特征。二值圖像操作只返回與二值圖像的形式或結(jié)構(gòu)有關(guān)的信息，如果希望對其他類型的圖像進行同樣的操作，則首先要將其轉(zhuǎn)換為二進制的圖像格式，可以通過調(diào)用MATLAB提供的im2bw（）來實現(xiàn)。

二值圖像經(jīng)常出現(xiàn)在數(shù)字圖像處理中作為圖像掩碼或者在圖像分割、二值化和dithering的結(jié)果中出現(xiàn)。一些輸入輸出設(shè)備，如激光打印機、傳真機、單色計算機顯示器等都可以處理二值圖像。

二值圖像經(jīng)常使用位圖格式存儲。

二值圖像可以解釋為二維整數(shù)格Z，圖像變形處理領(lǐng)域很大程度上就是受到這個觀點啟發(fā)。

圖像旋轉(zhuǎn)：

圖像旋轉(zhuǎn)是指圖像以某一點為中心旋轉(zhuǎn)一定的角度,形成一幅新的圖像的過程。當(dāng)然這個點通常就是圖像的中心。既然是按照中心旋轉(zhuǎn),自然會有這樣一個屬性：旋轉(zhuǎn)前和旋轉(zhuǎn)后的點離中心的位置不變.根據(jù)這個屬性,我們可以得到旋轉(zhuǎn)后的點的坐標(biāo)與原坐標(biāo)的對應(yīng)關(guān)系。由于原圖像的坐標(biāo)是以左上角為原點的,所以我們先把坐標(biāo)轉(zhuǎn)換為以圖像中心為原點。假設(shè)原圖像的寬為w,高為h,（x0,y0）為原坐標(biāo)內(nèi)的一點,轉(zhuǎn)換坐標(biāo)后的點為（x1,y1）。那么不難得到： x1 = x0-w/2;y1 =-y0 + h/2;在新的坐標(biāo)系下,假設(shè)（x0,y0）距離原點的距離為r,點與原點之間的連線與x軸的夾角為b,旋轉(zhuǎn)的角度為a,旋轉(zhuǎn)后的點為（x1,y1）

噪聲：

是電路或系統(tǒng)中不含信息量的電壓或電流。在工業(yè)與自然界中，存在著各種干擾源（噪聲源），如大功率電力電子器件的接入、大功率用電設(shè)備的開啟與斷開、雷擊閃電等都會使空間電場和磁場產(chǎn)生有序或無序的變化，這些都是干擾源（或噪聲源）。這些源產(chǎn)生的電磁波或尖峰脈沖通過磁、電耦合或是通過電源線等路徑進入放大電路，各種電氣設(shè)備，形成各種形式的干擾。

斑點噪聲：

斑點噪聲是SAR成像系統(tǒng)的一大特色，源自基本分辨單元內(nèi)地物的隨機散射，在圖像上表現(xiàn)為信號相關(guān)（如在空間上相關(guān)）的小斑點，它既降低了圖像的畫面質(zhì)量，又嚴(yán)重影響圖像的自動分割、分類、目標(biāo)檢測以及其它定量專題信息的提取。

SAR圖像斑點噪聲的去除一方面要抑制圖像均勻區(qū)域斑點噪聲，另一方面要保持圖像邊緣和紋理細節(jié)信息。SAR斑點噪聲的抑制可通過非相干多視處理，也可使用空間域濾波實現(xiàn)。非相干多視處理會降低圖像的地面分辨率。因此，涌現(xiàn)出了一系列空間域濾波方法，如均值濾波、中值濾波、Lee濾波、Kuan濾波、Frost濾波、Sigma濾波以及Gamma Map濾波等。但這類算法存在自身無法克服的矛盾：一方面為增強斑點去噪效果需選較大的濾波窗口，另一方面為保持圖像的實際分辨率要求所選的窗口較小。

高斯噪聲：

所謂高斯噪聲是指它的概率密度函數(shù)服從高斯分布（即正態(tài)分布）的一類噪聲。如果一個噪聲，它的幅度分布服從高斯分布，而它的功率譜密度又是均勻分布的，則稱它為高斯白噪聲。高斯白噪聲的二階矩不相關(guān)，一階矩為常數(shù)，是指先后信號在時間上的相關(guān)性。高斯白噪聲包括熱噪聲和散粒噪聲。

實驗中是通過MATLAB自帶的函數(shù)產(chǎn)生噪聲，各函數(shù)如下： J1=imnoise(I,'salt & pepper',0.05);%添加椒鹽噪聲

J2=imnoise(I,'gaussian',0,0.03);

%添加均值為0，方差為0.03的高斯噪聲。

六、Matlab優(yōu)勢

MATLAB是一個包含大量算法的集合。其可以快捷的實現(xiàn)用戶所需的各種計算功能。函數(shù)中所使用的算法都是科研和工程計算中的最新研究成果，而前經(jīng)過了各種優(yōu)化和差錯處理。在通常情況下，可以用它來代替底層編程語言，如C和C++。在計算要求相同的情況下，使用MATLAB的編程工作量會大大減少。MATLAB的這些函數(shù)集包括從最簡單最基本的函數(shù)到諸如矩陣，特征向量、快速傅立葉變換的復(fù)雜函數(shù)。函數(shù)所能解決的問題其大致包括矩陣運算和線性方程組的求解、微分方程及偏微分方程的組的求解、符號運算、傅立葉變換和數(shù)據(jù)的統(tǒng)計分析、工程中的優(yōu)化問題、稀疏矩陣運算、復(fù)數(shù)的各種運算、三角函數(shù)和其他初等數(shù)學(xué)運算、多維數(shù)組操作以及建模動態(tài)仿真等。圖形處理功能

圖形處理功能MATLAB自產(chǎn)生之日起就具有方便的數(shù)據(jù)可視化功能，以將向量和矩陣用圖形表現(xiàn)出來，并且可以對圖形進行標(biāo)注和打印。高層次的作圖包括二維和三維的可視化、圖象處理、動畫和表達式作圖。可用于科學(xué)計算和工程繪圖。新版本的MATLAB對整個圖形處理功能作了很大的改進和完善，使它不僅在一般數(shù)據(jù)可視化軟件都具有的功能（例如二維曲線和三維曲面的繪制和處理等）方面更加完善，而且對于一些其他軟件所沒有的功能（例如圖形的光照處理、色度處理以及四維數(shù)據(jù)的表現(xiàn)等），MATLAB同樣表現(xiàn)了出色的處理能力。同時對一些特殊的可視化要求，例如圖形對話等，MATLAB也有相應(yīng)的功能函數(shù)，保證了用戶不同層次的要求。另外新版本的MATLAB還著重在圖形用戶界面（GUI）的制作上作了很大的改善，對這方面有特殊要求的用戶也可以得到滿足

模塊集合工具箱

MATLAB對許多專門的領(lǐng)域都開發(fā)了功能強大的模塊集和工具箱。一般來說，它們都是由特定領(lǐng)域的專家開發(fā)的，用戶可以直接使用工具箱學(xué)習(xí)、應(yīng)用和評估不同的方法而不需要自己編寫代碼。目前，MATLAB已經(jīng)把工具箱延伸到了科學(xué)研究和工程應(yīng)用的諸多領(lǐng)域，諸如數(shù)據(jù)采集、數(shù)據(jù)庫接口、概率統(tǒng)計、樣條擬合、優(yōu)化算法、偏微分方程求解、神經(jīng)網(wǎng)絡(luò)、小波分析、信號處理、圖像處理、系統(tǒng)辨識、控制系統(tǒng)設(shè)計、LMI控制、魯棒控制、模型預(yù)測、模糊邏輯、金融分析、地圖工具、非線性控制設(shè)計、實時快速原型及半物理仿真、嵌入式系統(tǒng)開發(fā)、定點仿真、DSP與通訊、電力系統(tǒng)仿真等，都在工具箱（Toolbox）家族中有了自己的一席之地。

七、總結(jié)

運用matlab軟件對圖像進行處理，讓我鞏固了之前所學(xué)的知識，同時也在這次作業(yè)中更加了解到matlab語言在生活中的運用環(huán)境和掌握這門語言的重要性

八、存在問題

1.在進行圖像增強時要不要講圖像先進行平滑處理？ 2.如何增加這個算法的準(zhǔn)確度

3.在此次作業(yè)中，為何添加椒鹽噪聲時無法顯示

第五篇：數(shù)字圖像處理實驗報告

數(shù)字圖像處理

實驗報告

班級:通信103 學(xué)號:201027201 姓名:計富威 指導(dǎo)教師:孫潔

實驗一 MATLAB數(shù)字圖像處理初步

一、實驗?zāi)康呐c要求

1．熟悉及掌握在MATLAB中能夠處理哪些格式圖像。2．熟練掌握在MATLAB中如何讀取圖像。

3．掌握如何利用MATLAB來獲取圖像的大小、顏色、高度、寬度等等相關(guān)信息。

4．掌握如何在MATLAB中按照指定要求存儲一幅圖像的方法。5．圖像間如何轉(zhuǎn)化。

二、實驗內(nèi)容及步驟

1．利用imread()函數(shù)讀取一幅圖像，假設(shè)其名為”第一個.tif”，存入一個數(shù)組中；

>>I=imread('第一個.tif');2．利用whos命令提取該讀入圖像”第一個.tif”的基本信息； >>whos I 3．利用imshow()函數(shù)來顯示這幅圖像； >>imshow(I);

第一個.tif 4．利用imfinfo函數(shù)來獲取圖像文件的壓縮，顏色等等其他的詳細信息；

>>imfinfo('第一個.tif');5．利用imwrite()函數(shù)來壓縮這幅圖象，將其保存為一幅壓縮了像素的jpg文件,設(shè)為flower.jpg；語法：imwrite(原圖像，新圖像，‘quality’,q), q取0-100。

>>imwrite(I,'第一個.jpg','quality',50)6．同樣利用imwrite()函數(shù)將最初讀入的tif圖象另存為一幅bmp圖像，設(shè)為flower.bmp。>>imwrite(I,'第一個.bmp');7．用imread()讀入圖像：Lenna256.jpg 和camemaman.jpg； >>b=imread('lena256.bmp');>>c=imread('cameraman.tif');8．用imfinfo()獲取圖像Lenna256.jpg和camemaman.jpg 的大小； >>imfinfo('lena256.bmp');>>imfinfo('cameraman.tif');9．用figure,imshow()分別將Lenna256.jpg和camemaman.jpg顯示出來，觀察兩幅圖像的質(zhì)量。>>figure >>imshow(b);>>figure >>imshow(c);

（Lena256.jpg圖像截圖）

（cameraman.jpg圖像截圖）

10．用im2bw將一幅灰度圖像轉(zhuǎn)化為二值圖像，并且用imshow顯示出來觀察圖像的特征。>> d=im2bw(b);>>figure >>imshow(b);>>figure >>imshow(d);

（二值化截圖）

三、實驗總結(jié)

通過實驗MatLab軟件的基本使用有了基本的了解，學(xué)會了使用MatLab軟件來讀取一個特定格式的圖像，并通過相關(guān)的命令語句對圖像進行格式轉(zhuǎn)換、圖像壓縮、二值化等的處理，掌握了利用MATLAB來獲取圖像的大小、顏色、高度、寬度等等相關(guān)信息，掌握在MATLAB中如何通過imshow()語句來讀取圖像等等。

第二 圖像基本運算

一、實驗?zāi)康?/p>

1．了解圖像的算術(shù)運算在數(shù)字圖像處理中的初步應(yīng)用。2．體會圖像算術(shù)運算處理的過程和處理前后圖像的變化。

二、實驗原理

圖像的代數(shù)運算是圖像的標(biāo)準(zhǔn)算術(shù)操作的實現(xiàn)方法，是兩幅輸入圖像之間進行的點對點的加、減、乘、除運算后得到輸出圖像的過程。如果輸入圖像為A(x,y)和B(x,y)，輸出圖像為C(x,y)，則圖像的代數(shù)運算有如下四種形式：

C(x,y)= A(x,y)+ B(x,y)C(x,y)= A(x,y)-B(x,y)C(x,y)= A(x,y)* B(x,y)C(x,y)= A(x,y)/ B(x,y)

三、實驗步驟 1．圖像的加法運算

在MATLAB中，如果要進行兩幅圖像的加法，或者給一幅圖像加上一個常數(shù)，可以調(diào)用imadd函數(shù)來實現(xiàn)。imadd函數(shù)將某一幅輸入圖像的每一個像素值與另一幅圖像相應(yīng)的像素值相加，返回相應(yīng)的像素值之和作為輸出圖像。imadd函數(shù)的調(diào)用格式如下：

Z = imadd（X，Y）首先讀入兩幅圖像

>>a=imread('第二個原圖1.jpg');>>b=imread('第二個原圖2.jpg')通過一個加法操作：>> c=imadd(a,b);

給圖像的每一個像素加上一個常數(shù)可以使圖像的亮度增加。如截圖

第一張為原圖，第二張為亮度加50，第三張為亮度減50 2．圖像的減法運算

在MATLAB中，使用imsubtract函數(shù)可以將一幅圖像從另一幅圖像中減去，或者從一幅圖像中減去一個常數(shù)。imsubtract函數(shù)將一幅輸入圖像的像素值從另一幅輸入圖像相應(yīng)的像素值中減去，再將這個結(jié)果作為輸出圖像相應(yīng)的像素值。imsubtract函數(shù)的調(diào)用格式如下：

Z = imsubtract(X,Y)； 讀入一幅畫后通過減法 >>a3=imsubtract(a,50);

3.圖像的乘法運算

在MATLAB中，使用immultiply函數(shù)實現(xiàn)兩幅圖像的乘法。immultiply函數(shù)將兩幅圖像相應(yīng)的像素值進行元素對元素的乘法操作（MATLAB點乘），并將乘法的運算結(jié)果作為輸出圖形相應(yīng)的像素值。immulitply函數(shù)的調(diào)用格式如下：

Z = immulitply(X,Y)讀入一幅圖后通過乘法操作 >> a=imread('cameraman.tif');>> b=immultiply(a,1.5);

4．圖像的除法運算

在MATLAB中使用imdivide函數(shù)進行兩幅圖像的除法。imdivide函數(shù)對兩幅輸入圖像的所有相應(yīng)像素執(zhí)行元素對元素的除法操作（點除），并將得到的結(jié)果作為輸出圖像的相應(yīng)像素值。imdivide函數(shù)的調(diào)用格式如下：

Z = imdivide(X,Y)讀入一幅圖后通過除法操作

四、實驗總結(jié)

通過對圖像的四則運算了結(jié)圖像的不同變化過程，對軟件的進一步使用也有了更加深刻的認(rèn)識。

實驗三 圖像增強—空域濾波

一、實驗?zāi)康?/p>

進一步了解MatLab軟件/語言，學(xué)會使用MatLab對圖像作濾波處理，使學(xué)生有機會掌握濾波算法，體會濾波效果。

了解幾種不同濾波方式的使用和使用的場合，培養(yǎng)處理實際圖像的能力，并為課堂教學(xué)提供配套的實踐機會。

二、實驗設(shè)備與軟件

(1)IBM-PC計算機系統(tǒng)

(2)MatLab軟件/語言包括圖像處理工具箱(Image Processing Toolbox)

(3)實驗所需要的圖片

三、實驗內(nèi)容與步驟

a)調(diào)入并顯示原始圖像“原圖像.jpg”。>>I=imread('原圖像.jpg');b)利用imnoise命令在圖像“原圖像.jpg”上加入高斯(gaussian)噪聲

>>J = imnoise(I,'gauss',0.02);

%添加高斯噪聲 c)利用預(yù)定義函數(shù)fspecial命令產(chǎn)生平均(average)濾波器 d）分別采用3x3和5x5的模板，分別用平均濾波器以及中值濾波器，對加入噪聲的圖像進行處理并觀察不同噪聲水平下，上述濾波器處理的結(jié)果；

>>ave1=fspecial('average',3);

%產(chǎn)生3×3的均值模版 >>ave2=fspecial('average',5);

%產(chǎn)生5×5的均值模版 >>K = filter2(ave1,J)/255;

%均值濾波3×3 >>L = filter2(ave2,J)/255;

%均值濾波5×5 e）選擇不同大小的模板，對加入某一固定噪聲水平噪聲的圖像進行處理，觀察上述濾波器處理的結(jié)果。

>>M = medfilt2(J,[3 3]);

%中值濾波3×3模板 >>N = medfilt2(J,[4 4]);

%中值濾波4×4模板

f）利用imnoise命令在圖像Sample2-1.jpg 上加入椒鹽噪聲(salt & pepper)>>J = imnoise(I,'salt& pepper',0.02);

%添加椒鹽噪聲

四、實驗總結(jié)

椒鹽噪聲是由圖像傳感器，傳輸信道，解碼處理等產(chǎn)生的黑白相間的亮暗點噪聲。椒鹽噪聲往往由圖像切割引起，去除脈沖干擾及椒鹽噪聲最常用的算法是中值濾波。椒鹽噪聲是指兩種噪聲，一種是鹽噪聲，另一種是胡椒噪聲。鹽=白色，椒=黑色。前者是高灰度噪聲，后者屬于低灰度噪聲。一般兩種噪聲同時出現(xiàn)，呈現(xiàn)在圖像上就是黑白雜點。這點我們通過實驗結(jié)果可以明顯看到。中值濾波對于濾除圖像的椒鹽噪聲非常有效。

實驗四圖像分割

一、實驗?zāi)康?/p>

使用MatLab 軟件進行圖像的分割。使學(xué)生通過實驗體會一些主要的分割算子對圖像處理的效果，以及各種因素對分割效果的影響。

二、實驗要求

要求學(xué)生能夠自行評價各主要算子在無噪聲條件下和噪聲條件下的分割性能。能夠掌握分割條件(閾值等)的選擇。完成規(guī)定圖像的處理并要求正確評價處理結(jié)果，能夠從理論上作出合理的解釋。

三、實驗內(nèi)容與步驟

(1)使用Roberts 算子的圖像分割實驗，使用的原圖是cameraman.jpg,截圖如下

(2)使用Prewitt 算子的圖像分割實驗 截圖如下

(3)使用Sobel 算子的圖像分割實驗

(4)使用LoG(拉普拉斯-高斯)算子的圖像分割實驗

四、實驗結(jié)果

對Roberts算子、Prewitt 算子、Sobel 算子、LoG(拉普拉斯-高斯)算子的運算對圖像的結(jié)果有了基本的認(rèn)識，加深學(xué)習(xí)效果。

實驗五 形態(tài)學(xué)運算

1、實驗?zāi)康?/p>

學(xué)習(xí)常見的數(shù)學(xué)形態(tài)學(xué)運算基本方法，了解腐蝕、膨脹、開運算、閉運算取得的效果，培養(yǎng)處理實際圖像的能力，并為課堂教學(xué)提供配套的實踐機會。

2、實驗要求

利用MatLab工具箱中關(guān)于數(shù)學(xué)形態(tài)學(xué)運算的函數(shù)，計算本指導(dǎo)書中指定二值圖像進行處理。

3、實驗設(shè)備與軟件

1.LC-PC計算機系統(tǒng)

2.MatLab軟件/語言包括圖像處理工具箱(Image Processing Toolbox)3.實驗所需要的圖片

4、實驗內(nèi)容與步驟

1.調(diào)入并顯示圖像“原圖.jpg”； 2.調(diào)入并顯示圖像“原圖.jpg”；

3.選取合適的閾值，得到二值化圖像“原圖.jpg”； >>bw = im2bw(I,level);

%二值化 4.設(shè)置結(jié)構(gòu)元素；

5.對得到的二值圖像“原圖.jpg”進行腐蝕運算； >>BW2 = imerode(bw,SE1);

%腐蝕 6.對得到的二值圖像“原圖.jpg”進行膨脹運算； >>BW1 = imdilate(bw,SE);

%膨脹 7.對得到的二值圖像“原圖.jpg”進行開運算；

>>BW3 = bwmorph(bw, 'open');

%開運算 8.對得到的二值圖像“原圖.jpg”進行閉運算； >>BW4 = bwmorph(bw, 'close');

%閉運算 9.將兩種處理方法的結(jié)果作比較；

五、實驗總結(jié)

通過本次實驗，學(xué)習(xí)了常見的數(shù)學(xué)形態(tài)學(xué)運算基本方法，了解腐蝕、膨脹、開運算、閉運算取得的效果，培養(yǎng)處理實際圖像的能力，通過自己動手的實驗，對課本上的知識有了更加深刻的理解。