一種計(jì)算機(jī)層析成像設(shè)備和方法
【專利摘要】本發(fā)明公開了一種對(duì)被檢測(cè)對(duì)象的感興趣區(qū)域進(jìn)行CT成像的方法，包括步驟：獲取所述感興趣區(qū)域的CT投影數(shù)據(jù)；獲取B區(qū)域的CT投影數(shù)據(jù)；選擇一組覆蓋所述感興趣區(qū)域的PI線段，并且為所述PI線段組中的每條PI線段計(jì)算其上的重建圖像值；然后組合所述PI線段組中的所有PI線段上的重建圖像值而得到所述感興趣區(qū)域的圖像。本發(fā)明還公開了采用該方法的CT成像設(shè)備及其中的數(shù)據(jù)處理單元。由于只需要X射線束覆蓋感興趣區(qū)域和B區(qū)域即能夠精確重建獲得該感興趣區(qū)域的二維/三維斷層圖像，因此可以使用較小尺寸的探測(cè)器實(shí)現(xiàn)大尺寸物體任意位置感興趣區(qū)域的CT成像，這能夠大大降低CT掃描過程中的X射線輻射劑量。
【專利說(shuō)明】一種計(jì)算機(jī)層析成像設(shè)備和方法
[0001]本申請(qǐng)是申請(qǐng)日為2009年5月26日、申請(qǐng)?zhí)枮?00910085612.4以及發(fā)明名稱為“一種計(jì)算機(jī)層析成像設(shè)備和方法”的發(fā)明專利申請(qǐng)的分案申請(qǐng)。
【技術(shù)領(lǐng)域】
[0002]本發(fā)明涉及一種計(jì)算機(jī)層析成像(CT)設(shè)備，尤其是一種用于對(duì)感興趣(ROI)區(qū)域進(jìn)行CT成像的設(shè)備，以及在這些設(shè)備中使用的CT成像方法。
【背景技術(shù)】
[0003]自從1972年Hounsfield發(fā)明了第一臺(tái)CT機(jī)，CT技術(shù)給醫(yī)學(xué)診斷和工業(yè)無(wú)損檢測(cè)帶來(lái)了革命性的影響，CT已經(jīng)成為醫(yī)療、生物、航空航天、國(guó)防等行業(yè)重要的檢測(cè)手段之一。X射線錐束CT已經(jīng)在醫(yī)學(xué)臨床、安全檢查、無(wú)損檢測(cè)等領(lǐng)域得到了廣泛的應(yīng)用，特別是在醫(yī)學(xué)臨床診斷中，螺旋CT已經(jīng)成為不可或缺的檢查手段之一 [1]。
[0004]盡管目前CT技術(shù)已經(jīng)在工業(yè)、安檢、醫(yī)療等領(lǐng)域取得了巨大的成功，但由于工程應(yīng)用條件的復(fù)雜性和多樣性，對(duì)CT技術(shù)的進(jìn)一步發(fā)展提出了更高的要求。尤其是在工業(yè)應(yīng)用中，CT技術(shù)在大尺寸、高精度成像，醫(yī)療低劑量成像等方面存在較大困難。這主要由于:CT掃描視野(FOV)受X射線束寬度、探測(cè)器尺寸和掃描視角的限制，這使得對(duì)大物體掃描的投影數(shù)據(jù)可能存在探測(cè)器方向和掃描角度方向兩個(gè)方面的截?cái)啵欢F(xiàn)有的主流CT算法都是針對(duì)完整物體的全局重建方法，要求X射線束必須完全覆蓋物體斷層，對(duì)投影數(shù)據(jù)有截?cái)嗟那闆r難以處理。因此，目前對(duì)大物體或不規(guī)則物體進(jìn)行成像時(shí)，往往難以直接完成掃描，需要通過多次掃描后依靠數(shù)據(jù)重排等近似轉(zhuǎn)換方法才能重建出最終圖像，這對(duì)CT成像的速度和精度造成了負(fù)面的影響。
[0005]此外，CT設(shè)備中的探測(cè)器已經(jīng)成為限制CT設(shè)備的硬件成本的最關(guān)鍵環(huán)節(jié)，探測(cè)器的價(jià)格直接和探測(cè)器單元尺寸和數(shù)量成正比，而探測(cè)器價(jià)格的居高不下，這極大地限制了CT產(chǎn)品的成本空間。
[0006]另一方面，在醫(yī)療CT成像中，為了保證X射線投影數(shù)據(jù)不截?cái)啵壳癈T掃描使用的X射線束寬度必須覆蓋人體斷層的寬度方向，而真正感興趣往往只是人體的某個(gè)器官，這大大增加了人體在CT掃描過程中受到的不必要的輻射劑量，如果不改變目前CT設(shè)計(jì)的思路，很難降低輻射劑量。目前，醫(yī)療照射已經(jīng)成為全民最大的人為電離輻射來(lái)源，減少CT檢查的X射線劑量是關(guān)系到全體社會(huì)公眾及其后代健康的重大課題。[0007]另外，隨著人民生活水平的提高，人們對(duì)醫(yī)學(xué)診斷的要求越來(lái)越高，特別是對(duì)于某些人體特殊部位，例如:女性乳腺、耳蝸、牙齒等部位的醫(yī)學(xué)診斷需要很高的空間分辨率。對(duì)于這些情況，目前的主流全身螺旋CT機(jī)無(wú)法滿足正常醫(yī)療診斷的需求。隨著大面積平板探測(cè)器技術(shù)的飛速發(fā)展，目前醫(yī)用平板探測(cè)器技術(shù)已經(jīng)相當(dāng)成熟，如已經(jīng)應(yīng)用于高空間分辨率的X射線DR成像。目前基于非晶硒、非晶硅技術(shù)的平板探測(cè)器，其有效探測(cè)器面積超過了 500_X500mm，而探測(cè)器像素尺寸也達(dá)到了 0.1mm左右，而目前多層螺旋CT的探測(cè)器單位尺寸只有約0.5mm，利用平板探測(cè)器可以實(shí)現(xiàn)比現(xiàn)有螺旋CT空間分辨率高的多的CT圖像。但是，由于平板探測(cè)器數(shù)據(jù)量巨大，數(shù)據(jù)傳輸速度還無(wú)法滿足目前錐束CT成像的要求，目前還沒有真正應(yīng)用平板探測(cè)器進(jìn)行全身成像的CT設(shè)備。如果我們不是針對(duì)人體全身，而是利用小型平板探測(cè)器對(duì)感興趣的器官實(shí)現(xiàn)高精度感興趣區(qū)域CT成像，那么就能夠克服數(shù)據(jù)傳輸速度慢的技術(shù)瓶頸，使之成為可能。
[0008]針對(duì)CT系統(tǒng)的上述缺點(diǎn)和局限性，我們開始探索新的CT圖像重建方法和CT成像模式。事實(shí)上，很多工程應(yīng)用中并不要求對(duì)完整的物體進(jìn)行全局CT成像，只需要獲得某些感興趣區(qū)域(ROI)的物體圖像即可，特別是醫(yī)療臨床診斷中，只要能夠?qū)崿F(xiàn)對(duì)可疑病灶部位的成像即可[2]。
[0009]20世紀(jì)80年代初，人們已經(jīng)開始研究物體的局部CT成像問題，但受當(dāng)時(shí)CT重建理論所限，人們無(wú)法精確重建出物體局部的CT圖像，因此轉(zhuǎn)而尋求一種與物體斷層圖像相關(guān)的近似函數(shù)。1985年,Smith等提出了一種lambda tomography的局部重建算法,該算法利用局部投影數(shù)據(jù)重建出一個(gè)與ROI密度函數(shù)有相同奇異性的函數(shù)[3]。隨后，Katsevich提出了一種pseudolocal tomography的局部重建算法，這種算法通過重建密度函數(shù)的Hilbert變換的一部分來(lái)代替原函數(shù)M。但是，由于這些函數(shù)都無(wú)法代替真正的物體斷層圖像，難以滿足實(shí)際工程應(yīng)用的需求，極大地降低了局部CT成像對(duì)于實(shí)際工程應(yīng)用的意義。在隨后很長(zhǎng)的一段時(shí)間里，針對(duì)物體局部ROI的成像研究一直陷于停頓，無(wú)法得到解決。
[0010]最近幾年，CT重建理論獲得了較大的發(fā)展。2002年，Katsevich首先提出了一種基于濾波反投影(Filtered Back Projection,FBP)形式的錐束螺旋CT精確重建算法。該算法很好地解決了長(zhǎng)物體重建問題，在Z軸方向投影數(shù)據(jù)截?cái)嗟那闆r下，該算法仍然能夠精確地重建出被掃描部分的物體圖像，并且由于該算法是FBP形式的，在重建速度上和迭代類重建算法相比具有很大的優(yōu)勢(shì)，Katsevich為CT算法的發(fā)展提供了一個(gè)嶄新的思路[5_6]。2004 年,Zou 和 Pan 提出了一種反投影濾波(Back Projection Filtration,BPF)形式的螺旋CT精確重建算法，該方法 只需要理論上最少的投影數(shù)據(jù)就能夠精確重建出物體的斷層圖像[7_8]。至此，螺旋CT重建的基本理論問題得到了較好地解決。隨后，BPF算法被推廣應(yīng)用到平行束、扇束、錐束CT圖像重建中。Zou和Pan的BPF算法是基于PI線的重建算法，PI線是掃描軌跡上任意兩點(diǎn)連線上的一條線段，BPF算法要求每條PI線的兩個(gè)端點(diǎn)在物體支撐外，它的最大優(yōu)勢(shì)在于當(dāng)投影數(shù)據(jù)存在某些截?cái)鄷r(shí)仍然能夠精確重建出PI線上的圖像，這使得針對(duì)ROI而不是完整物體的CT重建成為可能。2006年，Defrise等在BPF算法的基礎(chǔ)上取得了更進(jìn)一步的成果，放寬了 PI線的限制，證明了當(dāng)PI線只有一個(gè)端點(diǎn)在物體支撐外時(shí)，由通過該P(yáng)I線的投影數(shù)據(jù)仍然能夠精確重建出PI線上圖像[9]。2007年，Wang等人進(jìn)一步證明了當(dāng)PI線完全在物體內(nèi)部時(shí)，如果已知該P(yáng)I線上一部分的圖像信息，就能夠通過截?cái)嗟耐队皵?shù)據(jù)精確重建出該P(yáng)I線上的物體圖像[1°]。但是，在實(shí)際CT工程應(yīng)用中，物體內(nèi)部PI線上的重建數(shù)值信息很難事先獲得，因此，Wang等人的方法在實(shí)際應(yīng)用中具有一定的局限性。
[0011]相關(guān)文獻(xiàn):
[1]李亮.CT投影變換理論及錐束重建方法研究[博士學(xué)位論文].北京:清華大學(xué)工程物理系，2007.[2]李亮，陳志強(qiáng)，康克軍，張麗，邢宇翔.感興趣區(qū)域CT圖像重建方法及模擬實(shí)驗(yàn)，CT理論與應(yīng)用研究，18: 1-7，2009.[3]AFaridanij E L Ritmanj K T Smith.“Local tomography”，SIAM Appl Math,52: 459-484，1992.[4]AKatsevichj A.G.Ramm.“Pseudolocal tomography”， SIAM Appl Math, 56:167-191，1996.[5]KatsevichA.A general scheme for constructing inversion algorithms forcone beam CT.1nt J Math Math Scij 2003， 21:1305-1321.[6]KatsevichA.An improved exact filtered backprojection algorithm forspiral computed tomography.Adv Appl Math, 2004， 32:681-697.[7]ZouYj Pan X.Exact image reconstruction on P1-1ines from minimum datain helical cone-beam CT.Phys Med Biol, 2004， 49:941-959.[8]ZouYj Pan X.1mage reconstruction on P1-1ines by use of filteredbackprojection in helical cone-beam CT.Phys Med Biol, 2004， 49:2717-2731.[9]MDefrisej F Nooj R Clackdoylej H Kud0.“Truncated Hilbert transformand image reconstruction from limited tomographic data”， Inverse Problems, 22:1037-1053， 2006.[10]YB Yej H Y Yuj Y Weij G Wang.“A General Local Reconstruction ApproachBased on a Truncated Hilbert Transform，，， International Journal of BiomedicalImaging, Volume 2007， Article ID 63634。
【發(fā)明內(nèi)容】

[0012]本發(fā)明的 申請(qǐng)人:發(fā)現(xiàn)了一個(gè)全新的CT精確重建數(shù)據(jù)完備性條件:
對(duì)于任意一點(diǎn)3 = (COS病，sm為)，如果投影數(shù)據(jù)同時(shí)滿足下述條件，則圖像函數(shù)/(F)可以被精確重建:
1.存在向量友=(-sin為，cos病)和兩條線段.￡/C L、Ln [L ,其中l(wèi)是通過毛點(diǎn)且平行于牙的在函數(shù)/(O支撐內(nèi)的線段，；
2.線段4至少有一個(gè)端點(diǎn)在函數(shù)支撐外，或線段^上的/(F)已知；
3.對(duì)于任意點(diǎn)IfUA，通過包含該點(diǎn)的微小鄰域內(nèi)任意點(diǎn)的任意角度癸e [爲(wèi)肩+ π]下的投影P (s,釣均被采集得到，其中s ?。
[0013]上述條件可以由圖1解釋，當(dāng)感興趣區(qū)域ROI位于物體支撐內(nèi)部時(shí)，F(xiàn)OV僅僅覆蓋ROI是不能精確重建ROI圖像的，需要增加額外的投影信息才能夠重建出該ROI的斷層圖像，分兩種情況:(a) (b)屬于一類情況，即增加一個(gè)額外的深色陰影區(qū)域，通過該陰影區(qū)域的投影數(shù)據(jù)也需要采集得到，并且該深色陰影區(qū)域的圖像已知；(c) (d)屬于另外一類情況，即增加額外的一個(gè)陰影區(qū)域B，該區(qū)域有一部分或全部在物體支撐外。只要滿足了上面兩種情況之一，都能夠?qū)崿F(xiàn)感興趣區(qū)域ROI的精確CT成像。對(duì)比圖1中的四種情況，我們可以發(fā)現(xiàn)對(duì)于后面兩種情況(c)(d)由于不需要已知該額外陰影區(qū)域的圖像信息，因此這兩種情況更方便在工程應(yīng)用中實(shí)現(xiàn)；特別是最后一種情況(d)由于陰影區(qū)域B位于物體支撐外，因而只需要測(cè)量到通過ROI和區(qū)域B的X射線投影即可精確重建ROI的斷層圖像，這就意味著該CT掃描方式對(duì)于醫(yī)療成像而言，被檢測(cè)者能夠接受更少的輻射劑量。被檢測(cè)者只需要接受掃描其可疑組織區(qū)域的X射線劑量和通過其身體外一個(gè)很小空氣區(qū)域的X射線劑量即可，相比目前的全身螺旋CT掃描，能夠在很大程度上減小被檢測(cè)者接受的輻射劑量。
[0014]本發(fā)明正是基于上述CT精確重建數(shù)據(jù)完備性條件，尤其是基于上述情況(d)而做出。
[0015]根據(jù)本發(fā)明的一個(gè)方面，提供了一種對(duì)被檢測(cè)對(duì)象的感興趣區(qū)域進(jìn)行CT成像的方法，包括步驟:獲取所述感興趣區(qū)域的CT投影數(shù)據(jù)；獲取B區(qū)域的CT投影數(shù)據(jù)，其中所述B區(qū)域的至少一部分位于所述被檢測(cè)對(duì)象的支撐之外，而且所述B區(qū)域被選為使得能夠選擇一組覆蓋所述感興趣區(qū)域的PI線段，其中每條經(jīng)過所述感興趣區(qū)域的PI線段均通過B區(qū)域；以及根據(jù)所述感興趣區(qū)域的CT投影數(shù)據(jù)和所述B區(qū)域的CT投影數(shù)據(jù)來(lái)重建所述感興趣區(qū)域的CT投影數(shù)據(jù)。
[0016]根據(jù)本發(fā)明的另一方面，提供了一種用于對(duì)被檢測(cè)對(duì)象的感興趣區(qū)域進(jìn)行CT成像的CT成像設(shè)備，包括:X射線發(fā)生器，用于產(chǎn)生用于掃描的X射線束；探測(cè)器裝置，用于探測(cè)透過被掃描區(qū)域的X設(shè)備以產(chǎn)生投影數(shù)據(jù)；被檢測(cè)對(duì)象承載裝置，用于承載所述被檢測(cè)對(duì)象進(jìn)出所述CT成像設(shè)備；主控制器，用于控制所述CT成像設(shè)備的操作，以使得所述X射線發(fā)生器所發(fā)出的X射線僅僅覆蓋所述感興趣區(qū)域來(lái)進(jìn)行CT掃描以獲得所述感興趣區(qū)域的投影數(shù)據(jù)，以及使得所述X射線發(fā)生器所發(fā)出的X射線僅僅覆蓋B區(qū)域來(lái)對(duì)所述B區(qū)域進(jìn)行CT掃描以獲得所述B區(qū)域的投影數(shù)據(jù)，所述B區(qū)域的至少一部分位于所述被檢測(cè)對(duì)象的支撐之外，而且所述B區(qū)域被選為使得能夠選擇一組覆蓋所述感興趣區(qū)域的PI線段，其中每條經(jīng)過所述感興趣區(qū)域的PI線段均通過B區(qū)域；以及數(shù)據(jù)處理單元，用于基于所獲得的感興趣區(qū)域的投影數(shù)據(jù)和B區(qū)域的投影數(shù)據(jù)來(lái)重建所述感興趣區(qū)域的圖像。
[0017]根據(jù)本發(fā)明提出的X射線CT成像設(shè)備和方法，只需要X光源發(fā)出的X射線束覆蓋感興趣區(qū)域即能夠精確重建獲得該感興趣區(qū)域的二維/三維斷層圖像，因此可以使用較小尺寸的探測(cè)器實(shí)現(xiàn)大尺寸物體任意位置感興趣區(qū)域的CT成像，提高了 CT掃描和圖像重建的速度，降低了 CT設(shè)備的硬件成本，同時(shí)能夠大大降低CT掃描過程中的X射線輻射劑量。因此，本發(fā)明具有很高的市場(chǎng)應(yīng)用潛`力。
【專利附圖】

【附圖說(shuō)明】
[0018]通過閱讀下文優(yōu)選實(shí)施方式的詳細(xì)描述，各種其他的優(yōu)點(diǎn)和益處對(duì)于本領(lǐng)域普通技術(shù)人員將變得清楚明了。附圖僅用于示出優(yōu)選實(shí)施方式的目的，而并不認(rèn)為是對(duì)本發(fā)明的限制。而且在整個(gè)附圖中，用相同的參考符號(hào)表示相同的部件。在附圖中:
圖1說(shuō)明了根據(jù)本發(fā)明的CT精確重建數(shù)據(jù)完備性條件；
圖2描述了根據(jù)本發(fā)明一個(gè)實(shí)施例的計(jì)算機(jī)層析(CT)成像方法的流程圖；
圖3描述了根據(jù)本發(fā)明一個(gè)實(shí)施例的、在圖2的CT成像方法中執(zhí)行的投影數(shù)據(jù)重建方法的流程圖；
圖4示出了根據(jù)本發(fā)明一個(gè)實(shí)施例的CT成像設(shè)備400的結(jié)構(gòu)框圖；
圖5示出了根據(jù)本發(fā)明一個(gè)實(shí)施例的CT成像設(shè)備400中的前準(zhǔn)直裝置420的結(jié)構(gòu)框
圖；
圖6示出了根據(jù)本發(fā)明一個(gè)實(shí)施例的CT成像設(shè)備400中的探測(cè)器位置控制裝置440的結(jié)構(gòu)框圖；以及
圖7示出了根據(jù)本發(fā)明一個(gè)實(shí)施例的CT成像設(shè)備進(jìn)行CT成像的數(shù)值模擬結(jié)果。
具體實(shí)施例
[0019]下面結(jié)合附圖和具體的實(shí)施方式對(duì)本發(fā)明作進(jìn)一步的描述。
[0020]如上面關(guān)于CT精確重建數(shù)據(jù)完備性條件的情況(d)所述，只需要通過掃描獲得了通過ROI和B區(qū)域的X射線投影即可精確重建ROI的斷層圖像。
[0021]圖2描述了根據(jù)本發(fā)明的一個(gè)實(shí)施例的計(jì)算機(jī)層析(CT)成像方法。在步驟S201中，確定要進(jìn)行CT成像的被檢測(cè)對(duì)象中的感興趣(ROI)區(qū)域，并且根據(jù)所確定的ROI區(qū)域來(lái)確定至少一部分位于被檢測(cè)對(duì)象支撐之外的B區(qū)域。B區(qū)域被選為使得能夠選擇一組覆蓋ROI區(qū)域的PI線段，其中保證每條經(jīng)過ROI區(qū)域的PI線均通過B區(qū)域(如圖1 (d)所示的情況)。B區(qū)域可以是一個(gè)任意大小和形狀的空間區(qū)域，考慮到實(shí)際探測(cè)器單元的尺寸和病人接受的輻射劑量?jī)蓚€(gè)因素，一般地，B區(qū)域可以選擇為直徑是探測(cè)器單元10倍的圓或球形空間區(qū)域。當(dāng)然，B區(qū)域可以不止一個(gè)。另外，為了進(jìn)一步提高圖像重建速度和精度，可以進(jìn)一步確定其它參考區(qū)域。
[0022]在確定了 ROI區(qū)域和B區(qū)域之后，在步驟S203中，將ROI區(qū)域調(diào)整到掃描視野范圍內(nèi)，以使得掃描X射線束僅僅覆蓋該ROI，并對(duì)該ROI區(qū)域進(jìn)行CT掃描以獲得ROI區(qū)域的CT投影數(shù)據(jù)。然后在步驟S205中，將B區(qū)域調(diào)整到掃描視野范圍內(nèi)，以使得X射線束僅僅覆蓋該B區(qū)域，并對(duì)該B區(qū)域進(jìn)行CT掃描以獲得B區(qū)域的CT投影數(shù)據(jù)。隨后在步驟S207中，基于在步驟S203獲得的ROI區(qū)域的CT投影數(shù)據(jù)和在步驟S205獲得的B區(qū)域的CT投影數(shù)據(jù)，根據(jù)下面如參考圖3所述的重建方法來(lái)重建出該ROI區(qū)域的CT圖像。在本發(fā)明的CT成像方法中，雖然需要進(jìn)行兩次CT掃描才能夠重建出ROI區(qū)域的CT圖像，但在每次掃描中，X射線束僅僅覆蓋到ROI和B區(qū)域，因此兩次掃描中被檢測(cè)對(duì)象受到的總X射線輻射劑量遠(yuǎn)遠(yuǎn)小于傳統(tǒng)CT掃描的福射劑量。
[0023]圖3具體描述了在圖2所示的步驟S207中執(zhí)行的重建方法。在步驟S301中，選擇一組能夠覆蓋ROI區(qū)域的PI線段，且保證每條經(jīng)過ROI區(qū)域的PI線均通過B區(qū)域。
[0024]該方法隨后進(jìn)入步驟S303-S307，其中對(duì)于在步驟S301中選定的PI線段組中的每條PI線，沿著該P(yáng)I線進(jìn)行圖像重建。然后在步驟S309中，對(duì)在所有PI線上的重建數(shù)據(jù)進(jìn)行組合而得到最終的ROI圖像。這種基于PI線的圖像重建過程在二維和三維重建中是一樣的，下面以二維圖像重建為例來(lái)說(shuō)明在步驟S303-S307中執(zhí)行的、基于一條PI線的圖像重建處理。
[0025]將沿著一條PI線的二維圖像定義為函數(shù)/(?) €i2(P)，則圖像重建處理就是尋找符合下面公式(I)所有5個(gè)約束的交集函數(shù)。
【權(quán)利要求】
1.一種對(duì)被檢測(cè)對(duì)象的感興趣區(qū)域進(jìn)行CT成像的方法，包括步驟:同時(shí)獲取所述感興趣區(qū)域的CT投影數(shù)據(jù)和B區(qū)域的CT投影數(shù)據(jù)，其中所述B區(qū)域的至少一部分位于所述被檢測(cè)對(duì)象的支撐之外，而且所述B區(qū)域被選為使得能夠選擇一組覆蓋所述感興趣區(qū)域的PI線段，其中每條經(jīng)過所述感興趣區(qū)域的PI線段均通過B區(qū)域；以及根據(jù)所述感興趣區(qū)域的CT投影數(shù)據(jù)和所述B區(qū)域的CT投影數(shù)據(jù)來(lái)重建所述感興趣區(qū)域的CT投影數(shù)據(jù)。
2.如權(quán)利要求1所述的CT成像方法，其中所述根據(jù)所述感興趣區(qū)域的CT投影數(shù)據(jù)和所述B區(qū)域的CT投影數(shù)據(jù)來(lái)重建所述感興趣區(qū)域的CT投影數(shù)據(jù)包括步驟: 對(duì)于所述PI線段組中的每條PI線段: 計(jì)算在所述PI線段上的一維Hilbert變換值，以及 對(duì)所計(jì)算的一維Hilbert變換值進(jìn)行有限Hilbert逆變換，以獲得在所述PI線段上的重建圖像值； 組合所述PI線段組中的所有PI線段上的重建圖像值而得到所述感興趣區(qū)域的圖像。
3.如權(quán)利要求2所述的CT成像方法，其中所述有限Hilbert逆變換步驟包括: 在重建圖像域和Hilbert變換空間域兩個(gè)域之間進(jìn)行凸集投影迭代來(lái)獲得滿足精度要求的所述在PI線段上的重建圖像。
4.如權(quán)利要求1所述的CT成像方法， 其中所述獲取所述感興趣區(qū)域的CT投影數(shù)據(jù)和所述獲取B區(qū)域的CT投影數(shù)據(jù)包括: 利用平行束、扇束或者錐束X射線來(lái)對(duì)所述感興趣區(qū)域和B區(qū)域進(jìn)行CT掃描，以獲得所述感興趣區(qū)域和B區(qū)域的CT投影數(shù)據(jù)。
【文檔編號(hào)】A61B6/06GK103584877SQ201310546294
【公開日】2014年2月19日 申請(qǐng)日期:2009年5月26日 優(yōu)先權(quán)日:2009年5月26日
【發(fā)明者】李亮, 陳志強(qiáng), 張麗, 康克軍, 李元景, 劉以農(nóng), 邢宇翔, 趙自然 申請(qǐng)人:清華大學(xué), 同方威視技術(shù)股份有限公司