【正文】 原圖像 (2) 新均值去噪圖像 (3)均值去噪圖像我們將以上圖像對(duì)比后，能觀察出新均值去噪法的去噪效果明顯比均值去噪法的去噪效果要好，而且新均值去噪法對(duì)圖像的模糊化程度要比均值去噪要小很多，所以新均值去噪法對(duì)汽車縱梁圖像的去噪還是很有效的。我使用的是迭代法閾值選取法，具體方法和步驟如下：1)選擇一個(gè)初始閾值T，例如，假定圖像像素的最大灰度值和最小灰度值分別為和，則初始閾值T可以選擇為。3)分別計(jì)算G1和G2包含的像素的灰度均值u1和u2 。5)重復(fù)步驟2)、3)、4)，直到連續(xù)兩次計(jì)算得到的T的差值滿足設(shè)定的范圍，從而完成閾值的自動(dòng)計(jì)算。 tmax=max(I(:))。 ok=true。 g2=Ith。 u2=mean(I(g2))。 ok=abs(ththnew)=1。end。 閾值計(jì)算完成以后，以閾值的百分之五為基準(zhǔn)，來確定圖像上的點(diǎn)的灰度最后取源圖像上的還是去噪后圖像的。 新均值去噪法的時(shí)間復(fù)雜度我們要評(píng)價(jià)一種方法的好壞，不僅要考慮處理結(jié)果的好壞，還要考慮該方法所消耗的成本。我們可以很容易的得到均值去噪的時(shí)間復(fù)雜度為n2,而新均值去噪法與均值去噪法的區(qū)別在于新思想部分，我們根據(jù)以上程序段可以分析到，新均值去法的時(shí)間復(fù)雜度也是n2,所以新均值去噪與均值去噪花費(fèi)的時(shí)間在一個(gè)級(jí)別，兩者消耗成本基本相同。在本次課題研究中我是利用Matlab的GUI模塊進(jìn)行設(shè)計(jì)的界面。 圖像輸入模塊，它是由一個(gè)Push Button模塊與Axes圖像顯示模塊組成。 加入噪聲模塊與去噪模塊，它們都是由一個(gè)Popup Menu模塊與Axes圖像顯示模塊組成。它們的使用是先將原圖片經(jīng)過加入噪聲模塊選擇一種噪聲加入，然后用去噪模塊中的各種去噪方法進(jìn)行去噪處理得出的去噪圖像顯示在去噪模塊中。 對(duì)汽車縱梁圖像噪聲處理模塊，它是由一個(gè)Popup Menu模塊與Axes圖像顯示模塊組成。 退出模塊 這個(gè)模塊比較簡(jiǎn)單，就是一個(gè)Push Button模塊構(gòu)成， 退出模塊這個(gè)模塊的作用主要是使系統(tǒng)安全退出，盡量避免非法退出，因?yàn)槟菢尤菀壮霈F(xiàn)一些意想不到的錯(cuò)誤。 界面總圖結(jié)束語經(jīng)過大四下學(xué)期的努力，我從最初的資料收集，對(duì)各種去噪算法的研究，設(shè)計(jì)最佳方法，界面設(shè)計(jì)，程序編寫到最后的調(diào)試，終于達(dá)到了本課題的要求，實(shí)現(xiàn)了對(duì)汽車縱梁圖像很好的去噪處理，能夠很好的將其噪聲去掉并且對(duì)圖像中的有用信息保留的最多。在研究的過程中尋找各種算法的優(yōu)缺點(diǎn)，然后經(jīng)過各種嘗試，最后在均值濾波的基礎(chǔ)上，將其改進(jìn)成為適合于汽車縱梁圖像噪聲處理的新方法，即新均值去噪法，能夠很好的對(duì)圖片進(jìn)行去噪處理。在設(shè)計(jì)過程中，我是以Matlab作為平臺(tái)的，先設(shè)計(jì)了一個(gè)界面，然后對(duì)界面中的各個(gè)模塊進(jìn)行程序編寫并且激活，這樣能夠很好的將我的畢業(yè)設(shè)計(jì)成果很好的展現(xiàn)給大家，就這樣我完成了整個(gè)課題研究。將新均值去噪與均值去噪進(jìn)行比較，它能在很好的處理各種噪聲的同時(shí)對(duì)圖像模糊程度降低，這樣也就達(dá)到了這個(gè)課題的研究意義。雖然我研究出來的新均值去噪法有很多優(yōu)點(diǎn)，但是我覺得它還有一點(diǎn)不足之處，那就是新的圖片要選定一個(gè)點(diǎn)是應(yīng)該選擇原圖像上的還是均值去噪圖像上的時(shí)，需要將原圖像與均值去噪圖像相減結(jié)果的絕對(duì)值與原圖像的閾值的百分之五進(jìn)行比較。當(dāng)這個(gè)原圖像閾值的百分比更適合時(shí)，我們處理出來的圖像的模糊度能降到更低，而且在這個(gè)同時(shí)圖像噪聲的處理效果也更好。致謝在畢業(yè)設(shè)計(jì)完成之際，感謝湖北汽車工業(yè)學(xué)院所有老師對(duì)我的關(guān)心和幫助。感謝指導(dǎo)老師周奎的悉心指導(dǎo)。在我的論文方面，周老師又盡可能的提出了寶貴的意見，在論文即將完成之際，工作繁忙的他也依然抽出寶貴的時(shí)間，給予悉心指導(dǎo)，提出很多修改意見，從而進(jìn)一步增強(qiáng)了本文的可讀性和嚴(yán)密性。同時(shí)我還要感謝一直關(guān)心和幫助我的同學(xué)，沒有你們的幫助我的課題也不會(huì)如此順利地完成。湖北汽車工業(yè)學(xué)院學(xué)報(bào)，.[3] 李任江,[J].長春工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2007年9月第28卷第3期 .[4] 李月晶,[J].長春理工大學(xué)學(xué)報(bào)(自然科學(xué)版) ,2008年6月第31卷第2期.[5] 肖生發(fā)，趙樹鵬．汽車構(gòu)造[M].北京：中國林業(yè)出版社；北京大學(xué)出版社，20068[6] 朱菊華，[J].計(jì)算機(jī)工程應(yīng)用，1001(3): 9395[7] 梁雯，[J].安徽大學(xué)學(xué)報(bào)(自然科學(xué)版)，1997(2):4247[8] 劉真，[J].北京廣播學(xué)院學(xué)報(bào)(自然科學(xué))1996(2):28一37[9] [M].南京:東南大學(xué)出版社，1999.[10] [M].北京：機(jī)械工業(yè)出版社，2010：115152[11] [M].北京：北京電子郵電大學(xué)出版社，2006：97106[12] [M].北京：中國電力出版社，：96140[13] 陸玲，王蕾，[M].北京：中國電力出版社，2007:5195[14] Rafael C. Gonzalez著, （MATLAB版）[M] .北京：電子工業(yè)出版社，2005.[15] (第二版)[M].西安：西安電子科技大學(xué)出版社，2008. [16] 陳楊，陳榮娟，郭穎輝．[M].西安：西安電子科技大學(xué)出版社，200210．[17] 姚敏．?dāng)?shù)字圖像處理[M].北京：機(jī)械工業(yè)出版社，20061．[18] Weiwei 基于DSP圖像處理的MATLAB輔助研究[J].河北：廊坊師范學(xué)院學(xué)報(bào),20065.[19] [J]。gui_State = struct(39。, mfilename, ... 39。, gui_Singleton, ... 39。, bishe_OpeningF, ... 39。, bishe_OutputF, ... 39。, [] , ... 39。, [])。amp。endif nargout [varargout{1:nargout}] = gui_mainf(gui_State, varargin{:})。endfunction bishe_OpeningF(hObject, eventdata, handles, varargin) = hObject。 function varargout = bishe_OutputF(hObject, eventdata, handles) varargout{1} = 。*.tif39。*.jpg39。*.bmp39。*.gif39。*.*39。選擇圖片39。if pathname == 0 return。global I。axes()。 function a_Callback(hObject, eventdata, handles)var=get(,39。)。axes()。switch var case 1 j=imnoise(I,39。)。 case 2 j=imnoise(I,39。)。 case 3 j=imnoise(I,39。 pepper39。imshow(j)。speckle39。imshow(j)。amp。BackgroundColor39。defaultUicontrolBackgroundColor39。BackgroundColor39。white39。end function b_Callback(hObject, eventdata, handles)var1=get(,39。)。axes()。if(var1==1)p=wiener2(j,[5 5])。 elseif(var1==2) p=medfilt2(j)。 elseif(var1==3) p=ordfilt2(j,5,ones(3,4))。 elseif(var1==4) j1=im2bw(j)。 j3=bwmorph(j2,39。)。close39。 imshow(p)。average39。 imshow(p)。 function b_CreateF(hObject, eventdata, handles)if ispc amp。 isequal(get(hObject,39。), get(0,39。)) set(hObject,39。,39。)。 function c_Callback(hObject, eventdata, handles)function s_CreateF(hObject, eventdata, handles)if ispc amp。 isequal(get(hObject,39。), get(0,39。)) set(hObject,39。,39。)。value39。 global I。 if(var2==1) q=wiener2(I,[5 5])。 elseif(var2==2) q=medfilt2(I)。 elseif(var2==3) q1=im2bw(I)。 q3=bwmorph(q2,39。)。close39。 imshow(q)。 imshow(q)。average39。 imshow(q)。 q1=filter2(fspecial(39。,3),I1)/255。 tmin=min(I(:))。 th=(tmin+tmax)/2。 while ok g1=I=th。 u1=mean(I(g1))。 thnew=(u1+u2)/2。 th=thnew。 th=floor(th)。 G=[]。 else G(s,t)=I(s,t)。