專利名稱：一種x射線源焦點(diǎn)漂移ct重建方法
技術(shù)領(lǐng)域：
本發(fā)明涉及一種X射線源焦點(diǎn)漂移CT重建方法，屬于X射線計(jì)算機(jī)層析成像(CT) 技術(shù)領(lǐng)域。
背景技術(shù)：
在X射線CT系統(tǒng)中，X射線源發(fā)出X射線，從不同角度穿過(guò)被檢測(cè)物體的某一區(qū)域， 放置于射線源對(duì)面的探測(cè)器在相應(yīng)角度接受，然后根據(jù)各角度射線不同程度的衰減，利 用一定的重建算法和計(jì)算機(jī)進(jìn)行運(yùn)算，重建出物體被掃描區(qū)域的射線線衰減系數(shù)分布映 射圖像，從而實(shí)現(xiàn)由投影重建圖像，無(wú)損地再現(xiàn)物體在該區(qū)域內(nèi)的介質(zhì)密度、成分和結(jié) 構(gòu)形態(tài)等特征。
現(xiàn)代CT系統(tǒng)一般采用濾波反投影(FBP)算法進(jìn)行圖像重建。FBP算法要求射線源焦 點(diǎn)、轉(zhuǎn)臺(tái)旋轉(zhuǎn)中心以及探測(cè)器的位置必須滿足以下四個(gè)條件(1)射線源焦點(diǎn)與旋轉(zhuǎn) 中心、探測(cè)器中心在一條直線上；(2)焦點(diǎn)與旋轉(zhuǎn)中心的連線必須垂直于探測(cè)器；(3) 在掃描過(guò)程中，焦點(diǎn)到旋轉(zhuǎn)中心的距離和旋轉(zhuǎn)中心到探測(cè)器中心的距離必須保持不變；
(4)在任意掃描角度下，必須同時(shí)滿足上述條件。在實(shí)際工業(yè)CT系統(tǒng)中，由于受到各 種因素的影響，系統(tǒng)往往不能同時(shí)滿足上述條件。如果采用經(jīng)典的FBP重建算法進(jìn)行重 建，容易產(chǎn)生變形和偽像，難以得到高質(zhì)量的CT圖像。
隨著研究的深入，在系統(tǒng)不能滿足上述條件時(shí)的圖像重建問(wèn)題得到了CT界的關(guān)注和 研究。陳明等，轉(zhuǎn)臺(tái)中心偏移與探測(cè)器偏轉(zhuǎn)造成的圖像偽影分析與校正，第十一屆中國(guó) 體視學(xué)與圖像分析學(xué)術(shù)會(huì)議論文集，2006: 489—492中，分析了在CT系統(tǒng)存在旋轉(zhuǎn)中心 偏移及探測(cè)器偏轉(zhuǎn)時(shí)采用傳統(tǒng)FBP算法重建圖像造成偽像的原因，給出了含有探測(cè)器偏 轉(zhuǎn)誤差的重建公式及數(shù)值驗(yàn)證結(jié)果；Ge Wang, Revisit to Derivative-Free Noncircular Fan-Beam Reconstruction, CT理論與應(yīng)用，2003， Vol. 12(2): 39-44中，分析了焦點(diǎn) 到旋轉(zhuǎn)中心距離的變化對(duì)圖像質(zhì)量的影響，并且給出了一種非圓軌跡下的扇束重建公 式。
綜上所述，盡管目前很多學(xué)者分析了探測(cè)器偏轉(zhuǎn)誤差、旋轉(zhuǎn)中心誤差、幾何距離變 化等導(dǎo)致系統(tǒng)不能滿足上述條件的因素，并推導(dǎo)了相應(yīng)的圖像重建公式，但尚未有研究 過(guò)X射線源焦點(diǎn)漂移這一導(dǎo)致系統(tǒng)不能滿足上述條件的因素，也未提出這種情況下的重建方法。
發(fā)明內(nèi)容
本發(fā)明的技術(shù)解決問(wèn)題是克服現(xiàn)有技術(shù)的不足，提供一種X射線源焦點(diǎn)漂移CT重 建方法，該方法實(shí)現(xiàn)了X射線源焦點(diǎn)漂移情況下的CT圖像重建，獲得了正確的CT圖像。
本發(fā)明的技術(shù)解決方案 一種X射線源焦點(diǎn)漂移CT圖像重建方法，包括CT投影數(shù)據(jù) 的加權(quán)、濾波和反投影重建三個(gè)步驟，CT投影數(shù)據(jù)的加權(quán)采用公式(1)進(jìn)行，CT投影 數(shù)據(jù)的濾波采用公式(2)和(3)進(jìn)行，CT圖像反投影重建采用公式(4) 、 (5)和(6) 進(jìn)行
Z)3-ZW(/ )(s-d⑨)
(1.
1
(Z)'+r.si ) ".r cos(- —爭(zhēng)d⑨ r sin(- — 6)
(1)
(2)
(3)
(4)
(5)
(6)
其中，A(S,^為加權(quán)后CT投影數(shù)據(jù)；^"釣為CT掃描形成的投影數(shù)據(jù)；^是線陣探
測(cè)器上每個(gè)探元的位置；^是旋轉(zhuǎn)角度；D為未發(fā)生射線源焦點(diǎn)偏移時(shí)，焦點(diǎn)到轉(zhuǎn)臺(tái)旋 轉(zhuǎn)中心距離；^")是射線源焦點(diǎn)偏移距離，是旋轉(zhuǎn)角度^的函數(shù)；"'(")是"(^)對(duì)"的 導(dǎo)數(shù)；/7'》,歷是濾波后CT投影數(shù)據(jù)；l川為濾波函數(shù)；h(>)是| P |的傅立葉變換；^"^) 為極坐標(biāo)形式的待重建的CT圖像；^為加權(quán)因子；"'為發(fā)生射線源焦點(diǎn)偏移時(shí)，焦點(diǎn)到
轉(zhuǎn)臺(tái)旋轉(zhuǎn)中心距離；f為發(fā)生射線源焦點(diǎn)偏移時(shí)，物體的實(shí)際旋轉(zhuǎn)角度；^是某旋轉(zhuǎn)角 度下，通過(guò)重建點(diǎn)的射線與探測(cè)器的交點(diǎn)。 下面進(jìn)一步描述上述公式的推導(dǎo)過(guò)程。
本發(fā)明的CT掃描結(jié)構(gòu)如圖1所示F和F'分別為射線源焦點(diǎn)的理想位置和實(shí)際位置；
d為焦點(diǎn)沿平行于探測(cè)器方向偏移的距離，是旋轉(zhuǎn)角度e的函數(shù)，隨e的不同而不同；D
和D'分別為F、 F'到旋轉(zhuǎn)中心O的距離；重建坐標(biāo)系為xoy; Ae是焦點(diǎn)偏離角度；F0、F'M垂直于探測(cè)器；a(r, 9 )為極坐標(biāo)形式的被重建圖像；A為任意一個(gè)重建點(diǎn)；重建點(diǎn) A到實(shí)際焦點(diǎn)F'的距離為L(zhǎng);過(guò)重建點(diǎn)A的射線與探測(cè)器相交于點(diǎn)s,，與F'M之間的夾角為 y0; s"為探測(cè)器上一個(gè)探元的位置，過(guò)s''的射線與F'M之間的夾角為〃'；旋轉(zhuǎn)中心0 到過(guò)探元s''的射線的垂直距離為Xr'';直線AC平行于F's''，旋轉(zhuǎn)中心O到AC的垂直距 離為Xr。
莊天戈，CT原理與算法，上海交通大學(xué)出版社，1992中，給出了平行束濾波反投影 的重建公式，如下
=去 ",0) /<力-力")化r" W
20
由圖l可得到以下幾何關(guān)系式:
sin(沖=,d(/ ) c。s(沖=-^_
廣綱2 >2+(^-,)2 sin(/") = ■', 『'■ . eos(/")=-
嗣、c
X八s".d".," ，必",;一"尋"一,V
Xr — Xr" = ￡. sinOO — y")=丄 sinOO) cos0") —丄.cos(>0) sin(;r") i = (D'+ r' sin(P'- 。) / cos(Ay5 + ;k0)
D2—
cos(A〃 + ;k0) = cos(A-). cos(;k0) - sin(A^). sin(沖=
(7)
(8)
(9)
^匚—(Z)'+r.sin(/ '—。)
h(s)為理想濾波函數(shù)- h(力-,乂i川) (12) 重置濾波函數(shù)變量-
(10)
(11)- Zr") = A( 利用變量代換-
Jl/ ieVD+(s"-dW) # (13)
d(/ ))2
根據(jù)雅可比行列式計(jì)算公式
。3-ZW(-)(5"-,)
d(/ ))2
s"-，、") (14)
l氺
W，
(1—-
，)
如2 )3 D-，)""-力/Z)
由(7) ~ (15)式可得 一)- f二'd，卿
(15)
1 Z)3—ZW ))
(1-
"'⑨
/入
CD'+r.sin( "-，(s廣柳
g. r cos(yg - + dQg)' r sin(;g - gj D + rsin(/ —0
(16) :)(17)
(18)
(19)
(20)
2^式(16)即為射線源焦點(diǎn)漂移情況下的CT菌像精確重建公式。
本發(fā)明與現(xiàn)有技術(shù)相比的優(yōu)點(diǎn)如下
(1) 現(xiàn)有方法不能實(shí)現(xiàn)X射線源焦點(diǎn)漂移情況下的CT圖像重建問(wèn)題，本發(fā)明解決 了這一問(wèn)題；
(2) 本發(fā)明重建過(guò)程僅包括加權(quán)、濾波和反投影重建三個(gè)步驟，簡(jiǎn)單、高效。
圖1為本發(fā)明考慮射線源焦點(diǎn)漂移的X射線CT圖像重建方法對(duì)應(yīng)的掃描幾何結(jié)構(gòu)
圖2為經(jīng)典的Sh印p-Logan仿真模型；圖3 (a)是射線源焦點(diǎn)漂移為^(") = 1()時(shí)，傳統(tǒng)FBP算法得到的CT重建圖像；
圖3 (b)是射線源焦點(diǎn)漂移為"(")=1()時(shí)，本發(fā)明重建方法得到的CT重建圖像；
圖3 (c)是射線源焦點(diǎn)漂移為"(")=1()時(shí)，傳統(tǒng)FBP算法得到的CT重建圖像與仿 真模型的差值圖像；
圖3 (d)是射線源焦點(diǎn)漂移為^灼=1()時(shí)，本發(fā)明重建方法得到的CT重建圖像與 仿真模型的差值圖像；
圖4")是射線源焦點(diǎn)漂移為^釣="/2 + 5()時(shí)，傳統(tǒng)FBP算法得到的CT重建圖像； 圖4 (b)是射線源焦點(diǎn)漂移為"(釣=^2 + 5()時(shí)，本發(fā)明重建方法得到的CT重建圖
像；
圖4(c)是射線源焦點(diǎn)漂移為^") = "/2 + 5()時(shí)，傳統(tǒng)FBP算法得到的CT重建圖像 與仿真模型的差值圖像；
圖4 (d)是射線源焦點(diǎn)漂移為"(釣="/2 + 5()時(shí)，本發(fā)明重建方法得到的CT重建圖 像與仿真模型的差值圖像；
圖5 (a)是射線源焦點(diǎn)漂移為"^^2^^sinW)+U時(shí)，傳統(tǒng)FBP算法得到的CT重 建圖像；
圖5 (b)是射線源焦點(diǎn)漂移為"(釣-2Mx(sin(^+0時(shí)，本發(fā)明重建方法得到的CT
圖5 (c)是射線源焦點(diǎn)漂移為^釣-2Mx(sin^)+i)時(shí)，傳統(tǒng)FBP算法得到的CT重 建圖像與仿真模型的差值圖像；
圖5 (d)是射線源焦點(diǎn)漂移為^釣-2W"sin^)+"時(shí)，本發(fā)明重建方法得到的CT 重建圖隼與仿真模型的差值圖像。
具體實(shí)施例方式
本發(fā)明針對(duì)射線源焦點(diǎn)發(fā)生漂移時(shí)的CT圖像重建問(wèn)題，根據(jù)發(fā)生X射線源焦點(diǎn)漂移 時(shí)的CT掃描幾何結(jié)構(gòu)，推導(dǎo)CT投影數(shù)據(jù)的加權(quán)、濾波和反投影重建公式，通過(guò)在傳統(tǒng)FBP 重建公式中引入X射線源焦點(diǎn)漂移距離參數(shù)，形成新的重建方法，實(shí)現(xiàn)正確CT重建。具 體實(shí)現(xiàn)公式參見(jiàn)前述的(1) - (6)。
為驗(yàn)證本發(fā)明，在計(jì)算機(jī)上進(jìn)行了仿真實(shí)驗(yàn)。仿真實(shí)驗(yàn)的具體步驟如下 (1)建立待重建的二維圖像模型。圖2顯示了本發(fā)明仿真實(shí)驗(yàn)采用的經(jīng)典 Sheep化og肌模型o(2) 仿真實(shí)驗(yàn)參數(shù)設(shè)置。D為1200， e在360度范圍內(nèi)以l度為增量取360個(gè)投影角 度，探測(cè)器通道個(gè)數(shù)s為256，重建圖像大小為256x256，射線源焦點(diǎn)漂移考慮三種情況: ) = 10、 af(/ ) = A/2 + 50、 i/(/ ) = 200x(sin(〃)+l)。
(3) 根據(jù)投影幾何關(guān)系，生成投影數(shù)據(jù)。
(4) 根據(jù)傳統(tǒng)FBP重建公式進(jìn)行投影數(shù)據(jù)加權(quán)、濾波和反投影計(jì)算，獲得CT圖像。
(5) 根據(jù)公式(16) (20)進(jìn)行投影數(shù)據(jù)加權(quán)、濾波和反投影計(jì)算，獲得CT圖像。
圖2是經(jīng)典Sh印p-Logan仿真模型；圖3 (a)是射線源焦點(diǎn)漂移為"々)=1G時(shí)，傳
統(tǒng)FBP算法得到的CT重建圖像；圖3 (b)是射線源焦點(diǎn)漂移為"(")=1()時(shí)，本發(fā)明重
建方法得到的CT重建圖像；圖3 (c)是射線源焦點(diǎn)漂移為"(")=1()時(shí)，傳統(tǒng)FBP算法
得到的CT重建圖像與仿真模型的差值圖像；圖3 (d)是射線源焦點(diǎn)漂移為^^) = 1()時(shí)， 本發(fā)明重建文壇法得到的CT重建圖像與仿真模型的差值圖像；圖4 (a)是射線源焦點(diǎn)
漂移為"(")=^/2 + 5()時(shí)，傳統(tǒng)FBP算法得到的CT重建圖像；圖4(b)是射線源焦點(diǎn)漂 移為^/ )-"/2 + 5()時(shí)，本發(fā)明重建方法得到的CT重建圖像；圖4(c)是射線源焦點(diǎn)漂 移為c/(")-"/2 + 50時(shí)，傳統(tǒng)FBP算法得到的CT重建圖像與仿真模型的差值圖像；圖4 (d)是射線源焦點(diǎn)漂移為^") = "/2 + 5()時(shí)，本發(fā)明重建方法得到的CT重建圖像與仿 真模型的差值圖像；圖5 (a)是射線源焦點(diǎn)漂移為"(々)-2^x(sin(々)+"時(shí)，傳統(tǒng)FBP 算法得到的CT重建圖像；圖5 (b)是射線源焦點(diǎn)漂移為"々)=2()()"sinW)+1hf，本發(fā)
明重建方法得到的CT重建圖像；圖5(c)是射線源焦點(diǎn)漂移為"(々)-2^x(sin(々)+0時(shí)， 傳統(tǒng)FBP算法得到的CT重建圖像與仿真模型的差值圖像；圖5 (d)是射線源焦點(diǎn)漂移 為fl^)-200x(sin(/ );H')時(shí)，本發(fā)明重建文壇法得到的CT重建圖像與仿真模型的差值圖 像。
由圖2、 3(a)、 4(a)、 5(a)可知，在發(fā)生射線源焦點(diǎn)漂移時(shí)，傳統(tǒng)FBP算法重建得到 的CT圖像，存在較大的變形和偽影。由圖2、 3(b) 、 4(b) 、 5(b)可知，在發(fā)生射線源焦 點(diǎn)漂移時(shí)，本發(fā)明重建方法能實(shí)現(xiàn)正確的CT重建，獲得正確的CT圖像。
權(quán)利要求
1、一種X射線源焦點(diǎn)漂移CT圖像重建方法，包括CT投影數(shù)據(jù)的加權(quán)、濾波和反投影重建三個(gè)步驟，其特征在于所述CT投影數(shù)據(jù)的加權(quán)采用公式(1)進(jìn)行，所述CT投影數(shù)據(jù)的濾波采用公式(2)和(3)進(jìn)行，所述CT圖像反投影重建采用公式(4)、(5)和(6)進(jìn)行，<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><msub> <mi>P</mi> <mi>e</mi></msub><mrow> <mo>(</mo> <mi>s</mi> <mo>,</mo> <mi>β</mi> <mo>)</mo></mrow><mo>=</mo><mi>p</mi><mrow> <mo>(</mo> <mi>s</mi> <mo>,</mo> <mi>β</mi> <mo>)</mo></mrow><mo>·</mo><mfrac> <mn>1</mn> <mrow><msup> <mi>D</mi> <mn>2</mn></msup><mo>+</mo><msup> <mrow><mo>(</mo><mi>d</mi><mrow> <mo>(</mo> <mi>β</mi> <mo>)</mo></mrow><mo>)</mo> </mrow> <mn>2</mn></msup> </mrow></mfrac><mo>·</mo><mfrac> <mrow><msup> <mi>D</mi> <mn>3</mn></msup><mo>-</mo><mi>Dd</mi><mrow> <mo>(</mo> <mi>β</mi> <mo>)</mo></mrow><mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>d</mi> <mrow><mo>(</mo><mi>β</mi><mo>)</mo> </mrow> <mo>)</mo></mrow> </mrow> <msqrt><msup> <mi>D</mi> <mn>2</mn></msup><mo>+</mo><msup> <mrow><mo>(</mo><mi>s</mi><mo>-</mo><mi>d</mi><mrow> <mo>(</mo> <mi>β</mi> <mo>)</mo></mrow><mo>)</mo> </mrow> <mn>2</mn></msup> </msqrt></mfrac><mo>·</mo><mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mfrac><mrow> <msup><mi>d</mi><mo>′</mo> </msup> <mrow><mo>(</mo><mi>β</mi><mo>)</mo> </mrow></mrow><mrow> <mi>D</mi> <mo>-</mo> <mi>d</mi> <mrow><mo>(</mo><mi>β</mi><mo>)</mo> </mrow> <mrow><mo>(</mo><mi>s</mi><mo>-</mo><mi>d</mi><mrow> <mo>(</mo> <mi>β</mi> <mo>)</mo></mrow><mo>)</mo> </mrow> <mo>/</mo> <mi>D</mi></mrow> </mfrac> <mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo></mrow> </mrow>]]></math></maths><maths id="math0002" num="0002" ><math><![CDATA[ <mrow><msub> <msup><mi>p</mi><mo>′</mo> </msup> <mi>e</mi></msub><mrow> <mo>(</mo> <mi>s</mi> <mo>,</mo> <mi>β</mi> <mo>)</mo></mrow><mo>=</mo><msub> <mi>p</mi> <mi>e</mi></msub><mrow> <mo>(</mo> <mi>s</mi> <mo>,</mo> <mi>β</mi> <mo>)</mo></mrow><mo>*</mo><mfrac> <mn>1</mn> <mn>2</mn></mfrac><mi>h</mi><mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo></mrow> </mrow>]]></math></maths>h(s)＝F-1(|ρ|)(3)<maths id="math0003" num="0003" ><math><![CDATA[ <mrow><mi>a</mi><mrow> <mo>(</mo> <mi>r</mi> <mo>,</mo> <mi>θ</mi> <mo>)</mo></mrow><mo>=</mo><munderover> <mo>∫</mo> <mn>0</mn> <mrow><mn>2</mn><mi>π</mi> </mrow></munderover><mfrac> <mn>1</mn> <msup><mi>U</mi><mn>2</mn> </msup></mfrac><msub> <msup><mi>p</mi><mo>′</mo> </msup> <mi>e</mi></msub><mrow> <mo>(</mo> <msub><mi>s</mi><mn>1</mn> </msub> <mo>,</mo> <mi>β</mi> <mo>)</mo></mrow><mi>dβ</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo></mrow> </mrow>]]></math></maths><maths id="math0004" num="0004" ><math><![CDATA[ <mrow><mi>U</mi><mo>=</mo><mfrac> <mrow><mo>(</mo><msup> <mi>D</mi> <mo>′</mo></msup><mo>+</mo><mi>r</mi><mo>·</mo><mi>sin</mi><mrow> <mo>(</mo> <msup><mi>β</mi><mo>′</mo> </msup> <mo>-</mo> <mi>θ</mi> <mo>)</mo></mrow><mo>)</mo> </mrow> <mrow><mi>D</mi><mo>-</mo><mfrac> <mrow><mi>d</mi><mrow> <mo>(</mo> <mi>β</mi> <mo>)</mo></mrow> </mrow> <mi>D</mi></mfrac><mrow> <mo>(</mo> <msub><mi>s</mi><mn>1</mn> </msub> <mo>-</mo> <mi>d</mi> <mrow><mo>(</mo><mi>β</mi><mo>)</mo> </mrow> <mo>)</mo></mrow> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo></mrow> </mrow>]]></math></maths><maths id="math0005" num="0005" ><math><![CDATA[ <mrow><msub> <mi>s</mi> <mn>1</mn></msub><mo>=</mo><mfrac> <mrow><mi>D</mi><mo>·</mo><mi>r</mi><mi>cos</mi><mrow> <mo>(</mo> <mi>β</mi> <mo>-</mo> <mi>θ</mi> <mo>)</mo></mrow><mo>+</mo><mi>d</mi><mrow> <mo>(</mo> <mi>β</mi> <mo>)</mo></mrow><mo>·</mo><mi>r</mi><mi>sin</mi><mrow> <mo>(</mo> <mi>β</mi> <mo>-</mo> <mi>θ</mi> <mo>)</mo></mrow> </mrow> <mrow><mi>D</mi><mo>+</mo><mi>r</mi><mi>sin</mi><mrow> <mo>(</mo> <mi>β</mi> <mo>-</mo> <mi>θ</mi> <mo>)</mo></mrow> </mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo></mrow> </mrow>]]></math></maths>其中，pe(s，β)為加權(quán)后CT投影數(shù)據(jù)；p(s，β)為CT掃描形成的投影數(shù)據(jù)；s是線陣探測(cè)器上每個(gè)探元的位置；β是旋轉(zhuǎn)角度D為未發(fā)生射線源焦點(diǎn)偏移時(shí)，焦點(diǎn)到轉(zhuǎn)臺(tái)旋轉(zhuǎn)中心距離；d(β)是射線源焦點(diǎn)偏移距離，是旋轉(zhuǎn)角度β的函數(shù)；d′(β)是d(β)對(duì)β的導(dǎo)數(shù)；p′e(s，β)是濾波后CT投影數(shù)據(jù)；|ρ|為濾波函數(shù)；h(s)是|ρ|的傅立葉變換；a(r，θ)為極坐標(biāo)形式的待重建的CT圖像；U為加權(quán)因子；D′為發(fā)生射線源焦點(diǎn)偏移時(shí)，焦點(diǎn)到轉(zhuǎn)臺(tái)旋轉(zhuǎn)中心距離；β′為發(fā)生射線源焦點(diǎn)偏移時(shí)，物體的實(shí)際旋轉(zhuǎn)角度；s1是某旋轉(zhuǎn)角度下，通過(guò)重建點(diǎn)的射線與探測(cè)器的交點(diǎn)。
全文摘要
本發(fā)明屬于X射線計(jì)算機(jī)層析成像(CT)技術(shù)領(lǐng)域，具體為一種X射線源焦點(diǎn)漂移CT重建方法。該方法基于濾波反投影重建原理，根據(jù)X射線源焦點(diǎn)漂移時(shí)的CT掃描幾何結(jié)構(gòu)，推導(dǎo)CT投影數(shù)據(jù)的加權(quán)、濾波和反投影計(jì)算公式，在重建中引入X射線源焦點(diǎn)漂移距離參數(shù)，從而形成一種新的重建算法。傳統(tǒng)CT重建算法對(duì)應(yīng)的CT掃描幾何結(jié)構(gòu)不同于發(fā)生X射線源焦點(diǎn)漂移時(shí)的結(jié)構(gòu)，因此，不能應(yīng)用傳統(tǒng)CT重建算法對(duì)發(fā)生X射線源焦點(diǎn)漂移的二維CT投影進(jìn)行圖像重建。本發(fā)明解決了射線源焦點(diǎn)漂移情況下的X射線CT重建問(wèn)題，重建過(guò)程簡(jiǎn)單、高效。
文檔編號(hào)G06T11/00GK101620736SQ20091009028
公開(kāi)日2010年1月6日 申請(qǐng)日期2009年8月4日 優(yōu)先權(quán)日2009年8月4日
發(fā)明者健 傅, 斌 李, 江柏紅 申請(qǐng)人:北京航空航天大學(xué)