專利名稱：人體骨骼未知體三維交互模型的建立方法及其應(yīng)用的制作方法
技術(shù)領(lǐng)域：
本發(fā)明屬于計算機成像技術(shù)領(lǐng)域，特別是醫(yī)學影像及手術(shù)模擬技術(shù)領(lǐng)域。其中包括一種人體骨骼未知體三維交互模型的建立方法和一種人體骨性三維交互模擬手術(shù)的實現(xiàn)方法。
背景技術(shù)：
在科學研究領(lǐng)域，研究的主要目的是理解自然的本質(zhì)。要達到這個目的，需經(jīng)過從觀察自然現(xiàn)象到模擬自然現(xiàn)象并分析模擬結(jié)果的過程，這依賴于科技進步發(fā)展。
具體到醫(yī)學方面，為了能夠準確認識疾病，為手術(shù)方案的制訂提供依據(jù)，計算機手術(shù)模擬應(yīng)運而生，是現(xiàn)代外科的重要發(fā)展方向，該領(lǐng)域的研究是計算機科學、數(shù)學、機械學、外科及生物醫(yī)學工程學等多學科高技術(shù)的結(jié)晶，已成為各國科技界研究的一個熱點。
日本學者Fujino于1989年第一次提出模擬手術(shù)的概念，即在計算機產(chǎn)生的三維實體模型上立體展示病變情況及手術(shù)過程。Murphy研究中提出了由CT數(shù)據(jù)來計劃全髖關(guān)節(jié)置換手術(shù)的概念。1990年Delp等人建立了人體下肢的圖形交互模型，進行下肢手術(shù)模擬的研究。近年來，隨著計算機的飛速發(fā)展，這個領(lǐng)域的研究不斷深入。1996年Koch利用有限元模型系統(tǒng)通過改變骨骼和軟組織的位置來模擬面部整形手術(shù)。1999年Christopher利用有限元模型進行髖關(guān)節(jié)穩(wěn)定性預(yù)測分析；Neumann，P.報道了利用三維切割交互模型進行頜面外科的手術(shù)設(shè)計的研究。國內(nèi)已有一些開發(fā)計算機手術(shù)模擬系統(tǒng)的報道，如計算機輔助上頜面手術(shù)模擬，計算機輔助放射治療。上述模擬手術(shù)的人體模型及目前國際上美國、韓國及我國建立的虛擬人，均屬于正常人已知體(即直接從人體或尸體上獲得數(shù)據(jù))的三維重建模型，與實際病人骨骼疾病本身具有著一定誤差。實際上，人體骨骼系統(tǒng)疾病的病理模型的具體數(shù)據(jù)并不能從人體或尸體上直接獲得，這些結(jié)構(gòu)屬于未知體。以兒童先天性髖脫位為例，不能直接獲得數(shù)據(jù)，需通過螺旋CT(SVCT)掃描獲得數(shù)據(jù)。
現(xiàn)有螺旋CT設(shè)備雖然能進行三維模型重建，但卻不能對模型進行任意切割，只能用于觀察，目前尚無法實現(xiàn)手術(shù)模擬。
計算機模擬手術(shù)的關(guān)鍵在于人體(也包括動物體等，為敘述方便，以下統(tǒng)稱人體)三維交互模擬的建立。通常三維模型的重建是利用三維數(shù)據(jù)場等值面生成的經(jīng)典算法Marching Cubes法，即移動立方體法。該算法由Lorensen等人于1987年提出，它對體數(shù)據(jù)中的體素(又稱體元)進行處理，生成三角面片集。體素是體數(shù)據(jù)的基本組成單位，每八個相鄰的采樣點所定義的立方體區(qū)域構(gòu)成一個體素。一個體素中可能包含多個三角面片?，F(xiàn)有手術(shù)模擬系統(tǒng)在進行三維重建時就是采用這種算法。為將重構(gòu)的三維模型用于手術(shù)模擬，最好要求重構(gòu)模型中的三角面片不但具有連接性，而且能夠快速定位，這樣才能滿足實時交互進行三維切割、移位、旋轉(zhuǎn)和拼合的需要。然而基于現(xiàn)有的這種Marching Cubes方法生成的三維重構(gòu)模型都是散亂的三角面片集，所得到的三角面片并非都是相連的，這種模型結(jié)構(gòu)在對模型進行切割、取舍等全局性模擬操作時帶來不便，使得現(xiàn)有此類系統(tǒng)特別復(fù)雜，通常需要在工作站上才能實現(xiàn)。
而且，現(xiàn)有的手術(shù)模擬系統(tǒng)一般是利用原始圖像數(shù)據(jù)直接進行三維重構(gòu)，這樣得到的三角面片集數(shù)量巨大，相應(yīng)增加了運算的復(fù)雜性，不但對設(shè)備要求更高，而且其反應(yīng)速度很慢。
而且，現(xiàn)有手術(shù)模擬系統(tǒng)并不能實現(xiàn)三維交互模型在任意位置的切割、平移、旋轉(zhuǎn)、拾取和拼合。

發(fā)明內(nèi)容
基于現(xiàn)有技術(shù)的不足，本發(fā)明所要解決的技術(shù)問題是，提供一種人體骨骼未知體三維交互模型的建立方法，能使建立過程中生成的三角面片相連，并可以實現(xiàn)快速查找和定位。
本發(fā)明所要解決的另一技術(shù)問題是，提供一種人體骨骼未知體三維交互模型的建立方法，其還可以簡化模型數(shù)據(jù)，提高系統(tǒng)反應(yīng)速度。
本發(fā)明所要解決的再一技術(shù)問題是，提供一種人體骨性三維交互模擬手術(shù)的實現(xiàn)方法，實現(xiàn)骨骼三維交互模型在任意位置的切割、平移、旋轉(zhuǎn)、拾取和拼合。
為解決上述技術(shù)問題，本發(fā)明提供了如下技術(shù)方案。
人體骨骼未知體三維交互模型的建立方法，包括一圖像數(shù)據(jù)獲取步驟、一模型數(shù)據(jù)提取步驟及一模型數(shù)據(jù)重建步驟，所述模型數(shù)據(jù)重建步驟包括一體素處理步驟，是利用移動立方體法(Marching Cubes法)對體數(shù)據(jù)的體素進行處理，生成三角面片集；一模型記錄步驟，是利用體數(shù)據(jù)的體素在模型中所處的行、層、分組位置與B+樹結(jié)構(gòu)的對應(yīng)關(guān)系，通過B+樹數(shù)據(jù)結(jié)構(gòu)方式，重新組織和記錄三維模型；以及一模型重建步驟，是利用B+樹結(jié)構(gòu)的數(shù)據(jù)搜索方式，根據(jù)B+樹結(jié)構(gòu)下記錄的三維模型的體素的連接關(guān)系，對體素及其中的三角面片定位，實現(xiàn)三維模型的重建。
其中，所述圖像數(shù)據(jù)獲取步驟是采用螺旋CT對包含未知體的實體進行掃描，獲得原始圖像數(shù)據(jù)。
其中，所述模型數(shù)據(jù)提取步驟是根據(jù)灰度等級對原始圖像數(shù)據(jù)進行選擇，對數(shù)據(jù)開窗進行處理，提取出其中的模型圖像數(shù)據(jù)，并去除圖像中的干擾數(shù)據(jù)。
其中，在所述模型數(shù)據(jù)提取步驟后，還包括一模型數(shù)據(jù)壓縮步驟，是根據(jù)原始圖像數(shù)據(jù)的密度，對模型圖像數(shù)據(jù)進行兩級重采樣，以簡化數(shù)據(jù)，實現(xiàn)數(shù)據(jù)壓縮。
其中，進行兩級重采樣所用的是領(lǐng)域平均濾波器。
所說的人體骨性三維交互模擬手術(shù)的實現(xiàn)方法包括如下步驟骨胳模型建立步驟，利用以上任一權(quán)利要求所述的未知體三維交互模型的建立方法，建立骨骼三維交互模型，其模型圖像數(shù)據(jù)以B+樹數(shù)據(jù)結(jié)構(gòu)形式組織記錄；骨刀模型建立步驟，根據(jù)骨刀曲面方程，確定骨刀模型所需物體空間大小，并進行體素處理步驟，是利用移動立方體法對體數(shù)據(jù)的體素進行處理，生成三角面片集；模型記錄步驟，是利用體數(shù)據(jù)的體素在模型中所處的行、層、分組位置與B+樹結(jié)構(gòu)的對應(yīng)關(guān)系，通過B+樹數(shù)據(jù)結(jié)構(gòu)方式，重新組織和記錄三維模型；模型重建步驟，是利用B+樹結(jié)構(gòu)的數(shù)據(jù)搜索方式，根據(jù)B+樹結(jié)構(gòu)下記錄的三維模型的體素的連接關(guān)系，對體素及其中的三角面片定位，實現(xiàn)三維模型的重建，其中，該骨刀模型的網(wǎng)格劃分方式及大小與所述未知體三維模型相同；以及手術(shù)模擬步驟，包含骨刀模型的平移和/或旋轉(zhuǎn)步驟、骨骼模型的切割步驟以及骨骼模型的拼合步驟或其組合。
其中，所述手術(shù)模擬步驟中的骨刀的平移和/或旋轉(zhuǎn)步驟，是根據(jù)平移和/或旋轉(zhuǎn)變換矩陣，生成骨刀曲面方程的新描述信息，再建立該新的骨刀模型，實現(xiàn)骨刀的平移或旋轉(zhuǎn)。
其中，所述手術(shù)模擬步驟中的未知體模型切割步驟，是利用骨刀模型的B+樹數(shù)據(jù)結(jié)構(gòu)中各層索引節(jié)點內(nèi)包含的層次包圍盒信息，通過對骨骼模型的B+樹數(shù)據(jù)結(jié)構(gòu)的不完全遍歷和回溯，定位骨骼模型與骨刀模型表面間的節(jié)點并分割骨骼模型。
其中的手術(shù)模擬步驟中的骨骼模型的拼合步驟，是根據(jù)變換矩陣變換未知體模型移動部分的所有面片的頂點的坐標和法向，再由新面片集生成新的B+樹，通過骨骼模型移動部分的平移和/或旋轉(zhuǎn)，實現(xiàn)骨骼模型的拼合。
所述骨刀模型可以包括平面骨刀模型和/或球面骨刀模型，其中，平面骨刀模型由以下參數(shù)確定局部坐標系點在全局坐標系中的坐標、刀平面方程和刀邊界形狀方程；球面骨刀模型由以下參數(shù)確定局部坐標系點在全局坐標系中的坐標、球面半徑和刀邊界形狀方程。
本發(fā)明具有如下有益效果本發(fā)明的人體骨骼未知體三維模型的建立方法由于采用了改進的移動立方體法，用B+樹數(shù)據(jù)結(jié)構(gòu)方式重新組織記錄三維模型，保證了三角面片的連接，可以實現(xiàn)快速查找和定位，建模效率高。并且，由于采用了螺旋CT獲取原始圖像掃描數(shù)據(jù)，有利于模型的真實再現(xiàn)，并且由于簡化了模型數(shù)據(jù)，可以進一步提高系統(tǒng)反應(yīng)速度。
本發(fā)明的人體骨性三維交互模擬手術(shù)的實現(xiàn)方法的骨骼和骨刀三維模型具有相同的網(wǎng)格劃分方式和大小，以及相同的B+樹結(jié)構(gòu)，實現(xiàn)了未知體交互模型在任意位置的切割、平移、旋轉(zhuǎn)、拾取和拼合。
通過該骨科三維交互模型，醫(yī)生可在虛擬環(huán)境中容易進行手術(shù)規(guī)劃和手術(shù)過程模擬，比較不同的手術(shù)方案和結(jié)果，以幫助制定最佳手術(shù)方案。


下面結(jié)合附圖，對本發(fā)明的具體實施方式
作詳細說明。
圖1為本發(fā)明所述人體骨骼未知體三維交互模型的重建方法的一種實施例的流程示意圖；圖2為先天性髖脫位患兒SVCT掃描數(shù)據(jù)經(jīng)二級重采樣再進行三維重構(gòu)的三角面片集簡化結(jié)果列表；圖3為本發(fā)明所述人體骨性三維交互模擬手術(shù)實現(xiàn)方法的一種實施例的流程示意圖。
具體實施例方式
參見圖1，所示為本發(fā)明所提出的人體骨骼未知體三維模型建立方法的一種實施方式的流程示意圖。該方法包括一圖像數(shù)據(jù)獲取步驟、一模型數(shù)據(jù)提取步驟、一模型數(shù)據(jù)壓縮步驟以及一模型數(shù)據(jù)重建步驟。
對于其中的圖像獲取步驟，由于未知體的三維重建不能通過在實體上直接測量得到體數(shù)據(jù)，所以原始圖像中體數(shù)據(jù)一般是利用CT掃描設(shè)備對包含未知體的實體進行掃描而獲得。
但是，傳統(tǒng)的CT由于接口問題，阻礙直接獲取原始掃描數(shù)據(jù)，需要依靠斷層膠片重建模型，其結(jié)果是數(shù)據(jù)丟失多，誤差大，失真明顯。因此，最好是利用螺旋CT。螺旋CT是連續(xù)不間斷地收集數(shù)據(jù)，克服了傳統(tǒng)CT的缺點，可直接提供高精度原始整體掃描數(shù)據(jù)，從而能更好地保證重建模型的精確性和真實性。
對于其中的模型數(shù)據(jù)提取步驟，是從通過圖像數(shù)據(jù)獲取步驟得到的原始圖像數(shù)據(jù)中選取適當灰度等級的數(shù)據(jù)開窗進行處理，提取出其中的模型圖像數(shù)據(jù)，并去除圖像中的干擾數(shù)據(jù)。
如果直接利用上述螺旋CT掃描數(shù)據(jù)進行三維重構(gòu)，那么得到的三角面片的量將非常大。例如，對于先天性髖脫位患兒的螺旋CT掃描數(shù)據(jù)直接進行重構(gòu)，得到的三角面片集大約含有290000個三角面片。如此龐大數(shù)量的三角面片將使三維重構(gòu)和響應(yīng)實時交互操作變得極為困難，對計算機的要求將非常高。為解決此問題，可對上述螺旋CT掃描數(shù)據(jù)進行簡化壓縮。
對高精度的螺旋CT掃描數(shù)據(jù)進行簡化可以采用三種方法。第一種方法是先求得精細模型，再對精細模型進行優(yōu)化和簡化，其計算量大。第二種方法是對重構(gòu)算法進行優(yōu)化，結(jié)果是重構(gòu)計算本身變得復(fù)雜，重構(gòu)時間長。第三種方法是先簡化原始數(shù)據(jù)，再生成簡化模型。其優(yōu)點是可以在不產(chǎn)生精細模型的情況下由簡化數(shù)據(jù)直接生成簡化模型，計算效率較高。在本發(fā)明的實施例中具體采用該第三種方法。
發(fā)明人在對上述先天性髖脫位患兒的螺旋CT掃描數(shù)據(jù)的研究中，分別采用三種濾波器進行數(shù)據(jù)重采樣實驗，以研究數(shù)據(jù)簡化效果。螺旋CT原始數(shù)據(jù)是規(guī)整網(wǎng)格數(shù)據(jù)，可直接選取重采樣點原始數(shù)據(jù)值作為簡化數(shù)據(jù)中該點的數(shù)據(jù)值。參見圖2，所示為簡化結(jié)果列表。從中可以看出，對掃描數(shù)據(jù)經(jīng)重采樣簡化再進行三維重構(gòu)可分別使產(chǎn)生的三角面片數(shù)據(jù)集壓縮到精細模型的20％(經(jīng)一級重采樣)和4％(經(jīng)二級重采樣)。研究結(jié)果表明，采用最簡便的領(lǐng)域平均濾波器進行兩級重采樣簡化數(shù)據(jù)生成的模型保留了原始數(shù)據(jù)中的信息，降低了計算量，取得了滿意的效果，可以滿足實時交互的需要。
因此，在本實施例中，所述模型數(shù)據(jù)壓縮步驟是根據(jù)原始圖像數(shù)據(jù)的密度，對模型圖像數(shù)據(jù)進行兩級重采樣，以簡化數(shù)據(jù)，實現(xiàn)數(shù)據(jù)壓縮，進行兩級重采樣所用的是領(lǐng)域平均濾波器。
對于其中的模型數(shù)據(jù)重建步驟，其具體包括一體素處理步驟、一模型記錄步驟及一模型重建步驟，其中體素處理步驟是利用移動立方體法對體數(shù)據(jù)的體素進行離散處理，生成三角面片集；模型記錄步驟是利用體數(shù)據(jù)的體素在模型中所處的行、層、分組位置與B+樹結(jié)構(gòu)的對應(yīng)關(guān)系，采用B+樹數(shù)據(jù)結(jié)構(gòu)方式，重新組織和記錄三維模型；模型重建步驟是利用B+樹結(jié)構(gòu)的數(shù)據(jù)搜索方式，根據(jù)B+樹結(jié)構(gòu)下記錄的三維模型的體素的連接關(guān)系，對體素及其中的三角面片定位，實現(xiàn)三維模型的重建。
可以看出，本發(fā)明對傳統(tǒng)移動立方體法作了改進，表現(xiàn)在，使用B+樹數(shù)據(jù)結(jié)構(gòu)分組存儲三角面片集。這樣，得到的整個髖關(guān)節(jié)三維模型為一個B+樹組成的森林樣數(shù)據(jù)結(jié)構(gòu)，B+樹是B樹的一種變形，其中所有的葉節(jié)點都相互連接，這樣就保證了三角面片之間的連接。B+樹提供了兩條搜索路徑，一條是由根至葉的查找，另一條是沿葉節(jié)點的查找，從而在模型重建時，明顯提高三角面片的定位速度，使三維重構(gòu)和分組可在對三維網(wǎng)格數(shù)據(jù)的一次掃描中完成。
參見圖3，所示為本發(fā)明骨科三維交互模擬手術(shù)的實現(xiàn)方法的流程示意圖。下面結(jié)合模擬小兒先天性髖脫位骨科手術(shù)過程的實例，對本方法進行詳細說明。
該方法包括骨骼模型建立步驟、骨刀模型建立步驟和手術(shù)模擬實施步驟。
其中，骨骼模型建立步驟是利用前述的未知體三維模型的建立方法，建立未知體三維模型，其模型圖像數(shù)據(jù)以B+樹數(shù)據(jù)結(jié)構(gòu)形式組織記錄。
選擇左側(cè)先天性髖脫位患兒1名，女，6歲，作為建模素材。采用PhilipTomoscan SR 7000型螺旋CT從髂骨翼至股骨近端連續(xù)掃描36層，螺旋層厚5mm，床進速度3mm/s，螺旋CT原始整體圖像數(shù)據(jù)通過HUB存儲于Dicom服務(wù)器。使用SGI圖形工作站，通過NFS共享方式，從Dicom服務(wù)器獲取螺旋CT原始圖像數(shù)據(jù)。
從上述原始圖像數(shù)據(jù)中選取適當灰度等級的數(shù)據(jù)開窗進行處理，提取出其中的模型圖像數(shù)據(jù)，并去除圖像中的干擾數(shù)據(jù)。
采用領(lǐng)域平均濾波器經(jīng)兩級重采樣，有效壓縮原始圖像數(shù)據(jù)。
采用前述改進的移動立方體法進行體素處理，并用B+樹方式分組存儲經(jīng)處理的體素及其所含三角面片信息。按照Marching Cubes方法可以發(fā)現(xiàn)其中隱含著三角面片的多層次包圍盒信息，由上向下依次為分組、層、行和體素，而用B+樹數(shù)據(jù)結(jié)構(gòu)記錄存儲三角面片符合移動立方體法處理過程的特點，這樣記錄的三維模型為一個B+數(shù)組成的森林樣數(shù)據(jù)結(jié)構(gòu)。
骨刀模型建立步驟，根據(jù)骨刀曲面方程，確定骨刀模型所需的物體空間大小，并將骨刀模型的物體空間劃分為與髖關(guān)節(jié)模型相同的網(wǎng)格空間；然后按層、行和體素順序，開始對每一邏輯立方體進行處理。當一個立方體不處于曲面邊界時，不處理該立方體，對處于邊界上的立方體則生成該立方體內(nèi)的局部數(shù)據(jù)場，并對該邏輯立方體內(nèi)生成的所有面片進行檢測，以保證該立方體內(nèi)屬于邊界外的面片不被包含入骨刀模型內(nèi)。
將骨刀模型離散為三角面片，使骨刀和髖關(guān)節(jié)模型物體空間的坐標原點重合，網(wǎng)格劃分和大小相同，根據(jù)曲面方程生成正規(guī)網(wǎng)格數(shù)據(jù)。再采用改進的移動立方體法進行體素處理，產(chǎn)生重構(gòu)骨刀模型的三角面片集，并將其組織到一棵B+樹中。
現(xiàn)分別以平面骨刀模型和球面骨刀模型為例，說明骨刀曲面方程的確定。平面骨刀模型，由局部坐標系點在全局坐標系中的坐標、刀平面方程和刀邊界形狀方程(4個平面方程)確定；球面骨刀模型由局部坐標系點在全局坐標系中的坐標、球面半徑和刀邊界形狀方程(2個平面方程)確定。
在手術(shù)模擬實施步驟中，包含骨刀模型的平移和/或旋轉(zhuǎn)步驟、未知體模型的切割步驟以及未知體模型的拼合步驟或其組合。實踐中的各種復(fù)雜術(shù)式都是通過這些基本操作或其組合實現(xiàn)的。
在手術(shù)模擬實施過程中，需要將骨刀模型移動到任意位置旋轉(zhuǎn)為任意角度后進行切割操作。平移和旋轉(zhuǎn)后骨刀模型的B+樹結(jié)構(gòu)必然徹底改變，而B+樹結(jié)構(gòu)又是實現(xiàn)前述快速切割算法的基礎(chǔ)。因此，如何處理作了平移和旋轉(zhuǎn)操作后的骨刀模型的B+樹是核心問題。一種思路是根據(jù)變換矩陣變換面片的所有頂點的坐標和法向，再由新面片集生成新的B+樹；另一種思路是根據(jù)變換矩陣生成整個骨刀模型的新描述信息，再由前述的骨刀模型生成方法生成新位置上的骨刀模型。按前一種算法，整個骨刀模型的每一個面片要被重新修改，而關(guān)鍵問題是變換過的面片很可能不再處于一個網(wǎng)格內(nèi)，因此必須將這樣的面片重新剖分。而按后一種算法，只要重新確定新位置上的骨刀模型經(jīng)過的網(wǎng)格，再處理這些網(wǎng)格就可以得到新模型。由于骨刀模型經(jīng)過的網(wǎng)格數(shù)較少，按第二種思路重新生成的面片比按第一種思路修改每一個面片的幾何構(gòu)造并進行剖分更為有效率。因此，本實施例中，采用了修改整個模型參數(shù)再重新生成骨刀模型的方式，先根據(jù)平移和/或旋轉(zhuǎn)變換矩陣，生成骨刀曲面方程的新描述信息，再建立該新的骨刀模型，實現(xiàn)骨刀的平移或旋轉(zhuǎn)。
骨刀模型移動到要求位置后，重新生成新的模型，以進行切割操作。此前我們已獲得了髖關(guān)節(jié)表面模型的B+樹表示，根據(jù)上面所述，骨刀模型是用與髖關(guān)節(jié)模型相同的B+樹結(jié)構(gòu)組織記錄三角面片集，兩模型的網(wǎng)格劃分方式及大小均相同。在為骨刀和髖關(guān)節(jié)模型建立了各自的B+樹表示后將兩模型放入同一物體空間，進行求交切割運算。
按照傳統(tǒng)的離散為三角面片的三維表面模型間的切割算法，需要對兩個表面模型逐個面片進行測試，以求得兩個模型的所有交點，再由交點開始利用面片的相鄰信息在被切割體的模型中跟蹤所有面片，確定每個面片應(yīng)屬于被分割后的哪一部分。為得到分割結(jié)果需要遍歷被切割模型的所有面片，對于復(fù)雜模型效率極低。
本發(fā)明實現(xiàn)切割的方法是在將骨刀和髖關(guān)節(jié)的三維表面模型放入同一物體空間，分別為兩個模型建立各自的B+樹表示后，利用B+樹中各層索引節(jié)點內(nèi)包含的層次包圍盒信息，通過對B+樹的不完全遍歷和回溯，只選擇子孫節(jié)點中包含交點的非葉節(jié)點進行擴展，達到快速定位表面模型間的交點并分割髖關(guān)節(jié)模型的目的。有效地提高了實際系統(tǒng)操作效率。
在其后的骨骼模型的拼合步驟中，是通過未知體模型移動部分的平移和/或旋轉(zhuǎn)，實現(xiàn)未知體模型的拼合。具體方式是根據(jù)變換矩陣變換未知體模型移動部分的所有面片的頂點的坐標和法向，再由得到新面片集生成新的B+樹，從而可以得到拼合后的骨骼模型。
本發(fā)明可具體應(yīng)用于微機環(huán)境，便于推廣使用，對提高整體醫(yī)療水平有重要意義，具有廣闊的應(yīng)用前景。
權(quán)利要求
1.人體骨骼未知體三維交互模型的建立方法，包括一圖像數(shù)據(jù)獲取步驟、一模型數(shù)據(jù)提取步驟及一模型數(shù)據(jù)重建步驟，其特征在于所述模型數(shù)據(jù)重建步驟包括一體素處理步驟，是利用移動立方體法對體數(shù)據(jù)的體素進行處理，生成三角面片集；一模型記錄步驟，是利用體數(shù)據(jù)的體素在模型中所處的行、層、分組位置與B+樹結(jié)構(gòu)的對應(yīng)關(guān)系，通過B+樹數(shù)據(jù)結(jié)構(gòu)方式，重新組織和記錄三維模型；一模型重建步驟，是利用B+樹結(jié)構(gòu)的數(shù)據(jù)搜索方式，根據(jù)B+樹結(jié)構(gòu)下記錄的三維模型的體素的連接關(guān)系，對體素及其中的三角面片定位，實現(xiàn)三維模型的重建。
2.如權(quán)利要求1所述的建立方法，其特征在于，所述圖像數(shù)據(jù)獲取步驟是采用螺旋CT對包含未知體的實體進行掃描，獲得原始圖像數(shù)據(jù)。
3.如權(quán)利要求2所述的建立方法，其特征在于，所述模型數(shù)據(jù)提取步驟是根據(jù)灰度等級對原始圖像數(shù)據(jù)進行選擇，對數(shù)據(jù)開窗進行處理，提取出其中的模型圖像數(shù)據(jù)，并去除圖像中的干擾數(shù)據(jù)。
4.如權(quán)利要求1所述的建立方法，其特征在于，在所述模型數(shù)據(jù)提取步驟后，還包括一模型數(shù)據(jù)壓縮步驟，是根據(jù)原始圖像數(shù)據(jù)的密度，對模型圖像數(shù)據(jù)進行兩級重采樣，以簡化數(shù)據(jù)，實現(xiàn)數(shù)據(jù)壓縮。
5.如權(quán)利要求5所述的建立方法，其特征在于，進行兩級重采樣所用的是領(lǐng)域平均濾波器。
6.人體骨性三維交互模擬手術(shù)的實現(xiàn)方法，其特征在于，該方法包括如下步驟骨胳模型建立步驟，利用以上任一權(quán)利要求所述的未知體三維交互模型的建立方法，建立骨骼三維交互模型，其模型圖像數(shù)據(jù)以B+樹數(shù)據(jù)結(jié)構(gòu)形式組織記錄；骨刀模型建立步驟，根據(jù)骨刀曲面方程，確定骨刀模型所需物體空間大小，并進行體素處理步驟，是利用移動立方體法對體數(shù)據(jù)的體素進行處理，生成三角面片集；模型記錄步驟，是利用體數(shù)據(jù)的體素在模型中所處的行、層、分組位置與B+樹結(jié)構(gòu)的對應(yīng)關(guān)系，通過B+樹數(shù)據(jù)結(jié)構(gòu)方式，重新組織和記錄三維模型；模型重建步驟，是利用B+樹結(jié)構(gòu)的數(shù)據(jù)搜索方式，根據(jù)B+樹結(jié)構(gòu)下記錄的三維模型的體素的連接關(guān)系，對體素及其中的三角面片定位，實現(xiàn)三維模型的重建，其中，該骨刀模型的網(wǎng)格劃分方式及大小與所述未知體三維模型相同；手術(shù)模擬步驟，包含骨刀模型的平移和/或旋轉(zhuǎn)步驟、骨骼模型的切割步驟以及骨骼模型的拼合步驟或其組合。
7.如權(quán)利要求6所述的實現(xiàn)方法，其特征在于，其中的手術(shù)模擬步驟中的骨刀的平移和/或旋轉(zhuǎn)步驟，是根據(jù)平移和/或旋轉(zhuǎn)變換矩陣，生成骨刀曲面方程的新描述信息，再建立該新的骨刀模型，實現(xiàn)骨刀的平移或旋轉(zhuǎn)。
8.如權(quán)利要求6所述的實現(xiàn)方法，其特征在于，其中的手術(shù)模擬步驟中的未知體模型切割步驟，是利用骨刀模型的B+樹數(shù)據(jù)結(jié)構(gòu)中各層索引節(jié)點內(nèi)包含的層次包圍盒信息，通過對骨骼模型的B+樹數(shù)據(jù)結(jié)構(gòu)的不完全遍歷和回溯，定位骨骼模型與骨刀模型表面間的節(jié)點并分割骨骼模型。
9.如權(quán)利要求6所述的實現(xiàn)方法，其特征在于，其中的手術(shù)模擬步驟中的骨骼模型的拼合步驟，是根據(jù)變換矩陣變換未知體模型移動部分的所有面片的頂點的坐標和法向，再由新面片集生成新的B+樹，通過骨骼模型移動部分的平移和/或旋轉(zhuǎn)，實現(xiàn)骨骼模型的拼合。
10.如權(quán)利要求6所述的實現(xiàn)方法，其特征在于，所述骨刀模型包括平面骨刀模型和/或球面骨刀模型，其中，平面骨刀模型由以下參數(shù)確定局部坐標系點在全局坐標系中的坐標、刀平面方程和刀邊界形狀方程；球面骨刀模型由以下參數(shù)確定局部坐標系點在全局坐標系中的坐標、球面半徑和刀邊界形狀方程。
全文摘要
本發(fā)明公開了一種人體骨骼未知體三維交互模型的建立方法，包括圖像數(shù)據(jù)獲取步驟、模型數(shù)據(jù)提取步驟及模型數(shù)據(jù)重建步驟，模型數(shù)據(jù)重建步驟包括體素處理步驟，利用移動立方體法對體數(shù)據(jù)的體素進行處理，生成三角面片集；模型記錄步驟，利用體數(shù)據(jù)的體素在模型中所處的行、層、分組位置與B+樹結(jié)構(gòu)的對應(yīng)關(guān)系，通過B+樹數(shù)據(jù)結(jié)構(gòu)方式，重新組織和記錄三維模型；模型重建步驟，利用B+樹結(jié)構(gòu)的數(shù)據(jù)搜索方式，根據(jù)B+樹結(jié)構(gòu)下記錄的三維模型的體素的連接關(guān)系，對體素及其中的三角面片定位，實現(xiàn)三維模型的重建。本發(fā)明還公開了一種應(yīng)用上述方法實現(xiàn)人體骨性三維交互模擬手術(shù)的方法，能實現(xiàn)實時交互的手術(shù)模擬，從而為制定最佳手術(shù)方案提供依據(jù)。
文檔編號G06T17/50GK1555030SQ20031012115
公開日2004年12月15日 申請日期2003年12月22日 優(yōu)先權(quán)日2003年12月22日
發(fā)明者李浩宇, 許忠信 申請人:李浩宇