CN1705457A - 基于两种不同x射线频谱的成像方法 - Google Patents

基于两种不同x射线频谱的成像方法 Download PDF

Info

Publication number
CN1705457A
CN1705457A CN200480001240.8A CN200480001240A CN1705457A CN 1705457 A CN1705457 A CN 1705457A CN 200480001240 A CN200480001240 A CN 200480001240A CN 1705457 A CN1705457 A CN 1705457A
Authority
CN
China
Prior art keywords
distribution
atomic number
values
contrast agent
ray
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.)
Pending
Application number
CN200480001240.8A
Other languages
English (en)
Inventor
比约恩·海斯曼
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 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 CN1705457A publication Critical patent/CN1705457A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/481Diagnostic techniques involving the use of contrast agents
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/40Arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/405Source units specially adapted to modify characteristics of the beam during the data acquisition process
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/482Diagnostic techniques involving multiple energy imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/504Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/58Testing, adjusting or calibrating thereof
    • A61B6/582Calibration
    • A61B6/583Calibration using calibration phantoms

Landscapes

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

Abstract

本发明涉及对检查对象、特别是患者(P)进行成像检查的方法。其中,首先给检查对象注射造影剂(KM)。然后至少确定两个X射线衰减值的空间分布,这些空间分布分别表示局部X射线衰减系数(μ(x,y))或与其线性相关的量(C)。该两个空间分布至少包括:基于第一X射线频谱确定的第一衰减值分布(μ1(x,y));基于与第一X射线频谱不同的第二X射线频谱确定的第二衰减值分布(μ2(x,y))。通过对两个衰减值空间分布的分析,确定一个或多个预定的原子序数值(Z;Z1,Z2,···)的空间分布,或确定包含关于造影剂(KM)在检查对象体内分布信息的、未预定义的检查对象中现有的原子序数值的空间分布(Z(x,y))。将原子序数空间分布用于图像中的造影剂(KM)显示。

Description

基于两种不同X射线频谱的成像方法
技术领域
本发明涉及一种用于为了医疗的目的对尤其是患者的检查对象进行成像检查的方法。该方法尤其适合于应用在断层造影方法或者应用在可以进行成像断层造影的检查装置、例如X射线计算机断层造影设备中。
背景技术
造影剂例如由DE 4433564A1、WO 00/16811或者DE 10002939C1公开。
放射摄影方法的结果,如计算机断层造影、乳腺造影、血管造影、X射线检查技术等类似方法的结果,首先是对X射线从X射线源至X射线检测器的路径上衰减的显示(投影图像)。这种衰减是由于放射路径上被透射的介质或物质造成的,因此,这种衰减也被认为是对放射路径上所有像素的衰减系数的线性积分。特别是在断层造影方法中,例如在X射线计算机断层造影中,借助于所谓的再现方法可以从投影的衰减数据计算出各像素的衰减系数(μ),并由此得到比纯粹对投影图像进行分析敏感得多的检查。
为了显示衰减分布,代之以使用衰减系数通常采用对水的衰减系数标准化的值,即所谓的CT数。CT数根据通过测量确定的当前衰减系数μ和参考衰减系数μH2O按照下式计算得出:
C = 1000 × μ - μ H 2 o μ H 2 o [ HU ] - - - ( 1 )
其中,CT数C以Hounsfield[HU]为单位。对于水来说,其值CH2O=0HU,而对于空气,该值CL=-1000HU。
由于两种表示相互可以转换或等同,以下所选用的衰减值或衰减系数的概念既指衰减系数μ又指CT值。此外,在本发明说明书范围内物质和组织的概念可以交换地使用。假定,在医疗上显示的检查的上下文中物质可以是解剖组织,以及反之在材料检验和可靠性检验中的组织可以被理解为检查对象的任意物质。
尽管基于局部衰减系数(μ)的图像的说服力明显地被提高,不过在个别情况下其在解释图像时会产生问题。即,局部提高的衰减值要么回溯到较高原子序数的物质、例如骨骼中的钙或造影剂中的碘,要么回溯到较高的柔软部分的密度、例如在肺节点(Lungenknoten)中那样。在位置 处的局部衰减系数μ取决于入射到该组织或物质中的X射线能量E和该局部组织或物质的密度ρ,如下式所示:
μ = μ ( E , r ‾ ) = ( μ ρ ) ( E , Z ) × ρ ( r → ) - - - ( 2 )
其中, 为取决于能量和物质的质量衰减系数,Z为(有效)原子序数。
因此,由其有效原子序数Z确定的、物质的取决于能量的X射线吸收,与由物质密度ρ造成的X射线吸收重叠。因此,不同的化学及物理组成的物质在X射线图像中可能具有相同的衰减值。反之,不能从X射线照片的衰减值中得出检查对象的物质组成。
为了解决该问题,要求有用于表示物质特征值的方法。例如在US4247774中,与计算机支持的断层造影方法结合公开了,采用不同的X射线频谱或不同的X射线量子能量来产生图像。这种方法通常称为双频谱CT。其利用的是以原子序数为条件的衰减系数μ的能量依赖性,也就是说,其基于这样的效应:较高原子序数的物质或组织对较低能量的X射线的吸收要远大于原子序数较低的物质或组织。反之,在较高X射线能量的情况下,衰减值得到平衡并且主要是物质密度的函数。这样,在双频谱CT中例如对在两个不同的X射线管电压下拍摄的图像中的区别进行计算。
只要没有另外给出,则在本说明书范围中不是按照严格的、针对元素的意义上使用原子序数的概念,而是表示组织(有关物质)的有效原子序数,其可以从参与组织或物质构成的元素的化学原子序数和原子量中计算出。
如果在X射线拍摄时还附加地采用所谓的基本物质分解方法,则还可以得到更具体的结论,该方法如W.Kalender等在下列文献中描述的那样:“Materialselektive Bildgebung und Dichtemessung mit der Zwei-Spektren-Methode,I.Grundlagen und Methodik”,W.Kalender,W.Bautz,D.Flesenberg,C.Suess und E.Klotz,Digit.Bilddiagn.7,1987,66-77,Georg Thieme出版社。在该方法中,利用较低能量和较高能量的X射线来测量检查对象的X射线衰减值,并将所获得的值分别与两种基本物质的相应参考值进行比较,如(对于骨骼物质的)钙和(对于柔软部分组织的)水进行比较。在此假设每个测量值都表示这两种基本物质的测量值的线性叠加。因此,对于图像显示的检查对象的每种元素,可以通过与基本物质的值比较而计算出骨骼部分和软组织部分,从而可以将原始拍摄转换为对这两种基本物质的显示。
因此,基本物质分解以及双频谱方法适合于分开或者区分在人或者动物组织内具有巨大差别的原子序数的、预定的解剖结构或物质类型。
在德国专利申请文献DE 10143131中公开了一种方法,其敏感性和说服力超过基本物质分解方法,并且例如可以进行高说服力的功能CT成像。该方法通过分析X射线设备的影响频谱的测量数据计算出平均密度 和有效原子序数
Figure A20048000124000062
的空间分布。由此得到极好的对比度,特别是针对检查对象的化学及物理组成。例如,在组织中原子序数分布的显示允许观察待检查对象的生物化学组成、与迄今密度均匀表示的器官进行对比、定量地确定人体组成部分(如碘等)以及例如将钙化所基于的原子序数分离出来。
发明内容
本发明要解决的技术问题是,提供一种在依赖于物质或原子序数的、基于X射线的成像中可以改善灵敏度或提高说服力的方法。
本发明的上述技术问题是通过一种用于对检查对象、特别是患者进行成像检查的方法解决的,其中,
a)为该检查对象注射造影剂,
b)然后至少确定两个X射线衰减值的空间分布,这些X射线衰减值分别表示局部X射线衰减系数(μ(x,y))或与其线性相关的量(C),其中,所述两个空间分布至少包括:
-基于第一X射线频谱确定的第一衰减值分布,
-基于与该第一X射线频谱不同的第二X射线频谱确定的第二衰减值分布,
c)通过对所述两个衰减值空间分布的分析,确定一个或多个预定原子序数值(Z;Z1,Z2,...)的空间分布,或者确定包含关于所注射的造影剂在检查对象体内分布的信息的、未预定义的检查对象中现有的原子序数值的空间分布(Z(x,y)),和
d)将所述原子序数空间分布用于造影剂成像显示。
本发明基于这样的思路:造影剂的使用可以改善在X射线计算机断层造影中的功能性成像。在此,造影剂迄今为止仅用来例如将血液在吸收中从其组织背景中消除。没有进行对物质或组织选择地分析。此外,本发明还基于这样的认知:通过给定可承受剂量或浓度的造影剂,可以实现借助于两种不同的X射线频谱测量原子序数的差。
在按照本发明的方法中,可以预定造影剂的原子序数值。该方法尤其可以与本文开始部分提到基本物质分解进行组合。
优选地,将原子序数空间分布作为二维或者三维场确定,其中每个场值是在由有关场表示的位置上的局部原子序数值。该方法尤其可以与本文开始部分提到德国专利申请文献DE 10143131的方法进行组合。该专利申请的公开内容、特别是其中的权利要求1至7明确地被引入到本专利申请中来。
此外,优选除了原子序数值分布之外,还确定另一个其场值分别反映局部密度值的二维或者三维场。
例如,可以将原子序数的空间分布这样用于成像:显示一幅图像,该图像仅显示了来自特定的(例如包括造影剂的原子序数的值的)原子序数区间或者超过一个特定原子序数值的数据。也可以将所测量的原子序数值转换为灰度值刻度或者彩色刻度,其中,可以将造影剂的原子序数值突出出来或者单独地进行染色,并且将该刻度成像地示出。可以为这种图像衬底或者叠加一幅常规的、非功能性的衰减图像。
按照一种特别优选的实施方式,将所确定的具有原子序数值的场和所确定的具有密度值的场用于计算造影剂的局部浓度或者局部数量。
在本发明的范围中,造影剂被理解为各种在给入检查对象中、特别是注射到患者体内之后导致吸收(即在X射线图像)中的对比度改善或者对比度增强的物质。其中包括常规的、如例如在灌注测量中注射入血管中使其在图像中突出的造影剂。不过,“造影剂”还可被理解为这样的物质,其例如按照钥匙-锁原理(Scluessel-Schloss-Prinzip)专门或有选择性地沉积或聚集在检查对象中的特定位置并因此允许对器官功能进行检查。后面提到的这种材料也可以是所谓的标记物或示踪物。这种标记物例如由具有对待检查的对象结构的高亲和力的生物大分子(如抗体、缩氨酸或糖分子)以及由在X射线图像中很易看出的(例如掺杂的)对比物质构成。大分子例如用作所谓的“新陈代谢标志”,其作用是,使得也被整体称为新陈代谢标记物的造影剂仅仅聚集在特定的区域、如肿瘤、炎症或其它病灶处。造影剂例如由本文开始部分提到的文献公开。
优选地,使用原子序数值大于20或者大于40的造影剂。造影剂尤其是原子序数值小于83或小于70的造影剂。
特别优选地,造影剂包括钆(Gadolinium)、碘(Jod)、镱(Ytterbium)、镝(Dysposium)、铁(Eisen)和/或铋(Wismut)。
按照另一种优选的实施方式,造影剂包含一种有机化合物,特别是例如糖的脂肪族碳水化合物和/或氨基酸或者缩氨酸。
造影剂可以构成为选择性地沉积在检查对象的特定位置上或特定的组织成分中。
在优选的实施方式中,在重量浓度为10-4至10-7的范围内、特别是在10-5至10-6的范围内添加所述造影剂(KM)。
在本说明书中所使用的概念“X射线频谱”具有比仅由设备的X射线源发射的X射线的频谱分布更多一层含义。在X射线检测器一侧也可以利用不同的效率转换不同的频谱成分并因此不同地加权。由此产生出的有效频谱分布在本说明书中也称为X射线频谱。
两个衰减值分布不必一定依次作为两幅图像利用不同的管电压来记录。因为每个X射线管发射出具有一定宽度的频谱,所以在对于选择频谱的实施方式中所属的接收单元也可以进一步或者完全同时地接收两个衰减值分布。为此,例如可以在射线路径上使用可以调整的滤波器和/或两个分开设置的X射线检测器阵列。
特别是利用选择量子能量的X射线检测器阵列来实现用于实施本发明方法的接收单元。
特别是,就使用在本文开始部分提到的德国专利申请文献DE 10143131中描述的方法而言,特别有利的是,确定密度和原子序数的第一衰减值分布的第一衰减值的第一函数关系,以及密度和原子序数的第二衰减值分布的、与第一衰减值对应的第二衰减值的至少第二函数关系,并且从第一函数关系与第二函数关系以及必要时其它函数关系的比较中,确定原子序数的空间分布以及可选择的密度的空间分布。
在此,优选地借助于对校准试样(Eichprobe)的参考测量或者以在物理模型基础上进行模拟的形式,至少对一个X射线频谱确定密度和原子序数的衰减值的函数关系。
按照另一种优选的实施方式,在确定密度和原子序数值对的基础上对第一衰减值分布和其它衰减值分布的每个对应衰减值这样将衰减值分布变换为密度分布和原子序数分布,使得该值对满足对于第一X射线频谱和至少另一个X射线频谱的、密度和原子序数的X射线吸收的特定函数关系。由此,对于图像元素的密度和原子序数可以简单地作为所记录的X射线吸收值分布的相互对应的X射线吸收值的函数关系的交集(Schnittmenge)来计算。
优选地,第一X射线频谱具有这样的量子能量,即,其相对于第二X射线频谱的量子能量通过光子效应有利于X射线吸收,使得在确定原子序数时得到高分辨率。
在本发明的一种优选的实施方式中,为了改变X射线频谱以便记录检查对象,至少改变X射线管的一个运行参数,其中,X射线源在第一运行状态发射第一X射线频谱,而在第二运行状态发射与第一X射线频谱不同的第二X射线频谱,从而可以在两个X射线频谱之间进行较快的转换。
此外,为了改变X射线频谱以便记录检查对象,有利的是改变检测器的特性,其中,X射线检测器将从X射线源接收的X射线的频谱的部分区域变换成相互独立的电信号,并且在此允许在不同的X射线频谱下对X射线吸收的分布同时进行记录。
附图说明
下面对照附图、根据实施方式对本发明作进一步的说明。图中:
图1示出按照本发明方法的一种实施方式的流程图,
图2示出根据等吸收曲线在不同组成的物质下相同衰减值μ的形成,
图3a举例示出作为根据图1方法的一部分的、用于确定等吸收曲线的计算方法的函数方案,
图3b举例示出作为根据图1方法的一部分的、将X射线衰减值变换为物质密度和原子序数值的流程图,以及
图4示出在两种不同的X射线频谱下一种组织类型的两条等吸收曲线。
具体实施方式
图1中作为流程图示意地示出了按照本发明方法的一种实施方式。在第一步骤1中,对患者P注入标记物或者造影剂KM,例如通过注射到血管中或者通过服用。然后,在一个仅示意示出的X射线计算机断层造影设备2中对患者P进行检查,并且是既对第一X射线频谱S1进行分析、也(同时或者随后)对第二X射线频谱S2进行分析(第二步骤3),两种X射线频谱通过对X射线计算机断层造影设备2的对应设置而加以选择。借助于在这样获得的原始数据的基础上进行图像再现(第三步骤5),对每个X射线频谱S1、S2产生一个衰减值分布,例如作为在具有坐标x和y的横向层析图像内的(线性)X射线衰减系数μ的分布μ1(x,y)以及μ2(x,y)。在第四步骤7中进行计算机支持的将X射线衰减系数的分布μ1(x,y)以及μ2(x,y)变换为原子序数分布Z(x,y)。原子序数分布Z(x,y)在第五步骤8中被用于在显示器9上显示造影剂KM的分布。
如果假定例如注射具有在106水分子上的一个Gd原子(对应于大约9ppm重量成分)的基于Gd的标记物,则产生效果Z=7.52。与水值Z=7.42相比,利用根据图3a和3b的方法可以证实该浓度。根据图3a和3b的方法可以作为第四步骤7的一部分来执行。该方法的细节在申请号为10143131的德国专利申请文献中有所描述,该文献明确地被引用。
如果要同时确定密度分布ρ(x,y),则可以在已知的(例如类似于水的)载体物质条件下定量地确定标记物质的浓度。
造影剂KM的原子序数应该具有与本底物质的原子序数、典型的为Z=7.42的水尽可能高的偏差。
作为对于根据图3a和3b的方法的替换,也可以使用本文开始部分描述的基本物质分解,以便确定在横向层面(x,y)中事先固定的两个或多个原子序数Z1,Z2,...的空间分布。这种分布同样可以用于在图像中显示造影剂KM。
为了解释根据图3a和3b的方法,下面首先解释:图2的等吸收曲线11连接了所有这样的值对(ρ,Z),它们在定义的X射线频谱下具有相同的衰减值μ以及C。图1的显示表明,不能单独从关于组织或物质的类型和组成的信息支持地推导出X射线图像衰减值。
X射线由不同的物质并且依赖于X射线的能量而不同强度地被衰减。这点可以归结于在不同物质下不同作用的衰减机制。
在本说明书的范围中,例如从参与构造的元素的原子序数Zi、其原子量Ai和其局部物质等效密度ρi中简化地计算特定组织类型的实际有效的原子序数Z作为原子序数:
Z = { Σ i ρ i A i ρ i Z i 4 Σ i ρ i A i ρ i } 1 3 - - - ( 3 )
对于纯钙得到ZCa=20,对于氢化钙大约为ZCaH2≌16.04,而对于水大约为ZH2O≌7.428。因此,可以通过原子序数Z很好地采集对象的化学或生物化学组成。
计算检查区域中的原子序数和密度分布的前提是,至少有两幅在拍摄几何上相同、但利用所使用的X射线的不同能量产生的该区域的X射线照片。在使用多于两幅利用不同的X射线能量记录的X射线照片的条件下,可以改善Z和ρ分辨率,但是由此也提高了射线负担。因此,在对患者进行检查的情况下并不总存在这种可能性。
将基于图像数据的衰减值变换为原子序数和物质或组织密度的分布图像的出发点是,了解对于X射线设备的每个X射线频谱的等吸收曲线。如已经提到的那样,在此X射线频谱不是被理解为由设备的X射线源发射的X射线的频谱分布的窄概念,而是一个扩展的概念,其还考虑在X射线检测器一侧对X射线管的发射频谱的不同频谱范围的不同加权。因此,所测量的衰减值由对X射线管发射的射线频谱的直接衰减以及所使用的X射线检测器的频谱效率给出。两个值都是针对设备的量,并且必须直接或间接地借助校准试样的衰减值来确定。它们是计算等吸收曲线的基础。
在图3a中示意性示出了用于建模或者用于计算一簇等吸收曲线的三种方法300,即,理论建模、试验确定以及利用通过试验确定的参数的曲线校准的理论建模。
原理上,要建立与在X射线照片中要求的用于覆盖X衰减的间隙的衰减值一样多的等吸收曲线。在此,不是对于每个在理论上出现的衰减值都计算一个等吸收曲线;而是可以根据需要通过插值或者其它适当的平均方法得到的未计算的等吸收曲线。
在图3a的流程图的左侧分支示出了理论建模的基本步骤。在步骤S302中,首先作为参数读入对于具有可用X射线管电压的设备特定的X射线发射频谱S(E)的数据。为此,X射线的频谱分布可以在准备阶段对每个X射线设备试验地测定,或者可以采用对专门X射线源类型特有的数据。在步骤S303确定检测器设备函数w(E)。为此,也可以在准备阶段对检测器装置进行精确地测量,或者采用对于检测器类型特有的数据、例如其频谱的技术指标。在步骤S304中,在物理模型的基础上以曲线簇Ci(ρ,Z)以及μi(ρ,Z)的形式计算等吸收曲线,该物理模型对于S(E)和w(E)的每个重要组合模拟了对于具有不同原子序数和在不同的物质密度下的物质的X射线衰减Ci以及μi
作为对于步骤S302至步骤S304的理论建模的替代,也可以试验地确定等吸收曲线的曲线簇。为此,在步骤S305在X射线设备中在S(E)和w(E)的不同重要组合的条件下测量具有不同的密度和平均原子序数的校准物质的X射线衰减。测量值构成了对于随后在步骤S306中对等吸收曲线的曲线簇Ci以及μi的支撑点。
作为另一种替代,可以利用试验确定的X射线衰减值对在理论基础上模型化的曲线族Ci以及μi进行校准。在步骤S307中,测量用于校准理论曲线簇所需的衰减值,如在上面步骤S305中描述的那样利用适当的校准物质或者模型在X射线设备中进行测量。与步骤S302至步骤S304的纯理论建模不同,在该方法中不以精确地了解X射线发射频谱S(E)和w(E)为前提条件,而是以在步骤S308中对等吸收曲线的曲线簇Ci以及μi建模的参数为前提。在步骤309中利用在步骤S307中试验地确定的校准值对曲线校准,最后定义了这些参数的、对于X射线设备的X射线发射频谱和检测器设备函数特定的值。
通过确定所需的X射线衰减值以及S(E)和w(E)的组合的等吸收曲线,实现了用于变换图像数据、在图像数据中表示X射线通过组织时的衰减值进行变换的前提条件其中,该衰减值表示在相应组织中原子序数以及物质密度的分布。
根据不同的任务表述也可以将三种方法混合地用于等吸收曲线的确定。例如,可以将由试验不能精确确定或者要在利用很大花费下才能确定或者根本不能确定的值,借助理论建模来完善或者以其精确性来精确表示。然后,在步骤S310中将利用不同方法得到的数据综合为一个统一的数据组,并且在步骤S311为图像变换做好准备。
图3b中示出了一种适合于本发明方法的变换方法320。该方法依赖于按照上述方法300之一确定的并且作为数据组在步骤S321准备好的等吸收曲线的曲线簇。
变换逐图像元素地进行。下面,假设对于基于两个在不同的X射线频谱、但在相同的拍摄几何下拍摄的X射线图像的X射线衰减值进行变换。这是实施按照本发明的变换的最小前提条件。但也可以采用在多于两个X射线的不同能量分布条件下的多于两个X射线照片。
在步骤S322中对待变换的图像元素进行选择,并且在随后的步骤S323中,从第一X射线照片中对该图像元素读入衰减值曲线簇C1以及μ1,和从第二X射线照片中读入C2以及μ2。在随后的步骤S324中,调用对于第一X射线照片使用的X射线频谱S1(E)和检测器设备函数w1(E),以及对于第二X射线照片对应的值S2(E)和w2(E)。这些值构成了用于随后选择与各个衰减值对应的等吸收曲线的参数。在此,频谱分布Si(E)和wi(E)也可以间接地确定,例如通过调用所使用的管电压U1以及U2或者X射线检测器的运行参数来确定。
在步骤S325中,从在步骤S321中提供的等吸收曲线数据组中选择出,在参数S1(E)和w1(E)下满足条件C1以及μ1的第一曲线和在参数S2(E)和w2(E)下满足条件C2以及μ2的第二曲线。在图4中示出了这样得到的第一等吸收曲线11和第二等吸收曲线41的例子。
在步骤S326中作为两条曲线11和41的交集计算出交点42。该曲线交点42例如可以通过局部的线性变换或者借助于迭代的交点求解得到。因为两条曲线11和41代表着同一图像元素的、因此对于待检查组织的同一部分区域的两种不同的衰减值,所以两种衰减值必须由同一物质或者组织类型引起。因此,曲线交点的坐标(ρ,Z)显现出了该图像元素所属的组织部分区域的物质密度和原子序数。
最后,在步骤S327中将这样确定的原子序数Z作为对应的图像元素值写入原子序数分布中,并在步骤S328中将类似确定的物质密度值ρ写入到密度分布中。将步骤S322至S328对于所有其余图像点进行重复,直到随后可以在步骤S329中输出图像。其中,可以跳过步骤S324,因为对于一幅图像的所有图像元素频谱分布Si(E)和wi(E)是相同的。

Claims (14)

1.一种用于对检查对象、特别是患者(P)进行成像检查的方法,其中,
a)为该检查对象注射造影剂(KM),
b)然后至少确定两个X射线衰减值的空间分布(μ1(x,y),μ2(x,y)),这些X射线衰减值分别表示局部X射线衰减系数(μ(x,y))或与其线性相关的量(C),其中,所述两个空间分布(μ1(x,y),μ2(x,y))至少包括:
-基于第一X射线频谱确定的第一衰减值分布(μ1(x,y)),
-基于与该第一X射线频谱不同的第二X射线频谱确定的第二衰减值分布(μ2(x,y)),
c)通过对所述两个衰减值空间分布(μ1(x,y),μ2(x,y))的分析,确定一个或多个预定原子序数值(Z;Z1,Z2,...)的空间分布,或者确定包含关于所注射的造影剂(KM)在检查对象体内分布的信息的、未预定义的检查对象中现有的原子序数值的空间分布(Z(x,y)),和
d)将所述原子序数空间分布(Z(x,y))用于造影剂(KM)成像显示。
2.根据权利要求1所述的方法,其中,预先定义所述造影剂(KM)的原子序数值。
3.根据权利要求1所述的方法,其中,将所述原子序数空间分布作为二维或者三维场确定,其中每个场值是在由有关场表示的位置(x,y)上的局部原子序数值(Z(x,y))。
4.根据权利要求3所述的方法,其中,除了所述原子序数分布之外,还确定另一个其场值分别反映局部密度值(ρ(x,y))的二维或者三维场。
5.根据权利要求4所述的方法,其中,将所确定的具有原子序数值(Z(x,y))的场和所确定的具有密度值(ρ(x,y))的场,用于计算造影剂局部浓度或者局部数量。
6.根据权利要求1至5中任一项所述的方法,其中,使用原子序数值大于20的造影剂(KM)。
7.根据权利要求6所述的方法,其中,使用原子序数值大于40的造影剂(KM)。
8.根据权利要求1至7中任一项所述的方法,其中,使用原子序数值小于83、特别是小于70的造影剂(KM)。
9.根据权利要求1至8中任一项所述的方法,其中,所述造影剂(KM)包含钆、碘、镱、镝、铁和/或铋。
10.根据权利要求1至9中任一项所述的方法,其中,所述造影剂(KM)包含一种有机化合物,特别是例如糖的脂肪族碳水化合物。
11.根据权利要求1至10中任一项所述的方法,其中,所述造影剂(KM)包含氨基酸或者缩氨酸。
12.根据权利要求1至11中任一项所述的方法,其中,将所述造影剂(KM)选择性地沉积在检查对象的特定位置上或特定的组织成分中。
13.根据权利要求1至12中任一项所述的方法,其中,在重量浓度为10-4至10-7的范围内、特别是在10-5至10-6的范围内添加所述造影剂(KM)。
14.根据权利要求1至13中任一项所述的方法,其中,确定密度和原子序数的第一衰减值分布的第一衰减值的第一函数关系(11),以及至少确定密度和原子序数的第二衰减值分布的、与第一衰减值对应的第二衰减值的第二函数关系(41),并且其中,从该第一函数关系(11)与第二函数关系(41)以及必要时其它函数关系的比较中,确定原子序数空间分布,并且可选地确定空间密度分布。
CN200480001240.8A 2003-03-14 2004-03-02 基于两种不同x射线频谱的成像方法 Pending CN1705457A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE10311628.1 2003-03-14
DE10311628A DE10311628B4 (de) 2003-03-14 2003-03-14 Bildgebungsverfahren

Publications (1)

Publication Number Publication Date
CN1705457A true CN1705457A (zh) 2005-12-07

Family

ID=32945926

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200480001240.8A Pending CN1705457A (zh) 2003-03-14 2004-03-02 基于两种不同x射线频谱的成像方法

Country Status (5)

Country Link
US (1) US7319739B2 (zh)
JP (1) JP4358852B2 (zh)
CN (1) CN1705457A (zh)
DE (1) DE10311628B4 (zh)
WO (1) WO2004080308A1 (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101015457B (zh) * 2005-12-23 2011-11-09 西门子公司 X射线装置的操作方法
CN102665557A (zh) * 2009-12-16 2012-09-12 佳能株式会社 X射线图像处理设备、x射线图像处理方法和计算机程序的存储介质
CN103654826A (zh) * 2012-08-31 2014-03-26 西门子公司 用于检测硅树脂植入物损伤的方法和计算机断层成像设备
CN109459452A (zh) * 2018-12-26 2019-03-12 上海联影医疗科技有限公司 有效原子序数计算方法、装置和存储介质
CN113281359A (zh) * 2021-06-22 2021-08-20 国家卫生健康委职业安全卫生研究中心(国家卫生健康委煤炭工业职业医学研究中心) 一种基于ct技术进行射线安检物性识别的方法和装置

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2003252103A1 (en) * 2002-07-23 2004-02-09 Ge Medical Systems Global Technology Company, Llc Methods and systems for detecting components of plaque
DE10347971B3 (de) 2003-10-15 2005-06-09 Siemens Ag Verfahren und Vorrichtung zur Bestimmung des Flüssigkeitstyps einer Flüssigkeitsansammlung in einem Objekt
DE102004008519B4 (de) 2004-02-20 2010-11-25 Siemens Ag Verfahren zur Visualisierung quantitativer Information in Datensätzen der medizinischen Bildgebung
US7209536B2 (en) * 2004-11-19 2007-04-24 General Electric Company CT colonography system
US7599465B2 (en) * 2004-11-19 2009-10-06 General Electric Company Detection of thrombi in CT using energy discrimination
DE102004060580A1 (de) * 2004-12-16 2006-06-29 Siemens Ag Verfahren zur Erzeugung einer computertomographischen Darstellung von Gewebestrukturen mit Hilfe einer Kontrastmittelapplikation
DE102005023165A1 (de) * 2005-05-19 2006-11-23 Siemens Ag Medizinisches Bildgebungssystem mit einem um einen Patienten verfahrbaren Teil und Kollisionsschutzverfahren
US20070025514A1 (en) * 2005-06-06 2007-02-01 Ruediger Lawaczeck X-ray arrangement for graphic display of an object under examination and use of the x-ray arrangement
DE102006009222B4 (de) * 2006-02-28 2008-02-28 Siemens Ag Verfahren und Vorrichtung zur Bestimmung der Konzentration einer Substanz in einem Körpermaterial mittels Mehr-Energie-Computertomographie
JP5105589B2 (ja) * 2007-07-11 2012-12-26 株式会社日立メディコ X線ct装置
US8290232B2 (en) * 2008-02-15 2012-10-16 Mayo Foundation For Medical Education And Research System and method for quantitative imaging of chemical composition to decompose more than two materials
US8218837B2 (en) * 2008-06-06 2012-07-10 General Electric Company Material composition detection from effective atomic number computation
CN101647706B (zh) * 2008-08-13 2012-05-30 清华大学 高能双能ct系统的图象重建方法
US9008372B2 (en) 2011-05-31 2015-04-14 Schlumberger Technology Corporation Method for determination of spatial distribution and concentration of contrast components in a porous and/or heterogeneous sample
US9513233B2 (en) * 2011-10-28 2016-12-06 The University Of Chicago Color x-ray histology for multi-stained biologic sample
JP2015520846A (ja) * 2012-05-01 2015-07-23 アナロジック コーポレイション Ct密度画像とスパースマルチエネルギーデータを用いてボクセルセットのz−有効値を決定する方法およびシステム{determinationofz−effectivevalueforsetofvoxelsusingctdensityimageandsparsemulti−energydata}
DE102012217301B4 (de) * 2012-09-25 2021-10-14 Bayer Pharma Aktiengesellschaft Kombination aus Kontrastmittel und Mammographie-CT-System mit vorgegebenem Energiebereich und Verfahren zur Erzeugung tomographischer Mammographie-CT-Aufnahmen durch diese Kombination

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4149081A (en) * 1976-11-29 1979-04-10 Varian Associates, Inc. Removal of spectral artifacts and utilization of spectral effects in computerized tomography
US4247774A (en) * 1978-06-26 1981-01-27 The United States Of America As Represented By The Department Of Health, Education And Welfare Simultaneous dual-energy computer assisted tomography
NL8003354A (nl) * 1980-06-09 1982-01-04 Philips Nv Stralingsonderzoekapparaat met beeldsubtractie.
US4709382A (en) * 1984-11-21 1987-11-24 Picker International, Inc. Imaging with focused curved radiation detectors
US4662379A (en) * 1984-12-20 1987-05-05 Stanford University Coronary artery imaging system using gated tomosynthesis
US5132998A (en) 1989-03-03 1992-07-21 Matsushita Electric Industrial Co., Ltd. Radiographic image processing method and photographic imaging apparatus therefor
DE4433564A1 (de) 1994-09-07 1996-04-04 Schering Ag Peptidische Kontrastmittel für die Röntgen-, Magnetresonanz- und Nukleardiagnostik, deren Herstellung und Verwendung
AU1366697A (en) * 1995-12-26 1997-07-28 Holomed Aps A method and system for generating an x-ray image
JP2000086538A (ja) * 1998-09-17 2000-03-28 Toshihiro Akaike Mri用造影剤
DE10002939C1 (de) * 2000-01-13 2001-09-20 Schering Ag Paramagnetische DOTA-Derivate, Verfahren zu ihrer Herstellung und ihre Verwendung für die Herstellung von pharmazeutischen Mitteln
DE10127267B4 (de) 2001-06-05 2008-01-03 Siemens Ag Bildgebendes medizinsches Röntgengerät
DE10143131B4 (de) 2001-09-03 2006-03-09 Siemens Ag Verfahren zur Ermittlung von Dichte- und Ordnungszahlverteilungen bei radiographischen Untersuchungsverfahren

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101015457B (zh) * 2005-12-23 2011-11-09 西门子公司 X射线装置的操作方法
CN102665557A (zh) * 2009-12-16 2012-09-12 佳能株式会社 X射线图像处理设备、x射线图像处理方法和计算机程序的存储介质
CN103654826A (zh) * 2012-08-31 2014-03-26 西门子公司 用于检测硅树脂植入物损伤的方法和计算机断层成像设备
US9211066B2 (en) 2012-08-31 2015-12-15 Siemens Aktiengesellschaft Method for detecting damage to silicone implants and computed tomography device
CN109459452A (zh) * 2018-12-26 2019-03-12 上海联影医疗科技有限公司 有效原子序数计算方法、装置和存储介质
CN109459452B (zh) * 2018-12-26 2021-10-22 上海联影医疗科技股份有限公司 有效原子序数计算方法、装置和存储介质
CN113281359A (zh) * 2021-06-22 2021-08-20 国家卫生健康委职业安全卫生研究中心(国家卫生健康委煤炭工业职业医学研究中心) 一种基于ct技术进行射线安检物性识别的方法和装置

Also Published As

Publication number Publication date
DE10311628A1 (de) 2004-10-07
WO2004080308A1 (de) 2004-09-23
JP4358852B2 (ja) 2009-11-04
JP2006520223A (ja) 2006-09-07
US20060269043A1 (en) 2006-11-30
DE10311628B4 (de) 2006-04-13
US7319739B2 (en) 2008-01-15

Similar Documents

Publication Publication Date Title
CN1705457A (zh) 基于两种不同x射线频谱的成像方法
US7457450B2 (en) Method for recording and evaluating image data with the aid of a tomography machine
JP4309760B2 (ja) X線撮影検査法における密度および原子番号の分布を求める方法
Clark et al. Micro-CT of rodents: state-of-the-art and future perspectives
US7856134B2 (en) Method for generating image data relating to a virtually prescribable X-ray tube voltage from first and second CT image data
US7778454B2 (en) Method and apparatus for determining the concentration of a substance in a body material by means of multi-energy computed tomography
US9211104B2 (en) Selection method for two contrast media for use in a dual-energy CT examination, contrast media combination and generation of CT images using a contrast media combination and different energy spectra
US8290232B2 (en) System and method for quantitative imaging of chemical composition to decompose more than two materials
US8364240B2 (en) Method and device for automatically differentiating types of kidney stones by means of computed tomography
EP3271898B1 (en) Image reconstruction based on energy-resolved image data from a photon-counting multi bin detector
JP6139207B2 (ja) X線撮像における造影剤推定のシステム及び方法
US20060067473A1 (en) Method and system for multi-energy tomosynthesis
CN101495039A (zh) 用于能谱ct的定量物质分解
Zaidi et al. Advances in attenuation correction techniques in PET
CN1614400A (zh) 位置分辨地确定检查对象内元素浓度的方法和装置
Sauter et al. Perfusion-ventilation CT via three-material differentiation in dual-layer CT: a feasibility study
US7197171B2 (en) Nuclear imaging
Rajendran MARS spectral CT technology for orthopaedic applications
EP3484368B1 (en) Identifying types of adipose tissue
Walter et al. Dual kVp material decomposition using flat-panel detectors
Matanaghi Bone-related applications of spectral photon counting CT
Valton et al. Evaluation of tomographic reconstruction methods for small animal microCT and microPET/CT
Bosmans et al. Physics of Medical Imaging
RU2254056C2 (ru) Способ получения рентгеновского изображения (варианты)
Xu et al. Chest tomosynthesis image enhancement by bone suppression using convolutional neural networks with synthetic data

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication