CN109688931A - 谱计算机断层扫描指纹识别 - Google Patents

谱计算机断层扫描指纹识别 Download PDF

Info

Publication number
CN109688931A
CN109688931A CN201780043267.0A CN201780043267A CN109688931A CN 109688931 A CN109688931 A CN 109688931A CN 201780043267 A CN201780043267 A CN 201780043267A CN 109688931 A CN109688931 A CN 109688931A
Authority
CN
China
Prior art keywords
cluster
histogram
dimensions histogram
projection
dimensions
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201780043267.0A
Other languages
English (en)
Other versions
CN109688931B (zh
Inventor
S·戈特曼
L·戈申
R·E·伦金斯基
M·A·刘易斯
T·C·泽斯比
X·段
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.)
Koninklijke Philips NV
University of Texas System
Original Assignee
Koninklijke Philips Electronics NV
University of Texas System
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 Koninklijke Philips Electronics NV, University of Texas System filed Critical Koninklijke Philips Electronics NV
Publication of CN109688931A publication Critical patent/CN109688931A/zh
Application granted granted Critical
Publication of CN109688931B publication Critical patent/CN109688931B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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/42Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4208Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
    • A61B6/4241Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
    • 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/46Arrangements for interfacing with the operator or the patient
    • A61B6/467Arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B6/469Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • 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/46Arrangements for interfacing with the operator or the patient
    • A61B6/461Displaying means of special interest
    • A61B6/463Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • 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/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Optics & Photonics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Human Computer Interaction (AREA)
  • Pulmonology (AREA)
  • Quality & Reliability (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

一种方法包括生成对象的谱投影数据,所述谱投影数据至少包括对应于第一能量范围的第一谱投影数据和对应于第二能量范围的第二谱投影数据,其中,所述第一能量范围和所述第二能量范围是不同的。所述方法还包括根据谱投影数据的组合来构建图像,并针对5能量范围并且根据谱投影数据来构建基础图像的集合。所述方法还包括从所述基础图像的集合构建多维直方图,其中,所述多维直方图包括至少两个轴,第一个轴对应于第一基础分量,并且第二个轴对应于第二基础分量,并且所述多维直方图包括聚类的集合,针对在谱投影数据中表示的每种材料10包括一个聚类。所述方法还包括在视觉上同时呈现所述图像和所述多维直方图。

Description

谱计算机断层扫描指纹识别
技术领域
以下总体涉及谱成像,更具体地涉及谱计算机断层摄影(CT)指纹识别。
背景技术
CT扫描器通常包括安装在可旋转机架上的X射线管,与一个或多个集成探测器相对。X射线管围绕定位于所述X射线管与一个或多个探测器之间的检查区域旋转,并且发射穿过检查区域和设置在检查区域中的对象和/或目标的辐射。一个或多个探测器探测穿过检查区域的辐射并生成指示检查区域和设置在其中的对象和/或目标的信号或投影数据。
重建投影数据以生成体积图像数据,其可用于生成一幅或多幅图像。得到的(一幅或多幅)图像包括以对应于相对辐射密度的灰度值表示的像素。这样的信息反映了被扫描的对象和/或目标的衰减特性,并且通常示出了诸如患者体内的解剖结构、无生命对象内的物理结构等等的结构。这些图像高度依赖于X射线源和光子探测器的特性。
探测到的辐射还包括谱信息,因为对象和/或目标对辐射的吸收取决于穿过其中的光子的能量。这样的谱信息提供额外的信息,例如指示组织的元素或材料成分(例如,原子序数)和/或对象和/或目标的材料的信息。然而,对于传统CT,投影数据不反映谱特性,因为它所代表的数据与在能量谱上积分的能量注量成比例。
针对谱(多能量)成像配置的CT扫描器利用谱特性。例如,对于双能量系统,可以生成反映被成像材料的固有特性的基础图像(例如,组织的每种成分的光电效应(PE)和康普顿散射(CS)行为)。
尽管这样的图像允许基于能量衰减特性来区分材料,但是它们不能很好地识别特定类型材料的存在或不存在,例如识别特定类型组织的正常和/或患病组织。
发明内容
本文中描述的各方面解决以上提及的问题和其他问题。
在一个方面中,一种方法包括生成对象的谱投影数据,至少包括对应于第一能量范围的第一谱投影数据和对应于第二能量范围的第二谱投影数据,其中所述第一能量范围和所述第二能量范围是不同的。所述方法还包括根据谱投影数据的组合来构建图像,并针对能量范围并且根据谱投影数据来构建基础图像的集合。所述方法还包括从所述基础图像的集合构建多维直方图,其中,所述多维直方图包括至少两个轴,第一个轴对应于第一基础分量,并且第二个轴对应于第二基础分量,并且所述多维直方图包括聚类的集合,针对在谱投影数据中表示的每种材料的包括一个聚类。所述方法还包括在视觉上同时呈现所述图像和所述多维直方图。
在另一方面中,一种系统包括存储有指令的存储器,包括谱CT指纹模块和被配置为至少执行与所述谱CT指纹模块相对应的指令的处理器。所述指令使所述处理器:生成对象的谱投影数据,至少包括对应于第一能量范围的第一谱投影数据和对应于第二能量范围的第二谱投影数据,其中,所述第一能量范围和所述第二能量范围是不同的;根据谱投影数据的组合来构建图像;针对所述能量范围并且根据所述谱投影数据来构建基础图像的集合;根据所述基础图像的集合构建多维直方图,其中,所述多维直方图包括至少两个轴,第一个轴对应于第一基础分量,并且第二个轴对应于第二基础分量,并且所述多维直方图包括聚类的集合,针对在谱投影数据中表示的每种材料所述聚类的集合包括一个聚类,并且在视觉上同时呈现所述图像和所述多维直方图。
在另一方面中,一种计算机可读介质被编码有计算机可执行指令,所述计算机可读指令在由计算机的处理器运行时,使所述处理器:生成对象的谱投影数据,至少包括对应于第一能量范围的第一谱投影数据和对应于第二能量范围的第二谱投影数据,其中,所述第一能量范围和所述第二能量范围是不同的;根据谱投影数据的组合来构建图像;针对所述能量范围并且根据所述谱投影数据来构建基础图像的集合;根据所述基础图像的集合构建多维直方图,其中,所述多维直方图包括至少两个轴,第一个轴对应于第一基础分量,并且第二个轴对应于第二基础分量,并且所述多维直方图包括聚类的集合,针对在谱投影数据中表示的每种材料所述聚类的集合包括一个聚类,并且在视觉上同时呈现所述图像和所述多维直方图。
在另一方面中,一种方法包括:利用成像系统来生成表示针对预定类型组织的已知正常组织的第一谱投影数据,至少包括对应于第一能量范围的第一谱投影数据和对应于第二能量范围的第二谱投影数据,其中,所述第一能量范围和所述第二能量范围是不同的。所述方法还包括利用成像系统来生成表示预定类型组织的异常组织的第二谱投影数据,至少包括对应于所述第一能量范围的第一谱投影数据和对应于所述第二能量范围的第二谱投影数据。所述方法还包括针对所述能量范围并且根据所述第一谱投影数据来构建基础图像的第一集合,并且针对所述能量范围并且根据所述第二谱投影数据来构建基础图像的第二集合。所述方法还包括根据基础图像的第一集合构建第一多维直方图,其中,所述第一多维直方图包括至少两个轴,第一个轴对应于第一基础分量,并且第二个轴对应于第二基础分量,根据基础图像的第二集合来构建第二多维直方图,其中,所述第二多维直方图包括至少两个轴,第一个轴对应于所述第一基础分量并且第二个轴对应于所述第二基础分量。所述方法还包括从所述第二多维直方图中识别不存在于所述第一多维直方图中的聚类,将识别出的聚类标记为表示已知的异常组织,并且将所述第一多维直方图保存为针对所述正常组织的谱成像指纹并且将所述第二多维直方图标记为针对异常组织的谱成像指纹。
附图说明
本发明可以采取各种部件和各部件的布置以及各种步骤和各步骤的安排的形式。附图仅出于图示优选的实施例的目的并且不应被解释为对本发明的限制。
图1示意性地图示了包括具有谱CT指纹模块的控制台的成像系统。
图2示意性地图示了成像系统的重建器的示例。
图3图形地示出了针对双能量投影数据的理论材料衰减分解(MAD)图。
图4描绘了2-D图像和相应的MAD图的视觉显示。
图5描绘了图4的2-D图像,其中识别出了它里面的ROI,并且更新的MAD图突出显示对应于ROI中的材料的聚类。
图6描绘了图4的MAD图,其中在它里面识别出了针对第一区域的ROI,并且更新的2-D图像突出显示与ROI中表示的材料相对应的像素或体素。
图7描绘了图4的MAD图,其中在它里面识别出了针对第二不同区域的ROI,并且更新的2-D图像突出显示与ROI中表示的材料相对应的像素或体素。
图8图示了用于创建与多能量投影数据相关的谱CT指纹的示例方法。
图9图示了用于创建与多能量投影数据相关的谱CT指纹的另一示例方法。
图10图示了用于创建与多能量投影数据相关的谱CT指纹的又一示例方法。
图11图示了用于结合多能量投影数据来识别异常的存在和/或类型的示例方法。
图12图示了用于结合多能量投影数据来识别异常的不存在的示例方法。
图13图示了用于结合多能量投影数据来识别异常的不存在的另一示例方法。
图14示出了用于结合多能量投影数据来可视化图像中感兴趣组织的示例方法。
图15示出了用于结合多能量投影数据来识别图像中感兴趣组织的示例方法。
图16示意性地示出了图1的变型,其中,谱CT指纹模块是远离成像系统的系统的一部分。
具体实施方式
图1示意性地示出了成像系统100,例如被配置用于谱(多能量)成像的计算机断层摄影(CT)扫描器。成像系统100包括大致固定的机架102和旋转机架104。旋转机架104由固定机架102可旋转地支撑,并围绕检查区域106绕纵向或z轴108旋转。对象支撑体110,例如卧榻,在检查区域中支撑目标或对象。对象支撑体110可以与执行成像流程协调地移动,以便相对于检查区域106引导对象或目标,以加载、扫描和/或卸载对象或目标。
辐射源112(例如X射线管)由旋转机架104可旋转地支撑。辐射源112与旋转机架104一起旋转,并且发射穿过检查区域106的多色辐射。在所示的实施例中,辐射源112是(标准)单X射线管,其配置为发射针对感兴趣的单个选定峰值发射电压(kVp)(即,所述kVp处的能谱)的辐射。在另一实例中,辐射源112被配置为在扫描期间在至少两个不同的发射电压(例如,80kVp,140kVp等)之间切换。在又一实例中,辐射源112包括在旋转机架104上在角度上偏置的两个或更多个X射线管,每个X射线管被配置为发射具有不同平均能量谱的辐射。皇家飞利浦有限公司的US8442184B2描述了具有kVp切换和多个X射线管的系统。
辐射敏感探测器阵列114在检查区域106对面与辐射源112相对,对向一角度弧。探测器阵列114包括沿z轴108方向相对于彼此布置的一行或多行探测器,并且探测穿过检查区域106的辐射。在所示实施例中,探测器阵列114包括能量分辨探测器,例如多层闪烁体/光传感器探测器。US8442184B2和US 7968853B2,两者均归属于皇家飞利浦有限公司,描述了一种多层探测器。利用能量分辨探测器,辐射源112包括标准,kVp切换和/或多个X射线管辐射源112。在另一实例中,探测器阵列114包括非能量分辨探测器,并且辐射源112包括kVp切换和/或多个X射线管辐射源112。探测器阵列114生成并输出指示N个不同能量中的每一个的信号(投影数据),其中N是等于或大于2的整数。
重建器116重建由探测器阵列114输出的信号(投影数据)。图2示意性地图示了重建器116的范例。在该示例中,重建器116接收针对N个不同能量中的每个通量的投影数据(N个能量投影数据)。谱分解器202将N个不同的能量分解成不同的基础投影数据,并且投影数据组合器204将N个不同的能量组合成不同的基础投影数据。例如,对于双能量(N=2),例如第一(较高)和第二(较低)能量配置,谱分解器202将第一能量投影数据和第二能量投影数据分解为光电效应基础投影数据和康普顿散射基础投影数据。对于同样的情况,投影数据组合器204组合第一能量投影数据和第二能量投影数据以产生等同于非谱图像数据的组合数据。
重建处理器206基于来自算法存储器208的重建算法来重建基础投影数据和/或组合投影数据。这可以包括重建一个或多个基础图像(例如,第一基础图像,第二基础图像,......,和第N基础图像)。例如,对于双能量示例,重建处理器206可以生成光电效应图像、康普顿散射图像、和/或组合光电效应/康普顿散射图像。这可以额外地包括从组合投影图像重建传统(非谱)图像。可以从该数据生成的其他图像包括但不限于低能量图像和高能量图像,单能量或单色图像(例如,140kVp),虚拟非对比(或未增强的)图像,有效Z(原子序数)图像,碘图等。重建器116输出重建的图像。
返回图1,操作者控制台118(计算系统)允许操作者控制系统100的操作。这包括选择成像采集协议(例如,多能量),选择重建算法(例如,多能量),调用扫描,调用可视化软件应用程序,与执行的可视化软件应用程序交互等。操作者控制台118包括输入/输出(I/O)120,其便于与至少一个输出设备122(例如显示监视器,胶片放映机(filmer)等)、(一个或多个)输入设备124(例如鼠标,键盘等)进行通信。操作者控制台118还包括至少一个处理器125(例如,中央处理单元或CPU,微处理器,控制器等)和计算机可读存储介质126(其排除瞬态介质),例如物理存储器和/或其他非瞬态存储器。计算机可读存储介质128存储计算机可读指令132和数据130。处理器125执行指令132。
所示的计算机可读指令130至少包括谱CT指纹模块134。在一个实例中,谱CT指纹模块134包括用于创建材料衰减分解(MAD)图(或直方图)的指令,其用作谱CT指纹。额外地,或者替代地,谱CT指纹模块128包括用于采用这样的MAD图的指令。如本文中所使用的,谱CT指纹将指示特定类型的组织/材料的谱CT数据映射到该类型的组织/材料。这样,谱CT指纹可用于构建视觉上强调特定类型的感兴趣组织/材料的图像。谱CT指纹还可以用于确定图像中特定类型的组织/材料的存在或不存在。谱CT指纹还可用于识别图像中的未知组织/材料。谱CT指纹是系统和能量不变的,并且可以结合任何CT成像系统使用。因此,它还可以用于确认CT成像系统被正确地校准。
图3图示了针对双能量的理论MAD图300。第一(或y-)轴302表示第一基础,第二(或x-)轴表示第二基础。在所示实施例中,第一轴302表示光电效应衰减系数,并且第二轴表示康普顿散射衰减系数。在一个变型中,系数值是无量纲的或具有其他单位,例如亨氏单位(HU)。可以从由重建器116生成的基础图像获得光电效应和康普顿散射衰减系数。MAD图300中表示的材料是出于解释目的并且不是非限制性的。其他材料,包括更少的或更多的,相似的或不同的,等等,可以在其他MAD图中表示。示例包括两个单能量图像,低能量和高能量图像对,等等。MAD图中的每个基础表示另一个维度。
在MAD图300中,原点306表示空气。第一聚类308提供了针对脂肪组织(脂肪)的谱CT指纹。第二聚类310提供了针对水的谱CT指纹。第三聚类312提供了针对软组织的谱CT指纹。第四聚类314提供了针对肺组织的谱CT指纹。第四聚类314在一定范围内延伸,例如,至少由于空气和肺组织的部分体积效应。第五聚类316提供了针对肺碘的谱CT指纹。第五聚类316在一定范围内延伸,例如,作为碘浓度水平的函数,在第一轴302上更高的水平由高的衰减系数值表示。第六聚类318提供针对钙的谱CT指纹。第六聚类318在一定范围内延伸,例如,至少由于骨和软组织的部分体积效应。纯钙在该范围内在320处表示。
回到图1,可以针对任何数量的组织/材料生成MAD图。这包括针对特定类型的已知正常组织和特定类型的已知异常组织生成MAD图。例如,可以根据已知具有健康肝脏(或其他组织)的患者群体的图像数据生成MAD图,并且可以根据已知具有特定肝脏(或其他组织)肿瘤或病变的患者群体的图像数据生成另一个MAD图。MAD图可以是针对整个身体、特定身体部位和/或特定器官的。在另一个示例中,可以根据通过扫描由特定材料(例如,骨)构造的“体模”和/或特定材料(例如,造影剂)的已知分布的图像数据生成MAD图。MAD图可以本地存储(例如,在存储器128中)和/或远程存储(例如,在健康护理实体内,在服务器上,在“云”中等)。
MAD图可以经由输出设备124的显示监视器,与图像同时地、单独地和/或与图像交替地可视地显示。图4-7图示了本文中描述的MAD图的示例可视化和使用。所有图4-7同时示出了输出设备124的显示监视器的图形用户界面(GUI)中的图像和MAD图。
图4示出了2-D图像402和相应的MAD图404。所述MAD图404是基于2-D图像402的全部像素或体素来创建的。类似于图3中的MAD图300,MAD图404是光电效应对康普顿散射衰减系数2-D直方图。在该示例中,使用预定水平和窗口衰减值设置来显示MAD图404,以在视觉上强调特定感兴趣组织的细节。聚类406表示类似空气的材料,聚类408表示组织/脂肪样材料。示出了针对区域412的放大视图410。
在图5中,在2-D图像402上识别感兴趣区域(ROI)502、504和506。ROI 502包括肺组织,ROI 504包括皮下脂肪组织,ROI 506包括肝脏组织。这些ROI可以被自动识别(通过用户批准、拒绝或修改),通过系统和用户操作半自动识别,或者由用户手动识别,例如,经由软件和诸如鼠标的输入设备124。所示的形状不是限制性的。本文中也预期其他形状,包括像正方形、矩形、圆形等的预定形状。
响应于对2-D图像402上的ROI 502、504和506的识别,MAP图404被更新为在视觉上强调对应于ROI 502、504和506的聚类406和408的部分。强调可以是颜色、灰度、轮廓、突出显示等。在图5中,通过放大视图410通过灰度显示强调,其中,第一突出显示508对应于ROI502(肺组织),第二突出显示510对应于ROI 504(脂肪组织)并且第三突出显示512对应于ROI 506(肝组织)。由于空气和肺组织的部分体积效应,ROI 508被模糊。
在图6中,在MAD图404上识别ROI 602、604和606。类似地,ROI602、604和606可以是半自动或手动识别的。在一个实例中,ROI 602、604和606中的一个或多个可以对应于要在图像402中被可视化的已知组织。在另一个示例中,ROI与未知组织相对应。针对空气状材料的聚类406(图4)的放大视图608是针对区域610的。
通过将MAD图404与感兴趣组织的MAD图比较,可以在MAD图404中识别区域610。由于MAD图提供组织的指纹,因此MAD图404中的与针对组织的MAD图中的区域相匹配的区域识别MAD图404中的针对该组织的区域。然后可以将ROI放置在匹配的区域上。在没有匹配的MAD图或者用户想要看到什么解剖结构对应于MAD图的特定区域的情况下,可以在该区域上放置ROI。在图6中,ROI 602、604和606对应于空气、床泡沫和肺。
响应于在MAD图404上识别ROI 602、604和606,更新2-D图像402以在视觉上强调2-D图像402的对应于ROI 602、604和606的部分。强调可以是颜色、灰阶、轮廓、突出显示等。在图6中,通过灰度显示强调,其中,第一突出显示612对应于ROI 602(空气),第二突出显示614对应于ROI 604(床泡沫)并且第三突出显示616对应于ROI 606(肺)。
如果用户希望在MAD图404中看到特定于与ROI相对应的材料的图像,则可以为该材料创建合适的材料分解,并且重建器116可以使用该材料分解来处理所接收的N能量投影数据。例如,MAD图可以用于通过识别针对MAD图中的硅的聚类来来识别对应于硅乳房植入物的像素或体素。可以识别所述聚类,并且然后可以生成硅基础图像,可以突出显示传统非谱图像中的硅等。
在图7中,在MAD图404上识别ROI 702、704和706。类似地,ROI702、704和706可以是半自动或手动识别的。在一个实例中,ROI 702、704和706中的一个或多个可以对应于要在图像402中被可视化的已知组织。针对组织/脂肪样材料的聚类408(图4)的放大视图708是针对区域710的。通过将MAD图404与感兴趣组织的MAD图比较,可以在MAD图404中识别区域710。
在图7中,ROI 702、704和706对应于脂肪组织、肝组织和肋骨软骨。响应于在MAD图404上识别ROI 702、704和706,用对应于ROI 702(脂肪组织)的第一突出显示712,对应于ROI704(肝组织)的第二突出显示714和对应于ROI 706(肋骨软骨)第三突出显示716来更新2-D图像402。同样,强调可以通过颜色、灰阶、轮廓、突出显示等来进行。
如图4-7中所示,并且如本文中所描述,MAD图中像素的位置取决于人体内组织和元素的已知特性(即,光电效应和康普顿散射的衰减系数)。这些位置占据整个MAD图的一小部分。在这些范围之外的任何像素都是由人体内的异物(即金属器具或硅树脂植入物)中的异物或材料引起的,或者是来自由于信号采集或图像重建中的误差而引起的伪影。这些伪影的精确图案和/或位置可以提供关于它们的起源和性质的线索。因此,MAD图可以提供对扫描器性能进行评估和改进两者的工具。
图8-14呈现了与本文描述的MAD图有关的各种方法。要理解,方法中的动作的顺序不是限制性的。这样,本文中预期其他顺序。另外,可以省略一个或多个动作和/或可以包括一个或多个额外的动作。
图8图示了用于创建与多能量投影数据相关的谱CT指纹的示例方法。
在802处,执行包括感兴趣材料(MOI)的对象或目标的谱扫描,创建多能量投影数据,投影数据的每个集合对应于不同的能量。
在804,利用多能量投影数据生成多个基础图像,每个基础图像对应于不同的能量。
在806处,从多个基础图像创建MAD图。
在808处,在MAD图中识别MAD图中的对应于感兴趣的材料的聚类。这可以使用来自文献和/或其他方面的已知光电效应和康普顿散射衰减系数,来自具有和没有感兴趣材料的先前扫描的数据,和/或其他来实现。
在810处,标记聚类以将聚类与感兴趣材料相关联。
在812处,将MAD图保存到例如计算机存储器以供使用。
图9图示了用于创建与多能量投影数据相关的谱CT指纹的另一示例方法。
在902处,执行除了图8中的感兴趣材料之外的包括感兴趣材料的对象或目标的谱扫描,从而创建多能量投影数据。
在904处,如本文所述地利用多能量投影数据生成多个基础图像。
在906处,从多个基础图像创建MAD图。
在908处,通过将MAD图与图8中存储的MAD图比较来在MAD图中对MAD图中的对应于不同感兴趣材料的聚类进行定位以定位MAD图中不在存储的MAD中的聚类。
在910处,标记聚类以将聚类与不同的感兴趣材料相关联。
在912处,保存MAD图。
图10图示了用于创建与多能量投影数据相关的谱CT指纹的另一示例方法。
在1002处,执行对包括没有感兴趣的材料的第一对象或目标的第一谱扫描,从而创建第一投影数据。
在1004处,执行对包括感兴趣的材料的第二对象或目标的第二谱扫描,从而创建第二投影数据。
动作1002和1004的投影数据可以是新采集的投影数据和/或从先前检查中挖掘的投影数据。此外,动作1004可以在动作1002之前。
在1006处,利用所述第一投影数据生成第一多个基础图像。
在1008处,利用所述第二投影数据生成第二多个基础图像。
在1010处,根据所述第一多个基础图像创建第一MAD图。
在1012处,根据所述第二多个基础图像创建第二MAD图。
在1014处,通过比较所述第一MAD图和所述第二MAD图来定位所述第一MAD图中的在所述第二MAD图中不存在的聚类,以定位所述第一MAD图中对应于所述感兴趣材料的聚类。
在1016处,标记聚类以将聚类与感兴趣材料相关联。
在1018处,保存所述第一MAD图。
图11图示了用于结合多能量投影数据来识别异常的存在和/或类型的示例方法。
在1102处,执行对象的谱扫描,创建投影数据。
在1104处,利用所述投影数据生成多个基础图像。
在1106处,根据所述多个基础图像创建MAD图。
在1108处,将所述MAD图与针对正常组织的MAD图进行比较。
在1110处,响应于MAD图包括在针对正常组织的MAD图中不存在的聚类而确定对象具有异常。
在1112处,将所述MAD图与针对异常的MAD图进行比较。
在1114处,响应于所述MAD图包括针对异常的与MAD图中表示异常的的聚类相匹配的聚类,确定所述对象具有异常。
在一个变型中,省略了动作1108和1110,并且仅将MAD图与针对异常的MAD图进行比较。在另一个变型中,省略了动作1112和1114,并且仅将MAD图与针对正常组织的MAD图进行比较。
图12图示了用于结合多能量投影数据来识别一种异常的不存在的示例方法。
在1202处,执行对象的谱扫描,创建投影数据。
在1204处,利用所述投影数据生成多个基础图像。
在1206处,根据所述多个基础图像创建MAD图。
在1208处,将所述MAD图与针对正常组织的MAD图进行比较。
在1210处,响应于MAD图与针对正常组织的的MAD图相匹配而确定对象不具有异常。
图13图示了用于结合多能量投影数据来识别一种异常的不存在的示例方法。
在1302处,执行对象的谱扫描,创建投影数据。
在1304处,利用所述投影数据生成多个基础图像。
在1306处,根据所述多个基础图像创建MAD图。
在1308处,将所述MAD图与针对异常的MAD图进行比较。
在1310处,响应于所述MAD图包括与针对异常的MAD图中表示异常的聚类相匹配的聚类,确定所述对象不具有异常。
图14示出了用于结合多能量投影数据来可视化图像中感兴趣组织的示例方法。
在1402处,执行对象的谱扫描,创建投影数据。
在1404处,利用所述投影数据生成多个基础图像。
在1406处,根据所述多个基础图像创建MAD图。
在1408处,识别MAD图中对应于感兴趣组织的聚类。这可以通过将MAD图与针对感兴趣组织和/或其他组织的MAD图进行比较来实现。
在1410处,针对聚类创建ROI。
在1412处,利用所述投影数据来生成图像。
在1414处,针对由聚类表示的材料来在图像中突出显示对应于感兴趣组织的像素或体素。
图15图示了用于结合多能量投影数据来识别图像中感兴趣组织的示例方法。
在1502处,执行对象的谱扫描,创建投影数据。
在1504处,利用所述投影数据生成多个基础图像。
在1506处,利用所述投影数据来生成图像。
在1508处,在图像中创建ROI。
在1510处,针对所述ROI创建MAD图。
在1512处,识别所述MAD图中的对应于所述ROI的聚类。
在1514处,通过所述聚类来识别所述ROI中的类型的材料。这可以通过将MAD图与针对不同类型的材料的MAD图进行比较和/或以其他方式来实现。
在1516处,利用识别出的材料来创建材料分解。
在1518处,使用所述材料分解来生成材料的谱图像。
本文中的动作可以通过被编码或录入在计算机可读存储介质上的计算机可读指令来实施,所述计算机可读指令当由(一个或多个)计算机处理器运行时,令所述(一个或多个)处理器执行所描述的动作。额外地或者替代地,所述计算机可读指令中的至少一个由信号、载波或其他瞬态介质所承载。
图16示意性地示出了图1的变型,其中与谱CT指纹模块132有关的操作者控制台118的部件是计算系统1602的一部分,计算系统1602是与成像系统100分开的装置。操作者控制台118包括诸如监视器或显示器的人类可读输出设备和诸如键盘和鼠标的输入设备,并且驻留在控制台118上的软件允许操作者与扫描器100交互。指令130包括用于重建体积图像数据的重建器。
计算系统160可以从成像系统100(控制台118和/或重建器116)、数据储存库1604、另一成像系统和/或其他设备接收投影。合适的数据存储库1604的范例是图片存档及通信(PACS)、放射信息系统(RIS)、医院信息系统(HIS)、电子医学数据(EMR)、数据库、服务器、成像系统、和/或其他数据存储库。
本文中讨论的MAD图也可用于质量保证。例如,由于针对特定材料的MAD图独立于X射线源和探测器阵列,因此可以通过扫描包括特定材料的体模来测试多个成像系统,针对每个成像系统生成针对该材料的MAD图,并将MAD图与该材料的已知有效参考MAD图进行比较。如果成像系统的MAD图偏离(一个或多个)参考MAD图,则例如可以重新校准成像系统。
己经参考优选实施例描述了本发明。通过阅读和理解前述的详细描述,本领域技术人员可以想到各种修改和变型。旨在将本发明解释为包括所有这样的修改和替代,只要它们落入权利要求书及其等价方案的范围之内。

Claims (23)

1.一种方法,包括:
利用成像系统生成对象的谱投影数据,所述谱投影数据至少包括对应于第一能量范围的第一谱投影数据和对应于第二能量范围的第二谱投影数据,其中,所述第一能量范围和所述第二能量范围是不同的;
根据所述谱投影数据的组合来构建图像;
针对所述能量范围并且根据所述谱投影数据来构建基础图像的集合;
根据所述基础图像的集合构建多维直方图,其中,所述多维直方图包括至少两个轴,第一个轴对应于第一基础分量,并且第二个轴对应于第二基础分量,并且所述多维直方图包括聚类的集合,针对在谱投影数据中表示的每种材料所述聚类的集合包括一个聚类;并且
在视觉上同时呈现所述图像和所述多维直方图。
2.根据权利要求1所述的方法,还包括:
将所显示的多维直方图与针对已知正常组织的多维直方图进行比较;并且
响应于在所显示的多维直方图中的聚类在针对已知正常组织的多维直方图中不存在,识别所显示的多维直方图中存在异常。
3.根据权利要求2所述的方法,还包括:
将所显示的多维直方图与针对已知异常组织的所述多维直方图进行比较;并且
通过将所述聚类与所述多维直方图中对应于所述已知异常组织的聚类进行匹配来识别所述已知异常组织的存在。
4.根据权利要求1至3中的任一项所述的方法,还包括:
将所显示的多维直方图与针对已知异常组织的多维直方图进行比较;并且
响应于所显示的多维直方图中的聚类在针对已知正常组织的所述多维直方图中不存在,识别所显示的多维直方图中存在所述己知异常组织。
5.根据权利要求1至3中的任一项所述的方法,还包括:
将所显示的多维直方图与针对已知正常组织的多维直方图进行比较;并且
识别所显示的多维直方图中的在针对已知正常组织的所述多维直方图中不存在的聚类。
6.根据权利要求1至5中的任一项所述的方法,还包括:
接收输入,所述输入识别针对所述多维直方图中的聚类的感兴趣区域;并且
更新所显示的图像以突出显示所述图像中的与所述多维直方图中的识别出的感兴趣区域中的材料相对应的一组像素或体素。
7.根据权利要求1至6中的任一项所述的方法,还包括:
接收识别所述图像中的感兴趣区域的输入;并且
更新所显示的多维直方图以突出显示与所述图像中的识别出的感兴趣区域中的材料相对应的所述聚类的集合中的聚类。
8.根据权利要求7所述的方法,还包括:
将突出显示的所述多维直方图与针对不同材料的多维直方图的库进行比较;
将突出显示的所述聚类与来自所述库的多维直方图中的聚类匹配;并且
基于与匹配的聚类相对应的所述多维直方图来识别所述材料的类型。
9.根据权利要求8所述的方法,还包括:
针对所述材料类型创建材料分解;并且
使用所述材料分解来构建所述材料类型的谱图像。
10.根据权利要求1所述的方法,还包括:
将显示的多维直方图与针对多种材料的多维直方图的库进行比较;并且
通过将所述聚类与来自多维直方图的所述库中的直方图中的一个中的聚类进行匹配来识别由所述聚类表示的材料的类型。
11.根据权利要求10所述的方法,其中,所述聚类对应于已知材料,还包括:
将所述聚类标记为表示所述已知材料;并且
将所显示的多维直方图保存到所述库中。
12.一种系统(100),包括
存储器(128),其存储指令(132),所述包括谱CT指纹模块(134);以及
处理器,其被配置为至少执行与所述谱CT指纹模块相对应的指令,其使得所述处理器:
生成对象的谱投影数据,所述谱投影数据至少包括对应于第一能量范围的第一谱投影数据和对应于第二能量范围的第二谱投影数据,其中,所述第一能量范围和所述第二能量范围是不同的;
根据所述谱投影数据的组合来构建图像;针对所述能量范围并且根据所述谱投影数据来构建基础图像的集合;
根据所述基础图像的集合构建多维直方图,其中,所述多维直方图包括至少两个轴,第一个轴对应于第一基础分量,并且第二个轴对应于第二基础分量,并且所述多维直方图包括聚类的集合,针对在谱投影数据中表示的每种材料所述聚类的集合包括一个聚类;并且
在视觉上同时呈现所述图像和所述多维直方图。
13.根据权利要求12所述的系统,其中,所述处理器:接收输入,所述输入识别针对所述多维直方图中的聚类的感兴趣区域;并且更新显示的图像以突出显示所述图像中的与所述多维直方图中的识别出的感兴趣区域中的材料相对应的一组像素或体素。
14.根据权利要求12至13中的任一项所述的系统,其中,所述处理器:接收输入,所述输入识别所述图像中的感兴趣区域;并且更新所显示的多维直方图以突出显示与所述图像中的识别出的感兴趣区域中的材料相对应的聚类的集合中的聚类。
15.根据权利要求14所述的系统,其中,所述处理器:将突出显示的所述多维直方图与针对不同材料的多维直方图的库进行比较;将突出显示的所述聚类与来自所述库的多维直方图中的聚类进行匹配;并且基于所述多维直方图来识别与匹配的聚类相对应的所述材料的类型。
16.根据权利要求15所述的系统,其中,所述处理器:创建针对材料类型的材料分解;并且使用所述材料分解来构建所述材料类型的谱图像。
17.根据权利要求12至16中的任一项所述的系统,其中,所述处理器:将所显示的多维直方图与针对已知正常组织的多维直方图进行比较;并且响应于所显示的多维直方图中的聚类在针对已知正常组织的所述多维直方图中不存在而识别所显示的多维直方图中存在异常。
18.根据权利要求17所述的系统,其中,所述处理器:将所显示的多维直方图与针对已知异常组织的所述多维直方图进行比较;并且通过将所述聚类与所述多维直方图中对应于所述已知异常组织的聚类进行匹配来识别所述已知异常组织的存在。
19.根据权利要求12至16中的任一项所述的系统,其中,所述处理器:将所显示的多维直方图与针对已知异常组织的多维直方图进行比较;并且响应于所显示的多维直方图中的聚类在针对已知正常组织的所述多维直方图中不存在而识别在所显示的多维直方图中存在所述己知异常组织。
20.根据权利要求12至16中的任一项所述的系统,其中,所述处理器:将所显示的多维直方图与针对已知正常组织的多维直方图进行比较;并且识别所显示的多维直方图中的在针对已知正常组织的所述多维直方图中不存在的聚类。
21.根据权利要求12所述的系统,其中,所述处理器:将所显示的多维直方图与针对多种材料的多维直方图的库进行比较;并且通过将所述聚类与来自多维直方图的所述库中的直方图中的一个中的聚类进行匹配来识别由所述聚类表示的材料的类型。
22.根据权利要求21所述的系统,其中,所述聚类对应于己知材料,并且所述处理器将所述聚类标记为表示所述已知材料并且将显示的多维直方图保存到所述库中。
23.一种编码有计算机可执行指令的计算机可读介质,所述计算机可读指令在由计算机的处理器运行时,使所述处理器:
生成对象的谱投影数据,所述谱投影数据至少包括对应于第一能量范围的第一谱投影数据和对应于第二能量范围的第二谱投影数据,其中,所述第一能量范围和所述第二能量范围是不同的;
根据所述谱投影数据的组合来构建图像;
针对所述能量范围并且根据所述谱投影数据来构建基础图像的集合;
根据所述基础图像的集合构建多维直方图,其中,所述多维直方图包括至少两个轴,第一个轴对应于第一基础分量,并且第二个轴对应于第二基础分量,并且所述多维直方图包括聚类的集合,针对在谱投影数据中表示的每种材料所述聚类的集合包括一个聚类;并且
在视觉上同时呈现所述图像和所述多维直方图。
CN201780043267.0A 2016-07-13 2017-07-12 谱计算机断层扫描指纹识别 Active CN109688931B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201662361702P 2016-07-13 2016-07-13
US62/361702 2016-07-13
PCT/EP2017/067644 WO2018011321A1 (en) 2016-07-13 2017-07-12 Spectral computed tomography fingerprinting

Publications (2)

Publication Number Publication Date
CN109688931A true CN109688931A (zh) 2019-04-26
CN109688931B CN109688931B (zh) 2024-01-12

Family

ID=59315651

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201780043267.0A Active CN109688931B (zh) 2016-07-13 2017-07-12 谱计算机断层扫描指纹识别

Country Status (5)

Country Link
US (1) US10820865B2 (zh)
EP (1) EP3484369B1 (zh)
JP (1) JP7258744B2 (zh)
CN (1) CN109688931B (zh)
WO (1) WO2018011321A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111340127A (zh) * 2020-03-04 2020-06-26 清华大学 基于材料聚类的能谱ct迭代材料分解方法和装置

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11786201B1 (en) * 2022-06-14 2023-10-17 Curvebeam Ai Limited Method and system for material decomposition in dual- or multiple-energy x-ray based imaging

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070189443A1 (en) * 2004-11-19 2007-08-16 Walter Deborah J Detection of thrombi in ct using energy discrimination
CN101166469A (zh) * 2005-04-26 2008-04-23 皇家飞利浦电子股份有限公司 用于光谱ct的双层探测器
CN101495039A (zh) * 2005-09-22 2009-07-29 皇家飞利浦电子股份有限公司 用于能谱ct的定量物质分解
US20100128844A1 (en) * 2008-11-26 2010-05-27 Brian Thomsen System and Method for Material Segmentation Utilizing Computed Tomography Scans
US20100328313A1 (en) * 2009-06-26 2010-12-30 Kabushiki Kaisha Toshiba Advanced clustering method for material separation in dual energy ct
CN102076263A (zh) * 2008-06-30 2011-05-25 皇家飞利浦电子股份有限公司 频谱ct
JP2011110245A (ja) * 2009-11-27 2011-06-09 Ge Medical Systems Global Technology Co Llc 画像表示装置、x線ct装置およびプログラム
US20120106816A1 (en) * 2010-10-29 2012-05-03 Bruno Kristiaan Bernard De Man System and method for material decomposition optimization in image domain
CN104700390A (zh) * 2013-12-09 2015-06-10 通用电气公司 双能ct造影剂增强扫描图像中的钙化部分识别方法
WO2016059527A2 (en) * 2014-10-13 2016-04-21 Koninklijke Philips N.V. Spectral imaging
CN105559813A (zh) * 2014-11-04 2016-05-11 株式会社东芝 医用图像诊断装置以及医用图像处理装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104700389B (zh) * 2013-12-09 2019-08-13 通用电气公司 双能ct扫描图像中的对象识别方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070189443A1 (en) * 2004-11-19 2007-08-16 Walter Deborah J Detection of thrombi in ct using energy discrimination
CN101166469A (zh) * 2005-04-26 2008-04-23 皇家飞利浦电子股份有限公司 用于光谱ct的双层探测器
CN101495039A (zh) * 2005-09-22 2009-07-29 皇家飞利浦电子股份有限公司 用于能谱ct的定量物质分解
CN102076263A (zh) * 2008-06-30 2011-05-25 皇家飞利浦电子股份有限公司 频谱ct
US20100128844A1 (en) * 2008-11-26 2010-05-27 Brian Thomsen System and Method for Material Segmentation Utilizing Computed Tomography Scans
US20100328313A1 (en) * 2009-06-26 2010-12-30 Kabushiki Kaisha Toshiba Advanced clustering method for material separation in dual energy ct
JP2011110245A (ja) * 2009-11-27 2011-06-09 Ge Medical Systems Global Technology Co Llc 画像表示装置、x線ct装置およびプログラム
US20120106816A1 (en) * 2010-10-29 2012-05-03 Bruno Kristiaan Bernard De Man System and method for material decomposition optimization in image domain
CN104700390A (zh) * 2013-12-09 2015-06-10 通用电气公司 双能ct造影剂增强扫描图像中的钙化部分识别方法
WO2016059527A2 (en) * 2014-10-13 2016-04-21 Koninklijke Philips N.V. Spectral imaging
CN105559813A (zh) * 2014-11-04 2016-05-11 株式会社东芝 医用图像诊断装置以及医用图像处理装置

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111340127A (zh) * 2020-03-04 2020-06-26 清华大学 基于材料聚类的能谱ct迭代材料分解方法和装置
CN111340127B (zh) * 2020-03-04 2023-01-10 清华大学 基于材料聚类的能谱ct迭代材料分解方法和装置

Also Published As

Publication number Publication date
US10820865B2 (en) 2020-11-03
EP3484369B1 (en) 2023-05-31
JP7258744B2 (ja) 2023-04-17
EP3484369A1 (en) 2019-05-22
US20190290226A1 (en) 2019-09-26
WO2018011321A1 (en) 2018-01-18
JP2019527885A (ja) 2019-10-03
CN109688931B (zh) 2024-01-12

Similar Documents

Publication Publication Date Title
RU2612575C2 (ru) Обработка изображений для спектральной компьютерной томографии
EP3213298B1 (en) Texture analysis map for image data
US11813101B2 (en) Image quality improved virtual non-contrast images generated by a spectral computed tomography (CT) scanner
CN105793894B (zh) 根据图像数据来进行骨骼分割
CN107209946A (zh) 图像数据分割和显示
KR20160119540A (ko) 단층 촬영 장치 및 그에 따른 단층 영상 처리 방법
CN109770933A (zh) 通过三维定位改进图像质量的系统和方法
CN109688931A (zh) 谱计算机断层扫描指纹识别
US20170181714A1 (en) Contrast agent dose simulation
JP7489387B2 (ja) 3d仮想内視鏡レンダリング
US11935162B2 (en) Protocol-dependent 2-D pre-scan projection image based on 3-D pre-scan volumetric image data
US11361415B2 (en) Material-selective adaptive blending of volumeiric image data
US20230148983A1 (en) Suppression of motion artifacts in computed tomography imaging
Pedrero Pérez Truncation artifact correction for micro-CT scanners

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