CN104427938A - 用于产生相位对比显示的方法和x射线*** - Google Patents
用于产生相位对比显示的方法和x射线*** Download PDFInfo
- Publication number
- CN104427938A CN104427938A CN201380036091.8A CN201380036091A CN104427938A CN 104427938 A CN104427938 A CN 104427938A CN 201380036091 A CN201380036091 A CN 201380036091A CN 104427938 A CN104427938 A CN 104427938A
- Authority
- CN
- China
- Prior art keywords
- phase
- electron density
- check object
- ray
- contrast
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 230000010363 phase shift Effects 0.000 claims abstract description 34
- 238000003384 imaging method Methods 0.000 claims description 21
- 239000000463 material Substances 0.000 claims description 18
- 230000000694 effects Effects 0.000 claims description 10
- 238000000354 decomposition reaction Methods 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 abstract description 13
- 230000001419 dependent effect Effects 0.000 abstract description 2
- 238000010521 absorption reaction Methods 0.000 description 20
- 238000005259 measurement Methods 0.000 description 11
- 230000009977 dual effect Effects 0.000 description 6
- 230000000875 corresponding effect Effects 0.000 description 4
- 230000005855 radiation Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 238000002083 X-ray spectrum Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005305 interferometry Methods 0.000 description 2
- 238000010587 phase diagram Methods 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000011630 iodine Substances 0.000 description 1
- 229910052740 iodine Inorganic materials 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 238000011005 laboratory method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000011573 trace mineral Substances 0.000 description 1
- 235000013619 trace mineral Nutrition 0.000 description 1
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/48—Diagnostic techniques
- A61B6/484—Diagnostic techniques involving phase contrast X-ray imaging
-
- 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/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
-
- 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/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- 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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/419—Imaging computed tomograph
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Public Health (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Optics & Photonics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Physiology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及一种用于产生检查对象(P)的相位对比显示的方法和X射线***,其中,通过利用至少两个不同的X射线能量谱(2,3;4,5)确定对于X射线辐射的能量相关的衰减值来确定在检查对象(P)中的电子密度分布,从先前确定的电子密度分布中确定相移值,并且从所计算的相移值中产生相位对比显示。
Description
技术领域
本发明涉及一种用于产生检查对象的相位对比显示的方法和一种用于执行该方法的X射线***。
背景技术
通过所谓的复数折射率表明了材料关于X射线光学性质的特征。虽然常规的具有固定光谱的X射线成像直接测量复数折射率的虚部,但是其不能访问对X射线辐射的相移进行描述的实部。存在以下观点,即可以将相位信息用于医学诊断,以便更好地分离软组织。
在过去研发了多种方法,其能够显示检查对象对穿透检查对象的电磁波(特别是特定能量的X射线)的相位的影响。通常将这样的显示称为相位对比拍摄或者断层成像的相位对比显示。关于这种已知技术的概况例如在Raupach R.,Flohr T.的出版物“Analytical evaluation of the signal and noisepropagation in X-ray differential phase-contrast computed tomography”;Phys.Med.Biol.2011,56:2219-2244以及其给出的其它参考目录中给出。在该方法中,连续尝试直接测量并图像显示在辐射穿过时出现的相移。
然而至今显示,至今所建议的方法虽然部分能够在实验室技术并使用高剂量下实现并且提供好的图像数据,但是在对于有生命的对象视作可承受的剂量负担范围中的实现似乎导致不充足的且非常高噪声的图像结果。
发明内容
因此,本发明要解决的技术问题是,找到一种用于基于透射的电磁辐射的相移值,图像地再现检查对象的方法,该方法在利用对于有生命的检查对象视作可承受的剂量负担的检查范围中提供尽可能低噪声的图像结果。
该技术问题通过独立权利要求的特征解决。本发明的优选扩展是从属权利要求的内容。
发明人认识到以下内容:
通过相位敏感的X射线方法(相位对比成像=PCI)可以测量相位信息。此外,已知除了评估信号衰减也评估X射线辐射的相位的多种可能性。然而,所有方法的共同点是,测量首先提供相位信息的空间导数,也就是微分信号。当然,通过积分能够从中重建绝对相位,但是具有以不好的方式损害噪声功率谱的结果:低频的噪声部分被放大。这特别是使CT重建下的投影的强度值的或吸收系数的量化意义和稳定性变差。在相同的信噪比(SNR)下,要考虑结构的更差的探测率。只有在非常高的空间分辨率下(以及与此关联的高的剂量值下),关于在相同剂量下可实现的SNR,微分测量相对于吸收才有优势。
然而,例如在CT***中,在不相应提高剂量以获得对于诊断最低需要的SNR情况下,不可能简单地提高空间分辨率。因此,只有当拥有具有显著改善的空间一致性的紧凑的X射线源时,相位信息的潜在增值才能按照与剂量无关的方式用于计算机断层成像。
另外,所需要的PCI设备相比常规的成像设备在技术上更复杂,机械上具有巨大挑战,并且因此明显更贵。在许多PCI设备下,直接测量相位显著延长测量时间,这主要原因是由于用于控制辐射的一致性的措施(例如焦点处的光栅(源光栅))而降低的X射线流,以及由于在干涉测量方法(例如“相位步进扫描”)下用于实际观察干涉的技术。
另外原则上已知,通过两个主要的物理过程,即通过光电效应和康普顿效应(μ=μPhoto+μCompton)确定(低于511keV的能量范围中的)X射线辐射的吸收,其中,康普顿效应基本上直接与观察材料的电子密度成比例(μCompton~Z/E~ρe),并且光电效应具有强的能量依赖性(μPhoto=Z3.8/E3~ρe(Z/E)3)。从利用分别不同的能量谱或分别不同的能量的至少两个吸收测量中可以确定各个效应在衰减中的分量,从而通过康普顿效应的部分能够确定被透射的材料的电子密度。
替代地,借助利用不同能量的至少两个衰减测量也进行将被检查的对象分成两个或更多的主要基本材料的材料分解。如果识别出从中产生的基本材料分量,则也可以从这样的衰减测量中确定在检查对象中存在的电子密度,因为对于各自材料的电子密度已知。
当然,也已知由材料的电子密度确定,关于透射时的相位改变,材料对于透射的电磁波的影响。因此,能够从材料中的电子密度的知识中确定期待的或出现的相移。原则上,这种方法相对于直接测量相移的优点在于,也能唯一确定超过π的相移。在直接相位对比测量方法中,多于π的相移不能被唯一识别,因为超过整数倍的π的相移会遗失关于相移几次超过π的信息。在此,仅测量两个驻波在+/-π范围中的相位差,而不是特定的波位置的实际运行时间差。
为此,建议具有以下步骤的方法,或者说建议执行以下流程的X射线设备:
-利用两个或多个X射线光谱或X射线能量测量吸收。这作为“双能量拍摄”被足够公知,并且能够以不同方式进行。一种优选变化是使用具有两个辐射器的CT(双源CT),在其中通过X射线光谱的专门预过滤能够优化光谱分离。
-确定在断层成像测量中的检查对象中的局部电子密度,或者确定在来自光谱吸收CT图像的投影数据或者说吸收投影数据中的电子密度线积分。为此,可以使用已知的方法,例如按照参与的吸收过程的展开或基本材料分解。在此,对于临床重要的组织能够达到精度<1%。
-通过使用检查对象中的电子密度与复数折射率的实部之间的物理关系,按照以下公式计算穿透检查对象的X射线的相位信息:
在此,NA表示阿伏伽德罗常数,re表示经典电子半径,ρ表示质量密度,A表示原子质量,Z表示核电荷数,f'表示特定于原子的修正系数,λ表示X射线辐射的波长,以及δ表示相移。
特定于原子的修正系数f'对于生物对象重要的元素处于范围f'/Z<~1%,对于微量元素仅处于范围0.1%,从而以高精确度简化为:
其中ρe表示电子密度。
对于化学键,为了计算相移δ应当根据化学计量成分和键的总密度适当进行加权。由此,能够以高精确度从电子密度中按照以下公式计算对于任意能量/光谱的实部或相位图(=δ图):
原则上能够将该计算既用于投影成像又用于断层成像。在投影成像的情况下,确定电子密度的线积分,从而借助公式(3)也确定相移δ的线积分。如果将该方法用于断层成像,则通过光谱吸收确定来确定局部电子密度,其通过公式(3)导出局部相移值δ。
以上描述的方法基本上由于克拉茂-克朗尼希关系而成立,该关系描述了在完全了解折射率的虚部的能量相关性下,则也已知实部作为能量的函数。因为这在通常情况下需要了解在所有能量情况下的吸收,所以在强的X射线辐射的情况下更适合:因为吸收基本上通过两个物理效应,即光电效应和康普顿散射,促成,所以在至少两个能量或能量谱下测量吸收足以。如果在所使用的X射线能量范围中还出现通过材料(例如具有K个边的碘)的附加的吸收部分,则可以利用第三能量或第三光谱进行测量,以便改善电子密度的计算精度,并由此改善相位信息。
这里描述的方法的主要优点在于,通过这里描述的方法计算的相位图具有与吸收图相同的噪声功率谱,这获得了普遍化的CT值的量化意义。在对于临床CT典型的空间分辨率下,在相同剂量下的SNR也比利用当前可用的紧凑型PCI组件的测量更好。
根据该认知,发明人建议一种用于产生检查对象的相位对比显示的方法,其中借助利用至少两个不同的X射线能量谱确定X射线辐射的能量相关的衰减值,首先确定检查对象中的电子密度分布,然后从先前确定的电子密度分布中计算相移值,最后从计算的相移值中产生相位对比显示。
利用该方法,在第一变化中可以从电子密度沿着在焦点与探测器之间的X射线的线积分中确定电子密度分布。也就是从投影能量相关的吸收拍摄中确定在各自辐射路径上的电子密度的投影的“面占有”,即沿着分别测量的X射线的累积的电子密度,并且从中确定总相移,必要时也可以超过π的界限。从中可以产生沿着测量的穿过检查对象的X射线积分的相移的投影显示作为相位对比显示。
相比只能确定在+/-π范围的相位差的直接测量的相位对比成像方法,该测量的优点在于,也能唯一确定在π范围之外的值。由此使大于π的射线相移在重建时不会导致计算错误,并且可以在没有这样的错误下从来自不同投影方向的多个投影的相位对比显示中重建断层成像的相位对比显示。
替代地,也可以首先进行吸收数据的重建,从而能够确定局部的电子密度和其在检查对象中的分布。由此确定在检查对象中的局部的电子密度值作为电子密度分布。然后,产生在检查对象中的局部相移值的断层成像显示以用于相位对比显示。
为了确定在检查对象中的电子密度分布,例如可以射线方面在投影的图像显示中或者体素方面在断层成像的图像显示中确定在测量的衰减值中的康普顿效应的部分。
按照另一种替换方案,也可以借助基本材料分解方法确定在检查对象中的电子密度分布。在这样的材料分解方法中,确定在检查对象中出现的两个已知典型材料的分密度。如果具有材料沿着每个测量射线的分密度或者具有在检查对象中的每个体素的分密度,则也能够容易地从考察的材料的已知材料特性中确定在那里具有的电子密度。
对于确定电子密度也有利的是,使用生物检查对象,优选患者作为检查对象。在生物检查对象中,也就是在临床重要的组织中,合理地仅出现其特定于原子的修正系数f'处于范围f'/Z<1%的元素(参见公式(1)),从而导出公式(2)的简化假设特别好的成立,并且由此以好的精确度描述按照公式(3)的从电子密度到相移的过度。
相应地,发明人也建议,为了从电子密度中确定相移,也使用公式
其中δ表示相移,NA表示阿伏伽德罗常数,re表示经典电子半径,ρe表示电子密度,以及λ表示X射线辐射的波长。
本发明的范围不仅包括以上描述的方法,而且还包括一种用于对检查对象进行成像的相位对比显示的X射线***,该***具有用于控制的计算机***,其中,在计算机***的存储器中存储了至少一个程序,该程序在运行中执行以上描述的方法的方法步骤。
这样的X射线***既可以是用于产生投影的X射线图像的***,也可以是用于产生断层成像的X射线图像的***。优选地,对于执行该方法可以使用关于机械上和电学上的构造都已知的双能量CT***,其在扫描检查对象时使用两个不同的、优选尽可能少重叠的X射线能量谱。替代地,当然也可以使用具有能量选择的探测器的CT***,利用该探测器能够有针对性地确定所选择的能量范围中的吸收行为。
附图说明
以下借助附图详细描述本发明,其中仅示出了对于理解本发明所需的特征。使用了以下附图标记:1:双能量CT***;2:第一X射线管;3:第一探测器;4:第二X射线管;5:第二探测器;6:机架壳体;8:患者卧榻;9:***轴;10:计算机***;P:患者;Prg1-Prgn:计算机程序。详细示出了:
图1示出了用于执行根据本发明的方法的双能量CT***;
图2示出了通过干涉方法以生物可承受的剂量对医学模体进行的相位对比CT拍摄,
图3示出了从图2中以如图2相同剂量对模体进行的吸收CT拍摄,
图4示出了通过干涉方法以相比图2高出10倍的分辨率和高出1000倍的剂量对模体进行的相位对比CT拍摄,
图5示出了以相比图2高出10倍的分辨率和高出1000倍的剂量对模体进行的吸收CT拍摄,
图6示出了图表,其显示对于相位对比CT所需的SNR,作为结构大小的函数
图7示出了通过干涉测量方法以根据当前医学CT检查典型的分辨率对模体进行的相位对比CT拍摄,
图8示出了从图7中通过根据本发明的方法以相应于图7的分辨率对模体进行的相位对比CT拍摄。
具体实施方式
图1示出了具有机架壳体6的双能量CT***1,在该壳体中,在未详细示出的机架上具有两个辐射器-探测器***2,3和4,5,其分别具有X射线管2或4以及对置的探测器3或5。利用这两个辐射器-探测器***从患者P产生利用不同X射线能量谱的CT拍摄,为进行检查借助可沿着***轴9移动的患者卧榻8将患者移动穿过在辐射器-探测器***之间的测量区。通过具有相应程序的计算机***10进行***控制。
根据本发明,在计算机***10的存储器中还存在程序Prg1-Prgn,其在运行中执行根据本发明的方法,通过从先前确定的吸收拍摄中确定患者中的局部电子密度,例如通过基本材料分解或确定由于康普顿效应的吸收部分。然后,从电子密度中计算期待的或者在测量中发生的在穿过患者的X射线辐射路径时的相移,并且将其作为断层成像的相位对比拍摄显示、打印和/或存储以用于其它应用。
要指出的是,借助在这里的成像CT***也可以利用两个不同能量或能量谱拍摄投影照片,例如以概览扫描的形式。利用该投影照片也可以为每个辐射或每个像素确定电子占据,从中能够确定在穿过检查对象的射线路径时的总相移,优选地也超过范围π。
如果已经将从多个能量拍摄的正弦图数据(Sinogrammdaten)转换成由射线方式的电子占据组成的数据组并且将其转换成相移信息,则可以从该相移信息中重建断层成像的相位对比拍摄。
为理解本发明,利用图2和图3对应示出了利用干涉测量方法拍摄的相位对比CT拍摄(图2)以及吸收CT拍摄(图3)。两个拍摄都以对于体内CT典型的相同分辨率和相同辐射剂量建立。在此容易看出,在图2中以干涉测量方式建立的相位对比拍摄具有明显更低的SNR。
图4和图5示出了与图2和图3对应的拍摄,然而其中,具有高10倍的分辨率和高1000倍的剂量。在此容易看出,在图4中以干涉测量方式建立的相位对比拍摄具有比图5中的吸收拍摄明显更高的SNR。
图6中的图表示出了根据结构大小(横坐标)对于相位对比CT拍摄需要的SNR(纵坐标),以便根据测试对象(例如在诊断成像中的病灶)的大小到达如同吸收CT拍摄相同的探测率。
最后,利用图7和图8示出了通过干涉测量方法建立的常规的相位对比CT拍摄(图7)和按照根据本发明的方法以相同剂量产生的相同模体的相位对比CT拍摄。显然可以看出,SNR和细节上都得到显著改善。
根据本发明的方法以常规的基于吸收的成像为基础确定相位信息。以这种方式可以绕开更换到相位敏感的PCI方法需要的复杂并昂贵的技术障碍以及风险。
此外,在CT典型的分辨率中可以期待所描述的方法具有更好的剂量效率。这一方面在于相位信息的更好的SNR本身,另一方面也在于在许多PCI方法中在患者身后多达50%的X射线量子被遗失而不能用于成像的事实,由此降低了在PCI方法中的剂量效率。
与PCI相反,从双能量CT拍摄确定的相位信息的噪声纹理(噪声功率谱)与传统的CT图像相同,由此易于医生理解。
虽然在细节上通过优选实施例详细图解和描述了本发明,但是本发明不限于公开的例子,并且技术人员可以从中推出其它变化,而不脱离本发明的保护范围。
Claims (13)
1.一种用于产生检查对象(P)的相位对比显示的方法,具有以下方法步骤:
1.1通过利用至少两个不同的X射线能量谱(2,3;4,5)确定对于X射线辐射的能量相关的衰减值,来确定在检查对象(P)中的电子密度分布,
1.2从先前确定的电子密度分布中确定相移值,并且
1.3从所计算的相移值中产生相位对比显示。
2.按照权利要求1所述的方法,其特征在于,从所述电子密度沿着在焦点与探测器之间的X射线的线积分中确定所述电子密度分布。
3.按照权利要求2所述的方法,其特征在于,产生沿着穿过检查对象(P)的进行测量的X射线积分的相移的投影显示,作为相位对比显示。
4.按照权利要求1至3中任一项所述的方法,其特征在于,从来自不同投影方向的多个投影的相位对比显示中重建至少一个断层成像的相位对比显示。
5.按照权利要求1所述的方法,其特征在于,确定在检查对象(P)中的局部电子密度值,作为所述电子密度分布。
6.按照权利要求5所述的方法,其特征在于,产生在检查对象(P)中的局部相移的断层成像显示,作为相位对比显示。
7.按照权利要求1至6中任一项所述的方法,其特征在于,通过确定在所测量的衰减值中的康普顿效应的部分,进行对在检查对象(P)中的电子密度分布的确定。
8.按照权利要求1至6中任一项所述的方法,其特征在于,借助基本材料分解方法进行对在检查对象(P)中的电子密度分布的确定。
9.按照权利要求1至8中任一项所述的方法,其特征在于,使用生物检查对象,优选患者,作为检查对象。
10.按照权利要求1至9中任一项所述的方法,其特征在于,为了从电子密度中确定相移,使用公式
其中,δ表示相移,NA表示阿伏伽德罗常数,re表示经典电子半径,ρe表示电子密度,以及λ表示X射线辐射的波长。
11.一种用于对检查对象(P)进行成像的相位对比显示的X射线***,具有用于控制的计算机***(10),其特征在于,在所述计算机***(10)的存储器中存储了至少一个程序(Prg1-Prgn),所述程序在运行中执行按照权利要求1至10中任一项所述的方法。
12.按照权利要求11所述的X射线***,其特征在于,所述X射线***被构造用于产生投影的X射线图像。
13.按照权利要求11所述的X射线***,其特征在于,所述X射线***被构造用于产生断层成像的X射线图像。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102012211146.8 | 2012-06-28 | ||
DE201210211146 DE102012211146A1 (de) | 2012-06-28 | 2012-06-28 | Verfahren und Röntgensystem zur Erzeugung einer Phasenkontrastdarstellung |
PCT/EP2013/060643 WO2014000996A1 (de) | 2012-06-28 | 2013-05-23 | Verfahren und röntgensystem zur erzeugung einer phasenkontrastdarstellung |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104427938A true CN104427938A (zh) | 2015-03-18 |
Family
ID=48579017
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380036091.8A Pending CN104427938A (zh) | 2012-06-28 | 2013-05-23 | 用于产生相位对比显示的方法和x射线*** |
Country Status (4)
Country | Link |
---|---|
US (1) | US20150117595A1 (zh) |
CN (1) | CN104427938A (zh) |
DE (1) | DE102012211146A1 (zh) |
WO (1) | WO2014000996A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107427690A (zh) * | 2015-03-20 | 2017-12-01 | 皇家飞利浦有限公司 | Mrcat图像中的不确定区域的后退解决方案 |
CN112577977A (zh) * | 2019-09-30 | 2021-03-30 | 中国科学院深圳先进技术研究院 | 相衬成像方法、装置、存储介质及电子设备 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3028036A1 (en) * | 2013-07-30 | 2016-06-08 | Koninklijke Philips N.V. | Monochromatic attenuation contrast image generation by using phase contrast ct |
DE102014202745B4 (de) * | 2014-02-14 | 2023-06-01 | Siemens Healthcare Gmbh | Untersuchungsvorrichtung und Verfahren zur kombinierten Röntgen- und Ultraschallabtastung |
JP6148415B1 (ja) * | 2014-06-16 | 2017-06-14 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | コンピュータ断層撮影(ct)ハイブリッドデータ収集 |
WO2016008956A1 (en) | 2014-07-17 | 2016-01-21 | Koninklijke Philips N.V. | Iterative reconstruction method for spectral, phase-contrast imaging |
DE102016209674B4 (de) * | 2016-06-02 | 2023-10-26 | Siemens Healthcare Gmbh | Ermittlung einer räumlichen Verteilung von Materialeigenschaftswerten auf Basis einer Einzelenergie-CT-Bildaufnahme mit Hilfe eines iterativen Optimierungsverfahrens |
WO2018091344A1 (en) * | 2016-11-16 | 2018-05-24 | Koninklijke Philips N.V. | Apparatus for generating multi energy data from phase contrast imaging data |
WO2019056309A1 (en) * | 2017-09-22 | 2019-03-28 | Shenzhen United Imaging Healthcare Co., Ltd. | METHOD AND SYSTEM FOR GENERATING A PHASE CONTRAST IMAGE |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1247687A (zh) * | 1996-12-24 | 2000-03-15 | X-射线技术股份有限公司 | 在相位衬比成象中的相位恢复 |
US6563906B2 (en) * | 2000-08-28 | 2003-05-13 | University Of New Brunswick | X-ray compton scattering density measurement at a point within an object |
WO2009058976A1 (en) * | 2007-10-30 | 2009-05-07 | Massachusetts Institute Of Technology | Phase-contrast x-ray imaging |
CN101900696A (zh) * | 2009-05-27 | 2010-12-01 | 清华大学 | 双能欠采样物质识别方法和*** |
CN101960298A (zh) * | 2008-03-12 | 2011-01-26 | 佳能株式会社 | X射线成像装置、x射线成像方法和x射线成像装置的控制方法 |
JP2011206188A (ja) * | 2010-03-29 | 2011-10-20 | Fujifilm Corp | 放射線撮影システム及び方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7286640B2 (en) * | 2004-04-09 | 2007-10-23 | Xradia, Inc. | Dual-band detector system for x-ray imaging of biological samples |
WO2007110795A2 (en) * | 2006-03-29 | 2007-10-04 | Philips Intellectual Property & Standards Gmbh | Effective dual-energy x-ray attenuation measurement |
US8995609B2 (en) * | 2011-08-02 | 2015-03-31 | Georgia Tech Research Corporation | X-ray compton scatter imaging on volumetric CT systems |
-
2012
- 2012-06-28 DE DE201210211146 patent/DE102012211146A1/de not_active Withdrawn
-
2013
- 2013-05-23 CN CN201380036091.8A patent/CN104427938A/zh active Pending
- 2013-05-23 US US14/408,314 patent/US20150117595A1/en not_active Abandoned
- 2013-05-23 WO PCT/EP2013/060643 patent/WO2014000996A1/de active Application Filing
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1247687A (zh) * | 1996-12-24 | 2000-03-15 | X-射线技术股份有限公司 | 在相位衬比成象中的相位恢复 |
US6563906B2 (en) * | 2000-08-28 | 2003-05-13 | University Of New Brunswick | X-ray compton scattering density measurement at a point within an object |
WO2009058976A1 (en) * | 2007-10-30 | 2009-05-07 | Massachusetts Institute Of Technology | Phase-contrast x-ray imaging |
CN101960298A (zh) * | 2008-03-12 | 2011-01-26 | 佳能株式会社 | X射线成像装置、x射线成像方法和x射线成像装置的控制方法 |
CN101900696A (zh) * | 2009-05-27 | 2010-12-01 | 清华大学 | 双能欠采样物质识别方法和*** |
JP2011206188A (ja) * | 2010-03-29 | 2011-10-20 | Fujifilm Corp | 放射線撮影システム及び方法 |
Non-Patent Citations (2)
Title |
---|
ANASTASIYA KHROMOVA: "A Study of Multiple Refractive Scattering of Monoenergetic X-Rays from Ensembles of Monodisperse Spheres", 《PHYSICS:MEDICAL PHYSICS》 * |
T.E.GUREYEV AND K.A.NUGENT: "Phase retrieval with the transport-of-intensity equation. II. Orthogonal series solution for nonuniform illumination", 《JOURNAL OF THE OPTICAL SOCIETY OF AMERICA》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107427690A (zh) * | 2015-03-20 | 2017-12-01 | 皇家飞利浦有限公司 | Mrcat图像中的不确定区域的后退解决方案 |
CN107427690B (zh) * | 2015-03-20 | 2020-09-01 | 皇家飞利浦有限公司 | Mrcat图像中的不确定区域的后退解决方案 |
CN112577977A (zh) * | 2019-09-30 | 2021-03-30 | 中国科学院深圳先进技术研究院 | 相衬成像方法、装置、存储介质及电子设备 |
CN112577977B (zh) * | 2019-09-30 | 2023-11-03 | 中国科学院深圳先进技术研究院 | 相衬成像方法、装置、存储介质及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
WO2014000996A1 (de) | 2014-01-03 |
US20150117595A1 (en) | 2015-04-30 |
DE102012211146A1 (de) | 2014-01-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104427938A (zh) | 用于产生相位对比显示的方法和x射线*** | |
JP5142540B2 (ja) | X線装置の焦点‐検出器システム | |
CN101011252B (zh) | 用于产生相位对比照片的x射线设备的焦点/检测器*** | |
US9842414B2 (en) | Monochromatic attenuation contrast image generation by using phase contrast CT | |
CN101873828B (zh) | 用于x射线相位对比成像的探测装置 | |
US7643866B2 (en) | Method for producing a computed tomography display of tissue structures by applying a contrast medium | |
Chabior et al. | Beam hardening effects in grating‐based x‐ray phase‐contrast imaging | |
CN1880949B (zh) | 计算吸收物质加权系数和改善图像对比度-噪声比的方法 | |
Garnett | A comprehensive review of dual-energy and multi-spectral computed tomography | |
Raupach et al. | Performance evaluation of x‐ray differential phase contrast computed tomography (PCT) with respect to medical imaging | |
US9155510B2 (en) | Systems and methods for generating x-ray phase contrast images using a conventional x-ray imaging system | |
Mackenzie et al. | Image simulation and a model of noise power spectra across a range of mammographic beam qualities | |
JP2007203064A (ja) | X線装置の焦点‐検出器装置 | |
Abadi et al. | Direct quantification of breast dose during coronary CT angiography and evaluation of dose reduction strategies | |
US11234663B2 (en) | Apparatus for generating multi energy data from phase contrast imaging data | |
Li et al. | Fundamental relationship between the noise properties of grating‐based differential phase contrast CT and absorption CT: Theoretical framework using a cascaded system model and experimental validation | |
Liu et al. | Synchrotron‐based high‐energy x‐ray phase sensitive microtomography for biomedical research | |
JP6148415B1 (ja) | コンピュータ断層撮影(ct)ハイブリッドデータ収集 | |
Fingerle et al. | Simulated cystic renal lesions: quantitative X-ray phase-contrast CT—an in vitro phantom study | |
JP4472492B2 (ja) | 被検体における元素濃度を位置分解して求める方法および装置 | |
Saito | Dual‐energy approach to contrast‐enhanced mammography using the balanced filter method: Spectral optimization and preliminary phantom measurement | |
Basharat et al. | Theoretical feasibility of dual‐energy radiography for structural and functional imaging of chronic obstructive pulmonary disease | |
Xu et al. | Intensity-based iterative reconstruction for helical grating interferometry breast CT with static grating configuration | |
CN111089869B (zh) | 多能探测器x射线相衬成像方法及***、储存介质及设备 | |
Patera et al. | X-ray grating interferometry design for the 4D GRAPH-X system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150318 |