【正文】 1991,752～756[11] Scheunders P. Multiscale edge representation applied to image fusion [J] .2002[12] Akerman A. Pyramidal techniques for multisensor fusion [C]. SPIE,1992,18（28）,124～131[13] Brown LG. A survey of image registration techniques [J]. Assoc Comput Mach,1992,24,325～376[14] Vajdic ., Katz ., etc. AIbased relational matching and multimodal medical image fusion: Generalized 3D approaches [J]. SPIE,1994, 1565～1573[15] Li H., Dellerck R., De Cuyper B., etc. Object recognition in brain CTscans: Knowledgebased fusion of data from multiple feature extractors [J]. IEEE Trans Med Image . 1995,14（2）,212～229[16] Toet A, Ruyven J .Merging thermal and visual images by a contract pyramid [J]. Optical Engineering,1989,28（7）,789～792[17] Whittington G., Spracklen T. The application of neural network model to sensor date fusion. Proc SPIE, 1990, 276287[18] Daubechies I. Orthonomal bases of pactly supported wavelets. Communication on Pure and Applied Mathemathics,1988,41:909～996[19] 唐向宏,李齊良. 時頻分析與小波變換. 北京:科學(xué)出版社,2008:118～120[20] 霍宏濤. 數(shù)字圖像處理. 北京:北京理工大學(xué)出版社,2002:194～195[21] S. Arivazhagan, etc. A modified statistical approach for image fusion using wavelet ,2009（1）,137～144[22] 陳煜,田裕鵬. 圖像融合中的圖像配準(zhǔn)方法. 激光與紅外,2004,34（1）, 17～26[23] 謝曉鋒,楊之廉等. 微粒群算法綜述. 控制與決策,2003,18（2）,129～ 130 [24] 王進,季薇,鄭寶玉. 粒子種群優(yōu)化（PSO）算法的性能研究. 南京郵電學(xué)院學(xué)報,2005,25（4）,25～31[25] Eberhart R，Kennedy J，A new optimizer using particle swarm theory, Proc 6th Int Symposium on Micro Machine and Human Science, N agoya, 1995,39～43[26] JiannDer lee, ChungHsien Huang, etc. Improved accuracy of brain MRI/SPECT registration using a twocluster SPECT normalization algorithm and a binative similarity measure: application to the evaluation of Parkinson’s disease, Ann Nucl Med,2007,20～29[27] 張煜,劉哲星,李樹祥等. 醫(yī)學(xué)圖像信息融合技術(shù)的發(fā)展. 國外醫(yī)學(xué)生物醫(yī)學(xué)工程分冊,2000,23（4）,202～205[28] Pinz A, Bartl R. Information fusion in image understanding [J]. IEEE Trans Med Image,1992,12（4）,234～255[30]劉貴喜,楊萬海. 基于小波分解的圖像融合方法及性能評價. 西安電子科技大學(xué)自動化學(xué)報,2002,28（6）,27～45[31]張汗靈,楊帆. 基于互信息和混合優(yōu)化算法的多模醫(yī)學(xué)圖像配準(zhǔn). 湖南大學(xué)學(xué)報（自然科學(xué)版）,33（1）,2006,117～120[32]朱齊亮,樊要玲等. 基于互信息的醫(yī)學(xué)圖像配準(zhǔn)中的優(yōu)化算法. 科技信息, 2007（32）,368[33]Lau HY, Kagawa K,Lee WR, etc. CTMRI image fusion for 3D conformal prostate radiotherapy use in patients with altered pelvic anatomy [J]. The British Journal of Radiology,1996,69,1165～1170[34]張孫瑋,呂伯昇,張迅. 科技論文寫作入門. 北京:化學(xué)工業(yè)出版社,2007:171～179致 謝 本文的完成得到眾多老師和親人朋友的支持與幫助，借此機會： 深深感謝我的導(dǎo)師侯海燕老師對我設(shè)計上的悉心指導(dǎo)，從論文設(shè)計工作的進行，到畢業(yè)論文的撰寫和修改，傾注了侯老師的大量心血。她對工作嚴謹認真，一絲不茍的態(tài)度，更是我學(xué)習(xí)的榜樣。 感謝同組的同學(xué)對我的幫助，他們使我得到許多啟示。這是一個團結(jié)向上的集體，在大學(xué)度過這四年給我留下終生難忘的回憶。很榮幸結(jié)識各位師長和同學(xué)，正是這種愉悅的氛圍，讓我的大學(xué)生涯充滿快樂。 最后，我要特別感謝我的父母，我所有的成績都是他們的支持和鼓勵下才能取得的，他們的支持和鼓勵，使我一直不畏前路的任何困難，勇往直前，奮斗不息。附 錄%讀入圖像% Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global im[filename,pathname] =...uigetfile({39。*.jpg39。39。*.bmp39。39。*.gif39。 39。*.tif39。},39。選擇圖像39。)。str=[pathname filename]。im=imread(str)。axes()。imshow(im)。%圖像分割% Executes on button press in pushbutton3.function pushbutton3_Callback(hObject, eventdata, handles)% hObject handle to pushbutton3 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)global im[counts,x] = imhist(im)。J1=imhist(im)。J2=medfilt2(im)。J3=imadjust(J2,[ ],[0 1])。figure,imshow(J3)。J5=imhist(J3)。% figure,plot(J5)。n=4。wname=39。sym839。[xd,cxd,lxd]=wden(J5,39。rigrsure39。,39。s39。,39。mln39。,n,wname)。% figure,plot(xd)。[C,L]=wavedec(xd,5,39。db439。)。 %得到xd信號的小波系數(shù)cA5=appcoef(C,L,39。db439。,5)。 %得到xd信號的近似系數(shù)cD5=detcoef(C,L,5)。cD4=detcoef(C,L,4)。cD3=detcoef(C,L,3)。cD2=detcoef(C,L,2)。cD1=detcoef(C,L,1)。A5=wrcoef(39。a39。,C,L,39。db439。,5)。%重構(gòu)第五層小波近似分量% figure,plot(A5)。A6=diff(A5)。% figure。% plot(A6)k=1。 for i=5:250 if ((A6(i4)0)amp。amp。(A6(i3)0)amp。amp。(A6(i2)0)amp。amp。(A6(i1)0)amp。amp。(A6(i)0)amp。amp。(A6(i+1)0)amp。amp。(A6(i+2)0)amp。amp。(A6(i+3)0)amp。amp。(A6(i+4)0)) k=k+1。 th(k)=i。 endend[r,s]=size(J3)M=length(th)。th(1)=0。th(k+1)=256。J6=J3。for i=1:r for j=1:s for k=1:M1 if ((J3(i,j)=th(k))amp。amp。(J3(i,j)th(k+1))) J6(i,j)=ceil((th(k+1)th(k))/2+th(k))。 end end endendaxes()。imshow(J6)。%退出% Executes on button press in pushbutton4.function pushbutton4_Callback(hObject, eventdata, handles)% hObject handle to pushbutton4 (see GCBO)% eventdata reserved to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)close(