CN102138804A - 在计算开支最小化条件下降低焦外辐射效应的ct数据处理 - Google Patents
在计算开支最小化条件下降低焦外辐射效应的ct数据处理 Download PDFInfo
- Publication number
- CN102138804A CN102138804A CN201010580280XA CN201010580280A CN102138804A CN 102138804 A CN102138804 A CN 102138804A CN 201010580280X A CN201010580280X A CN 201010580280XA CN 201010580280 A CN201010580280 A CN 201010580280A CN 102138804 A CN102138804 A CN 102138804A
- Authority
- CN
- China
- Prior art keywords
- measurement data
- att
- detecting element
- efs
- data
- 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.)
- Granted
Links
Images
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/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
- A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
- A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
Abstract
本发明涉及一种用于从测量数据(p)中重建检查对象的图像数据(f)的方法,其中,所述测量数据(p)是在计算机断层扫描系统(C1)的发射焦点辐射和焦外辐射的辐射源(C2、C4)与检查对象之间相对旋转运动时由检测器(C3、C5)采集的。检查投影的测量数据(p)的彼此的差别,以及进行测量数据校正,以降低焦外辐射的效应。在此,仅关于所选取的测量数据(p)进行所述校正,而测量数据(p)的选取则取决于检查结果。在使用得到校正的测量数据(p)情况下进行图像数据(f)的重建。
Description
技术领域
本发明涉及一种用于从测量数据中重建检查对象的图像数据的方法,其中,所述测量数据是在计算机断层扫描系统的发射焦点辐射和焦外辐射的辐射源与检查对象之间相对旋转运动时由检测器检测的。
背景技术
利用CT系统扫描检查对象的方法众所周知。在此方面,例如使用圆扫描、具有进给的顺序圆扫描或者螺旋扫描。在这些扫描中,利用至少一个X射线源和至少一个相对的检测器从不同的拍摄角度拍摄检查对象的吸收数据并将这样收集的吸收数据或投影借助相应的重建方法计算出检查对象的剖面图像。
为从计算机断层扫描仪(CT机)的X射线-CT-数据组中,也就是从所检测的投影中重建计算机断层扫描的图像,目前作为标准方法使用所谓的滤波反投影(Filtered Back Projection;FBP)。
重建的CT图像的对比度和清晰度此外取决于焦点的尺寸,也就是说,发射X辐射的X射线管阳极的那个范围。通常X射线管既发射焦点射线,也发射焦外辐射,也就是(可能明显)超出固有的有效焦点以外形成的辐射。焦外辐射扩大了X射线管的辐射面积并因此使图像的对比度和清晰度变差。
发明内容
本发明的目的在于,提供一种用于重建CT图像的方法,其中,需要考虑,X射线管既发射焦点辐射,也发射焦外辐射。在这种情况下特别是注意降低计算开支。此外提供一种相应的控制和计算单元、CT系统、计算机程序和计算机程序制品。
依据本发明的方法涉及从测量数据中重建检查对象的图像数据。在此,所述测量数据是在计算机断层扫描系统的发射焦点辐射和焦外辐射的辐射源与检查对象之间相对旋转运动时由检测器采集的。检查投影的测量数据的彼此的差别。此外进行测量数据的校正,以降低焦外辐射的效应。在这种情况下,仅关于所选取的测量数据进行所述校正,而测量数据的选取则取决于检查结果。在使用得到校正的测量数据情况下进行图像数据的重建。
在数据采集后存在大量投影的测量数据。从这些测量数据中要重建检查对象的图像。为此通常使用至少180°投影角度范围的测量数据。在图像重建之前进行EFS校正(EFS:焦外辐射)。为此处理一个、多个或者最好所有投影的测量数据。最好在求对数之前对按照扇形或者锥形射线几何形状的强度数据进行EFS校正。
关于投影,不是校正所有测量数据,而是仅校正所选取的测量数据。根据对测量数据的彼此差别的检查来选取需要校正的这些测量数据。根据这些差别来校正投影的某些测量数据,而不校正其他测量数据。所要校正的测量数据的这种选取在典型的临床数据组情况下明显减少图像重建总共所需的计算时间。
为对所选取的测量数据进行EFS校正使用本身公知的方法。此方面的例子是利用专用的卷积核(Faltungskern)对强度数据进行卷积,该卷积核精确或者良好近似地产生EFS份额的“去卷积”。EFS校正改变测量数据,使得它们类似于没有或者具有少量EFS所获得的那些测量数据。
最后对于所观察的投影存在新的测量数据组,其中对某些测量数据进行了EFS校正,而对其他测量数据没有进行EFS校正。这种新的测量数据组用于重建检查对象的图像。该图像与在使用原始测量数据的情况下获得的图像相比得到改进:对比度和清晰度得到提高,因为基于EFS的伪影不存在或者很少存在。
在本发明的进一步构成中,仅对与其他测量数据具有最小差别的测量数据进行校正。为此可以确定阈值,超过该阈值则说明,需要对相应的测量数据进行EFS校正。所述其他测量数据(关于其确定差别,也就是各个所观察的测量数据与其进行比较)可以按照不同方式来确定。
例如,检测器可以具有至少一个包括多个检测元件的检测行,在检查时将一个检测行的检测元件的测量数据与相同检测行的处于确定距离上的另一个检测元件的测量数据进行比较。这种距离可以意味着使用直接相邻的检测元件;但距离最好更大。对所有检测元件最好使用相同的距离。
例如,检测器可以具有至少一个包括多个检测元件的检测行,在检查时将一个检测行的检测元件的测量数据与相同检测行的处于确定距离上的另外两个检测元件的测量数据进行比较。具有优点的是,另外两个检测元件与所观察的检测元件的两侧对称设置。但两侧上的距离也可以不同。在这里同样适用的是,可以使用直接相邻的检测元件;但距离最好更大。对所有检测元件最好使用相同的距离。此外,另外两个检测元件可以处于所观察的检测元件的相同侧上。
例如,检测器可以具有至少一个包括多个检测元件的检测行,并在检查时将一个检测行的检测元件的测量数据与相同检测行的处于确定距离上的另外四个检测元件的测量数据进行比较。具有优点的是,另外四个检测元件与所观察的检测元件的两侧对称设置,也就是说,四个检测元件的两个设置在检测元件的左侧并且四个检测元件的两个设置在右侧,其中,到左侧的距离与到右侧的距离相同。但检测元件两侧上的距离也可以不同。在这里也适用的是,可以使用直接相邻的检测元件;但距离最好更大。对所有检测元件最好使用相同的距离。此外,另外四个检测元件可以处于所观察的检测元件的相同侧上。
最好对检测行的所有检测元件进行比较。如果一个检测元件设置在行的边缘上,那么不能保持确定的距离,因为不存在这种距离的检测元件,所以对该检测元件要么不进行比较,要么选择其他距离。如果检测器具有多行,那么最好对每个检测行的所有检测元件进行比较。
关于为比较而使用的唯一的检测行内部的检测元件的确定所介绍的方法也可以关于在不同检测行内的检测元件被采用。
在本发明的构成中,测量数据校正与一个检测元件相关地进行,方法是,通过卷积将该检测元件的测量数据与相同检测行的一定数量的相邻检测元件的测量数据相关联。在这种相互联系中特别有益的是,其他检测元件与所观察的检测元件具有的确定距离相当于相邻检测元件数量的四分之一或者至少大约四分之一。
依据本发明的进一步构成,这样进行选取,使测量数据的校正显而易见仅在测量数据对比度大的区域内起作用。为此目的,注意力特别是转向投影的测量值之间的大的差别上。这样可以有针对性地使受到EFS特别不利影响的那些图像区得到改进。
依据本发明的优选构成,为选取测量数据,使用显示需要进行校正的那些测量数据的函数。为此例如适合使用数字函数。具有优点的是,对于在函数从需要进行校正的测量数据向不需要进行校正的测量数据过渡的情况下的测量数据,计算已校正和未校正测量数据的混合。这一点相当于叠化测量数据并避免函数过渡上的伪影。
依据本发明的控制和计算单元用于从CT系统的测量数据中重建检查对象的图像数据。该控制和计算单元包括用于储存程序代码的程序存储器,其中,程序存储器内-必要时除了别的之外-存在适用于实施上述类型方法的程序代码。依据本发明的CT系统包括这种控制和计算单元。此外,该系统可以包括例如用于采集测量数据所需的其他组成部分。
依据本发明的计算机程序具有程序代码装置,当在计算机上执行计算机程序时,所述程序代码装置适用于实施上述类型的方法。
依据本发明的计算机程序制品包括储存在计算机可读数据载体上的程序代码装置,其适用于在计算机上执行计算机程序的情况下实施上述类型的方法。
附图说明
下面借助实施例对本发明进行详细说明。其中:
图1示出具有图像重建组成部分的计算机断层扫描系统的实施例的第一示意图;
图2示出具有图像重建组成部分的计算机断层扫描系统的实施例的第二示意图;
图3示出具有焦外辐射的X射线管的示意组成部分;
图4示出具有焦点辐射和焦外辐射的CT数据采集;
图5示出不同检测元件的CT测量数据(衰减数据)的曲线和具有从中计算的对比度值的曲线;
图6A示出具有图5平方对比度值的曲线和具有差值的曲线;
图6B示出用于显示活化区的曲线和具有叠化系数的曲线。
具体实施方式
图1首先示意示出具有图像重建装置C21的第一计算机断层扫描系统C1。处于机架外壳C6内的是这里未示出的封闭机架,其上设置第一X射线管C2及相对的检测器C3。可选的是,在这里所示的CT系统中,设置第二X射线管C4及相对的检测器C5,从而通过附加可供使用的辐射器/检测器组合可以达到更高的时间分辨率,或者在使用不同的X射线能量光谱时,辐射器/检测器系统内也可以进行“双能量(Dual-Energy)”检查。
CT系统C1还具有病床C8,检查时患者可以在病床上沿着也称为z轴的系统轴线C9移动到测量场内,其中,扫描本身可以作为没有患者进给的纯圆形扫描,仅在所关心的检查范围内进行。在这种情况下,X射线源C2或C4分别环绕患者旋转。在此,检测器C3或C5与X射线源C2或C4相对地平行地同时运转,以采集投影测量数据,这些投影测量数据然后被用于重建剖面图像。作为在各个扫描之间逐步地移动患者通过检查区的顺序扫描的替换,不言而喻,也存在螺旋扫描的可能性,其中在采用X射线辐射的循环扫描期间连续沿者系统轴线C9移动患者通过X射线管C2或C4与检测器C3或C5之间的检查区。通过患者沿轴线C9的运动和X射线源C2或C4的同时循环,在测量期间在X射线源C2或C4相对于患者的螺旋扫描时形成螺旋轨迹。这种轨迹也可以通过如下来实现,即在患者不运动的情况下机架沿轴线C9移动。
CT系统C1通过具有在存储器中存在计算机程序代码Prg1-Prgn的控制和计算单元C10控制。从控制和计算单元C10可以通过控制接口24传递采集控制信号AS,以便依据确定的测量协议控制CT系统C1。
由检测器C3或C5采集的投影测量数据p(下面也称为原始数据)通过原始数据接口C23传送到控制和计算单元C10。然后这些原始数据p(在需要时在适当的预处理后)在图像重建组成部分C21内被进一步处理。图像重建组成部分C21在该实施例中在控制和计算单元C10内以处理器上的软件方式,例如以一个或者多个计算机程序代码Prg1-Prgn的方式实现。由图像重建组成部分C21重建的图像数据f然后被存储在控制和计算单元C10的存储器C22内和/或者以常见的方式在控制和计算单元C10的屏幕上被输出。它们也可以通过图1中未示出的接口被输入到与计算机断层扫描系统C1连接的网络内,例如放射信息系统(RIS)并存储在那里接入的大容量存储器内或者作为图像输出。
控制和计算单元C10附加地还可以执行EKG函数,其中,导线C12用于在患者与控制和计算单元C10之间传导EKG势能。图1中所示的CT系统C1附加地还具有造影剂注射器C11,通过其可以附加地将造影剂注入患者的血液循环内,从而可以更好显示患者的血管,特别是跳动的心脏的心室。此外还可以实施灌注测量,对于灌注测量同样适用所提出方法。
图2示出一个C弧系统,其中与图1的CT系统不同,外壳C6携带C弧C7,其上一方面固定了X射线管C2、另一方面固定了相对的检测器C3。C弧C7为扫描同样环绕系统轴线C9回转,从而可以从大量的扫描角度进行扫描并可以从大量的投影角度中测定相应的投影数据p。图2的C弧系统C1与图1中的CT系统一样具有图1所述类型的控制和计算单元C10。
本发明可以在图1和2所示的两个系统中使用。此外,原则上也可以用于其他CT系统,例如具有形成完整一圈检测器的CT系统。
图3示出X射线管的示意图,正如可以在图1和2的系统中使用的那样。通过利用在阴极K与阳极A之间施加的电压来加速从热阴极K发出的电子e-,产生由X射线管发射的X射线辐射。在快速的电子e-进入阳极材料(例如钨)内时,形成X射线辐射。这种辐射主要相当于电子e-的减速辐射。
重建的图像的清晰度主要取决于X射线管阳极A上焦斑(Brennfleck)的尺寸。有效焦斑Fok,也就是发射大部分X射线辐射的阳极A的范围称为有效焦点。在诊断的X射线管情况下焦点尺寸通常在0.3mm到2mm之间。根据X射线管的结构的不同,在本身的有效焦点Fok的外面在数厘米的范围上发出X射线辐射,其因此使图像的对比度显示不清晰、变差。
这种寄生的X射线辐射称为焦外辐射,简称EFS(英语:EFR,extra focalradiation,也写做off-focus radiation)。EFS的形成做如下解释:一部分高速击中阳极A的电子e-要么弹性地从阳极A反射,要么它们在阳极A内引发重新离开阳极表面的次级电子。这些散射的初级或者次级电子e- 散射的能量与初级电子e-的能量相比减少约20%。通过阳极A的电场的吸引,电子e- 散射再次击中阳极A。由这些电子e- 散射产生的X射线辐射为焦外辐射。由于电子e- 散射前面的能量损耗,EFS平均比焦点的X射线辐射更柔。在此,散射电子e- 散射的击中部位一般情况下远离本来的焦点Fok。电子e- 散射扩大了发射区和由此的成像的辐射源,使有效焦点Fok扩展。这一点在图3中通过有效焦点Fok旁边的区域Δ表示。由X射线管发射的全部辐射的EFS比例根据X射线管的结构典型地最高约为10%。
如果没有将EFS全部析出,那么它成为扫描检查对象所使用的X射线辐射的组成部分。图4示意示出通过检查对象O的投影值的采集。检查对象O内部空间上的衰减分布或密度分布为f(x)。这一点适用于通过从所采集的投影值(也就是从测量数据)中重建来测定。然后f(x)可以作为灰度值图像显示。h(t)表示阳极上的发射分布;它包括焦点辐射和焦外辐射。发射分布h(t)即说明从阳极的哪个点发出多少X射线发射。
阳极的延伸-简化为一维-采用t表示。η D表示确定的检测器像素或元件。检测器具有一个或者多个检测行,其中,每行包括大量并排设置的检测元件。ζF(t)、ζF(t′)和ζF(t″)为阳极的位置t、t′和t″向检测器像素η D的X射线。X射线通过检查对象O的轨迹沿线性参数s分布。角度α是投影角度,该角度在CT系统的辐射源/接收器对环绕检查对象O旋转时改变。即可以通过α的数据来识别确定的投影。在每个投影中,由每个检测元件采集测量值。可以看出,无论是由焦点发出的射线ζF(t′),还是由焦点外面发出的射线ζF(t)和ζF(t″),均为检测器像素η D的测量结果提供份额。
EFS的存在使CT机的调制传输函数变差。借助调制传输函数可以将重建的图像内部的对比度与成像的对象内部的对比度进行比较。调制传输函数是本振频率的函数,其作为每个长度单位的线对(Linienpaar)的数量给出。通过可能与有效焦点距离相当远的EFS,调制传输函数特别是在本振频率小的情况下变差。这样使图像的对比度变差。特别是低对比度的对象组成部分在高对比度的对象组成部分的附近更差地成像。例如在拍摄头颅时会出现这种不希望的效应。在头颅边缘的附近,例如在靠近外头盖骨的大脑区内,EFS可能使对比度发生明显变化。
为消除EFS的消极影响,在重建检查对象的图像之前可以进行测量数据的EFS去卷积(英语:EFS deconvolution)。为此使用去卷积函数,其考虑哪些射线来自焦点和哪些射线是EFS。为此需要认识X射线源的辐射特性(Abstrahlcharakteristik)以及所检查的对象的(至少近似的)模型。通过EFS去卷积从测量数据中“计算出”EFS射线。
进行EFS去卷积,方法是,在使用去卷积函数的情况下,为每个投影重新计算(也就是校正)所有测量数据,也就是每个检测行的每个检测元件的测量数据。关于确定的检测元件,这一点意味着:去卷积函数具有确定的宽度,例如49;与此相应观察一系列的测量数据,这些测量数据由所观察的检测元件一侧上、中间的所观察的检测元件上的检测元件的24,和所观察的检测元件另一侧上的检测元件的24个测量数据组成。长度49的这一系列测量数据现在利用去卷积函数卷积。该计算的结果是各个检测元件的EFS校正测量值。对于各个投影的所有检测元件并且对于所有投影进行相应的计算。因此产生一个用于图像重建的得到修正的测量数据组。由此重建的图像消除了通过EFS引起的伪影。
但卷积运算的计算非常费时。与此相应,通过测量数据的去卷积明显提高计算CT图像的持续时间。因此通常仅对于所选取的图像进行EFS去卷积,即对于这种图像进行:其在对象边缘附近的低对比度显示对诊断是特别重要的。
为取得测量数据的快速EFS去卷积,按如下所述进行:
仅在如下图像区内应用EFS去卷积:在该图像区中,EFS由于伪影或质量损失而变得明显。上面已经详细介绍过,由于EFS导致的质量损失主要在对比度非常大的对象组成部分的附近出现。与此相应,EFS去卷积被限制在高对比度的区域上。这一点在图5中示意示出。
图5中横坐标上标出信道编号N。在这种情况下,这是一行检测元件的编号。图5的曲线ATT示出对数的衰减值。为获得这些数值,对各个测量值(也就是由检测元件接收的X射线强度)标准化、求对数和改变符号。曲线ATT因此相当于测量数据。导致图5的曲线ATT的拍摄是头颅扫描。右侧和左侧上的两个波峰来自头部外壳,在拍摄期间头颅由头部外壳保持。陡峭的上升和下降由如下的X射线引起:该X射线纵向地透射头皮和外头盖骨。圆形顶峰来自既穿过头盖骨也穿过脑内部的X射线。
可以看出,衰减分布中出现两个大的突变:大致是在信道编号270和信道编号460处。在这种情况下这是从空气向头颅的过渡。
现在为每个投影和每个信道如下检查测量数据组,在哪些区域内出现大的对比度。高对比度相当于彼此相邻的检测元件的测量值的大的值差。这种检查的结果作为曲线CON在图5中示出。数值CON可以按照不同的方式计算。一种具有优点的可能性如下:
S(i)是求对数的测量值,也就是曲线ATT的单值。参数i在此方面表示信道编号。为所有信道i按照下列凭经验推导的公式计算对比度CON:
CON(i)=(|S(i)-S(i+w)|+|S(i)-S(i-w)|)2
即计算在当前观察的信道的值S(i)与向一侧位移w个信道的信道的值S(i+w)之间的差,并计算在当前观察的信道的值S(i)与向另一侧位移w个信道的信道的值S(i-w)之间的差。
可以根据需要选取参数w。特别合适的是,w大约相当于EFS去卷积函数的滤波宽度的四分之一。在图5和6的具有EFS去卷积函数的49个滤波元件的例子中w选择为12。
为计算CON(i),将数值S(i)仅与两个处于确定的相邻区域内的数值S(i+w)和S(i-w)进行比较就足够了。分别仅使用所观察的信道的右侧和左侧的一个比较值就足够了,因为需要观察缓慢变化的结构,而不是局部噪声。与此相应,曲线CON的计算仅需很少时间。特别地,它-与每个单个信道i相关一与EFS去卷积相比计算强度小得多,该EFS去卷积不仅涉及另外两个检测元件的测量值,而且涉及大量其他检测元件的。
作为分别仅采用在右侧和左侧上的一个数值用于比较的例子的替换,也可以将每侧上的两个数值用于比较。在这种情况下,对比度CON的公式为:
CON(i)=(|S(i)-S(i+w)|+|S(i)-S(i-w)|+|S(i)-S(i+v)|+|S(i)-S(i-v)|)2
在这种情况下,v是与w不同的第二参数。
将数值CON(i)分别与可选择的(例如对所有信道N相同的)参数c0进行比较。从中确定图6B的曲线ACT。正如借助信道编号N所看到的那样,图6A和6B的曲线仅涉及图5所示区域的一部分。该部分包括图5曲线ATT的两个左侧波峰。
在各个数值CON(i)超过极限值c0的情况下,曲线ACT具有数值1,如果不是这种情况则具有数值0。数值0在这种情况下意味着,在各个检测元件的附近不存在大的测量值差别,而数值1则表示这种差别的存在。可以看出,曲线ACT仅在曲线ATT的陡峭上升的区域内具有数值1。
从曲线ACT可以知道,为哪些检测元件进行EFS去卷积:仅为曲线ACT具有数值1的那些检测元件进行EFS去卷积。因此这些测量值依据S(i)=>∑(i)通过去卷积的测量值取代,其中,∑(i)是EFS校正的数值。关于单个检测元件的EFS校正的计算在这种情况下可以按照传统的本身公知的方式进行。曲线ACT具有数值0的其他检测元件的测量值在无EFS校正的情况下被用于图像重建。即此时存在修改的EFS校正数据组,图像重建可以该数据组为依据。
在该经修改的数据组内不做改变地接收了绝大多数测量值,而对于所选择的测量值进行EFS校正。即,为进行测量数据的EFS校正,不是像传统那样校正所有测量值,而是仅校正在其相邻区域内确定了大对比度或大测量值差别的那些测量数据。这样做的优点是,用于计算新EFS校正的测量数据组所需的计算时间极大减少。该减少以哪个倍数进行,基本上取决于参数w和c0的选择以及所扫描的对象。正如已经提到的那样,在具有49个滤波元件的卷积核的图5和6的例子中选择为12。对于头颅检查的情况来说,c0例如可以选择为0.75,这在典型的扫描时使EFS去卷积所需的计算时间减少约70%-80%。
使用参数CON,其使得测量值仅在如下区域内通过EFS校正而改变,在该区域中存在大的对比度值。取决于参数CON来作出关于实施EFS校正的决定,是一种合理的方法,正如借助图6A所看到的那样。该图以缩放的方式示出图5的曲线CON。曲线DIF表示与原始的测量值相比的改变,其是在完全EFS校正时获得的。(曲线DIF的确定对图像重建来说没必要;它的进行只是为了演示该方法的合理性。)即为确定曲线DIF,对所有测量值S(i)进行EFS校正。从原始测量值中减去这些校正的测量值。这种比较的结果是曲线DIF。它因此表明在完全EFS校正的情况下哪些检测元件中出现变化。显而易见,曲线CON主要在曲线DIF与0明显不同的位置上明显不等于0。这一点意味着,仅对于值CON2与0明显不同的那些信道进行EFS校正就足够。
为了在与曲线ACT的突变相应的边界上不引起伪影,在这些边界上可以进行EFS校正的和未EFS校正的测量值的混合。通过这种叠化在边界上产生软过渡。为此计算叠化函数F(i),例如图6B的曲线FAD中所示。
然后通过下列公式进行叠化:
∑(i)=>F(i)*∑(i)+(1-F(i))*S(i)
也就是说,使用相关信道i的叠化函数F(i)的数值,以便从各个信道EFS的EFS校正的数据∑(i)和未EFS校正的数据S(i)中计算混合,也就是加权总和。
一个检测元件的测量值与其他检测元件的比较和取决于此的EFS校正借助单个检测行进行了说明。作为对此的替换或者补充,比较和EFS校正也可以在检测器的维中与之垂直地进行;在这种情况下,将检测元件与另一个检测行的检测元件进行比较并在该方向上使用EFS校正。该方向在图1和2中采用C9标注。该方法类似也可以在检测器的两个维中使用;在这种情况下,为EFS校正进行二维卷积。
前面借助实施例对本发明进行了说明。不言而喻,可以进行大量的变化和修改,而不偏离本发明的框架。
Claims (16)
1.用于从测量数据(p、ATT)中重建检查对象(O)图像数据(f)的方法,其中,所述测量数据(p、ATT)是在计算机断层扫描系统(C1)的发射焦点辐射和焦外辐射的辐射源(C2、C4)与检查对象(O)之间相对旋转运动时由检测器(C3、C5)采集的,
检查投影的测量数据(p、ATT)的彼此的差别(CON),
进行测量数据校正,以降低焦外辐射的效应,
其中,所述校正仅关于所选取的测量数据(p、ATT)进行,而测量数据(p、ATT)的选取则取决于检查结果,
在使用得到校正的测量数据(p、ATT)情况下进行图像数据(f)的重建。
2.按权利要求1所述的方法,其中,
仅对与其他测量数据(p、ATT)具有最小差别的测量数据(p、ATT)进行所述校正。
3.按权利要求1或2所述的方法,其中,
所述检测器(C3、C5)具有至少一个包括多个检测元件(N)的检测行,并在检查时将一个检测行的检测元件(N)的测量数据(p、ATT)与相同检测行的处于确定的距离上的另一个检测元件(N)的测量数据(p、ATT)进行比较。
4.按权利要求1或2所述的方法,其中,
所述检测器(C3、C5)具有至少一个包括多个检测元件(N)的检测行,并在检查时将一个检测行的检测元件(N)的测量数据(p、ATT)与相同检测行的处于确定的距离上的另外两个检测元件(N)的测量数据(p、ATT)进行比较。
5.按权利要求1或2所述的方法,其中,
所述检测器(C3、C5)具有至少一个包括多个检测元件(N)的检测行,并在检查时将一个检测行的检测元件(N)的测量数据(p、ATT)与相同检测行的处于确定的距离上的另外四个检测元件(N)的测量数据(p、ATT)进行比较。
6.按权利要求3-5之一所述的方法,其中,所述比较对检测行的所有检测元件进行。
7.按权利要求1-6之一所述的方法,其中,
所述测量数据校正关于一个检测元件(N)进行,方法是通过卷积将该检测元件(N)的测量数据(p、ATT)与相同检测行的一定数量的相邻检测元件(N)的测量数据(p、ATT)相关联。
8.按权利要求7和权利要求3-6之一所述的方法,其中,所述确定的距离大致相当于相邻检测元件(N)数量的四分之一。
9.按权利要求1-8之一所述的方法,其中,所述选取这样进行,使测量数据校正仅在测量数据(p、ATT)的对比度大的区域内起作用。
10.按权利要求1-9之一所述的方法,其中,
对于测量数据(p、ATT)的选取,使用显示需要进行校正的那些测量数据(p、ATT)的函数(ACT)。
11.按权利要求10所述的方法,其中,
对于在所述函数(ACT)从需要进行校正的测量数据(p、ATT)向不需要进行校正的测量数据(p、ATT)过渡的情况下的测量数据,计算已校正的和未校正的测量数据(p、ATT)的混合。
12.按权利要求1-11之一所述的方法,其中,
所述检测器(C3、C5)具有多个包括大量检测元件(N)的检测行,并在检查时将一个检测行的一个检测元件(N)的测量数据(p、ATT)与一个或多个其他检测行的一个或者多个其他检测元件(N)的测量数据(p、ATT)进行比较。
13.一种控制和计算单元(C10),用于从CT系统(C1)的测量数据(p、ATT)中重建检查对象(O)的图像数据(f),包括用于储存程序代码(Prg1-Prgn)的程序存储器,其中,在所述程序存储器中存在实施按权利要求1-12之一所述方法的程序代码(Prg1-Prgn)。
14.一种CT系统(C1),具有按权利要求13所述的控制和计算单元(C10)。
15.一种具有程序代码装置(Prg1-Prgn)的计算机程序,以便在计算机上执行计算机程序的情况下实施按权利要求1-12之一所述的方法。
16.一种计算机程序制品,包括储存在计算机可读数据载体上的计算机程序的程序代码装置(Prg1-Prgn),以便在计算机上执行计算机程序的情况下实施按权利要求1-12之一所述的方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102009048073.0 | 2009-10-01 | ||
DE200910048073 DE102009048073B3 (de) | 2009-10-01 | 2009-10-01 | CT-Datenbearbeitungsverfahren, CT-System, Computerprogramm und computerlesbarer Datenträger zur Reduzierung von Extrafokalstrahlungs-Effekten unter Minimierung des Rechenaufwandes |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102138804A true CN102138804A (zh) | 2011-08-03 |
CN102138804B CN102138804B (zh) | 2014-11-12 |
Family
ID=44312633
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010580280.XA Expired - Fee Related CN102138804B (zh) | 2009-10-01 | 2010-09-30 | 在计算开支最小化条件下降低焦外辐射效应的ct数据处理 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN102138804B (zh) |
DE (1) | DE102009048073B3 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106667515A (zh) * | 2012-11-05 | 2017-05-17 | 上海联影医疗科技有限公司 | Ct扫描仪散焦校正方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102011006579A1 (de) * | 2011-03-31 | 2012-10-04 | Siemens Aktiengesellschaft | Verfahren zur Erzeugung von Bilddaten eines Untersuchungsobjekts, Projektionsdatenverarbeitungseinrichtung, Röntgensystem und Computerprogramm |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5140520A (en) * | 1986-12-12 | 1992-08-18 | Yokogawa Medical Systems, Limited | Convolution processing method for x-ray tomographic apparatus |
US5757951A (en) * | 1994-12-20 | 1998-05-26 | Picker International, Inc. | Correction of off-focal radiation |
JP2002538704A (ja) * | 1999-03-01 | 2002-11-12 | ミラダ ソリューションズ リミテッド | X線画像処理方法及び装置 |
US6628744B1 (en) * | 2002-09-26 | 2003-09-30 | Koninklijke Philips Electronics N.V. | Off-focal radiation correction in CT |
DE10222702A1 (de) * | 2002-05-22 | 2003-12-24 | Siemens Ag | Verfahren zur Korrektur der Extrafokalstrahlung einer Röntgenröhre in der Computertomographie |
-
2009
- 2009-10-01 DE DE200910048073 patent/DE102009048073B3/de not_active Expired - Fee Related
-
2010
- 2010-09-30 CN CN201010580280.XA patent/CN102138804B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5140520A (en) * | 1986-12-12 | 1992-08-18 | Yokogawa Medical Systems, Limited | Convolution processing method for x-ray tomographic apparatus |
US5757951A (en) * | 1994-12-20 | 1998-05-26 | Picker International, Inc. | Correction of off-focal radiation |
JP2002538704A (ja) * | 1999-03-01 | 2002-11-12 | ミラダ ソリューションズ リミテッド | X線画像処理方法及び装置 |
DE10222702A1 (de) * | 2002-05-22 | 2003-12-24 | Siemens Ag | Verfahren zur Korrektur der Extrafokalstrahlung einer Röntgenröhre in der Computertomographie |
US6628744B1 (en) * | 2002-09-26 | 2003-09-30 | Koninklijke Philips Electronics N.V. | Off-focal radiation correction in CT |
Non-Patent Citations (1)
Title |
---|
王淑云等: "医用诊断X线管焦点辐射对焦点MTF的影响", 《泰山医学院学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106667515A (zh) * | 2012-11-05 | 2017-05-17 | 上海联影医疗科技有限公司 | Ct扫描仪散焦校正方法 |
CN106725567A (zh) * | 2012-11-05 | 2017-05-31 | 上海联影医疗科技有限公司 | 电子计算机x射线断层扫描仪 |
CN106725568A (zh) * | 2012-11-05 | 2017-05-31 | 上海联影医疗科技有限公司 | Ct扫描仪散焦强度测量方法 |
CN106667515B (zh) * | 2012-11-05 | 2019-11-19 | 上海联影医疗科技有限公司 | Ct扫描仪散焦校正方法 |
CN106725568B (zh) * | 2012-11-05 | 2019-12-20 | 上海联影医疗科技有限公司 | Ct扫描仪散焦强度测量方法 |
Also Published As
Publication number | Publication date |
---|---|
DE102009048073B3 (de) | 2011-07-07 |
CN102138804B (zh) | 2014-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6985482B2 (ja) | X線コンピュータ断層撮影装置、スキャン計画設定支援装置、医用画像診断システム、制御方法及び制御プログラム | |
EP1959835B1 (en) | Systems and methods for scanning and data acquisition in computed tomography (ct) applications | |
US7403588B2 (en) | Method for production of tomographic slice images of an examination object using at least two X-ray beams at offset angles, and a computed tomography scanner for carrying out this method | |
US9107642B2 (en) | Method and apparatus for tomographic X-ray imaging and source configuration | |
US8483361B2 (en) | Anode target for an x-ray tube and method for controlling the x-ray tube | |
US8208600B2 (en) | X-ray generating apparatus and X-ray CT apparatus using the same | |
US8983024B2 (en) | Tetrahedron beam computed tomography with multiple detectors and/or source arrays | |
US20050195935A1 (en) | X-ray CT system and X-ray apparatus | |
US7460635B2 (en) | X-ray CT apparatus, method of controlling the same, and program | |
US10709408B2 (en) | Medical image diagnosis apparatus and control method | |
US20160000396A1 (en) | X-ray ct apparatus and image processing apparatus | |
JP2004180715A (ja) | X線コンピュータ断層撮影装置 | |
JP2007144172A (ja) | モーション・アーティファクト補正と共にct画像再構成を実行する方法及びシステム | |
JP2009125250A (ja) | X線ct装置 | |
JP6595154B2 (ja) | X線ctを用いた画像診断装置 | |
US7949087B2 (en) | Radiography apparatus | |
JP6615439B2 (ja) | X線ct装置 | |
JP2021010727A (ja) | X線ctシステム及び医用処理装置 | |
CN102138804B (zh) | 在计算开支最小化条件下降低焦外辐射效应的ct数据处理 | |
JP2010142478A (ja) | X線ct装置 | |
CN211049378U (zh) | 一种大覆盖范围高速螺旋扫描ct机 | |
CN110881993A (zh) | X射线ct装置、医用图像处理装置以及x射线ct系统 | |
JP5823178B2 (ja) | X線ct装置 | |
JP5203750B2 (ja) | 心電同期スキャン方法及びx線コンピュータ断層撮影装置 | |
JP2011136002A (ja) | X線ct装置 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141112 Termination date: 20190930 |