CN102804228A - 功能图像数据增强和/或增强器 - Google Patents
功能图像数据增强和/或增强器 Download PDFInfo
- Publication number
- CN102804228A CN102804228A CN2011800144924A CN201180014492A CN102804228A CN 102804228 A CN102804228 A CN 102804228A CN 2011800144924 A CN2011800144924 A CN 2011800144924A CN 201180014492 A CN201180014492 A CN 201180014492A CN 102804228 A CN102804228 A CN 102804228A
- Authority
- CN
- China
- Prior art keywords
- image data
- functional image
- functional
- anatomical
- gradient
- 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
- 239000003623 enhancer Substances 0.000 title 1
- 238000000034 method Methods 0.000 claims abstract description 33
- 230000002708 enhancing effect Effects 0.000 claims abstract description 24
- 230000006870 function Effects 0.000 claims description 63
- 238000012545 processing Methods 0.000 claims description 54
- 238000009792 diffusion process Methods 0.000 claims description 43
- 238000009499 grossing Methods 0.000 claims description 25
- 238000005070 sampling Methods 0.000 claims description 9
- 238000012952 Resampling Methods 0.000 claims description 8
- 238000002203 pretreatment Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 6
- 230000006978 adaptation Effects 0.000 claims description 4
- 238000012805 post-processing Methods 0.000 claims description 3
- 238000005728 strengthening Methods 0.000 claims 2
- 238000003384 imaging method Methods 0.000 description 34
- 230000008569 process Effects 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 210000003484 anatomy Anatomy 0.000 description 6
- 230000003252 repetitive effect Effects 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 4
- 239000004744 fabric Substances 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 230000000877 morphologic effect Effects 0.000 description 3
- 239000000700 radioactive tracer Substances 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000002902 bimodal effect Effects 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 239000002872 contrast media Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000035790 physiological processes and functions Effects 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000004308 accommodation Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000003707 image sharpening Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003331 infrared imaging Methods 0.000 description 1
- 238000012905 input function Methods 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 239000006249 magnetic particle Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012633 nuclear imaging Methods 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 238000001959 radiotherapy Methods 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- 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
- 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/73—Deblurring; Sharpening
-
- 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
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
一种用于增强功能图像数据的方法,其包括:获得功能图像数据;获得对应于所述功能图像数据的解剖学图像数据;以及通过在所述功能图像数据和解剖学图像数据的基础上使所述功能图像数据扩散而生成增强的功能图像数据。
Description
技术领域
以下总体涉及功能成像数据的增强,将借助对多模态成像系统的特定应用来对其进行描述,其中,所述多模态成像系统的一个成像系统主要提供功能信息,所述多模态成像系统的另一个成像系统主要提供解剖学信息。这里也可以设想单独的功能成像系统和解剖学成像系统。
背景技术
就例如PET-CT、SPECT-CT和PET-MRI等双模态成像系统而言,PET或SPECT成像数据主要提供功能信息,而CT或MRI成像数据则主要提供解剖学信息。除了常用的几何配准、辐射衰减校正和融合可视化之外,在两个信息集之间几乎不采用算法交互。例如,曾经有人演示过一些尝试,试图利用更加高级的交互作用进行患者运动校正以及功能图像的协同增强。
就协同增强而言,有文献公开了试图提高PET空间解析度,增强图像对比度,校正分体积效应,降低图像噪声,以及向功能图像添加可能出现在解剖学图像中的精细结构的技术。这些技术主要以小波分解、隐马尔可夫树和正则化迭代解卷积这些数学技术为基础。这些方案中的一些,例如,隐马尔可夫树非常复杂,而且难以优化。此外,所公布的结果只表明了部分改善,其可能不足以实现可靠的临床使用。
所公开的其他针对双模态成像的相关技术采用独特的定靶CT造影剂材料和PET放射示踪剂的相关联的生理分布。借助这种技术,CT图像还将提供功能信息。其他的一些技术则尝试采用解剖学数据校正功能核成像数据中的运动伪像。还有其他的一些技术对肺结节的PET图像施行校正,其前提是能够在CT图像上准确地检测到结节。
从上文来看,存在一种有待解决的需求,即,需要新的并且不明显的使用功能成像数据和解剖学成像数据的系统和/或方法。
发明内容
本申请的各个方面解决了上述和其他问题。
根据一个方面,一种用于增强功能图像数据的方法包括:获得功能图像数据;获得对应于所述功能图像数据的解剖学图像数据;以及通过在所述功能图像数据和解剖学图像数据的基础上扩散所述功能图像数据而生成增强的功能图像数据。
根据另一方面,一种功能图像数据增强器包括处理器,该处理器通过对功能图像数据应用非线性扩散算法而生成增强的功能图像数据。
根据另一方面,一种计算机可读存储介质编码有指令,所述指令在通过计算机执行时使所述计算机的处理器执行下述步骤:生成经调适的功能图像数据,其中,基于所述功能图像数据和所述解剖学图像数据中的结构之间的相关性对功能图像数据中的结构进行调适,使之适应于解剖学图像数据中的结构。
本发明可以采取各种部件和部件布置的形式,以及各种步骤和步骤布置的形式。附图的用途仅在于对优选实施例进行图示,而不应理解为对本发明构成限制。
附图说明
图1示出了一种包括用于采集功能图像数据的成像系统和用于采集解剖学图像数据的成像系统的示范性多模态成像系统以及被配置为基于所述功能图像数据和解剖学图像数据增强功能图像数据的计算部件。
图2示出了一种包括功能图像数据增强处理部件的示范性功能图像数据增强器。
图3示出了一种示范性功能图像数据增强处理部件。
图4示出了一种示范性功能图像数据增强工作流程。
图5示出了一种功能图像数据增强的示范性方法,其包括对功能图像数据进行调适和平滑处理。
图6示出了一种用于对功能图像数据进行调适和平滑处理的示范性方法。
具体实施方式
下文涉及基于功能图像数据和解剖学图像数据来增强功能图像数据。在一个例子中,其包括采用相关非线性扩散方案增强功能图像数据,在所述方案中,通过某种方式改变功能图像数据的结构,使它们的流类似于解剖学图像数据的结构。一般而言,只要是在所述功能和解剖学图像数据中的对应形态相关的情况下都可以存在这种处理,并且两种图像数据集之间的调适的程度是可选择的。
在一些功能成像模态中,通过成像系统能够检测到的适当的示踪剂、标记、造影剂材料等的摄入水平来表示生理功能。然而,在其他功能成像模态中,通过例如在测量身体的电磁信号的MEG和EEG中直接测量生物学信号来表示功能信息。应当认识到,此处采用的摄入一词通常描述任何相关的功能图像信息。
应当认识到所述功能和解剖学图像数据可以来自于多模态成像系统,例如,PET-CT、SPECT-CT、PET-MRI和具有光学成像的临床前CT,并且/或者可以来自于单独的成像系统。其他适当的成像模态包括US、x射线放射照相术、功能MRI、磁性粒子成像(MPI)、红外成像、脑磁波描记(MEG)、脑电波描记(EEG)和/或其他成像系统。
图1示出了一种示范性系统100。所述系统包括至少一个成像系统102。所示出的成像系统102包括至少两个子成像系统,即,至少生成功能图像数据的第一成像系统104和至少生成解剖学图像数据的第二成像系统106。在另一实施例中,第一成像系统104反过来或者另外又生成解剖学图像数据,并且/或者第二成像系统106反过来或者另外又生成功能图像数据。在又一实施例中,所述第一和第二成像系统104和106是分开的单独的系统。
系统100还包括计算部件108,例如,工作站、台式计算机、膝上型电脑、手持计算设备等。所述计算系统包括一个或多个处理器110和存储器112,例如,所述存储器是采用数据114和/或指令116对其进行编码的计算机可读介质,在通过处理器110执行时,所述数据和/或指令将使计算系统108执行各种操作。至少采用功能图像数据增强指令118对所示出的存储器112进行了编码。输入/输出120提供了用于接收诸如算法参数的输入和/或将诸如受到增强的功能图像数据传送至(例如)一个或多个远程装置122的接口。
适当的远程装置122的例子包括但不限于显示器、计算机、数据库、图片存档及通信系统(PACS)和/或其他装置。受到增强的功能图像数据能够通过联系初始功能图像数据、解剖学图像数据和/或其他图像数据来进行可视化和/或做进一步处理。或者或此外,所述功能图像数据增强指令118能够通过第一或第二成像系统104或106、另一成像系统、PACS系统、其他计算部件和/或其他系统的至少其中之一来存储和/或实施。
以下进行更详细的描述,在一个例子中,所述功能图像数据增强指令118使所述功能图像数据中的结构适合于所述解剖学图像数据中的结构,例如,至少在对应的模态相关的地方是这样。其结果是,相对于省略了功能图像数据增强指令的配置,获得了功能信息和解剖学信息之间的更好的调适。此外,其可以在(例如)背景摄入未发生变化的地方以更高的信噪比(SNR)实现功能图像数据的图像锐化。
因而,能够实现功能图像数据中的形态结构的更好的描绘,与此同时获得更加准确的功能图像值。就作为功能模态的一个代表实例的核医学而言,这样做可以得到改进的临床应用,例如,对小的转移的检测和功能评估、勾勒肿瘤轮廓以用于放射治疗规划和后续措施以及对放射示踪剂摄入的定量评估。其还可以促进扫描时间和/或放射示踪剂剂量的降低。
图2示出了基于解剖学图像数据而增强功能图像数据的功能图像数据增强器202。功能图像数据增强器202能够通过处理器110(图1)执行功能图像数据增强指令118(图1)的和/或其他方式来执行。所示出的功能图像数据增强器202接收功能图像数据、解剖学图像数据以及任选的一个或多个参数。所述功能图像数据和/或解剖学图像数据可以包括重构的图像数据。在一个例子中,对所述重构的图像数据应用诸如衰减校正、空间配准等的各种处理。
所述功能图像数据增强器202包括对所述功能和/或解剖学图像数据进行预处理的相关增强前处理部件204。例如,在一个非限制性例子中,所述相关增强前处理部件204对所述功能和/或解剖学图像数据进行重新采样。例如,在一个例子中,其可以包括对功能图像数据(或解剖学图像数据)重新采样,以使其适合所述解剖学图像数据(或所述功能图像数据)相同的体素大小和位置(例如,沿x、y和/或z方向)。在一个例子中,对具有相对较为粗略的采样网格的图像数据进行预处理,从而使其图像数据适合具有相对较为精细的采样网格的图像的采样网格。在另一个例子中,对两种图像数据集都重新采样。在又一个例子中,对两种图像数据集都不重新采样。
在另一非限制性实施例中,相关增强前处理部件204对所述图像数据集的至少其中之一,例如,对所述功能图像数据集应用锐化算法。适当的锐化算法包括图像解卷积技术,例如,Lucy-Ri chardson解卷积或者其他基于MLEM的解卷积。可以通过2D-轴体层的锐化和/或1D轴向锐化或者完全的3D锐化来应用锐化。也可以设想其他算法。在一个例子中,所述锐化在只增加适度的噪声的情况下使功能图像数据显著锐化。可以将所述锐化算法等同地应用于整个功能图像数据体或其预定部分。
在另一实施例中,省略相关增强前处理部件204。
所述功能图像数据增强器202还包括相关增强处理部件206。在图示的实施例中,相关增强处理部件206基于解剖学图像数据对经过预处理的功能图像数据应用相关非线性扩散算法。在一个例子中,其可以包括按照非线性的方式对经过预处理的功能图像数据进行平滑处理、边缘增强或形变处理中的一者或多者,其中,所述功能图像数据的梯度和所述解剖学图像数据的梯度控制所述平滑处理/边缘增强/形变。
出于(例如)增强算法的内部迭代的目的,可以按照一种不同的方式或者按照类似的方式使解剖学图像数据平滑。所述相关非线性扩散算法可以是(例如)具有预定次数的迭代的单迭代或多迭代算法。方程1示出了这样的算法的非限制性例子:
方程1:
采用方程1,在图像数据处理过程中,生理功能信息既作为解剖学信息的空间梯度的函数又作为功能信息的空间梯度的函数而以迭代的方式发生变化。所述迭代变化通过这样一种方式发生,即,使得所述功能信息的形态结构倾向于适应解剖学信息的相关形态结构。这种适配趋势随着结构相关性的增大而增大。
在图示的实施例中,将各种输入参数用于图像数据增强。例如,一个适当的参数包括锐度参数。这样的参数可以界定功能图像数据锐化的水平以及对具有满足预定摄入范围的摄入的小的图案的相关强度校正。另一个适当的参数包括调适参数。这一参数可以确定功能图像数据到解剖学图像数据的取决于相关性的形态调适的水平。
因而,这一参数控制赋予解剖学形态的相对重要性(或权重)。另一个适当的参数包括背景参数。这一参数能够用来界定功能图像数据内的被看作是不应受到增强的背景摄入的强度范围。这里也可以设想额外的和/或替代的参数。在一个非限制性实施例中,对于至少一个所述参数而言,用户可以将参数设定为,或者可以采用默认文件将参数设定为开启或关闭、低、中、高或其他。
功能图像数据增强器202还包括相关增强后处理部件208。其中,相关增强前处理部件204对所述图像数据进行上采样和/或尺寸重设,可以采用相关增强后处理部件208在初始图像数据体素尺寸和/或位置的基础上对图像数据进行下采样和/或尺寸重设。在一个例子中,采用反走样算法对所述图像数据进行重新采样和/或尺寸重设。
功能图像数据增强器202的输出至少包括增强的功能图像数据。如文中讨论的,可以采用初始功能图像数据、解剖学图像数据和/或其他图像数据使增强的功能图像数据可视化和/或对其做进一步处理。
应当认识到,如果局部功能图像数据形态与对应的解剖学图像数据形态相似或者与之相关,那么使功能图像数据增强的概率就高。然而,如果不存在这样的局部相关性,那么也不会使所述局部功能图像数据劣化。在缺少相关性的情况下,例如,在表现功能的生理特征不具有对应的显著解剖学形态的情况下,这一特点可能是有利的,而且其还能够缓解由患者运动或者功能形态与解剖学形态之间的不准确空间配准导致的伪像。
此外,所述增强不需要任何分割或复杂的形态分析,因为通过图像梯度反映了形态相关性。因而,可以按照相同的方式对整个体积应用所述增强而不必识别具体的结构。此外,内部和/或用户参数未强加锐利的边界或阈值,从而允许灵活直观的调整。
图3示出了图2的相关增强处理部件206的示范性实施例。主处理块302包括调适处理块304和平滑处理块306,以及处理块参数设置器308。所述处理块参数设置器308包括指示所述主块302、调适处理块304和/或平滑处理块306的迭代次数的参数。
在图示的实施例中,处理块参数设置器308伴随主块302的每次迭代更新所述调适处理块304和/或平滑处理块306的一个或多个参数。在一个例子中,所述处理块参数设置器308可以相对于主块302的前一次迭代提高主块302的每次迭代中的所述块304的迭代的次数。处理块参数设置器308还可以相对于302的前一次迭代更新302的每次迭代当中的304的非线性扩散算子中的非线性项,以降低扩散通量具有正斜率的梯度范围。这一受到调节的参数更新可以实现随梯度幅度的降低而执行相对更强的扩散。其可以促进针对锐化和针对平滑的调适作用的均衡。其还可以避免增强过程中的伪(spurious)分段效应。
调适处理块304包括解剖学图像数据(AID)梯度确定器310、功能图像数据(FID)扩散器312和解剖学图像数据(AID)扩散器314。调适处理块304调适功能图像数据形态,使之适合解剖学图像数据形态。调适处理块304进行预定次数的迭代运算,例如,所述次数是由处理块参数设置器308确定的,所述迭代次数随着主块302的每次迭代而增大。
AID梯度确定器310确定解剖学图像数据的梯度。采用所述梯度控制功能图像数据和解剖学图像数据二者的扩散。功能图像数据扩散器312基于所述梯度使所述功能图像数据扩散。在一个非限制性实例中,在与所述解剖学图像数据的扩散相关的情况下使所述功能图像数据扩散。所述解剖学图像数据扩散器314基于所述梯度使所述解剖学图像数据扩散。在一个非限制性实例中,还通过非线性方式使所述解剖学图像平滑。调适处理块304的输出是中间增强功能图像数据。
平滑处理块306包括功能图像数据(FID)梯度确定器316和功能图像数据(FID)扩散器318。在图示的实施例中,平滑处理块306基于一种非线性方案使所述中间功能图像数据平滑或使其模糊。在所示的例子中,使平滑处理块306进行预定次数的迭代运算,例如,所述次数是由处理块参数设置器308确定的,所述迭代次数随着主块302的每次迭代而增大。
处理块参数设置器308还可以相对于302的前一次迭代而更新302的每次迭代当中的306的非线性扩散算子中的非线性项,以降低扩散通量具有正斜率的梯度范围。这一受到调节的参数更新可以实现:(例如)根据调适块304的扩散处理,随着梯度幅度的降低而执行相对更强的扩散。
在一个例子中,所述平滑处理有助于在不引入过多的噪声和/或对解剖学图像数据的看起来过于像伪像的调适的情况下创建所述功能图像数据的更加逼真的外观。FID梯度确定器316确定调适块304输出的功能图像数据(中间增强功能图像数据)的梯度。采用所述梯度控制受到调适的功能解剖学图像数据的扩散。功能图像数据扩散器318基于所述梯度使所述增强的功能图像数据扩散。
主块302通常进行多次迭代,对于主块302的每次迭代,调适处理块304和/或平滑处理块306都进行多次迭代。在所有的迭代之后,主块302的输出是增强的功能图像数据。
如上文指出的,计算部件108(图1)应用相关非线性扩散算法来增强功能图像数据。下文提供了适当的相关非线性扩散算法的非限制性例子。
在方程2中示出了所述非线性扩散方案的一般微分方程:
方程2:
其中,I是图像像素值的阵列,t与迭代相关,g的定义可以如方程3中所示:
方程3:
其中,λ是变量,Cm通过方程4获得:
方程4:
Cm=root(1-exp(-x)-m·x·exp(-x))。
就一种方案而言,相关非线性扩散使得一幅图像中的变化与第二幅图像中的梯度相关。作为方程5的函数,其实现了第一幅图像到第二幅图像的调适:
方程5:
为了使功能图像数据适合于解剖学图像数据的图案,将功能图像数据用于I1,将解剖学图像数据用于I2。
表1示出了用于所述扩散的示范性伪码。
表1
第1行中的表达式定义了输入功能图像数据和解剖学图像数据。第2行开始了所述主块302的迭代(图3)。第3行开始所述调适处理块304的内部相关扩散迭代。第4行涉及解剖学图像数据的作为其自身的梯度的函数的扩散。第5行涉及功能图像数据的作为解剖学图像数据梯度的函数的扩散。第7行将中间增强功能图像定义为向所述平滑处理块306的输入。
第8行开始所述平滑处理块306的内部扩散迭代。第9行涉及功能图像数据的作为功能图像数据梯度的函数的扩散。第11行指示了用于块302的下一外部迭代的两图像值。第12行描述了更新内部迭代的次数和非线性参数(λ)的趋势。第14行指示了新的功能图像数据的最终结果。(第4、5、9行内的)因数s为增量因数(例如,常数s=0.2)。
图4示出了功能图像数据增强的示范性工作流程。在402中,诸如医生、放射科医师、技师等的经授权的用户选择功能图像数据和解剖学图像数据,例如,PET图像数据和对应的CT图像数据。若不然也可以自动选择输入数据。
在404中,获得一个或多个处理参数。如文中讨论的,适当的参数包括锐度、调适和背景。锐度参数决定功能图像数据的锐化水平,调适参数定义功能图像数据到解剖学图像数据的取决于相关性的形态调适的水平,背景参数定义功能图像数据中被看作是背景摄入的强度的粗略范围。如文中所述,这些参数可以包括默认的和/或用户定义的参数。
在406中,基于所述解剖学图像数据和所述参数增强所述功能图像数据。增强的功能图像数据可以包括将所述图像数据标识为增强的功能图像数据的标记。在一个例子中,其有助于区分初始图像数据与增强的图像数据。
在408中,对增强的功能图像数据进行显示和/或否则对其进行处理。可以联系初始功能图像数据、解剖学图像数据和/或其他图像数据对所述增强的功能图像数据进行显示和/或否则对其进行处理。
图5示出了例如与图4的操作406有关的功能图像数据增强的示范性方法。
在502中,对功能图像数据和/或解剖学图像数据进行重新采样和/或尺寸重设,以使其适合相同的体素大小和位置。
在504中,对经处理的功能图像数据的整体或者预定子集(例如,等同地)应用锐化算法。
在506中,对经锐化的功能图像数据应用相关非线性扩散,以生成增强的功能图像数据。在一个例子中,其包括对功能图像数据进行调适和平滑处理,例如,如下文结合图6所述,或者采取其他方式。
在508中,对经处理的功能图像数据进行下采样和/或将其尺寸重设至初始功能图像数据体素大小和/或位置。
在510中,对增强的功能图像数据加以显示和/或对其加以利用,例如,单独利用或者联系初始功能图像数据、解剖学图像数据和/或其他数据对其加以利用。
图6示出了例如结合图5的操作506对功能图像数据进行调适和平滑处理的示范性方法。
在602中,确定解剖学图像数据的梯度。
在604中,采用所述梯度使所述功能图像数据扩散(使之平滑和/或对其进行边缘增强和/或使之发生形变)。
在606中,判断是否针要对经扩散的功能图像数据重复操作602和604。
如果是,那么在608中,采用所述梯度使所述解剖学图像数据扩散,并利用经扩散的功能图像数据和经扩散的解剖学图像数据重复操作602-606。
如果否,那么在610中确定经扩散的功能图像数据的梯度。
在612中,采用所述梯度使经扩散的功能图像数据扩散。
在614中,判断是否要针对经扩散的功能图像数据重复操作610和612。同样地,能够预置迭代次数。
如果是,那么采用来自操作614的输出功能图像数据重复操作610-614。
如果否,那么在616中,采用来自上一步骤612的经处理的功能图像数据、来自上一步骤608的经处理的解剖学图像数据和更新后的处理参数来判断是否要重复操作602-614。
在618中输出增强的功能图像数据。
上述操作可以通过计算机可读指令来执行,在通过计算机处理器执行所述指令时,所述指令将使所述处理器执行文中描述的操作。在这种情况下,将所述指令存储在与相关计算机关联和/或能被其访问的诸如内存的计算机可读存储介质中。
已经参考优选实施例描述了本发明。在阅读并理解了前述详细说明后,本领域技术人员可以想到修改和变化。这意味着,应当将本发明构造为包括所有此类落在所附权利要求及其等同形式的范围内的修改和变化。
Claims (27)
1.一种用于增强功能图像数据的方法,包括:
获得功能图像数据;
获得对应于所述功能图像数据的解剖学图像数据;以及
通过基于所述功能图像数据和所述解剖学图像数据使所述功能图像数据扩散而生成增强的功能图像数据。
2.根据权利要求1所述的方法,生成所述增强的功能图像数据包括:
基于所述功能图像数据和所述解剖学图像数据中的结构之间的相关性对所述功能图像数据中的结构进行调适,使之适合于所述解剖学图像数据中的结构。
3.根据权利要求1或2所述的方法,还包括:
采用相关非线性扩散算法生成所述增强的功能图像数据。
4.根据权利要求3所述的方法,还包括:
以迭代的方式应用所述算法,其中,对于相继的迭代而言,相对于前一迭代以更高的强度对低梯度区域应用所述算法。
5.根据权利要求3或4所述的方法,其中,对于相继的迭代而言,相对于前一迭代,降低了用于区分低梯度区域和高梯度区域的梯度值。
6.根据权利要求1所述的方法,还包括:
在增强所述功能图像数据之前使所述功能图像数据锐化。
7.根据权利要求6所述的方法,其中,对所述功能图像数据均匀地应用锐化。
8.根据权利要求1、3、4和5中的任何一项所述的方法,还包括:
调适功能图像数据形态,使之适合于解剖学图像数据形态。
9.根据权利要求8所述的方法,其中,调适包括:
a)确定用于所述解剖学图像数据的梯度;
b)基于所述梯度对所述解剖学图像数据进行平滑处理;以及
c)基于所述梯度和经平滑处理的解剖学图像数据,对所述功能图像数据进行平滑处理和/或边缘增强和/或使之发生形变。
10.根据权利要求9所述的方法,还包括:
采用经平滑处理的解剖学图像数据和经平滑处理的/边缘增强的/形变的功能图像数据对操作a)-c)进行预定次数的重复。
11.根据权利要求10所述的方法,还包括:
a)确定中间增强功能图像数据的梯度;以及
b)基于所述功能图像数据梯度对所述中间增强功能图像数据进行平滑处理。
12.根据权利要求11所述的方法,还包括:
采用经平滑处理的功能图像数据对a)-b)进行预定次数的重复。
13.根据权利要求12所述的方法,还包括:
在增强所述功能图像数据之前对所述功能图像数据或所述解剖学图像数据的至少其中之一进行上采样。
14.根据权利要求13所述的方法,还包括:
基于对初始功能图像数据的采样对所述增强的功能图像数据进行下采样。
15.一种功能图像数据增强器(202),包括:
通过对功能图像数据应用非线性扩散算法而生成增强的功能图像数据的处理器(206)。
16.根据权利要求15所述的增强器,还包括:
在应用所述相关非线性扩散算法之前使功能图像数据锐化的相关增强前处理部件(204)。
17.根据权利要求15或16所述的增强器,还包括:
调适处理块(304),其基于所述功能图像数据和所述解剖学图像数据中的结构之间的相关性对功能图像数据中的结构进行调适,使之适合于解剖学图像数据中的结构。
18.根据权利要求17所述的增强器,其中,所述调适处理块(304)以迭代的方式应用所述算法,对于相继的迭代而言,相对于前一迭代以更高的强度对低梯度区域应用所述算法。
19.根据权利要求18所述的增强器,其中,对于相继的迭代而言,相对于前一迭代,降低了用于区分低梯度区域和高梯度区域的梯度值。
20.根据权利要求17到19中的任何一项所述的增强器,还包括:
对经调适的功能图像数据进行平滑处理以生成增强的功能图像数据的平滑处理块(306)。
21.根据权利要求15到20中的任何一项所述的增强器,其中,所述处理器(206)基于一个或多个处理参数生成增强的功能图像数据,所述参数包括选自由下述参数构成的集合中的至少一个参数:锐度参数、调适参数和背景参数。
22.根据权利要求21所述的增强器,其中,所述锐度参数界定功能图像数据锐化的水平以及具有满足预定范围的功能指示水平的小的图案的相关强度校正。
23.根据权利要求21到22中的任何一项所述的增强器,其中,所述调适参数决定着功能图像数据到解剖学图像数据的取决于相关性的形态调适的水平。
24.根据权利要求21到23中的任何一项所述的增强器,其中,所述背景参数界定所述功能图像数据中的被看作是不应被增强的背景功能指示水平的强度范围。
25.根据权利要求15到24中的任何一项所述的增强器,还包括:
在应用所述相关非线性扩散算法之前对所述功能图像数据进行重新采样和/或尺寸重设的相关增强前处理部件(204);以及
在应用所述相关非线性扩散算法之后对所述功能图像数据进行重新采样和/或尺寸重设的相关增强后处理部件(208)。
26.一种采用编码有指令的计算机可读存储介质,所述指令在通过计算机的处理器执行时使所述计算机执行下述步骤:
生成经调适的功能图像数据,其中,基于所述功能图像数据和所述解剖学图像数据中的结构之间的相关性对功能图像数据中的结构进行调适,使之适合于解剖学图像数据中的结构。
27.根据权利要求26所述的计算机可读存储介质,其中,所述指令在通过所述处理器执行时使所述计算机执行下述步骤:
基于所述经调适的功能图像数据的梯度对所述经调适的功能图像数据进行平滑处理。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US31505310P | 2010-03-18 | 2010-03-18 | |
US61/315,053 | 2010-03-18 | ||
PCT/IB2011/050664 WO2011114243A1 (en) | 2010-03-18 | 2011-02-17 | Functional image data enhancement and/or enhancer |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102804228A true CN102804228A (zh) | 2012-11-28 |
CN102804228B CN102804228B (zh) | 2015-08-19 |
Family
ID=44065678
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201180014492.4A Expired - Fee Related CN102804228B (zh) | 2010-03-18 | 2011-02-17 | 功能图像数据增强和/或增强器 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8989464B2 (zh) |
EP (1) | EP2548172A1 (zh) |
JP (1) | JP5815573B2 (zh) |
CN (1) | CN102804228B (zh) |
WO (1) | WO2011114243A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104410929A (zh) * | 2014-12-19 | 2015-03-11 | 三星电子(中国)研发中心 | 字幕图像的处理方法和装置 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012056880A (ja) * | 2010-09-08 | 2012-03-22 | Idemitsu Kosan Co Ltd | インドロカルバゾール化合物、有機エレクトロルミネッセンス素子用材料、及びそれを用いた有機エレクトロルミネッセンス素子 |
US8942445B2 (en) * | 2012-09-14 | 2015-01-27 | General Electric Company | Method and system for correction of lung density variation in positron emission tomography using magnetic resonance imaging |
KR101461099B1 (ko) * | 2012-11-09 | 2014-11-13 | 삼성전자주식회사 | 자기공명영상장치 및 기능적 영상획득방법 |
WO2019090023A1 (en) * | 2017-11-03 | 2019-05-09 | General Electric Company | System and method for interactive representation learning transfer through deep learning of feature ontologies |
CN112614137B (zh) * | 2019-10-06 | 2022-08-30 | 山东九州信泰信息科技股份有限公司 | 基于对象数量的参数调控系统 |
US11309072B2 (en) | 2020-04-21 | 2022-04-19 | GE Precision Healthcare LLC | Systems and methods for functional imaging |
KR102342940B1 (ko) * | 2021-06-23 | 2021-12-27 | 주식회사 포바이포 | 선예도 개선을 위한 심층 기울기 사전정보기반 l0 스무딩 방법 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040260170A1 (en) * | 2003-06-20 | 2004-12-23 | Confirma, Inc. | System and method for adaptive medical image registration |
CN101178814A (zh) * | 2007-11-30 | 2008-05-14 | 华南理工大学 | 一种融合解剖与功能成像信息数据场的半透明体绘制方法 |
CN101310687A (zh) * | 2007-05-22 | 2008-11-26 | 西门子公司 | 数据分析的方法 |
Family Cites Families (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9003A (en) * | 1852-06-08 | Improvement in wheel-cultivators | ||
US6021A (en) * | 1849-01-09 | Cast-iron cab-wheel | ||
EP0731956A4 (en) * | 1993-11-30 | 1997-04-23 | Arch Dev Corp | AUTOMATIC METHOD AND SYSTEM FOR ALIGNING AND CORRELATING TWO DIFFERENT MODALITIES |
US6740883B1 (en) | 1998-08-14 | 2004-05-25 | Robert Z. Stodilka | Application of scatter and attenuation correction to emission tomography images using inferred anatomy from atlas |
US6878941B2 (en) | 2002-04-09 | 2005-04-12 | Elgems Ltd. | Gamma camera and CT system |
US6490476B1 (en) * | 1999-10-14 | 2002-12-03 | Cti Pet Systems, Inc. | Combined PET and X-ray CT tomograph and method for using same |
GC0000235A (en) * | 2000-08-09 | 2006-03-29 | Shell Int Research | Processing an image |
US6728424B1 (en) * | 2000-09-15 | 2004-04-27 | Koninklijke Philips Electronics, N.V. | Imaging registration system and method using likelihood maximization |
EP1502237A2 (en) * | 2002-04-03 | 2005-02-02 | Segami S.A.R.L. | Image registration process |
DE10231061A1 (de) | 2002-07-10 | 2004-01-22 | Philips Intellectual Property & Standards Gmbh | Verfahren und System zur Verbesserung des Informationsgehaltes in einem Bild |
GB2391125B (en) * | 2002-07-19 | 2005-11-30 | Mirada Solutions Ltd | Registration of multi-modality data in imaging |
WO2004089218A1 (ja) * | 2003-04-04 | 2004-10-21 | Hitachi Medical Corporation | 機能画像の表示方法及び装置 |
US20050031176A1 (en) | 2003-08-08 | 2005-02-10 | Hertel Sarah R. | Method and apparatus of multi-modality image fusion |
EP1828977B1 (en) | 2004-11-17 | 2013-09-11 | Koninklijke Philips Electronics N.V. | Restoration of the nuclear medicine 2d planar image by iterative constrained deconvolution |
JP4969445B2 (ja) * | 2005-06-15 | 2012-07-04 | 株式会社日立メディコ | 磁気共鳴イメージング装置 |
US7348564B2 (en) | 2005-12-12 | 2008-03-25 | General Electric Company | Multi modality imaging methods and apparatus |
AU2007244705A1 (en) | 2006-04-27 | 2007-11-08 | Barnes-Jewish Hospital | Detection and imaging of target tissue |
US8160676B2 (en) * | 2006-09-08 | 2012-04-17 | Medtronic, Inc. | Method for planning a surgical procedure |
US20080095414A1 (en) | 2006-09-12 | 2008-04-24 | Vladimir Desh | Correction of functional nuclear imaging data for motion artifacts using anatomical data |
JP5462865B2 (ja) | 2008-05-15 | 2014-04-02 | コーニンクレッカ フィリップス エヌ ヴェ | 不完全な解剖学的画像を補償するための非減衰補正pet放出画像の使用 |
-
2011
- 2011-02-17 JP JP2012557629A patent/JP5815573B2/ja not_active Expired - Fee Related
- 2011-02-17 US US13/583,934 patent/US8989464B2/en not_active Expired - Fee Related
- 2011-02-17 CN CN201180014492.4A patent/CN102804228B/zh not_active Expired - Fee Related
- 2011-02-17 WO PCT/IB2011/050664 patent/WO2011114243A1/en active Application Filing
- 2011-02-17 EP EP11711657A patent/EP2548172A1/en not_active Withdrawn
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040260170A1 (en) * | 2003-06-20 | 2004-12-23 | Confirma, Inc. | System and method for adaptive medical image registration |
CN101310687A (zh) * | 2007-05-22 | 2008-11-26 | 西门子公司 | 数据分析的方法 |
CN101178814A (zh) * | 2007-11-30 | 2008-05-14 | 华南理工大学 | 一种融合解剖与功能成像信息数据场的半透明体绘制方法 |
Non-Patent Citations (1)
Title |
---|
杜丽娅: "基于B样条的医学图像弹性配准", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104410929A (zh) * | 2014-12-19 | 2015-03-11 | 三星电子(中国)研发中心 | 字幕图像的处理方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
US8989464B2 (en) | 2015-03-24 |
JP2013521919A (ja) | 2013-06-13 |
JP5815573B2 (ja) | 2015-11-17 |
EP2548172A1 (en) | 2013-01-23 |
US20130004038A1 (en) | 2013-01-03 |
CN102804228B (zh) | 2015-08-19 |
WO2011114243A1 (en) | 2011-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Würfl et al. | Deep learning computed tomography: Learning projection-domain weights from image domain in limited angle problems | |
Wang et al. | A review on medical imaging synthesis using deep learning and its clinical applications | |
CN102804228B (zh) | 功能图像数据增强和/或增强器 | |
CN110337673B (zh) | 基于深度学习的在断层摄影重建中估计数据的方法和系统 | |
Yokota et al. | Dynamic PET image reconstruction using nonnegative matrix factorization incorporated with deep image prior | |
US11200711B2 (en) | Smart filtering for PET imaging including automatic selection of filter parameters based on patient, imaging device, and/or medical context information | |
Naqa et al. | Deblurring of breathing motion artifacts in thoracic PET images by deconvolution methods | |
Le Pogam et al. | Evaluation of a 3D local multiresolution algorithm for the correction of partial volume effects in positron emission tomography | |
Joemai et al. | Metal artifact reduction for CT: Development, implementation, and clinical comparison of a generic and a scanner‐specific technique | |
CN103959329A (zh) | 图像域去噪 | |
CN104252714A (zh) | 时变数据的重建 | |
Riblett et al. | Data‐driven respiratory motion compensation for four‐dimensional cone‐beam computed tomography (4D‐CBCT) using groupwise deformable registration | |
Pei et al. | Non-rigid craniofacial 2D-3D registration using CNN-based regression | |
da Costa-Luis et al. | Micro-networks for robust MR-guided low count PET imaging | |
Kazantsev et al. | An anatomically driven anisotropic diffusion filtering method for 3D SPECT reconstruction | |
Zhang et al. | PET image reconstruction using a cascading back-projection neural network | |
Xu et al. | Statistical iterative reconstruction to improve image quality for digital breast tomosynthesis | |
Li et al. | Multienergy cone-beam computed tomography reconstruction with a spatial spectral nonlocal means algorithm | |
CN106251295A (zh) | 考虑轮廓显著性数据的结果图像的重建 | |
Abascal et al. | A novel prior-and motion-based compressed sensing method for small-animal respiratory gated CT | |
Jiang et al. | Fast four‐dimensional cone‐beam computed tomography reconstruction using deformable convolutional networks | |
Dutta et al. | Deep learning framework to synthesize high-count preclinical PET images from low-count preclinical PET images | |
Chillaron et al. | Evaluation of image filters for their integration with LSQR computerized tomography reconstruction method | |
Daum | Model-constrained non-rigid registration in medicine | |
Us et al. | Combining dual-tree complex wavelets and multiresolution in iterative CT reconstruction with application to metal artifact reduction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150819 Termination date: 20170217 |
|
CF01 | Termination of patent right due to non-payment of annual fee |