TWI517093B - Computer tomography reconstruction method - Google Patents
Computer tomography reconstruction method Download PDFInfo
- Publication number
- TWI517093B TWI517093B TW102136721A TW102136721A TWI517093B TW I517093 B TWI517093 B TW I517093B TW 102136721 A TW102136721 A TW 102136721A TW 102136721 A TW102136721 A TW 102136721A TW I517093 B TWI517093 B TW I517093B
- Authority
- TW
- Taiwan
- Prior art keywords
- image
- projection
- data
- dimensional
- difference value
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
Description
本發明係關於一種電腦斷層影像重建方法,特別是指一種透過二維系統矩陣進行投影、比對、反投影以及修正之方式,以降低運算量及儲存空間的電腦斷層影像重建方法。
目前電腦斷層重建影像中使用系統矩陣之重建法,雖已普遍廣泛地被使用在疊代式重建法;然而,隨著電腦科技技術以及影像擷取技術之提升,影像偵檢器之解析度越來越高,伴隨之資料量也越來越龐大,需要消耗相當大的運算及儲存資源。雖然一般儀器公司之開發採用多部影像工作站串聯以達成其運算目的,但是在技術上仍需儲存多組影像之數據,極佔儲存空間與運算之記憶體容量。
請參考美國公開專利第US6850585B2號和第US20120098832號報告,其係揭露一種使用多角度二維平面造影以呈現三維影像的影像重建方法,係針對偵檢系統掃描一待測物所產生之至少一條射線穿越待測物後所量測之一量測值進行影像重建,該影像重建方法包括有下列步驟:首先以複數個分別具有數值之立體像素對該待測物建立一影像空間;然後,對沿該射線方向進行投影並與該量測計數值比較以得到一校正計數值;接著,據該校正計數值與該反應線上每一個立體像素值,進行反投影;最後,根據該校正計數值與每一個立體像素值依其對應關係進行
修正以得到修正後的數值,並以更新後的值取代原數值,並視使用者需求重複疊代方式數次進行影像重建。而上述所知的技術中,疊代前需根據每一個影像的畫素進行三維影像的處理,其運算量較高。
本發明的目的,係在於提供一種電腦斷層影像重建方法,係利用不同角度的二維影像以及虛擬的二維系統矩陣,以圓球體空間為基礎,進行投影、比對、反投影以及修正等步驟進行處理,進而在預定差異值內獲得所欲獲得的斷層影像輸出,藉此以降低所需運算量與儲存空間。
為達上述目的,本發明係提供一種電腦斷層影像重建方法,其步驟係包括:收集資訊步驟:以一影像系統在每隔一預定角度之一位置環繞一物體,以該物體為中心收集該物體之不同方向之相對應的一投影資料,並將該投影資料記錄成該影像系統的一幾何參數,各方向的該幾何參數係建立成一組三維影像重建空間資訊;投影步驟:建構該影像系統之虛擬數據模型並對其立體特性簡化成一簡化數據模型,並虛擬化一組感興趣區域之一數據資料當作一預測資料,最後將該感興趣空間之數據資料透過該簡化數據模型進行模擬投影;比對步驟:將該模擬投影所得到的資訊與該影像系統所收集各角度之投影影像資訊進行差異比對,計算一差異值;反投影步驟:該差異值再透過該簡化數據模型進行反投影之運算;修正步驟:根據該反投影後之差異值修正原始三維影像重建空間的資訊;
循環步驟:確認該差異值是否小於一預定差異值,若否的話,則回到該投影步驟,若是的話,則下一組三維影像重建空間資訊開始依序進行該投影步驟、該比對步驟、該反投影步驟以及該修正步驟,直至各方向的該組三維影像重建空間資訊依序完成該投影步驟、該比對步驟、該反投影步驟以及該修正步驟,且該組三維影像重建空間的數據與前一次疊代之差異值小於該預定差異值
影像輸出步驟:輸出最終之經來自各方向所得到的模擬投影之數據所形成之三維影像重建影像。
其中,該感興趣區域係為一圓球體空間。
其中,該簡化數據模型係為一截圓影像資料。
其中,該影像系統係一影像偵檢器及相對應之一能量源,該物體係位在該影像偵檢器及該能量源之間。
其中,該能量源係為光源、輻射射源或非輻射射源。
其中,該比對步驟系採用以一投影資料中心為軸進行旋轉之方式對模擬投影之資訊以及該組三維影像重建空間資訊進行比對。
其中,該影像輸出步驟中所輸出之影像係為三維或多維斷層影像。
10‧‧‧物體
100‧‧‧能量源
13‧‧‧截圓區域資料
14‧‧‧數據
20‧‧‧影像偵檢器
21‧‧‧投影資料
22‧‧‧投影區域
23‧‧‧投影資料
30‧‧‧位置
31‧‧‧投影資料
32‧‧‧區域
33‧‧‧區域
θ‧‧‧角度
Φ‧‧‧角度
步驟S1~S7‧‧‧本發明之電腦斷層影像重建方法
第1圖係表示本發明電腦斷層影像重建方法中物體透過影像系統呈現投影資料之示意圖。
第2圖係表示本發明電腦斷層影像重建方法中簡化成二維幾何架構的示意圖。
第3圖係表示本發明電腦斷層影像重建方法中進行投影、反投影以及修正之示意圖。
第4圖係表示本發明電腦斷層影像重建方法在另一組旋轉角度θ後進行投影比對之示意圖。
第5圖係表示本發明電腦斷層影像重建方法在不同角度Φ時的示意圖。
第6圖係表示本發明電腦斷層影像重建方法的流程圖。
第7圖係表示一組感興趣區域(圓球體空間)之數據資料作為預定資料並經過二維系統矩陣所投影出之初始投影影像。
第8圖係為本發明之影像偵檢器所收集到的投影資料(泛指其他非起始之投影資料)。
第9圖係表示本發明所收取到的投影資料根據其物理特性轉換成物體衰減圖(minus log map)。
第10圖係表示本發明其中一三維弦波圖(sinogram)。
第11圖係表示透過本發明多方投影資料收集所輸出的影像(projenogram)。
第1圖係表示本發明電腦斷層影像重建方法中物體透過影像系統呈現投影資料之示意圖。
請參考第1圖,影像系統係可包括一能量源100及一影像偵檢器20;物體10係位在能量源100與影像偵檢器20之間,能量源100係可為光
源、輻射射源或者是非輻射射源等,一般係為X射線,而物體10係可為能量源100之吸收體或衰減體;以X射線之能量源為例,能量源100朝物體10投射,會在影像偵檢器20上形成一投影資料21,假設物體10係為一球體空間。
第2圖係表示本發明電腦斷層影像重建方法中簡化成二維幾何架構的示意圖。
請參考第2圖所示,將物體10的三維架構透過簡化且可重複使用之二維幾何架構來呈現,而物體10之重建空間係利用一二維系統矩陣在此平面上(二維)形成一截圓區域資料13。
第3圖係表示本發明電腦斷層影像重建方法中進行投影、反投影以及修正之示意圖。
請參考第3圖所示,透過截圓區域13進行投影,可在投影區域22產生一組虛擬投影資料,並與原投影資料在投影區域22進行比對,再將其差異進行反投影,並針對截圓區域資料13進行修正。
第4圖係表示本發明電腦斷層影像重建方法在另一組旋轉角度θ後進行投影比對之示意圖。
在上述流程完成後,係以相同的二維系統矩陣重複套用,繼續進行另一組以投影資料中心(即y軸)旋轉角度θ後之數據14,並進行投影與其對應之投影資料23進行比對,重複疊代修正之步驟,依此類推,直到完成投影資料21上所有的資料比對。
第5圖係表示本發明電腦斷層影像重建方法在不同角度Φ時的示意圖。
在不同角度Φ時,影像偵檢器20的位置30與其投影資料31上
對應之區域32及33,係根據前述重複之步驟以同樣之二維系統進行不同方位之資料比對,直到所有資料之誤差值小於一預定誤差值,即可輸出數據以顯示三維或多維斷層影像。
第6圖係表示本發明電腦斷層影像重建方法的流程圖。
請參考第6圖所示,本發明的電腦斷層影像重建方法之步驟如下。
收集資訊步驟S1:以影像系統在每隔一預定角度之一位置環繞一物體10,以物體10為中心收集物體10之不同方向之相對應的一投影資料21,並將投影資料21記錄成影像系統的一幾何參數,各方向的幾何參數係建立成一組三維影像重建空間資訊。
投影步驟S2:建構影像系統之虛擬數據模型並對其立體特性簡化成一簡化數據模型(即二維系統矩陣),並虛擬化一組感興趣區域之一數據資料當作一預測資料(如第7圖所示),將感興趣空間之數據資料透過簡化數據模型(即二維系統矩陣)進行模擬投影(如第8圖所示)。
比對步驟S3:將模擬投影所得到的資訊與影像系統所收集該組三維影像重建空間資訊進行差異比對,計算一差異值。
反投影步驟S4:上述的差異值再透過簡化數據模型(即二維系統矩陣)進行反投影之疊代運算。
修正步驟S5:根據上述的差異值修正一差異影像。
循環步驟S6:確認上述的差異值是否小於一預定差異值,若否的話,則回到投影步驟,若是的話,則下一組三維影像重建空間資訊開始依序進行投影步驟、比對步驟、反投影步驟以及修正步驟,直至各方
向的該組三維影像重建空間資訊依序完成投影步驟、比對步驟、反投影步驟以及修正步驟,且各方向的修正步驟中的差異值小於預定差異值。
影像輸出步驟S7:輸出最終之所有各方向所得到的模擬投影之數據所形成之影像。
其中,感興趣區域係為一圓球體空間;簡化數據模型係為一截圓影像資料13;比對步驟係採用以投影資料中心(如y軸)為軸進行旋轉之方式對模擬投影之資訊以及該組三維影像重建空間資訊進行比對。
第9圖係表示本發明所收取到的投影資料根據其物理特性轉換成物體衰減圖(minus log map)。
第10圖係表示本發明其中一三維弦波圖(sinogram),係指來自所有方向之物體衰減圖,亦即本發明中依照不同角度Φ依序堆疊而成的圖式。
第11圖係指透過本發明多方投影資料收集所輸出的影像(projenogram)。
請參考第11圖所示,其係將物體衰減圖之座標系(u,v)轉換成極座標(R,θ)並將其繪出即可;每一次經二維系統矩陣運算後的數值以及在影像偵檢器20上的投影便可透過此資料收集進行比對。
此外,若將此資料依照所對應之方向回填至影像偵檢器20所在之平面的話,係可還原原始影像之投影,如第9圖所示。
因此,根據上述的方法,係可有效地減少傳統三維或多維影像之重建時所建立龐大系統矩陣運算,本發明使用時僅需二維截平面之陣列資料(二維系統矩陣)來建立部分的系統矩陣以取代傳統建置三維或多維
影像之系統矩陣,透過其球座標之對稱特性減少運算時記憶體用量與儲存之空間量,達到計算量縮減之效果。
當前述係針對本發明之各實施例時,本發明之其他或進一步的實施例可被設計出而無須違反其基本範圍,且其基本範圍係由下列的申請專利範圍所界定。
雖然本發明以相關的較佳實施例進行解釋,但是這並不構成對本發明的限制。應說明的是,本領域的技術人員根據本發明的思想能夠構造出很多其他類似實施例,這些均在本發明的保護範圍之中。
步驟S1~S7‧‧‧本發明之電腦斷層影像重建方法
Claims (7)
- 一種電腦斷層影像重建方法,其步驟係包括:收集資訊步驟:以一影像系統在每隔一預定角度之一位置環繞一物體,以該物體為中心收集該物體之不同方向之相對應的一投影資料,並將該投影資料記錄成該影像系統的一幾何參數,各方向的該幾何參數係建立成一組三維影像重建空間資訊;投影步驟:建構該影像系統之虛擬數據模型並對其立體特性簡化成一簡化數據模型,並虛擬化一組感興趣區域之一數據資料當作一預測資料,最後將該感興趣空間之數據資料透過該簡化數據模型進行模擬投影;比對步驟:將該模擬投影所得到的資訊與該影像系統所收集各角度之投影影像資訊進行差異比對,計算一差異值;反投影步驟:該差異值再透過該簡化數據模型進行反投影之運算;修正步驟:根據該反投影後之差異值修正原始三維影像重建空間的資訊;循環步驟:確認該差異值是否小於一預定差異值,若否的話,則回到該投影步驟,若是的話,則下一組三維影像重建空間資訊開始依序進行該投影步驟、該比對步驟、該反投影步驟以及該修正步驟,直至各方向的該組三維影像重建空間資訊依序完成該投影步驟、該比對步驟、該反投影步驟以及該修正步驟,且該組三維影像重建空間的數據與前一次疊代之差異值小於該預定差異值;以及影像輸出步驟:輸出最終之經來自各方向所得到的模擬投影之數據所形成之三維影像重建影像。
- 依據申請專利範圍第1項所述的方法,其中,該感興趣區域係為一圓球體空間。
- 依據申請專利範圍第2項所述的方法,其中,該簡化數據模型係為一截圓影像資料。
- 依據申請專利範圍第1項所述的方法,其中,該影像系統係一影像偵檢器及相對應之一能量源,該物體係位在該影像偵檢器及該能量源之間。
- 依據申請專利範圍第4項所述的方法,其中,該能量源係為光源、輻射射源或非輻射射源。
- 依據申請專利範圍第1項所述的方法,其中,該比對步驟系採用以一投影資料中心為軸進行旋轉之方式對模擬投影之資訊以及該組三維影像重建空間資訊進行比對。
- 依據申請專利範圍第1項所述的方法,其中,該影像輸出步驟中所輸出之影像係為三維或多維斷層影像。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW102136721A TWI517093B (zh) | 2013-10-11 | 2013-10-11 | Computer tomography reconstruction method |
US14/285,045 US9858690B2 (en) | 2013-10-11 | 2014-05-22 | Computed tomography (CT) image reconstruction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW102136721A TWI517093B (zh) | 2013-10-11 | 2013-10-11 | Computer tomography reconstruction method |
Publications (2)
Publication Number | Publication Date |
---|---|
TW201514919A TW201514919A (zh) | 2015-04-16 |
TWI517093B true TWI517093B (zh) | 2016-01-11 |
Family
ID=52809662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW102136721A TWI517093B (zh) | 2013-10-11 | 2013-10-11 | Computer tomography reconstruction method |
Country Status (2)
Country | Link |
---|---|
US (1) | US9858690B2 (zh) |
TW (1) | TWI517093B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI676024B (zh) * | 2017-04-17 | 2019-11-01 | 美商雷神公司 | 用以組合彩色3d影像之系統及方法 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108882897B (zh) * | 2015-10-16 | 2022-01-25 | 瓦里安医疗系统公司 | 一种对患者进行成像的方法和系统 |
CN105551001B (zh) | 2015-12-11 | 2019-01-15 | 沈阳东软医疗系统有限公司 | 一种图像重建方法、装置及设备 |
CN105574904B (zh) | 2015-12-11 | 2019-01-11 | 沈阳东软医疗系统有限公司 | 一种图像重建方法、装置及设备 |
TWI639414B (zh) * | 2017-11-17 | 2018-11-01 | 財團法人國家同步輻射研究中心 | 電腦斷層掃描影像的對位方法 |
CN109223016B (zh) * | 2018-07-18 | 2022-06-07 | 江苏一影医疗设备有限公司 | Ct成像方法及成像设备 |
CN113782174B (zh) * | 2021-11-12 | 2022-02-18 | 极限人工智能有限公司 | Ct机的校验方法、装置、电子设备及可读存储介质 |
Family Cites Families (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5253171A (en) | 1990-09-21 | 1993-10-12 | General Electric Company | Parallel processing method and apparatus based on the algebra reconstruction technique for reconstructing a three-dimensional computerized tomography (CT) image from cone beam projection data |
US6751283B2 (en) | 2002-08-06 | 2004-06-15 | Koninklijke Philips Electronics, N.V. | Reconstruction method for tilted-gantry computed tomography |
AU2003262726A1 (en) | 2002-08-21 | 2004-03-11 | Breakaway Imaging, Llc | Apparatus and method for reconstruction of volumetric images in a divergent scanning computed tomography system |
US7787669B2 (en) | 2003-01-30 | 2010-08-31 | Koninklijke Philips Electronics N.V. | Reconstruction of local patient doses in computed tomography |
WO2004072905A1 (en) | 2003-02-14 | 2004-08-26 | Koninklijke Philips Electronics N.V. | System and method for helical cone-beam computed tomography with exact reconstruction |
DE602004004877T2 (de) | 2003-02-14 | 2007-11-08 | Koninklijke Philips Electronics N.V. | System und verfahren für exakte bildrekonstruktion für spiralkegelstrahl-compu tertomographie mit redundanten daten |
US6850585B2 (en) * | 2003-03-05 | 2005-02-01 | Ge Medical Systems Global Technology Company, Llc | Progressive updating approach for volumetric CT image reconstruction |
US7535988B2 (en) | 2003-05-27 | 2009-05-19 | Clean Earth Technologies, Llc | Method for fast image reconstruction with compact radiation source and detector arrangement using computerized tomography |
JP2007530126A (ja) | 2004-03-25 | 2007-11-01 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | ゲート式コーンビームctの動き補償を用いた再構成 |
JP2008510509A (ja) | 2004-08-20 | 2008-04-10 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 実際測定値及び仮想測定値からの対象物画像の再構成のためのコンピュータ断層撮影方法及びコンピュータ断層撮影装置 |
RU2009110163A (ru) | 2006-08-22 | 2010-09-27 | Конинклейке Филипс Электроникс, Н.В. (Nl) | Реконструкция компьютерной томографии по двум наклонным окружностям |
US8126108B2 (en) | 2007-04-18 | 2012-02-28 | Agency For Science, Technology And Research | Method and apparatus for reorientated resconstruction of computed tomography images of planar objects |
WO2010028027A1 (en) | 2008-09-03 | 2010-03-11 | Mayo Foundation For Medical Education And Research | Method for reconstruction in dual energy, dual source helical computed tomography |
US8326054B2 (en) | 2008-11-26 | 2012-12-04 | Wisconsin Alumni Research Foundation | Method for prior image constrained image reconstruction in cardiac cone beam computed tomography |
CN102362292A (zh) | 2009-03-26 | 2012-02-22 | 皇家飞利浦电子股份有限公司 | 用于计算机断层摄影图像重建的方法和设备 |
WO2011011684A2 (en) | 2009-07-24 | 2011-01-27 | The Research Foundation Of State University Of New York | Selection of optimal views for computed tomography reconstruction |
WO2011100723A2 (en) | 2010-02-12 | 2011-08-18 | The Regents Of The University Of California | Graphics processing unit-based fast cone beam computed tomography reconstruction |
EP2577608B1 (en) | 2010-05-27 | 2014-07-16 | Koninklijke Philips N.V. | Improved reconstruction for cone-beam computed tomography imaging with off-center flat panel detector |
US9536314B2 (en) * | 2010-10-20 | 2017-01-03 | Siemens Medical Solutions Usa, Inc. | Image reconstruction |
CN102456227B (zh) | 2010-10-28 | 2015-05-27 | 清华大学 | Ct图像重建方法及装置 |
-
2013
- 2013-10-11 TW TW102136721A patent/TWI517093B/zh active
-
2014
- 2014-05-22 US US14/285,045 patent/US9858690B2/en active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI676024B (zh) * | 2017-04-17 | 2019-11-01 | 美商雷神公司 | 用以組合彩色3d影像之系統及方法 |
Also Published As
Publication number | Publication date |
---|---|
US20150103968A1 (en) | 2015-04-16 |
TW201514919A (zh) | 2015-04-16 |
US9858690B2 (en) | 2018-01-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
TWI517093B (zh) | Computer tomography reconstruction method | |
US8000513B2 (en) | System and method for 3D time of flight PET forward projection based on an exact axial inverse rebinning relation in fourier space | |
US8983162B2 (en) | Method and apparatus for estimating monte-carlo simulation gamma-ray scattering in positron emission tomography using graphics processing unit | |
Scheins et al. | Analytical calculation of volumes-of-intersection for iterative, fully 3-D PET reconstruction | |
WO2020237873A1 (zh) | 基于神经网络的螺旋ct图像重建方法和设备及存储介质 | |
JP2009534631A (ja) | ダーティアイソトープのpet再構成 | |
Pang et al. | Accelerating simultaneous algebraic reconstruction technique with motion compensation using CUDA-enabled GPU | |
US7769217B2 (en) | Fast iterative 3D PET image reconstruction using a set of 2D linogram transformations | |
Matej et al. | Image-based modeling of PSF deformation with application to limited angle PET data | |
Feng et al. | Use of three-dimensional Gaussian interpolation in the projector/backprojector pair of iterative reconstruction for compensation of known rigid-body motion in SPECT | |
CN103908278A (zh) | 图像重建方法及装置、医疗成像系统 | |
US20060098857A1 (en) | Method for faster iterative reconstruction for converging collimation spect with depth dependent collimator response modeling | |
KR101283266B1 (ko) | Gpu를 이용한 양전자 방출 단층 촬영 영상에서의 몬테카를로 시뮬레이션 감마선 산란 추정 방법 및 장치 | |
JP2010204755A (ja) | 画像処理装置、画像再構成システム、画像処理方法およびプログラム | |
TWI509564B (zh) | 三維成像的投影方法 | |
JP7273272B2 (ja) | 角度オフセットによる断層画像データの取得方法、取得装置、および制御プログラム | |
Miao | Comparative studies of different system models for iterative CT image reconstruction | |
Li et al. | Parallel computing approach for efficient 3-D X-ray-simulated image reconstruction | |
CN105832358B (zh) | 一种基于系统校准的旋转双平板pet系统的成像方法 | |
Liu et al. | Cooperative data fusion of transmission and surface scan for improving limited-angle computed tomography reconstruction | |
Nguyen et al. | GPU-accelerated iterative reconstruction from Compton scattered data using a matched pair of conic projector and backprojector | |
Lékó et al. | Scale invariance in projection selection using binary tomography | |
Blackston et al. | Progress Update on Iterative Reconstruction of Neutron Tomographic Images, 2019 | |
US9279887B2 (en) | Self-similar, transversely shift-invariant point response function in single photon emission computed tomography | |
Fan et al. | A block-wise approximate parallel implementation for ART algorithm on CUDA-enabled GPU |