CN111000578A - 图像重建方法、装置、ct设备及ct系统 - Google Patents
图像重建方法、装置、ct设备及ct系统 Download PDFInfo
- Publication number
- CN111000578A CN111000578A CN201911359736.7A CN201911359736A CN111000578A CN 111000578 A CN111000578 A CN 111000578A CN 201911359736 A CN201911359736 A CN 201911359736A CN 111000578 A CN111000578 A CN 111000578A
- Authority
- CN
- China
- Prior art keywords
- sampling
- data
- axial
- target
- raw data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000002591 computed tomography Methods 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000005070 sampling Methods 0.000 claims abstract description 336
- 230000006870 function Effects 0.000 claims description 25
- 238000013519 translation Methods 0.000 claims description 24
- 238000012545 processing Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 10
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- 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]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (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系统扫描过程中,探测穿过扫描对象的X射线并转换为电信号;
所述处理电路,用于将所述电信号转换成脉冲信号,采集脉冲信号的能量信息;
所述CT设备,用于:
获取对受检对象进行多排CT扫描得到的原始轴向采样生数据;
对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据;
根据所述目标轴向采样生数据进行图像重建,得到目标重建图像。
本发明实施例提供的技术方案可以包括以下有益效果:
本发明实施例,通过获取对受检对象进行多排CT扫描得到的原始轴向采样生数据,对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据,根据所述目标轴向采样生数据进行图像重建,得到目标重建图像,利用对原始轴向采样生数据的上采样,能够获得数倍于原始轴向采样生数据的轴向采样生数据,以软件方式提高了轴向采样率,去除了风车伪影,提高了重建图像的图像质量,并且不需要使用具有支持高轴向采样率的硬件的价格高昂的设备,降低了使用成本。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本说明书。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本说明书的实施例,并与说明书一起用于解释本说明书的原理。
图1是具有风车伪影的CT图像的示例图。
图2是本发明实施例提供的图像重建方法的流程示例图。
图3是本发明实施例提供的上采样示意图。
图4是本发明实施例提供的图像重建装置的功能方块图。
图5是本发明实施例提供的CT设备的一个硬件结构图。
图6是根据本发明实施例提供的图像重建方法获得的CT图像的示例图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本发明实施例的一些方面相一致的装置和方法的例子。
在本发明实施例使用的术语是仅仅出于描述特定本发明实施例的目的,而非旨在限制本发明实施例。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。还应当理解,本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解,尽管在本发明实施例可能采用术语第一、第二、第三等来描述各种信息,但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如,在不脱离本发明实施例范围的情况下,第一信息也可以被称为第二信息,类似地,第二信息也可以被称为第一信息。取决于语境,如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
风车伪影的形成原因是轴向采样率不足,导致重建得到的CT交替出现白黑条,严重影响了图像质量。图1是具有风车伪影的CT图像的示例图。如图1所示,在图像区域存在非常明显的明暗交替的条纹,这些条纹即为风车伪影。
目前去除风车伪影的方法主要分为两种,一种方法是从硬件出发来提高硬件的轴向采样率从而去除风车伪影,如Z飞焦点技术,这种方法需要硬件的支持,设备价格高,因此使用成本较高;另一种方法是从软件出发,通过算法来去除风车伪影。目前的主要算法为对图像进行后处理来去除风车伪影,这种图像后处理算法复杂度较高,运行效率较低。
本发明实施例提供的图像重建方法是通过软件算法来去除风车伪影的方法。本发明实施例通过对多排CT扫描得到的原始轴向采样生数据进行上采样,能够获得数倍于原始轴向采样生数据的轴向采样生数据,以软件方式提高了轴向采样率,从而去除了风车伪影,提高了重建图像的图像质量,并且算法复杂度低,运行效率高。
下面通过实施例对本发明的图像重建方法进行详细说明。
图2是本发明实施例提供的图像重建方法的流程示例图。如图2所示,本实施例中,图像重建方法可以包括:
S201,获取对受检对象进行多排CT扫描得到的原始轴向采样生数据。
S202,对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据。
S203,根据所述目标轴向采样生数据进行图像重建,得到目标重建图像。
本实施例中,对受检对象进行多排CT扫描采集到的生数据是在轴向采样的生数据。原始轴向采样生数据中的每一组生数据包括采样值和采样位置,其中采样位置为采样点的轴向位置。
假设原始轴向采样生数据的轴向采样数(也即物理切片数)为N,轴向起始采样位置为z0,轴向采样间隔为Δz,螺距为p,物理切片厚度为ΔT。则原始轴向采样生数据包括采样值P0、P1……PN-1,以及采样值P0、P1……PN-1对应的采样位置z0、z0+Δz……z0+(N-1)Δz。
本实施例中,上采样是指将原始轴向采样生数据的采样数量翻倍。例如假设原始轴向采样生数据对应的数据矩阵大小为672*32,则可以通过上采样的方式将该数据矩阵扩大到672*64、672*128等等,如此,即可将轴向采样生数据的数据量扩大至原始轴向采样生数据的数据量的数倍。这相当于提高了轴向采样率,因此能够去除风车伪影,提高重建图像的图像质量。
本实施例在不需要设备具有支持高轴向采样率的硬件,因此使用成本低。
并且,本实施例算法复杂度低,运行效率高。
步骤S203中,由于目标轴向采样生数据数倍于原始轴向采样生数据,因此根据目标轴向采样生数据数重建得到的重建图像,相比于根据原始轴向采样生数据重建得到的重建图像具有更少的风车伪影,提高了图像质量。
在一个示例性的实现过程中,步骤S202可以包括:
根据上采样倍数M,确定M个平移偏移量,M为自然数,且M≥2;
分别按照所述M个平移偏移量中的每个平移偏移量,将所述原始轴向采样生数据在轴向进行平移,得到M组中间采样数据,所有中间采样中数据中任意两个相邻采样数据对应的轴向位置之差相等;
根据所述M组中间采样数据,得到目标轴向采样生数据。
举例说明。以前述的原始轴向采样生数据(采样值P0、P1……PN-1,采样值P0、P1……PN-1对应的采样位置z0、z0+Δz……z0+(N-1)Δz)作为输入生数据InSample。请参见图3。图3是本发明实施例提供的上采样示意图。
将生数据InSample向左平移δ,得到第一组中间采样数据InnerSample1,将生数据InSample向右平移δ,得到第二组中间采样数据InnerSample2。
InnerSample1中,采样值保持不变,仍然为P0、P1……PN-1,但采样值P0、P1……PN-1对应的轴向位置(即采样位置)变为z0-δ、z0+Δz-δ……z0+(N-1)Δz-δ,为了方便计算,将这些采样位置分别记为
InnerSample2中,采样值保持不变,仍然为P0、P1……PN-1,但采样值P0、P1……PN-1对应的轴向位置(即采样位置)变为z0+δ、z0+Δz+δ……z0+(N-1)Δz+δ,为了方便计算,将这些采样位置分别记为
在一个示例性的实现过程中,根据所述M组中间采样数据,得到目标轴向采样生数据,可以包括:
对于所述M组中间采样数据中的每一中间采样数据,获取该中间采样数据在所在中间采样数据组的指定邻域内的采样数据,作为基础采样数据;
根据该中间采样数据的轴向位置和各个基础采样数据的轴向位置,确定所述各个基础采样数据的权重值;
根据所述各个基础采样数据的权重值和所述各个基础采样数据的采样值,获得该中间采样数据对应的目标轴向采样生数据的目标采样值;
将该中间采样数据的轴向位置确定为该中间采样数据对应的目标轴向采样生数据的目标采样位置。
在一个示例性的实现过程中,根据该中间采样数据的轴向位置和各个基础采样数据的轴向位置,确定所述各个基础采样数据的权重值,包括:
对于每个基础采样数据,获取该中间采样数据的轴向位置和该基础采样数据的轴向位置的差值;
将所述差值作为预设权重函数的自变量的值,计算所述预设权重函数的函数值,以所述函数值作为该基础采样数据的权重值。
举例说明。请参见图3。输出生数据OutSample中的采样值可以根据如下的公式(1)计算得到:
其中,其中n=0,1,2……N-1,c为权值归一化参数,W为权值函数,W可以根据如下的公式(2)计算得到:
W(x)=cos(πx) (2)
c可以根据如下的公式(3)计算得到:
W-1、W0、W1根据上述的公式(2)计算得到。
请参见图3,对于InnerSample1中的中间采样数据P1(采样位置)来说,基础采样数据为P0(采样位置)、P2(采样位置),中间采样数据P1对应的目标轴向采样生数据为图3中的则可以根据基础采样数据为P0(采样位置)、P1(采样位置)、P2(采样位置)利用公式(1)计算得到采样值
在一个示例性的实现过程中,根据上采样倍数M,确定M个平移偏移量之前,可以包括:
获取上采样倍数M。
在一个示例性的实现过程中,获取上采样倍数M,可以包括:
读取系统中的默认上采样倍数值,以所述默认上采样倍数值作为上采样倍数M。
本实施例中,在系统中设置默认上采样倍数值,图像重建过程中直接根据系统中的默认上采样倍数值对原始轴向采样生数据进行上采样,速度快,处理效率高。
其中,默认上采样倍数值可以由开发人员根据经验值设定。例如假设2倍采样就可以达到符合要求的去除风车伪影的效果,则可以将默认上采样倍数值设置为2。
在一个示例性的实现过程中,获取上采样倍数M,可以包括:
接收用户输入的上采样倍数值,以所述输入的上采样倍数值作为上采样倍数M。
本实施例中,由用户来输入上采样倍数M,可以灵活性地适应不同用户的不同要求。例如,用户对于去除风车伪影的效果要求较低时,可以输入较小的值作为上采样倍数M;当用户对于去除风车伪影的效果要求较高时,可以输入较大的值作为上采样倍数M。
图6是根据本发明实施例提供的图像重建方法获得的CT图像的示例图。对比图1和图6,可知本发明实施例提供的图像重建方法能够有效去除风车伪影,提高图像质量。
本发明实施例提供的图像重建方法,通过获取对受检对象进行多排CT扫描得到的原始轴向采样生数据,对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据,根据所述目标轴向采样生数据进行图像重建,得到目标重建图像,利用对原始轴向采样生数据的上采样,能够获得数倍于原始轴向采样生数据的轴向采样生数据,以软件方式提高了轴向采样率,去除了风车伪影,提高了重建图像的图像质量,并且不需要使用具有支持高轴向采样率的硬件的价格高昂的设备,降低了使用成本。
基于上述的方法实施例,本发明实施例还提供了相应的装置、设备及存储介质实施例。
图4是本发明实施例提供的图像重建装置的功能方块图。如图4所示,本实施例中,图像重建装置可以包括:
获取模块410,用于获取对受检对象进行多排CT扫描得到的原始轴向采样生数据;
上采样模块420,用于对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据;
重建模块430,用于根据所述目标轴向采样生数据进行图像重建,得到目标重建图像。
在一个示例性的实现过程中,上采样模块420可以具体用于:
根据上采样倍数M,确定M个平移偏移量,M为自然数,且M≥2;
分别按照所述M个平移偏移量中的每个平移偏移量,将所述原始轴向采样生数据在轴向进行平移,得到M组中间采样数据,所有中间采样中数据中任意两个相邻采样数据对应的轴向位置之差相等;
根据所述M组中间采样数据,得到目标轴向采样生数据。
在一个示例性的实现过程中,根据所述M组中间采样数据,得到目标轴向采样生数据,包括:
对于所述M组中间采样数据中的每一中间采样数据,获取该中间采样数据在所在中间采样数据组的指定邻域内的采样数据,作为基础采样数据;
根据该中间采样数据的轴向位置和各个基础采样数据的轴向位置,确定所述各个基础采样数据的权重值;
根据所述各个基础采样数据的权重值和所述各个基础采样数据的采样值,获得该中间采样数据对应的目标轴向采样生数据的目标采样值;
将该中间采样数据的轴向位置确定为该中间采样数据对应的目标轴向采样生数据的目标采样位置。
在一个示例性的实现过程中,根据该中间采样数据的轴向位置和各个基础采样数据的轴向位置,确定所述各个基础采样数据的权重值,包括:
对于每个基础采样数据,获取该中间采样数据的轴向位置和该基础采样数据的轴向位置的差值;
将所述差值作为预设权重函数的自变量的值,计算所述预设权重函数的函数值,以所述函数值作为该基础采样数据的权重值。
在一个示例性的实现过程中,根据上采样倍数M,确定M个平移偏移量之前,包括:
获取上采样倍数M。
在一个示例性的实现过程中,获取上采样倍数M,包括:
读取系统中的默认上采样倍数值,以所述默认上采样倍数值作为上采样倍数M。
在一个示例性的实现过程中,获取上采样倍数M,包括:
接收用户输入的上采样倍数值,以所述输入的上采样倍数值作为上采样倍数M。
本发明实施例还提供了一种CT设备。图5是本发明实施例提供的CT设备的一个硬件结构图。如图5所示,CT设备包括:内部总线501,以及通过内部总线连接的存储器502,处理器503和外部接口504,其中,所述外部接口,用于连接CT系统的探测器,所述探测器包括多个探测器室及相应的处理电路;
所述存储器502,用于存储图像重建逻辑对应的机器可读指令;
所述处理器503,用于读取存储器502上的机器可读指令,并执行所述指令以实现如下操作:
获取对受检对象进行多排CT扫描得到的原始轴向采样生数据;
对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据;
根据所述目标轴向采样生数据进行图像重建,得到目标重建图像。
在一个示例性的实现过程中,对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据,包括:
根据上采样倍数M,确定M个平移偏移量,M为自然数,且M≥2;
分别按照所述M个平移偏移量中的每个平移偏移量,将所述原始轴向采样生数据在轴向进行平移,得到M组中间采样数据,所有中间采样中数据中任意两个相邻采样数据对应的轴向位置之差相等;
根据所述M组中间采样数据,得到目标轴向采样生数据。
在一个示例性的实现过程中,根据所述M组中间采样数据,得到目标轴向采样生数据,包括:
对于所述M组中间采样数据中的每一中间采样数据,获取该中间采样数据在所在中间采样数据组的指定邻域内的采样数据,作为基础采样数据;
根据该中间采样数据的轴向位置和各个基础采样数据的轴向位置,确定所述各个基础采样数据的权重值;
根据所述各个基础采样数据的权重值和所述各个基础采样数据的采样值,获得该中间采样数据对应的目标轴向采样生数据的目标采样值;
将该中间采样数据的轴向位置确定为该中间采样数据对应的目标轴向采样生数据的目标采样位置。
在一个示例性的实现过程中,根据该中间采样数据的轴向位置和各个基础采样数据的轴向位置,确定所述各个基础采样数据的权重值,包括:
对于每个基础采样数据,获取该中间采样数据的轴向位置和该基础采样数据的轴向位置的差值;
将所述差值作为预设权重函数的自变量的值,计算所述预设权重函数的函数值,以所述函数值作为该基础采样数据的权重值。
在一个示例性的实现过程中,根据上采样倍数M,确定M个平移偏移量之前,包括:
获取上采样倍数M。
在一个示例性的实现过程中,获取上采样倍数M,包括:
读取系统中的默认上采样倍数值,以所述默认上采样倍数值作为上采样倍数M。
在一个示例性的实现过程中,获取上采样倍数M,包括:
接收用户输入的上采样倍数值,以所述输入的上采样倍数值作为上采样倍数M。
本发明实施例还提供一种CT系统,包括探测器、扫描床和CT设备,所述探测器包括多个探测器室及相应的处理电路;其中:
所述探测器室,用于在所述CT系统扫描过程中,探测穿过扫描对象的X射线并转换为电信号;
所述处理电路,用于将所述电信号转换成脉冲信号,采集脉冲信号的能量信息;
所述CT设备,用于:
获取对受检对象进行多排CT扫描得到的原始轴向采样生数据;
对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据;
根据所述目标轴向采样生数据进行图像重建,得到目标重建图像。
在一个示例性的实现过程中,对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据,包括:
根据上采样倍数M,确定M个平移偏移量,M为自然数,且M≥2;
分别按照所述M个平移偏移量中的每个平移偏移量,将所述原始轴向采样生数据在轴向进行平移,得到M组中间采样数据,所有中间采样中数据中任意两个相邻采样数据对应的轴向位置之差相等;
根据所述M组中间采样数据,得到目标轴向采样生数据。
在一个示例性的实现过程中,根据所述M组中间采样数据,得到目标轴向采样生数据,包括:
对于所述M组中间采样数据中的每一中间采样数据,获取该中间采样数据在所在中间采样数据组的指定邻域内的采样数据,作为基础采样数据;
根据该中间采样数据的轴向位置和各个基础采样数据的轴向位置,确定所述各个基础采样数据的权重值;
根据所述各个基础采样数据的权重值和所述各个基础采样数据的采样值,获得该中间采样数据对应的目标轴向采样生数据的目标采样值;
将该中间采样数据的轴向位置确定为该中间采样数据对应的目标轴向采样生数据的目标采样位置。
在一个示例性的实现过程中,根据该中间采样数据的轴向位置和各个基础采样数据的轴向位置,确定所述各个基础采样数据的权重值,包括:
对于每个基础采样数据,获取该中间采样数据的轴向位置和该基础采样数据的轴向位置的差值;
将所述差值作为预设权重函数的自变量的值,计算所述预设权重函数的函数值,以所述函数值作为该基础采样数据的权重值。
在一个示例性的实现过程中,根据上采样倍数M,确定M个平移偏移量之前,包括:
获取上采样倍数M。
在一个示例性的实现过程中,获取上采样倍数M,包括:
读取系统中的默认上采样倍数值,以所述默认上采样倍数值作为上采样倍数M。
在一个示例性的实现过程中,获取上采样倍数M,包括:
接收用户输入的上采样倍数值,以所述输入的上采样倍数值作为上采样倍数M。
本发明实施例还提供一种计算机可读存储介质,其上存储有计算机程序,其中,所述程序被处理器执行时实现如下操作:
获取对受检对象进行多排CT扫描得到的原始轴向采样生数据;
对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据;
根据所述目标轴向采样生数据进行图像重建,得到目标重建图像。
在一个示例性的实现过程中,对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据,包括:
根据上采样倍数M,确定M个平移偏移量,M为自然数,且M≥2;
分别按照所述M个平移偏移量中的每个平移偏移量,将所述原始轴向采样生数据在轴向进行平移,得到M组中间采样数据,所有中间采样中数据中任意两个相邻采样数据对应的轴向位置之差相等;
根据所述M组中间采样数据,得到目标轴向采样生数据。
在一个示例性的实现过程中,根据所述M组中间采样数据,得到目标轴向采样生数据,包括:
对于所述M组中间采样数据中的每一中间采样数据,获取该中间采样数据在所在中间采样数据组的指定邻域内的采样数据,作为基础采样数据;
根据该中间采样数据的轴向位置和各个基础采样数据的轴向位置,确定所述各个基础采样数据的权重值;
根据所述各个基础采样数据的权重值和所述各个基础采样数据的采样值,获得该中间采样数据对应的目标轴向采样生数据的目标采样值;
将该中间采样数据的轴向位置确定为该中间采样数据对应的目标轴向采样生数据的目标采样位置。
在一个示例性的实现过程中,根据该中间采样数据的轴向位置和各个基础采样数据的轴向位置,确定所述各个基础采样数据的权重值,包括:
对于每个基础采样数据,获取该中间采样数据的轴向位置和该基础采样数据的轴向位置的差值;
将所述差值作为预设权重函数的自变量的值,计算所述预设权重函数的函数值,以所述函数值作为该基础采样数据的权重值。
在一个示例性的实现过程中,根据上采样倍数M,确定M个平移偏移量之前,包括:
获取上采样倍数M。
在一个示例性的实现过程中,获取上采样倍数M,包括:
读取系统中的默认上采样倍数值,以所述默认上采样倍数值作为上采样倍数M。
在一个示例性的实现过程中,获取上采样倍数M,包括:
接收用户输入的上采样倍数值,以所述输入的上采样倍数值作为上采样倍数M。
对于装置和设备实施例而言,由于其基本对应于方法实施例,所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的模块可以是或者也可以不是物理上分开的,作为模块显示的部件可以是或者也可以不是物理模块,即可以位于一个地方,或者也可以分布到多个网络模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本说明书方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下,在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外,在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。
本领域技术人员在考虑说明书及实践这里申请的发明后,将容易想到本说明书的其它实施方案。本说明书旨在涵盖本说明书的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本说明书的一般性原理并包括本说明书未申请的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本说明书的真正范围和精神由下面的权利要求指出。
应当理解的是,本说明书并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。本说明书的范围仅由所附的权利要求来限制。
以上所述仅为本说明书的较佳实施例而已,并不用以限制本说明书,凡在本说明书的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本说明书保护的范围之内。
Claims (10)
1.一种图像重建方法,其特征在于,包括:
获取对受检对象进行多排CT扫描得到的原始轴向采样生数据;
对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据;
根据所述目标轴向采样生数据进行图像重建,得到目标重建图像。
2.根据权利要求1所述的方法,其特征在于,对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据,包括:
根据上采样倍数M,确定M个平移偏移量,M为自然数,且M≥2;
分别按照所述M个平移偏移量中的每个平移偏移量,将所述原始轴向采样生数据在轴向进行平移,得到M组中间采样数据,所有中间采样中数据中任意两个相邻采样数据对应的轴向位置之差相等;
根据所述M组中间采样数据,得到目标轴向采样生数据。
3.根据权利要求2所述的方法,其特征在于,根据所述M组中间采样数据,得到目标轴向采样生数据,包括:
对于所述M组中间采样数据中的每一中间采样数据,获取该中间采样数据在所在中间采样数据组的指定邻域内的采样数据,作为基础采样数据;
根据该中间采样数据的轴向位置和各个基础采样数据的轴向位置,确定所述各个基础采样数据的权重值;
根据所述各个基础采样数据的权重值和所述各个基础采样数据的采样值,获得该中间采样数据对应的目标轴向采样生数据的目标采样值;
将该中间采样数据的轴向位置确定为该中间采样数据对应的目标轴向采样生数据的目标采样位置。
4.根据权利要求3所述的方法,其特征在于,根据该中间采样数据的轴向位置和各个基础采样数据的轴向位置,确定所述各个基础采样数据的权重值,包括:
对于每个基础采样数据,获取该中间采样数据的轴向位置和该基础采样数据的轴向位置的差值;
将所述差值作为预设权重函数的自变量的值,计算所述预设权重函数的函数值,以所述函数值作为该基础采样数据的权重值。
5.根据权利要求2所述的方法,其特征在于,根据上采样倍数M,确定M个平移偏移量之前,包括:
获取上采样倍数M。
6.根据权利要求5所述的方法,其特征在于,获取上采样倍数M,包括:
读取系统中的默认上采样倍数值,以所述默认上采样倍数值作为上采样倍数M。
7.根据权利要求5所述的方法,其特征在于,获取上采样倍数M,包括:
接收用户输入的上采样倍数值,以所述输入的上采样倍数值作为上采样倍数M。
8.一种图像重建装置,其特征在于,包括:
获取模块,用于获取对受检对象进行多排CT扫描得到的原始轴向采样生数据;
上采样模块,用于对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据;
重建模块,用于根据所述目标轴向采样生数据进行图像重建,得到目标重建图像。
9.一种CT设备,其特征在于,包括:内部总线,以及通过内部总线连接的存储器、处理器和外部接口;其中,所述外部接口,用于连接CT系统的探测器,所述探测器包括多个探测器室及相应的处理电路;
所述存储器,用于存储图像重建逻辑对应的机器可读指令;
所述处理器,用于读取所述存储器上的所述机器可读指令,并执行如下操作:
获取对受检对象进行多排CT扫描得到的原始轴向采样生数据;
对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据;
根据所述目标轴向采样生数据进行图像重建,得到目标重建图像。
10.一种CT系统,其特征在于,包括探测器、扫描床和CT设备,所述探测器包括多个探测器室及相应的处理电路;其中:
所述探测器室,用于在所述CT系统扫描过程中,探测穿过扫描对象的X射线并转换为电信号;
所述处理电路,用于将所述电信号转换成脉冲信号,采集脉冲信号的能量信息;
所述CT设备,用于:
获取对受检对象进行多排CT扫描得到的原始轴向采样生数据;
对所述原始轴向采样生数据进行上采样,得到目标轴向采样生数据;
根据所述目标轴向采样生数据进行图像重建,得到目标重建图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911359736.7A CN111000578B (zh) | 2019-12-25 | 2019-12-25 | 图像重建方法、装置、ct设备及ct系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911359736.7A CN111000578B (zh) | 2019-12-25 | 2019-12-25 | 图像重建方法、装置、ct设备及ct系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111000578A true CN111000578A (zh) | 2020-04-14 |
CN111000578B CN111000578B (zh) | 2023-05-02 |
Family
ID=70118606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911359736.7A Active CN111000578B (zh) | 2019-12-25 | 2019-12-25 | 图像重建方法、装置、ct设备及ct系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111000578B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101379709A (zh) * | 2006-07-19 | 2009-03-04 | 高通股份有限公司 | 带有偏移量的∑-△调制 |
CN103054599A (zh) * | 2011-10-18 | 2013-04-24 | 株式会社东芝 | X射线ct装置及x射线ct装置的动作方法 |
US20150071533A1 (en) * | 2013-09-12 | 2015-03-12 | At&T Intellectual Property I, L.P. | Guided image upsampling using bitmap tracing |
CN108022216A (zh) * | 2017-04-14 | 2018-05-11 | 上海联影医疗科技有限公司 | 一种降低mpr伪影的系统和方法 |
CN108460740A (zh) * | 2018-03-06 | 2018-08-28 | 赛诺威盛科技(北京)有限公司 | Ct螺旋重建图像伪影去除方法 |
CN110047138A (zh) * | 2019-04-24 | 2019-07-23 | 复旦大学 | 一种磁共振薄层图像重建方法 |
-
2019
- 2019-12-25 CN CN201911359736.7A patent/CN111000578B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101379709A (zh) * | 2006-07-19 | 2009-03-04 | 高通股份有限公司 | 带有偏移量的∑-△调制 |
CN103054599A (zh) * | 2011-10-18 | 2013-04-24 | 株式会社东芝 | X射线ct装置及x射线ct装置的动作方法 |
US20150071533A1 (en) * | 2013-09-12 | 2015-03-12 | At&T Intellectual Property I, L.P. | Guided image upsampling using bitmap tracing |
CN108022216A (zh) * | 2017-04-14 | 2018-05-11 | 上海联影医疗科技有限公司 | 一种降低mpr伪影的系统和方法 |
CN108460740A (zh) * | 2018-03-06 | 2018-08-28 | 赛诺威盛科技(北京)有限公司 | Ct螺旋重建图像伪影去除方法 |
CN110047138A (zh) * | 2019-04-24 | 2019-07-23 | 复旦大学 | 一种磁共振薄层图像重建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111000578B (zh) | 2023-05-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110461228B (zh) | 使用多对比度和深度学习改进医学图像的质量 | |
CN109785243B (zh) | 基于对抗生成网络未配准低剂量ct的去噪方法、计算机 | |
CN111192228B (zh) | 图像处理方法、装置、ct设备及ct系统 | |
CN110070612B (zh) | 一种基于生成对抗网络的ct图像层间插值方法 | |
CN108776318B (zh) | 磁共振多对比度图像重建 | |
KR101351583B1 (ko) | 의료 영상 이미징 방법, 그에 따른 의료 진단 장치 및 그에 따른 기록 매체 | |
CN111110260B (zh) | 一种图像重建方法、装置及终端设备 | |
CN104700440B (zh) | 磁共振部分k空间图像重建方法 | |
JP2013506520A (ja) | コーンビームct画像におけるアーチファクト低減のための方法 | |
CN116612131A (zh) | 一种基于ADC-UNet模型的心脏MRI结构分割方法 | |
JP2021146220A (ja) | フィルタ生成方法及びフィルタ生成システム | |
CN111369465B (zh) | Ct动态图像增强方法及装置 | |
CN110246200B (zh) | 磁共振心脏电影成像方法、装置及磁共振扫描仪 | |
CN111127581A (zh) | 图像重建方法、装置、ct设备及ct系统 | |
CN111462273B (zh) | 图像处理方法、装置、ct设备及ct系统 | |
CN102890817A (zh) | 医学图像的处理方法及系统 | |
CN111000578A (zh) | 图像重建方法、装置、ct设备及ct系统 | |
WO2021051049A1 (en) | Few-view ct image reconstruction system | |
CN116071270A (zh) | 基于可变形卷积生成对抗网络的电子数据生成方法及系统 | |
WO2009004571A1 (en) | Method and apparatus for image reconstruction | |
CN111311531B (zh) | 图像增强方法、装置、控制台设备及医学成像系统 | |
Bao et al. | Random multi-channel image synthesis for multiplexed immunofluorescence imaging | |
JPWO2022163402A5 (zh) | ||
CN114494014A (zh) | 一种磁共振图像超分辨率重建方法及装置 | |
Lim et al. | Automatic Artifact Detection Algorithm in Fetal MRI |
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 |