CN111201452B - 体积图像数据的材料选择性自适应混合 - Google Patents
体积图像数据的材料选择性自适应混合 Download PDFInfo
- Publication number
- CN111201452B CN111201452B CN201880065817.3A CN201880065817A CN111201452B CN 111201452 B CN111201452 B CN 111201452B CN 201880065817 A CN201880065817 A CN 201880065817A CN 111201452 B CN111201452 B CN 111201452B
- Authority
- CN
- China
- Prior art keywords
- image data
- voxels
- volume image
- spectral
- values
- 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.)
- Active
Links
- 239000000463 material Substances 0.000 title claims abstract description 69
- 238000002156 mixing Methods 0.000 title claims abstract description 24
- 230000003044 adaptive effect Effects 0.000 title description 3
- 230000003595 spectral effect Effects 0.000 claims abstract description 131
- 239000000203 mixture Substances 0.000 claims abstract description 16
- 238000003384 imaging method Methods 0.000 claims abstract description 14
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical group [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 claims description 22
- 229910052791 calcium Inorganic materials 0.000 claims description 22
- 239000011575 calcium Substances 0.000 claims description 22
- 238000000034 method Methods 0.000 claims description 18
- 125000001475 halogen functional group Chemical group 0.000 claims description 17
- 239000002872 contrast media Substances 0.000 claims description 8
- 210000004872 soft tissue Anatomy 0.000 claims description 7
- PNDPGZBMCMUPRI-UHFFFAOYSA-N iodine Chemical group II PNDPGZBMCMUPRI-UHFFFAOYSA-N 0.000 claims description 6
- 230000005855 radiation Effects 0.000 description 21
- 238000002591 computed tomography Methods 0.000 description 12
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 10
- 229910052740 iodine Inorganic materials 0.000 description 10
- 239000011630 iodine Substances 0.000 description 10
- 230000011218 segmentation Effects 0.000 description 10
- 238000001228 spectrum Methods 0.000 description 7
- 208000004434 Calcinosis Diseases 0.000 description 6
- 208000031481 Pathologic Constriction Diseases 0.000 description 6
- 230000002308 calcification Effects 0.000 description 6
- 210000001519 tissue Anatomy 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000036262 stenosis Effects 0.000 description 5
- 208000037804 stenosis Diseases 0.000 description 5
- 210000004204 blood vessel Anatomy 0.000 description 4
- 230000007704 transition Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- 206010060965 Arterial stenosis Diseases 0.000 description 2
- 230000002730 additional effect Effects 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- LEHOTFFKMJEONL-UHFFFAOYSA-N Uric Acid Chemical compound N1C(=O)NC(=O)C2=C1NC(=O)N2 LEHOTFFKMJEONL-UHFFFAOYSA-N 0.000 description 1
- TVWHNULVHGKJHS-UHFFFAOYSA-N Uric acid Natural products N1C(=O)NC(=O)C2NC(=O)NC21 TVWHNULVHGKJHS-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- -1 iodine, metals Chemical class 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 229940116269 uric acid Drugs 0.000 description 1
Classifications
-
- 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
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
- G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
-
- 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]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- 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
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20004—Adaptive image processing
- G06T2207/20012—Locally adaptive
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Medical Informatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Pulmonology (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
一种成像系统(500)包括被配置为产生投影数据的数据采集系统(515)和具有重建算法(518)和至少一个混合算法(524)的至少一个存储器设备。所述成像系统还包括重建器(516)，所述重建器被配置为利用所述重建算法来重建所述投影数据并且至少生成对应于第一基础材料成份的第一谱体积图像数据和对应于第二基础材料成份的第二谱体积图像数据，并且利用所述至少一个混合算法将所述第一谱体积图像数据与所述第二谱体积图像数据进行混合以产生混合体积图像数据。
Description
技术领域
以下总体上涉及成像并且更具体地涉及体积图像数据的材料选择性自适应混合并且具体应用于计算机断层摄影(CT)(包括谱和非谱CT)描述。
背景技术
计算机断层摄影(CT)扫描器通常包括安装在与一行或多行探测器相对的可旋转机架上的X射线管。X射线管围绕位于X射线管与一行或多行探测器之间的检查区域旋转并且发射穿过检查区域和被布置在检查区域中的受检者和/或对象的辐射。一行或多行探测器探测穿过检查区域的辐射并且生成指示检查区域和被布置在其中的受检者和/或对象的信号或投影数据。
投影数据被重建以借助于可以被用于生成一幅或多幅图像的计算机生成体积图像数据。(一幅或多幅)所得图像包括根据对应于相对放射密度的灰度级值表示的像素。这样的信息反映扫描的受检者和/或对象的衰减特性，并且通常示出诸如患者内的解剖结构、无生物对象内的物理结构等的结构。这些图像取决于X射线源和光子探测器的性质。
在检查的对象或者身体部分上给予的辐射还包括谱信息，因为由受检者和/或对象造成的辐射的吸收取决于贯穿其的光子的能量。这样的谱信息提供额外信息，诸如指示组织的元素或者材料组成(例如，原子数)和/或受检者和/或对象的材料的信息。然而，投影数据未反映谱特性，因为其表示的数据与在能量谱上积分的能量通量成比例。
被配置用于谱(多能量)成像的CT扫描器利用谱特性。例如，利用双能量系统，可以生成反映所成像的材料的固有性质(组织组成)的基础图像(例如，组织的每个分量的光电效应(PE)和康普顿散射(CS)行为等)。尽管这样的图像允许基于能量衰减特性的材料的鉴别，但是在一些实例中，存在图像质量改进与退化之间的折中。
例如，在通过钙化动脉斑块的动脉狭窄的非谱CT血管造影成像中，钙可以引起钙化表现为比其实际大小更大的“晕状”伪影。晕状伪影在较高的X射线能量处减少。然而，碘化的造影剂与周围组织之间的CT值差异(“对比”)被降低，并且该对比差异已经被用于分割动脉一用于进一步分析，例如，利用计算机断层摄影血流储备分数(CT-FFR)模拟，其已经被用于确定例如病变是否是流量限制的并且保证任何介入等。
图1-4图示了“晕状”伪影。图1示出了较低能量多平面重新格式化(MPR)图像中的钙化102和血管分割104，并且图2示出了较高能量MPR图像中的钙化202和血管分割204。图3A和3B是钙化102和202和分割104和204的子部分的放大视图。图3A中的放大(由于“晕状”)钙化102具有高度302，其大于图3B中的较高能量图像中的钙化102的高度304。因此，图3A中的血管306的管腔的大小小于图3B中的较高能量图像中的血管308的管腔的对应的大小，从而指示与图3B中的308相比较图3A中的更大的狭窄。
尽管图3B的较高能量图像中的狭窄更准确地反映实际狭窄几何结构，但是由于“晕状”伪影减少，碘化的造影剂与周围组织之间的减小的对比能够导致更不准确的血管分割。范例示出在图4A和4B中，其是分割104和204的子部分的放大视图。从图4A和4B，至少由于降低的对比，图4B的较高能量图像的血管壁比图4A中的血管壁更不规则。
发明内容
本文所描述的方法解决上文提到的问题和其他问题。
在一个方面中，一种成像系统包括被配置为产生投影数据的数据采集系统。所述成像系统还包括至少一个存储器设备。所述存储器至少包括重建算法和至少一个混合算法。所述成像系统还包括重建器，所述重建器被配置为利用所述重建算法来重建所述投影数据并且至少生成对应于第一基础材料成份的第一谱体积图像数据和对应于第二基础材料成份的第二谱体积图像数据，并且利用所述至少一个混合算法将所述第一谱体积图像数据与所述第二谱体积图像数据进行混合以产生混合体积图像数据。
在另一方面中，一种方法包括：执行对比增强扫描，其产生投影数据；并且重建所述投影数据以产生非谱体积图像数据。所述方法还包括重建所述投影数据以产生针对第一基础材料的第一谱体积图像数据，并且重建所述投影数据以产生针对第二基础材料的第二谱体积图像数据。所述方法还包括将所述非谱体积图像数据、所述第一谱体积图像数据和所述第二谱体积图像数据中的至少两个体积图像数据进行混合，从而产生混合体积图像数据。
在另一方面中，一种计算机可读介质编码有计算机可执行指令，所述指令当由计算机的处理器运行时使所述处理器：接收根据对比增强扫描生成的投影数据；重建所述投影数据以产生非谱体积图像数据；重建所述投影数据以产生针对第一基础材料的第一谱体积图像数据；重建所述投影数据以产生针对第二基础材料的第二谱体积图像数据；并且将所述非谱体积图像数据、所述第一谱体积图像数据和所述第二谱体积图像数据中的至少两个体积图像数据进行混合，从而产生混合体积图像数据。
附图说明
本发明可以采取各种部件和部件布置的形式，并且可以采取各种步骤和各步骤安排的形式。附图仅出于图示优选的实施例的目的而不应被解释为对本发明的限制。
图1示出了具有钙、晕状伪影和分割血管的较低能量对比增强图像的现有技术范例。
图2示出了具有钙、比图1更小的晕状伪影和分割血管的较高能量对比增强图像的现有技术范例。
图3A示出了图1的晕状伪影的放大视图。
图3B示出了图2的降低的晕状伪影的放大视图。
图4A示出了图1中的血管的正常部分的分割的放大视图。
图4B示出了图2中的血管的正常部分的分割的放大视图。
图5示意性地图示了具有重建器、重建算法和混合算法的成像系统。
图6图示了根据本文中的实施例的范例方法。
图7图示了根据本文中的实施例的范例方法。
图8示出了常规(120kVp)图像。
图9示出了虚拟70keV单能图像。
图10示出了虚拟120keV单能图像。
图11示出了混合70/120keV图像。
具体实施方式
图5示意性地图示了成像系统500，诸如计算机断层摄影(CT)扫描器。成像系统500包括总体固定机架502和旋转机架504。机架504的旋转部分由固定部分502支撑并且关于纵轴或z轴508围绕检查区域106旋转。受检者支撑体510(诸如卧榻)支撑检查区域中的对象或受检者。受检者支撑体510与执行成像流程协调可移动，从而相对于检查区域506引导受检者或对象，以用于加载、扫描和/或卸载受检者或对象。
诸如X射线管的辐射源512由旋转机架504可旋转地支撑。辐射源512随着旋转机架504旋转并且发射穿过检查区域506的X射线辐射。在所图示的实施例中，辐射源512是单个X射线管，其被配置为发射用于单个选定感兴趣峰值发射电压(kVP)的宽带辐射(即，该kVp处的能量谱)。在另一实例中，辐射源512被配置为在扫描期间在至少两个不同发射电压(例如，70keV、100keV等)之间切换。在又一实例中，辐射源512包括在旋转机架504上角偏移的两个或更多个X射线管，每个被配置为发射具有不同平均能量谱的辐射。US 8442184 B2描述了具有kVp切换和多个X射线管的系统，并且通过引用将其整体并入本文。
辐射谱敏感探测器阵列514跨检查区域506与辐射源512相对围成一角弧。探测器阵列514包括一排或多排探测器，其沿着z轴508方向相对于彼此布置并且探测穿过检查区域506的辐射。在所图示的实施例中，探测器阵列214包括能量分辨探测器，诸如多层闪烁体/光电传感器探测器(例如，US 7968853 B2，通过引用将其整体并入本文)和/或光子计数(直接转换)探测器(例如，WO2009072056A2，通过引用将其整体并入本文)。利用能量分辨探测器，辐射源512包括宽带kVp切换和/或多个X射线管辐射源512。在另一实例中，探测器阵列514包括非能量分辨探测器，并且辐射源512包括kVp切换和/或多个X射线管辐射源512。探测器阵列514生成指示不同能量的投影数据。
辐射源512和辐射谱敏感探测器阵列514可以被统称为数据采集的一部分(515)。
重建器516利用多个不同的重建算法518重建该相同组投影数据，包括(一个或多个)谱重建算法520和(一个或多个)非谱重建算法522。(一个或多个)非谱重建算法522产生常规宽带(非谱)体积图像数据，例如，通过组合谱投影数据和重建组合的体积图像数据。(一个或多个)谱重建算法520产生基础体积图像数据，例如，第一基础体积图像数据、第二基础体积图像数据、...、第N基础体积图像数据。例如，对于双能量，重建器516可以生成光电效应和康普顿散射体积图像数据集、单能/单色体积图像数据集(例如，70keV和100keV)、钙和碘体积图像数据集、骨和软组织体积图像数据集等。其他数据集包括有效Z(原子数)、k边缘等体积图像数据集。
如下面更详细地描述的，重建器516被配置为将体积图像数据混合在一起以基于(一个或多个)混合算法524产生混合体积图像数据。例如，重建器516可以选择性地组合两个或更多个虚拟单能和/或材料选择性体积图像数据集的不同区域。这可以减轻特定能量改进视觉特性同时使另一视觉特性退化的实例。重建器516包括至少一个处理器(例如，中央处理单元或者CPU、微处理器等)和计算机可读存储介质(其不包括瞬态介质)，诸如物理存储器、存储器设备和/或其他非瞬态存储器。计算机可读存储介质存储计算机可读指令和数据(例如，算法518和/或524)。至少一个处理器被配置为运行指令。至少一个处理器还可以运行由瞬态介质携带的指令，诸如信号、载波和/或其他瞬态介质。
操作者控制台526允许操作者控制系统500的操作。这包括选择图像采集协议(例如，多能量)、选择重建算法(例如，多能量)、选择混合算法、调用扫描、调用可视化软件应用、与运行可视化软件应用交互等。成像采集协议可以备选地自动选择，例如，经由人工智能和/或以其他方式。操作者控制台526包括(一个或多个)输出设备(诸如显示监测器、录像机等)和(一个或多个)输入设备(诸如鼠标、键盘等)。(一个或多个)重建算法518和/或(一个或多个)混合算法524可以存储在操作者控制台526、重建器516、位于远程位置的分离的计算机工作站的计算机可读存储介质中和/或其他存储介质中。
如上文简要讨论的，在一个实施例中，重建器516被配置为基于(一个或多个)混合算法524选择性地混合不同能量体积图像数据的不同区域以产生混合体积图像数据。例如，在一个实例中，重建器516采用重建算法，其产生非谱体积图像数据、针对第一感兴趣材料的第一谱体积图像数据和针对第二感兴趣材料的第二谱体积图像数据，全部根据来自相同扫描的相同投影数据的子集或全部。
对于体积图像数据，重建器516创建识别包含感兴趣材料的体素的图。例如，对于第一谱体积图像数据，产生第一图，其在体素基础上指代哪些体素主要包含第一材料。例如，主要具有第一材料的体素可以各自给定图中的一(1)的值，并且其他体素可以各自给定图中的零(0)的值。同样针对(一个或多个)其他感兴趣材料的体积图像数据完成。这些图然后被用于混合体积图像数据。
例如，对于主要包含第一材料的体素，如由对应的第一图所确定的，非谱体积图像数据的体素值向下加权并且第一谱体积图像数据的体素值向上加权，并且组合加权的体素值(例如，求和和标准化、平均等)。加权可以呈现从零(0)到一(1)的值。例如，向上加权和向下加权可以是1.0和0.0、0.9和0.1、0.8和0.2、...、0.5和0.5、0.1和0.9、和/或其他权重。
同样地，对于主要包含第二材料的体素，如由对应的第二图所确定的，非谱体积图像数据的体素值向下加权并且第二谱体积图像数据的体素值向上加权，并且组合加权的体素值。之间的体素可以被混合以用于平滑转换。这可以包括使用平滑函数，诸如S形、双曲正切、和/或其他平滑函数。通常，平滑函数将是单调递增(或者递减)函数。
对于不具有材料的体素，使用非谱体积图像数据的体素值，或者用于非谱体积图像数据的体素值的加权是一并且用于谱体积图像数据的体素值的加权是零。同样地，这些体素和针对第一或第二材料加权的体素之间的体素的值被混合以用于其之间的平滑转换。
在一个实例中，混合性质的该位置特异性调整允许取决于其材料成份而优化图像对比。
在选定协议是产生要针对通过钙化动脉斑块的动脉狭窄评价的体积图像数据的对比增强扫描的血管造影协议的情况下，第一材料和第二材料可以分别是钙和碘。在该实例中，重建器516重建钙基础、较高能量(例如，100keV)，和/或强调钙的其他谱体积图像数据，以及碘基础、较低能量(例如，70keV)，和/或强调碘的其他谱体积图像数据。在该实例中，第一图和第二图指示主要钙和主要碘体素。
对于主要包含钙的体素，非谱体积图像数据的体素值向下加权并且钙谱体积图像数据的体素值向上加权。对于主要包含碘的体素，非谱体积图像数据的体素值向下加权并且碘谱体积图像数据的体素值向上加权。不包含任一者的体素的值根据非谱体素图像数据确定。钙与碘之间、钙与非谱之间和碘与非谱之间的体素的值被混合以用于其之间的平滑转换。
在该范例中，混合体积图像数据减少关于较高能量图像数据的晕状伪影，同时维持碘化造影剂与具有非谱(或较低能量)图像数据的周围组织之间的CT值差异。即，在混合体积图像数据中，钙谱体积图像数据在更大的程度上(与其他体积图像数据有关)贡献于表示钙的体素，碘谱体积图像数据在更大的程度上(与其他体积图像数据有关)贡献于表示碘的体素，并且剩余体素表示其他(例如，非谱)体积图像数据或者之间的转换和/或具有调节值的体素。
该方法很适于自动和/或手动识别狭窄和分段血管的应用。一个这样的应用是血流储备分数(FFR)模拟，诸如CT-FFR。在一个实例中，相对于使用非谱体积图像数据、钙体积图像数据、碘体积图像数据或其他谱体积图像数据，根据所得的混合体积图像数据，狭窄的识别以及狭窄之前和之后的血管的壁的直径可以更准确地确定。在心脏应用中，这样的混合可以与用于一个或多个心脏相位的数据一起利用。
前述内容为钙和碘提供特定范例。其他感兴趣材料可以包括尿酸、金、骨、碘、金属和/或其他材料。通常，重建器516可以产生并且混合针对任何两个(或更多个)基础感兴趣材料的体积图像数据，以产生优化用于包含在其中的所有材料(或其预定子集)的图像质量的单个混合体积图像数据集，其中，个体一组体积图像数据可以强调一种材料对另一材料的损害。
在上文中，非谱体积图像数据和谱体积图像数据混合。在变型中，仅谱体积图像数据集混合。在另一变型中，仅一组谱体积图像数据与非谱体积图像数据混合。例如，在一个实例中，仅钙谱体积图像数据与非谱体积图像数据混合以减少晕状，同时维持来自非谱体积图像数据的碘和软组织对比差异。
图6图示了根据本文所描述的(一个或多个)实施例的范例方法。
应意识到，方法中的动作的排序是非限制性的。这样一来，在本文中预期其他排序。另外，可以省略一个或多个动作和/或可以包括一个或多个额外动作。
在602处，执行对比增强扫描。
在604处，重建非谱体积图像数据。
在606处，重建强调第一感兴趣材料的第一谱体积图像数据。
在608处，重建强调第二感兴趣材料的第二谱体积图像数据。
在610处，非谱体积图像数据、第一谱体积图像数据和/或第二谱体积图像数据被混合在一起以产生混合体积图像数据，如本文和/或以其他方式所描述的。
图7图示了根据本文所描述的(一个或多个)实施例的范例方法。
应意识到，方法中的动作的排序是非限制性的。这样一来，在本文中预期其他排序。另外，可以省略一个或多个动作和/或可以包括一个或多个额外动作。
在702处，执行对比增强扫描。
在704处，重建强调第一感兴趣材料的第一谱体积图像数据。
在706处，重建强调第二感兴趣材料的第二谱体积图像数据。
在708处，第一谱体积图像数据和第二谱体积图像数据被混合在一起以产生混合体积图像数据，如本文和/或以其他方式所描述的。
以上可以通过编码或嵌入在计算机可读存储介质上的计算机可读指令来实施，其当由(一个或多个)处理器运行时，使(一个或多个)处理器执行所描述的动作。额外地或者备选地，计算机可读指令中的至少一个由信号、载波或其他瞬态介质携带，所述介质不是计算机可读存储介质。
图8-11图示了本文所描述的方法的非限制性范例。图8示出了常规(非谱)120kVp图像800，图9示出了虚拟70keV单能图像900，图10示出了虚拟120keV单能图像1000，并且图11示出了根据图9的虚拟70keV单能图像和图10的虚拟120keV单能图像构建的混合图像1100。
如本文所讨论的，较高keV图像1000将具有相对于常规图像800的晕状伪影802和较低keV图像900的晕状伪影902的较少晕状伪影1002(例如，较小直径)。此外，较低keV图像900将具有相对于常规图像800的血管管腔分割804的更准确的血管管腔分割904和由于增加的造影剂/软组织对比的较高keV图像1000的血管管腔分割1004。
在该范例中，混合图像1100包括较高keV图像1000的部分和较低keV图像900的部分。特别地，包括表示钙的像素的区域1102、1104、1106和1108根据具有较低晕状伪影的较高keV图像1000构建，并且包括表示血管管腔而非钙的像素的区域1110、1112、1114和1118根据具有更准确的血管管腔分割的较低keV图像900构建。
应理解，该范例是非限制性的。混合的概念也可以应用于常规和较高keV图像以利用来自较高keV图像的钙的降低的晕状。通常，用于至少对比增强斑块研究的单能图像提供造影剂/软组织对比与晕状伪影之间的折中。通过选择性地混合不同单能图像的部分，混合图像可以利用增强造影剂/软组织对比和减少的晕状伪影两者构建，其比非谱图像或任何个体单能图像更准确。
尽管在附图和前面的描述中已经详细图示和描述了本发明，但是这些图示和描述应被视为说明性或示范性的而非限制性的；本发明不限于所公开的实施例。通过研究附图、公开内容以及权利要求书，本领域技术人员在实践所要求保护的本发明时能够理解和实现对所公开的实施例的其他变型。
在权利要求书中，“包括”一词不排除其他元件或步骤，并且词语“一”或“一个”不排除多个。单个元件或其他单元可以履行权利要求书中所记载的若干项目的功能。尽管在互不相同的从属权利要求中记载了特定措施，但是这并不指示不能有利地使用这些措施的组合。
计算机程序可以被存储/分布在合适的介质上，例如与其他硬件一起提供或作为其他硬件的部分提供的光学存储介质或固态介质，但计算机程序也可以以其他形式来分布，例如经由因特网或者其他有线或无线电信系统分布。权利要求书中的任何附图标记不应被解读为对范围的限制。
Claims (15)
1.一种成像系统(500)，包括：
数据采集系统(515)，其被配置为产生投影数据；
至少一个存储器设备，其具有重建算法(518)和至少一个混合算法(524)；以及
重建器(516)，其被配置为：
利用所述重建算法来重建所述投影数据并且至少生成对应于第一基础材料成份的第一谱体积图像数据和对应于第二基础材料成份的第二谱体积图像数据；
创建针对所述第一谱体积图像数据的第一图，所述第一图识别哪些体素主要包含所述第一基础材料成份并且哪些体素不主要包含所述第一基础材料成份；
创建针对所述第二谱体积图像数据的第二图，所述第二图识别哪些体素主要包含所述第二基础材料成份并且哪些体素不主要包含所述第二基础材料成份；并且
利用所述至少一个混合算法并且基于所述第一图和所述第二图将所述第一谱体积图像数据与所述第二谱体积图像数据进行混合以产生混合体积图像数据，所述混合体积图像数据拥有具有与由所述第一图识别的主要包含所述第一基础材料成份的第一谱体积图像数据的体素的值与来自所述第二谱体积图像数据的体素的值的加权组合相对应的值的体素，以及具有与由所述第二图识别的主要包含所述第二基础材料成份的第二谱体积图像数据的体素的值与来自所述第一谱体积图像数据的体素的值的加权组合相对应的值的体素。
2.根据权利要求1所述的系统，其中，所述重建器还被配置为重建所述投影数据并且生成非谱体积图像数据并且将所述非谱体积图像数据与所述第一谱体积图像数据或者所述第二谱体积图像数据中的至少一个谱体积图像数据进行混合以产生所述混合体积图像数据，所述混合体积图像数据拥有具有来自仅所述非谱体积图像数据的体素的值的体素，以及具有与由所述第一图识别的主要包含所述第一基础材料成份的第一谱体积图像数据的体素的值或者由所述第二图识别的主要包含所述第二基础材料成份的第二谱体积图像数据的体素的值与来自所述非谱体积图像数据的体素的值的加权组合相对应的值的体素。
3.根据权利要求1所述的系统，其中，所述重建器还被配置为重建所述投影数据并且生成非谱体积图像数据并且将所述非谱体积图像数据与所述第一谱体积图像数据和所述第二谱体积图像数据进行混合以产生所述混合体积图像数据，所述混合体积图像数据拥有：具有来自仅所述非谱体积图像数据的体素的值的体素；具有与由所述第一图识别的主要包含所述第一基础材料成份的第一谱体积图像数据的体素的值与来自所述非谱体积图像数据的体素的值的加权组合相对应的值的体素；以及具有与由所述第二图识别的主要包含所述第二基础材料成份的第二谱体积图像数据的体素的值与来自所述非谱体积图像数据的体素的值的加权组合相对应的值的体素。
4.根据权利要求1至3中的任一项所述的系统，其中，所述重建器被配置为通过在逐体素基础上组合体素值进行混合。
5.根据权利要求4所述的系统，其中，所述重建器被配置为通过对强调感兴趣材料成份的体积图像数据的体素的值进行向上加权并且对其他材料的体素的值进行向下加权来混合相同位置处的体素。
6.根据权利要求5所述的系统，其中，每个体素贡献于最终体素值。
7.根据权利要求5所述的系统，其中，仅一个体素贡献于最终体素值。
8.根据权利要求5至7中的任一项所述的系统，其中，所述重建器被配置为将其之间的体素的值进行平滑。
9.根据权利要求5至7中的任一项所述的系统，其中，所述重建器被配置为将邻近于其的体素的值进行平滑。
10.根据权利要求1至9中的任一项所述的系统，其中，所述第一基础材料成份是钙并且所述第二基础材料成份是碘，并且所述体积图像数据被混合以降低晕状伪影，同时维持造影剂与软组织之间的对比。
11.一种方法，包括：
存储重建算法(518)和至少一个混合算法(524)；
产生投影图像；
利用所述重建算法来重建所述投影数据并且至少生成对应于第一基础材料成份的第一谱体积图像数据和对应于第二基础材料成份的第二谱体积图像数据；
创建针对所述第一谱体积图像数据的第一图，所述第一图识别哪些体素主要包含所述第一基础材料成份并且哪些体素不主要包含所述第一基础材料成份；
创建针对所述第二谱体积图像数据的第二图，所述第二图识别哪些体素主要包含所述第二基础材料成份并且哪些体素不主要包含所述第二基础材料成份；并且
利用所述至少一个混合算法并且基于所述第一图和所述第二图将所述第一谱体积图像数据与所述第二谱体积图像数据进行混合以产生混合体积图像数据，所述混合体积图像数据拥有具有与由所述第一图识别的主要包含所述第一基础材料成份的第一谱体积图像数据的体素的值与来自所述第二谱体积图像数据的体素的值的加权组合相对应的值的体素，以及具有与由所述第二图识别的主要包含所述第二基础材料成份的第二谱体积图像数据的体素的值与来自所述第一谱体积图像数据的体素的值的加权组合相对应的值的体素。
12.根据权利要求11所述的方法，还包括：
重建所述投影数据并且生成非谱体积图像数据；并且
将所述非谱体积图像数据与所述第一谱体积图像数据或者所述第二谱体积图像数据中的至少一个谱体积图像数据进行混合以产生所述混合体积图像数据，所述混合体积图像数据拥有具有来自仅所述非谱体积图像数据的体素的值的体素，以及具有与由所述第一图识别的主要包含所述第一基础材料成份的第一谱体积图像数据的体素的值或者由所述第二图识别的主要包含所述第二基础材料成份的第二谱体积图像数据的体素的值与来自所述非谱体积图像数据的体素的值的加权组合相对应的值的体素。
13.根据权利要求11所述的方法，还包括：
重建所述投影数据并且生成非谱体积图像数据；并且
将所述非谱体积图像数据与所述第一谱体积图像数据和所述第二谱体积图像数据进行混合以产生所述混合体积图像数据，所述混合体积图像数据拥有：具有来自仅所述非谱体积图像数据的体素的值的体素；具有与由所述第一图识别的主要包含所述第一基础材料成份的第一谱体积图像数据的体素的值与来自所述非谱体积图像数据的体素的值的加权组合相对应的值的体素；以及具有与由所述第二图识别的主要包含所述第二基础材料成份的第二谱体积图像数据的体素的值与来自所述非谱体积图像数据的体素的值的加权组合相对应的值的体素。
14.根据权利要求11至13中任一项所述的方法，还包括：
通过向上加权强调感兴趣材料成份的体积图像数据的体素的值并且向下加权其他材料的体素的值来混合相同位置处的体素。
15.根据权利要求11至14中的任一项所述的方法，其中，所述第一基础材料成份是钙并且所述第二基础材料成份是碘，并且所述体积图像数据被混合以降低晕状伪影，同时维持造影剂与软组织之间的对比。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US201762569756P | 2017-10-09 | 2017-10-09 | |
| US62/569756 | 2017-10-09 | ||
| PCT/EP2018/076055 WO2019072554A1 (en) | 2017-10-09 | 2018-09-26 | SELECTIVE ADAPTIVE MIXING OF VOLUMETRIC IMAGE DATA MATERIALS |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111201452A CN111201452A (zh) | 2020-05-26 |
| CN111201452B true CN111201452B (zh) | 2023-09-12 |
Family
ID=63708373
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201880065817.3A Active CN111201452B (zh) | 2017-10-09 | 2018-09-26 | 体积图像数据的材料选择性自适应混合 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US11361415B2 (zh) |
| EP (1) | EP3695251A1 (zh) |
| JP (1) | JP7265541B2 (zh) |
| CN (1) | CN111201452B (zh) |
| WO (1) | WO2019072554A1 (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102021213438A1 (de) | 2021-11-29 | 2023-06-01 | Siemens Healthcare Gmbh | Verfahren zur Auswertung eines angiographischen Dualenergie-Computertomographiedatensatzes, Auswertungseinrichtung, Computerprogramm und elektronisch lesbare Datenträger |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101313332A (zh) * | 2005-09-06 | 2008-11-26 | 皇家飞利浦电子股份有限公司 | 具有多个能窗的计算机断层摄影中的数据处理和分析 |
| CN103649990A (zh) * | 2011-07-15 | 2014-03-19 | 皇家飞利浦有限公司 | 用于谱ct的图像处理 |
Family Cites Families (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7356115B2 (en) | 2002-12-04 | 2008-04-08 | Varian Medical Systems Technology, Inc. | Radiation scanning units including a movable platform |
| US7968853B2 (en) | 2005-04-26 | 2011-06-28 | Koninklijke Philips Electronics N.V. | Double decker detector for spectral CT |
| DE102006015452A1 (de) * | 2006-03-31 | 2007-10-11 | Siemens Ag | Verfahren und Vorrichtung zur Detektion von chemischen Anomalien und/oder Auffälligkeiten in Weichgewebe eines Objektbereiches |
| US8233683B2 (en) * | 2007-08-24 | 2012-07-31 | Siemens Aktiengesellschaft | Methods for non-linear image blending, adjustment and display |
| WO2009072056A2 (en) | 2007-12-04 | 2009-06-11 | Koninklijke Philips Electronics N. V. | Monolithically integrated crystalline direct-conversion semiconductor detector for detecting incident x-radiation at ultra-fine pitch and method for manufacturing such an x-ray semiconductor detector |
| JP5329103B2 (ja) * | 2008-02-01 | 2013-10-30 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像処理装置及びｘ線ｃｔ装置 |
| BRPI0910206A2 (pt) | 2008-06-30 | 2015-09-29 | Koninkl Philips Electronics Nv | sistema de formação de imagem de tomografia computadorizada e método |
| JP5666781B2 (ja) * | 2009-02-27 | 2015-02-12 | 株式会社東芝 | Ｘ線撮影装置及びｘ線画像処理方法 |
| JP6073616B2 (ja) * | 2011-09-28 | 2017-02-01 | 東芝メディカルシステムズ株式会社 | Ｘ線ｃｔ装置、画像処理装置及びプログラム |
| JP2013103023A (ja) * | 2011-11-15 | 2013-05-30 | Otoichi:Kk | 画像診断支援装置および画像診断支援方法 |
| US9271688B2 (en) * | 2012-03-28 | 2016-03-01 | General Electric Company | System and method for contrast agent estimation in X-ray imaging |
| DE102012214472B4 (de) * | 2012-08-14 | 2020-10-01 | Siemens Healthcare Gmbh | Verfahren zur Ermittlung von Dualenergie-Bilddatensätzen und eine Röntgeneinrichtung dazu |
| US9042512B2 (en) * | 2012-11-13 | 2015-05-26 | General Electric Company | Multi-sector computed tomography image acquisition |
| US9483851B2 (en) * | 2014-07-14 | 2016-11-01 | General Electric Company | Systems and methods for filtering for image generation |
| JP6313168B2 (ja) * | 2014-09-02 | 2018-04-18 | キヤノンメディカルシステムズ株式会社 | Ｘ線ｃｔ装置、画像処理装置及び画像処理プログラム |
| JP6738320B2 (ja) | 2014-09-05 | 2020-08-12 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | スペクトル画像データの視覚化 |
| US20170014069A1 (en) | 2015-07-14 | 2017-01-19 | Koninklijke Philips N. V. | Spectral imaging based fluid volume map |
-
2018
- 2018-09-26 CN CN201880065817.3A patent/CN111201452B/zh active Active
- 2018-09-26 WO PCT/EP2018/076055 patent/WO2019072554A1/en unknown
- 2018-09-26 JP JP2020517152A patent/JP7265541B2/ja active Active
- 2018-09-26 US US16/753,825 patent/US11361415B2/en active Active
- 2018-09-26 EP EP18779344.3A patent/EP3695251A1/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101313332A (zh) * | 2005-09-06 | 2008-11-26 | 皇家飞利浦电子股份有限公司 | 具有多个能窗的计算机断层摄影中的数据处理和分析 |
| CN103649990A (zh) * | 2011-07-15 | 2014-03-19 | 皇家飞利浦有限公司 | 用于谱ct的图像处理 |
Also Published As
| Publication number | Publication date |
|---|---|
| EP3695251A1 (en) | 2020-08-19 |
| CN111201452A (zh) | 2020-05-26 |
| JP7265541B2 (ja) | 2023-04-26 |
| US20210407056A1 (en) | 2021-12-30 |
| WO2019072554A1 (en) | 2019-04-18 |
| US11361415B2 (en) | 2022-06-14 |
| JP2020536602A (ja) | 2020-12-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10147168B2 (en) | Spectral CT | |
| Boas et al. | CT artifacts: causes and reduction techniques | |
| US9498179B1 (en) | Methods and systems for metal artifact reduction in spectral CT imaging | |
| US11813101B2 (en) | Image quality improved virtual non-contrast images generated by a spectral computed tomography (CT) scanner | |
| US20190142357A1 (en) | Systems and methods for improving image quality with three-dimensional scout | |
| CN111201452B (zh) | 体积图像数据的材料选择性自适应混合 | |
| CN110123352B (zh) | 成像参数的基于谱拓扑图的确定 | |
| US10383589B2 (en) | Direct monochromatic image generation for spectral computed tomography | |
| JP7258744B2 (ja) | スペクトルコンピュータ断層撮影フィンガープリンティング | |
| EP3968859B1 (en) | Protocol-dependent 2-d pre-scan projection image based on 3-d pre-scan volumetric image data | |
| US11875431B2 (en) | Method for providing automatic adaptive energy setting for ct virtual momochromatic imaging | |
| US20200342639A1 (en) | Single ct backprojector with one geometry calculation per voxel for multiple different types of projection data | |
| US20200294653A1 (en) | Spectral (mulit-energy) imaging visualization | |
| Tang | MARS Spectral CT: Image quality performance parameters using the Medipix3. 0 detector |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |