CN106056642A - 多谱式计算机断层成像 - Google Patents
多谱式计算机断层成像 Download PDFInfo
- Publication number
- CN106056642A CN106056642A CN201610220577.2A CN201610220577A CN106056642A CN 106056642 A CN106056642 A CN 106056642A CN 201610220577 A CN201610220577 A CN 201610220577A CN 106056642 A CN106056642 A CN 106056642A
- Authority
- CN
- China
- Prior art keywords
- ext
- spectrum
- kext
- value
- picture point
- 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
- 238000013170 computed tomography imaging Methods 0.000 title claims abstract description 16
- 238000001228 spectrum Methods 0.000 claims abstract description 82
- 238000000034 method Methods 0.000 claims abstract description 46
- 238000005259 measurement Methods 0.000 claims description 41
- 238000013016 damping Methods 0.000 claims description 25
- 230000015572 biosynthetic process Effects 0.000 claims description 14
- 239000000463 material Substances 0.000 claims description 12
- 239000000203 mixture Substances 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 8
- 239000003086 colorant Substances 0.000 claims description 4
- 230000007423 decrease Effects 0.000 claims description 2
- 238000003384 imaging method Methods 0.000 claims description 2
- 238000003325 tomography Methods 0.000 claims 1
- 230000000875 corresponding effect Effects 0.000 description 21
- 230000001276 controlling effect Effects 0.000 description 9
- 230000003595 spectral effect Effects 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000007689 inspection Methods 0.000 description 5
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 239000000193 iodinated contrast media Substances 0.000 description 2
- 229910052740 iodine Inorganic materials 0.000 description 2
- 239000011630 iodine Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- ODINCKMPIJJUCX-UHFFFAOYSA-N Calcium oxide Chemical compound [Ca]=O ODINCKMPIJJUCX-UHFFFAOYSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 230000002308 calcification Effects 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- 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/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- 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]
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/08—Volume rendering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/408—Dual energy
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Computer Graphics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明描述了一种多谱式成像方法(100,200)、优选一种CT成像方法。在多谱式成像方法(100,200)中,采集检查对象(O)的待拍摄区域(VOI)的谱分辨的投影测量数据(PMD(sn(f))),其与多个预先确定的不同的分谱(sn(f))相关联。重建谱分辨的图像数据(BD(sn(f)),其针对待拍摄区域(VOI)的每个图像点(xi)包括多个衰减值(μn(xi,sn(f))),其中,所述衰减值(μn(xi,sn(f)))分别与预先确定的不同的分谱(sn(f))中的一个相关联。此外,基于所述多个衰减值(μn(xi,sn(f)))针对每个图像点(xi)确定极限衰减值(μext(xi,sn_ext(f))),并且确定代表性图像数据集合(BDR),其中,对每个图像点(xi)分配所确定的极限衰减值(μext(xi,sn_ext(f)))。此外,描述了一种图像数据确定装置(40)。另外,描述了一种计算机断层成像系统(1)。
Description
技术领域
本发明涉及一种多谱式成像方法。此外,本发明涉及一种图像数据确定装置。另外,本发明还涉及一种计算机断层成像系统。
背景技术
借助现代的成像方法,经常产生二维或者三维的图像数据,其可以用于将拍摄的检查对象可视化,此外还可以用于其它应用。
这些成像方法经常基于X射线的采集,其中,产生所谓的投影测量数据。例如,可以借助计算机断层成像系统(CT系统)来获取投影测量数据。在CT系统中,通常布置在机架上的X光源与相对地布置的X光检测器的组合围绕检查对象(下面不失一般性将其称为患者)所在的测量空间运行。在此,旋转中心(也称为“对称中心”)与所谓的系统轴z一致。在一次或更多次旋转期间,用X光源的X射线照射患者,其中,借助相对的X光检测器采集投影测量数据或X光投影数据。
所产生的投影测量数据、也简称为投影数据,特别地与X光检测器的结构有关。X光检测器通常具有通常以规则的像素阵列的形式布置的多个检测单元。检测单元分别针对射向检测单元的X射线产生检测信号,其在特定时刻关于X射线的强度和谱分布被分析,以便得出关于检查对象的结论并且产生投影测量数据。
传统上,CT图像数据具有带有灰度值的图像点,灰度值分别对应于待拍摄的区域在每个时刻的测量的衰减值。下面在提到图像点时,由此例如意为二维像素和三维体素。测量的衰减值是待拍摄的区域中的对应的对象点对所入射的X射线已经历的总衰减有贡献的量值。该衰减值以亨斯菲尔德单位(Hounsfield-Unit)给出。在此,值-1000HU对应于由于空气而产生的衰减,0HU相当于与水的衰减对应的衰减。通常,图像数据作为灰度值编码的二维图像来显示。
此外,为了考虑待拍摄的体积区域的第三个维度,存在不同的显示可能性。例如,并排显示多个二维图像。在此,每个二维图像表示三维体积的一个层。也可以先后显示二维图像,从而可以利用鼠标或者其它进行控制的输入装置翻阅所拍摄的体积。为了形成各个层图像,例如可以对与一个层相关联的三维图像数据材料进行重整。在此,例如在层方向上、即在与层平面垂直的方向上对图像数据求平均值,并且在与该层相关联的二维图像中显示所计算的平均值。这种重整也称为多平面重整(Multi-Planar-Reformatierung,MPR)。替换地,也可以在与层平面垂直的方向上进行最大值或最小值的计算,并且在与该层相关联的二维图像中显示所计算的最大值或最小值。在显示最大值的情况下,结合这种过程,也称为最大强度值的投影(maximumintensity projection=MIP)。对于显示最小值的情况,结合这种过程,也称为最小强度值的投影(minimum intensity projection=MinIP)。
在CT系统的一系列应用中,从涉及同一对象的独立的测量中采集多个数据集合或投影数据集合。这些数据集合例如在采用利用多个能量阈值来记录、即所谓的多能量扫描的情况下被记录。在多能量扫描中,利用一个或更多个能量阈值采集量子计数检测器的数据,其中,被能量阈值分离的各个能量区域与不同的数据集合相关联。
在这种情况下,对于待拍摄的区域的每个点,存在多个衰减值,其分别与不同的数据集合中的一个相关联,并且合起来描述入射的X射线的衰减值的谱相关性。原则上,谱相关性可以通过针对不同的谱或谱间隔给出不同的衰减值来表征,或者也可以通过给出各个测量点的不同基材的分量来表征。在此,基材与分谱相关联,所述分谱一起给出图像数据的谱相关性的总谱。也就是说,图像数据的谱相关性对应于在测量中产生的数据的第四个维度。然而,很难直观地显示四维的测量数据。
传统上,将这些四维数据分割为分别与一个谱分量对应的三维体积数据,并且单独显示根据谱分量分离的三维数据。然而,四维测量数据的这种显示非常不清楚并且复杂。
发明内容
由此,本发明要解决的技术问题是,在多谱式成像中产生能够清楚地显示、同时包括重要的谱图像信息的图像数据集合。
上述技术问题通过根据本发明的多谱式成像方法、根据本发明的图像数据确定装置和根据本发明的计算机断层成像系统来解决。
在优选可以构造为多谱式CT成像方法的根据本发明的多谱式成像方法中,采集检查对象的待拍摄区域的谱分辨的投影测量数据。投影测量数据与多个预先确定的不同的分谱相关联。分谱例如可以对应于通过X光检测器的能量阈值定义的频率间隔。分谱也可以通过与各个基材相关联的谱分量来定义。此外,在根据本发明的多谱式成像方法中,根据所采集的投影测量数据重建谱分辨的图像数据。这些图像数据包括针对待拍摄区域的每个图像点的多个衰减值,其中,图像数据分别与预先确定的不同的分谱中的一个相关联。也就是说,每个图像点不仅与一个衰减值相关联,而且其分别取决于所述分谱而与多个衰减值相关联。此外,在根据本发明的方法中,基于多个衰减值针对每个图像点确定极限衰减值。极限衰减值例如可以是最小值或者最大值,也就是说,与一个图像点相关联的谱衰减值中的、显示最小衰减或最大衰减的值的谱衰减值。最后,产生代表性的图像数据集合,其中,对每个图像点分配分别关于各个图像点确定的极限衰减值。这种图像数据集合提供关于在待拍摄的区域中存在的待拍摄结构的压缩信息。例如,如果相应地选择衰减值的最大值作为极限值,则产生具有良好的信噪比的图像。此外,可以假设,一个图像点的相应地最大的衰减值表示该位置处的主要材料以及还最好地再现结构特征。
特别是用于多谱式CT成像的根据本发明的图像数据确定装置具有用于从检查对象的待拍摄区域中采集谱分辨的投影测量数据的输入接口,其中,投影测量数据与多个预先确定的不同的分谱相关联。根据本发明的图像数据确定装置还具有用于重建谱分辨的图像数据的重建单元,所述谱分辨的图像数据包括针对待拍摄区域的每个图像点的与预先确定的分谱相关联的多个衰减值。用于确定每个图像点的极限衰减值的衰减值确定单元也是根据本发明的图像数据确定装置的一部分。此外,根据本发明的图像数据确定装置包括用于产生代表性图像数据集合的图像数据集合产生单元,其中,对每个图像点分配所确定的极限衰减值。
根据本发明的计算机断层成像系统具有根据本发明的图像数据确定装置。
根据本发明的图像数据确定装置的主要部件大部分可以以软件部件的形式构造。这特别涉及重建单元、衰减值确定单元和图像数据集合产生单元。然而,原则上,这些部件部分地、特别是在涉及特别快速的计算时也可以以软件支持的硬件、例如FPGA等的形式来实现。所需的接口例如在仅涉及来自其它软件部件的数据的传输时同样可以作为软件接口来构造。但是其也可以作为由合适的软件控制的、按照硬件构建的接口来构造。
特别地,根据本发明的图像数据确定装置可以是CT系统的用户终端或者控制装置的一部分。
在很大程度上按照软件实现具有如下优点:还可以以简单的方式通过软件更新来升级迄今为止使用的控制装置,以便以根据本发明的方式工作。在这方面,上述技术问题还通过对应的具有计算机程序的计算机程序产品来解决,其能够以程序段直接加载到计算机断层成像系统的控制装置的存储装置中,以便在控制装置中执行程序时,执行根据本发明的方法的所有步骤。除了计算机程序之外,这种计算机程序产品在必要时还可以包括附加组成部分、例如文档和/或附加部件、还有硬件部件、例如用于使用软件的硬件密钥(软件狗等)。
为了向控制装置传输和/或为了存储在控制装置上或中,可以使用计算机可读介质、例如存储棒、硬盘或其它便携式或者固定地构建的数据载体,其存储有控制装置的运算单元可读取并且可执行的计算机程序的程序段。为此,运算单元例如可以具有一个或更多个协作的微处理器等。
从属权利要求以及下面的描述相应地包含本发明的特别有利的构造和展开。在此,特别地也可以将一种权利要求类别的权利要求与另一种权利要求类别的从属权利要求类似的展开。此外,在本发明的范围内,还可以将不同的实施例和权利要求的不同特征与新的实施例组合。
在根据本发明的一个构造中,将分谱与不同的频率间隔相关联。如已经提及的那样,检测器利用多个能量阈值按照频率间隔分离地采集X射线。根据采集的计数速率,对这些单个的频率间隔分配不同的衰减值。所采集的衰减值形成谱,谱也可以分解成分别与基材相关联的分谱,其中,在分解中确定的系数表示关于在该图像点出现的基材的份额的信息。
在根据本发明的方法的一个有利构造中,极限衰减值是最大值。其意为对共同像素点分配的多个衰减值中的最大衰减值。
替换地,极限衰减值也可以是最小值。其意为对共同图像点分配的多个衰减值中的最小衰减值。
作为最小衰减值的显示的示例,在利用含碘造影剂进行拍摄时,通过关于所有阈值或分谱选择“最小”衰减值可以获得基本上“无碘”的图像,从而该显示近似于本征扫描(natives Scan)。
在选择“最大”衰减值时,对于使用含碘造影剂的情况,获得最大碘对比,由此获得最佳对比图像。
替换地,利用“最小”衰减值的显示适合于对血管中的钙的显示,因为由于石灰或骨头固有的高对比度,作为效果可能出现在文献中描述的盛开(Blooming),其在这种显示中同样减小,并且能够更好地判断钙化的真实程度。
在根据本发明的方法的一个特别有效的变形方案中,在产生代表性图像数据集合时,附加地对每个图像点分配关于与所确定的极限衰减值相关联的极限分谱的信息。换句话说,将代表性图像数据集合与关于极限衰减值的谱来源的信息相关联。
在根据本发明的方法的一个替换变形方案中,对每个分谱分配谱值。谱值应当理解为根据其值能够推断相应的衰减值的谱来源的变量。该谱值例如可以与具有特定KV值的检测器的能量阈值相关联。其也可以与特定基材相关联。
特别优选的是,关于与所确定的极限衰减值相关联的极限分谱的附加信息包括与相应地相关联的极限分谱相关联的极限谱值。在该变形方案中,极限谱值给出关于所确定的极限衰减值的谱来源的信息。
在根据本发明的方法的一个能够特别简单地实现的构造中,代表性图像数据集合被构造为作为灰度级值显示与每个图像点相关联的极限衰减值。
在根据本发明的方法的一个特别有效的构造中,代表性图像数据集合被构造为,利用其强度取决于相应的极限衰减值的大小的不同颜色,根据相关联的极限谱值,显示与每个图像点相关联的极限衰减值。
在根据本发明的方法的一个替换构造中,针对每个图像点确定附加的第二高或第二小的谱衰减值或者另外的附加谱衰减值。此外,在产生代表性图像数据集合的步骤中,附加地将每个图像点与所确定的附加谱衰减值和其相关联的分谱和/或与分谱对应的谱值相关联。
代表性图像数据集合例如可以以如下方式显示:利用强度与相应的谱衰减值的大小有关的不同的颜色,与相应的相关联的谱值有关地显示与共同的图像点相关联的极限谱衰减值和附加谱衰减值,并且与该共同的图像点相关联的极限谱衰减值和附加谱衰减值共同作为颜色的颜色混合来显示,所述颜色和与该共同的图像点相关联的极限谱衰减值和附加谱衰减值相关联。换句话说,在该构造中,通过颜色混合,通过不同的颜色显示一个图像点中的衰减的谱分量,所述颜色共同得到颜色混合。在此,各个颜色的强度由相关联的衰减值的大小确定。以这种方式,图像数据状态得到一个附加维度,该附加维度说明关于每个图像点中的谱分量的信息。
优选,与各个谱衰减值相关联的颜色值包括基本色。
还可以修改根据本发明的方法,使得将所产生的代表性图像数据集合划分为作为二维图像数据显示的多个层,其中,二维显示作为空间MIP或MinIP或MPR显示来进行。
附图说明
下面,参考附图根据实施例再一次对本发明进行详细说明。
其中:
图1示出了图示根据本发明的第一实施例的多谱式成像方法的流程图,
图2示出了图示根据本发明的第二实施例的多谱式成像方法的流程图,
图3示出了图示根据本发明的第三实施例的多谱式成像方法的流程图,
图4示出了根据本发明的实施例的图像数据确定装置的框图,
图5示出了根据本发明的实施例的计算机断层成像系统的示意图。
具体实施方式
在图1中示出了根据本发明的第一实施例的多谱式成像方法100。首先,在方法100的步骤1.I中,拍摄检查对象O(参见图5)的待拍摄区域VOI的谱分辨的投影测量数据PMD(sn(f))。为此,向CT系统1(参见图5)的测量单元传输采集控制信号并且采集投影测量数据PMD。投影测量数据例如可以借助具有多个能量阈值的CT检测器来拍摄。这些投影测量数据例如包括与不同的分谱sn(f)相关联的多个独立的投影测量数据集合PMDn。这些分谱sn(f)例如可以是简单地通过至少一个阈值划分的分谱。也就是说,在这种情况下,每个分谱sn(f)与总谱的一个频段相关联。但是,分谱也可以分解成与所谓的基材相关联的分谱,其叠加得到所采集的投影测量数据的总谱。在这种情况下,代替将投影测量数据集合PMDn与频段相关联,可以将投影测量数据集合与各个基材相关联。这两种类型的投影测量数据PMD的谱分解都应当包含在分解为投影测量数据集合PMDn的谱分解的形成中。
现在,在步骤1.II中,重建谱分辨的图像数据BD(sn(f))。在这种情境中,重建应当理解为通常的例如借助滤波反投影将独立的投影测量数据集合PMDn重建为图像数据的图像数据重建。可以将图像数据BD划分为谱分辨的图像数据集合BDn,其包括相应地与独立的投影测量数据集合PMDn中的一个相关联的、与位置有关的、即分别与一个图像点xi相关联的衰减值(μn(xi,sn(f)))。也就是说,换句话说,在划分为N个分谱时(例如当检测器包括N-1个能量阈值时),每个图像点xi与N个在谱方面不同的衰减值(μn(xi,sn(f)))的集合相关联。
现在,在步骤1.III中,关于每个图像点xi,确定在谱方面不同的衰减值(μn(xi,sn(f)))的极值(μext(xi,sn_ext(f)))。例如,这可以是谱方面不同的衰减值μn(xi,sn(f))的集合的衰减值,该衰减值是其中具有最大值的。
在步骤1.IV中,产生代表性图像数据集合BDR,其作为衰减值包括在步骤1.III中确定的极值μext(xi,sn_ext(f))。
在步骤1.V中,作为灰度值GSW以图形表示所确定的图像数据,其各自的灰度级对应于分别相关联的极值μext(xi,sn_ext(f))。所选择的图形表示也可以直观地表示为谱方向上的MIP或MinIP。
在图2中示出了根据本发明的第二实施例的方法200。在方法200中,除了已经在第一实施例中选择的作为谱方向上的MIP或MinIP的表示之外,还可以以图形表示关于相应的衰减值的谱来源的信息。以这种方式,图像数据包含一种类型的附加的谱维度。
步骤2.I至2.IV对应于步骤1.I至1.IV,因此参考结合图1对此的说明。在步骤2.V中,除了极限衰减值之外,还将关于极值的谱来源的信息sn_ext加入图像数据BDR中。此外,该谱信息sn_ext由谱值vsn_ext表征。在步骤2.VI中,将谱值vsn_ext分别与色值FW相关联,其中,色值FW例如对应于色谱中的颜色C,其与投影测量数据集合PMDn中的一个对应或相关联,所述投影测量数据集合与不同的分谱sn(f)相关联。也就是说,利用该颜色C表示图像数据的来源,其中,将与分别与一个分谱或谱段相关联的投影测量数据集合PMDn中的一个的关联理解为来源。在步骤2.VII中,与所确定的色值FW和代表性图像数据BDR对应地,作为彩色图像数据C以图形表示所确定的图像数据。在表示中,颜色的强度对应于相应地相关联的谱极值μext(xi,sn_ext(f)),并且颜色对应于相应的谱来源。也可以作为谱方向上的MIP或MinIP直观地表示所选择的图形表示。
在图3中示出了根据本发明的第三实施例的多谱式成像方法300。步骤3.I至3.II对应于步骤1.I至1.II或2.I至2.II,因此关于其描述参考结合图1的说明。然而,与在前面描述的实施例100、200中不同,现在,在步骤3.III中,在产生代表性图像数据集合BDR时,不仅确定衰减值的数据集合的极值μext(xi,sn_ext(f)),还确定分别与具有其它分谱sn_2ext(f),sn_kext(f)的其它投影测量数据集合PMDn相关联的另外的附加的谱衰减值μ2ext(xi,sn_2ext(f)),μkext(xi,sn_kext(f))。在步骤3.IV中,附加地将每个图像点xi与所确定的附加的谱衰减值μ2ext(xi,sn_2ext(f)),μkext(xi,sn_kext(f))相关联。在步骤3.V中,对各个图像点xi分配其相关联的分谱sn_2ext(f),sn_kext(f)以及与分谱sn_2ext(f),sn_kext(f)对应的谱值vs2_kext,vsn_kext。在步骤3.VI中,与所确定的谱值VSn_ext,VS2_kext,VSn_kext有关地对各个图像点分配色值FW。因为现在对每个图像点xi分配了多个色值FW,因此对于每个图像点xi形成了与所确定的色值FW对应的各个颜色C的颜色混合CM。在步骤3.VII中,与所确定的代表性图像数据BDR和相关联的颜色混合CM对应地以图形表示所确定的图像数据。在表示时,颜色混合CM中的颜色C的强度对应于分别相关联的衰减值μ2ext(xi,sn_2ext(f)),μkext(xi,sn_kext(f)),μext(xi,sn_ext(f))的大小,并且颜色混合CM中的颜色C对应于相关联的衰减值(μ2ext(xi,sn_2ext(f)),μkext(xi,sn_kext(f)),μext(xi,sn_ext(f))的相应的谱来源。所描述的方法优选自动执行。
对此,在图4中示出了根据本发明的实施例的图像数据确定装置40。图像数据确定装置40例如可以是在图5中示出的CT系统1的控制装置20的一部分。图像数据确定装置40包括从CT系统1(参见图5)的测量装置(参见图5)接收投影数据PMD的输入接口41。投影数据PMD是检查对象O的待拍摄区域VOI的谱分辨的投影测量数据PMD(sn(f)),其与多个预先确定的不同的分谱sn(f)相关联。投影测量数据PMD(sn(f))被传输到重建谱分辨的图像数据BD(sn(f))的重建单元42,所述谱分辨的图像数据针对待拍摄区域VOI的每个图像点xi,包括多个与预先确定的分谱sn(f)相关联的衰减值μn(xi,sn(f))。重建的图像数据BD(sn(f))被传输到针对每个图像点xi确定至少一个极限衰减值μext(xi,sn_ext(f))的衰减值确定单元43。随后,向图像数据集合产生单元44传输极限衰减值μext(xi,sn_ext(f)),图像数据集合产生单元44由此产生代表性图像数据集合BDR,其中,对每个图像点xi分配所确定的极限衰减值μext(xi,sn_ext(f))。
随后,代表性图像数据集合BDR传输到输出接口45。输出接口45将代表性图像数据集合BDR输出到其它单元、例如存储单元或图像显示单元、例如终端或屏幕。
在图5中示意性地示出了根据本发明的实施例的具有根据本发明的图像数据确定装置40的计算机断层成像系统(CT系统)1。
在此,CT系统1主要由扫描器10构成,其中,在机架11上,具有检测器16和与检测器16相对的X光源15的投影数据获取单元5围绕测量空间12旋转。患者支承装置3或患者台3位于扫描器10前面,患者台3的上部2可以与位于其上的患者O一起移入扫描器10,以便使患者O相对于检测器系统16移动通过测量空间12。扫描器10和患者台3由控制装置20控制,从控制装置20经由通常的控制接口23获得采集控制信号AS,以便以传统方式根据预先给定的测量协议控制整个系统。通过沿着与纵向通过测量空间12的系统轴z对应的z方向移动患者O以及X光源15同时旋转,针对X光源15相对于患者O在测量期间得到螺旋轨迹。在此,检测器16总是相对于X光源15一起并行地运行,以便采集投影测量数据PMD,然后将其用于重建体积和/或层图像数据。
也可以执行连续测量方法,其中,驶向z方向上的固定位置,然后在一次旋转、一次部分旋转或者多次旋转期间,在相关z位置采集所需的投影测量数据PMD,以便重建该z位置的片段图像,或者以便根据多个z位置的投影数据重建体积图像数据。
根据本发明的方法原则上也可用于例如具有多个X光源和/或检测器和/或具有形成一个完整的环的检测器的其它CT系统。
将由检测器16获取的投影测量数据PMD(下面也称为原始数据)经由在本实施例中是图像数据确定装置40的输入接口41的原始数据接口,传输到控制装置20或包含于其中的图像数据确定装置40。然后,在图像数据确定装置40中以上面描述的方式对该原始数据进行处理。在本实施例中,图像数据确定装置40在控制装置20中在处理器上以软件的形式实现。
在图像数据确定装置40中进行处理之后,将所确定的代表性图像数据BDR输出到CT系统的控制装置20的存储单元22以及例如输出单元。
最后,应当再一次指出,前面描述的方法和装置仅仅是本发明的优选实施例,本领域技术人员可以对本发明进行变形,而不脱离由权利要求给出的本发明的范围。本发明不局限于应用于医疗领域,相反本发明原则上还可以应用于为了其它目的、例如为了进行材料检测等而进行的CT图像拍摄。为了完整起见,还应当指出,不定冠词“一”或“一个”的使用不是排他的,相关特征也可以以复数的形式存在。类似地,术语“单元”不排除其由在必要时也可能分布在空间中的多个部件构成。
Claims (15)
1.一种优选用于多谱式CT成像的多谱式成像方法(100,200),具有以下步骤:
-采集检查对象(O)的待成像区域(VOI)的谱分辨的投影测量数据(PMD(sn(f))),其与多个预先确定的不同的分谱(sn(f))相关联,
-重建谱分辨的图像数据(BD(sn(f)),其针对待成像区域(VOI)的每个图像点(xi)包括多个衰减值(μn(xi,sn(f))),其中,所述衰减值(μn(xi,sn(f)))分别与预先确定的不同的分谱(sn(f))中的一个相关联,
-基于所述多个衰减值(μn(xi,sn(f)))针对每个图像点(xi)确定极限衰减值(μext(xi,sn_ext(f))),
-产生代表性图像数据集合(BDR),其中,对每个图像点(xi)分配所确定的极限衰减值(μext(xi,sn_ext(f)))。
2.根据权利要求1所述的方法(100,200),其中,所述分谱(sn(f))与不同的频率间隔(If)或者不同的基材(BMn)相关联。
3.根据权利要求1或2所述的方法(100,200),其中,所述极限衰减值(μext(xi,sn_ext(f)))是最大值或最小值。
4.根据权利要求1至3中任一项所述的方法(100,200),其中,在产生代表性图像数据集合(BDR)时,附加地对每个图像点(xi)分配关于与所确定的极限衰减值(μext(xi,sn_ext(f)))相关联的极限分谱(sn_ext(f))的信息。
5.根据权利要求1至4中任一项所述的方法(100,200),其中,对分谱(sn(f))中的每一个分配谱值(VSn)。
6.根据权利要求4或5所述的方法(100,200),其中,关于与所确定的极限衰减值(μext(xi,sn_ext(f)))相关联的极限分谱(sn_ext(f))的附加信息包括与分别相关联的极限分谱(sn_ext(f))相关联的极限谱值(VSn_ext)。
7.根据权利要求1至6中任一项所述的方法(100,200),其中,以作为灰度级值显示与每个图像点(xi)相关联的极限衰减值(μext(xi,sn_ext(f)))的方式,显示代表性图像数据集合(BDR)。
8.根据权利要求4至7中任一项所述的方法(100,200),其中,以与相关联的极限谱值(VSn_ext)有关地利用不同的颜色(C)显示与每个图像点(xi)相关联的极限衰减值(μext(xi,sn_ext(f)))的方式,显示代表性图像数据集合(BDR),其中,颜色(C)的强度(I)与相应的极限衰减值(μext(xi,sn_ext(f)))的大小有关。
9.根据权利要求4至8中任一项所述的方法(200),其中,针对每个图像点(xi)确定附加的第二高或第二小的谱衰减值(μ2ext(xi,sn_2ext(f)))或者另外的附加谱衰减值(μkext(xi,sn_kext(f))),并且在产生代表性图像数据集合(BDR)的步骤中,附加地将每个图像点(xi)与所确定的附加谱衰减值(μ2ext(xi,sn_2ext(f)),μkext(xi,sn_kext(f)))和其相关联的分谱(sn_2ext(f),sn_kext(f))和/或与分谱(sn_2ext(f),sn_kext(f))对应的谱值(VS2_kext,VSn_kext)相关联。
10.根据权利要求9所述的方法(200),其中,代表性图像数据集合(BDR)以如下方式被显示:与分别相关联的谱值(VSn_ext,VSn_2ext,VSn_kext)有关地利用不同的颜色(C)显示与共同的图像点(xi)相关联的极限谱衰减值和附加谱衰减值(μext(xi,sn_ext(f)),μkext(xi,sn_kext(f))),其中,不同的颜色(C)的强度(I)与相应的谱衰减值(μext(xi,sn_ext(f)),μkext(xi,sn_2ext(f)),μkext(xi,sn_kext(f)))的大小有关,并且与该共同的图像点(xi)相关联的极限谱衰减值和附加谱衰减值(μext(xi,sn_ext(f)),μkext(xi,sn_kext(f)))共同作为颜色(C)的颜色混合(CM)来显示,所述颜色和与该共同的图像点(xi)相关联的极限谱衰减值和附加谱衰减值(μext(xi,sn_ext(f)),μkext(xi,sn_kext(f)))相关联。
11.根据权利要求1至10中任一项所述的方法(100,200),其中,将所产生的代表性图像数据集合(BDR)划分为多个层,所述多个层作为二维图像数据被显示,其中,二维显示作为空间上的MIP或MinIP或MPR显示来进行。
12.一种特别用于多谱式CT成像的图像数据确定装置(40),具有:
-输入接口(41),用于从检查对象(O)的待成像区域(VOI)中采集谱分辨的投影测量数据(PMD(sn(f))),所述投影测量数据与多个预先确定的不同的分谱(sn(f))相关联,
-重建单元(42),用于重建谱分辨的图像数据BD(sn(f)),所述谱分辨的图像数据针对待成像区域(VOI)的每个图像点(xi)包括多个与预先确定的分谱(sn(f))相关联的衰减值(μn(xi,sn(f))),
-衰减值确定单元(43),用于针对每个图像点(xi)确定极限衰减值(μext(xi,sn_ext(f))),
-图像数据集合产生单元(44),用于产生代表性图像数据集合(BDR),其中,对每个图像点(xi)分配所确定的极限衰减值(μext(xi,sn_ext(f)))。
13.一种计算机断层成像系统(1),包括具有根据权利要求12所述的图像数据确定装置(40)的控制装置(20)。
14.一种具有计算机程序的计算机程序产品,其能够以程序段直接加载到计算机断层成像系统(1)的控制装置的存储装置中,以便在计算机断层成像系统(1)的控制装置中执行所述计算机程序时,执行根据权利要求1至11中任一项所述的方法的所有步骤。
15.一种计算机可读介质,其上存储有运算单元能读取并且能执行的程序段,以便在由所述运算单元执行所述计算机段时,执行根据权利要求1至11中任一项所述的方法的所有步骤。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102015206630.4A DE102015206630B4 (de) | 2015-04-14 | 2015-04-14 | Multispektrale CT-Bildgebung |
DE102015206630.4 | 2015-04-14 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106056642A true CN106056642A (zh) | 2016-10-26 |
CN106056642B CN106056642B (zh) | 2019-09-06 |
Family
ID=57043203
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610220577.2A Expired - Fee Related CN106056642B (zh) | 2015-04-14 | 2016-04-11 | 多谱式计算机断层成像 |
Country Status (3)
Country | Link |
---|---|
US (1) | US10290122B2 (zh) |
CN (1) | CN106056642B (zh) |
DE (1) | DE102015206630B4 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110123352A (zh) * | 2018-02-02 | 2019-08-16 | 西门子保健有限责任公司 | 成像参数的基于谱拓扑图的确定 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102015206630B4 (de) * | 2015-04-14 | 2022-05-05 | Siemens Healthcare Gmbh | Multispektrale CT-Bildgebung |
DE102017213479A1 (de) | 2017-08-03 | 2019-02-07 | Siemens Healthcare Gmbh | Computertomographische Aufnahme mit verschiedenen Energieschwellensätzen |
CN111278362B (zh) | 2017-09-22 | 2023-12-05 | 芝加哥大学 | 用于低剂量多谱x射线断层扫描的系统和方法 |
US11158095B2 (en) * | 2018-08-24 | 2021-10-26 | General Electric Company | System and method for reducing artifact bloom in a reconstructed object |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101308102A (zh) * | 2008-07-16 | 2008-11-19 | 中北大学 | 一种计算机断层扫描成像装置和方法 |
CN101566590A (zh) * | 2009-04-20 | 2009-10-28 | 西北工业大学 | 面阵探测器射线数字成像中的散射强度分布获取方法 |
US20120183118A1 (en) * | 2011-01-14 | 2012-07-19 | The Cleveland Clinic Foundation | Automated parameter selection for a tomographic imaging device |
CN102853791A (zh) * | 2012-02-23 | 2013-01-02 | 浙江大学 | 扫描超声波显微镜同时测量薄层材料厚度、声速、密度和衰减的方法 |
US20130022259A1 (en) * | 2011-07-22 | 2013-01-24 | Jae-Hyun Kwon | Apparatus for generating diagnosis image, medical imaging system, and method for processing image |
CN102902052A (zh) * | 2011-07-27 | 2013-01-30 | 莱卡微系统Cms有限责任公司 | 显微镜照明方法及显微镜 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6081577A (en) * | 1998-07-24 | 2000-06-27 | Wake Forest University | Method and system for creating task-dependent three-dimensional images |
US6453068B1 (en) * | 1999-09-17 | 2002-09-17 | Xerox Corporation | Luminance enhancement with overshoot reduction control based on chrominance information |
DE102006015454A1 (de) * | 2006-03-31 | 2007-10-18 | Siemens Ag | Verfahren und Vorrichtung zur automatischen Differenzierung von Nierensteintypen mittels Computertomographie |
DE102006015451A1 (de) * | 2006-03-31 | 2007-10-11 | Siemens Ag | Verfahren und Vorrichtung zur automatischen Differenzierung von Knochen oder anderen kalziumhaltigen Materialien und Kontrastmittel in Weichteilgewebe |
WO2008026150A1 (en) * | 2006-08-28 | 2008-03-06 | Koninklijke Philips Electronics N.V. | Method and apparatus for image enhancement |
AU2008330068B8 (en) * | 2007-11-27 | 2013-11-21 | Exxonmobil Upstream Research Company | Method for determining the properties of hydrocarbon reservoirs from geophysical data |
US8244013B2 (en) * | 2008-09-10 | 2012-08-14 | Siemens Medical Solutions Usa, Inc. | Medical image data processing and interventional instrument identification system |
US8885905B2 (en) * | 2009-12-03 | 2014-11-11 | Cedars-Sinai Medical Center | Method and system for plaque characterization |
WO2011067873A1 (ja) * | 2009-12-04 | 2011-06-09 | パナソニック株式会社 | 二酸化炭素還元方法、並びに、それに用いる二酸化炭素還元触媒および二酸化炭素還元装置 |
DE102010011911B4 (de) * | 2010-03-18 | 2019-09-26 | Siemens Healthcare Gmbh | Tomosyntheseverfahren mit einer iterativen Maximum-A-Posteriori-Rekonstruktion |
KR101689867B1 (ko) | 2010-09-15 | 2016-12-27 | 삼성전자주식회사 | 영상을 처리하는 방법, 이를 수행하는 영상처리장치 및 의료영상시스템 |
US8803490B2 (en) * | 2011-06-30 | 2014-08-12 | Broadcom Corporation | Current-mode active termination |
US9295443B2 (en) * | 2012-02-22 | 2016-03-29 | Koninklijke Philips N.V. | Method and system for reducing localized artifacts in imaging data |
JP6242631B2 (ja) * | 2012-08-30 | 2017-12-06 | 東芝メディカルシステムズ株式会社 | 医用画像処理装置及びx線コンピュータ断層撮影装置 |
US9401012B2 (en) * | 2012-10-18 | 2016-07-26 | Athentech Technologies Inc. | Method for correcting purple distortion in digital images and a computing device employing same |
US10282820B2 (en) * | 2013-02-21 | 2019-05-07 | Koninklijke Philips N.V. | Structure propagation restoration for spectral CT |
WO2016103839A1 (ja) * | 2014-12-22 | 2016-06-30 | オリンパス株式会社 | 超音波診断装置、超音波診断装置の作動方法および超音波診断装置の作動プログラム |
DE102015206630B4 (de) * | 2015-04-14 | 2022-05-05 | Siemens Healthcare Gmbh | Multispektrale CT-Bildgebung |
-
2015
- 2015-04-14 DE DE102015206630.4A patent/DE102015206630B4/de not_active Expired - Fee Related
-
2016
- 2016-04-07 US US15/092,784 patent/US10290122B2/en not_active Expired - Fee Related
- 2016-04-11 CN CN201610220577.2A patent/CN106056642B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101308102A (zh) * | 2008-07-16 | 2008-11-19 | 中北大学 | 一种计算机断层扫描成像装置和方法 |
CN101566590A (zh) * | 2009-04-20 | 2009-10-28 | 西北工业大学 | 面阵探测器射线数字成像中的散射强度分布获取方法 |
US20120183118A1 (en) * | 2011-01-14 | 2012-07-19 | The Cleveland Clinic Foundation | Automated parameter selection for a tomographic imaging device |
US20130022259A1 (en) * | 2011-07-22 | 2013-01-24 | Jae-Hyun Kwon | Apparatus for generating diagnosis image, medical imaging system, and method for processing image |
CN102902052A (zh) * | 2011-07-27 | 2013-01-30 | 莱卡微系统Cms有限责任公司 | 显微镜照明方法及显微镜 |
CN102853791A (zh) * | 2012-02-23 | 2013-01-02 | 浙江大学 | 扫描超声波显微镜同时测量薄层材料厚度、声速、密度和衰减的方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110123352A (zh) * | 2018-02-02 | 2019-08-16 | 西门子保健有限责任公司 | 成像参数的基于谱拓扑图的确定 |
CN110123352B (zh) * | 2018-02-02 | 2023-06-06 | 西门子保健有限责任公司 | 成像参数的基于谱拓扑图的确定 |
Also Published As
Publication number | Publication date |
---|---|
DE102015206630B4 (de) | 2022-05-05 |
CN106056642B (zh) | 2019-09-06 |
US20160307340A1 (en) | 2016-10-20 |
DE102015206630A1 (de) | 2016-10-20 |
US10290122B2 (en) | 2019-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5400326B2 (ja) | トモシンセシス画像を表示するための方法 | |
US9401019B2 (en) | Imaging tomosynthesis system, in particular mammography system | |
CN109035355B (zh) | 用于pet图像重建的系统和方法 | |
JP5036534B2 (ja) | 標的病変における変化の精密な測定評価のための医療用撮像システム | |
CN106056642A (zh) | 多谱式计算机断层成像 | |
US9140803B2 (en) | Acquisition protocol assessment apparatus | |
JP5442530B2 (ja) | 画像処理装置、画像表示装置およびプログラム並びにx線ct装置 | |
CN106232007B (zh) | X射线ct装置和处理装置 | |
US9993216B2 (en) | Creating a resultant image for a specifiable, virtual x-ray quanta energy distribution | |
KR101579111B1 (ko) | 영상 진단 방법, 장치 및 기록매체 | |
CN105377141B (zh) | X射线ct装置 | |
CN108511043B (zh) | 基于数值模拟的x-ct虚拟数据采集及图像重建方法及系统 | |
CN109381805A (zh) | 确定涉及多个组织区域的局部的组织功能的功能参数 | |
WO2013094483A1 (ja) | 医用画像診断装置及び医用画像診断装置を用いた位相決定方法 | |
CN105934199A (zh) | 用于处理体腔的医学图像的设备和方法 | |
CN107530042A (zh) | X射线ct装置和迭代修正参数决定方法 | |
EP3220826B1 (en) | Method and apparatus for processing medical image | |
CN108601570B (zh) | 断层摄影图像处理设备和方法以及与方法有关的记录介质 | |
US20170202532A1 (en) | Data processing method, data processing device, and x-ray ct apparatus | |
CN109419526A (zh) | 用于数字乳房断层合成中的运动评估和校正的方法和系统 | |
US10346957B2 (en) | Method for image improvement of image data from a dental image generation system | |
CN108430332A (zh) | 断层扫描设备及其用于重建断层扫描图像的方法 | |
US20230334732A1 (en) | Image rendering method for tomographic image data | |
US20140049542A1 (en) | Layer Display of Volume Data | |
JP7486614B2 (ja) | スペクトル暗視野イメージング |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190906 |