CN1210000C - 在利用仅覆盖一半视场的减小尺寸的检测器的计算机x-射线断层成像***中应用的方法和装置 - Google Patents
在利用仅覆盖一半视场的减小尺寸的检测器的计算机x-射线断层成像***中应用的方法和装置 Download PDFInfo
- Publication number
- CN1210000C CN1210000C CN 00800601 CN00800601A CN1210000C CN 1210000 C CN1210000 C CN 1210000C CN 00800601 CN00800601 CN 00800601 CN 00800601 A CN00800601 A CN 00800601A CN 1210000 C CN1210000 C CN 1210000C
- Authority
- CN
- China
- Prior art keywords
- projection
- detector
- data
- ray
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 claims description 51
- 238000002591 computed tomography Methods 0.000 claims description 44
- 230000000007 visual effect Effects 0.000 claims description 22
- 230000005855 radiation Effects 0.000 claims description 18
- 230000007704 transition Effects 0.000 claims description 18
- 238000003384 imaging method Methods 0.000 claims description 17
- 238000002601 radiography Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims 1
- 238000009499 grossing Methods 0.000 abstract description 5
- 238000005303 weighing Methods 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 238000013170 computed tomography imaging Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000003491 array Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000010408 sweeping Methods 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010010 raising Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Pulmonology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Facsimile Scanning Arrangements (AREA)
- Image Input (AREA)
- Image Processing (AREA)
Abstract
描述了一种获得对象的投影数据的计算机X-射线断层成像(CT)***和方法,其包括X-射线源和检测器。该检测器相对于与CT***的旋转中心在检测器上的投影相对应的中心位置移动了其宽度的一半。依据本发明的方法,对于每个投影视图,选择表示最靠近成像***的旋转中心在检测器上的投影的位置的检测器元件的值Va。然后,对于所选择的检测器元件,检测器元件值Vb是从在相同方向所重构的数据的正向投影或相反的方向中估计的。然后选择能够消除Va和Vb之差的平滑函数。应用加权的平滑函数以消除Va和Vb之差。当将真实的投影数据和所估计的投影数据相结合时应用该过程来减小幅值差,由此产生平滑的过渡区。
Description
本申请要求分别以在1999年4月15日和1999年11月19日申请的申请号为No.60/129,398和60/166,500的临时申请的申请日为优先权,在此以引用的方式将它们全部结合在本申请中。
本发明涉及一种应用在体积型计算机X-射线断层成像(VCT)***中应用的方法和装置,该***应用一种尺寸减小的面积检测器,这种面积检测器仅覆盖一半的视场,由此降低这种面积检测器的尺寸和成本但不增加或基本不增加假象。
计算机X-射线断层成像(CT)是这样的一种技术:一般包含对患者进行X-射线辐射、采集患者的部分身体的数字X-射线投影数据以及处理和背式投影数字X-射线投影数据以产生图象,然后将该图象显示在CT***的监视器上。CT***通常包括台架、工作台、X-射线管、X-射线检测器阵列、计算机和显示检测器。计算机给台架的控制器发送指令以使台架使X-射线管和/或检测器阵列以特定的转速转动。
在第三代CT***中,在部分由检测器阵列和X-射线管组成的台架和患者的身体之间产生相对转动。由于产生这种相对转动,计算机控制通过X-射线管和检测器阵列执行的数据采集过程以采集数字X-射线照相。然后计算机进行处理并通过执行重构算法背向投影数字X-射线照相数据,并在显示监视器上显示所重构的CT图像。
如今所应用的许多CT***都利用在台架中的单行检测器,这种单行的检测器通常称为检测器元件的线性阵列。更先进的CT***应用两至四个线性检测器阵列以构成多行检测器。虽然这两者检测器结构都可以用于螺旋扫描方案,但是由于通过增加检测器阵列的螺旋间距多行检测器能够在更少的时间中扫描患者特定的轴线区域,所以它有利于患者扫描。螺旋间距通常定义为在台架旋转一圈中支撑患者的工作台的位移与检测器间距之比。例如,一个螺旋间距是指在CT***的CT台架的旋转一圈中将患者工作台移动等于检测器间距的量。
通常,线性检测器或多行检测器阵列覆盖由X-射线源发射的X-射线扇形束的整个视场。换句话说,通过检测器阵列吸收穿过或照射所扫描的对象的面积的X-射线,该对象可能是或不是患者。
在CT成像***中,比较理想的是并且在某些情况下也是必需的是减小检测器阵列的尺寸。例如,在新近发展的CT技术中应用包括许多行线性检测器阵列的面积检测器阵列进行CT数据采集。当前,仍然还没有能够覆盖整个成像的视场或患者范围的检测器面板。此外,一些应用线性检测器阵列的***支持对于所扫描的患者的很大的视场。理想的是在这种情况下也减小检测器阵列的大小和成本。
用于克服这些局限性的一种方法是将相对较小的检测器阵列平移其宽度的一半。例如,假设为覆盖患者的所需的视场的检测器阵列的最初的尺寸应该是80厘米。可以应用等于最初检测器的宽度一半的更小的检测器,即在这种情况为40厘米。这种检测器偏移它的一半的宽度(在这种情况为20厘米)以使它大致覆盖CT成像***的视场的一半。在本实例中,通过宽度等于它的最初宽度值的一半的检测器获得了在患者上的相同的视场。
还能够增加具有固定宽度的检测器的***的视场。通常,CT成像***的中心旋转的投影与检测器面板的中心对准。在CT成像中的旋转中心是X-射线源和检测器阵列绕其旋转的点的物理位置。然而,通过使检测器相对于它的原始位置偏移它的一半宽度能够增加这种***的视场(FOV)。虽然检测器仍然测量穿过成像***的物理旋转中心(即,在ISO)的X-射线的投影数据,但是成像***的旋转中心的投影是在已经平移的线性检测器、多行检测器或面积检测器的边沿附近。这种结构反过来又有效地使原始成像***结构的视场加倍,这就能够较大地增加成像***的视场。将检测器移动它的一半宽度的***结构通常称为半检测器移动。
在扇形束CT***中,在这种CT***中X射线源是辐射具有孔径张角的X-射线的点,该具有孔径张角的X-射线仅辐射检测器面板并类似于扇形,需要采集CT台架的整个旋转的一部分旋转的投影数据。具体地说,需要在台架绕患者的圆形轨道旋转180°加上该扇形角的角度间距的同时采集投影数据。在再一次测量中,扇形角X-射线的孔径张角的度量,具有该孔径张角的X-射线仅辐射在成像***的轴向平面中的检测器阵列。可以清楚地看到,由于不需要在台架绕患者旋转的整个360°的过程中测量投影,一些投影数据必定是冗余的。
在CT***的半检测器移动结构中,在台架的整个360°旋转中采集数据。在每个视角或台架上,仅测量一半的投影数据。应用来自台架的其它视图的数据来完成在给定的视角上的投影数据。在本领域中完成这种处理的方法是公知的。然而,当将所测量的覆盖成像***的视场的一半的投影数据与从其它的台架的视图中产生的数据相结合时,所得到的投影数据并不能与投影数据的中心附近相匹配。如果没有降低或消除这些不匹配的话,它将在所重构的图像中产生不希望的假象。
当前应用减小由在视场中的投影数据的不连续性引起的假象的一种技术是利用加权函数来平滑在过渡区中的数据的不连续性。这种技术要求检测器具有额外的检测器元件,这些额外的检测器元件延伸通过成像***的旋转中心在检测器上的投影。由于台架绕患者旋转360°,在两方向上稍稍移动延伸通过旋转中心在检测器上的投影的检测器面板区称为过渡区。实际的数据是通过检测器在一半过渡区中测量的,以及从台架的可替换视图中产生第二半过渡区中的数据。对在过渡区内所测量的和所综合的数据乘以用于平滑不连续性的加权系数。通常,较大的过渡区产生更好的图像质量,但是由于这种***的结构的视场稍小于比应用半检测器移动结构所产生的视场,所以它还导致了更高的***成本。可替换的是,需要更多的延伸过旋转中心在检测器上的投影的检测器元件,这造成了更高的***成本。
需要改善对所测量的数据和在过渡区内所产生的数据的完整性,以能够实现或要求检测器阵列的完整视场,以便单个检测器元件数可以减少。
在应用半检测器移动结构的多行或体积型CT***中,在成像***的一半视场中测量投影数据,同时另一半投影数据必需从反向的射线中产生或合成(即,从台架的角度位置中采集射线,该射线提供了所成像的材料的的线性衰减系数的类似的线积分数据)。不幸的是,如果检测器是它的最初的宽度的两倍并且没有偏移,则在CT台架的其它的投影角上测量的投影数据与已经测量的射线方向不具有相同的方向。这是因为从源发射的X-射线具有径向角度(即,射线在轴向平面之中的方向)和仰角(即,射线在轴向平面之外的方向)。因此,人们需要一种VCT***,这种VCT***应用在半检测器移动结构中的面积检测器并实现其优点,由此克服前述的困难。
本发明提供一种获得对象的投影数据的计算机X-射线断层成像(CT)***,该CT***部分包括X-射线源和检测器。检测器相对于与CT***的旋转中心在检测器上的投影相对应的中心位置移动它的宽度的一半。依据本发明的方法,对于每个投影视图,选择表示连接源和检测器元件的射线穿过对象的线性衰减的线积分的所处理的检测器元件值Va;这个信号的位置表示来自最靠近成像***的旋转中心在检测器上的投影的位置的检测器元件的信号。然后,对于所选择的检测器元件,检测器元件值Vb是从在相同方向所重构的数据的正向投影或相反的方向中估计的。然后选择能够消除Va和Vb之差的平滑函数。然后应用该平滑函数来消除a和Vb之差。
附图1所示为本发明的CT***方块图。
附图2所示为依据本发明的方法应用的检测器偏移。
附图3所示为依据优选的实施例说明本发明的方法的方块图。
在描述本发明的方法和装置之前,参考附图1总体上讨论一下本发明的VCT***。附图1所示为适合于实施本发明的方法和装置的体积型CT扫描***的方块图。虽然可以理解的是本发明并不限于对任何特定对象的成像,但仍然结合体积型CT扫描***在重构患者的解剖学特征的图像的应用中进行讨论。正如本领域的熟练技术人员将理解的是,本发明还可以用于工业过程。此外,本发明并不限于医用CT设备,而是包括工业***,在这些工业***中X-射线源和检测器结构都保持固定同时在扫描时间中绕对象旋转。
在体积型CT扫描***中,台架绕对象比如病人旋转,并采集投影数据。计算机1控制体积型CT扫描***的运行。当在此称台架的旋转时,该术语是指X-射线管2的旋转和/或检测器3的旋转,可取的是该检测器3是一种较高分辨率的面积检测器。台架包括X-射线管2和面积检测器3。控制器4A和4B受体积型CT扫描***计算机1控制并分别连接到X-射线管2和检测器3。控制器4A和4B使适当的旋转运动传递到X-射线管2和/或检测器3。并不是每个控制器都需要。可以使用单一控制器部件使台架旋转。还应该指出的是,为实施本发明的方法计算机1控制图像扫描时间、图像分辨率和/或轴向覆盖区的变化。
当对检测器3进行采集时计算机1通过给数据采集***6发指令并控制台架的速度来控制数据采集过程。此外,计算机1指令数据采集***6构造通过面积检测器3获得的射线照相的分辨率,由此能够改变***的分辨率。如图所示数据采集***6包括读出电子***。
面积检测器3包括检测器元件阵列(未示)。每个检测器元件测量与其相关的强度值,该强度值与辐照到检测器元件上的X-射线能量的大小相关。当本发明的装置和方法并入到体积型CT扫描***中,产生了一种新的体积型CT扫描***。因此,本发明还提供一种新颖的体积型CT扫描***。
还应该指出的是本发明并不限于任何特定的计算机来执行本发明的数据采集和处理的任务。如这里所使用的术语“计算机”是指任何能够执行计算并需要完成本发明的任务的机器。因此,用于实现本发明的控制算法10的计算机可以是能够执行所需的任务的任何计算机。
关于本发明,一种通过实施变型的数据平滑方案的方法消除了需要应用附加的检测器元件来覆盖过渡区的需要。此外,如果已经应用覆盖整个视场的较大的检测器阵列来采集数据,则一种变型的方法应用迭代算法来估计已经测量的投影数据。该跌代算法包括如下的步骤:重构CT数据、估计在所测量的投影数据和通过正向投影程序产生的数据之间的误差、消除和/或减小该误差以及重复该过程直到满足停止准则。该停止准则可以是对测量在所重构的数据在两种所跌代的方法之间的差的度量的阈值。这种度量的一个实例是在本领域中十分公知的最小二乘方准则。正向投影程序是一种通过计算沿着将X-射线源的位置连接到检测器的射线在对象或患者的CT重构中的线性衰减系数的线积分来模拟数据实验采集的方法。
由于在过渡区中的误差可以以类似的方式处理,下文将在相同的条件中讨论这两种方法。
本发明的一种技术通过正向投影从前面的重复步骤中所获得的重构的数据或通过对从一组反向射线中获得的冗余检测器数据进行插值来形成X-射线投影数组{Pa},该反向射线在{Pa}的相反的方向中形成了另一投影数据组。如前面所提到的,这种正向投影技术是这样的一种方法,即射线从假想的X-射线源发射;这些射线对着每个检测器元件穿过所重构的体积。沿着该射线,对所重构的值沿着射线的线性衰减值求和并表示为线性衰减系数的线积分。在本领域中十分公知的几种技术都可以用于实施正向投影处理。
正向投影所重构的数据的跌代技术(表示为IFPT)通常适合于产生与较大的锥形角相对应的投影数据(即,当用于VGT***中时),而对冗余的投影数据进行插值的技术(表示为RPDT)更适合于对更靠近中平面(即,零度锥角的平面)的投影数据。与扇形角类似,锥形角是指在与扇形角方向正交的方向上从X-射线源发射的X-射线的角度范围。在应用IFPT或RPDT所获得的估计检测器值和原始值(如果实际已经测量了该数据则将获得该值)之间的差可能使CT***的旋转中心在检测器上的投影的图像产生畸变。为简洁起见,CT***的旋转中心在检测器上的投影的位置称为中心射线位置。
RPDT技术对于多行检测器或面积检测器的中央的行有效,而IFPT在对于RPDT方法无效的锥角上的多行或面积检测器中的行都有效。
为降低在所重构的图像中的畸变或假象,已经研究出了一种应用平滑函数的方法,参考附图3该平滑函数表述如下:
1.对于每个投影视图,选择最接近在检测器上的中心射线位置21的已知的检测器元件,在此以后称为Va。
2.对于相同的检测器元件,得出一种估计值22(即,对来自替换的视图中的投影数据的插值(RPDT)或对所重构的数据进行正向投影(IFPT)),在此以后称为Vb。
3.产生适当的平滑函数23。
4.通过在Va和Vb之间的差值的一半加权平滑函数,并根据需要将该函数加到所测量的数据和合成的数据中/从所测量的数据和合成的数据中减去该函数。平滑函数降低在中间射线位置上的在Va和Vb之间的差并逐渐平滑在成像***24的视场的中心附近的区域中的这种差。
这可以从下面的讨论中得出,使d=Va-Vb,其中d是在视场的中心的投影数据中的不连续量。作为实例,可以用来逐步平滑该差的可能的平滑函数是一种如下定义的指数函数:
V=0.5de-ax 0
(等式1)
这里X0是距离与检测器元件值Va相对应的检测器位置的距离的绝对值,a是控制与平滑函数相关的曲线的斜率的系数。将指数函数加到位于中心射线位置(对应于成像***中的旋转中心在检测器上的投影的检测器位置)的一侧的投影值中(从位于中心射线位置(对应于成像***中的旋转中心在检测器上的投影的检测器位置)的一侧的投影值中减去该指数函数)以降低(提高)所估计的值,以及从在中心射线位置的替换侧的投影值中减去(加到在中心射线位置的替换侧的投影值中)以降低(增加)更高(更低)的所测量的值。换句话说,它提供了一种方法,当将真实投影数据和估计投影数据结合在一起时这种方法降低了在中心射线的投影数据的不一致性,因此该方法在数据中提供了一种光滑的过渡区,该过渡区降低或消除了假象25。
应用上面所描述的方法,可以消除检测器在中心射线位置上的重叠量,由此提供一种成本效率合算的方案设计。
目前,在已有的技术中如何在半检测器移动结构的体积型CT(VCT)***中应用面积检测器并不清楚。前文已经描述了VCT***和通过产生过渡区消除在成像***的视场中的投影数据的不连续性的不同的公知的技术,下文将描述本发明的另一方面。因此,将概括上述讨论以得到一般性的实施例。
应用变量fθ和fθ’来分别表示在源角度θ处所得到的正向和反向射线(相同的角度方位但以相反的方向通过的射线)的信号强度,这里
fθ(n)=0 对于N2<n<N (等式2)
fθ’(n)=0 对于1<n<N2 (等式3)
应用变量n表示在检测器阵列中的位置。例如,如果检测器阵列具有1000个检测器,1<n<1000。N是水平检测器元件的数量,N2是中心射线在检测器上的位置。理想的fθ(N2)应该精确地等于fθ’(N2),因为两者都穿过对象的相同的部分。但是由于下面的原因这永远不可能:
(a)每个射线的实际形状是从源发出并在检测器终止的空的四面体。没有完全相同的射线通过对象的相同的部分,除非对象是完全均匀的并圆形对称。
(b)在扫描周期中对象/患者的运动都可能在每个正向/反向射线对中引入附加的误差。
(c)迄今为止还没有研制出完美的有效的插值方案。这就是说插值过程可能引入误差。
假设d(θ)是在角度源位置θ的fθ(N2)和fθ’(N2)之差,如果d(θ)是完全随机的,则可能由所重构的图像感应的误差会被与其它的CT随机误差相连的量子噪声所掩盖。然而,如果误差是某种***的误差,它将在所重构的图像中引入明显的假象。由于这个原因,必需在过渡区应用平滑处理。换句话说,研究出一种平滑函数来使在检测器元件N2所表示的中心射线位置周围的fθ和fθ’的幅值误差更小。
应用W和W’分别表示fθ和fθ’的的平滑函数。当对W和W’求导时,必需考虑一定的规则,对于本领域的熟练技术人员来说这些规则都是可以理解的。此外,正如本领域的熟练技术人员将会理解的是,除了在此所特别说明的函数外,不同的平滑函数可以用于此目的。例如
W(n)+W’(n)=1 对于所有的n (等式4)
这里Δn设定W和W’的平滑范围,δ是微分算子。应该指出的是,对于这种类型的应用常规的平滑函数还通常称作尾翼函数(feathering function)。应用W和W’在正向和反向射线之间寻找过渡区。应该指出的是,为使平滑函数有效,Δn必需是比零大的整数。实际上,Δn越大,平滑函数的效果越好。然而,Δn增加得太多可能要求附加的检测器元件来延伸到在中心射线位置的检测器元件N2之外。因此,应该选择Δn足够大但是又不大到要求对过渡区增加附加的检测器元件。这样,对fθ和fθ’作如下的限制:
fθ(n)=0 对于N2+Δn<n<N (等式6)
fθ’(n)=0 对于1<n<N2-Δn (等式7)
在滤波背向-投影方法中所应用的实际检测器信号是Wfθ和W’fθ’。还应该指出的是,由于更宽的过渡区易于消除在正向和反向射线之间的许多不匹配的误差,所以对于每个半视场(FOV)的投影数据不需要叠加反向射线。换句话说,每半个FOV数据(加上附加的Δn检测器值)填充以零以获得长度为N的检测器数据,在此之后进行常规的滤波背向投影程序。不包含插值过程。
起因是当开始增加在CT***中即在面积检测器中的检测器的行数时附加检测器元件(Δn乘以行数)变得更大。因此,需要通过设计一种使Δn最小的方法和装置来改善常规的方法。
在本发明中这种方法可以将Δn减小到1,而计算机模拟表明常规的平滑方法要求Δn大约20才能实现相当的假象水平。在VCT应用中这种优点更有意义,在VCT中在面积检测器的过渡区中所需的检测器元件的数目可能比在线性阵列中所需的元件还多三次幂的数量级。
在正向和反向射线中可能有***误差,通过插值或正向投影获得在反向射线中***误差。通过求fθ(N2)和fθ’(N2)之差测量幅值误差。使d(θ)=fθ(N2)-fθ’(N2)。将d(θ)看作fθ(N2)和fθ’(N2)之间的幅值误差,在此θ是X-射线源的角度位置。
我们的方法是应用例如在等式1中所述的指数函数来消除幅值误差,在等式1中a是控制指数函数的平滑性的控制系数。将正向和反向射线函数fθ(n)和fθ’(n)分别转换成投影数据gθ(n),1<n<N:
gθ(n)=fθ(n)-pθ(N2-n) (等式8)
gθ’(n)=fθ’(n)+pθ(n-N2) (等式9)
使等式8和等式9的gθ(N2)=gθ’(N2)。分量pθ(n)表示平滑函数;这种函数的一种实例描述在等式1中。
对于每个投影图像,实施下面的过程:
1.获得原始的半FOV投影数据,称其为fθ(n),并依据等式2补零。
2.获得反向射线fθ’(n)的数组,并依据等式3进行补零。
3.依据等式8和等式9基于幅值误差d(θ)(这里d(θ)=fθ(N2)-fθ’(N2))应用平滑函数。
4.产生gθ(n)以形成N-检测器数组。
5.对hθ(n)应用常规的滤波背向投影。
6.对所有的投影角度重复步骤1至5。
作为本发明的一种方法,对于在那些可以通过对在不同的角度之外的其它投影数据进行插值来获得反向射线的任何CT扫描器来说上面的过程都是有效的。在如下的2D扇形束中比较理想:另一半FOV数据总是可以大致地从冗余的扇形束投影数据中计算出。当将这种方法扩展到3D的VCT中时,当使用圆形轨道时仅在中平面(仰角为零的X-射线平面)上能够精确地插值。
我们的结果表明当将相同的方法应用到应用圆形轨道的VCT中时,对于在±1.5度的锥形角度(或者在此称为仰角)这种方法比常规的平滑方法(应用超过中心射线的位置的20个附加的检测器)更优越。如果应用完整的检测器,由于反向射线的角度方位与已经测量的数据差别极大,所以对于较大的锥形角度这种结果就不正确。为理解这种概念,考虑X-射线穿过CT扫描器的旋转轴线并偏离检测器的中平面;因此仰角不为零。现在,转动射线源和检测器点180°。无疑,该射线并不穿过患者模型的相同的区域。为补救这种情况,可以应用迭代的方法来提高图像的质量。该过程如下:
1.依据上述的步骤1至6获得初始的3D图像。从上述的讨论中可知,在所重构的图像中极有可能存在假象。
2.在该方法的第二次迭代中应用正向投影方法获得每个半FOV投影数据组的“反向射线”并依据上文所述的依据步骤3至步骤6将它们与从来自另一半FOV投影中所测量的投影数据相结合形成一个完整的投影数据组。在投影数据组中逐行地实施平滑函数。
3.继续到本方法的步骤2直到重构收敛,即图像的质量不再改变/提高。
应该注意的是,结合一定的实施例已经讨论了本发明。然而,本发明并不限于这些实施例。例如,这里所讨论的平滑函数并不意味着排除可以用于该目的的该类型的其它函数。此外,这些方案并非限于一种扫描方案,即它们还可以应用到轴向扫描(在扫描周期中患者的工作台并不移动)和螺旋扫描方案中。在本领域的熟练技术人员会理解这些方式,在这些方式中应用这些概念并类推以实现对特定领域的应用很有用的其它面积检测器扫描方案。
Claims (9)
1.一种获得对象的投影数据的体积型计算机X-射线断层成像即VCT***,该VCT***包括:
X-射线源,该X-射线源对对象投影X-射线;
检测器,该检测器相对于与VCT***的旋转中心在检测器上的投影相对应的中心位置偏移了其宽度的一半,该检测器接收从X-射线源投影的X-射线并产生响应辐照在其上的X-射线的电信号;
从检测器读取电信号并将该电信号转换为数字信号的数据采集***;以及
能够执行重构算法的计算机,该计算机从数据采集***部分接收电信号,其中当计算机执行重构算法以处理所说的数字信号时,计算机重构图像。
2.权利要求1所述的VCT***,其中该检测器包括多个检测器元件,且其中,对于多个投影视图中的每一个,通过计算机选择与检测器元件相关的检测器元件值Va,该检测器元件值Va表示连接该X-射线源和具有该检测器元件值的检测器元件的射线穿过的对象的线性衰减系数的线积分,以及其中该计算机估计与检测器元件值Va相对的检测器元件的值Vb,其中对于每个投影视图,Va对应于所测量的投影数据,而Vb对应于估计的投影数据,其中一旦获得每个投影视图的值Va和Vb,计算机应用平滑函数来消除在Va和Vb之间的差。
3.权利要求2所述的VCT***,其中该检测器是具有相对较高分辨率的面积检测器。
4.权利要求3所述的VCT***,其中在应用平滑函数之前,对该平滑函数进行加权,其中将经加权的平滑函数应用到所测量的和所估计的投影数据中以减小在检测器的中心射线位置上的值Va和Vb之间的差,并消除在VCT***的视场的中心附近的区域中的值Va和Vb之间的差。
5.权利要求4所述的VCT***,其中通过将Va和Vb之差的一半加入到相关的所测量的和所估计的投影数据中来对平滑函数进行加权。
6.权利要求4所述的VCT***,其中通过从相关的所测量的和所估计的投影数据中减去Va和Vb之差的一半来对平滑函数进行加权。
7.一种应用计算X-射线机断层成像(CT)***获得对象的投影数据的方法,该方法包括如下的步骤:
从X-射线源向对象投影X-射线;
在检测器的投影视图上接收CT***在检测器上的投影的X-射线,该检测器相对于与VCT***的旋转中心在检测器上的投影相对应的中心位置偏移了其宽度的一半,该检测器包括许多检测器元件,该检测器元件产生响应在辐射在其上的X-射线的电信号;
数字化该电信号;
对于每个投影视图,选择最靠近CT***的中心射线位置的检测器元件的检测器元件值Va;
对于所选择的检测器元件值,经过从反方向采集的投影数据的插值或从在相同的方向中重构的数据的正向投影中估计与检测器元件值Va相关的检测器元件的值Vb;
选择能够消除在值Va和Vb之间的差的平滑函数,投影视图的值Va对应于所测量的投影数据,而投影视图的值Vb对应于估计的投影数据;
对平滑函数进行加权;以及
当结合所测量的投影数据和估计的投影数据时应用加权的平滑函数来减小差值以产生平滑的过渡区。
8.如权利要求7所述的方法,其中通过将Va和Vb之差的一半加入到相关的所测量的和所估计的投影数据中来对平滑函数进行加权。
9.如权利要求7所述的方法,其中通过从相关的所测量的和所估计的投影数据中减去Va和Vb之差的一半来对平滑函数进行加权。
Applications Claiming Priority (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12939899P | 1999-04-15 | 1999-04-15 | |
US60/129,398 | 1999-04-15 | ||
US60/129398 | 1999-04-15 | ||
US16650099P | 1999-11-19 | 1999-11-19 | |
US60/166500 | 1999-11-19 | ||
US60/166,500 | 1999-11-19 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1300201A CN1300201A (zh) | 2001-06-20 |
CN1210000C true CN1210000C (zh) | 2005-07-13 |
Family
ID=26827535
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 00800601 Expired - Fee Related CN1210000C (zh) | 1999-04-15 | 2000-04-14 | 在利用仅覆盖一半视场的减小尺寸的检测器的计算机x-射线断层成像***中应用的方法和装置 |
Country Status (5)
Country | Link |
---|---|
JP (2) | JP2002541896A (zh) |
CN (1) | CN1210000C (zh) |
AU (1) | AU4352400A (zh) |
DE (1) | DE10081367T1 (zh) |
WO (1) | WO2000062674A1 (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2820045A1 (fr) | 2001-01-29 | 2002-08-02 | Chabunda Christophe Mwanza | Dispositif bilina-imatron de stereoradiotherapie par acceleration et collision des particules des 2 faisceaux d'irriadiation synchrones associe a un dispositif de verification instantannee de delivrance des radiations |
US7099428B2 (en) | 2002-06-25 | 2006-08-29 | The Regents Of The University Of Michigan | High spatial resolution X-ray computed tomography (CT) system |
GB2422759B (en) * | 2004-08-05 | 2008-07-16 | Elekta Ab | Rotatable X-ray scan apparatus with cone beam offset |
JP4820561B2 (ja) * | 2005-03-14 | 2011-11-24 | 株式会社東芝 | 核医学診断装置 |
WO2007020318A2 (en) * | 2005-08-17 | 2007-02-22 | Palodex Group Oy | X-ray imaging apparatus and x-ray imaging method for eccentric ct scanning |
CN100435733C (zh) * | 2005-12-31 | 2008-11-26 | 清华大学 | X-ct扫描*** |
US7945012B2 (en) | 2006-08-17 | 2011-05-17 | Koninklijke Philips Electronics N.V. | Computed tomography image acquisition |
WO2008092009A2 (en) | 2007-01-24 | 2008-07-31 | Imaging Sciences International Llc | Adjustable scanner |
DE602008005584D1 (de) | 2007-06-11 | 2011-04-28 | Morgan Electro Ceramics Bv | Bildgebungssystem und bildgebungsverfahren zur abbildung eines interessensbereiches |
CN102124496B (zh) * | 2008-05-09 | 2013-05-29 | 皇家飞利浦电子股份有限公司 | 用于产生活动对象的图像的设备 |
CN102325499B (zh) | 2009-01-21 | 2014-07-16 | 皇家飞利浦电子股份有限公司 | 用于大视场成像及运动伪影的探测和补偿的方法和装置 |
CN104545976B (zh) * | 2014-12-30 | 2017-04-19 | 上海优益基医疗器械有限公司 | 计算机体层摄影方法和装置 |
CN110121300B (zh) * | 2016-12-21 | 2023-06-16 | 皇家飞利浦有限公司 | 缩放式射线照相重建 |
CN107714072B (zh) * | 2017-11-20 | 2019-09-20 | 中国科学院高能物理研究所 | 缺失数据的补偿方法、计算机断层扫描成像方法及*** |
US10426424B2 (en) | 2017-11-21 | 2019-10-01 | General Electric Company | System and method for generating and performing imaging protocol simulations |
CN117257340A (zh) * | 2023-11-21 | 2023-12-22 | 北京朗视仪器股份有限公司 | 一种异形探测器、医疗影像设备以及图像补全方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0824674B2 (ja) * | 1986-03-31 | 1996-03-13 | 株式会社東芝 | X線ct装置 |
US5265142A (en) * | 1992-05-08 | 1993-11-23 | General Electric Company | Image reconstruction technique for a computer tomography system |
JP3637074B2 (ja) * | 1992-12-15 | 2005-04-06 | 株式会社東芝 | ヘリカルスキャン方式のコンピュータ断層撮影装置 |
US5400379A (en) * | 1994-02-25 | 1995-03-21 | General Electric Company | Multi-slice x-ray CT using a detector mask |
US5848117A (en) * | 1996-11-27 | 1998-12-08 | Analogic Corporation | Apparatus and method for computed tomography scanning using halfscan reconstruction with asymmetric detector system |
JPH10290798A (ja) * | 1997-04-17 | 1998-11-04 | Ge Yokogawa Medical Syst Ltd | 投影データ測定方法および装置並びにx線ct装置 |
DE19721535C2 (de) * | 1997-05-22 | 2001-09-06 | Siemens Ag | Röntgen-Computertomograph zur Erzeugung von Röntgenschattenbildern |
-
2000
- 2000-04-14 JP JP2000611814A patent/JP2002541896A/ja active Pending
- 2000-04-14 CN CN 00800601 patent/CN1210000C/zh not_active Expired - Fee Related
- 2000-04-14 DE DE10081367T patent/DE10081367T1/de not_active Withdrawn
- 2000-04-14 AU AU43524/00A patent/AU4352400A/en not_active Abandoned
- 2000-04-14 WO PCT/US2000/010170 patent/WO2000062674A1/en active Application Filing
-
2010
- 2010-11-15 JP JP2010254398A patent/JP5194095B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
AU4352400A (en) | 2000-11-02 |
JP5194095B2 (ja) | 2013-05-08 |
JP2011031070A (ja) | 2011-02-17 |
CN1300201A (zh) | 2001-06-20 |
JP2002541896A (ja) | 2002-12-10 |
DE10081367T1 (de) | 2001-06-28 |
WO2000062674A1 (en) | 2000-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1210000C (zh) | 在利用仅覆盖一半视场的减小尺寸的检测器的计算机x-射线断层成像***中应用的方法和装置 | |
CN105596022B (zh) | 锥束ct动态成像的方法和设备 | |
CN1292712C (zh) | 多模态成像方法和装置 | |
US7062009B2 (en) | Helical interpolation for an asymmetric multi-slice scanner | |
CN1182544C (zh) | 数字化无胶片x-射线投影成象***和方法 | |
JP6139092B2 (ja) | X線ct装置およびシステム | |
CN1217625C (zh) | 用于具有任意节距的多切面ct成像的部分扫描加权 | |
CN101036164A (zh) | 采用与体素相关插值的图像重建 | |
WO1995022115A1 (en) | Reconstruction of images from cone beam data | |
CN1502309A (zh) | 便于减少伪影的方法和装置 | |
WO1999001065A1 (en) | Iterative cone-beam ct reconstruction | |
JP2002320612A (ja) | 高品質画像再構成のためのコンピュータ断層撮影(ct)重み付け法 | |
US7529335B2 (en) | Voxel-driven spiral reconstruction for cone-beam computer tomography | |
JP2010259653A (ja) | X線コンピュータ断層撮影装置 | |
CN1605322A (zh) | 使扫描图像中的模糊最小化的方法和设备 | |
US8798350B2 (en) | Method and system for reconstruction algorithm in cone beam CT with differentiation in one direction on detector | |
CN101056584A (zh) | 用于检查周期性运动的对象的ct方法 | |
EP1927018B1 (en) | Improved rebinning for computed tomography imaging | |
CN1418599A (zh) | 荧光检查计算机断层方法 | |
JP4121198B2 (ja) | コンピュータ断層撮影システム用の検出器 | |
IL130718A (en) | Methods and apparatus for helical multi-frame image reconstruction in a computed tomography fluoro system | |
JP2002034970A (ja) | マルチ・スライスct走査の螺旋再構成の方法及び装置 | |
US20090116721A1 (en) | Method and system for multi-detector-row computed tomography utilizing projection data up-sampling with shift | |
US7020317B2 (en) | Computed tomography apparatus | |
JP2003225230A (ja) | コンピュータ断層撮影装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20050713 Termination date: 20150414 |
|
EXPY | Termination of patent right or utility model |