CN107978003B - 一种ct图像金属伪影的处理方法及装置 - Google Patents

一种ct图像金属伪影的处理方法及装置 Download PDF

Info

Publication number
CN107978003B
CN107978003B CN201711034156.1A CN201711034156A CN107978003B CN 107978003 B CN107978003 B CN 107978003B CN 201711034156 A CN201711034156 A CN 201711034156A CN 107978003 B CN107978003 B CN 107978003B
Authority
CN
China
Prior art keywords
image
metal
original image
region
original
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
Application number
CN201711034156.1A
Other languages
English (en)
Other versions
CN107978003A (zh
Inventor
楼珊珊
任亮
逄岭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Neusoft Medical Systems Co Ltd
Original Assignee
Neusoft Medical Systems Co Ltd
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 Neusoft Medical Systems Co Ltd filed Critical Neusoft Medical Systems Co Ltd
Priority to CN201711034156.1A priority Critical patent/CN107978003B/zh
Publication of CN107978003A publication Critical patent/CN107978003A/zh
Priority to US16/172,894 priority patent/US10783676B2/en
Application granted granted Critical
Publication of CN107978003B publication Critical patent/CN107978003B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/20Image enhancement or restoration using local operators
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • 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]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20224Image subtraction

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Abstract

本发明实施例公开了一种CT图像金属伪影的处理方法,在去除CT图像的金属伪影时,根据金属介入前第一CT扫描的第一原始图像和金属介入后第二CT扫描的第二原始图像共同确定第二原始图像中的金属区域的第三图像,以获得金属区域对应的金属通道位置。之后,建立第二原始图像的模型图像,将原始生数据中所述金属通道位置处的数据替换为第二原始图像的模型图像对应的模型生数据中金属通道位置处的数据,并获得替换后的原始生数据对应的第四图像。最后,将第三图像融合至第四图像中,以获得去除金属伪影后的第五图像。这样,可以准确地获得金属介入后第二CT扫描的第二原始图像中的金属区域,从而提高金属区域的提取精度,提高去除金属伪影的CT图像的质量。

Description

一种CT图像金属伪影的处理方法及装置
技术领域
本发明涉及医学影像的图像处理技术领域,尤其涉及一种CT图像金属伪影的处理方法及装置。
背景技术
CT(Computed Tomography,电子计算机断层扫描)是医学领域中常用的扫描设备,其通过X射线对被扫描对象进行扫描,X射线经过被扫描对象衰减后,被检测器检测并转换为数据信号,该数据信号被称为原始生数据,通过原始生数据重建出CT图像。
在临床应用时,被扫描对象通常是患者,在CT扫描时会有金属介入患者的情形,例如需要在患者体内置入手术探针等情形,在被扫描的物体中包含金属时,由于金属对X射线的高衰减特性,导致获得的数据信号存在误差,这会导致重建出来图像存在金属伪影,金属伪影为从金属部位发射出的充满图像的放射状细纹,会降低CT图像的质量,需要将金属伪影去除掉。
目前,主要是基于滤波反投影重建方法,将重建的CT图像上的金属伪影去除,在该方法中,对通过原始生数据重建的CT图像进行阈值分割,提取出金属图像,而后通过修正原始生数据,再次重新建像获得去除金属伪影的图像。而由于不同患者的身体情况并不相同,通过单一的阈值分割的方法提取的金属图像,会导致去除金属伪影后金属图像部分的不准确性,影响CT图像的质量。
发明内容
本发明提供了一种CT图像金属伪影的处理方法及装置,以准确地确定金属区域的位置,提高金属区域提取的准确性,进而提高去除金属伪影后图像的准确性,保证CT图像的质量。
本发明提供了一种CT图像金属伪影的处理方法,所述方法包括:
在预设的扫描重建设置条件下,获得金属介入前第一扫描的第一原始图像;以及,获得金属介入后第二扫描的原始生数据和第二原始图像;
根据所述第一原始图像和所述第二原始图像,获得所述第二原始图像中金属区域的第三图像;
根据所述第三图像,获得所述金属区域对应的金属通道位置;
将所述第二原始图像中所述金属区域的像素值设置为预设值,获得模型图像,所述预设值不同于所述第二原始图像中其他各物质的像素值,并获得所述模型图像对应的模型生数据;
将所述第二扫描的原始生数据中所述金属通道位置处的数据替换为所述模型生数据中对应位置处的数据;
获得替换后的原始生数据对应的第四图像;
将所述第三图像和所述第四图像融合,获得目标图像。
可选地,所述根据所述第一原始图像和所述第二原始图像,获得所述第二原始图像中金属区域的第三图像,包括:
以所述第一原始图像为背景图像,获得所述第二原始图像与所述第一原始图像的差异图像;
通过像素过滤,保留所述差异图像中的金属区域,以获得第三图像。
可选地,所述将所述第二原始图像中金属区域的像素值设置为预设值,获得模型图像,包括:
对所述第一原始图像和所述第二原始图像进行像素对此,在所述第二原始图像中确定出金属区域;
将所述金属区域的像素值设置为预设值。
可选地,所述获得模型图像,还包括:
将所述第二原始图像中的其他区域区分为不同的物质区域,所述其他区域为第二原始图像中金属区域之外的部分;
将所述其他区域中各物质区域的各像素值分别置换为与各物质对应的特定像素值;或者,若所述物质区域中包括骨头区域,则将所述其他区域中所述骨头区域之外的各物质区域的各像素值分别置换为与各物质对应的特定像素值。
可选地,所述将所述第二原始图像中的其他区域区分为不同的物质区域,包括:
通过m个像素阈值将所述其他区域划分为n个不同的物质区域,其中,n=m+1,m≥1。
可选地,所述特定像素值为各物质区域像素值的均匀值。
一种CT图像金属伪影的处理装置,包括:
图像及数据获取单元,用于在预设的扫描重建设置条件下,获得金属介入前第一CT扫描的第一原始图像;以及,获得金属介入后第二CT扫描的原始生数据和第二原始图像;
金属图像获取单元,根据所述第一原始图像和所述第二原始图像,获得所述第二原始图像中金属区域的第三图像;
金属通道位置确定单元,用于根据所述第三图像,获得所述金属区域对应的金属通道位置;
模型图像获取单元,用于将所述第二原始图像中所述金属区域的像素值设置为预设值,获得模型图像,所述预设值不同于所述第二原始图像中其他各物质的像素值,并获得所述模型图像对应的模型生数据;
金属区域数据替换单元,用于将所述第二扫描的原始生数据中所述金属通道位置处的数据替换为所述模型生数据中对应位置处的数据;
替换图像获取单元,用于获得替换后的原始生数据对应的第四图像;
目标图像获取单元,用于将所述第三图像和所述第四图像融合,获得目标图像。
可选地,所述金属图像获取单元包括:
差异图像获取单元,用于以所述第一原始图像为背景图像,获得所述第二原始图像与所述第一原始图像的差异图像;
过滤单元,用于通过像素过滤,保留所述差异图像中的金属区域,以获得第三图像。
可选地,所述模型图像获取单元包括:
金属区域设置单元,用于通过对所述第一原始图像和所述第二原始图像的进行像素对比,在所述第二原始图像中确定出金属区域,并将所述金属区域的像素值设置为预设值。
可选地,所述模型图像获取单元还包括:
区域划分单元,用于将所述第二原始图像中的其他区域区分为不同的物质区域,所述其他区域为第二原始图像中金属区域之外的部分;
像素置换单元,用于将所述其他区域中各物质区域的各像素值分别置换为与各物质对应的特定像素值;或者,若所述物质区域中包括骨头区域,则将所述其他区域中所述骨头区域之外的各物质区域的各像素值分别置换为与各物质对应的特定像素值。
可选地,所述区域划分单元中,通过m个像素阈值将所述其他区域划分为n个不同的物质区域,其中,n=m+1,m≥1。
可选地,所述特定像素值为各物质区域像素值的均匀值。
本发明实施例提供的一种CT图像金属伪影的处理方法,在去除CT图像的金属伪影时,根据金属介入前第一CT扫描的第一原始图像和金属介入后第二CT扫描的第二原始图像共同确定第二原始图像中的金属区域的第三图像,以获得金属区域对应的金属通道位置。之后,建立第二原始图像的模型图像,将原始生数据中所述金属通道位置处的数据替换为第二原始图像的模型图像对应的模型生数据中金属通道位置处的数据,并获得替换后的原始生数据对应的第四图像。最后,将第三图像融合至第四图像中,以获得去除金属伪影后的目标图像。这样,通过将金属介入前第一CT扫描的第一原始图像和金属介入后的第二原始图像,可以准确地获得金属区域以及对应的金属通道位置,从而提高金属区域的提取精度,将原始生数据中的金属通道位置对应的数据替换之后,可以有效地去除金属伪影,进而提高去除金属伪影的CT图像的质量。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的CT图像金属伪影的处理方法的应用场景示意图;
图2为本发明实施例的CT图像金属伪影的处理方法的流程示意图;
图3为本发明实施例的CT图像金属伪影的处理装置的结构框图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
正如背景技术中的描述,基于滤波反投影重建方法去除CT图像上的金属伪影,是根据金属区域和除金属区域的其他物质区域灰度等级的不同而设置阈值,从而将CT图像上的各像素点进行划分,以区分出金属区域和其他物质区域。但是,由于不同患者的身体情况的不相同,导致金属区域和其他物质区域灰度差异不明显,从而在设置阈值时,可能会出现将不属于金属区域的其他物质区域划分进来或将属于金属区域的部分区域遗漏的情况,造成提取的金属区域不准确,进而导致去除金属伪影后金属图像部分的不准确性,影响CT图像的质量。为此,本发明旨在寻求一种新的方法,可以有效地提高提取金属区域的准确性,进而提高去除金属伪影后图像的准确性,保证CT图像的质量。
举例说明,本发明实施例可以应用到如图1所示的场景。在所述场景下,用户101可以与用户终端102进行交互,也即,用户101可以操作用户终端102,用户终端102可以响应用户101的操作执行相应的动作以及向用户101呈现图像。
具体地,用户可以预先在用户终端102设置扫描重建设置条件,用户终端102可以在该条件下获得金属介入前第一CT扫描的第一原始图像,以及获得金属介入后第二CT扫描的原始生数据和第二原始图像。然后,用户终端102可以根据第一原始图像和第二原始图像,可以通过图像处理获得第二原始图像中的金属区域的第三图像。之后,用户终端102可以根据第三图像,获得金属区域对应的金属通道位置。再后,用户终端102可以将所述第二原始图像中所述金属区域的像素值设置为预设值,获得模型图像,所述预设值不同于所述第二原始图像中其他各物质的像素值,并获得所述模型图像对应的模型生数据。最后,用户终端102可以将所述第三图像融合至所述第四图像中,以获得目标图像,并将目标图像呈现给用户,目标图像为去除伪影后的CT图像,以便用户101根据用户终端102呈现的目标图像对患者进行诊断。
需要注意的是,上述应用场景仅是为了便于理解本发明而示出,本发明的实施方式在此方面不受任何限制。相反,本发明的实施方式可以应用于适用的任何场景。
下面结合附图,详细说明本申请实施例如何处理CT图像金属伪影的。
参见图2,图2为本申请实施例提供的一种CT图像金属伪影的处理方法的流程图,该方法包括:
S201、在预设的扫描重建设置条件下,获得金属介入前第一CT扫描的第一原始图像;以及,获得金属介入后第二CT扫描的原始生数据和第二原始图像。
在临床应用时,患者在被CT进行扫描时,有时会由于金属介入患者而使重建出来的CT图像存在金属伪影,影响CT图像质量,影响医务人员根据CT图像对患者进行诊断。在本实施例中,为了去除CT图像中的金属伪影,首先需要确定金属区域,由于金属介入前后,利用CT扫描得到的图像会有所差异,因此,为了确定金属区域以及后续伪影的去除,用户终端需要获得金属介入前第一CT扫描的第一原始图像,以及获得金属介入后第二CT扫描的原始生数据和第二原始图像。
需要说明的是,为了使用户终端获得的金属介入前第一CT扫描的第一原始图像和金属介入后第二CT扫描的原始生数据和第二原始图像具有可比性,第一CT扫描和第二CT扫描需要在相同的条件下进行,即预设的第一CT扫描和第二CT扫描的扫描重建设置条件应该是相同的。其中,扫描重建设置条件可以包括通道的数量、层数、投影数量等,通过对扫描重建设置条件的预先设置,使两次扫描在同一扫描重建设置条件下进行。
例如,在患者体内置入手术探针的情形下,在通过手术向患者体内植入探针之前,在预设的扫描重建设置条件下,利用CT对患者进行扫描,用户终端得到患者的第一原始图像,在通过手术向患者体内植入探针之后,在与上述预设的扫描重建设置条件相同的扫描重建设置条件下,利用CT对患者进行扫描,用户终端得到患者的原始生数据和第二原始图像,以便在后续过程进行利用。
S202、根据所述第一原始图像和所述第二原始图像,获得所述第二原始图像中的金属区域的第三图像。
用户终端可以对第一原始图像和第二原始图像进行图像处理,从而获得第二原始图像中的金属区域的第三图像。
由于第二原始图像与第一原始图像相比,其区别在于第二原始图像中包括金属区域以及由于金属的存在而出现的金属伪影,因此,用户终端在对第一原始图像和第二原始图像进行图像处理时,可以以所述第一原始图像为背景图像,获得第二原始图像与背景图像的差异图像,该差异图像即包括了金属区域和由于金属的存在而出现的金属伪影;然后,用户终端可以利用金属区域与金属伪影像素值的不同,通过像素过滤的方法,将差异图像中的金属伪影过滤掉,保留差异图像中的金属区域,以获得只包括金属区域的第三图像。
S203、根据所述第三图像,获得所述金属区域对应的金属通道位置。
在图像处理中,图像经过正投影可以得到与该图像对应的生数据,因此,用户终端可以对第三图像进行正投影,得到第三图像对应的生数据,即金属区域的生数据,从而根据金属区域的生数据与金属通道位置的对应关系,获得金属区域在原始生数据中对应的通道位置,即金属通道位置,该金属通道位置涵盖整个金属区域,可以具体为金属区域数据对应的通道起始位置和通道终止位置。
S204、将所述第二原始图像中所述金属区域的像素值设置为预设值,获得模型图像,所述预设值不同于所述第二原始图像中其他各物质的像素值,并获得所述模型图像对应的模型生数据。
该步骤中,对第二原始图像进行处理,获得模型图像,在一些实施例中,可以仅将第二原始图像中金属区域的像素值设置为预设值,预设值不同于所述第二原始图像中其他各物质的像素值,来获得模型图像,这样,模型图像中金属区域的像素值被均一化,相应地,该模型图像对应的生数据中金属通道位置的数据被理想化。
在另外一些实施例中,可以将第二原始图像中金属区域的像素值设置为预设值,同时,将金属区域之外的其他区域进行像素处理,其他区域也即非金属区域,可以将其他区域中各物质区域的各像素值分别置换为与各物质对应的特定像素值;或者,若所述物质区域中包括骨头区域,则将其他区域中所述骨头区域之外的各物质区域的各像素值分别置换为与各物质对应的特定像素值。这样,可以消除非金属物质的噪声,尤其是非骨头物质的噪声,提高模型图像的质量,进而提高模型图像对应的生数据的质量。
以下将具体描述模型图像的获取方法,具体的,
在步骤S401,将所述第二原始图像中所述金属区域的像素值设置为预设值,所述预设值不同于所述第二原始图像中其他各物质的像素值。
在建立第二原始图像的模型图像时,用户终端可以通过在第二原始图像中确定出金属区域,从而将第二原始图像分为金属区域和非金属区域,将金属区域设置某个特定的像素值表示。其中,由于第二原始图像中包括金属区域,而第一原始图像(即背景图像)不包含金属区域,因此,用户终端可以通过对所述第一原始图像和所述第二原始图像进行像素对比,在第二原始图像中确定出金属区域,然后用户终端将所确定出的金属区域的像素值设置为预设值,从而建立第二原始图像的模型图像中的金属区域。其中,所述预设值例如可以是水对应的像素值,或其他的区别于第二原始图像中各物质的像素值。
在步骤S402,将第二原始图像中的其他区域区分为不同的物质区域,其他区域为第二原始图像中金属区域之外的部分。
其他区域也即确定金属区域之后,第二原始图像中的非金属区域,在建立第二原始图像的模型图像时,用户终端可以将非金属区域进行区分,即可以将金属区域之外的部分区分出不同的物质区域,将各个物质区域采用与各个物质对应的像素值进行表示。具体地,用户终端可以通过像素阈值的设置,将所述第二原始图像中非金属区域区分为不同的物质区域。
其中,所述像素阈值可以根据第二原始图像非金属区域包括的各个物质进行设置。一般来说,设置像素阈值的数量m可以依据第二原始图像中非金属区域包含的物质种类数量n而确定,例如,第二原始图像中非金属区域包括的物质的种类数n为3,那么,可以设置n-1=2个像素阈值进行区分。每个像素阈值的预设值可以根据第二原始图像中非金属区域各物质的像素分布范围确定,该预设值可以根据经验值来确定,通过n-1个像素阈值,将非金属区域的像素值划分为m个区域,每个区域代表一种物质,从而,将非金属区域划分为m种物质对应的m个物质区域。
在一些情况下,所述第二原始图像中非金属区域可能包括空气区域、软组织区域和骨头区域,用户终端需要对空气区域、软组织区域和骨头区域这三种物质区域进行区分。依据上述的像素阈值设定方法,可以设置两个像素阈值,即第一阈值和第二阈值,用户终端可以通过所述第一阈值区分软组织区域和空气区域,通过所述第二阈值区分软组织区域和骨头区域。在步骤S403,进行其他区域像素值的置换。
在一些实施例中,可以将其他区域中各物质区域的各像素值分别置换为与各物质对应的特定像素值。在这些实施例中,将其他区域中的各物质区域的像素值进行置换,从而得到模型图像,该模型图像具有更低的噪声,使得模型图像对应的生数据中金属通道位置的数据被理想化,以便于后续利用模型图像对应的生数据进行原始生数据的替换,进而获得去除伪影的图像。
在另一些实施例中,若物质区域中包括骨头区域,可以将其他区域中骨头区域之外的各物质区域的各像素值分别置换为与各物质对应的特定像素值。由于骨头物质的特殊性,骨头区域可以保持原有的像素值,并不进行像素值的置换,来获得模型图像,该模型图像同样具有较低的噪声。
在利用像素阈值将第二原始图像中非金属区域区分为不同的物质区域之后,由于即使在同一物质区域中各像素值也可能有所不同,为了使同一物质区域对应的像素值均匀,可以将非金属区域的各个物质区域中各像素值分别被置换为与各物质对应的特定像素值,例如,将所述软组织区域和空气区域中各像素值分别置换为软组织对应的特定像素值和空气对应的特定像素值。
以上述第二原始图像中非金属区域包括空气区域、软组织区域和骨头区域的情况为例,当利用第一阈值和第二阈值区分出空气区域、软组织区域和骨头区域这三个区域后,用户终端可以将空气区域的各像素值置换为空气区域中空气对应的特定像素值,将软组织区域的各像素值置换为软组织区域中软组织对应的特定像素值。进一步地地,还可以将骨头区域的各像素值置换为骨头区域中骨头对应的特定像素值。
其中,所述特定像素值可以为各区域像素值的均匀值。
该均匀值可以为各物质区域的像素值的均值,
综上,用户终端可以将空气区域的各像素值置换为空气区域像素值的均值,将软组织区域的各像素值置换为软组织区域像素值的均值,将骨头区域域的各像素值置换为骨头区域像素值的均值或保留骨头区域的像素值。上述特定像素值除了使用各区域像素值的均值外,还可以使用各区域像素值的正态分布的中央值或者各区域像素值比较集中靠近的一个值等,也就是说,特定像素值为代表各物质区域的像素值的均匀值,该均匀值能够表征该物质区域中的物质在图像中体现的真实像素值,各区域的均值、正态分布值或集中分布值等都可以体现各区域中的物质,这样,在进行替换之后,各区域中的像素值即能准确表征该物质同时该区域图像为均匀图像,从而,获得较为理想的模型图像。
在步骤S404,获得所述模型图像对应的模型生数据。
通过将模型数据进行正投影,获得模型图像对应的生数据,记做模型生数据。
由于该模型图像按照不同物质区域进行均匀值替换,替换后的各物质区域都是均匀的,模型图像不存在伪影的影响而造成非均匀的部分,因此,该模型图像近似为一种不包含伪影的理想图像,那么,基于该模型图像进行正投影得到的模型生数据也是不包含伪影数据的生数据。
S205、将所述第二扫描的原始生数据中所述金属通道位置处的数据替换为所述模型生数据中对应位置处的数据,并获得替换后的原始生数据对应的第四图像。
第二扫描的原始生数据中包括了第二原始图像的金属及伪影对应的数据,所述第二原始图像的金属及伪影对应的数据位于金属通道位置处,如果可以将原始生数据中第二原始图像中金属及伪影对应的生数据去除,将得到不包括金属和金属伪影的第二原始图像。因此,为了实现得到不包括金属和金属伪影的图像,在本实施例中,用户终端可以根据S203获得的第二原始图像生数据中金属通道的位置,将第二原始图像生数据中所述金属通道位置处的数据替换为所述模型生数据中所述金属通道位置处的数据,由于模型图像中的金属区域的像素值被设置为预设值,是去除了金属影响后的数据,因此替换后的原始生数据是去除了金属影响后的生数据。
在图像处理中,生数据经过反投影可以得到与该生数据对应的图像,因此,用户终端可以对替换后的原始生数据进行反投影,获得替换后的原始生数据对应的第四图像,该第四图像是去除了金属和金属伪影后的图像。
S206、将所述第三图像融合至所述第四图像中,以获得目标图像。
其中,第三图像为第二原始图像中的金属区域的图像,第四图像为去除了金属和金属伪影后的图像,用户终端可以通过图像处理的方法将第三图像融合至所述第四图像中,具体地,例如用户终端可以是将第四图像中与第三图像中金属区域相同位置的像素值替换为第三图像中金属区域的像素值,从而得到包括金属区域且去除了金属伪影的目标图像,实现了去除CT图像中金属伪影的目的。
本发明实施例提供的一种CT图像金属伪影的处理方法,在去除CT图像的金属伪影时,根据金属介入前第一CT扫描的第一原始图像和金属介入后第二CT扫描的第二原始图像共同确定第二原始图像中的金属区域的第三图像,并获得金属区域对应的金属通道位置。之后,建立第二原始图像的模型图像,将第二原始图像生数据中所述金属通道位置处的数据替换为第二原始图像的模型图像对应的模型生数据中对应位置处的数据,并获得替换后的原始生数据对应的第四图像。最后,将第三图像融合至第四图像中,以获得去除金属伪影后的目标图像。这样,通过将金属介入前第一CT扫描的第一原始图像和金属介入后的第二原始图像,可以准确地获得金属区域以及对应的金属通道位置,从而提高金属区域的提取精度,将原始生数据中的金属通道位置对应的数据替换之后,可以有效地去除金属伪影,进而提高去除金属伪影的CT图像的质量。
基于上述实施例提供的一种CT图像金属伪影的处理方法,本发明实施例提供了一种CT图像金属伪影的处理装置。参见图3,示出了一种CT图像金属伪影的处理装置,该装置包括:
图像及数据获取单元301,用于在预设的扫描重建设置条件下,获得金属介入前第一CT扫描的第一原始图像;以及,获得金属介入后第二CT扫描的原始生数据和第二原始图像;
金属图像获取单元302,根据所述第一原始图像和所述第二原始图像,获得所述第二原始图像中金属区域的第三图像;
金属通道位置确定单元303,用于根据所述第三图像,获得所述金属区域对应的金属通道位置;
模型图像获取单元304,用于将所述第二原始图像中所述金属区域的像素值设置为预设值,获得模型图像,所述预设值不同于所述第二原始图像中其他各物质的像素值,并获得所述模型图像对应的模型生数据;
金属区域数据替换单元305,用于将所述第二扫描的原始生数据中所述金属通道位置处的数据替换为所述模型生数据中对应位置处的数据;
替换图像获取单元306,用于获得替换后的原始生数据对应的第四图像;
目标图像获取单元307,用于将所述第三图像和所述第四图像融合,获得目标图像。
进一步地,所述金属图像获取单元302包括:
差异图像获取单元,用于以所述第一原始图像为背景图像,获得所述第二原始图像与所述第一原始图像的差异图像;
过滤单元,用于通过像素过滤,保留所述差异图像中的金属区域,以获得第三图像。
进一步地,所述模型图像获取单元304包括:
金属区域设置单元,用于通过对所述第一原始图像和所述第二原始图像的进行像素对比,在所述第二原始图像中确定出金属区域,并将所述金属区域的像素值设置为预设值。
进一步地,所述模型图像获取单元304还包括:
区域划分单元,用于将所述第二原始图像中的其他区域区分为不同的物质区域,所述其他区域为第二原始图像中金属区域之外的部分;
像素置换单元,用于将所述其他区域中各物质区域的各像素值分别置换为与各物质对应的特定像素值;或者,若所述物质区域中包括骨头区域,则将所述其他区域中所述骨头区域之外的各物质区域的各像素值分别置换为与各物质对应的特定像素值。
进一步地,所述区域划分单元中,通过m个像素阈值将所述其他区域划分为n个不同的物质区域,其中,n=m+1,m≥1。
进一步地,所述特定像素值为各物质区域像素值的均匀值。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述得比较简单,相关之处参见方法实施例的部分说明即可。以上所描述的系统实施例仅仅是示意性的,其中所述作为分离部件说明的模块或单元可以是或者也可以不是物理上分开的,作为模块或单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

Claims (12)

1.一种CT图像金属伪影的处理方法,其特征在于,所述方法包括:
在预设的扫描重建设置条件下,获得金属介入前第一CT扫描的第一原始图像;以及,获得金属介入后第二CT扫描的原始生数据和第二原始图像;
根据所述第一原始图像和所述第二原始图像,获得所述第二原始图像中金属区域的第三图像;
根据所述第三图像,获得所述金属区域对应的金属通道位置,所述金属通道位置处的数据包括所述金属区域及金属伪影对应的生数据;
将所述第二原始图像中所述金属区域的像素值设置为预设值,获得模型图像,所述预设值不同于所述第二原始图像中其他各物质的像素值,并获得所述模型图像对应的模型生数据;
将所述第二CT扫描的原始生数据中所述金属通道位置处的数据替换为所述模型生数据中对应位置处的数据;
获得替换后的原始生数据对应的第四图像;
将所述第三图像和所述第四图像融合,获得目标图像。
2.根据权利要求1所述的处理方法,其特征在于,所述根据所述第一原始图像和所述第二原始图像,获得所述第二原始图像中金属区域的第三图像,包括:
以所述第一原始图像为背景图像,获得所述第二原始图像与所述第一原始图像的差异图像;
通过像素过滤,保留所述差异图像中的金属区域,以获得第三图像。
3.根据权利要求1所述的处理方法,其特征在于,所述将所述第二原始图像中金属区域的像素值设置为预设值,获得模型图像,包括:
对所述第一原始图像和所述第二原始图像进行像素对此,在所述第二原始图像中确定出金属区域;
将所述金属区域的像素值设置为预设值。
4.根据权利要求1或3所述的处理方法,其特征在于,所述获得模型图像,还包括:
将所述第二原始图像中的其他区域区分为不同的物质区域,所述其他区域为第二原始图像中金属区域之外的部分;
将所述其他区域中各物质区域的各像素值分别置换为与各物质对应的特定像素值;或者,若所述物质区域中包括骨头区域,则将所述其他区域中所述骨头区域之外的各物质区域的各像素值分别置换为与各物质对应的特定像素值。
5.根据权利要求4所述的处理方法,其特征在于,所述将所述第二原始图像中的其他区域区分为不同的物质区域,包括:
通过m个像素阈值将所述其他区域划分为n个不同的物质区域,其中,n=m+1,m≥1。
6.根据权利要求4所述的处理方法,其特征在于,所述特定像素值为各物质区域像素值的均匀值。
7.一种CT图像金属伪影的处理装置,其特征在于,包括:
图像及数据获取单元,用于在预设的扫描重建设置条件下,获得金属介入前第一CT扫描的第一原始图像;以及,获得金属介入后第二CT扫描的原始生数据和第二原始图像;
金属图像获取单元,根据所述第一原始图像和所述第二原始图像,获得所述第二原始图像中金属区域的第三图像;
金属通道位置确定单元,用于根据所述第三图像,获得所述金属区域对应的金属通道位置,所述金属通道位置处的数据包括所述金属区域及金属伪影对应的生数据;
模型图像获取单元,用于将所述第二原始图像中所述金属区域的像素值设置为预设值,获得模型图像,所述预设值不同于所述第二原始图像中其他各物质的像素值,并获得所述模型图像对应的模型生数据;
金属区域数据替换单元,用于将所述第二CT扫描的原始生数据中所述金属通道位置处的数据替换为所述模型生数据中对应位置处的数据;
替换图像获取单元,用于获得替换后的原始生数据对应的第四图像;
目标图像获取单元,用于将所述第三图像和所述第四图像融合,获得目标图像。
8.根据权利要求7所述的装置,其特征在于,所述金属图像获取单元包括:
差异图像获取单元,用于以所述第一原始图像为背景图像,获得所述第二原始图像与所述第一原始图像的差异图像;
过滤单元,用于通过像素过滤,保留所述差异图像中的金属区域,以获得第三图像。
9.根据权利要求7所述的装置,其特征在于,所述模型图像获取单元包括:
金属区域设置单元,用于通过对所述第一原始图像和所述第二原始图像的进行像素对比,在所述第二原始图像中确定出金属区域,并将所述金属区域的像素值设置为预设值。
10.根据权利要求7或9所述的装置,其特征在于,所述模型图像获取单元还包括:
区域划分单元,用于将所述第二原始图像中的其他区域区分为不同的物质区域,所述其他区域为第二原始图像中金属区域之外的部分;
像素置换单元,用于将所述其他区域中各物质区域的各像素值分别置换为与各物质对应的特定像素值;或者,若所述物质区域中包括骨头区域,则将所述其他区域中所述骨头区域之外的各物质区域的各像素值分别置换为与各物质对应的特定像素值。
11.根据权利要求10所述的装置,其特征在于,所述区域划分单元中,通过m个像素阈值将所述其他区域划分为n个不同的物质区域,其中,n=m+1,m≥1。
12.根据权利要求10所述的装置,其特征在于,所述特定像素值为各物质区域像素值的均匀值。
CN201711034156.1A 2017-10-30 2017-10-30 一种ct图像金属伪影的处理方法及装置 Active CN107978003B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201711034156.1A CN107978003B (zh) 2017-10-30 2017-10-30 一种ct图像金属伪影的处理方法及装置
US16/172,894 US10783676B2 (en) 2017-10-30 2018-10-29 Processing metal artifacts in CT image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711034156.1A CN107978003B (zh) 2017-10-30 2017-10-30 一种ct图像金属伪影的处理方法及装置

Publications (2)

Publication Number Publication Date
CN107978003A CN107978003A (zh) 2018-05-01
CN107978003B true CN107978003B (zh) 2021-03-30

Family

ID=62012769

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711034156.1A Active CN107978003B (zh) 2017-10-30 2017-10-30 一种ct图像金属伪影的处理方法及装置

Country Status (2)

Country Link
US (1) US10783676B2 (zh)
CN (1) CN107978003B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110547821B (zh) * 2018-05-31 2024-04-26 通用电气公司 追踪目标介入物的ct系统、方法及相关计算机可读介质
CN110097517B (zh) * 2019-04-28 2022-12-27 东软医疗系统股份有限公司 去除图像伪影的方法及装置
CN113554563B (zh) * 2021-07-23 2024-05-14 上海友脉科技有限责任公司 一种医学图像处理方法、介质及电子设备
CN113520441B (zh) * 2021-08-03 2022-11-29 浙江大学 消除ct高阻射物伪影干扰的组织显像方法及系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102567958A (zh) * 2010-12-31 2012-07-11 东软飞利浦医疗设备系统有限责任公司 一种去除ct图像中金属伪影的图像后处理方法
CN103971387A (zh) * 2013-01-29 2014-08-06 上海联影医疗科技有限公司 Ct图像重建方法
WO2015168147A1 (en) * 2014-04-29 2015-11-05 Carl Zeiss X-ray Microscopy, Inc. Segmentation and spectrum based metal artifact reduction method and system
CN105225208A (zh) * 2015-10-14 2016-01-06 上海联影医疗科技有限公司 一种计算机断层成像金属伪影校正方法及装置
US9349198B2 (en) * 2013-07-26 2016-05-24 General Electric Company Robust artifact reduction in image reconstruction
CN105701778A (zh) * 2016-01-11 2016-06-22 赛诺威盛科技(北京)有限公司 一种从ct图像中去除金属伪影的方法
CN106683143A (zh) * 2016-01-21 2017-05-17 上海联影医疗科技有限公司 图像金属伪影校正方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7835785B2 (en) * 2005-10-04 2010-11-16 Ascension Technology Corporation DC magnetic-based position and orientation monitoring system for tracking medical instruments
US9173617B2 (en) * 2011-10-19 2015-11-03 Mayo Foundation For Medical Education And Research Method for controlling radiation dose and intravenous contrast dose in computed tomography imaging
CN104644203B (zh) * 2014-09-02 2018-01-23 沈阳东软医疗系统有限公司 一种剂量调制扫描方法和装置
CN104504655A (zh) * 2014-12-04 2015-04-08 沈阳东软医疗系统有限公司 一种ct金属伪影处理方法和装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102567958A (zh) * 2010-12-31 2012-07-11 东软飞利浦医疗设备系统有限责任公司 一种去除ct图像中金属伪影的图像后处理方法
CN103971387A (zh) * 2013-01-29 2014-08-06 上海联影医疗科技有限公司 Ct图像重建方法
US9349198B2 (en) * 2013-07-26 2016-05-24 General Electric Company Robust artifact reduction in image reconstruction
WO2015168147A1 (en) * 2014-04-29 2015-11-05 Carl Zeiss X-ray Microscopy, Inc. Segmentation and spectrum based metal artifact reduction method and system
CN105225208A (zh) * 2015-10-14 2016-01-06 上海联影医疗科技有限公司 一种计算机断层成像金属伪影校正方法及装置
CN105701778A (zh) * 2016-01-11 2016-06-22 赛诺威盛科技(北京)有限公司 一种从ct图像中去除金属伪影的方法
CN106683143A (zh) * 2016-01-21 2017-05-17 上海联影医疗科技有限公司 图像金属伪影校正方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"CT金属伪影校正研究";李铭;《中国博士学位论文全文数据库(电子期刊)信息科技辑》;20150915;I138-63 *
"金属伪影校正在工业CT中的应用研究";陈云斌等;《中国体视学与图像分析》;20170331;第22卷(第3期);第307-313页 *

Also Published As

Publication number Publication date
CN107978003A (zh) 2018-05-01
US20190130610A1 (en) 2019-05-02
US10783676B2 (en) 2020-09-22

Similar Documents

Publication Publication Date Title
CN107978003B (zh) 一种ct图像金属伪影的处理方法及装置
CN106709917B (zh) 神经网络模型训练方法、装置及系统
US9437017B2 (en) Method and apparatus for metal artifact elimination in a medical image
JP5275668B2 (ja) 画像処理装置及び画像処理方法
JP5207649B2 (ja) 骨または他のカルシウム含有物質と軟部組織内の造影剤との自動差別化方法および装置
CN108158597B (zh) 确定原始x射线能量数据的方法、装置及ct设备
US10339655B2 (en) Automated image evaluation in x-ray imaging
US7103135B2 (en) Method and apparatus for generating an improved image of natural tissue in reconstructing body images from 3D-measurements
EP2891455B1 (en) Radiation tomographic image-generating apparatus and radiation tomographic image-generating method
US7602951B2 (en) Method and system for providing defective cell correction in a medical imaging device
JP2005296605A (ja) 放射線写真画像を診断関連領域と診断非関連領域とにセグメント化する方法
CN110610498A (zh) 乳腺钼靶图像处理方法、系统、存储介质及设备
US10950015B2 (en) Medical image processing method, and computer readable storage medium
EP2989983B1 (en) Image processing device
US8503746B2 (en) Medical image processing device, medical image processing method, and medical image processing program
KR101467380B1 (ko) 의료 영상 화질 개선 방법 및 그 장치
Rovaris et al. Segmentation methods for micro CT images: a comparative study using human bone samples
CN110866898A (zh) 一种乳腺图像处理方法、设备及存储介质
US10872445B2 (en) Apparatus for tomosynthesis image reconstruction
CN107330895B (zh) 医学图像的分割方法及装置
CN111803107B (zh) 一种金属检测方法、装置、设备及存储介质
CN110811662A (zh) 一种扫描剂量调制的方法、装置、设备及存储介质
KR102477991B1 (ko) 의료 영상 처리 방법 및 장치
JP2006263055A (ja) X線画像処理システム及びx線画像処理方法
CN112562030A (zh) 图像重建方法、装置和电子设备

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
CB02 Change of applicant information

Address after: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province

Applicant after: DongSoft Medical System Co., Ltd.

Address before: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province

Applicant before: Dongruan Medical Systems Co., Ltd., Shenyang

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant