專利名稱:用于處理醫(yī)用圖像的方法和設(shè)備的制作方法
技術(shù)領(lǐng)域:
本發(fā)明涉及用于處理醫(yī)用圖像的方法和對(duì)應(yīng)的設(shè)備。它尤其適于但不限于處理磁共振圖像,例如人腦的磁共振圖像。
背景技術(shù):
在神經(jīng)成像領(lǐng)域,從腦部的磁共振圖像(MRI)自動(dòng)抽取特征是一種越來越漸重要的方法。因?yàn)槟X部解剖學(xué)在受試者之間變化很顯著,并且在衰老或疾病進(jìn)展期間可以發(fā)生顯著的改變,所以近年來找到在特征抽取期間處理解剖學(xué)差異的適當(dāng)方式已經(jīng)獲得了日漸增長(zhǎng)的關(guān)注?;趫D冊(cè)的方法是用于處理此可變性的最流行的方法之一。在本文的語(yǔ)境中,“圖冊(cè)(atlas)”是一種數(shù)據(jù)集(其可以是3D圖像、2D圖像、任何維的圖像或一組圖像),其具有注釋或標(biāo)記以便識(shí)別圖像內(nèi)的點(diǎn)、區(qū)域或結(jié)構(gòu)?;趫D冊(cè)的方法假定圖冊(cè)能夠以概率或統(tǒng)計(jì)學(xué)方式編碼解剖學(xué)可變性。當(dāng)構(gòu)建代表性的圖冊(cè)時(shí),重要的是使所有圖像與不偏向總體中任何特定子群的模板配準(zhǔn)。Avants和Gee (2004)以及Joshi等Q004)已經(jīng)提出了兩種使用用于形狀平均和圖冊(cè)構(gòu)建的大形變微分同胚(diffeomorphic)背景的方法。用于共配準(zhǔn)圖像的無模板方法為空間圖像歸一化形成確定的框架。與尋找單個(gè)代表性平均圖冊(cè)的方式相背,兩種更近期的方法描述了識(shí)別圖像數(shù)據(jù)集中不同總體的模式的方式(Blezek和Miller,2007 ;Sabuncu等,2008)。為設(shè)計(jì)取決于受試者信息的可變圖冊(cè),近年來已經(jīng)將多種方法應(yīng)用于表征隨時(shí)間變化和疾病進(jìn)展期間腦部形狀的解剖學(xué)變化的問題。Davis等Q007)描述了總體形狀回歸的方法,其中使核回歸適應(yīng)于微分同胚的流形并且將其用于獲得年齡依賴性圖冊(cè)。Ericsson等Q008)提出了構(gòu)建患者特異性圖冊(cè)的方法,其中不同的平均腦圖冊(cè)被構(gòu)建在根據(jù)元信息如性別、年齡或臨床因素的小形變背景中。從磁共振(MR)腦圖像數(shù)據(jù)中抽取特征或生物標(biāo)志的方法通常開始于自動(dòng)分割目的區(qū)域。非常流行的分割方法是使用標(biāo)記傳播,標(biāo)記傳播通過對(duì)兩種圖像進(jìn)行比對(duì)來將標(biāo)記由圖冊(cè)圖像轉(zhuǎn)換到未見的目標(biāo)圖像。典型但不是必要地,圖冊(cè)被手動(dòng)標(biāo)記。Bajcsy等(1983)以及最近的Gee等(199 和Collins等(19卯)提出了使用此方法的早期工作。標(biāo)記傳播的精確度很大程度上取決于基礎(chǔ)的圖像比對(duì)的精確度。為了克服對(duì)單一分割的依賴,Warfield等Q004)提出了 STAPLE方法,該方法對(duì)分割集合計(jì)算真實(shí)分割的概率估計(jì)值算。Rohlfing等Q004)證明了多分類器框架的提高的穩(wěn)健性和精確度,在所述框架中從多圖冊(cè)傳播的標(biāo)記在分類器融合步驟中結(jié)合從而獲得對(duì)目標(biāo)圖像的最終分割。Heckemarm等Q006)成功地將標(biāo)記傳播與分類器融合結(jié)合用于分割腦部MR圖像中的大量結(jié)構(gòu)。由于解剖學(xué)變化具有寬的范圍,在多圖冊(cè)分割中對(duì)圖冊(cè)的選擇成為重要的問題。對(duì)給定目標(biāo)的合適圖冊(cè)的選擇有助于保證圖冊(cè)-目標(biāo)配準(zhǔn)和隨后的分割盡可能的精確?!龅萉007)描述了通過合并圖冊(cè)選擇來改善單個(gè)圖冊(cè)中分割結(jié)果的不同方法。Aljabar等Q009)研究了在多圖冊(cè)分割期間最優(yōu)圖冊(cè)選擇的不同的相似性度量。Van Rikxoort等(2008)提出了這樣的方法,其中在圖像的不同子窗口中分別進(jìn)行圖冊(cè)組合直到滿足收斂性判據(jù)。這些方法顯示有意義的是為每個(gè)目標(biāo)圖像單獨(dú)選擇合適的圖冊(cè)。盡管MR腦部圖像數(shù)量的增加是有用的,但是高質(zhì)量手動(dòng)圖冊(cè)的產(chǎn)生是勞動(dòng)密集并且成本高的任務(wù)(參見例如Hammers等Q003))。這意味著圖冊(cè)的數(shù)目通常相對(duì)有限,并且在大多數(shù)情況中,限于特定的總體(例如年輕人、健康受試者)。即使使用選擇方式,這也會(huì)限制圖冊(cè)數(shù)據(jù)庫(kù)的適用性。為克服這一點(diǎn),Tang等Q009)設(shè)法通過利用自單一模板圖像和練習(xí)圖像之間的轉(zhuǎn)換學(xué)習(xí)的形變的PCA模型來產(chǎn)生多種圖冊(cè)圖像。通過用取樣自該模型的若干轉(zhuǎn)換對(duì)初始模板進(jìn)行轉(zhuǎn)換來產(chǎn)生可能的圖冊(cè)。假設(shè)是,通過發(fā)現(xiàn)未見圖像的合適圖冊(cè),可以容易地獲得對(duì)此模板快速且精確的配準(zhǔn)。然而,具有比練習(xí)數(shù)據(jù)更高的變異水平的測(cè)試數(shù)據(jù)將對(duì)此方式提出重大挑戰(zhàn)。此外,對(duì)高度可變練習(xí)數(shù)據(jù)集的使用可能導(dǎo)致無代表性的PCA模型以及不同圖像和單個(gè)模板間的配準(zhǔn)誤差增加的可能性。此限制使得此方式僅適用于可以容易地獲得從所有練習(xí)圖像至單個(gè)初始模板的良好配準(zhǔn)的情況。基于圖冊(cè)的分割受益于對(duì)與目標(biāo)圖像相似的圖冊(cè)的選擇OVu等,2007 ;Aljabar等,2009)。然而,實(shí)際上,初始圖冊(cè)可能僅表示目標(biāo)圖像總體的特定子群。因此需要能夠?qū)⑾鄬?duì)少量的圖冊(cè)傳播到顯示大量解剖學(xué)變化的大的并且不同組的MR腦部圖像。相對(duì)于直接多圖冊(cè)傳播,將自動(dòng)標(biāo)記的腦部圖像用于標(biāo)記未見圖像的現(xiàn)有技術(shù)不能提高分割精確度。在(Heckemarm等,2006)中,當(dāng)多個(gè)相對(duì)一致的圖冊(cè)傳播到被當(dāng)做用于分割未見圖像的單個(gè)圖冊(cè)使用的隨機(jī)選擇的中間圖像時(shí),相比直接多圖冊(cè)傳播和融合的0:84,所得的采用手工描繪(delineation)的平均Dice重疊是080。在另一個(gè)實(shí)驗(yàn)中,單個(gè)圖冊(cè)被傳播到之后進(jìn)一步用于多圖冊(cè)分割的隨機(jī)選擇的中間對(duì)象,產(chǎn)生最好0:78的與手工描繪的Dice重疊。更多的背景技術(shù)由US2007/0053589 AUUS 2008/0154118 Al和WO 2009/093146Al提供,所有上述文獻(xiàn)都公開了用于分割圖像數(shù)據(jù)的方法。發(fā)明概述根據(jù)本發(fā)明的第一方面,提供了在所附權(quán)利要求的權(quán)利要求1中所述的方法。因此提供了由計(jì)算機(jī)處理器進(jìn)行的處理醫(yī)用圖像的方法,并且所述方法包括以下步驟(a)獲得包含一個(gè)或多個(gè)圖像的一個(gè)或多個(gè)圖冊(cè),在所述一個(gè)或多個(gè)圖像中已經(jīng)用標(biāo)記數(shù)據(jù)標(biāo)記了一個(gè)或多個(gè)解剖學(xué)特征;(b)獲得多個(gè)未標(biāo)記的圖像;(c)將所述標(biāo)記的和未標(biāo)記的圖像進(jìn)行比較并且選擇與一個(gè)或多個(gè)所述標(biāo)記的圖像最相似的一個(gè)或多個(gè)未標(biāo)記的圖像;(d)使標(biāo)記數(shù)據(jù)從一個(gè)或多個(gè)最接近的標(biāo)記圖像傳播到那些選定的一個(gè)或多個(gè)圖像中的每個(gè)圖像,由此標(biāo)記所述選定的一個(gè)或多個(gè)圖像中的每個(gè)圖像的對(duì)應(yīng)的一個(gè)或多個(gè)解剖學(xué)特征并且使所述選定的一個(gè)或多個(gè)圖像成為標(biāo)記的一個(gè)或多個(gè)圖像;以及(e)從步驟(c)開始反復(fù)地重復(fù),由此標(biāo)記所述未標(biāo)記的圖像中的其他圖像。術(shù)語(yǔ)“標(biāo)記的”應(yīng)當(dāng)被廣義地解釋從而包括對(duì)解剖學(xué)特征的任何一種描繪、分割或注釋。類似地,術(shù)語(yǔ)“標(biāo)記數(shù)據(jù)”應(yīng)當(dāng)被廣義地解釋從而包括使解剖學(xué)特征在醫(yī)學(xué)圖像上得到描繪、分割或注釋的任何一種編碼。通過使標(biāo)記數(shù)據(jù)反復(fù)地從最接近的標(biāo)記的圖像傳播到未標(biāo)記的圖像,可以使用結(jié)構(gòu)上相似的圖冊(cè)分割每個(gè)未標(biāo)記的圖像。結(jié)果,可以使標(biāo)記的圖像和未標(biāo)記的圖像之間相對(duì)大的差異分解為比較相似的初始未標(biāo)記圖像之間許多的差異,通過上述方式標(biāo)記數(shù)據(jù)得到傳播,并且使得配準(zhǔn)誤差得到降低。在從屬權(quán)利要求中限定了優(yōu)選的、任選的特征。因此,優(yōu)選地,比較標(biāo)記和未標(biāo)記的圖像的步驟包括將圖像嵌入低維坐標(biāo)系統(tǒng)中。這使得標(biāo)記的和未標(biāo)記的圖像以計(jì)算上有效的方式進(jìn)行比較并且使得所述差異以計(jì)算上有效的方式進(jìn)行定量評(píng)價(jià)。在某些實(shí)施方案中,低維坐標(biāo)系統(tǒng)可以是二維坐標(biāo)空間,因此使對(duì)圖像之間差異的分析和處理進(jìn)一步簡(jiǎn)化。優(yōu)選地,比較標(biāo)記的和未標(biāo)記的圖像的步驟包括通過比較一組圖像中每對(duì)圖像的一個(gè)或多個(gè)各自的解剖學(xué)特征來限定一組成對(duì)的相似性度量。尤其優(yōu)選地,盡管本領(lǐng)域技術(shù)人員將理解存在可以實(shí)現(xiàn)此步驟的其他方式,此步驟還包括對(duì)成對(duì)的相似性度量進(jìn)行譜分析操作。成對(duì)的相似性度量可以表示一對(duì)圖像之間的強(qiáng)度相似性,和/或一對(duì)圖像之間形
變的量。優(yōu)選地,傳播標(biāo)記數(shù)據(jù)的步驟包括基于分類器融合技術(shù)使標(biāo)記數(shù)據(jù)從多個(gè)最接近的標(biāo)記的圖像傳播。這使得一個(gè)或多個(gè)選定的圖像以更高的精度進(jìn)行標(biāo)記。優(yōu)選地,該方法在步驟(d)后以及步驟(e)前還包括對(duì)新近傳播的標(biāo)記數(shù)據(jù)進(jìn)行基于強(qiáng)度的精修操作的步驟,以便進(jìn)一步使標(biāo)記過程期間的配準(zhǔn)誤差的累積最小化。圖像可以是不同受試者的。備選地,至少一些圖像可以是相同受試者的但是在不同的時(shí)間點(diǎn)獲取的圖像,由此使得受試者自身變化(intra-subject variance)得到識(shí)別和研究。圖像可以是磁共振圖像,或本領(lǐng)域技術(shù)人員熟悉的其他醫(yī)用圖像。該方法還可以包括對(duì)表示病況出現(xiàn)與否的解剖學(xué)特征進(jìn)行標(biāo)記并且使用該特征得到該病況的生物標(biāo)志。在生物標(biāo)志的基礎(chǔ)上,該方法還可以包括將受試者分配以診斷分類,和/或量化受試者對(duì)治療的反應(yīng),和/或選擇受試者的治療。根據(jù)本發(fā)明的第二方面,提供了成像設(shè)備,所述成像設(shè)備被配置成執(zhí)行根據(jù)本發(fā)明的第一實(shí)施方案的方法。該成像設(shè)備可以是醫(yī)用掃描儀,如MRI掃描儀或其他類型。根據(jù)本發(fā)明的第三方面,提供了圖像處理設(shè)備,所述圖像處理設(shè)備被配置成執(zhí)行根據(jù)本發(fā)明的第一實(shí)施方案的方法。根據(jù)本發(fā)明的第四方面,提供了計(jì)算機(jī)系統(tǒng),所述計(jì)算機(jī)系統(tǒng)被配置成執(zhí)行根據(jù)本發(fā)明的第一實(shí)施方案的方法。根據(jù)本發(fā)明的第五方面,提供了計(jì)算機(jī)程序,所述計(jì)算機(jī)程序包括用于執(zhí)行根據(jù)本發(fā)明的第一實(shí)施方案的方法的編碼的指令。根據(jù)本發(fā)明的第六方面,提供了編碼根據(jù)本發(fā)明第五實(shí)施方案的計(jì)算機(jī)程序的計(jì)算機(jī)可讀介質(zhì)或物理載波信號(hào)。附圖簡(jiǎn)述現(xiàn)在將僅通過實(shí)例并且參照附圖來描述本發(fā)明的實(shí)施方案,在附圖中
圖1圖示使用我們的新方法的圖冊(cè)傳播過程;圖2圖示結(jié)果,該結(jié)果顯示對(duì)在四個(gè)受試者組(健康的年輕人、老年人對(duì)照、MCI、AD)之中所選擇的不同特征維度的辨別能力;圖3圖示MNI152腦部圖冊(cè),其顯示用于評(píng)價(jià)成對(duì)圖像相似性的海馬周圍的目的區(qū)域;圖4圖示坐標(biāo),該坐標(biāo)嵌入有基于健康受試者的30個(gè)圖冊(cè)和來自老年癡呆癥患者和年齡匹配的對(duì)照受試者的796個(gè)圖像;圖5圖示對(duì)橫切切片上右側(cè)海馬的分割結(jié)果的比較;圖6圖示隨著距原始的圖冊(cè)集的距離增加分割精確度的變化情況,用于評(píng)估的每個(gè)圖像亞組由一個(gè)條形圖表示;圖7圖示手動(dòng)和自動(dòng)分割的平均海馬體積;以及圖8是Bland-Altman圖,其顯示基于對(duì)海馬的手動(dòng)和自動(dòng)分割的體積度量之間的一致。優(yōu)選實(shí)施方案描述這些實(shí)施方案表示為申請(qǐng)人所知的實(shí)施本發(fā)明的最佳方式。然而,它們不是可以實(shí)現(xiàn)此目的的唯一方式。首先,本實(shí)施方案采取用于處理醫(yī)用(或其他)圖像的方法或算法的形式。所述方法或算法可以被集成到能夠通過計(jì)算機(jī)處理器執(zhí)行的計(jì)算機(jī)程序或指令代碼集中。計(jì)算機(jī)處理器可以是常規(guī)(足夠高的性能)計(jì)算機(jī)的處理器,或者其他一些圖像處理設(shè)備或計(jì)算機(jī)系統(tǒng)的處理器。備選地,計(jì)算機(jī)處理器可以集成至一臺(tái)醫(yī)用成像設(shè)備如MRI掃描儀中,或者與一臺(tái)醫(yī)用成像設(shè)備如MRI掃描儀通信??梢栽谟?jì)算機(jī)可讀介質(zhì)或數(shù)據(jù)載體如CD-ROM、DVD或固態(tài)存儲(chǔ)裝置上提供計(jì)算機(jī)程序或指令代碼集。備選地,它可以作為數(shù)據(jù)信號(hào)從相連的計(jì)算機(jī)或通過局域網(wǎng)或廣域網(wǎng)如因特網(wǎng)下載。作為進(jìn)一步的備選方案,可以將計(jì)算機(jī)程序或指令代碼集硬編碼在被配置成執(zhí)行它的計(jì)算機(jī)處理器(或與其相關(guān)聯(lián)的存儲(chǔ)器)中。初步概述我們的方法以獲得一個(gè)或多個(gè)預(yù)先存在的圖冊(cè)開始,在所述一個(gè)或多個(gè)預(yù)先存在的圖冊(cè)中一組數(shù)字圖像已經(jīng)被標(biāo)記或注釋。從例如MRI掃描儀或另一臺(tái)醫(yī)用成像設(shè)備同樣獲得其上傳播有標(biāo)記或注釋的一組圖像。要研究的圖像可以是腦部圖像。備選地它們可以是人(或動(dòng)物)身體的其他部分的圖像,如膝蓋的圖像-例如為了診斷骨關(guān)節(jié)炎。在圖1中描繪了使用我們的新方法的圖冊(cè)傳播和分割過程,其顯示5個(gè)步驟。首先,在步驟(1)中,將所有標(biāo)記的圖像(即圖冊(cè))和未標(biāo)記的圖像嵌入低維流形(low-dimensional manifold)中。在步驟O)中,選擇N個(gè)與標(biāo)記的圖像最接近的未標(biāo)記的圖像用于分割。然后,在步驟⑶中,使M個(gè)最接近的標(biāo)記的圖像與選定的圖像(圖示了一個(gè)選定的圖像的實(shí)例)中的每個(gè)圖像配準(zhǔn)。在步驟中,強(qiáng)度精修用來獲得每個(gè)選定的圖像的標(biāo)記圖。然后,在步驟(5)中,重復(fù)步驟(幻_(4)直到更多的圖像(并且優(yōu)選地所有圖像)被標(biāo)記。如前所述,基于圖冊(cè)的分割受益于對(duì)與目標(biāo)圖像相似的圖冊(cè)的選擇。我們的方法提供這樣的框架,其中這通過首先將所有圖像嵌入低維坐標(biāo)系統(tǒng)得到保證,所述低維坐標(biāo)系統(tǒng)提供圖像之間的距離度量并且允許識(shí)別圖像的鄰近區(qū)。在由坐標(biāo)系統(tǒng)嵌入學(xué)習(xí)得到的流形中,傳播框架可以得到識(shí)別并且標(biāo)記的圖冊(cè)可以以逐步的方式得到傳播,即從初始圖冊(cè)開始到整個(gè)總體被分割為止。使用在其鄰近區(qū)內(nèi)的圖冊(cè)來分割各個(gè)圖像,這意味著相異圖像之間的形變被分解為比較相似的圖像之間若干小的形變并且配準(zhǔn)誤差得到降低。為了進(jìn)一步使配準(zhǔn)誤差的累積最小化,在每個(gè)標(biāo)記傳播步驟后完成對(duì)分割的基于強(qiáng)度的精修。一旦被分割,圖像又可以在隨后的分割步驟中被用作圖冊(cè)。在總體中的所有圖像被分割后,它們表示大的圖冊(cè)數(shù)據(jù)庫(kù),從所述圖冊(cè)數(shù)據(jù)庫(kù)中可以選擇用于分割未見圖像的合適亞組。通過對(duì)它們的成對(duì)相似性應(yīng)用譜分析步驟來獲得其中嵌入有圖像的坐標(biāo)系統(tǒng)。當(dāng)對(duì)特定目標(biāo)圖像進(jìn)行標(biāo)記的圖冊(cè)的傳播和融合時(shí),使它們提供的信息與基于目標(biāo)圖像強(qiáng)度的模型結(jié)合從而產(chǎn)生最終分割。因此,為了通過具有高水平的受試者自身變化的圖像數(shù)據(jù)集傳播初始集合的圖冊(cè),學(xué)習(xí)得到表示該數(shù)據(jù)集的流形,其中局部鄰近區(qū)內(nèi)的圖像彼此相似。該流形由所有圖像的嵌入坐標(biāo)表示。通過將其中各個(gè)頂點(diǎn)表示圖像的完整圖應(yīng)用譜分析步驟來獲得此嵌入,并且圖像間的所有成對(duì)相似性被用于限定該圖中的邊權(quán)。成對(duì)相似性可以測(cè)量為圖像間的強(qiáng)度相似性或圖像間的形變的量或兩者的組合。在隨后的步驟中,在新限定的坐標(biāo)系統(tǒng)內(nèi)傳播圖冊(cè)。在第一步中,初始集合的圖冊(cè)被傳播到它們的局部鄰近區(qū)中的多個(gè)圖像并且用于標(biāo)記它們。以此方式標(biāo)記的圖像自身變成圖冊(cè),并且在隨后的步驟中,在整個(gè)數(shù)據(jù)集內(nèi)進(jìn)一步傳播。以此方式,各個(gè)圖像通過使用緊鄰它的多個(gè)圖冊(cè)進(jìn)行標(biāo)記,這具有使配準(zhǔn)誤差降低的好處。在此技術(shù)的擴(kuò)展中,在不同時(shí)間從相同受試者獲得的一個(gè)或多個(gè)掃描(所謂的“縱向”掃描)可以得到標(biāo)記。在將多個(gè)圖冊(cè)傳播到各個(gè)基線掃描后,從多個(gè)圖冊(cè)獲得的空間先驗(yàn)(spatial prior)不僅可以被用于分割基線掃描(如最初進(jìn)行的那樣)而且還可以被用于分割縱向掃描。因此,此擴(kuò)展技術(shù)允許對(duì)不同時(shí)間點(diǎn)(例如第0天、第3天、第15天等)進(jìn)行同時(shí)分割,這又允許測(cè)量時(shí)間點(diǎn)之間的差異。因此,在相繼的時(shí)間點(diǎn)從基線圖像獲取的受試者圖像可以被同時(shí)分割并被用于識(shí)別受試者自身變化(即在單個(gè)受試者內(nèi)在不同時(shí)間點(diǎn)的解剖學(xué)結(jié)構(gòu)的差異)。圖形構(gòu)建和流形嵌入為確定中間圖冊(cè)傳播步驟,將所有圖像嵌入到通過坐標(biāo)系統(tǒng)表示的流形中,所述坐標(biāo)系統(tǒng)通過應(yīng)用譜分析步驟獲得。譜分析技術(shù)具有的優(yōu)點(diǎn)是基于對(duì)數(shù)據(jù)項(xiàng)如圖像之間的成對(duì)的相似性度量產(chǎn)生特征坐標(biāo)。這與需要數(shù)據(jù)項(xiàng)之間的距離度量的方法如多維量表(MDS)形成對(duì)比。在譜分析步驟后,學(xué)習(xí)得到的坐標(biāo)系統(tǒng)中的兩個(gè)圖像之間的距離不僅依賴于它們之間最初的成對(duì)相似性,而且還依賴于各個(gè)圖像與總體中的剩余部分所具有的所有成對(duì)相似性。這使得坐標(biāo)系統(tǒng)嵌入中的距離成為比可能對(duì)噪音敏感的單個(gè)成對(duì)的相似性度量更穩(wěn)健的接近性度量??梢栽赩ox Luxburg(2007)中發(fā)現(xiàn)對(duì)譜分析方法的詳細(xì)介紹并且在Chimg(1997)中可以獲得進(jìn)一步的細(xì)節(jié)。譜分析步驟應(yīng)用于完全的、加權(quán)的和非直接的圖形G = (V,E),數(shù)據(jù)集中的各個(gè)圖像由一個(gè)頂點(diǎn)Vi表示。兩個(gè)頂點(diǎn)Vi和Vj之間的非負(fù)權(quán)Wij被相應(yīng)圖像之間的相似性Sij限定。在本發(fā)明中,使用基于強(qiáng)度的相似性。通過收集每個(gè)圖像對(duì)的邊權(quán)Wij = Sij來獲得G
的權(quán)矩陣W并且對(duì)角矩陣T含有各個(gè)頂點(diǎn)的度數(shù)和《=Σ^; °
J來源于譜分析步驟的特征數(shù)據(jù)的維度可以由使用者選擇。在我們的工作中,我們依次測(cè)試了特征數(shù)據(jù)的各個(gè)維度并且評(píng)估了在四個(gè)受試者組(年輕人、AD、MCI和老年人對(duì)照受試者)之間的辨別能力。基于每個(gè)聚類(cluster)的質(zhì)心使用平均聚類間距離測(cè)量各個(gè)特征維度的辨別能力。對(duì)于被研究的群組,當(dāng)使用二維特征時(shí)辨別能力最大并且在之后降低(參見圖2、。我們因此選擇使用2D譜特征作為其中嵌入數(shù)據(jù)的坐標(biāo)空間。圖像相似性在我們的方法的優(yōu)選實(shí)施方案中,我們使用圖像對(duì)Ii和Ij之間基于強(qiáng)度的相似性。此相似性基于歸一化的交互信息(NMI) (Studholme等,1999),其與圖像I的熵H(I)以及兩個(gè)圖像的相關(guān)熵Hdi ;Ij) 一起被定義為
權(quán)利要求
1.一種處理醫(yī)用圖像的方法,所述方法由計(jì)算機(jī)處理器執(zhí)行并且包括以下步驟(a)獲得包含一個(gè)或多個(gè)圖像的一個(gè)或多個(gè)圖冊(cè),在所述一個(gè)或多個(gè)圖像中已經(jīng)用標(biāo)記數(shù)據(jù)標(biāo)記了一個(gè)或多個(gè)解剖學(xué)特征;(b)獲得多個(gè)未標(biāo)記的圖像;(c)將所述標(biāo)記的和未標(biāo)記的圖像進(jìn)行比較并且選擇與一個(gè)或多個(gè)所述標(biāo)記的圖像最相似的一個(gè)或多個(gè)未標(biāo)記的圖像;(d)使標(biāo)記數(shù)據(jù)從一個(gè)或多個(gè)最接近的標(biāo)記圖像傳播到那些選定的一個(gè)或多個(gè)圖像中的每個(gè)圖像,由此標(biāo)記所述選定的一個(gè)或多個(gè)圖像中的每個(gè)圖像的對(duì)應(yīng)的一個(gè)或多個(gè)解剖學(xué)特征并且使所述選定的一個(gè)或多個(gè)圖像成為標(biāo)記的一個(gè)或多個(gè)圖像;以及(e)從步驟(c)開始反復(fù)地重復(fù),由此標(biāo)記所述未標(biāo)記的圖像中的其他圖像。
2.權(quán)利要求1所述的方法,其中所述將所述標(biāo)記的和未標(biāo)記的圖像進(jìn)行比較的步驟包括將所述圖像嵌入到低維坐標(biāo)系統(tǒng)中。
3.權(quán)利要求1或權(quán)利要求2所述的方法,其中所述低維坐標(biāo)系統(tǒng)是二維坐標(biāo)空間。
4.任一前述權(quán)利要求所述的方法,其中所述將所述標(biāo)記的和未標(biāo)記的圖像進(jìn)行比較的步驟包括通過比較所述圖像的組中每對(duì)圖像的一個(gè)或多個(gè)各自的解剖學(xué)特征來限定一組成對(duì)的相似性度量。
5.權(quán)利要求4所述的方法,其中所述將所述標(biāo)記的和未標(biāo)記的圖像進(jìn)行比較的步驟還包括對(duì)所述成對(duì)的相似性度量進(jìn)行譜分析操作。
6.權(quán)利要求4或權(quán)利要求5所述的方法,其中所述成對(duì)的相似性度量表示一對(duì)圖像之間的強(qiáng)度相似性。
7.權(quán)利要求4、權(quán)利要求5或權(quán)利要求6所述的方法,其中所述成對(duì)的相似性度量表示一對(duì)圖像之間形變的量。
8.任一前述權(quán)利要求所述的方法,其中所述傳播標(biāo)記數(shù)據(jù)的步驟包括基于分類器融合技術(shù)使標(biāo)記數(shù)據(jù)從多個(gè)所述最接近的標(biāo)記的圖像傳播。
9.任一前述權(quán)利要求所述的方法,所述方法在步驟(d)后及步驟(e)前還包括對(duì)新近傳播的標(biāo)記數(shù)據(jù)進(jìn)行基于強(qiáng)度的精修操作的步驟。
10.任一前述權(quán)利要求所述的方法,其中所述圖像是不同受試者的圖像。
11.權(quán)利要求1至9中任一項(xiàng)所述的方法,其中至少一些所述圖像是相同受試者的但在不同的時(shí)間點(diǎn)獲取的圖像。
12.任一前述權(quán)利要求所述的方法,其中所述圖像是磁共振圖像。
13.任一前述權(quán)利要求所述的方法,所述方法還包括對(duì)表示病況出現(xiàn)與否的解剖學(xué)特征進(jìn)行標(biāo)記并且使用所述特征得到所述病況的生物標(biāo)志。
14.一種成像設(shè)備,所述成像設(shè)備被配置成執(zhí)行任一前述權(quán)利要求所述的方法。
15.權(quán)利要求14所述的成像設(shè)備,所述成像設(shè)備是醫(yī)用掃描儀。
16.權(quán)利要求15所述的成像設(shè)備,所述成像設(shè)備是MRI掃描儀。
17.一種圖像處理設(shè)備,所述圖像處理設(shè)備被配置成執(zhí)行權(quán)利要求1至13中任一項(xiàng)所述的方法。
18.一種計(jì)算機(jī)系統(tǒng),所述計(jì)算機(jī)系統(tǒng)被配置成執(zhí)行權(quán)利要求1至13中任一項(xiàng)所述的方法。
19.一種計(jì)算機(jī)程序,所述計(jì)算機(jī)程序包括當(dāng)在計(jì)算機(jī)處理器上運(yùn)行時(shí)用于執(zhí)行權(quán)利要求1至13中任一項(xiàng)所述的方法的編碼的指令。
20.編碼權(quán)利要求19所述的計(jì)算機(jī)程序的計(jì)算機(jī)可讀介質(zhì)或物理載波信號(hào)。
21.權(quán)利要求13所述的方法,所述方法還包括在所述生物標(biāo)志的基礎(chǔ)上給受試者分配以診斷分類。
22.權(quán)利要求13所述的方法,所述方法還包括在所述生物標(biāo)志的基礎(chǔ)上量化受試者對(duì)治療的反應(yīng)。
23.權(quán)利要求13所述的方法,所述方法還包括在所述生物標(biāo)志的基礎(chǔ)上選擇受試者的治療。
24.一種處理醫(yī)用圖像的方法,所述方法基本上如在本文中參照附圖的任一組合所描述的以及如在附圖的任一組合中所圖示的那樣。
25.一種成像設(shè)備、圖像處理設(shè)備或計(jì)算機(jī)系統(tǒng),所述成像設(shè)備、圖像處理設(shè)備或計(jì)算機(jī)系統(tǒng)基本上如在本文中參照附圖的任一組合所描述的那樣以及如在附圖的任一組合中所圖示的那樣。
26.一種計(jì)算機(jī)程序,所述計(jì)算機(jī)程序基本上如在本文中參照附圖的任一組合所描述的那樣以及如在附圖的任一組合中所圖示的那樣。
全文摘要
一種用于處理醫(yī)用圖像的方法,所述方法由計(jì)算機(jī)處理器執(zhí)行并且包括以下步驟(a)獲得包含一個(gè)或多個(gè)圖像的一個(gè)或多個(gè)圖冊(cè),在所述一個(gè)或多個(gè)圖像中已經(jīng)用標(biāo)記數(shù)據(jù)標(biāo)記了一個(gè)或多個(gè)解剖學(xué)特征;(b)獲得多個(gè)未標(biāo)記的圖像;(c)將標(biāo)記的和未標(biāo)記的圖像進(jìn)行比較并且選擇與一個(gè)或多個(gè)標(biāo)記的圖像最相似的一個(gè)或多個(gè)未標(biāo)記的圖像;(d)使標(biāo)記數(shù)據(jù)從一個(gè)或多個(gè)最接近的標(biāo)記圖像傳播到那些選定的一個(gè)或多個(gè)圖像中的每個(gè)圖像,由此標(biāo)記選定的一個(gè)或多個(gè)圖像中的每個(gè)圖像的對(duì)應(yīng)的一個(gè)或多個(gè)解剖學(xué)特征并且使選定的一個(gè)或多個(gè)圖像成為標(biāo)記的一個(gè)或多個(gè)圖像;以及(e)從步驟(c)開始反復(fù)地重復(fù),由此標(biāo)記未標(biāo)記的圖像中的其他圖像。
文檔編號(hào)G06K9/46GK102598054SQ201080043906
公開日2012年7月18日 申請(qǐng)日期2010年9月30日 優(yōu)先權(quán)日2009年9月30日
發(fā)明者伊恩·保羅·丹尼爾·呂凱, 保羅·阿爾雅巴, 羅賓·沃爾茨 申請(qǐng)人:帝國(guó)創(chuàng)新有限公司