【正文】 eshgrid(x,y)。% 計(jì)算傳遞函數(shù) end result(i,j)=h*g(i,j)。n=fix(N/2)。)。subplot(3,2,3)。)。z=log(abs(result))。end end subplot(3,2,4)。n=fix(N/2)。[M,N]=size(g)。%取幅度 mesh(z)。x=0:1:255。f=double(J1)。imshow(J1)。J1=imread(39。)。%新建圖像顯示窗口 % [n,m]=size(c)。title(39。subplot(3,2,3)。mesh(z)。對(duì)該圖進(jìn)行低通濾波 選取D=190z=log(abs(result))。[n,m]=size(i)。mesh(z)。z=log(abs(i))。)。subplot(3,2,1)。clear。五、思考題分析為什么圖像通過(guò)低通濾波器后變得模糊？為什么通過(guò)高通濾波器后得到銳化結(jié)果？答：圖像的精細(xì)結(jié)構(gòu)及突變部分主要由高頻成分起作用，故經(jīng)低通濾波后圖像的精細(xì)結(jié)構(gòu)消失，變得模糊；經(jīng)高通濾波后圖像得到銳化。)。%顯示均衡化后的圖像 subplot(2,2,4)。imshow(J)。title(39。原圖39。J=histeq(Im)。%計(jì)算圖像的歸一化直方圖 subplot(1,3,1),imshow(I)； subplot(1,3,2),imhist(I,64)； subplot(1,3,3),plot(J)；（2）直方圖均衡化clear all。)。六、實(shí)驗(yàn)心得體會(huì)本次實(shí)驗(yàn)中，因?yàn)槌鯇W(xué)這個(gè)軟件，我學(xué)習(xí)到了在程序中關(guān)于圖像的運(yùn)用，以及也復(fù)習(xí)了課本上的許多知識(shí)，加深了對(duì)直方圖均衡化的理解。子圖直方圖之和為整圖的直方圖。五、思考題（1）、灰度直方圖可以反映一幅圖像的哪些特征？ 答：表征了圖像的一維信息。當(dāng)時(shí)，此時(shí)維納濾波就變成 了逆濾波；當(dāng)時(shí)，表明維納濾波避免了逆濾波中 出現(xiàn)的對(duì)噪聲過(guò)多的放大作用,也就是說(shuō)圖像維納復(fù)原比逆濾波復(fù)原效果好。圖像退化是由于成像系統(tǒng)受各種因素的影響，導(dǎo)致了圖像質(zhì)量的降低，稱(chēng)之為圖像退化。圖像增強(qiáng)和圖像復(fù)原之間有何區(qū)別？圖像增強(qiáng)：利用一定的技術(shù)手段，不用考慮圖像是否失真（即原 始圖像在變換后可能會(huì)失真）而且不用分析圖像降質(zhì)的原因。十一、實(shí)驗(yàn)步驟以圖像中值濾波操作為例說(shuō)明實(shí)驗(yàn)的具體步驟，其它實(shí)驗(yàn)項(xiàng)目的步驟與此類(lèi)似。關(guān)于圖像復(fù)原的詳細(xì)原理可以參考相關(guān)書(shū)籍，這里不再贅述。和194。圖像增強(qiáng)可以用空域法和頻域法分別實(shí)現(xiàn)，空域法主要是在空間域中對(duì)圖像象素灰度值直接進(jìn)行運(yùn)算處理，一般包括中值濾波、模板平滑和梯度銳化等，空域法可以用下式來(lái)描述：g(x,y)＝f(x,y)*h(x,y)其中f(x,y)是處理前圖像，g(x,y)表示處理后圖像，h(x,y)為空間運(yùn)算函數(shù)。進(jìn)行一次鏡像變換，順（逆）時(shí)針旋轉(zhuǎn)兩次，再以與第一次相反的方向鏡像變換?；叶染馐窃鰪?qiáng)圖像的有效方法之一。它的操作和閾值變換類(lèi)似。步驟14和15是使用圖像格式轉(zhuǎn)換工具的方法；1在桌面上雙擊圖像格式轉(zhuǎn)換工具Jpg2bmp的圖標(biāo)，進(jìn)入轉(zhuǎn)換工具界面，如圖所示；1按照界面提示，把JPG格式的圖像文件轉(zhuǎn)換成8位BMP圖像。至于FFT的編程實(shí)現(xiàn)，系統(tǒng)采用的方法是：首先編制一個(gè)一維FFT程序模塊，然后調(diào)用該模塊對(duì)圖像數(shù)據(jù)的列進(jìn)行一維FFT，再對(duì)行進(jìn)行一維FFT，最后計(jì)算并顯示幅度譜。下面以圖像點(diǎn)運(yùn)算中的閾值變換為例給出編程實(shí)現(xiàn)的程序流程圖，如下：數(shù)字圖象處理實(shí)驗(yàn)指導(dǎo)書(shū)實(shí)驗(yàn)一 圖象變換實(shí)驗(yàn)圖像正交變換的基本原理及編程實(shí)現(xiàn)步驟 數(shù)字圖像的處理方法主要有空域法和頻域法，點(diǎn)運(yùn)算和幾何變換屬于空域法。點(diǎn)運(yùn)算又被稱(chēng)為對(duì)比度增強(qiáng)、對(duì)比度拉伸或灰度變換。點(diǎn)運(yùn)算可以看作是“從象素到象素”的復(fù)制操作，而這種復(fù)制操作是通過(guò)灰度變換函數(shù)實(shí)現(xiàn)的。感覺(jué)做matlab的實(shí)驗(yàn)和平時(shí)編寫(xiě)C代碼的時(shí)候有很多地方是相同的，都會(huì)經(jīng)歷：編碼，debug，正確這三個(gè)過(guò)程，這個(gè)過(guò)程本身就具有挑戰(zhàn)性和趣味性。%f的值控制在[0,1]之間%彩色圖像平滑處理fen=imsubtract(Ha,fb)。%對(duì)亮度分量進(jìn)行平滑操作 %將H、S、I三個(gè)分量重新組合起來(lái) H1=cat(3,H,S,I_filtered)。,25)。S=h(:,:,2)。Ha=im2double(Ha)。%對(duì)飽和度分量進(jìn)行平滑操作 I_filtered=imfilter(I,w,39。)。average39。%直接對(duì)原圖像的RGB三個(gè)分量分別進(jìn)行平滑操作 H=f(:,:,1)。程序代碼如下：f=imread(39。map(2,:)=。[x,map]=gray2ind(f1,128)。)。clear。,f,25,m,C)。E251=f。%轉(zhuǎn)換到double類(lèi)型 [C,m]=covmatrix(I)。%顯示取得選定區(qū)域[M,N,K]=size(g)。%取得選定區(qū)域圖像的綠色分量G=immultiply(mask,f(:,:,2))。39。%顯示彩色梯度對(duì)應(yīng)的角度圖像figure,imshow(PPG,[ ])。)。clear。%拉普拉斯算子處理 figure,imshow(g)。1 1 1]。D:PictureFig0604(a)(iris).tif39。%將H、S、I三個(gè)分量重新組合起來(lái)f=hsi2rgb(H1)。I_filtered=imfilter(I,w,39。I=h(:,:,3)。figure,imshow(fc)。%f2的圖像色彩偏紅，拖動(dòng)映射曲線，%調(diào)整映射參數(shù)，使圖像的色彩看起來(lái)比較正常。space39。figure,imshow(f2)。,f)。D:PictureFig0614(a)(Chalk Original).tif39。figure,imshow(hsi_image(:,:,1))。figure,imshow(ntsc_image(:,:,1))。D:PictureFig0604(a)(iris).tif39。%f轉(zhuǎn)換為灰度圖像g1=dither(g)。dither39。nodither39。D:PictureFig0604(a)(iris).tif39。第一篇：數(shù)字圖像處理——彩色圖像實(shí)驗(yàn)報(bào)告(1)對(duì)彩色圖像的表達(dá)和顯示* * * * * * * * * * * *顯示彩色立方體* * * * * * * * * * * * *rgbcube(0,0,10)。%從對(duì)角線觀察彩色立方體%* * * * * * * * * *索引圖像的顯示和轉(zhuǎn)換* * * * * * * * * f=imread(39。%f是RGB真彩圖像%rgb圖像轉(zhuǎn)換成8色索引圖像，不采用抖動(dòng)方式[X1,map1]=rgb2ind(f,8,39。%采用抖動(dòng)方式轉(zhuǎn)換到8色索引圖像[X2,map2]=rgb2ind(f,8,39。%顯示效果要好一些 g=rgb2gray(f)。%顯示抖動(dòng)處理后的二值圖像程序運(yùn)行結(jié)果：*彩色立方體原圖不采用抖動(dòng)方式轉(zhuǎn)換到8色索引圖像采用抖動(dòng)方式轉(zhuǎn)換到8色索引圖像灰度圖像抖動(dòng)處理后的二值圖像(2)彩色空間轉(zhuǎn)換f=imread(39。%轉(zhuǎn)換到NTSC彩色空間%f是RGB真彩圖像ntsc_image=rgb2ntsc(f)。%顯示色差信息%轉(zhuǎn)換到HIS彩色空間hsi_image=rgb2hsi(f)。%顯示亮度信息程序運(yùn)行結(jié)果：原圖轉(zhuǎn)換到NTSC彩色空間顯示亮度信息顯示色差信息顯示色差信息轉(zhuǎn)換到HIS彩色空間顯示色差信息顯示飽和度信顯示亮度信息(3)彩色變換f=imread(39。image39。)。,f2,39。)。)。S=h(:,:,2)。,25)。%對(duì)亮度分量進(jìn)行平滑操作 H1=cat(3,H,S,I_filtered)。%彩色圖像平滑處理 %* * * **** * * * * * * * * * * * * 彩色圖像銳化 * * * * * * * * * * * * * * * * * * * * * *fb=imread(39。18 1。)。%顯示差值圖像程序運(yùn)行結(jié)果：原圖彩色圖像平滑處理拉普拉斯算子處理銳化處理(5)彩色圖像分割%* * * * * * * * * * * * * * * * * 彩色邊緣檢測(cè)* * * * * * * * * * * *close all。39。%顯示彩色梯度圖像 figure,imshow(A,[ ])。程序運(yùn)行結(jié)果：%顯示差值圖像彩色梯度圖像 彩色梯度對(duì)應(yīng)的角度圖像分量梯度圖像兩種梯度差值圖像%* * * * * * * * * * * * * * * * * 彩色圖像分割 * * * * * * * * * * * * * * * * * * * * * * * * *f=imread(39。%取得選定區(qū)域圖像的紅色分量%手工選定圖中紅色區(qū)域R=immultiply(mask,f(:,:,1))。%取得選定區(qū)域的RGB圖像 figure,imshow(g)。%找到mask中的像素在重排后的圖像g中位置 I=double(I(idx,1:3))。,f,25,m)。mahalanobis39。程序代碼如下%思路：把原圖想轉(zhuǎn)換為128位索引圖close all。D:39。%將f轉(zhuǎn)換成灰度圖像 figure,imshow(f1)。map(: ,1)=。程序運(yùn)行結(jié)果：原RGB圖像灰度圖像索引圖 修改索引表后的圖像(2)（4）是將圖像轉(zhuǎn)換到HIS空間后再對(duì)亮度分量進(jìn)行平滑操作，試直接對(duì)原圖像的RGB三個(gè)分量分別進(jìn)行平滑操作，并計(jì)算兩種操作方法的圖像差值。figure,imshow(f)。%色度分量%飽和度分量%亮度分量%2525的方形平滑算子 w=fspecial(39。replicate39。)。%對(duì)亮度分量進(jìn)行平滑操作 %將H、S、I三個(gè)分量重新組合起來(lái)Ha=cat(3,H_filtered,S_filtered,I_filtered)。%轉(zhuǎn)換到HIS彩色空間H=h(:,:,1)。average39。)。figure,imshow(fb)。(1)通過(guò)這個(gè)實(shí)驗(yàn)我了解了RGB圖像以及索引圖像在MATLAB中的存儲(chǔ)和詢(xún)問(wèn)方式(2)了解RGB空間到其他彩色空間的轉(zhuǎn)換方式和轉(zhuǎn)換方法(3)知道了一些彩色圖像的空間濾波方法(4)了解在RGB向量空間中檢測(cè)彩色邊緣的方法(5)了解在RGB向量空間中對(duì)圖像進(jìn)行分割的方法遇到的問(wèn)題：在做思考題(2)的時(shí)候，計(jì)算兩種平滑處理的得到圖像的差值的中程序出現(xiàn)了問(wèn)題，經(jīng)過(guò)看出錯(cuò)信息，認(rèn)真檢查程序后，發(fā)現(xiàn)是因?yàn)橹苯悠交幚淼膱D像是unit8類(lèi)型的，而轉(zhuǎn)換到HIS空間后，得到的圖像是double類(lèi)型，二者的類(lèi)型不同，導(dǎo)致使用imsubtract函數(shù)的時(shí)候出現(xiàn)了問(wèn)題。圖像點(diǎn)運(yùn)算是一種簡(jiǎn)單而重要的處理技術(shù)，它能讓用戶(hù)改變圖像數(shù)據(jù)占據(jù)的灰度范圍。另外，點(diǎn)運(yùn)算處理將改變圖像的灰度直方圖分布。實(shí)驗(yàn)系統(tǒng)提供了圖像灰度直方圖、點(diǎn)運(yùn)算和幾何變換相關(guān)內(nèi)容的文字說(shuō)明，用戶(hù)在操作過(guò)程中可以參考。正交變換實(shí)驗(yàn)的重點(diǎn)是快速傅立葉變換（FFT），其原理過(guò)于復(fù)雜，可以參考有關(guān)書(shū)籍，這里不再贅述。打開(kāi)計(jì)算機(jī)，在系統(tǒng)桌面上雙擊“數(shù)字圖像處理實(shí)驗(yàn)教學(xué)軟件系統(tǒng)”的可執(zhí)行文件“圖象處理”的圖標(biāo)，進(jìn)入實(shí)驗(yàn)系統(tǒng)；執(zhí)行文件→打開(kāi)，在OPEN對(duì)話(huà)框中選擇待處理的圖像，按【OK】后系統(tǒng)顯示出圖像；執(zhí)行查看→圖像基本信息，將顯示圖像基本信息對(duì)話(huà)框，如圖所示；數(shù)字圖象處理實(shí)驗(yàn)指導(dǎo)書(shū)實(shí)驗(yàn)一 圖象變換實(shí)驗(yàn)執(zhí)行查看→灰度直方圖，查看圖像的灰度直方圖，如圖所示；執(zhí)行圖像變換→正交變換→傅立葉變換，查看圖像的頻率域分布情況，如圖所示；數(shù)字