CN102160085A - 包括针对散射辐射的校正单元的成像装置 - Google Patents
包括针对散射辐射的校正单元的成像装置 Download PDFInfo
- Publication number
- CN102160085A CN102160085A CN2009801360618A CN200980136061A CN102160085A CN 102160085 A CN102160085 A CN 102160085A CN 2009801360618 A CN2009801360618 A CN 2009801360618A CN 200980136061 A CN200980136061 A CN 200980136061A CN 102160085 A CN102160085 A CN 102160085A
- Authority
- CN
- China
- Prior art keywords
- value
- attenuating elements
- radiation
- scattering
- probe
- 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
- 230000005855 radiation Effects 0.000 title claims abstract description 162
- 238000003384 imaging method Methods 0.000 title claims abstract description 66
- 238000012937 correction Methods 0.000 title abstract description 7
- 239000000523 sample Substances 0.000 claims description 129
- 238000005259 measurement Methods 0.000 claims description 22
- 238000004590 computer program Methods 0.000 claims description 21
- 238000000034 method Methods 0.000 claims description 20
- 238000012360 testing method Methods 0.000 claims description 16
- 230000015572 biosynthetic process Effects 0.000 claims description 13
- 238000003325 tomography Methods 0.000 claims description 11
- 230000001915 proofreading effect Effects 0.000 claims description 10
- 230000002238 attenuated effect Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 24
- 230000000694 effects Effects 0.000 description 7
- 230000001419 dependent effect Effects 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 101000632261 Homo sapiens Semaphorin-3A Proteins 0.000 description 2
- 101000650811 Homo sapiens Semaphorin-3D Proteins 0.000 description 2
- 102100027974 Semaphorin-3A Human genes 0.000 description 2
- 102100027746 Semaphorin-3D Human genes 0.000 description 2
- 230000002411 adverse Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000003760 hair shine Effects 0.000 description 2
- 238000002600 positron emission tomography Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012633 nuclear imaging Methods 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000001959 radiotherapy Methods 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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 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
-
- 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4291—Arrangements for detecting radiation specially adapted for radiation diagnosis the detector being combined with a grid or grating
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及一种用于产生对象的感兴趣区域的图像的成像装置。所述成像装置包括:辐射源(2),其用于发射辐射(4);以及探测器(6),其用于测量穿过所述感兴趣区域之后的辐射(4)并用于根据所测量的辐射(4)产生测量的探测值。所述成像装置还包括:衰减元件,其用于在所述辐射(4)穿过所述感兴趣区域之前衰减所述辐射(4);以及衰减元件散射值提供单元(12),其用于提供衰减元件散射值,该值取决于所述衰减元件引起的辐射(4)的散射。探测值校正单元(17)基于所提供的衰减元件散射值校正所述测量的探测值,并且重建单元(18)根据经校正的探测值重建感兴趣区域的图像。
Description
技术领域
本发明涉及一种用于产生对象的感兴趣区域的图像的成像装置、成像方法和计算机程序。本发明还涉及一种用于校正由成像装置采集的探测值的校正装置、校正方法和校正计算机程序。
背景技术
US 7,065,234B2公开了一种计算断层摄影***,其包括:辐射源,其发射辐射;以及探测器,其用于探测穿过感兴趣区域之后的辐射,以用于根据所探测的辐射产生探测值。针对散射效应校正探测值,并使用经校正的探测值重建感兴趣区域的图像。这一经重建的图像仍然包括由散射效应引起的图像伪影,由此降低了经重建的图像的质量。
发明内容
本发明的目的在于提供一种用于产生对象的感兴趣区域的图像的成像装置、成像方法和成像计算机程序,其中改善了感兴趣区域的重建图像的质量。本发明的另一目的是提供一种用于校正探测值的校正装置、校正方法和校正计算机程序,其校正探测值,从而改进利用经校正的探测值重建的图像的质量。
在本发明的第一方面中,提供了一种用于产生对象的感兴趣区域的图像的成像装置,其中所述成像装置包括:
-辐射源,其用于发射辐射,
-探测器,其用于测量穿过所述感兴趣区域之后的辐射并且根据所测量的辐射产生测量的探测值,
-衰减元件,其用于在辐射穿过所述感兴趣区域之前衰减所述辐射,
-衰减元件散射值提供单元,其用于提供衰减元件散射值,该值取决于所述衰减元件引起的辐射的散射,
-探测值校正单元,其用于基于所提供的衰减元件散射值校正所述测量的探测值,
-重建单元,其用于根据经校正的探测值重建感兴趣区域的图像。
探测器测量的探测值受到辐射在衰减元件中散射的不利影响。于是,通过针对衰减元件中的辐射散射校正探测值,即,通过提供取决于衰减元件引起的辐射散射的衰减元件散射值并且通过基于所确定的衰减元件散射值校正探测值,改善了探测值的质量,从而改善了根据经校正的探测值重建的图像的质量。
成像装置优选是计算断层摄影装置,其中X射线源发射X射线辐射并相对于对象的感兴趣区域旋转,以从不同角度方向照射感兴趣区域。探测器探测从不同方向照射感兴趣区域的,穿过衰减元件和感兴趣区域之后的辐射并且产生取决于所探测的辐射的探测值。计算断层摄影装置的衰减元件优选为楔形物。可以将楔形物称为蝴蝶结或射束成形器。在计算断层摄影装置中使用楔形物是为了降低在横向对象区域中通过对象的辐射量,对象优选为患者。具体而言,楔形物在要成像的这种对象(具体而言,患者)薄的地方是厚的,而楔形物在要成像的这种对象预计厚的地方是薄的。衰减元件散射值提供单元提供取决于楔形物中X射线辐射的散射的衰减元件散射值,探测值校正单元基于所确定的衰减元件散射值校正测量的探测值。重建单元优选适于使用反投影算法或另一种计算断层摄影重建算法根据经校正的探测值重建感兴趣区域的图像。
术语“感兴趣区域的”包括整个对象或仅包括对象的一部分。
辐射优选形成像锥形射束或扇形射束的射束,其中对象的尺寸和布置优选使得,如果辐射源适于在测量探测值的同时绕旋转轴旋转,则对象在垂直于旋转轴的方向上完全位于射束之内。
应当指出,“基于所提供的衰减元件散射值校正测量的探测值”这一表述也包括进一步处理衰减元件散射值,例如以便考虑对象的衰减,其中优选基于这些进一步处理的衰减元件散射值校正测量的探测值。
优选成像装置还包括用于确定对象散射值的对象散射确定单元,对象散射值取决于对象对辐射的散射,其中探测值校正单元适于基于所确定的对象散射值校正测量的探测值。由于探测值也可能受到辐射在对象中的散射的不利影响,因此考虑这些额外的散射效应,即,确定取决于辐射在对象中的散射的对象散射值以及还基于所确定的对象散射值校正探测值,进一步提高了探测值的质量,从而提高了利用经校正的探测值重建的感兴趣区域的图像的质量。
优选通过已知方法,具体而言,通过例如在WO 2007/148263A1中公开的基于与对象特有卷积核的卷积的独立方法,确定对象散射值并基于所确定的对象散射值校正测量的探测值,在此通过引用将该文献并入。
进一步优选地,该成像装置包括对象所处的检查区以及用于准直辐射源发射的辐射的准直器,准直器位于衰减元件和检查区之间,其中辐射源、探测器和准直器适于针对辐射的不同准直测量校准值,并且其中衰减元件散射值提供单元适于基于校准值确定衰减元件散射值。优选地,衰减元件散射值提供单元适于通过将一函数拟合到校准值来确定衰减元件散射值,所述函数考虑了校准值的衰减元件散射部分和未散射部分并且取决于不同的准直。进一步优选地,该函数对于描述校准值的衰减元件散射部分随准直变大而变大具有平方根相关性,对于描述校准值的未散射部分具有常数相关性,其中衰减元件散射值取决于衰减元件散射部分。由于探测值中由辐射在衰减元件中的散射引起的部分具有取决于辐射准直的平方根特性,因此可以通过将具有平方根相关性的函数拟合到针对辐射的不同准直测量的校准值,而容易且高精度地确定这一部分,从而确定取决于这一部分的衰减元件散射值。
进一步优选地,
-探测器包括若干探测器元件,
-辐射源、探测器和准直器适于为探测器元件的子组针对辐射的不同准直测量校准值,
-衰减元件散射值提供单元适于针对探测器元件的子组确定衰减元件散射值,
-成像装置还包括衰减元件散射值比提供单元,其用于提供探测器元件的子组的衰减元件散射值与子组之外的探测器元件的衰减元件散射值的比值,
-所述衰减元件散射值提供单元适于通过将针对子组的衰减元件散射值乘以所提供的比值来确定针对子组之外的探测器元件的衰减元件散射值。这允许测量校准值并基于仅针对探测器元件的子组的校准值确定衰减元件散射值,并且简单地通过将探测器元件的子组的衰减元件散射值乘以所提供的比值来确定针对其他探测器元件的衰减元件散射值,由此减少确定衰减元件散射值所需的工作量和时间。
进一步优选地,探测器包括二维探测表面,其中探测器元件的子组是所述探测表面上的一条线上的探测器。这条线优选是中心线,具体而言,如果成像装置为计算断层摄影装置,这条线优选相对于平行于旋转轴的方向位于中心。
优选调整辐射源、准直器和探测器,从而使得对于用于确定校准值的每种准直,探测器元件的子组,优选地探测器元件的线,进一步优选探地测器元件的中心线,被完全照射。
衰减元件散射值比提供单元优选为模拟单元,其中通过对成像装置建模,具体而言,通过对辐射源、探测器和衰减元件建模,并通过执行Monte-Carlo模拟来对衰减元件散射值的产生建模。如果成像装置为计算断层摄影装置,则衰减元件散射值比提供单元优选适于针对垂直于旋转轴的探测器元件的中心线与也垂直于旋转轴的另一条探测器元件的线的每种组合提供衰减元件散射值比,具体而言,衰减元件散射值的单一比值。
在实施例中,衰减元件散射值比提供单元是已经存储了衰减元件散射值比的存储单元,或者它可以例如通过使用上述模拟计算这些比值,其中优选地,对于位于平行于旋转轴的同一条线上且位于垂直于旋转轴的探测器元件的中心线上以及也垂直于旋转轴的另一条探测器元件的线上的每一对探测器元件,确定衰减元件散射值比,并且其中对针对中心线和单根另一条线确定的衰减元件散射值比求平均值,以针对垂直于旋转轴的中心线和也垂直于旋转轴的另一条单根线的组合确定单一衰减元件散射值比。这减少了衰减元件散射值比提供单元必须要提供的衰减元件散射值比的量。
进一步优选地,该成像装置还包括:
-空气扫描值提供单元,其用于提供空气扫描值,以及
-空气扫描值校正单元,其用于基于所述衰减元件散射值校正空气扫描值,
其中所述探测值校正单元适于基于经校正的空气扫描值校正探测值。
成像装置例如是计算断层摄影装置,并且其从不同方向照射感兴趣区域,其中优选以不同角范围,例如10度,划分不同方向。优选通过以下方式确定空气扫描值:在感兴趣区域中没有对象且辐射源绕感兴趣区域转动的同时采集探测值,并对属于相同角度范围的探测值求平均值来确定相应角度范围的空气扫描值。于是,优选针对不同采集条件,针对每个角度范围确定空气扫描值。例如,如果辐射源为X射线管,则由辐射源的管电流、准直、旋转速度等定义采集条件。由于探测值校正单元使用已经基于衰减元件散射值校正的空气扫描值来校正探测值,所以进一步提高了探测值的质量,从而提供了利用经校正的探测值重建的图像的质量。空气扫描值提供单元优选是已经存储了空气扫描值的存储单元。在另一实施例中,空气扫描值提供单元例如通过如上所述对感兴趣区域中没有对象时采集的探测值求平均值来确定空气扫描值。
进一步优选地,探测值校正单元适于针对已经在衰减元件中被散射并被对象(如果在衰减元件和探测器之间存在对象)衰减的辐射的衰减,调整衰减元件散射值,并基于所调整的衰减元件散射值校正所测量的探测值。优选基于测量的探测值执行这种调整。由于在用于测量对象的探测值的测量中,已经被衰减元件散射的辐射被对象衰减,因此必须要调整衰减元件散射值,从而考虑对象的这种衰减,所述衰减元件散射值考虑到衰减元件引起的辐射散射而未考虑对象造成的散射辐射的衰减。也可以将这定义为衰减元件散射值通过对象的传播。分别针对对象的衰减或衰减元件散射值通过对象的传播对衰减元件散射值进行这种调整进一步改善了经校正的探测值的质量,从而改善了利用经校正的探测值重建的图像的质量。
进一步优选地,通过将相应衰减元件散射值乘以取决于对象引起的衰减的相应对象衰减值执行衰减元件散射值的调整。对象衰减值优选是由对象引起的相应衰减或对象引起的这种衰减的低通滤波值。这种低通滤波器优选包括低通滤波器系数,根据***几何结构和位于探测器的探测器元件间的所用抗散射网格的接收角度调整低通滤波器系数。调整这些低通滤波器系数,从而使得由对象引起的滤波后衰减,在这种情况下即对象衰减值,反映了从探测器像素看到的且在抗散射网格的整个接收角度上积分的由对象引起的平均衰减。因此,进一步优选地,低通滤波器是二维的,且针对每种探测器尺度单独地选择低通系数。在一个优选实施例中,将低通滤波器实现为二项式低通滤波器。
进一步优选地,
-成像装置包括:
a)空气扫描值提供单元,其用于提供空气扫描值,以及
b)空气扫描值校正单元,其用于基于衰减元件散射值校正空气扫描值,
-探测值校正单元适于迭代地校正测量的探测值,其中
a)在迭代步骤中,通过从除以相应经校正的空气扫描值的相应测量的探测值减去乘以相应对象衰减值的相应衰减元件散射值来计算经校正的探测值,所述相应对象衰减值取决于由所述对象(31)引起的衰减,
b)一开始,相应对象衰减值取决于相应测量的探测值除以相应经校正的空气扫描值,在进一步迭代步骤中,对象衰减值取决于计算的相应经校正的探测值,即在前面的迭代步骤中计算出的所计算的相应探测值。
进一步优选地,探测值校正装置,尤其是在迭代步骤中,还适于也从除以相应经校正的空气扫描值的相应测量的探测值减去除以相应经校正的空气扫描值的相应对象散射值,以计算经校正的探测值。
应当理解,将衰减元件散射值、对象散射值、经校正的空气扫描值和对象衰减值分配给相应探测值。例如,将衰减的衰减元件散射值、对象散射值、经校正的空气扫描值和对象衰减值分配给对应于探测器的相同探测器元件和相同采集条件,即例如辐射源的相同角度位置、相同准直、相同管电流(如果辐射源为X射线管)等的探测值。优选将相应的衰减元件散射值、相应的对象散射值、相应的经校正的空气扫描值和相应的对象衰减值分配给同一探测值。
成像装置优选适于针对探测器的探测器元件和采集条件的每种组合确定衰减元件散射值。进一步优选地,成像装置针对这些组合的每种确定对象散射值、对象衰减值和经校正的空气扫描值。
在实施例中,衰减元件散射值提供单元是已经存储了衰减元件散射值的存储单元。衰减元件散射值提供单元还可适于例如通过Monte-Carlo模拟的成像装置的模拟,来确定衰减元件散射值。
在本发明的另一方面中,提供了一种用于校正测量的探测值的校正装置,为所述校正装置提供根据测量的辐射产生的探测值,所述辐射是由辐射源发射的并已穿过衰减元件,所述衰减元件用于在辐射穿过感兴趣区域之前衰减所述辐射,所述校正装置包括:
-衰减元件散射值提供单元,其用于提供取决于衰减元件引起的辐射的散射的衰减元件散射值,
-探测值校正单元,其用于基于所提供的衰减元件散射值校正测量的探测值。
优选地,校正装置还包括用于根据经校正的探测值重建感兴趣区域的图像的重建单元。
在本发明的另一方面中,提供了一种用于产生对象的感兴趣区域的图像的成像方法,其中该成像方法包括以下步骤:
-由辐射源发射辐射,
-在辐射穿过感兴趣区域之前,由衰减元件衰减所述辐射,
-探测器测量穿过感兴趣区域之后的辐射并根据所测量的辐射产生测量的探测值,
-衰减元件散射值提供单元提供取决于衰减元件引起的辐射的散射的衰减元件散射值,
-探测值校正单元基于所提供的衰减元件散射值校正测量的探测值,
-重建单元根据经校正的探测值重建感兴趣区域的图像。
在本发明的另一方面中,提供了一种用于校正测量的探测值的校正方法,其中该校正方法包括以下步骤:
-提供根据测量的辐射产生的测量的探测值,所述辐射是由辐射源发射的并且已经穿过衰减元件,所述衰减元件用于在辐射穿过感兴趣区域之前衰减所述辐射,
-衰减元件散射值提供单元提供取决于衰减元件引起的辐射的散射的衰减元件散射值,
-探测值校正单元基于所提供的衰减元件散射值校正所测量的探测值。
优选地,校正方法还包括重建单元根据经校正的探测值重建感兴趣区域的图像的步骤。
在本发明的另一方面中,提供了一种用于产生对象的感兴趣区域的图像的成像计算机程序,其中该成像计算机程序包括用于当成像计算机程序在控制根据权利要求1所述的成像装置的计算机上运行时,令所述成像装置执行根据权利要求12所述的成像方法的步骤的程序代码段。
在本发明的另一方面中,提供了一种用于校正探测值的校正计算机程序,其中校正计算机程序包括用于当校正计算机程序在控制根据权利要求11所述的校正装置的计算机上运行时,令所述校正装置执行根据权利要求13所述的校正方法的步骤的程序代码段。
应当理解,根据权利要求1所述的成像装置、根据权利要求11所述的校正装置、根据权利要求12所述的成像方法、根据权利要求13所述的校正方法、根据权利要求14所述的成像计算机程序和根据权利要求15所述的校正计算机程序具有与从属权利要求中限定的类似和/或相同的优选实施例。
应当理解,本发明的优选实施例还可以是从属权利要求与相应的独立权利要求的任意组合。
附图说明
参考下文描述的实施例,本发明的这些和其他方面将显而易见并得到阐述。在下述附图中:
图1示意性和示范性地示出了用于产生对象的感兴趣区域的图像的成像装置的实施例;
图2示意性和示范性地示出了成像装置的若干元件的布置;
图3示意性和示范性地示出了成像装置的探测器的探测表面;
图4示意性和示范性地示出了校准值与不同准直的相关性;
图5示出了流程图,其图示出了用于产生对象的感兴趣区域的图像的成像方法的实施例;
图6示出了流程图,其图示出了衰减元件散射值的确定;
图7示意性和示范性地示出了衰减元件散射值与不同准直的相关性;以及
图8根据不同准直和沿探测器探测表面上的中心线定位的不同探测器元件示意性和示范性地示出了衰减元件散射值。
具体实施方式
图1示意性和示范性地示出了用于对对象的感兴趣区域进行成像的成像装置。在这一实施例中,成像装置为计算断层摄影装置。该计算断层摄影装置包括能够围绕旋转轴R旋转的扫描架1,所述旋转轴R平行于z方向延伸。将辐射源2安装到扫描架1上,在这一实施例中,所述辐射源2为X射线管。辐射源2设有准直器3,在这一实施例中,所述准直器从辐射源2产生的辐射形成锥形辐射束4。所述辐射穿过诸如患者的对象(图1中未示出)和优选在检查区5中位于对象内的感兴趣区域,在这一实施例中,所述检查区是圆柱形的。在穿过检查区5之后,辐射4入射到探测器6上,探测器6包括二维探测表面。探测器6安装在扫描架1上。
计算断层摄影装置包括两个电机7、8。优选通过电机7以恒定的但是可调整的角速度驱动扫描架。电机8被设置为使诸如患者的对象平行于旋转轴R或z轴的方向发生位移,所述对象被安置在检查区5内的患者台上。例如,由控制单元9控制这些电机7、8,从而使得辐射源2和检查区5,并且因此检查区5内的感兴趣区域沿着螺旋形轨迹彼此相对移动。但是,也可能不使对象或检查区5移动,而只旋转辐射源2,即,辐射源2相对于对象或检查区5沿着圆形轨迹移动。此外,在另一实施例中,准直器3可以适于形成另一种射束形状,具体而言为扇形射束,探测器6可以包括探测表面,其根据另一种射束形状,具体而言,扇形射束成形。
在辐射源2和检查区5的相对运动期间,探测装置6根据入射在探测器6的探测表面上的辐射产生探测值。
成像装置,在这一实施例中即计算断层摄影装置,还包括衰减元件30,其用于在辐射4穿过对象31,从而穿过感兴趣区域之前对其进行衰减。图2中示范性和示意性地示出了衰减元件30。图2中所示的衰减元件30是楔形的,从而使得其在区域32中具有更大的厚度,其中辐射束4穿过该区域32并且仅穿过对象(应当是患者)的较薄部分,并且具有较薄区域33,其沿穿过对象31较厚区域的辐射束定位。衰减元件30的形状优选像具有中心凹陷的长方体,在垂直于旋转轴R或z方向的相交平面中为椭圆形。
参考图2,探测器6包括二维探测表面34,探测表面34包括若干探测器元件和抗散射网格。图2中示意性和示范性地示出了分配到探测器元件35的抗散射网格的叶片36。虚线37表示角度范围,从该角度范围,辐射可以进入叶片36之间的空间并被探测器元件35探测到。图2还图示了衰减元件30和对象31中的若干散射事件。
射束42在衰减元件30中的散射位置43处被散射,由此产生散射辐射束44,其在对象31中的散射位置45处也被散射,由此产生被探测器元件35探测到的另一散射束40。射束46在衰减元件30中的散射位置47处被散射,由此产生散射束41,散射束41穿过对象31而未被再次散射,其也被探测器元件35探测到。辐射束48穿过衰减元件30而在衰减元件30中未被散射,但这一辐射束48在对象31中的散射位置49处被散射,由此产生散射束39,其也被探测器元件35探测到。最后,辐射束38穿过衰减元件30和对象31,均未被散射,这一辐射束可以被视为初级辐射,其也被探测器元件35探测到。于是,探测器元件35探测到的探测值不仅包括初级辐射38,而且还包括散射辐射39、40和41。在下文中,忽略了两次散射辐射40,因为它对探测值的贡献小到可以忽略。
应当指出,附图仅仅是示范性和示意性的,并未按比例绘制。例如,在图2中,与总的探测器表面34相比,探测器元件35过大。
在这一实施例中,探测值是投影值并是由探测器6探测到的,将探测值提供给校正设备10,用于针对散射辐射校正探测值。具体而言,校正设备适于减少或消除探测值中的散射部分39、40和41,以便获得优选仅取决于初级辐射38的经校正的探测值。可以将校正设备10视为校正装置。
校正装置10包括衰减元件散射值提供单元12,其用于提供取决于衰减元件30引起的辐射散射的衰减元件散射值。为了确定这些衰减元件散射值,辐射源2、探测器6和准直器3适于针对辐射4的不同准直情况测量校准值,其中衰减元件散射值提供单元12适于基于校准值确定衰减元件散射值。准直器3位于检查区5和衰减元件30之间,其中准直器3适于提供准直,平行于z方向或旋转轴R准直是不同的,即,不同的准直对应于在平行于z方向或旋转轴R的方向上的不同延伸。优选地,准直器3包括两个叶片,它们形成狭缝,其中不同准直对应于在平行于z方向或旋转轴R的方向上不同的狭缝宽度。
优选地,在没有对象31时,针对探测器6探测表面34上的探测器元件50的中心线51确定校准值。这一中心线51相对于旋转轴R位于探测表面34的中心,在图3中示意性且示范性地示出了这一中心线。中心线51垂直于旋转轴R。改变准直器3的准直并测量中心线51的探测值,其中针对不同准直测量的这些探测值为校准值。选择不同的准直,从而使得中心线51始终完全被辐射源2发射的辐射4的初级辐射照射。
图4示意性和示范性地示出了分别由中心线51的探测器元件以两种校准coll 1、coll 2测量的两个校准值M(coll 1)、M(coll 2)。衰减元件散射值提供单元12适于通过将一函数拟合到校准值来确定衰减元件散射值,该函数考虑了校准值的衰减元件散射部分和未散射部分并且取决于不同的准直。在这一实施例中,该函数对于描述校准值的衰减元件散射部分随准直变大而变大具有平方根相关性,对于描述校准值的未散射部分具有常数相关性,其中衰减元件散射值取决于衰减元件散射部分。被拟合到校准值的函数优选由以下方程定义:
其中
在方程(1)中,Mi,j(coll)表示探测值,在这种情况下为校准值,是利用i表示的探测器元件和j表示的辐射源角位置以准直coll测量的。Pi,j这一项表示针对i表示的探测器元件和j表示的辐射源位置拟合的探测值(在这种情况下为校准值)的未散射部分。探测值(在这种情况下即校准值)的衰减元件散射部分由Si,j(coll)表示,其具有如方程(2)所示的平方根特性。bi,j这一项表示针对衰减元件散射部分Si,j(coll)的平方根特性的拟合参数。
在这一实施例中,由以下方程定义针对i表示的探测器元件和j表示的辐射源位置的衰减元件散射值SPRi,j(coll):
也可以将衰减元件散射值SPRi,j(coll)视为衰减元件散射特征。
在实施例中,并非根据辐射源的角位置确定衰减元件散射值。在这种情况下,可以针对单个辐射源位置确定校准值,或者可以对针对不同辐射源位置且针对相同准直和相同探测器元件确定的衰减元件散射值求平均值。
在已经针对探测器6的探测表面34的中心线51确定衰减元件散射值之后,衰减元件散射值比提供单元13提供中心线51的衰减元件散射值与中心线51之外的探测器元件的衰减元件散射值的比值,其中所述中心线51形成探测表面34上探测器元件50的子组。在这一实施例中,衰减元件散射值比提供单元13适于针对垂直于旋转轴R且不是中心线51的线和中心线51的每种组合提供衰减元件散射值比,具体而言,提供单一衰减元件散射值比。例如,如果已经针对中心线51的第一探测器元件53确定了衰减元件散射值,则将这一衰减元件散射值乘以已针对中心线51和第二探测器元件54所在的外部线52提供的比值,以确定第二探测器元件54的衰减元件散射值。通过这种方式,对于中心线51之外的所有探测器元件,通过将已针对中心线51的探测器元件确定的相应衰减元件散射值乘以衰减元件散射值比提供单元13提供的相应比值来确定衰减元件散射值。在另一实施例中,衰减元件散射值比提供单元13适于针对沿平行于旋转轴R的线定位的探测器元件和中心线51上位于平行于旋转轴R的同一条线上的相应探测器元件的每种组合提供衰减元件散射值比,即,在这一另一实施例中,并不是通过将为中心线51确定的对应衰减元件散射值乘以同一衰减元件散射值比来确定线52的所有衰减元件散射值。
衰减元件散射值比提供单元13优选为存储衰减元件散射值比的存储单元,并且/或者衰减元件散射值比提供单元13可以适于通过如下方式确定衰减元件散射值比:模拟探测器6探测的探测值的衰减元件散射部分的产生,从而通过对探测过程建模,即,例如通过对辐射源、衰减元件和探测器建模,并利用Monte-Carlo模拟来模拟衰减元件散射值的产生,其中通过计算模拟的衰减元件散射值的相应比值来确定该比值。
在另一实施例中,衰减元件散射值提供单元可以适于在不使用校准值的情况下,利用整个成像***的Monte Carlo模拟来确定衰减元件散射值。
校正设备10还包括对象散射确定单元14,其用于确定取决于对象31对辐射的散射的对象散射值。优选通过基于如WO 2007/148263A1中公开的,与对象特有卷积核的卷积的独立方法来确定这些对象散射值在此通过引用将该文献并入。
校正设备10还包括用于提供空气扫描值的空气扫描值提供单元15和用于基于衰减元件散射值校正空气扫描值的空气扫描值校正单元16。优选根据以下方程执行对经校正的空气扫描值的确定:
空气扫描值提供单元15优选是存储了在无对象时测量的且优选预先计算出的空气扫描值的存储单元。必须要理解,未经校正的空气扫描值仍然包含来自衰减元件的散射辐射的贡献。
校正装置10还包括用于基于所确定的衰减元件散射值校正探测值的探测值校正单元17。在这一实施例中,探测值校正单元17适于基于所确定的衰减元件散射值、所确定的对象散射值和所确定的经校正的空气扫描值来校正探测值。探测值校正单元17优选适于根据以下方程校正探测值
其中
且
在方程(6)中,表示探测值在衰减元件30中被散射且被对象31衰减的部分。这一项表示仅由对象31造成的辐射衰减,而表示对象造成的辐射的低通滤波的衰减。在这一实施例中,被视为对象的衰减值。在另一实施例中,可以省略低通滤波器,可以将视为对象衰减值。
低通滤波器优选包括低通滤波器系数,根据***几何结构和位于探测器的探测器元件之间的所用抗散射网格的接收角度调节低通滤波器系数。调节这些低通滤波器系数,从而使得由对象引起的经滤波的衰减,在这种情况下即对象衰减值,反映了从探测器像素看到的且在抗散射网格的整个接收角度上积分的由对象引起的平均衰减。因此,进一步优选地,低通滤波器是二维的,且针对每种探测器尺度单独地选择低通系数。在一个优选实施例中,将低通滤波器实现为二项式低通滤波器。
作为范例,可以通过以下方式找到适当的低通系数:假设探测器中心处的探测器像素中心中有虚拟辐射源,针对若干离散位置计算抗散射网格导致的虚拟衰减,所述离散位置分布在焦斑位置周围,而这些位置的间距在等角点处投影的有效间距应当与在连接焦斑和相邻探测器像素的相邻线之间找到的有效间距相同。可以通过如下方式找到用于适当二维二项式低通的系数:选择用于每个分立尺度的系数,从而使得两个尺度的卷积密切反映上述一般低通的特征。
应当指出,在这一实施例中,为了校正测量的探测值,需要对象衰减值即仅由对象造成的衰减不过,方程(5)需要的方程(6)中使用的项等价于项,该项仅由方程(5)计算,因此在评估方程(5)期间不容易得到。因此,优选实施例是使用迭代过程,其中,对于迭代指数k=1的第一次迭代,利用测量的探测值除以经校正的空气扫描值对进行初始化,在每次进一步迭代中,由上次迭代的校正值更新在一次迭代到下一次时观察到小于预定义阈值时,可以终止迭代。在另一实施例中,将迭代次数固定为预定义次数,例如3次。
在方程(4)到(6)中,假设变量对应于相同的准直。因此在这些方程中未示出对准直的相关性。
校正设备10还包括用于根据经校正的探测值重建感兴趣区域的图像的重建单元18。重建单元18优选适于执行反投影算法以重建感兴趣区域的图像。在另一实施例中,重建单元18可以适于使用另一种重建算法,例如拉冬反演。
向显示器单元11提供重建单元18重建的图像以显示重建图像。
校正设备优选也受到控制单元9的控制。
在下文中,将参考图5中所示的流程图描述用于产生对象的感兴趣区域的图像的成像方法。
在这一实施例中,用于产生对象的感兴趣区域的图像的成像方法是计算断层摄影方法。在步骤101中,辐射源2绕着旋转轴R旋转,而对象或检查区5不移动,即辐射源2沿着绕对象或检查区5的圆形轨迹移动。在另一实施例中,辐射源2相对于对象沿着另一轨迹,例如螺旋形轨迹移动。辐射源2发射至少在感兴趣区域的中穿过衰减元件和对象的辐射。探测器6探测到已穿过衰减元件和对象的辐射,产生测量的探测值,即投影值。将测量的探测值传输到校正装置10。
在步骤102中,提供取决于衰减元件30引起的辐射4散射的衰减元件散射值,在步骤103中,空气扫描值提供单元15提供空气扫描值。可以彼此独立地执行步骤101、102和103。在步骤104中,空气扫描值校正单元16基于衰减元件散射值校正空气扫描值。可以在执行步骤101之前或之后执行这一步骤104。优选地,在步骤101之前执行步骤104。
在步骤105中,对象散射确定单元14根据在步骤101中探测的测量的探测值确定对象散射值,在步骤106中,探测值校正单元17基于所确定的对象散射值、所提供的衰减元件散射值和经校正的空气扫描值校正探测值。
在步骤107中重建单元18使用经校正的探测值重建感兴趣区域的图像,在步骤108中,在显示器单元11上示出重建图像。
在下文中将参考图6所示的流程图描述确定衰减元件散射值的实施例。
在步骤201中,由具有不同准直的辐射照射探测器6的探测表面34,其中针对每种准直,读出中心线51的探测值。这些探测值被视为校准值。对于校准值的每次测量,用准直的辐射完全照射中心线51。
在步骤202中,衰减元件散射值提供单元12基于校准值确定衰减元件散射值,衰减元件散射值取决于衰减元件30引起的辐射4的散射。在这一实施例中,将具有平方根特性的上述拟合函数拟合到校准值,并使用方程(3)确定针对中心线51的衰减元件散射值。
在步骤203,衰减元件散射值比提供单元13提供中心线51的衰减元件散射值与探测表面34上其他线中的探测器元件的衰减元件散射值的衰减元件散射值比。此外,将这些衰减元件散射值比乘以针对中心线51确定的衰减元件散射值,以确定也用于探测表面34其他线的衰减元件散射值。
如果在一实施例中,校准值也用作空气扫描值或空气扫描值用作校准值,可以分别省略步骤201或步骤103。
也可以将校正设备10视为独立的***,为其提供例如在步骤101中采集的测量的探测值,其中校正设备至少包括:衰减元件散射值提供单元,其用于提供取决于衰减元件引起的辐射的散射的衰减元件散射值;以及探测值校正单元,其用于基于所提供的衰减元件散射值校正探测值。也可以被视为校正装置的这一校正设备优选还包括用于根据经校正的探测值重建感兴趣区域的图像的重建单元。一种对应的校正方法至少包括提供探测值、提供衰减元件散射值以及基于所提供的衰减元件散射值校正探测值的步骤。
在步骤201中采集校准值的同时,可以执行若干采集,具体而言,可以执行在辐射源的不同位置处的若干采集,以便产生冗余校准值,具体而言,优选对冗余校准值进行平均,以便使量子噪声保持小值。
当在步骤202中确定衰减元件散射值期间,具体而言,利用具有平方根特性的拟合函数确定衰减元件散射值期间,优选仅使用针对较大准直探测的,即针对大叶片开口探测的校准值,具体而言,仅使用对应于大于预定义准直值的准直的校准值。优选选择这一预定义准直值,从而使得校准值的质量足以用于确定可靠的衰减元件散射值。具体而言,优选仅使用利用与探测表面34上垂直于旋转轴R的八条或更多条线探测器元件对应的准直测量的校准值。
图7示意性和示范性地示出了针对探测器元件和固定辐射源位置并针对不同准直coll确定的衰减元件散射值SPR。在这一实施例中,由准直器3的狭缝在旋转轴R方向上的宽度定义准直。十字表示测量的校准值,线条表示拟合函数的被拟合平方根部分。
图8示意性和示范性地示出了针对探测器的探测表面的中心线,根据中心线中不同探测器元件d并根据不同准直确定的衰减元件散射值SPR,其中不同准直是由准直器3的狭缝在平行于旋转轴的方向上的宽度定义的。术语“x sl”,即例如“4 sl”表示选择准直,从而使得仅照射垂直于旋转轴R的x条线。
上述方法允许通过经每幅视图中看到的衰减传播衰减元件散射估计在像患者那样的对象存在时探测的衰减元件散射效应。因此将所确定的衰减元件散射值乘以根据初级辐射的衰减的衰减,所述初级辐射衰减即是当前投影,即当前辐射源位置的未被散射的辐射。为了考虑与仅源自辐射源焦点的初级辐射相比散射辐射入射角度的更宽范围,因此为了考虑衰减的不同路径,优选利用低通滤波器通过在初级衰减特征的局部区域上求平均值来计算衰减。由此,基于抗散射网格的角度传输函数指定所需低通滤波器的程度。
探测值校正单元还可以适于对探测值执行射束硬化校正,以进一步提高探测值的质量。如果成像装置为计算断层摄影装置,可以组合散射校正和射束硬化校正,以提供改进的计算断层摄影图像链,该图像链以物理上适当的方式应对散射并将散射效应与成像链中其他校正完全分开。
尽管在上述实施例中,成像装置优选是计算断层摄影扫描器,但在其他实施例中,成像装置可以是另一种成像模态,例如,另一种X射线成像装置或使用另一种辐射的成像装置,像诸如正电子发射断层摄影装置或单光子发射计算断层摄影装置的核成像装置。具体而言,该成像装置可以是基于例如C型臂的平面探测器锥形束计算断层摄影装置、辐射治疗中的板上成像机或核医学中的混合式锥形束计算断层摄影装置。
在成像装置中成像的对象优选是患者。不过,它也可以是像技术对象的另一个对象,例如在机场必须要检查的行李。
通过研究附图、说明书和权利要求,本领域技术人员能够在实践所主张的本发明的过程中理解并实施针对所公开的实施例的其他变型。
在权利要求中,“包括”一词不排除其他元件或步骤,不定冠词“一”不排除多个。
单个单元或设备可以实现权利要求中引用的若干项的功能。在互不相同的从属权利要求中引用某些措施这一事实不表示不能有利地采用这些措施的组合。
由一个或几个单元或设备执行的计算和确定,例如衰减元件散射值的计算和确定以及图像重建,可以由任意其他数量的单元或设备执行。例如,可以由单个单元或由任意其他数量的不同单元执行用于确定衰减元件散射值、用于校正空气扫描值和用于重建图像的计算。可以将根据上述成像方法和/或校正方法的计算和确定和/或成像装置和/或校正装置的控制实现为计算机程序的程序代码段和/或专用硬件。
可以将计算机程序存储/分布在适当的介质当中,例如,所述介质可以是光存储介质或者与其他硬件一起提供的或者作为其他硬件的部分的固态介质,但是,也可以使所述计算机程序通过其他形式分布,例如,通过因特网或者其他有线或无线电信***。
权利要求中的任何附图标记不应被解释为限制范围。
Claims (15)
1.一种用于产生对象(31)的感兴趣区域的图像的成像装置,所述成像装置包括:
-辐射源(2),其用于发射辐射(4),
-探测器(6),其用于在所述辐射(4)穿过所述感兴趣区域之后测量所述辐射(4)并且根据所测量的辐射(4)产生测量的探测值,
-衰减元件(30),其用于在所述辐射(4)穿过所述感兴趣区域之前衰减所述辐射(4),
-衰减元件散射值提供单元(12),其用于提供取决于所述衰减元件(30)引起的所述辐射(4)的散射的衰减元件散射值,
-探测值校正单元(17),其用于基于所提供的衰减元件散射值校正所测量的探测值,
-重建单元(18),其用于根据经校正的探测值重建所述感兴趣区域的图像。
2.根据权利要求1所述的成像装置,其中,所述成像装置包括要在其中放置所述对象(31)的检查区(5)以及用于准直所述辐射源(2)发射的所述辐射的准直器(3),所述准直器(3)位于所述衰减元件(30)和所述检查区(5)之间,其中,所述辐射源(2)、所述探测器(6)和所述准直器(3)适于针对所述辐射(4)的不同的准直测量校准值,并且其中,所述衰减元件散射值提供单元(12)适于基于所述校准值确定所述衰减元件散射值。
3.根据权利要求2所述的成像装置,其中,所述衰减元件散射值提供单元(12)适于通过将一函数拟合到所述校准值来确定所述衰减元件散射值,所述函数考虑了所述校准值的衰减元件散射部分和未散射部分并且取决于所述不同的准直。
4.根据权利要求3所述的成像装置,其中,所述函数对于描述所述校准值的衰减元件散射部分随准直变大而变大具有平方根相关性,并且对于描述所述校准值的未散射部分具有常数相关性,其中,所述衰减元件散射值取决于所述衰减元件散射部分。
5.根据权利要求2所述的成像装置,其中,
-所述探测器(6)包括若干探测器元件(50),
-所述辐射源(2)、所述探测器(6)和所述准直器(3)适于为所述探测器元件的子组(51)针对所述辐射(4)的不同的准直测量校准值,
-所述衰减元件散射值提供单元(12)适于针对所述探测器元件的子组(51)确定所述衰减元件散射值,
-所述成像装置还包括衰减元件散射值比提供单元(13),其用于提供探测器元件的所述子组(51)的衰减元件散射值与所述子组之外的探测器元件的衰减元件散射值的比值,
-所述衰减元件散射值提供单元(12)适于通过将针对所述子组(51)的所述衰减元件散射值乘以所提供的比值来确定针对所述子组之外的探测器元件的衰减元件散射值。
6.根据权利要求5所述的成像装置,其中,所述探测器(6)包括二维探测器表面,其中,探测器元件的所述子组(51)是所述探测器表面上的探测器元件的线。
7.根据权利要求1所述的成像装置,其中,所述成像装置还包括:
-空气扫描值提供单元(15),其用于提供空气扫描值,以及
-空气扫描值校正单元(16),其用于基于所述衰减元件散射值校正所述空气扫描值,
其中,所述探测值校正单元(17)适于基于经校正的空气扫描值校正所测量的探测值。
8.根据权利要求1所述的成像装置,其中,所述探测值校正单元(17)适于针对已经在所述衰减元件(30)中被散射并且被当在所述衰减元件(30)和所述探测器(6)之间存在对象(31)时的所述对象(31)衰减的辐射的衰减,调整所述衰减元件散射值,并且基于所调整的衰减元件散射值校正所测量的探测值。
9.根据权利要求8所述的成像装置,其中,通过将相应的衰减元件散射值乘以取决于所述对象(31)引起的衰减的相应的对象衰减值来执行所述衰减元件散射值的调整。
10.根据权利要求1所述的成像装置,其中,
-所述成像装置包括:
a)空气扫描值提供单元(15),其用于提供空气扫描值,以及
b)空气扫描值校正单元(16),其用于基于所述衰减元件散射值校正所述空气扫描值,
-所述探测值校正单元适于迭代地校正所测量的探测值,其中
a)在迭代步骤中,通过从除以相应的经校正的空气扫描值的相应的所测量的探测值中减去乘以相应的对象衰减值的相应的衰减元件散射值来计算所述经校正的探测值,所述相应的对象衰减值取决于由所述对象(31)引起的衰减,
b)一开始,所述相应的对象衰减值取决于所述相应的所测量的探测值除以所述相应的经校正的空气扫描值,并且在进一步迭代步骤中,所述对象衰减值取决于计算的相应的经校正的探测值。
11.一种用于校正测量的探测值的校正装置,为所述校正装置提供根据测量的辐射产生的探测值,所述辐射是由辐射源(2)发射的并且穿过衰减元件(30),所述衰减元件用于在所述辐射(4)穿过感兴趣区域之前衰减所述辐射(4),所述校正装置包括:
-衰减元件散射值提供单元(12),其用于提供取决于所述衰减元件(30)引起的所述辐射的散射的衰减元件散射值,
-探测值校正单元(17),其用于基于所提供的衰减元件散射值校正所测量的探测值。
12.一种用于对对象(31)的感兴趣区域进行成像的成像方法,所述成像方法包括:
-由辐射源(2)发射辐射(4),
-衰减元件(30)在所述辐射(4)穿过所述感兴趣区域之前衰减所述辐射(4),
-探测器(6)在所述辐射(4)穿过所述感兴趣区域之后测量所述辐射(4)并根据所测量的辐射(4)产生测量的探测值,
-衰减元件散射值提供单元(12)提供取决于所述衰减元件(30)引起的所述辐射(4)的散射的衰减元件散射值,
-探测值校正单元(17)基于所提供的衰减元件散射值校正所测量的探测值,
-重建单元(18)根据经校正的探测值重建所述感兴趣区域的图像。
13.一种用于校正测量的探测值的校正方法,所述校正方法包括以下步骤:
-提供根据测量的辐射产生的测量的探测值,所述辐射是由辐射源(2)发射的并且穿过衰减元件(30),所述衰减元件用于在所述辐射穿过感兴趣区域之前衰减所述辐射,
-衰减元件散射值提供单元(12)提供取决于所述衰减元件(30)引起的所述辐射的散射的衰减元件散射值,
-探测值校正单元(17)基于所提供的衰减元件散射值校正所述测量的探测值。
14.一种用于对对象(31)的感兴趣区域进行成像的成像计算机程序,所述成像计算机程序包括用于当所述成像计算机程序在控制根据权利要求1所述的成像装置的计算机上运行时,令所述成像装置执行根据权利要求12所述的成像方法的步骤的程序代码段。
15.一种用于校正测量的探测值的校正计算机程序,所述校正计算机程序包括用于当所述校正计算机程序在控制根据权利要求11所述的校正装置的计算机上运行时,令所述校正装置执行根据权利要求13所述的校正方法的步骤的程序代码段。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP08164401 | 2008-09-16 | ||
EP08164401.5 | 2008-09-16 | ||
PCT/IB2009/053936 WO2010032163A1 (en) | 2008-09-16 | 2009-09-09 | Imaging apparatus including correction unit for scattered radiation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102160085A true CN102160085A (zh) | 2011-08-17 |
CN102160085B CN102160085B (zh) | 2013-11-20 |
Family
ID=41514965
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009801360618A Active CN102160085B (zh) | 2008-09-16 | 2009-09-09 | 成像装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US8290116B2 (zh) |
CN (1) | CN102160085B (zh) |
WO (1) | WO2010032163A1 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105764422A (zh) * | 2014-10-01 | 2016-07-13 | 皇家飞利浦有限公司 | 成像装置和方法 |
CN107111888A (zh) * | 2015-01-07 | 2017-08-29 | 皇家飞利浦有限公司 | 用于校正能量依赖投影值的校正装置 |
CN108877892A (zh) * | 2017-05-12 | 2018-11-23 | 瓦里安医疗***公司 | 自动估算和减少计算机断层成像扫描中的散射 |
CN108931539A (zh) * | 2017-05-25 | 2018-12-04 | 北京君和信达科技有限公司 | 探测器自检方法、装置、介质及辐射型检查*** |
CN111936887A (zh) * | 2018-01-19 | 2020-11-13 | 皇家飞利浦有限公司 | 用于重建对象的感兴趣区域的谱图像的设备、***、方法和计算机程序 |
CN113453623A (zh) * | 2019-12-05 | 2021-09-28 | 皇家飞利浦有限公司 | 用于dax成像的全视场散射估计 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8712124B2 (en) * | 2011-06-21 | 2014-04-29 | General Electric Company | Artifact removal in nuclear images |
US10485498B2 (en) | 2016-01-29 | 2019-11-26 | Koninklijke Philips N.V. | Cone beam computed tomography projection values providing system and method |
WO2018082088A1 (zh) * | 2016-11-07 | 2018-05-11 | 深圳先进技术研究院 | 一种用于锥束ct图像散射修正的阻挡光栅优化方法及装置 |
WO2019128731A1 (en) * | 2017-12-29 | 2019-07-04 | Shenzhen United Imaging Healthcare Co., Ltd. | Systems and methods for scatter correction of image |
CN111096761B (zh) * | 2018-10-29 | 2024-03-08 | 上海西门子医疗器械有限公司 | 修正楔形滤波器散射的方法、装置和相关设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5394454A (en) * | 1992-05-09 | 1995-02-28 | U.S. Philips Corporation | Filter method for an x-ray system, and device for carrying out such a filter method |
CN1475000A (zh) * | 2000-10-13 | 2004-02-11 | ����˹�ش�ѧ | 使用圆加多弧轨迹进行圆锥射束体计算机断层摄影的***和方法 |
CN1818973A (zh) * | 2006-03-02 | 2006-08-16 | 西安交通大学 | 基于hl一致性条件的ct投影数据射束硬化效应校正方法 |
CN1954779A (zh) * | 2005-10-10 | 2007-05-02 | 西门子公司 | 用于ct***的辐射校正的方法 |
US20070165772A1 (en) * | 2006-01-14 | 2007-07-19 | General Electric Company | Methods and apparatus for scatter correction |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6183983A (ja) * | 1984-09-29 | 1986-04-28 | Shimadzu Corp | ポジトロンctの補正方法 |
US4809172A (en) | 1987-03-11 | 1989-02-28 | Kabushiki Kaisha Toshiba | Method for obtaining image data with a tomographic apparatus |
JPH0614911A (ja) * | 1992-07-02 | 1994-01-25 | Toshiba Corp | X線診断方法及びその装置 |
GB0116825D0 (en) * | 2001-07-10 | 2001-08-29 | Koninl Philips Electronics Nv | Determination of material parameters |
US6618466B1 (en) | 2002-02-21 | 2003-09-09 | University Of Rochester | Apparatus and method for x-ray scatter reduction and correction for fan beam CT and cone beam volume CT |
JP3942178B2 (ja) * | 2003-07-29 | 2007-07-11 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ctシステム |
US7477725B2 (en) * | 2003-11-11 | 2009-01-13 | Koninklijke Philips Electronics N.V. | Computed examination of an object by using coherent-scattered radiation |
US7065234B2 (en) | 2004-02-23 | 2006-06-20 | General Electric Company | Scatter and beam hardening correction in computed tomography applications |
DE102004060609A1 (de) * | 2004-12-16 | 2006-06-29 | Yxlon International Security Gmbh | Verfahren zum Messen des Impulsübertragungsspektrums von elastisch gestreuten Röntgenquanten |
JP2008525142A (ja) | 2004-12-29 | 2008-07-17 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | X線投影のアーチファクト補正のための装置及び方法 |
US7535987B2 (en) * | 2005-06-30 | 2009-05-19 | Kabushiki Kaisha Toshiba | X-ray CT apparatus |
US7463712B2 (en) * | 2006-05-18 | 2008-12-09 | The Board Of Trustees Of The Leland Stanford Junior University | Scatter correction for x-ray imaging using modulation of primary x-ray spatial spectrum |
WO2007148263A1 (en) | 2006-06-22 | 2007-12-27 | Philips Intellectual Property & Standards Gmbh | Method and system for error compensation |
-
2009
- 2009-09-09 CN CN2009801360618A patent/CN102160085B/zh active Active
- 2009-09-09 WO PCT/IB2009/053936 patent/WO2010032163A1/en active Application Filing
- 2009-09-09 US US13/063,486 patent/US8290116B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5394454A (en) * | 1992-05-09 | 1995-02-28 | U.S. Philips Corporation | Filter method for an x-ray system, and device for carrying out such a filter method |
CN1475000A (zh) * | 2000-10-13 | 2004-02-11 | ����˹�ش�ѧ | 使用圆加多弧轨迹进行圆锥射束体计算机断层摄影的***和方法 |
CN1954779A (zh) * | 2005-10-10 | 2007-05-02 | 西门子公司 | 用于ct***的辐射校正的方法 |
US20070165772A1 (en) * | 2006-01-14 | 2007-07-19 | General Electric Company | Methods and apparatus for scatter correction |
CN1818973A (zh) * | 2006-03-02 | 2006-08-16 | 西安交通大学 | 基于hl一致性条件的ct投影数据射束硬化效应校正方法 |
Non-Patent Citations (3)
Title |
---|
HABIB ZAIDI.MOHAMMAD REZA AY: "Current status and new horizons in Monte Carlo simulation of X-ray CT scanners", 《MED BIO ENG COMPUT》 * |
MOHAMMAD REZA AY AND HABIB ZAIDI: "Development and validation of MCNP4C-based Monte Carlo simulator for fan-and cone-beam x-ray CT", 《PHYSICS IN MEDICINE AND BIOLOGY》 * |
WOJCIECH ZBIJEWSKI AND FREEK J.BEEKMAN: "Efficient Monte Carlo based scatter artifact reduction in cone-beam micro-CT", 《IEEE TRANSACTIONS ON MEDICAL IMAGING》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105764422A (zh) * | 2014-10-01 | 2016-07-13 | 皇家飞利浦有限公司 | 成像装置和方法 |
CN105764422B (zh) * | 2014-10-01 | 2019-01-29 | 皇家飞利浦有限公司 | 成像装置和方法 |
CN107111888A (zh) * | 2015-01-07 | 2017-08-29 | 皇家飞利浦有限公司 | 用于校正能量依赖投影值的校正装置 |
CN108877892A (zh) * | 2017-05-12 | 2018-11-23 | 瓦里安医疗***公司 | 自动估算和减少计算机断层成像扫描中的散射 |
CN108877892B (zh) * | 2017-05-12 | 2022-03-22 | 瓦里安医疗***公司 | 自动估算和减少计算机断层成像扫描中的散射 |
CN108931539A (zh) * | 2017-05-25 | 2018-12-04 | 北京君和信达科技有限公司 | 探测器自检方法、装置、介质及辐射型检查*** |
CN111936887A (zh) * | 2018-01-19 | 2020-11-13 | 皇家飞利浦有限公司 | 用于重建对象的感兴趣区域的谱图像的设备、***、方法和计算机程序 |
CN113453623A (zh) * | 2019-12-05 | 2021-09-28 | 皇家飞利浦有限公司 | 用于dax成像的全视场散射估计 |
Also Published As
Publication number | Publication date |
---|---|
US8290116B2 (en) | 2012-10-16 |
CN102160085B (zh) | 2013-11-20 |
WO2010032163A1 (en) | 2010-03-25 |
US20110164722A1 (en) | 2011-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102160085B (zh) | 成像装置 | |
CN103400401B (zh) | 改善图像分辨率的***和方法 | |
CN101810488B (zh) | 探测值校正设备 | |
CN103329168B (zh) | 针对spect/ct***的迭代锥形射束ct重建的截断补偿 | |
US8000513B2 (en) | System and method for 3D time of flight PET forward projection based on an exact axial inverse rebinning relation in fourier space | |
US6339652B1 (en) | Source-assisted attenuation correction for emission computed tomography | |
US8008625B2 (en) | Method and apparatus for high-sensitivity single-photon emission computed tomography | |
CN101262819B (zh) | Ct散射的直接测量及校正 | |
CN101489486B (zh) | 基于散射辐射分数的x-射线探测器增益校准 | |
US20030138074A1 (en) | Method for correcting stray radiation in an x-ray computed tomograph scanner | |
US8712134B2 (en) | Method and system for expanding axial coverage in iterative reconstruction in computer tomography (CT) | |
CN104334081A (zh) | X射线ct成像器的运动层分解校准 | |
CN100358471C (zh) | 放射计算机断层摄像装置及断层图像数据生成方法 | |
CN103596502A (zh) | 针对断层摄影成像***的自适应校准 | |
CN102177430A (zh) | 使用经移位的几何结构改进ct图像采集的方法和设备 | |
US7737972B2 (en) | Systems and methods for digital volumetric laminar tomography | |
CN103229211A (zh) | 前向投影设备 | |
US8389943B2 (en) | Modeling of the point-spread-function in single-pinhole and multi-pinhole spect reconstruction | |
US8351572B1 (en) | Method and system to reconstruct treatment dose to a patient from integrated exit-transit images of radiation fields taken during treatment | |
US6921902B2 (en) | Scatter correction device for radiative tomographic scanner | |
Ordoñez et al. | Fast in situ image reconstruction for proton radiography | |
US7385200B2 (en) | Re-binning method for nuclear medicine imaging devices | |
CN101291627A (zh) | 先进的csct检测器形状 | |
Zeng et al. | Image reconstruction algorithm for a rotating slat collimator | |
EP3731760B1 (en) | Processing pipeline for prompt particle image reconstruction |
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 |