專利名稱：一種醫(yī)學(xué)圖像目標(biāo)區(qū)域勾畫方法
技術(shù)領(lǐng)域：
本發(fā)明涉及圖像處理，特別是利用計算機對4D圖像中的目標(biāo)區(qū)域進行勾畫的方法。
背景技術(shù)：
目標(biāo)靶區(qū)勾畫是放療技術(shù)中極為重要的步驟，該步的準(zhǔn)確度直接決定了放療計劃的精確程度。目前臨床中都是由醫(yī)師手動對目標(biāo)進行勾畫來制定放療計劃。而4D-CT圖像實際是由一組各個不同時相的常規(guī)CT圖像組成的，圖像數(shù)目眾多，通常約為1000 2000 張。巨大的圖像數(shù)目極大地加重了醫(yī)師勾畫的工作量，這種重復(fù)勾畫不僅浪費了寶貴的人力和精力，而且容易引起醫(yī)師疲勞而產(chǎn)生人為失誤等。實現(xiàn)目標(biāo)靶區(qū)邊界的自動勾畫，最直觀的思路應(yīng)該是直接對每個呼吸時相的圖像采用圖像分割技術(shù)來自動提取出目標(biāo)靶區(qū)。區(qū)域生長法是目前圖像分割中比較常用的方法，需要首先選取出一個或多個像素點作為種子點，然后逐一對種子像素周圍鄰域內(nèi)的像素根據(jù)事先確定的生長準(zhǔn)則進行判定，若滿足，則將該像素合并到種子像素集中，將這些新像素當(dāng)作新的種子像素繼續(xù)進行上面操作過程，直到再沒有滿足條件的像素點可被包括進來。然而種子點的選取一般都是人工手動完成的，而對4D-CT的大量圖像來說工作量仍然很大。由于呼吸運動屬于半周期運動，具有一定的規(guī)律性，而由此引起的目標(biāo)運動和變形也應(yīng)服從一定的統(tǒng)計規(guī)律。因此，也有學(xué)者提出在一些先驗知識的基礎(chǔ)上建立一個概率模型來描述目標(biāo)靶區(qū)的運動和變形，從而實現(xiàn)對其輪廓的自動跟蹤。常用的概率模型有貝葉斯濾波模型(bayesian filter)、卡爾曼濾波模型(kalman filter)、粒子濾波模型 (particle filter)等。然而建立一個概率模型需要大量的訓(xùn)練數(shù)據(jù)，而呼吸運動較大的個體差異性使得模型建立難度進一步加大。并且，當(dāng)出現(xiàn)不規(guī)則呼吸運動時，很難再用經(jīng)典的統(tǒng)計模型來描述人體內(nèi)的器官運動和變形。近年來，基于圖像配準(zhǔn)的輪廓映射方法成為目標(biāo)跟蹤和分割問題的比較流行的解決方法，其思路是首先在參考圖像中手動勾畫目標(biāo)輪廓，然后將輪廓點映射到其他待處理圖像上的對應(yīng)位置。該類方法的關(guān)鍵在于能夠得到準(zhǔn)確的點點間的映射關(guān)系，該映射關(guān)系通常都是采用變形配準(zhǔn)來實現(xiàn)的。其缺點在于直接對整幅圖像進行變形配準(zhǔn)，因此計算量很大，而輪廓映射的準(zhǔn)確度也會受到離目標(biāo)輪廓較遠(yuǎn)區(qū)域的影響。

發(fā)明內(nèi)容
本發(fā)明的目的是提供一種醫(yī)學(xué)圖像目標(biāo)區(qū)域勾畫方法，以快速、準(zhǔn)確地從醫(yī)學(xué)圖像中勾畫出目標(biāo)區(qū)域，大大減輕工作人員在圖像分割中的工作量。本發(fā)明的總體思路是手動勾畫出參考時相上的目標(biāo)輪廓，而在其他時相上以該目標(biāo)輪廓為參考，通過尋找與輪廓控制區(qū)域相似的圖像小塊來實現(xiàn)目標(biāo)靶區(qū)的自動跟蹤， 最后通過在跟蹤得到的每兩個相鄰輪廓控制區(qū)域中心之間進行非線性插值實現(xiàn)對其他時相靶區(qū)邊界的自動勾畫。進一步，可在輪廓跟蹤過程中，通過引入一個約束項來保證跟蹤到的邊界的完整性和連貫性。本發(fā)明的技術(shù)方案如下一種醫(yī)學(xué)圖像目標(biāo)區(qū)域勾畫方法，包括以下步驟(1)對圖像進行去噪及增強邊界特征的預(yù)處理；(2)在參考時相上勾畫出目標(biāo)輪廓，并根據(jù)該目標(biāo)輪廓選取若干輪廓控制區(qū)域；(3)分別提取所述若干輪廓控制區(qū)域的圖像紋理特征；(4) 根據(jù)圖像紋理特征，在待處理時相上搜索，跟蹤參考時相上的每個輪廓控制區(qū)域在待處理時相上的對應(yīng)位置；(5)在待處理時相上跟蹤得到的每相鄰兩個輪廓控制區(qū)域的中心之間進行插值處理，從而完成待處理時相上目標(biāo)區(qū)域的自動勾畫。一種較佳方案中，步驟(1)采用各向異性濾波算法對圖像進行預(yù)處理。一種較佳方案中，步驟(2)通過以下方法選取輪廓控制區(qū)域把目標(biāo)輪廓作為一個多邊形，分別以多邊形的每個頂點為中心構(gòu)建輪廓控制區(qū)域。更好的一種方案是，當(dāng)多邊形的兩個相鄰頂點之間的距離大于設(shè)定的閾值A(chǔ)時，在該兩個相鄰頂點之間的目標(biāo)輪廓上再選取至少一個像素點作為中心構(gòu)建輪廓控制區(qū)域。一種較佳方案中，步驟(3)采用SObel算子提取輪廓控制區(qū)域的圖像紋理特征。一種較佳方案中，步驟(4)根據(jù)醫(yī)學(xué)圖像對應(yīng)的器官運動幅度的先驗知識，限定待處理時相上每個輪廓控制區(qū)域的搜索范圍，該搜索范圍的中心位置與參考時相上的相應(yīng)輪廓控制區(qū)域的中心位置對應(yīng)。步驟(4)中，可以通過以下方案跟蹤參考時相上的每個輪廓控制區(qū)域在待處理時
相上的對應(yīng)位置通過公式
權(quán)利要求
1.一種醫(yī)學(xué)圖像目標(biāo)區(qū)域勾畫方法，其特征在于，包括以下步驟(1)對圖像進行去噪及增強邊界特征的預(yù)處理；(2)在參考時相上勾畫出目標(biāo)輪廓，并根據(jù)該目標(biāo)輪廓選取若干輪廓控制區(qū)域；(3)分別提取所述若干輪廓控制區(qū)域的圖像紋理特征；(4)根據(jù)圖像紋理特征，在待處理時相上搜索，跟蹤參考時相上的每個輪廓控制區(qū)域在待處理時相上的對應(yīng)位置；(5)在待處理時相上跟蹤得到的每相鄰兩個輪廓控制區(qū)域的中心之間進行插值處理，從而完成待處理時相上目標(biāo)區(qū)域的自動勾畫。
2.根據(jù)權(quán)利要求1所述的目標(biāo)區(qū)域勾畫方法，其特征在于步驟(1)中采用各向異性濾波算法對圖像進行預(yù)處理。
3.根據(jù)權(quán)利要求1所述的目標(biāo)區(qū)域勾畫方法，其特征在于步驟(2)中通過以下方法選取輪廓控制區(qū)域把目標(biāo)輪廓作為一個多邊形，分別以多邊形的每個頂點為中心構(gòu)建輪廓控制區(qū)域。
4.根據(jù)權(quán)利要求3所述的目標(biāo)區(qū)域勾畫方法，其特征在于步驟(2)中，當(dāng)多邊形的兩個相鄰頂點之間的距離大于設(shè)定的閾值A(chǔ)時，在該兩個相鄰頂點之間的目標(biāo)輪廓上再選取至少一個像素點作為中心構(gòu)建輪廓控制區(qū)域。
5.根據(jù)權(quán)利要求1所述的目標(biāo)區(qū)域勾畫方法，其特征在于步驟(3)中采用sobel算子提取輪廓控制區(qū)域的圖像紋理特征。
6.根據(jù)權(quán)利要求1所述的目標(biāo)區(qū)域勾畫方法，其特征在于步驟(4)中，根據(jù)醫(yī)學(xué)圖像對應(yīng)的器官運動幅度的先驗知識，限定待處理時相上每個輪廓控制區(qū)域的搜索范圍，該搜索范圍的中心位置與參考時相上的相應(yīng)輪廓控制區(qū)域的中心位置對應(yīng)。
7.根據(jù)權(quán)利要求6所述的目標(biāo)區(qū)域勾畫方法，其特征在于步驟(4)中，通過公式^ = ^^( .,,- )2計算區(qū)域特征向量之間的歐式距離來衡量圖像的相似性，其中，表示參考時相上第i個輪廓控制區(qū)域特征向量的第A個元素，&〃表示待處理時相上第i個輪廓控制區(qū)域的對應(yīng)搜索范圍內(nèi)的第J個被搜索區(qū)域的第A個元素，#表示特征向量的長度，SLj越小，參考時相上的對應(yīng)輪廓控制區(qū)域與待處理時相上的被搜索區(qū)域之間的相似性越大，將相似性最大的被搜索區(qū)域作為參考時相上的對應(yīng)輪廓控制區(qū)域在待處理時相上的對應(yīng)位置。
8.根據(jù)權(quán)利要求6所述的目標(biāo)區(qū)域勾畫方法，其特征在于步驟(4)中，通過公式計算區(qū)域特征向量之間的歐式距離來衡量圖像的相似性，其中，表示參考時相上第i個輪廓控制區(qū)域特征向量的第A個元素，&〃表示待處理時相上第i個輪廓控制區(qū)域的對應(yīng)搜索范圍內(nèi)的第J個被搜索區(qū)域的第A個元素，#表示特征向量的長度，SLj越小，參考時相上的對應(yīng)輪廓控制區(qū)域與待處理時相上的被搜索區(qū)域之間的相似性越大，同時通過公式Cu =Ici-C^計算相鄰輪廓控制區(qū)域相關(guān)性的變化程度，其中，Cri表示參考時相上第i個輪廓控制區(qū)域與待處理時相上對應(yīng)的相鄰輪廓控制區(qū)域之間的相關(guān)系數(shù)，以u表示待處理時相上搜索范圍內(nèi)第i個輪廓控制區(qū)域的第J·個被搜索區(qū)域與參考時相上對應(yīng)相鄰輪廓控制區(qū)域之間的相關(guān)系數(shù)，然后通過公式= ^gmmj-Bij- = argmm/5^ +^Cij)確定參考時相上的對應(yīng)輪廓控制區(qū)域在待處理時相上的對應(yīng)位置，其中μ為相似性項SLj與約束項CLj之間的相對權(quán)重。
全文摘要
一種醫(yī)學(xué)圖像目標(biāo)區(qū)域勾畫方法，包括以下步驟(1)預(yù)處理；(2)在參考時相上勾畫出目標(biāo)輪廓，并根據(jù)該目標(biāo)輪廓選取若干輪廓控制區(qū)域；(3)分別提取所述若干輪廓控制區(qū)域的圖像紋理特征；(4)根據(jù)圖像紋理特征，跟蹤參考時相上的每個輪廓控制區(qū)域在待處理時相上的對應(yīng)位置；(5)在待處理時相上跟蹤得到的每相鄰兩個輪廓控制區(qū)域的中心之間進行插值處理，從而完成待處理時相上目標(biāo)區(qū)域的自動勾畫。本發(fā)明勾畫方法能夠根據(jù)參考時相上的目標(biāo)輪廓和輪廓控制區(qū)域，快速、準(zhǔn)確地自動勾畫出其它時相上的目標(biāo)區(qū)域，能夠大大減輕工作人員在圖像分割中的工作量。
文檔編號G06T7/00GK102509286SQ20111030208
公開日2012年6月20日 申請日期2011年9月28日 優(yōu)先權(quán)日2011年9月28日
發(fā)明者李哲, 田珍, 袁克虹 申請人:清華大學(xué)深圳研究生院