本發(fā)明涉及一種spect成像方式,特別涉及一種在不改變錐形束spect檢測器準(zhǔn)直器孔徑大小的前提下提高成像分辨率的方法,屬于醫(yī)學(xué)圖像重建領(lǐng)域。
背景技術(shù):
單光子發(fā)射計算機斷層成像(singlephotonemissioncomputedtomography,spect)是核醫(yī)學(xué)的一項重要技術(shù),目前在臨床上有廣泛的應(yīng)用。在單光子發(fā)射斷層成像中,放射性示蹤劑被注入到病人體內(nèi),根據(jù)示蹤劑在人體內(nèi)分布,重建其在人體內(nèi)的密度分布圖像,該圖像能夠直觀反映人體器官新陳代謝水平的高低,從而實現(xiàn)功能成像的目的,對于疾病的早期診斷具有重要的意義。
在spect中,放射性示蹤劑通過輻射γ光子,穿過人體后被檢測器檢測,最終重建得到斷層圖像。在spect中,檢測器的準(zhǔn)直器孔徑是限制圖像空間分辨率和靈敏度的關(guān)鍵因素。準(zhǔn)直器孔徑越小,圖像的空間分辨率越高,但接收到的γ光子數(shù)量越少,靈敏度越低,成像質(zhì)量會變差。準(zhǔn)直器孔徑越大,接收到的γ光子數(shù)量越多,靈敏度越高,但圖像的空間分辨率越低,成像質(zhì)量也會變差。
超分辨率重建是從多幅低分辨率圖像重建出高分辨率圖像。作為一種有效的提高圖像分辨率的方法,它已經(jīng)被廣泛應(yīng)用于汽車牌照、遙感衛(wèi)星,醫(yī)學(xué)成像等領(lǐng)域。近年來超分辨率重建被應(yīng)用到磁共振、ct、pet、spect成像領(lǐng)域,以提高圖像分辨率,但這些方法主要應(yīng)用在圖像域或平行及扇形檢測器投影重建領(lǐng)域。
錐形投影spect系統(tǒng)與傳統(tǒng)的平行束及扇形束spect相比,具有高分辨率的特點,有利于對小目標(biāo)進行檢測及重建,錐形束投影重建和超分辨率重建相結(jié)合,將會極大提高spect的成像分辨率。
技術(shù)實現(xiàn)要素:
為了解決單光子發(fā)射計算機斷層成像(singlephotonemissioncomputedtomography,spect)分辨率低的問題,本發(fā)明公開的一種提高錐形束spect成像分辨率的方法要解決的技術(shù)問題是提供一種提高成像分辨率的方法,在不改變檢測器中準(zhǔn)直器孔徑大小的前提下,提高spect圖像的分辨率。
本發(fā)明是通過下述技術(shù)方案實現(xiàn)的。
本發(fā)明公開的一種提高錐形束spect成像分辨率的方法,包括如下步驟:在不改變檢測器中準(zhǔn)直器孔徑大小的前提下進行低分辨率投影圖采集;由于多幅低分辨率圖像都對每個角度同一目標(biāo)進行成像,但相互之間又包含不同的投影數(shù)據(jù)信息,通過對采集到的多幅低分辨率投影圖像的數(shù)據(jù)信息進行處理,得到每個角度下的高分辨率投影數(shù)據(jù);由此能夠從多幅低分辨率投影圖像中重建出高分辨率spect圖像,即提高重建spect圖像的分辨率。
所述的重建出高分辨率spect圖像方法可選用解析法或迭代法。
本發(fā)明公開的一種提高錐形束spect成像分辨率的方法,包括如下步驟:
步驟一、低分辨率投影圖采集。在不改變檢測器中準(zhǔn)直器孔徑大小的前提下,檢測器繞物體進行旋轉(zhuǎn),每次旋轉(zhuǎn)
步驟一中單個旋轉(zhuǎn)角度測量每組m個投影數(shù)據(jù)的具體采集方法通過方法1或方法2實現(xiàn)。
方法1:檢測器在其所在的平面內(nèi)繞檢測器中心旋轉(zhuǎn),共需要旋轉(zhuǎn)m次,每次旋轉(zhuǎn)后采集單個投影數(shù)據(jù),并記錄每次旋轉(zhuǎn)的角度。
所述的檢測器在其所在的平面內(nèi)繞檢測器中心旋轉(zhuǎn)的旋轉(zhuǎn)角度為固定角度均勻旋轉(zhuǎn)或不固定角度非均勻旋轉(zhuǎn),優(yōu)選在設(shè)定的角度θ范圍內(nèi)均勻旋轉(zhuǎn)m個角度,即每次旋轉(zhuǎn)
方法2:檢測器中心每次需平移預(yù)設(shè)的距離,共需要平移m次,每次平移后采集單個投影數(shù)據(jù),并記錄每次平移的距離。
每次平移的距離優(yōu)選檢測器中心每次平移高分辨率圖像一個像素大小的距離。
步驟二、對步驟一中檢測器繞物體進行旋轉(zhuǎn),每次旋轉(zhuǎn)
步驟三、根據(jù)步驟二得到的高分辨率投影數(shù)據(jù)重建出高分辨率spect圖像。
由于步驟一中多幅低分辨率投影圖像都對同一目標(biāo)進行成像,但相互之間又包含不同的投影數(shù)據(jù)信息,由此能夠得到高分辨率投影圖像,相應(yīng)的重建圖像的分辨率也會相應(yīng)的提高。由此能夠從多幅低分辨率投影圖像中重建出高分辨率spect圖像,因此能夠在不改變檢測器中準(zhǔn)直器孔徑大小的前提下,提高spect圖像的分辨率。
步驟三中重建高分辨率spect圖像方法選用解析法或迭代法。
步驟一中所述的檢測器繞物體旋轉(zhuǎn)的次數(shù)n,以及在每個角度下分別在其所在位置平面測得的投影數(shù)據(jù)個數(shù)m與實際低分辨率圖像大小和需獲得的高分辨率圖像大小有關(guān)。若低分辨率檢測器分辨率為ls×ls個像素,目標(biāo)獲得的高分辨率像素為hs×hs個像素,則m、n、hs、ls可由公式(1)、(2)、(3)給出。
n≥hs(1)
m≥de×de(3)
步驟二中所述的對每個角度下采集到的m幅低分辨率圖像進行處理,根據(jù)步驟一中檢測器繞物體進行旋轉(zhuǎn),單個旋轉(zhuǎn)角度測量每組m個投影數(shù)據(jù)的兩種采集方法,有針對兩種采集方法的不同實現(xiàn)步驟,即分別針對方法1、方法2采集到的m幅低分辨率圖像進行處理的方法。
針對方法1中對檢測器繞物體旋轉(zhuǎn)后采集到的m幅低分辨率圖像進行處理的方法包括如下步驟:
記
步驟2.1:以檢測器所在的平面建立坐標(biāo)系,坐標(biāo)系的中心為檢測器平面中心。對phi按建立的坐標(biāo)系進行旋轉(zhuǎn)操作,共需要旋轉(zhuǎn)m次。設(shè)旋轉(zhuǎn)第i個角度后的高分辨率投影為ph(i+1),則:
其中(x,y)為phi中像素點的坐標(biāo),(x',y')為ph(i+1)中像素點的坐標(biāo),γ為ph(i+1)相對于phi旋轉(zhuǎn)的角度。根據(jù)公式(4)算得的ph(i+1)中像素點的坐標(biāo)(x',y'),由雙線性插值公式求得ph(i+1)中每個點的像素值。
步驟2.2:將每個低分辨率區(qū)域?qū)?yīng)的高分辨率區(qū)域塊的所有像素點按照公式(5)疊加,即對ph(i+1)進行降采樣得到pl(i+1)。
步驟2.3:將經(jīng)步驟2.2降采樣得到的低分辨率圖像與直接采集到的低分辨率圖像按照公式(6)進行比較,求得真實低分辨率投影與目前低分辨率率投影的差值δpl(i+1)。
步驟2.4:將經(jīng)步驟2.3得到的對應(yīng)像素點的差值乘以每次調(diào)節(jié)的步長作為調(diào)節(jié)權(quán)重,得到對高分辨率圖像的調(diào)節(jié)權(quán)重如公式(7)所示。
x'=x×de+i(i=0,1,2…(de-1))
y'=y(tǒng)×de+i(i=0,1,2…(de-1))
求得ph(i+1)中每個點的調(diào)節(jié)權(quán)重δph(i+1),其中step為每次的調(diào)節(jié)步長,調(diào)節(jié)步長一般在0到1之間取值。
步驟2.5:由步驟2.4得到的調(diào)節(jié)權(quán)重δph(i+1)按照公式(8)更新ph(i+1)。
ph(i+1)=ph(i+1)+δph(i+1)(8)
步驟2.6;回到步驟2.1進行下一次旋轉(zhuǎn)角度調(diào)節(jié)或進行下一次循環(huán),直到達到要求的迭代結(jié)束條件,得到滿足投影分辨率要求的高分辨率投影圖像,結(jié)束迭代過程,即得到檢測器繞物體旋轉(zhuǎn)時單個角度下的高分辨率投影數(shù)據(jù)。
所述的迭代結(jié)束條件為:計算前后兩次迭代過程高分辨率圖像的差值的二范數(shù),并與結(jié)束迭代閾值條件進行比較,從而確定是否結(jié)束迭代。所述的迭代閾值條件根據(jù)需要獲得的高分辨率投影精度而定。
針對方法2中對檢測器繞物體旋轉(zhuǎn)后采集到的m幅低分辨率圖像進行處理的方法包括如下步驟:
記
步驟2.7:以檢測器所在的平面建立坐標(biāo)系,坐標(biāo)系的中心為檢測器平面中心。對phi按建立的坐標(biāo)系進行平移操作,共需要平移m次。設(shè)平移第i次后的高分辨率投影為ph(i+1),則:
(x',y')=(x,y)+(δx,δy)(9)
其中(δx,δy)為ph(i+1)中心相對于phi中心的位移。(x',y')和(x,y)分別為ph(i+1)和phi上的點。由公式(10)求得ph(i+1)中每個點的像素值。
ph(i+1)(x',y')=phi(x0,y0)(10)
步驟2.8:按步驟2.2-步驟2.5對高分辨率圖像進行調(diào)節(jié)。
步驟2.9:回到步驟2.7進行下一次平移調(diào)節(jié)或進行下一次循環(huán),直到達到要求的迭代結(jié)束條件,得到滿足投影分辨率要求的高分辨率投影圖像,結(jié)束迭代過程,即得到檢測器繞物體旋轉(zhuǎn)時單個角度下的高分辨率投影數(shù)據(jù)。
所述的迭代結(jié)束條件同步驟2.6里面的迭代結(jié)束條件,即所述的迭代結(jié)束條件為:計算前后兩次迭代過程高分辨率圖像的差值的二范數(shù),并與結(jié)束迭代閾值條件進行比較,從而確定是否結(jié)束迭代。所述的迭代閾值條件根據(jù)需要獲得的高分辨率投影精度而定。
有益效果:
本發(fā)明公開的一種提高錐形束spect成像分辨率的方法,由于多幅低分辨率投影圖像都對同一目標(biāo)進行成像,但相互之間又包含不同的信息,由此能夠根據(jù)步驟二的處理方法從多幅低分辨率投影圖像來獲取一幅高分辨率投影圖像,相應(yīng)的重建圖像的分辨率也會相應(yīng)的提高。因此能夠在不改變檢測器中準(zhǔn)直器孔徑大小的前提下,提高spect圖像的分辨率。
附圖說明
圖1為錐形spect檢測器模型示意圖。
圖2為錐形檢測器繞物體旋轉(zhuǎn)采集示意圖。
圖3為檢測器所在平面坐標(biāo)系示意圖,坐標(biāo)系的中心為檢測器中心。
圖4為方法一中檢測器平面繞其自身平面中心旋轉(zhuǎn)示意圖。
圖5為方法二中檢測器中心每次移動位置示意圖。
圖6為理想高分辨率投影圖像與實際低分辨率投影圖像對應(yīng)關(guān)系示意圖。
圖7為檢測器繞物體旋轉(zhuǎn)時單個角度下由多幅低分辨率投影圖像獲得高分辨率投影圖像流程圖。
圖8為由低分辨率投影圖像獲得高分辨率投影圖像時單次迭代調(diào)整過程示意圖。
圖9為二維線性插值坐標(biāo)示意圖。
圖10為不同層原始圖像、真實模擬高分辨率投影圖fdk重建圖像、由低分辨率投影旋轉(zhuǎn)移動獲得高分辨率投影經(jīng)fdk重建圖像、由低分辨率投影平行移動獲得高分辨率投影經(jīng)fdk重建圖像以及低分辨率fdk重建圖像。其中檢測器中心位于物體44層,檢測器焦點到檢測器中心的距離為380個像素,焦點到物體中心的距離為310個像素。
具體實施方式
下面結(jié)合附圖和具體實施例對本發(fā)明方法的實施方式做詳細說明。
實施例1:
檢測器繞物體旋轉(zhuǎn)的同時,其自身也需在其所在平面進行移動。本實例以低分辨率檢測器分辨率為32×32個像素,目標(biāo)獲得的高分辨率像素為128×128,n=128,m=16,de=4,為例進行說明,但是并不說明該發(fā)明僅限于這些條件,本發(fā)明中所涉及的范圍僅由權(quán)利說明書限定。
本實施例公開的一種針對錐形束spect超分辨率成像方法,其具體實施步驟如下:
步驟一、低分辨率投影圖采集。在不改變檢測器中準(zhǔn)直器孔徑大小的前提下,檢測器繞物體進行旋轉(zhuǎn),每次旋轉(zhuǎn)
在每個角度下,分別在檢測器所在位置平面測得一組投影數(shù)據(jù),每組投影數(shù)據(jù)包含16個投影數(shù)據(jù),單個旋轉(zhuǎn)角度測量每組16個投影數(shù)據(jù)結(jié)束后,檢測器繞物體繼續(xù)旋轉(zhuǎn)
方法1:檢測器在其所在的平面內(nèi)繞檢測器中心旋轉(zhuǎn),在
方法2:檢測器中心每次平移高分辨率圖像一個像素大小的距離,共需要平移16次,每次平移后采集單個投影數(shù)據(jù),并記錄每次平移的距離。檢測器中心移動位置如圖5所示。圖5中的小黑點為檢測器每次移動后的中心位置。
步驟二、對步驟一中檢測器繞物體進行旋轉(zhuǎn),每次旋轉(zhuǎn)
步驟二中所述的每個角度下采集到的16幅低分辨率投影圖像進行處理,根據(jù)步驟一中檢測器繞物體進行旋轉(zhuǎn),單個旋轉(zhuǎn)角度測量每組16個投影數(shù)據(jù)的兩種采集方法,有針對兩種采集方法的不同實現(xiàn)步驟,即分別針對方法1、方法2采集到的16幅低分辨率圖像進行處理的方法。
針對方法1中對檢測器繞物體旋轉(zhuǎn)后采集到的16幅低分辨率投影圖像進行處理的方法包括如下步驟:
記
圖6為理想高分辨率投影圖像與實際低分辨率投影圖像之間關(guān)系示意圖。在實際采集過程中,通過檢測器平面繞其自身旋轉(zhuǎn),可以獲得圖中右邊每個角度下的低分辨率圖像。圖中標(biāo)粗黑色方框代表高分辨率像素點與低分辨率像素點之間的對應(yīng)關(guān)系。
圖7為由16幅32×32像素的低分辨率投影圖像獲得1幅128×128像素高分辨率投影圖像流程圖。
流程圖中初始高分辨率投影圖像獲得說明:
初始高分辨率圖像獲得選用下述兩種方法之一:(1),對高分辨率圖像隨機賦值,此種方法在最終迭代后會有大量的隨機噪聲。(2),由采集到的低分辨率圖像插值成高分辨率圖像大小,并將此圖像作為原始高分辨率圖像,此種方法在最終迭代后獲得高分辨率圖像較好,且收斂速度較快。
單次迭代說明:根據(jù)低分辨率圖像采集時每次旋轉(zhuǎn)的角度大小,高分辨率圖像繞其中心每次旋轉(zhuǎn)相同的角度,然后進行降采樣。此處降采樣的方法為:將每個低分辨率區(qū)域?qū)?yīng)的高分辨率區(qū)域塊的所有像素疊加(如圖8所示:將高分辨率圖像每個加粗黑色方框內(nèi)的所有像素值相加,作為對應(yīng)低分辨率圖像中的一個像素值)。將經(jīng)降采樣得到的低分辨率圖像與采集到的低分辨率圖像進行比較,將對應(yīng)像素點的差值乘以每次調(diào)節(jié)的步長作為調(diào)節(jié)權(quán)重,對高分辨率圖像進行調(diào)節(jié)。
對應(yīng)過程如圖8所示:(a)為高分辨率圖像,(b)為由高分辨率圖像經(jīng)降采樣得到的圖像,(c)為實際采集到的低分辨率圖像,(d)為由(b)和(c)相減的差值得到的圖像,將得到的(d)作為調(diào)節(jié)依據(jù),對原始高分辨率圖像(a)進行調(diào)節(jié),最終得到調(diào)節(jié)后的高分辨率圖像(e)。然后獲得的(e)進入下一次調(diào)節(jié)過程。
迭代結(jié)束條件:迭代結(jié)束條件的表示有兩種方式,這兩種方式實質(zhì)上是一樣的,(1),通過計算圖8中(a)與(e)的差值的二范數(shù)并與結(jié)束迭代閾值條件進行比較,從而確定是否結(jié)束迭代。(2),計算圖8中(d)矩陣的二范數(shù)并與結(jié)束迭代閾值條件進行比較,從而確定是否結(jié)束迭代。
步驟2.1:對phi按圖(3)所示建立坐標(biāo)系進行旋轉(zhuǎn)操作,共需要旋轉(zhuǎn)16次。設(shè)旋轉(zhuǎn)第i個角度后的高分辨率投影為ph(i+1),則根據(jù)公式(4)算得的ph(i+1)中像素點的坐標(biāo)(x',y'),由線性插值公式(11)求得ph(i+1)中每個點的像素值。圖9為二維線性插值坐標(biāo)示意圖。
步驟2.2:將每個低分辨率區(qū)域?qū)?yīng)的高分辨率區(qū)域塊的所有像素按照公式(5)疊加,即對ph(i+1)進行降采樣得到pl(i+1)。
步驟2.3:將經(jīng)步驟2.2降采樣得到的低分辨率圖像與直接采集到的低分辨率圖像按照公式(6)進行比較,求得真實低分辨率投影與目前低分辨率率投影的差值δpl(i+1)。
步驟2.4:將經(jīng)步驟2.3得到的對應(yīng)像素點的差值乘以每次調(diào)節(jié)的步長作為調(diào)節(jié)權(quán)重,得到對高分辨率圖像的調(diào)節(jié)權(quán)重如公式(7)所示。求得ph(i+1)中每個點的調(diào)節(jié)權(quán)重δph(i+1),其中step為每次的調(diào)節(jié)步長,step=1。
步驟2.5:由步驟2.4得到的調(diào)節(jié)權(quán)重δph(i+1)按照公式(8)更新ph(i+1)。
步驟2.6:回到步驟2.1進行下一次旋轉(zhuǎn)角度調(diào)節(jié)或進行下一次循環(huán),直到達到要求的迭代結(jié)束條件,得到滿足投影分辨率要求的高分辨率投影圖像,結(jié)束迭代過程,即得到檢測器繞物體旋轉(zhuǎn)時單個角度下的高分辨率投影數(shù)據(jù)。
所述的迭代結(jié)束條件為:計算前后兩次迭代過程高分辨率圖像的差值的二范數(shù),并與結(jié)束迭代閾值條件進行比較,從而確定是否結(jié)束迭代。所述的迭代閾值條件根據(jù)需要獲得的高分辨率投影精度而定。
針對方法2中對檢測器繞物體旋轉(zhuǎn)后采集到的16幅低分辨率投影圖像進行處理的方法包括如下步驟:
記
步驟2.7:對phi按圖(4)所示的坐標(biāo)系進行平移操作,共需要平移16次。設(shè)平移第i次后的高分辨率投影為ph(i+1),則根據(jù)公式(9),(10)求得ph(i+1)中每個點的像素值。
步驟2.8:按步驟2.2-步驟2.5對高分辨率圖像進行調(diào)節(jié)。
步驟2.9:回到步驟2.7進行下一次平移調(diào)節(jié)或進行下一次循環(huán),直到達到要求的迭代結(jié)束條件,得到滿足投影分辨率要求的高分辨率投影圖像,結(jié)束迭代過程,即得到檢測器繞物體旋轉(zhuǎn)時單個角度下的高分辨率投影數(shù)據(jù)。
所述的迭代結(jié)束條件同步驟2.6里面的迭代結(jié)束條件。
步驟三、根據(jù)步驟二得到的高分辨率投影數(shù)據(jù)根據(jù)fdk重建算法重建出高分辨率spect圖像。圖10為利用fdk重建算法對投影圖像重建所得結(jié)果。
以上結(jié)合具體實施例對本發(fā)明的技術(shù)方案和具體實施方式作了說明,但這些說明不能被理解為限制了本發(fā)明的范圍,這些僅是舉例說明,可以對這些實施方式做出多種變更或修改,而不背離本發(fā)明的原理和實質(zhì)。本發(fā)明的保護范圍由隨附的權(quán)利要求書限定,任何在本發(fā)明權(quán)利要求基礎(chǔ)上的改動都是本發(fā)明的保護范圍。