【導讀】DICOM標準是醫(yī)學影像信息學領域的國際通用標準。在臨床應用中常需要獲。門的應用軟件來解決這類具體問題。本文研究了應用matlab對醫(yī)學骨科圖像進行處理以獲。取相關信息的方法。首先研究了對圖像進行邊緣提取的常用方法經(jīng)典算子法和基于數(shù)學形。接著研究了區(qū)域質(zhì)心，圓弧的圓心和半徑以及兩點之間的距離等量化信息的測量方法。最后應用matlab平臺設計。獨立進行研究工作所取得的成果，成果不存在知識產(chǎn)權爭議?；蜃珜戇^的作品成果。對本文的研究做出重要貢獻的個人和集體均已在文中以。本聲明的法律后果由本人承擔。本人完全了解濱州學院關于收集、保存、使用畢業(yè)設計（論文）的規(guī)定。文）；同意學校在不以營利為目的的前提下，建立目錄檢索與閱覽服務系統(tǒng)，公布設計（論文）的部分或全部內(nèi)容，允許他人依法合理使用。

　　

【正文】 所示，則： 22( ) ( )i i id X A Y B? ? ? ? (4 6)? 標記質(zhì)心 尋找連通區(qū)域 標記連通區(qū)域 廈門大學通信工程系 20xx 級本科畢業(yè)論文 31 點 到圓邊緣的距離的平方與半徑的平方的差為： 2 2 2 2i i i i i id R X Y a X b Y c? ? ? ? ? ? ? ? (4 7)? ( , )iiXY id ( , )AB 圖 47 令 為的平方和： 22 2 2( , , ) ( )i i i i iQ a b c X Y a X b Y c? ??? ? ? ? ? ????? (4 8)? 即求參數(shù) ,使得 的值最小。對 求偏導 得并令其為 0 得： 222222( , , ) 2 ( ) 0( , , ) 2 ( ) 0( , , ) 2 ( ) 0i i i i ii i i i ii i i iQ a b c X Y a X b Y c XaQ a b c X Y a X b Y c YbQ a b c X Y a X b Y cc?? ? ? ? ? ? ?? ???? ? ? ? ? ? ?? ???? ? ? ? ? ? ?? ????? (49)? 解此方程組可得 的極值，由于平方差大于 0，因此存在極小值，無極大值，比較所有的極值點即可得到極小值。將 取極小值時的 代入式 45 即得所求。 在編寫程序時 將最小二乘法編成一個函數(shù) circfit()， 在編寫回調(diào)函數(shù)時調(diào)用 。 (3) 任意兩點間距離的測量 要測量圖像中任意兩點間的距離相對簡單，首先獲取起點和終點的坐標值，然后用兩點間的距離公式即可求出結(jié)果。 R 廈門大學通信工程系 20xx 級本科畢業(yè)論文 32 第 5 章 全文總結(jié)以及展望 本文研究內(nèi)容總結(jié) 本文主要研究了醫(yī)學圖像 DICOM 以及應用 MATLAB 對醫(yī)學骨科圖像進行處理 ，提取目標圖像 的邊緣 ， 并在此基礎上 實現(xiàn)對區(qū)域 重心 ， 圓心 和半徑及 兩點間距離等 信息的 測量 。 最后應用 Matlab 平臺設計出一個 骨科圖像信息處理系統(tǒng)。 本文的特色和優(yōu)點 本文研究的內(nèi)容有以下幾個特色和優(yōu)點： ◆針對骨科圖像的 具體 特點 ， 結(jié)合圖像預處理 ， 經(jīng)典邊緣算子 ，數(shù)學 形態(tài)學 等方法較好的 提取了目標圖像的 邊緣； ◆充分利用 MATLAB 強大的數(shù)學計算功能 ， 在提取出邊緣的基礎上實現(xiàn)了對 圖像相關信息的 定量測量 ； ◆利用 GUIDE 設計出一個處理以上信息的 信息處理系統(tǒng)， 使研究的成果能夠更好的展現(xiàn)以及應用于實踐 。 在研究過程中 遇到的問題 在研究本課題中遇到 過 很多問題 ，主要的 有以下幾個 方面 : ◆所用的算法難以適應眾多的圖像 使用單一的算法進行邊緣提取時往往會出現(xiàn)這樣的情況 ： 一種算子對提取這張圖像效果較好 ， 但應用于另處一張圖像時效果就沒有其它的算法好 。 為了使所用的算法具有一定的普遍適應性 ， 將多種方法結(jié)合起 ， 最后通過實踐發(fā)現(xiàn)將圖像歸一化 ， 經(jīng)典算子和形態(tài)學 廈門大學通信工程系 20xx 級本科畢業(yè)論文 33 方法結(jié)合起來 ， 效果不錯 。 ◆ 對于存在斷點的邊緣 難以提取重心 因為質(zhì) 心的提取主要是依靠尋找連通區(qū)域 ，因此對于有斷點的邊緣以提取質(zhì)主， 最后通過形態(tài)學處理將斷點連接起來 ， 使問題得以解決 。 ◆ 將單個算法程序 綜合到 GUIDE 中的轉(zhuǎn)換 由于剛開始是一個問題一個問題的解決 ， 到 最 后應用 GUIDE 進行綜合編程時有許多地方不同 ， 因此走了很多彎路 ， 后來通過老師的指導和更深入的研究 ， 使問題得以解決 。 研究課題的展望 在研究 的過程中 我 發(fā)現(xiàn)還有以下幾點可以做進一步的研究 : ◆實現(xiàn)應用 MATLAB 對 DICOM 醫(yī)學圖像的直接處理 。 由于 MATLAB 還不能對 DICOM 進行直接的處理 ， 而是要先將其轉(zhuǎn)化成 MATLAB 可以處理的圖像格式 ， 因此可以嘗試應用 MATLAB將其直接轉(zhuǎn)換成常規(guī)圖片格式 。 ◆ 研究更 具普遍適應性和穩(wěn)健性的邊緣提取算法 。 ◆設計出功能更強大 ， 測量 信息更多更精確， 更具實踐性和操作性的系統(tǒng) 。 廈門大學通信工程系 20xx 級本科畢業(yè)論文 34 致 謝語 首先感謝我的指導老師趙彩丹老師，她為我提供了良好的學習環(huán)境并且在完成畢業(yè)設計過程中，從 確定題目、解決一個個問題， 論文修改直到完成，這每一步都離不開她所給予 的 悉心指點和幫助。在此對趙老師致以誠摯的感謝！同時我還要感謝 108 實驗室的學長學姐以及我身邊的朋友們，感謝他們在完成畢業(yè)設計過程中曾給予的幫助和支持。通過畢業(yè)設計，我得到了極大鍛煉，分析 、 研究 和 解決問題的能力得到提高，這一切都離不開大家的幫助，謝謝！ 廈門大學通信工程系 20xx 級本科畢業(yè)論文 35 參 考文獻 [1] 高向軍，田聯(lián)房，王立非等 .利用 MATLAB 實現(xiàn)醫(yī)學圖像處理與分析 [J].計算機應用 ,20xx,25(6):1. [2] 秦襄培 .MATLAB 圖像處理與界面編程寶典 [M].北京 :電子工業(yè)出版社 ,20xx. [3] 孫亦博 ,朱玲利 .DIOCM 醫(yī)學圖像與 BMP 格式的轉(zhuǎn)換研究 [J].現(xiàn)代計算機， . [4] 姜炳旭，劉杰，孫可 .Sobel 邊緣檢測的細化 [J].沈陽師范大學學報（自然科學版）， 20xx， 28（ 4）： 504. [5] 王娟，周金芝 .基于 Matlab 的形態(tài)學圖像處理研究 [J].計算機論壇， 20xx. [6] 李晨丹，徐進 .指紋圖像預處理和特征提取的 Matlab 實現(xiàn) [J].計算機工程與科學， 20xx， 31（ 7）： 62. [7] 徐國旺，廖明潮 .擬合圓的幾種方法 [J].武漢工業(yè)學院學報， 20xx. 廈門大學通信工程系 20xx 級本科畢業(yè)論文 36 附錄 (程序代碼 )： 下面是骨科圖像信息處理系統(tǒng)實現(xiàn)的 matlab 代碼及相關注釋： function varargout = bishe(varargin) gui_Singleton = 1。 gui_State = struct(39。gui_Name39。, mfilename, ... 39。gui_Singleton39。, gui_Singleton, ... 39。gui_OpeningF39。, @bishe_OpeningF, ... 39。gui_OutputF39。, @bishe_OutputF, ... 39。gui_LayoutF39。, [] , ... 39。gui_Callback39。, [])。 if nargin amp。amp。 ischar(varargin{1}) = str2func(varargin{1})。 end if nargout [varargout{1:nargout}] = gui_mainf(gui_State, varargin{:})。 else gui_mainf(gui_State, varargin{:})。 廈門大學通信工程系 20xx 級本科畢業(yè)論文 37 end function bishe_OpeningF(hObject, eventdata, handles, varargin) = hObject。 guidata(hObject, handles)。 function varargout = bishe_OutputF(hObject, eventdata, handles) varargout{1} = 。 %讀取圖片 function pushbutton1_Callback(hObject, eventdata, handles) [fname,fpath]=uigetfile({39。*.dcm。*.jpg39。},39。File Selector39。)。 =fname。 [pathstr, name, ext, versn] = fileparts(fname)。 if ext==39。.dcm39。 =diread([fpath fname])。 elseif ext==39。.jpg39。 =imread([fpath fname])。 else warndlg(39。請輸入 DICOM 或 JPEG 圖片文件 39。,39。Warning!39。)。 end axes()。 imshow(,[])。 廈門大學通信工程系 20xx 級本科畢業(yè)論文 38 guidata(hObject,handles)。 %選取區(qū)域 function pushbutton2_Callback(hObject, eventdata, handles) rect = round(getrect())。 =(rect(2):rect(2)+rect(4),rect(1):rect(1)+rect(3))。 axes()。 imshow(,[])。 guidata(hObject,handles)。 %提取邊緣 function pushbutton3_Callback(hObject, eventdata, handles) I=。 I1=double(I)。 M=mean2(I)。 VAR=std2(I)。 M0=。 V0=。 [row,col]=size(I1)。 for i=1:row。 廈門大學通信工程系 20xx 級本科畢業(yè)論文 39 for j=1:col。 if(I1(i,j)M) Nor(i,j)=M0+(V0*(I1(i,j)M)^2/VAR)^(1/2)。 else Nor(i,j)=M0(V0*(I1(i,j)M)^2/VAR)^(1/2)。 end end end T=50。 X0=I1。 [row,col]=size(X0)。 g=zeros(row,col)。 X=Nor。 for i=2:row1 for j=2:col1 Dx=(X(i+1,j1)X(i1,j1)+2*(X(i+1,j)X(i1,j))+X(i+1,j+1)X(i1,j+1))^2。 Dy=(X(i1,j+1)X(i1,j1)+2*(X(i,j+1)X(i,j1))+X(i+1,j+1)X(i+1,j1))^2。 g(i,j)=round(sqrt(Dx+Dy))。 end end 廈門大學通信工程系 20xx 級本科畢業(yè)論文 40 for i=1:row for j=1:col if g(i,j)T X(i,j)=255。