CN111192228B - 图像处理方法、装置、ct设备及ct*** - Google Patents
图像处理方法、装置、ct设备及ct*** Download PDFInfo
- Publication number
- CN111192228B CN111192228B CN202010003011.0A CN202010003011A CN111192228B CN 111192228 B CN111192228 B CN 111192228B CN 202010003011 A CN202010003011 A CN 202010003011A CN 111192228 B CN111192228 B CN 111192228B
- Authority
- CN
- China
- Prior art keywords
- target
- image
- reconstructed image
- original reconstructed
- value
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 17
- 210000004872 soft tissue Anatomy 0.000 claims abstract description 45
- 238000002591 computed tomography Methods 0.000 claims description 99
- 238000004088 simulation Methods 0.000 claims description 42
- 238000000034 method Methods 0.000 claims description 29
- 230000008569 process Effects 0.000 claims description 14
- 238000010586 diagram Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 210000002216 heart Anatomy 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 210000003205 muscle Anatomy 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/77—Retouching; Inpainting; Scratch removal
-
- 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]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明实施例提供一种图像处理方法、装置、CT设备及CT***。本发明实施例通过获取待处理的目标原始重建图像,目标原始重建图像为含有锥角伪影的CT图像,将目标原始重建图像中空气区域的像素值全部设置为第一值,将原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像,根据目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像,从目标原始重建图像中减去所述目标误差图像,得到目标图像,基于图像域获取能够更准确地代表锥角伪影的误差图像,误差小,从而能够更准确地去除目标原始重建图像中的锥角伪影,提高图像质量。
Description
技术领域
本发明涉及医学图像处理技术领域,尤其涉及一种图像处理方法、装置、CT设备及CT***。
背景技术
CT(Computed Tomography,电子计算机断层扫描)成像是一种安全、准确的无创影像学技术,在医学领域应用广泛。然而,在某些扫描类型的CT重建图像中,存在比较多的锥角伪影,例如圆周锥束CT扫描、螺旋锥束CT扫描等。在这些CT扫描中,锥角越大,锥角伪影越严重。锥角伪影的存在严重影响了图像质量,给医生使用CT图像辅助诊断带来很大干扰。
发明内容
为克服相关技术中存在的问题,本发明提供了一种图像处理方法、装置、CT设备及CT***,减小图像中的锥角伪影,提高图像质量。
根据本发明实施例的第一方面,提供一种图像处理方法,包括:
获取待处理的目标原始重建图像;所述目标原始重建图像为含有锥角伪影的CT图像;
将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像;
根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像;
从所述目标原始重建图像中减去所述目标误差图像,得到目标图像。
根据本发明实施例的第二方面,提供一种图像重建装置,包括:
获取模块,用于获取待处理的目标原始重建图像;所述目标原始重建图像为含有锥角伪影的CT图像;
设置模块,用于将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像;
确定模块,用于根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像;
相减模块,用于从所述目标原始重建图像中减去所述目标误差图像,得到目标图像。
根据本发明实施例的第三方面,提供一种CT设备,包括:内部总线,以及通过内部总线连接的存储器、处理器和外部接口;其中,所述外部接口,用于连接CT***的探测器,所述探测器包括多个探测器室及相应的处理电路;
所述存储器,用于存储图像处理逻辑对应的机器可读指令;
所述处理器,用于读取所述存储器上的所述机器可读指令,并执行如下操作:
获取待处理的目标原始重建图像;所述目标原始重建图像为含有锥角伪影的CT图像;
将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像;
根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像;
从所述目标原始重建图像中减去所述目标误差图像,得到目标图像。
根据本发明实施例的第四方面,提供一种CT***,包括探测器、扫描床和CT设备,所述探测器包括多个探测器室及相应的处理电路;其中:
所述探测器室,用于在所述CT***扫描过程中,探测穿过扫描对象的X射线并转换为电信号;
所述处理电路,用于将所述电信号转换成脉冲信号,采集脉冲信号的能量信息;
所述CT设备,用于:
获取待处理的目标原始重建图像;所述目标原始重建图像为含有锥角伪影的CT图像;
将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像;
根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像;
从所述目标原始重建图像中减去所述目标误差图像,得到目标图像。
本发明实施例提供的技术方案可以包括以下有益效果:
本发明实施例,通过获取待处理的目标原始重建图像,目标原始重建图像为含有锥角伪影的CT图像,将目标原始重建图像中空气区域的像素值全部设置为第一值,将原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像,根据目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像,从目标原始重建图像中减去所述目标误差图像,得到目标图像,基于图像域获取能够更准确地代表锥角伪影的误差图像,误差小,从而能够更准确地去除目标原始重建图像中的锥角伪影,提高图像质量。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本说明书。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本说明书的实施例,并与说明书一起用于解释本说明书的原理。
图1是本发明实施例提供的图像处理方法的一个流程示例图。
图2是本发明实施例提供的图像处理装置的功能方块图。
图3是本发明实施例提供的CT设备的一个硬件结构图。
图4是本发明实施例提供的图像处理方法的又一个流程示例图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本发明实施例的一些方面相一致的装置和方法的例子。
在本发明实施例使用的术语是仅仅出于描述特定本发明实施例的目的,而非旨在限制本发明实施例。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。还应当理解,本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解,尽管在本发明实施例可能采用术语第一、第二、第三等来描述各种信息,但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如,在不脱离本发明实施例范围的情况下,第一信息也可以被称为第二信息,类似地,第二信息也可以被称为第一信息。取决于语境,如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
在螺旋锥束CT扫描、圆周锥束CT扫描等CT扫描中,使用解析类重建算法,例如近似解析重建算法,对CT扫描生数据进行图像重建会造成严重的伪影。并且,随着锥角的增大,伪影会变得更加严重,因此这类伪影被称为锥角伪影。
相关技术中,基于正弦域生数据处理来抑制或者减少锥角伪影。该方法的原理是:对于使用解析类重建算法重建得到的原始重建图像进行前向正投影(Forward Project),得到正弦域模拟投影;用正弦域模拟投影减去正弦域生数据投影得到正弦域误差投影;对正弦域误差投影进行重建得到误差图像;用原始重建图像减去误差图像得到修正后的锥角伪影减少的图像。
这种技术中,用于对原始重建图像进行修正的误差图像是对正弦域中模拟投影和生数据投影的偏差进行重建得到的。如果采集的生数据投影质量较差或者数量较少,则会引入较大误差。因此这种技术去除锥角伪影的效果较差,图像质量较低。
本发明实施例提供一种图像处理方法,基于图像域去除锥角伪影,能够获得较好的去除锥角伪影的效果,提高图像质量。
本发明实施例提供的图像处理方法可以适用于所有包含锥角伪影的CT图像。
下面通过实施例对本发明的图像处理方法进行详细说明。
图1是本发明实施例提供的图像处理方法的一个流程示例图。如图1所示,本实施例中,图像处理方法可以包括:
S101,获取待处理的目标原始重建图像,目标原始重建图像为含有锥角伪影的CT图像。
S102,将目标原始重建图像中空气区域的像素值全部设置为第一值,将原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像。
S103,根据目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像。
S104,从目标原始重建图像中减去所述目标误差图像,得到目标图像。
本实施例中,目标原始重建图像可以是对CT扫描生数据利用解析类重建算法重建得到的图像。例如,目标原始重建图像可以是对螺旋锥束CT扫描采集的生数据使用近似解析重建算法得到的CT重建图像。
需要说明的是,本实施例中的目标原始重建图像不限于对螺旋锥束CT扫描采集的生数据使用解析类重建算法得到的CT重建图像、对圆周锥束CT扫描采集的生数据使用解析类重建算法得到的CT重建图像,还可以是其他的含有锥角误差的CT图像。
本实施例中,目标原始重建图像可以是对人体进行CT扫描后重建得到的图像。
在对人体进行CT扫描后得到的CT图像中,通常包括三个部分,即:空气区域、软组织区域、骨区域。其中,软组织可以是心脏、肝脏、肌肉等。
每个区域的CT值都在一定的范围之内。因此,在一个示例中,可以根据CT值判断CT图像中的任一个像素点属于哪个区域。
锥角伪影通常存在于CT图像的空气区域和软组织区域,而骨区域不存在伪影。通过将目标原始重建图像中空气区域的像素值全部设置为第一值,将原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像,使得目标基图像中不包含锥角伪影。
需要说明的是,第一值和第二值是不同的。
在一个示例中,第一值为所述目标原始重建图像中空气区域的CT值取值范围内的CT值;所述第二值为所述目标原始重建图像中软组织区域的CT值取值范围内的CT值。
例如,假设目标原始重建图像中空气区域的CT值的最大值为最小值为/>则第一值可以是区间/>中的任意一个值。
同样,假设目标原始重建图像中软组织区域的CT值的最大值为最小值为则第一值可以是区间/>中的任意一个值。
在一个示例中,所述第一值可以为所述目标原始重建图像中空气区域的最大CT值和最小CT值的平均值;所述第二值可以为所述目标原始重建图像中软组织区域的最大CT值和最小CT值的平均值。
例如,假设目标原始重建图像中空气区域的CT值取值范围是目标原始重建图像中软组织区域的CT值取值范围是/>则第一值可以为/>第二值可以为/>
本实施例中,根据目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像,由于目标基图像中不含有锥角伪影,因此根据目标基图像确定的目标误差图像能够得到更加准确地表示锥角伪影,从而有助于准确地去除目标原始重建图像中的锥角伪影。
本实施例中,通过将目标原始重建图像中与目标误差图像作差得到目标图像,由于目标误差图像能够得到更加准确地表示锥角伪影,因此目标图像中的锥角伪影更少,去除效果更好。
在一个示例性的实现过程中,在步骤S102之前,还可以包括:
确定原始重建图像中的空气区域和软组织区域。
例如,可以根据空气的CT值范围判断图像中的任一个像素点是否属于空气区域。同理,可以根据软组织的CT值范围判断图像中的任一个像素点是否属于软组织区域。
图4是本发明实施例提供的图像处理方法的又一个流程示例图。图4示出了根据目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像的过程。如图4所示,在一个示例性的实现过程中,步骤S103可以包括:S1031,对所述目标基图像进行前向投影,得到目标模拟投影;
S1032,对所述目标模拟投影进行图像重建,得到目标模拟重建图像;
S1033,从所述目标模拟重建图像中减去所述目标基图像,得到所述目标原始重建图像中的锥角伪影对应的目标误差图像。
本实施例中,目标模拟投影为正弦域模拟投影。
对所述目标基图像进行前向投影,得到目标模拟投影,可以包括:
将所述目标基图像与预设的前向投影算子相乘,得到目标模拟投影。
其中,前向投影算子可以是一个函数。
其中,对所述目标模拟投影进行图像重建采用的重建算法为解析类重建算法,例如近似解析重建算法。目标模拟重建图像中含有因解析类重建导致的锥角伪影。
由于目标模拟重建图像中含有锥角伪影,而目标基图像中不含有锥角伪影,因此二者作差能够得到比较准确地表示锥角伪影的误差图像。
下面通过一个示例对本发明实施例的图像处理方法作进一步详细说明。
假设图像ffbp为待处理的目标原始重建图像。图像ffbp是通过对人体进行螺旋CT扫描,并使用近似解析重建算法重建得到的CT图像,图像ffbp中空气区域的CT值取值范围是目标原始重建图像中软组织区域的CT值取值范围是/>
则对图像ffbp的处理过程如下:
第一步,获取图像ffbp。
第二步,对图像ffbp进行图像分割,以确定图像ffbp中的空气区域和软组织区域。
第三步,将图像ffbp中空气区域的所有像素值设置为将图像ffbp中软组织区域的所有像素值设置为/>得到图像ffbp对应的基图像ft,图像ft中不含锥角伪影。
第四步,利用前向投影算子p,根据下面的公式(1)对图像ft进行前向投影,得到基图像ft对应的模拟投影ps。
ps=p×ft (1)
第五步,对模拟投影ps使用近似解析重建算法进行重建,得到模拟重建图像fs,图像fs含有因近似解析导致的锥角伪影。
第六步,根据如下的公式(2)计算得到图像ffbp对应的误差图像fbias。fbias刻画了锥角伪影。
fbias=fs-ft (2)
第七步,根据如下的公式(3)计算得到图像ffbp对应的去除锥角伪影的目标图像f*。
f*=ffbp-fbias (3)
如果图像f*符合要求,则停止,以图像f*作为最终的目标图像;如果图像f*不符合要求,则用图像f*代替图像ffbp,重复上述的第二步至第七步。
本发明实施例提供的图像处理方法,通过获取待处理的目标原始重建图像,目标原始重建图像为含有锥角伪影的CT图像,将目标原始重建图像中空气区域的像素值全部设置为第一值,将原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像,根据目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像,从目标原始重建图像中减去所述目标误差图像,得到目标图像,基于图像域获取能够更准确地代表锥角伪影的误差图像,误差小,从而能够更准确地去除目标原始重建图像中的锥角伪影,提高图像质量。
基于上述的方法实施例,本发明实施例还提供了相应的装置、设备及存储介质实施例。
图2是本发明实施例提供的图像重建装置的功能方块图。如图2所示,本实施例中,图像重建装置可以包括:
获取模块210,用于获取待处理的目标原始重建图像;所述目标原始重建图像为含有锥角伪影的CT图像;
设置模块220,用于将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像;
确定模块230,用于根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像;
相减模块240,用于从所述目标原始重建图像中减去所述目标误差图像,得到目标图像。
在一个示例性的实现过程中,确定模块230可以具体用于:
对所述目标基图像进行前向投影,得到目标模拟投影;
对所述目标模拟投影进行图像重建,得到目标模拟重建图像;
从所述目标基图像中减去所述目标模拟重建图像,得到所述目标原始重建图像中的锥角伪影对应的目标误差图像。
在一个示例性的实现过程中,所述第一值为所述目标原始重建图像中空气区域的CT值取值范围内的CT值;所述第二值为所述目标原始重建图像中软组织区域的CT值取值范围内的CT值。
在一个示例性的实现过程中,所述装置还可以包括:
区域确定模块,用于确定所述原始重建图像中的空气区域和软组织区域。
在一个示例性的实现过程中,所述目标原始重建图像是对CT扫描生数据利用解析类重建算法重建得到的图像。
本发明实施例还提供了一种CT设备。图3是本发明实施例提供的CT设备的一个硬件结构图。如图3所示,CT设备包括:内部总线301,以及通过内部总线连接的存储器302,处理器303和外部接口304,其中,所述外部接口,用于连接CT***的探测器,所述探测器包括多个探测器室及相应的处理电路;
所述存储器302,用于存储图像处理逻辑对应的机器可读指令;
所述处理器303,用于读取存储器302上的机器可读指令,并执行所述指令以实现如下操作:
获取待处理的目标原始重建图像;所述目标原始重建图像为含有锥角伪影的CT图像;
将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像;
根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像;
从所述目标原始重建图像中减去所述目标误差图像,得到目标图像。
在一个示例性的实现过程中,根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像,包括:
对所述目标基图像进行前向投影,得到目标模拟投影;
对所述目标模拟投影进行图像重建,得到目标模拟重建图像;
从所述目标模拟重建图像中减去所述目标基图像,得到所述目标原始重建图像中的锥角伪影对应的目标误差图像。
在一个示例性的实现过程中,所述第一值为所述目标原始重建图像中空气区域的CT值取值范围内的CT值;所述第二值为所述目标原始重建图像中软组织区域的CT值取值范围内的CT值。
在一个示例性的实现过程中,将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像之前,还包括:
确定所述原始重建图像中的空气区域和软组织区域。
在一个示例性的实现过程中,所述目标原始重建图像是对CT扫描生数据利用解析类重建算法重建得到的图像。
本发明实施例还提供一种CT***,包括探测器、扫描床和CT设备,所述探测器包括多个探测器室及相应的处理电路;其中:
所述探测器室,用于在所述CT***扫描过程中,探测穿过扫描对象的X射线并转换为电信号;
所述处理电路,用于将所述电信号转换成脉冲信号,采集脉冲信号的能量信息;
所述CT设备,用于:
获取待处理的目标原始重建图像;所述目标原始重建图像为含有锥角伪影的CT图像;
将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像;
根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像;
从所述目标原始重建图像中减去所述目标误差图像,得到目标图像。
在一个示例性的实现过程中,根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像,包括:
对所述目标基图像进行前向投影,得到目标模拟投影;
对所述目标模拟投影进行图像重建,得到目标模拟重建图像;
从所述目标模拟重建图像中减去所述目标基图像,得到所述目标原始重建图像中的锥角伪影对应的目标误差图像。
在一个示例性的实现过程中,所述第一值为所述目标原始重建图像中空气区域的CT值取值范围内的CT值;所述第二值为所述目标原始重建图像中软组织区域的CT值取值范围内的CT值。
在一个示例性的实现过程中,将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像之前,还包括:
确定所述原始重建图像中的空气区域和软组织区域。
在一个示例性的实现过程中,所述目标原始重建图像是对CT扫描生数据利用解析类重建算法重建得到的图像。
本发明实施例还提供一种计算机可读存储介质,其上存储有计算机程序,其中,所述程序被处理器执行时实现如下操作:
获取待处理的目标原始重建图像;所述目标原始重建图像为含有锥角伪影的CT图像;
将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像;
根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像;
从所述目标原始重建图像中减去所述目标误差图像,得到目标图像。
在一个示例性的实现过程中,根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像,包括:
对所述目标基图像进行前向投影,得到目标模拟投影;
对所述目标模拟投影进行图像重建,得到目标模拟重建图像;
从所述目标模拟重建图像中减去所述目标基图像,得到所述目标原始重建图像中的锥角伪影对应的目标误差图像。
在一个示例性的实现过程中,所述第一值为所述目标原始重建图像中空气区域的CT值取值范围内的CT值;所述第二值为所述目标原始重建图像中软组织区域的CT值取值范围内的CT值。
在一个示例性的实现过程中,将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像之前,还包括:
确定所述原始重建图像中的空气区域和软组织区域。
在一个示例性的实现过程中,所述目标原始重建图像是对CT扫描生数据利用解析类重建算法重建得到的图像。
对于装置和设备实施例而言,由于其基本对应于方法实施例,所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的模块可以是或者也可以不是物理上分开的,作为模块显示的部件可以是或者也可以不是物理模块,即可以位于一个地方,或者也可以分布到多个网络模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本说明书方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
本领域技术人员在考虑说明书及实践这里申请的发明后,将容易想到本说明书的其它实施方案。本说明书旨在涵盖本说明书的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本说明书的一般性原理并包括本说明书未申请的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本说明书的真正范围和精神由下面的权利要求指出。
应当理解的是,本说明书并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。本说明书的范围仅由所附的权利要求来限制。
以上所述仅为本说明书的较佳实施例而已,并不用以限制本说明书,凡在本说明书的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本说明书保护的范围之内。
Claims (8)
1.一种图像处理方法,其特征在于,包括:
获取待处理的目标原始重建图像;所述目标原始重建图像为含有锥角伪影的CT图像;
将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像;
根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像,包括:对所述目标基图像进行前向投影,得到目标模拟投影;对所述目标模拟投影进行图像重建,得到目标模拟重建图像;从所述目标模拟重建图像中减去所述目标基图像,得到所述目标原始重建图像中的锥角伪影对应的目标误差图像;
从所述目标原始重建图像中减去所述目标误差图像,得到目标图像。
2.根据权利要求1所述的方法,其特征在于,所述第一值为所述目标原始重建图像中空气区域的CT值取值范围内的CT值;所述第二值为所述目标原始重建图像中软组织区域的CT值取值范围内的CT值。
3.根据权利要求2所述的方法,其特征在于,将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像之前,还包括:
确定所述原始重建图像中的空气区域和软组织区域。
4.根据权利要求1所述的方法,其特征在于,所述目标原始重建图像是对CT扫描生数据利用解析类重建算法重建得到的图像。
5.一种图像处理装置,其特征在于,包括:
获取模块,用于获取待处理的目标原始重建图像;所述目标原始重建图像为含有锥角伪影的CT图像;
设置模块,用于将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像;
确定模块,用于根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像,具体用于:对所述目标基图像进行前向投影,得到目标模拟投影;对所述目标模拟投影进行图像重建,得到目标模拟重建图像;从所述目标基图像中减去所述目标模拟重建图像,得到所述目标原始重建图像中的锥角伪影对应的目标误差图像;
相减模块,用于从所述目标原始重建图像中减去所述目标误差图像,得到目标图像。
6.根据权利要求5所述的装置,其特征在于,所述第一值为所述目标原始重建图像中空气区域的CT值取值范围内的CT值;所述第二值为所述目标原始重建图像中软组织区域的CT值取值范围内的CT值。
7.一种CT设备,其特征在于,包括:内部总线,以及通过内部总线连接的存储器、处理器和外部接口;其中,所述外部接口,用于连接CT***的探测器,所述探测器包括多个探测器室及相应的处理电路;
所述存储器,用于存储图像处理逻辑对应的机器可读指令;
所述处理器,用于读取所述存储器上的所述机器可读指令,并执行如下操作:
获取待处理的目标原始重建图像;所述目标原始重建图像为含有锥角伪影的CT图像;
将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像;
根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像,包括:对所述目标基图像进行前向投影,得到目标模拟投影;对所述目标模拟投影进行图像重建,得到目标模拟重建图像;从所述目标模拟重建图像中减去所述目标基图像,得到所述目标原始重建图像中的锥角伪影对应的目标误差图像;
从所述目标原始重建图像中减去所述目标误差图像,得到目标图像。
8.一种CT***,其特征在于,包括探测器、扫描床和CT设备,所述探测器包括多个探测器室及相应的处理电路;其中:
所述探测器室,用于在所述CT***扫描过程中,探测穿过扫描对象的X射线并转换为电信号;
所述处理电路,用于将所述电信号转换成脉冲信号,采集脉冲信号的能量信息;
所述CT设备,用于:
获取待处理的目标原始重建图像;所述目标原始重建图像为含有锥角伪影的CT图像;
将所述目标原始重建图像中空气区域的像素值全部设置为第一值,将所述原始重建图像中软组织区域的像素值全部设置为第二值,得到目标基图像;
根据所述目标基图像确定所述目标原始重建图像中的锥角伪影对应的目标误差图像,包括:对所述目标基图像进行前向投影,得到目标模拟投影;对所述目标模拟投影进行图像重建,得到目标模拟重建图像;从所述目标模拟重建图像中减去所述目标基图像,得到所述目标原始重建图像中的锥角伪影对应的目标误差图像;
从所述目标原始重建图像中减去所述目标误差图像,得到目标图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010003011.0A CN111192228B (zh) | 2020-01-02 | 2020-01-02 | 图像处理方法、装置、ct设备及ct*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010003011.0A CN111192228B (zh) | 2020-01-02 | 2020-01-02 | 图像处理方法、装置、ct设备及ct*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111192228A CN111192228A (zh) | 2020-05-22 |
CN111192228B true CN111192228B (zh) | 2024-03-08 |
Family
ID=70708094
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010003011.0A Active CN111192228B (zh) | 2020-01-02 | 2020-01-02 | 图像处理方法、装置、ct设备及ct*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111192228B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111798534B (zh) * | 2020-07-17 | 2024-03-08 | 东软医疗***股份有限公司 | 图像重建方法、装置、控制台设备及ct*** |
CN112184850B (zh) * | 2020-09-30 | 2024-01-05 | 沈阳先进医疗设备技术孵化中心有限公司 | 图像处理方法、装置、控制台设备及ct*** |
CN112244884B (zh) * | 2020-10-27 | 2023-08-29 | 沈阳先进医疗设备技术孵化中心有限公司 | 骨图像获取方法、装置、控制台设备及ct*** |
CN112258596A (zh) * | 2020-10-27 | 2021-01-22 | 沈阳先进医疗设备技术孵化中心有限公司 | 图像生成方法、装置、控制台设备及ct*** |
CN112562030A (zh) * | 2020-12-08 | 2021-03-26 | 赛诺威盛科技(北京)有限公司 | 图像重建方法、装置和电子设备 |
CN113240762A (zh) * | 2021-04-06 | 2021-08-10 | 沈阳先进医疗设备技术孵化中心有限公司 | 图像处理方法、装置及电子设备 |
CN113284205B (zh) * | 2021-04-23 | 2024-01-02 | 沈阳先进医疗设备技术孵化中心有限公司 | Ct迭代重建方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103186883A (zh) * | 2011-12-29 | 2013-07-03 | 上海联影医疗科技有限公司 | 一种ct图像重建中骨硬化伪影的校正方法 |
WO2018126434A1 (zh) * | 2017-01-06 | 2018-07-12 | 深圳先进技术研究院 | Ct图像阴影校正方法、装置及电子设备 |
CN110097517A (zh) * | 2019-04-28 | 2019-08-06 | 东软医疗***股份有限公司 | 去除图像伪影的方法及装置 |
CN110335325A (zh) * | 2019-06-27 | 2019-10-15 | 深圳安科高技术股份有限公司 | 一种ct图像重建方法及其*** |
CN110400359A (zh) * | 2019-07-23 | 2019-11-01 | 上海联影医疗科技有限公司 | 图像伪影的消除方法、装置、设备及存储介质 |
-
2020
- 2020-01-02 CN CN202010003011.0A patent/CN111192228B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103186883A (zh) * | 2011-12-29 | 2013-07-03 | 上海联影医疗科技有限公司 | 一种ct图像重建中骨硬化伪影的校正方法 |
WO2018126434A1 (zh) * | 2017-01-06 | 2018-07-12 | 深圳先进技术研究院 | Ct图像阴影校正方法、装置及电子设备 |
CN110097517A (zh) * | 2019-04-28 | 2019-08-06 | 东软医疗***股份有限公司 | 去除图像伪影的方法及装置 |
CN110335325A (zh) * | 2019-06-27 | 2019-10-15 | 深圳安科高技术股份有限公司 | 一种ct图像重建方法及其*** |
CN110400359A (zh) * | 2019-07-23 | 2019-11-01 | 上海联影医疗科技有限公司 | 图像伪影的消除方法、装置、设备及存储介质 |
Non-Patent Citations (1)
Title |
---|
何畅 ; 王剑钢 ; 杨旭 ; .锥形束CT图像去噪与重建算法的研究进展.中国医疗设备.2013,(06),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN111192228A (zh) | 2020-05-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111192228B (zh) | 图像处理方法、装置、ct设备及ct*** | |
CN111915696B (zh) | 三维图像数据辅助的低剂量扫描数据重建方法及电子介质 | |
CN110751702B (zh) | 图像重建方法、***、装置及存储介质 | |
CN106530236B (zh) | 一种医学图像处理方法及*** | |
CN107115119B (zh) | Pet图像衰减系数的获取方法、衰减校正的方法及*** | |
CN108876730B (zh) | 校正运动伪影的方法、装置及设备和存储介质 | |
US8658979B2 (en) | Nuclear image reconstruction | |
CN112017131B (zh) | Ct图像金属伪影去除方法、装置及计算机可读存储介质 | |
CN111627083B (zh) | 骨硬化伪影校正方法、装置、计算机设备和可读存储介质 | |
CN113706409A (zh) | 基于人工智能的cbct图像增强方法、装置和存储介质 | |
CN114972115A (zh) | 一种图像降噪方法、装置、设备及介质 | |
CN111612867B (zh) | 运动伪影校正方法、装置、计算机设备和可读存储介质 | |
CN111311531B (zh) | 图像增强方法、装置、控制台设备及医学成像*** | |
JP7238134B2 (ja) | Pet画像化中の自動的な動き補正 | |
CN107569251B (zh) | 医学成像方法和***及非暂态计算机可读存储介质 | |
CN113313649B (zh) | 图像重建方法及装置 | |
CN112258596A (zh) | 图像生成方法、装置、控制台设备及ct*** | |
CN112244884B (zh) | 骨图像获取方法、装置、控制台设备及ct*** | |
CN110327066B (zh) | 心脏运动信号获取方法、装置、计算机设备和存储介质 | |
CN110706338B (zh) | 图像重建方法、装置、ct设备及ct*** | |
CN113192155A (zh) | 螺旋ct锥束扫描图像重建方法、扫描***及存储介质 | |
CN111754436A (zh) | 医学图像伪影校正的加速方法、计算机设备和存储介质 | |
CN111260575B (zh) | 数据处理方法、装置、ct设备及ct*** | |
CN111932475B (zh) | 滤波方法、装置、ct设备及ct*** | |
CN111915516B (zh) | 运动补偿方法、装置、ct设备及ct*** |
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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20240204 Address after: 110167 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Applicant after: Shenyang Neusoft Medical Systems Co.,Ltd. Country or region after: China Address before: Room 336, 177-1, Chuangxin Road, Hunnan New District, Shenyang City, Liaoning Province Applicant before: Shenyang advanced medical equipment Technology Incubation Center Co.,Ltd. Country or region before: China |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |