【正文】 2)。 %建立2列0矩陣為存放中心點坐標for i=1:m for j=1:n p1=a0(i)。 %找到骨架灰色區(qū)域m=length(a0)。skel39。 %轉化為黑色為1白色為0 為了后面find函數(shù)尋找1 find函數(shù)是尋找矩陣中為1的坐標并記錄 endendlk=edge(j0,39。 J0=imread (t)。十、附錄附錄一：半徑求解：clcjieguo=zeros (100,4)。如胰腺及其周圍血管的三維重建和顯示及臨床意義，此研究實現(xiàn)了胰腺及其周圍血管的三維重建和顯示,探討出一種內臟器官和血管的重建方式對斷層影像解剖學教學、臨床影像學診斷和治療方案的設計具有明確的輔助作用。如果知道象素圓的生成算法，我們就可以運用本模型的思想，將模型推廣到實數(shù)領域，從而使模型的精度提高。對于不規(guī)則的三維空間物體，我們的算法不適應，但可在其基礎上改進。七、模型評價與改進算法的精確度較高，具有很好的穩(wěn)定性和可信度，對于半徑均勻的血管有很好的適應性，在模型算法的過程中，我們在保證精確度的同時，使算法盡量簡單，從而提高編程的運行速度；運用函數(shù)確定輪廓和骨架的位置，直觀且容易理解；在確定半徑時，我們采用平均值法，以減少計算過程中產(chǎn)生的誤差；對于中軸線參數(shù)方程的確定，7階多項式擬合方式，使得曲線較為光滑，視覺效果較好。根據(jù)中軸線的參數(shù)方程及100張切片的最大內切圓的半徑，運用編程（見附錄三）我們得到切片疊加的原圖形和血管重建后的三維立體圖形見下圖圖8： 圖7 疊加后血管的還原圖 圖8重建后的血管三維圖像六、模型檢驗 運用求得的滾動球的半徑和中軸線，用球心沿中軸線運動的方法產(chǎn)生一簇球面，其包絡面生成一段新的血管。根據(jù)圖4分析，圖形是呈盤旋形狀的，隨著變大，在軸變化趨勢為：陡—平—陡。對于階多項式擬合中的最高階次的取值遵循兩個原則：1） 偏差平方和盡量小，這樣可以提高數(shù)據(jù)的擬合度；2）擬合最高次數(shù)不能太高，最好控制在9以內，當次數(shù)過高時，會增大誤差，自變量若有稍微的變動時，應變量的變化會非常大。骨架看作點的集合，計算骨架上每一點到輪廓的距離：取最小值，即該點垂直于外邊緣，也就是存在著以該點為圓心的圓相切于外邊緣。找到所有內點的半徑后, 其中半徑最大的內點即為所要找的球心, 對應的內切圓即為最大內切圓。具體位置見下圖2： 圖1 輪廓 圖2 骨架接下來以其中一張圖片為例進行詳細研究，運用函數(shù)確定輪廓和骨架的位置，并用二維矩陣來存儲輪廓線及骨架，這樣可以減少冗雜的信息量而又突出形狀特點，具體方法如下：在確定骨架時，逐次去掉切片圖像輪廓線上的點，同時要保證區(qū)域的連通性，余下的部分就是物體的骨架。所以，所以將原來單色圖像格式運用讀取函數(shù)轉化為二維的矩陣，其中1代表黑色像素點，0代表白色像素點。具體過程：由于血管是由半徑固定的球滾動包絡形成，所以過球心的截圓的半徑也就是血管的半徑，題中給出了血管被切片機平行切割后的100張切片圖像，這100張切片圖像是由無窮多個球被截的圓疊加而成的，又由于管道中軸線與每張切片有且只有一個交點，即球心。根據(jù)擬合出的方程，畫出中軸線在平面的投影，擬合曲線及立體圖形。所以每張切片的最大內切圓的圓心位于血管的中軸線上，該圓的半徑等于血管半徑，即問題就可以轉化為求切片上最大內切圓的半徑。對于問題一，要求解管道的半徑。問題三：根據(jù)中軸線及半徑對血管進行三維重建。為簡化起見，假設：管道中軸線與每張切片有且只有一個交點；球半徑固定；切片間距以及圖象象素的尺寸均為1。假設某些血管可視為一類特殊的管道，該管道的表面是由球心沿著某一曲線（稱為中軸線）的球滾動包絡而成。 斷面可用于了解生物組織、器官等的形態(tài)。作者簽名： 日期： 年 月 日導師簽名： 日期： 年 月 日指導教師評閱書指導教師評價：一、撰寫（設計）過程學生在論文（設計）過程中的治學態(tài)度、工作精神□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格學生掌握專業(yè)知識、技能的扎實程度□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格學生綜合運用所學知識和專業(yè)技能分析和解決問題的能力□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格研究方法的科學性；技術線路的可行性；設計方案的合理性□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格完成畢業(yè)論文（設計）期間的出勤情況□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格二、論文（設計）質量論文（設計）的整體結構是否符合撰寫規(guī)范？□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格是否完成指定的論文（設計）任務（包括裝訂及附件）？□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格三、論文（設計）水平論文（設計）的理論意義或對解決實際問題的指導意義□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格論文的觀念是否有新意？設計是否有創(chuàng)意？□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格論文（設計說明書）所體現(xiàn)的整體水平□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格建議成績：□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格（在所選等級前的□內畫“√”）指導教師： （簽名） 單位： （蓋章）年 月 日評閱教師評閱書評閱教師評價：一、論文（設計）質量論文（設計）的整體結構是否符合撰寫規(guī)范？□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格是否完成指定的論文（設計）任務（包括裝訂及附件）？□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格二、論文（設計）水平論文（設計）的理論意義或對解決實際問題的指導意義□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格論文的觀念是否有新意？設計是否有創(chuàng)意？□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格論文（設計說明書）所體現(xiàn)的整體水平□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格建議成績：□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格（在所選等級前的□內畫“√”）評閱教師： （簽名） 單位： （蓋章）年 月 日教研室（或答辯小組）及教學系意見教研室（或答辯小組）評價：一、答辯過程畢業(yè)論文（設計）的基本要點和見解的敘述情況□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格對答辯問題的反應、理解、表達情況□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格學生答辯過程中的精神狀態(tài)□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格二、論文（設計）質量論文（設計）的整體結構是否符合撰寫規(guī)范？□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格是否完成指定的論文（設計）任務（包括裝訂及附件）？□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格三、論文（設計）水平論文（設計）的理論意義或對解決實際問題的指導意義□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格論文的觀念是否有新意？設計是否有創(chuàng)意？□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格論文（設計說明書）所體現(xiàn)的整體水平□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格評定成績：□ 優(yōu) □ 良 □ 中 □ 及格 □ 不及格（在所選等級前的□內畫“√”）教研室主任（或答辯小組組長）： （簽名）年 月 日教學系意見：系主任： （簽名）年 月 日一、問題重述生物體的外部結構具有繁雜多樣性，可以通過肉眼觀察，但若想了解去內部錯綜復雜的結構，就需要借助一定的輔助工具，人們常利用的是分解和