CN102440801A - 校准具有至少一个焦点-探测器组合的ct 系统的方法 - Google Patents

校准具有至少一个焦点-探测器组合的ct 系统的方法 Download PDF

Info

Publication number
CN102440801A
CN102440801A CN2011103093517A CN201110309351A CN102440801A CN 102440801 A CN102440801 A CN 102440801A CN 2011103093517 A CN2011103093517 A CN 2011103093517A CN 201110309351 A CN201110309351 A CN 201110309351A CN 102440801 A CN102440801 A CN 102440801A
Authority
CN
China
Prior art keywords
detector
die body
focus
value
described method
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
Application number
CN2011103093517A
Other languages
English (en)
Other versions
CN102440801B (zh
Inventor
S.卡普勒
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Medical Ag
Original Assignee
Siemens AG
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Siemens AG filed Critical Siemens AG
Publication of CN102440801A publication Critical patent/CN102440801A/zh
Application granted granted Critical
Publication of CN102440801B publication Critical patent/CN102440801B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
    • A61B6/4241—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
    • 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/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • 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/58—Testing, adjusting or calibrating thereof
    • A61B6/582—Calibration
    • A61B6/583—Calibration using calibration phantoms
    • 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/58—Testing, adjusting or calibrating thereof
    • A61B6/582—Calibration
    • A61B6/585—Calibration of detector units
    • 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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4007—Arrangements for generating radiation specially adapted for radiation diagnosis characterised by using a plurality of source units
    • A61B6/4014—Arrangements for generating radiation specially adapted for radiation diagnosis characterised by using a plurality of source units arranged in multiple source-detector units
    • 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
    • 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/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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Pulmonology (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明涉及一种用于校准具有至少一个焦点-探测器组合(2，3)的CT系统(C1)的方法，该焦点-探测器组合具有量子计数的探测器(D1)，该量子计数的探测器具有多个探测器元件，其中，该焦点-探测器组合(2，3)被围绕测量区域(11)和其中布置的系统轴(9)可旋转地设置，并且从焦点(F1)向探测器发送X射线束(S1)，该X射线束具有在能量范围上的X射线能量谱，其中，即使在量子计数的探测器的瘫痪范围中也将所获得的CT扫描的实际衰减值与理论上的单能量的额定衰减值比较并且对于每个探测器元件确定在额定衰减值和实际衰减值之间的变换函数，以及由此对所测量的衰减值进行校准。

Description

校准具有至少一个焦点-探测器组合的CT 系统的方法
技术领域
本发明涉及一种用于借助对至少一个模体的扫描来校准具有至少一个焦点-探测器组合的CT系统的方法，该焦点-探测器组合具有量子计数的探测器，该量子计数的探测器具有多个探测器元件，其中该焦点-探测器组合围绕测量区域和其中布置的系统轴可旋转地设置，并且从焦点向探测器发送X射线束，该X射线束具有在能量范围上的X射线能量谱，其中，将由此所确定的实际衰减值与额定衰减值进行比较，并且在随后对患者的扫描中对所确定的衰减值进行校正，以及利用校正的衰减值进行CT图像数据的重建并且输出和/或继续处理CT图像数据。
背景技术
本文开头提到种类的校准方法在计算机断层造影系统(CT)方面是一般公知的。然而，在该公知的校准方法中，探测器不是在落入探测器的饱和范围中的测量范围中被校准、而是在如下区域中被校准，在该区域中，假定在击中探测器的射线强度和信号应答之间是线性关系。量子计数的探测器位于X射线诊断地利用的射线强度的范围，然而在CT检查中至少部分地位于强烈非线性的区域中，其中此处已经出现探测器中的瘫痪迹象。
发明内容
本发明要解决的技术问题是，找到一种改进的校准方法，该方法既可以在非饱和的区域中也可以在出现瘫痪的区域中提供令人满意的结果。
发明人认识到以下：
由于在临床计算机断层造影(＝CT)中出现的高的X射线量子流和由此所产生的探测器信号的瘫痪，该瘫痪导致在探测器信号和射线强度之间的非线性，在将来的量子计数的CT系统中必须实现相当复杂的校准流程。探测器信号的非线性，确切说是其计数率在此受到各种影响，诸如：
·到达探测器上的X射线谱；
·到达探测器上的X射线流；
·有效的像素大小；
·计数器的触发阈值的精确值；
·由探测器所提供的、通过感应产生的脉冲的特性和有效持续时间；
·以及电子电路的信号形状特性。
目前作为对于临床CT探测器的候选尚待解决的直接转换器-探测器材料，诸如CdTe、CZT等，就其电荷收集效率(Ladungssammlungseffizienz)和脉冲特性来说具有突出的相当弱的均匀性。这意味着，这样的探测器的未校准的信号是很不均匀的并且大多要求对各个探测器像素或探测器通道进行开销大的个别校准。
校准方法(如在具有能量集成的经典的CT系统所使用的那样)只有在没有超过一定的临界量子流时才可以应用于量子计数的系统。然而，该临界量子流通常低于在临床CT中出现的量子流。目前没有公知用于校准具有量子计数的临床CT系统的方法。
按照本发明，建议了一种方法，该方法既包括探测器信号的线性化(Linearisierung)也包括射束硬化校正。可以理解的是，在本发明的范围内也可以分开地应用两个子方面。在事先调整了单个的探测器像素或探测器通道的计数器的触发阈值之后，应用以下描述的按照本发明的方法。
相应于按照本发明的方法，将不同大小的特定材料、典型地是水的模体，例如圆柱体、椭圆、直角平行六面体等，优选是偏心地这样定位在CT系统的测量区域中并且扫描，使得每个探测器通道在一个完整运转期间覆盖(überstreicht)在最大和最小衰减之间的尽可能完整的区域。此外，进行开敞扫描(Freiluft-Scans)。对模体的所述扫描优选在恒定的管电流(mA)的情况下进行，校准具有对于该管电流的随后的有效性。开敞扫描既可以在相同的也可以在不同的mA值的情况下进行。管电压在所有的扫描中应该选择为相同。在此，这里所建议的方法的主要核心在于在以下描述的对这里所获得的数据的分析并且特别地在于对对应的额定衰减值的确定。
在此，首先将模体的所测量的强度数据I借助来自于开敞扫描的强度数据I0对于每个通道分开地换算为衰减数据A：
A＝-ln(I)+ln(I0)+c，
其中，c是可自由选择的参数。将该数据A与基于单色X射线谱的额定衰减值a比较。该额定衰减值可以按照后面将详细示出的不同方式获得。随后，对于每个通道分开地应用以下流程：
1.将所有模体扫描的所有投影的所测量的衰减数据A映射到矢量A中。
2.另一个矢量a包含相应于A的单个分量的额定衰减数据。
3.在散布图(＝Scatter-Plot)中，相对额定值a绘制A的测量的值。使用拟合或回归算法，以按照解析关系映射数据：a(A，p)，其中p是匹配的参数。有意义地进行对于曲线变化优化的、借助多项式进行的逐段的匹配。
最后，函数a(A，p)代表探测器通道的校准规则A→a(A，p)，借助该校准规则保证了探测器数据的均匀性和线性化(包括射束硬化校正)。
额定衰减值的确定可以按照以下不同的方式进行：
-从正弦图拟合(Sinogramm-Fit)中确定额定衰减数据，作为理论的衰减数据。借助于所测量的正弦图，从模体扫描中确定模体与对称中心的距离，以及模体的方位角。然后，解析地计算基于单能量的X射线谱的额定衰减值。
-重建衰减数据A。通过在断层造影的图像数据中的阈值形成和分割，利用然后对模体与对称中心的距离以及方位角的确定，来识别模体。然后解析地计算基于单能量的X射线谱的额定衰减值。
-重建衰减数据A。通过在断层造影的图像数据中的阈值形成和分割识别模体。修改图像数据，例如将作为模体分割的像素置为0HU并且位于模体外部的像素置于-1000HU。在考虑精确的扫描器几何的条件下进行图像数据的前向投影，其中来自于前向投影的数据代表了基于单能量的X射线谱的额定衰减值。
该方法的其他版本例如是：
-迭代地执行用于借助拟合或前向投影来改进位置确定的上述方法；
-多次运行管电流的调制，以便对于在一个校准过程(Kalibrierungs-durchgang)内部不同的管电流产生校准规则；
-在机架不旋转的静态扫描中使用离心地布置的圆柱形模体或者横梁形模体；
-在旋转扫描中，在对管电流而不是离心模体进行调制的情况下进行开敞扫描(Freiluft-Scans)，以便仅进行探测器信号的线性化而不进行射束硬化校正；
-借助有理函数或指数函数进行逐段的匹配；
-借助线性内插或(B)样条内插逐段地内插平均的数据点。
利用上面描述的方法，示出了即使在最高的管电流和探测器信号部分瘫痪的情况下，就探测器信号的线性度和射束硬化校正方面，具有量子计数的探测器的临床CT系统的有效和高质量的校准。
此外，还可以将按照本发明的方法应用于构造为双源CT系统并且既具有常规的探测器系统又具有量子计数的探测器系统的CT系统中。如果常规的探测器已经是校准的，则可以将该常规的探测器用来改善对量子计数的探测器的数据的校准的质量。在此的基础是常规的探测器系统的有效的校准以及通过公知的对齐方法(Alignment-Verfahren)确定的、常规的探测器系统的数据网格(Datenraster)到量子计数的探测器系统的数据网格的几何变换规则。
在这些前提条件下，可以通过如下来优化上面描述的对额定衰减值的确定，即，将关于线性度和射线硬化完全校准的、集成的探测器系统的投影数据变换到量子计数的探测器系统的网格(Raster)并且直接用作额定衰减数据。
此外，在考虑探测器网格(Detektorraster)的变换的条件下，从集成的探测器系统的数据的正弦图拟合，确定对于量子计数的探测器系统的额定衰减值，作为理论的衰减数据。
替换地，可以通过在集成的探测器系统的断层造影图像数据中的阈值形成和分割来确定模体的位置，特别是通过对模体与CT系统的对称中心的距离的确定和对方位角的确定。然后，基于单能量的X射线谱解析地利用这样已知的模体几何形状和位置来计算对于量子计数的探测器系统的额定衰减值。
此外，可以通过在集成的探测器系统的断层造影图像数据中的阈值形成来定位并且分割模体，方法是，将模体内部的断层造影图像数据置为第一预定的值，例如0HU，并且将模体外部的图像数据置为第二预定的值，例如1000HU。然后，在关于量子计数的探测器系统考虑精确的扫描器几何形状的条件下进行图像数据的前向投影。来自于前向投影的数据此时直接代表对于基于单能量的X射线谱的量子计数的探测器系统的额定衰减值。
如果集成的探测器系统具有带有比量子计数的探测器系统更大的视野的探测器，则在重建时利用来自于集成的探测器系统的数据来补充量子计数的探测器系统的数据，以便保证上述方法的更好的收敛性。
因此，该方法描述了一种对具有量子计数的探测器的临床CT系统的有效和高质量的校准，特别是在既具有量子计数又具有常规的探测器的混合双源CT系统的情况下还有对量子计数的探测器的校准。在此，按照本发明的校准方法既涉及探测器信号的线性度又涉及射束硬化校正，即使在最高的管电流和探测器信号的部分瘫痪的情况下。
相应于上述发明思路，发明人由此建议一种用于校准具有至少一个焦点-探测器组合的CT系统的方法，该焦点-探测器组合具有量子计数的探测器，该探测器具有多个探测器元件，其中，焦点-探测器组合围绕测量区域和其中布置的系统轴可旋转地布置，并且从焦点向探测器发送X射线束，该X射线束具有在能量区域上的X射线能量谱，其中，按照本发明的方法具有以下方法步骤：
-至少一次性地将至少一个模体这样布置在测量区域中，使得围绕系统轴一次运转的焦点-探测器组合的每个测量射束至少一次穿过最大衰减的区域并且一次穿过最小衰减的区域，
-从环绕地(umlaufend)布置的多个投影角至少一次性地扫描至少一个模体并且关于多个剂量功率值(Dosisleistungswerten)进行没有模体的环绕扫描，
-对每个探测器元件，形成一个在空间上彼此相同的射线的、从量子计数的探测器的测量值中所计算的实际衰减值和额定衰减值之间的散布图，其中额定衰减值涉及预定的奇异的(
Figure BDA0000098341430000051
)射线能量，
-对于每个探测器元件，确定在额定衰减值和实际衰减值之间的至少一个变换函数，
-利用具有量子计数的探测器的至少一个焦点-探测器组合扫描患者并输出探测器测量值，利用探测器测量值计算衰减值，对每个探测器元件借助至少一个变换函数校正计算的衰减值，利用校正了的衰减值重建CT图像数据并且输出和/或进一步处理CT图像数据。
在此要提到的是，用来进行模体的扫描和空气测量的多个剂量功率值，优选地应该位于在以后在临床操作中实际上发生的X射线管的剂量功率调整的范围中和后面的在探测器上变得有效的射线强度中。
优选地，可以通过散布图的值的至少一个曲线匹配计算变换函数，其中曲线匹配也可以逐段地在散布图中进行。
按照一种实施变形，可以直接从模体在测量区域中的预先已知的位置和模体材料的预先已知的特定于能量的吸收系数，包括焦点-探测器系统的预先已知的空间结构，计算额定衰减值。
此外，还可以从正弦图拟合中确定额定衰减值，作为理论上的衰减数据，方法是：
-借助所测量的正弦图从模体扫描中确定模体与对称中心的距离以及模体的方位角，并且
-基于这样所确定的定位和模体的已知的几何形状，在假定单能量的X射线谱的条件下解析地计算额定衰减值。
替换地，存在如下可能性：通过从至少一个投影角通过至少一个焦点-探测器组合对模体的先前的扫描确定至少一个模体的位置。这简化了模体在测量区域中的定位，因为这里定位本身可以不太严格和精确地进行。为了确定位置，在此可以使用一种用于分割和识别在断层造影图像数据中的至少一个模体的基于阈值的方法。
此外，还可以通过如下确定至少一个模体的额定衰减值，即，利用CT扫描并且进行至少一个模体的重建并且确定CT图像数据，在CT图像数据中将具有高于预先给出的阈值的CT值的所有像素/体素利用预先给出的恒定的CT值来代替并且从这些改变了的CT图像数据中通过前向投影计算额定衰减值。
在所有上述方法变形中具有优势的是，将所使用的不同大小的模体相对于系统轴离心地布置。模体可以构造为球形的或圆柱形的或椭圆形的。
还提出，对至少一个模体利用不同的剂量功率进行多次扫描，以便在一次过程(Durchgang)之内产生对于不同的剂量功率值的校准数据。
此外，优选地还可以对至少一个模体利用在至少一个焦点-探测器组合的一次旋转期间改变的剂量功率进行至少一次扫描和/或在没有模体的情况下利用不同的剂量功率进行多次扫描。此外，还可以在没有模体的情况下利用在至少一个焦点-探测器组合的一次旋转期间改变的剂量功率进行至少一次扫描。
还提出，为了确定至少一个变换函数，特别是逐段地执行以下提到的函数拟合(Funktionsfit)中的一个：
-样条内插，
-B样条内插，
-拟合或匹配，优选为指数函数的最小平方拟合，
-拟合或匹配，有理函数的最小平方拟合。
还可以关于其中同时使用至少两个焦点-探测器系统的CT系统执行上面描述的校准方法，其中设置具有集成的探测器元件的至少一个常规的探测器。特别地在这种情况下，为了确定额定衰减值至少还可以使用常规的探测器的探测器数据，其中至少一个常规的探测器应该是已经校准了的。
在此还具有优势的是，具有常规的探测器的至少一个焦点-探测器系统和具有量子计数的探测器的至少一个焦点-探测器系统以相同的X射线谱工作。
同样具有优势的是，具有常规的探测器的至少一个焦点-探测器系统和具有量子计数的探测器的至少一个焦点-探测器系统在利用不同的几何网格(＝几何的数据网格)扫描的情况下工作，并且为了校准将数据从常规的探测器的网格传输到量子计数的探测器的网格或者反过来。例如可以通过内插来将来自于常规的探测器系统的几何网格的数据传输到量子计数的探测器的网格。
附图说明
以下借助附图结合优选实施例详细描述本发明，其中仅示出对于理解本发明所需的特征。使用以下附图标记：1：CT设备；2：第一X射线管；3：第一检测器；4：第二X射线管(可选的)；5：第二检测器(可选的)；6：支架壳体；7：模体；8：患者卧榻；9：系统轴；10：计算机；11：测量区域；12：散布图；a：额定衰减值；A：实际衰减值；D1：量子计数的探测器；D2：集成的探测器；F1，F2：焦点；P1至P3：模体；Prg1至Prgn：计算机程序；S1，S2：射线束；Sp1，Sp2：正弦图。
具体地：
图1示出了CT系统，
图2示出了具有三个模体变形的CT系统的焦点-探测器系统，
图3示出了对于两个不同大小的圆柱形模体(上边的小，下边的大)具有测量的实际衰减值(左边)和额定衰减值(右边)的CT系统的正弦图，
图4示出了根据相对的额定和实际衰减值的散布图，
图5示出了模体的重建的CT截面图(左边)和从中通过阈值分析确定的并且与模体的已知几何形状匹配的截面图(右边)，并且
图6示出了具有集成的和量子计数的探测器的混合CT系统的两个焦点-探测器系统。
具体实施方式
图1示出了待校准的CT系统1的例子，其中位于可移动的患者卧榻8上的模体7沿着系统轴9被移动通过机架壳体6中的开口，而此处未详细示出的机架围绕模体7旋转地运动。位于机架上的是至少一个由第一X射线管2和对置的量子计数的探测器3组成的管探测器系统。可选地，还可以在相同的机架上布置具有第二X射线管4和第二探测器5的第二管探测器系统。该第二探测器5可以是第二量子计数的探测器或者可以使用常规的集成的探测器，其相应于上面描述的方法可以简化量子计数的探测器的校准。对模体7的扫描通过计算机系统10来控制，其中探测器数据的分析和校准计算也可以通过该计算机系统10来进行。
相应地，该计算机系统10还具有程序Prg1至Prgn，利用这些程序可以在CT系统1的运行中执行按照本发明的方法。当然还可以通过这些程序进行随后对患者的扫描、图像重建和图像输出或图像数据的进一步处理，其中由探测器系统接收的测量数据被换算为吸收数据并且其相应于找到的变换函数为了校准的目的而被校正。
图2示出了类似于图1的CT系统的简单的焦点-探测器系统，包括三个具有横截面中不同直径的相对于系统轴9离心地布置的圆柱形模体P1至P3。焦点-探测器系统包括焦点F1，成束的射线扇形或射线锥从该焦点延伸到相对布置的量子计数的探测器D1。用虚线示出了最大可用的测量区域(＝FOV＝Field Of View，视野)11。焦点-探测器系统实施为可绕系统轴9旋转，从而在对模体P1至P3的依次执行的旋转扫描时从实际衰减值中建立正弦图，如在图3中举例示出的。
在图3的左列示出了由测量的实际衰减值得到的两个这样的正弦图SP1和SP2，其是由具有图2的模体P1和P2的CT系统产生的。同时在右边与之并列布置地绘出了具有额定衰减值的对应的正弦图，其来源于根据模体材料的代表性射线能量的已知的吸收值进行的理论计算。作为代表性射线能量，例如可以使用基于对正弦图SP1和SP2的测量所使用的射线谱而加权平均的射线能量。
图4示出了利用这些额定和实际衰减值的数据所产生的散布图12。由此相对于在横坐标上的利用给出的X射线能量谱记录的相应的实际衰减值A，在纵坐标上绘出了单色的额定衰减值a。在此要指出的是，在此处示出的测量值的具体情况下，在量子计数的探测器中的射线强度是借助所谓的在背负式(Pickaback)触发电路记录的，该电路已经补偿了在高的测量区域中的测量情况下发生的瘫痪现象。该电路和该测量方法在文件号为DE 10 2009 055 807.1的专利申请文件中有详细描述。该文献的公开内容应该全部合并到本申请书中。
在图4中示出的在衰减的额定和实际值之间的相关性(Korrelation)此时借助函数匹配(Funktionsanpassungen)、必要时还利用关于射线强度逐段的函数匹配按照简单的公式来表达(niedergelegt)。然而还存在如下可能性，即，建立在测量的区域上具有平均的实际衰减值对额定衰减值的“查询(look up)”表，而不是函数匹配，必要时内插中间值，并且利用该“查询”表进行校准。
图5左边示出了模体的重建后的CT截面图，而在右边示出了从中通过阈值分析确定的并且与模体的已知几何形状(此处是圆形)匹配的截面图。如果模体的位置不是要通过对模体在CT的测量区域中的实际的测量来进行，而是要对于额定衰减值的进一步的理论计算借助计算机断层造影拍摄来确定其位置，则可以使用该方法。
此外，在图6中还示出混合CT系统的双焦点-探测器系统的实施。第一探测器系统包括将其射线束S1投射到第一量子计数的探测器D1上的第一焦点F1。第二焦点-探测器系统包括将其射线束S2投射到第二集成的探测器D2上的第二焦点F2。两个焦点-探测器系统在该例子中以90°可旋转地围绕系统轴9偏移地布置。在本例中，集成的探测器还具有探测器元件的比量子计数的探测器更窄网眼的(engmaschigeres)网格，其还具有更窄的扇形角。
如上所述，现在可以借助已经校准的集成的探测器D2和对模体P1至P3的扫描进行对量子计数的探测器D1的校准，其中如上所述必须考虑两个探测器的不同的测量网格。
总之，通过本发明因此描述了用于校准CT系统的方法，其中，即使在量子计数的探测器的瘫痪范围中也将所获得的CT扫描的实际衰减值与理论上的单能量的额定衰减值进行比较，并且对于每个探测器元件确定在额定衰减值和实际衰减值之间的变换函数，以及由此对所测量的衰减值进行校准。
可以理解的是，前面所提到的本发明特征不仅可以按照分别给出的组合，而且可以按照其他组合或独立使用，而不脱离本发明的范围。

Claims (22)

1.一种用于校准具有至少一个焦点-探测器组合(2，3)的CT系统(C1)的方法，该焦点-探测器组合具有量子计数的探测器(D1)，该量子计数的探测器具有多个探测器元件，其中，该焦点-探测器组合(2，3)围绕测量区域(11)和其中布置的系统轴(9)可旋转地设置，并且从焦点(F1)向探测器发送X射线束(S1)，该X射线束具有在能量范围上的X射线能量谱，该方法具有以下方法步骤：
1.1.至少一次性地将至少一个模体(7，P1，P2，P3)这样布置在测量区域(11)中，使得围绕系统轴(9)一次运转的焦点-探测器组合(2，3)的每个测量射束至少一次穿过最大衰减的区域并且一次穿过最小衰减的区域，
1.2.从环绕地布置的多个投影角至少一次扫描至少一个模体(7，P1，P2，P3)，并且关于多个剂量功率值进行没有模体的环绕扫描，
1.3.对每个探测器元件和剂量功率值，形成在空间上彼此相同的射线的、从量子计数的探测器的测量值中所计算的实际衰减值和额定衰减值之间的散布图(12)，其中，所述额定衰减值涉及预定的奇异的射线能量，
1.4.对于每个探测器元件确定在额定衰减值和实际衰减值之间的至少一个变换函数，
1.5.利用具有量子计数的探测器(D1)的至少一个焦点-探测器组合(2，3)扫描患者并输出探测器测量值，利用所述探测器测量值计算衰减值，对每个探测器元件借助至少一个变换函数校正所计算的衰减值，利用校正了的衰减值重建CT图像数据并且输出和/或进一步处理CT图像数据。
2.根据权利要求1所述的方法，其特征在于，通过散布图(12)的值的至少一个曲线匹配计算所述变换函数。
3.根据权利要求2所述的方法，其特征在于，在散布图中逐段地进行所述曲线匹配。
4.根据权利要求1至3中任一项所述的方法，其特征在于，根据模体在测量区域中的预先已知的位置和模体材料的预先已知的特定于能量的吸收系数，包括焦点-探测器系统(2，3)的预先已知的空间结构，计算所述额定衰减值。
5.根据权利要求1至3中任一项所述的方法，其特征在于，从正弦图拟合中确定额定衰减值，作为理论上的衰减数据，方法是：
-借助所测量的正弦图(SP1，SP2)从模体扫描中确定模体(7，P1，P2，P3)与对称中心的距离以及模体(7，P1，P2，P3)的方位角，并且
-基于这样所确定的定位和模体(7，P1，P2，P3)的已知的几何形状，在假定单能量的X射线谱的条件下解析地计算所述额定衰减值。
6.根据权利要求1至3中任一项所述的方法，其特征在于，通过从至少一个投影角通过至少一个焦点-探测器组合(2，3)对模体(7，P1，P2，P3)的先前的扫描，确定至少一个模体(7，P1，P2，P3)的位置。
7.根据权利要求1至3中任一项所述的方法，其特征在于，使用用于确定所述至少一个模体(7，P1，P2，P3)在断层造影图像数据中的位置的基于阈值的方法，并且利用这样所确定的位置信息计算所述额定衰减值。
8.根据权利要求1至3中任一项所述的方法，其特征在于，通过如下确定至少一个模体(7，P1，P2，P3)的额定衰减值：利用CT扫描并且进行至少一个模体的重建并且确定CT图像数据，在所述CT图像数据中将具有高于预先给出的阈值的CT值的所有像素/体素利用预先给出的恒定的CT值来代替，并且从这些改变了的CT图像数据中通过前向投影计算所述额定衰减值。
9.根据权利要求1至8中任一项所述的方法，其特征在于，使用与系统轴(9)离心地布置的模体(7，P1，P2，P3)作为模体(7，P1，P2，P3)。
10.根据权利要求1至9中任一项所述的方法，其特征在于，还使用至少一个横梁形模体作为模体(7，P1，P2，P3)。
11.根据权利要求1至10中任一项所述的方法，其特征在于，还使用至少一个圆柱形模体模体作为模体(7，P1，P2，P3)。
12.根据权利要求1至11中任一项所述的方法，其特征在于，对至少一个模体(7，P1，P2，P3)利用不同的剂量功率进行多次扫描。
13.根据权利要求1至12中任一项所述的方法，其特征在于，对至少一个模体(7，P1，P2，P3)利用在至少一个焦点-探测器组合(2，3)的一次旋转期间改变的剂量功率进行至少一次扫描。
14.根据权利要求1至13中任一项所述的方法，其特征在于，在没有模体(7，P1，P2，P3)的情况下利用不同的剂量功率进行多次扫描。
15.根据权利要求1至14中任一项所述的方法，其特征在于，在没有模体的情况下利用在至少一个焦点-探测器组合(2，3)的一次旋转期间改变的剂量功率进行至少一次扫描。
16.根据权利要求1至15中任一项所述的方法，其特征在于，为了确定至少一个变换函数，执行以下提到的函数拟合中的一个、特别是逐段地执行：
-样条内插，
-B样条内插，
-指数函数的最小平方拟合，
-有理函数的最小平方拟合。
17.根据权利要求1至16中任一项所述的方法，其特征在于，在CT系统(1)中同时使用至少两个焦点-探测器系统(2，3；4，5)，其中设置具有集成的探测器元件的至少一个常规的探测器(D2)。
18.根据权利要求1至17中任一项所述的方法，其特征在于，为了确定额定衰减值，至少还使用所述常规的探测器(D2)的探测器数据。
19.根据权利要求17或18所述的方法，其特征在于，所述至少一个常规的探测器(D2)是已经校准了的。
20.根据权利要求17至19中任一项所述的方法，其特征在于，具有常规的探测器(D2)的至少一个焦点-探测器系统(4，5)和具有量子计数的探测器(D1)的至少一个焦点探-测器系统(2，3)以相同的X射线谱工作。
21.根据权利要求17至20中任一项所述的方法，其特征在于，具有常规的探测器(D2)的至少一个焦点-探测器系统(4，5)和具有量子计数的探测器(D1)的至少一个焦点-探测器系统(2，3)在利用不同的几何网格(＝几何的数据网格)扫描的情况下工作，并且为了校准将数据从常规的探测器(D2)的网格传输到量子计数的探测器(D1)的网格。
22.根据权利要求21所述的方法，其特征在于，将来自于常规的探测器系统(D2)的几何网格的数据通过内插传输到量子计数的探测器(D1)的网格。
CN201110309351.7A 2010-10-13 2011-10-13 校准具有至少一个焦点-探测器组合的ct 系统的方法 Active CN102440801B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102010042388A DE102010042388A1 (de) 2010-10-13 2010-10-13 Verfahren zur Kalibrierung eines CT-Systems mit zumindest einer Fokus-Detektor-Kombination mit einem quantenzählenden Detektor
DE102010042388.2 2010-10-13

Publications (2)

Publication Number Publication Date
CN102440801A true CN102440801A (zh) 2012-05-09
CN102440801B CN102440801B (zh) 2014-05-14

Family

ID=45895604

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110309351.7A Active CN102440801B (zh) 2010-10-13 2011-10-13 校准具有至少一个焦点-探测器组合的ct 系统的方法

Country Status (3)

Country Link
US (1) US8628241B2 (zh)
CN (1) CN102440801B (zh)
DE (1) DE102010042388A1 (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103340643A (zh) * 2013-06-09 2013-10-09 成都军区昆明总医院 应用多功能标准体模对ct关键技术指标进行检测的方法
CN104111469A (zh) * 2013-04-17 2014-10-22 上海联影医疗科技有限公司 模拟x射线穿透模体组件的衰减强度的方法
CN104603606A (zh) * 2012-09-24 2015-05-06 西门子公司 用于确定通过待检查的对象引起的x射线辐射的衰减的方法和装置
CN107693037A (zh) * 2013-09-18 2018-02-16 上海联影医疗科技有限公司 Pet扫描装置及其时间偏移校正方法
CN109416405A (zh) * 2016-06-29 2019-03-01 皇家飞利浦有限公司 光子计数计算机断层摄影
CN113796879A (zh) * 2021-09-27 2021-12-17 北京万东医疗科技股份有限公司 一种球管出射能谱验证方法、装置、电子设备及存储介质

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140072108A1 (en) * 2010-07-16 2014-03-13 David P. Rohler Methods and apparatus for extended low contrast detectability for radiographic imaging systems
DE102012205290A1 (de) * 2012-03-30 2013-10-02 Siemens Aktiengesellschaft Verfahren und CT-System zur Artefaktreduktion bei computertomographischen Aufnahmen mit zählendem Detektor
US20130301799A1 (en) * 2012-05-14 2013-11-14 Samsung Electronics Co., Ltd. X-ray imaging apparatus and control method therefor
DE102012216272A1 (de) * 2012-09-13 2014-03-13 Siemens Aktiengesellschaft Röntgenfokusjustage
JP6691734B2 (ja) * 2013-12-25 2020-05-13 キヤノンメディカルシステムズ株式会社 医用画像処理装置、ｘ線診断装置及び医用画像処理プログラム
KR101664137B1 (ko) * 2014-12-30 2016-10-10 삼성전자주식회사 검출기 어셈블리, 이를 포함하는 컴퓨터 단층 촬영 장치 및 그 제어방법
US10098603B2 (en) * 2016-04-15 2018-10-16 Toshiba Medical Systems Corporation Method for estimation and correction of grid pattern due to scatter
JP6878590B2 (ja) * 2016-12-16 2021-05-26 プリズマティック、センサーズ、アクチボラグＰｒｉｓｍａｔｉｃ Ｓｅｎｓｏｒｓ Ａｂ スペクトルコンピュータ断層撮影データからの従来のコンピュータ断層撮影画像の再生
CN109685877B (zh) * 2018-12-27 2022-11-25 重庆大学 一种基于自适应投影图像特征区域匹配的微纳ct焦点漂移校正方法
CN109682843B (zh) * 2019-02-13 2021-07-06 重庆交通大学 一种对ct系统的参数标定方法
US11490874B2 (en) 2021-03-19 2022-11-08 Canon Medical Systems Corporation Self calibration method and apparatus for correcting offset angle in a photon counting computed tomography system
CN114041811A (zh) * 2021-11-03 2022-02-15 海南大学 模体和模体在扫描设备中的应用方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5822393A (en) * 1997-04-01 1998-10-13 Siemens Aktiengesellschaft Method for adaptively modulating the power level of an x-ray tube of a computer tomography (CT) system
US20040022364A1 (en) * 2002-05-31 2004-02-05 Karl Stierstorfer Method for determining correction coefficients for detector channels for a computed tomograph
US20060109950A1 (en) * 2004-11-24 2006-05-25 Arenson Jerome S Method and system of ct data correction
CN101061958A (zh) * 2006-04-28 2007-10-31 西门子公司 对ct系统进行散射校正的方法以及ct系统
CN101061957A (zh) * 2006-04-28 2007-10-31 西门子公司 在x射线ct中进行散射校正的方法和应用该方法的ct
CN101262819A (zh) * 2005-09-13 2008-09-10 皇家飞利浦电子股份有限公司 Ct散射的直接测量及校正
JP2009165518A (ja) * 2008-01-11 2009-07-30 Ge Medical Systems Global Technology Co Llc Ｘ線ｃｔ装置
US20100128844A1 (en) * 2008-11-26 2010-05-27 Brian Thomsen System and Method for Material Segmentation Utilizing Computed Tomography Scans
US20100195804A1 (en) * 2009-02-02 2010-08-05 Arineta Ltd. Method for calibration of a ct scanner

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL7710052A (nl) * 1977-09-14 1979-03-16 Philips Nv Inrichting voor computer-tomografie.
JPS5594241A (en) * 1979-01-11 1980-07-17 Hitachi Medical Corp Xxray transverse layer device
FR2632749B1 (fr) * 1988-06-10 1990-08-24 Gen Electric Cgr Procede et systeme d'etalonnage d'un scanner a rayons x en utilisant un ou plusieurs etalons circulaires excentres
US5473663A (en) * 1994-09-12 1995-12-05 General Electric Company Method for evaluating the performance of detectors in a computed tomography system
US5473656A (en) * 1994-09-15 1995-12-05 General Electric Company Computed tomography system with correction for z-axis detector non-uniformity
US5734691A (en) * 1996-12-23 1998-03-31 General Electric Company Detector z-axis gain non-uniformity correction in a computed tomography system
US5774519A (en) * 1997-01-30 1998-06-30 Analogic Corporation Method of and apparatus for calibration of CT scanners
US7136516B2 (en) * 2002-01-25 2006-11-14 General Electric Company Method and system for segmenting magnetic resonance images
US6848827B2 (en) * 2003-05-13 2005-02-01 General Electric Company Method and apparatus for calibrating detector spectral response
US7086780B2 (en) * 2004-05-20 2006-08-08 General Electric Company Methods for spectrally calibrating CT imaging apparatus detectors
DE102004062857A1 (de) 2004-12-27 2006-07-13 Siemens Ag Verfahren zur Ermittlung wenigstens eines Skalierungsfaktors für Messwerte eines Computertomographiegerätes
DE102006025401B4 (de) 2006-05-31 2015-12-24 Siemens Aktiengesellschaft Erfassungsvorrichtung und Verfahren zum Erfassen eines Objekts mittels einer Projektion durch das Objekt hindurch mittels Röntgenstrahlen und einem kontrast- und rauschoptimierten quantenzählenden Detektor
US7922390B2 (en) * 2007-03-29 2011-04-12 Varian Medical Systems, Inc. Method and apparatus to facilitate forming detector-level calibration information for a high energy-based scanner
DE102009055807B4 (de) 2009-11-26 2016-11-24 Siemens Healthcare Gmbh Schaltungsanordnung zur Zählung von Röntgenquanten einer Röntgenstrahlung mittels quantenzählender Detektoren sowie anwendungsspezifische integrierte Schaltung und Strahler-Detektor-System
US8315352B2 (en) * 2010-09-16 2012-11-20 General Electric Company System and method of spectral calibration and basis material decomposition for X-ray CT systems

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5822393A (en) * 1997-04-01 1998-10-13 Siemens Aktiengesellschaft Method for adaptively modulating the power level of an x-ray tube of a computer tomography (CT) system
US20040022364A1 (en) * 2002-05-31 2004-02-05 Karl Stierstorfer Method for determining correction coefficients for detector channels for a computed tomograph
US20060109950A1 (en) * 2004-11-24 2006-05-25 Arenson Jerome S Method and system of ct data correction
CN101262819A (zh) * 2005-09-13 2008-09-10 皇家飞利浦电子股份有限公司 Ct散射的直接测量及校正
CN101061958A (zh) * 2006-04-28 2007-10-31 西门子公司 对ct系统进行散射校正的方法以及ct系统
CN101061957A (zh) * 2006-04-28 2007-10-31 西门子公司 在x射线ct中进行散射校正的方法和应用该方法的ct
JP2009165518A (ja) * 2008-01-11 2009-07-30 Ge Medical Systems Global Technology Co Llc Ｘ線ｃｔ装置
US20100128844A1 (en) * 2008-11-26 2010-05-27 Brian Thomsen System and Method for Material Segmentation Utilizing Computed Tomography Scans
US20100195804A1 (en) * 2009-02-02 2010-08-05 Arineta Ltd. Method for calibration of a ct scanner

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104603606A (zh) * 2012-09-24 2015-05-06 西门子公司 用于确定通过待检查的对象引起的x射线辐射的衰减的方法和装置
US9664801B2 (en) 2012-09-24 2017-05-30 Siemens Aktiengesellschaft Method and device for determining the x-ray radiation attenuation caused by the object to be examined
CN104111469A (zh) * 2013-04-17 2014-10-22 上海联影医疗科技有限公司 模拟x射线穿透模体组件的衰减强度的方法
CN104111469B (zh) * 2013-04-17 2015-06-10 上海联影医疗科技有限公司 模拟x射线穿透模体组件的衰减强度的方法
CN103340643A (zh) * 2013-06-09 2013-10-09 成都军区昆明总医院 应用多功能标准体模对ct关键技术指标进行检测的方法
CN103340643B (zh) * 2013-06-09 2014-12-10 成都军区昆明总医院 应用多功能标准体模对ct关键技术指标进行检测的方法
CN107693037A (zh) * 2013-09-18 2018-02-16 上海联影医疗科技有限公司 Pet扫描装置及其时间偏移校正方法
CN107693037B (zh) * 2013-09-18 2021-04-30 上海联影医疗科技股份有限公司 Pet扫描装置及其时间偏移校正方法
CN109416405A (zh) * 2016-06-29 2019-03-01 皇家飞利浦有限公司 光子计数计算机断层摄影
CN113796879A (zh) * 2021-09-27 2021-12-17 北京万东医疗科技股份有限公司 一种球管出射能谱验证方法、装置、电子设备及存储介质
CN113796879B (zh) * 2021-09-27 2024-03-08 北京万东医疗科技股份有限公司 一种球管出射能谱验证方法、装置、电子设备及存储介质

Also Published As

Publication number Publication date
DE102010042388A1 (de) 2012-04-19
CN102440801B (zh) 2014-05-14
US8628241B2 (en) 2014-01-14
US20120093282A1 (en) 2012-04-19

Similar Documents

Publication Publication Date Title
CN102440801B (zh) 校准具有至少一个焦点-探测器组合的ct 系统的方法
CN101138501B (zh) 用于生成目标的多光谱图像的方法和系统
JP5290145B2 (ja) ダーティアイソトープのｐｅｔ再構成
Garayoa et al. A study on image quality provided by a kilovoltage cone‐beam computed tomography
CN102955164B (zh) 用于使多通道量子计数射线探测器的阈值均匀化的方法
CN103065340B (zh) 计算机断层摄影（ct）中用于扩张迭代近似重建的轴向范围的方法以及系统
EP1941457B1 (en) Imaging method and device with dual reading scanner
US9076255B2 (en) Method and system for reconstruction of tomographic images
CN101109719A (zh) 辐射成像装置和辐射成像方法
US20100119034A1 (en) Computed tomography device and method using circular-pixel position-adaptive interpolation
US20150146844A1 (en) Method and system for spectral computed tomography (ct) with inner ring geometry
JP2009529394A (ja) ラドンデータから画像関数を再構成する方法
CN103654833B (zh) Ct探测器偏转角的确定方法和装置
JP2008526284A (ja) エネルギー入力ビームをコリメートするための方法及び装置
WO2013192600A1 (en) Methods and systems for estimating scatter
CN1847883A (zh) 放射线诊断装置
CN102764137A (zh) 一种静态ct扫描仪及其散射x光子校正方法
CN103690185A (zh) 用于多层x射线计算机断层扫描系统的可调式光子探测系统
US20120163557A1 (en) Methods and phantoms for calibrating an imaging system
US9629602B2 (en) System and method for ultra-high resolution tomographic imaging
JP5265353B2 (ja) 断層撮像（ｉｍａｇｉｎｇｔｏｍｏｇｒａｐｈｙ）のための方法及び装置
CN109146800B (zh) 锥束计算机断层成像图像校正方法和系统
Thanasupsombat et al. A Simple Scatter Reduction Method in Cone‐Beam Computed Tomography for Dental and Maxillofacial Applications Based on Monte Carlo Simulation
AU2019201125B2 (en) Geometric alignment, sample motion correction, and intensity normalization of computed tomography projections using pi-line optimization
US11375964B2 (en) Acquisition method, acquisition device, and control program for tomographic image data by means of angular offset

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
TR01 Transfer of patent right

Effective date of registration: 20220126

Address after: Erlangen

Patentee after: Siemens Healthineers AG

Address before: Munich, Germany

Patentee before: SIEMENS AG

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240903

Address after: German Phu F Haim

Patentee after: Siemens Medical AG

Country or region after: Germany

Address before: Erlangen

Patentee before: Siemens Healthineers AG

Country or region before: Germany

TR01 Transfer of patent right