【正文】 在論文的撰寫過程中老師們給予我很大的幫助，幫助解決了不少的難點，使得論文能夠及時完成，這里一并表示真誠的感謝。最后，我要特別感謝我的導(dǎo)師趙達(dá)睿老師、和研究生助教熊偉麗老師。感謝老師四年來對我孜孜不倦的教誨，對我成長的關(guān)心和愛護(hù)。最后，我要感謝我的父母對我的關(guān)系和理解，如果沒有他們在我的學(xué)習(xí)生涯中的無私奉獻(xiàn)和默默支持，我將無法順利完成今天的學(xué)業(yè)。從他身上，我學(xué)到了許多能受益終生的東西。這期間凝聚了很多人的心血，在此我表示由衷的感謝。首先非常感謝學(xué)校開設(shè)這個課題，為本人日后從事計算機(jī)方面的工作提供了經(jīng)驗，奠定了基礎(chǔ)。畢業(yè)設(shè)計（論文）使用授權(quán)聲明本人完全了解濱州學(xué)院關(guān)于收集、保存、使用畢業(yè)設(shè)計（論文）的規(guī)定。對本文的研究做出重要貢獻(xiàn)的個人和集體均已在文中以明確方式標(biāo)明。本人完全意識到本聲明的法律結(jié)果由本人承擔(dān)。 end endend%求百分比baifenbi=chonghegs/nbiao。for i=1:512 for j=1:512 if(biaozhunjuzheng(i,j)==0) nbiao=nbiao+1。 %通過矩陣的或運算得到填充后的新圖imshow(you|zuo)%原圖的黑點的個數(shù)biaozhunjuzheng=imread(strcat(int2str(pn),39。%再填充右邊界的左半部分for i=1:512 for j=512:1:2 if(juzheng(i,j)==0amp。計算新舊切片的重合度function baifenbi=baifenbi(pnjj,pn) %輸出擬合圖與原切片圖的重合度%填充新圖juzheng=pnjj。endfigure(1)。u=round(u)。lk=edge(pnj,39。)axis([0 512 0 512])。end%根據(jù)新的切平面的輪廓坐標(biāo)得到新輪廓的圖像矩陣plot(v(:,1),u(:,1),39。[real(r(1))+256 real(r(2))+256]]。 [g,h]=solve(m,n)。v=[]。 %中心軸線方程關(guān)于z的導(dǎo)數(shù)即[a(k),b(k),1]為z在k處的切線的方向向量q(k)=px(1)*c.^7+px(2)*c.^6+px(3)*c.^5+px(4)*c.^4+px(5)*c.^3+px(6)*c.^2+px(7)*c+px(8)。px=[,]。w=zeros(991)。 else b0(i,j)=1。title(39。)。X軸39。red39。%?。?！存儲角度for i=1:100 x=*sin(theta).*sin(phi)+cenn(i,1)。cenn(:,3)=z。 y1=polyval(py,z)。%通過matlab擬合工具箱發(fā)現(xiàn)zx平面7次擬合已經(jīng)足夠，擬合效果足夠。%讀入中心坐標(biāo)x=b(:,1)。)。Z軸39。ylabel(39。,2)。% 繪制三維螺旋線figureplot3(x, y, t,39。,jieguo)%導(dǎo)入excle方便使用數(shù)據(jù)附錄二：最大內(nèi)切圓半徑及圓心坐標(biāo)最大內(nèi)切圓圓心坐標(biāo)最大內(nèi)切圓圓心坐標(biāo)切片序號最大內(nèi)切圓半徑橫坐標(biāo)縱坐標(biāo)豎坐標(biāo)切片序號最大內(nèi)切圓半徑橫坐標(biāo)縱坐標(biāo)豎坐標(biāo)401506341704616371411496642714616372421486843724616373431486844736515874441437845746815775451378846756515876461378847768115277471161154877811527848115116497881152794911511650791351188050114117518013611781511141175281136117825211311853821371168353112119548313811584541111205584138115855511112056851391148656631515786139114875775145588713911488588114259881401138959511566089140113906051156619017267916131162629117267926231162629217267936331162629317267946431162629418243956531162629518724966631162629618724976731162629718724986831162629818724996931162629918818100附錄三：擬合后三維直線圖：% 繪制螺旋線% 繪制原理：隨著時間的延續(xù)或z坐標(biāo)的升高，t=0:99。%Y軸坐標(biāo)jieguo(k+1,3)=k+1。 %找到其中最大的半徑，并把半徑和圓心坐標(biāo)存儲%%%%%%%%%%%%%%%%%%%%%%輸出結(jié)果x=a0(zdxh)256。 % 通過循環(huán)求各個骨架上的各個點到輪廓的距離 end [zx,zxxh]=min(jl(i,:))。 q1=b0(i)。 %輪廓灰色區(qū)域個數(shù)n=length(x0)。,inf)。sobel39。%將bmp格式轉(zhuǎn)化為01矩陣，黑色為0，白色為1for i=1:1:512。%存儲以后生成的結(jié)果for k=0:99 t=strcat (int2str (k),39?；跀?shù)字虛擬人體血管的三維重建對教學(xué)和臨床的應(yīng)用都有一定實際意義。此模型具有一定的推廣意義，尤其是在醫(yī)學(xué)觀測領(lǐng)域，醫(yī)生在不能破壞人體構(gòu)造前提下,只能通過X 光技術(shù)獲得斷層圖像信息，然而僅憑一維的斷層圖像,醫(yī)生是很難了解其內(nèi)部的復(fù)雜結(jié)構(gòu)的，因此，序列圖像的三維重現(xiàn)技術(shù)卻使斷層圖像得到立體還原，使人體內(nèi)部構(gòu)造一目了然。八、模型推廣如果考慮中軸線與切片的交點不止一個，比如有兩個交點，仍然可以利用我們的離散模型去求解，因為不含軸心的切片中一定不含有最大圓。模型不能夠?qū)Π霃阶兓那闆r適應(yīng)，這是因為我們的理論根據(jù)就是等徑管道的幾何特性，另外，由于實際中圖像邊界上的點是連續(xù)的，只有位置，沒有大小，圖像在轉(zhuǎn)換為格式時，像素表示的邊界是離散的，產(chǎn)生了系統(tǒng)誤差。用原來100張的平面截新的血管，生成新的100張橫斷面圖像，抽樣比較新、舊100張橫截面圖像之間的差別。結(jié)合圖3和圖5分析隨著的不斷變大，圖形盤旋上升的趨勢為：平—陡—平。根據(jù)這兩個原則確定方程的最高階次，此時擬合出的中軸線較為光滑。最后取最小距離的最大值，即是最大內(nèi)切圓的半徑：（4）100張圖片的半徑完成骨架上任一點對輪廓線像素點的遍歷求距離, 并比較記錄的所有最小距離值, 最小距離的最大值即為該切片圖像的最大內(nèi)切圓的半徑，通過循環(huán)100次可以得到這100張圖片的最大內(nèi)切圓的半徑。若一張切片出現(xiàn)多個圓心時，鑒于所給切片是格式的圖像, 是對原血管的切片的連續(xù)圖形離散化而得到的近似圖形, 則實際計算時有可能出現(xiàn)在一個切片中同時找到多個最大內(nèi)切圓, 即內(nèi)切圓不惟一。運用邊緣檢測大幅度地減少了數(shù)據(jù)量，并且剔除了可以認(rèn)為不相關(guān)的信息，保留了圖像重要的結(jié)構(gòu)屬性。對于一個像素的圖像，每一個像素都有自己的一個確定的坐標(biāo)。所以在每張切片上總存在著以交點為圓心的圓，此時半徑最大，這最大半徑同時也是管道的半徑。對于問題三，基于問題一、二的求解結(jié)果，根據(jù)中軸線的