CN103918004A - 金属伪影校正算法的自适应应用 - Google Patents
金属伪影校正算法的自适应应用 Download PDFInfo
- Publication number
- CN103918004A CN103918004A CN201280054611.3A CN201280054611A CN103918004A CN 103918004 A CN103918004 A CN 103918004A CN 201280054611 A CN201280054611 A CN 201280054611A CN 103918004 A CN103918004 A CN 103918004A
- Authority
- CN
- China
- Prior art keywords
- image
- neighborhood
- corrector
- corrective action
- initial
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 32
- 239000002184 metal Substances 0.000 title description 10
- 230000003044 adaptive effect Effects 0.000 title description 3
- 230000009471 action Effects 0.000 claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 21
- 238000012545 processing Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 15
- 230000003321 amplification Effects 0.000 claims description 5
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 5
- 230000001419 dependent effect Effects 0.000 claims description 4
- 230000006870 function Effects 0.000 description 12
- 238000002591 computed tomography Methods 0.000 description 7
- 239000011159 matrix material Substances 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003702 image correction Methods 0.000 description 3
- 230000002238 attenuated effect Effects 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 239000007943 implant Substances 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 210000004872 soft tissue Anatomy 0.000 description 2
- 230000003936 working memory Effects 0.000 description 2
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000002195 synergetic effect Effects 0.000 description 1
- 230000007704 transition Effects 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/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/77—Retouching; Inpainting; Scratch removal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- 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]
-
- 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/10116—X-ray image
-
- 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/20—Special algorithmic details
- G06T2207/20004—Adaptive image processing
- G06T2207/20012—Locally adaptive
-
- 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/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
-
- 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)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
一种针对图像伪影来校正图像的装置和方法。初始图像由图像伪影校正器(190)校正。如此校正的样本校正图像与所述初始图像进行比较以获得关于校正动作的信息。之后所述校正动作由控制器(140)自适应地重新应用以获得改进的经校正的图像,由此确保先前存在的伪影被去除,并且新伪影的产生得以避免。
Description
技术领域
本发明涉及一种用于针对图像伪影校正图像的装置,一种针对图像伪影校正图像的方法,一种用于针对图像伪影校正图像的医学图像处理系统,一种计算机程序产品,以及一种计算机可读介质。
背景技术
金属植入物或其他高密度部分经常会导致CT图像中的图像伪影。在大多数情况下减少图像伪影的金属伪影减少(MAR)算法是已知的。一种MAR算法可从US7991243中获知。
然而,仍然存在这样的问题,即使目前已知的最佳算法也偶然地产生新伪影。
发明内容
因此,需要改进图像伪影减少算法。
本发明的目的是通过独立权利要求的主题解决的,其中进一步的实施例被包含在从属权利要求中。
应当注意,下面描述的各方面也适用于针对图像伪影校正图像的方法、针对图像伪影校正图像的医学图像处理系统、计算机程序产品以及计算机可读介质。
根据一个方面,提供一种用于针对伪影校正图像的装置。
所述装置包括用于接收初始图像的输入单元,用于处理所述图像以进行校正的处理单元,以及用于输出所述初始图像的经校正的版本的输出单元。
所述处理单元包括校正采样器,所述校正采样器被配置为将所述初始图像转发到伪影校正器并从那里接收经校正的图像的样本。样本图像是由所述伪影校正器对所述初始图像应用的校正动作的结果。
所述处理单元还包括比较器,所述比较器被配置为比较所述初始图像和经校正的样本图像以便建立表示所述伪影校正器对所述初始图像的所述校正动作的校正器图像。
所述处理单元还包括伪影校正控制器,所述伪影校正控制器被配置为将各自的校正动作以加权方式自适应地重新应用于所述初始图像中的多个图像点中的任意一个。针对每个图像点的用于对校正动作进行加权的权重或权重因子与该图像点周围的邻域中的图像信息能由所述校正器图像中的对应邻域中的对应图像信息补偿的程度相关。
之后,被如此自适应地重新校正的初始图像由所述输出单元输出为所述经校正的图像。
换言之,所述校正动作由所述控制器自适应地重新应用以获得改进的最终经校正的图像,由此确保先前存在的伪影被去除,并且新伪影的产生得以避免。
所述装置允许应用现有的伪影校正器来实施已知的金属伪影减少(MAR)算法。
所述装置能够作为“附加件”添加到现有MAR系统上。所述装置是“元图像伪影校正器”,因为其会检查由现有MAR系统输出的校正结果,并且之后通过使用校正器图像作为“路线图”自适应地重新应用校正动作来改进校正结果。所述校正器图像记录现有MAR系统对用于初始图像的校正动作的估计。
所述装置以“双”自适应的方式(即在定量和空间意义上)起作用:控制器确定在何处,即在初始图像中的哪些图像点处,要应用估计的校正动作,以及在该区域要应用多少校正动作,即衰减或放大多少。
所述装置以两阶段的方式工作:在第一,样本轮次中,伪影校正器被全局应用于初始图像以获得样本经校正的图像,并因此收集关于该算法对初始图像的校正动作的信息。
校正动作被记录在校正器图像中。使用校正器图像,记录于其中的校正动作在第二,最终轮次中被自适应且选择性地应用,但这次校正动作仅被局部应用于这样的图像区域,在所述图像区域中,“校正图像”显示也存在于原始的初始图像中的特征。
所述控制器在邻域方向起作用。以图像平面的先前定义的区域(其可以包括对于图像帧周围的边界区域而言安全的整个图像平面)中的每个图像点为中心,定义了子区域或邻域(“补片”)。之后该补片被映射到校正器图像以定义与初始图像邻域相对应的校正器图像邻域,以便为每个图像点获得一对邻域。
补偿程度(以下称为“程度”)测量两个邻域中的图像信息彼此补偿或相互抵消的适当性。该程度或与该程度相关联的校正权重指示图像信息在该对邻域上为镜像的或相似的范围。如此,该程度是各自的初始图像邻域和校正器图像邻域对中的每个的属性。因为邻域是彼此对应的,所以任一邻域都能够认为“具有”该程度。
邻域的高补偿程度或能力导致校正动作的维持或甚至放大(理想情况下如果伪影被完全补偿或抵消则为非衰减的,即权重或衰减因子接近“1”),而低补偿程度或能力导致衰减或完全消除,其权重因子的绝对值小于1或甚至接近零。
换言之,补偿程度越高,各自的中心图像点处的校正动作的放大越高,并且补偿程度越低,各自的中心图像点处的校正动作的衰减越高。
图像伪影是由这样的局部像素模式定义的图像特征,所述局部像素模式预期还在校正图像中显示,尽管在相反的像素或体素强度中,使得该伪影能够在将校正图像自适应地应用到初始图像之后被补偿或抵消。
所有MAR算法都能够按照这样的方式被制定出:即存在遭受金属伪影的初始图像,而算法创建了将被添加到初始图像的校正器图像。所述装置因此可以用于不明确遵循此形式而是直接创建经校正的图像的MAR算法。
在一个实施例中,计算经校正得到图像与原始图像之间的差异,并且将该差异用作校正器图像。
所述装置实现了这样的理念:将校正图像“局部”应用于其补偿伪影的区域中但避免其将创建新伪影的应用。
所述装置允许通过使用补偿程度的概念在这两种情况之间进行区分。根据一个实施例,权重是使用针对由每个中心初始图像点的邻域对形成的组合图像信息的熵度量来计算的。
根据一个备选实施例,计算邻域对中的每个中的图像信息之间的统计相关系数。
根据一个实施例,所述控制器被配置为调整针对校正器图像邻域的先前设置的默认尺寸,直到校正器图像邻域中的图像信息的熵超过预定阈值。
动态调整(以像素数或体素数量来测量的)邻域尺寸允许牵制算法的运行时间,因为运行时间会随邻域尺寸增大而增加。选择阈值熵值允许平衡运行时间的大小:由于小邻域中缺少图像信息或结构,对于过小的邻域尺寸,熵成本函数可能会变得相当平坦。然而,过大的邻域在计算上是过量的,这进一步阻碍了存在于小区域中的伪影的适当校正。
所述装置通过考虑邻域中的图像结构来选择合适的邻域尺寸,其中,“结构”是由熵函数测量的,其中,图像的熵优选被定义为在该邻域中的灰度值的归一化直方图的熵。如果校正图像在邻域中是相当平坦的,则像素灰度值直方图是高度峰化的,并具有高的熵。使用熵作为结构度量,这意味着大邻域应该被选择。另一方面,如果校正图像在邻域中具有细小条纹,则直方图显示出几个峰并具有较低的熵。因此,较小的邻域能够被选择。
应当理解,所述装置也可以用于不同于MAR的伪影校正算法。伪影不一定由“金属”部分造成,而可能由任何其他高辐射衰减部分引起。所述装置可以用于任何图像伪影校正器,不管是什么导致了图像伪影,或不管现有伪影校正器潜在的具体算法是什么。
在一个备选的装置外实施例中,所述控制器控制在远程伪影控制器处产生最终经校正的图像。
定义
放大校正动作包括将校正动作维持为在单个给定的图像点由MAR所提供的校正动作,在这种情况中权重至少等于单位量。
衰减校正动作包括消除在单个给定的图像点处的校正权重,在这种情况下权重大约为零。
图像应被广泛地解释为保存数值数据项目的至少2维阵列、矩阵或类似的数据结构,每个数值数据项目都可由至少二维的坐标i,j寻址。
图像信息或特征是由横跨构成图像平面中的补片或区域的像素的具体像素值分布给定的具体像素值模式。
附图说明
现在将参考以下附图描述本发明的示范性实施例,在附图中:
图1示意性示出根据本发明的一个实施例的用于图像校正的装置的框图;
图2示意性示出根据本发明的一个实施例的图1的装置的操作;
图3示出根据本发明的一个实施例的校正图像的方法的流程图。
具体实施方式
在图1的右侧示出用于针对伪影校正图像的装置的框图。
所述装置包括输入单元或接口器件105、处理单元110和输出单元或接口器件150。
输入105被配置为访问数据库系统180以从其中检索具有图像伪影的初始数字图像。
数据库180可以被布置成医疗设施中的PACS并且保存医学图像数据,例如2D计算机断层摄影(CT)投影图像或重建的2D或3D横截面图像(通常也被称为切片图像或“切片”)。这些图像具有适当的数字格式,例如DICOM格式。
所述装置被布置成经由接口器件105和适当的计算机通信网络连接到数据库180。
在图1的左侧示出了实施已知的金属伪影减少(MAR)算法的伪影校正器模块MAR190。图1示意性示出的该装置的实施例被布置成现有MAR190系统的“附加件”。MAR190可以访问与该图像相关的额外数据,例如初始图像所基于的投影数据,或者访问关于高度衰减部分的现有信息,例如整形外科植入物的CAD信息。
然而,在另一个实施例中,所述装置可以包括本地MAR作为处理单元110的部件。
所述装置的处理单元110包括校正采样器120、比较器130和伪影校正控制器140。
在一个实施例中,装置部件105、150、120、130和140运行为处理单元110上的软件例程。然而,在备选实施例中也可以预见到所述装置的分布式架构,其中,这些组件被连接在适当的通信网络中。所述部件还可以被布置成专用FPGA或硬连线的独立芯片。所述部件可以在适当的科学计算平台(例如或)中被编程并且之后可以被编译成例如维持在库中且当被处理器110调用时被链接的C++或C例程。
从广义上讲,用于针对伪影校正图像的该装置经由输入单元105连接到数据库180并从那里检索具有图像伪影的初始图像。该图像伪影可能是在拍摄对象的CT图像时由高密度部分(例如存在于对象中的金属部分)导致的结果。
该图像可以是从利用CT系统采集的一组投影数据重建的切片图像。
感兴趣对象可以是人体或动物体的一部分,并且高密度部分可以包括金属部分或其他高密度部分,例如内嵌在患者的软组织中的骨骼结构。图像伪影(例如条纹、失真和阴影)可能由于该部分对CT运行中所用的辐射施加了比周围软组织的平均衰减更高的衰减而存在于所采集的初始图像中。
在读入用于图像校正的初始图像后,校正采样器120将初始图像转发到MAR190并请求MAR190处理该初始图像以输出经MAR校正的样本图像,所述经MAR校正的样本图像之后在校正采样器120处被接收回。
之后初始图像及其经MAR校正的样本图像被转发到比较器130。比较器130比较这两幅图像并产生记录MAR190对初始图像的校正动作的校正器图像。
校正器图像表示对初始图像的MAR校正的针对每个图像点i,j的校正动作。之后校正器图像被转发到伪影校正控制器140以通过使用校正器图像并结合初始图像来直接获得最终经校正图像。
之后控制器140将校正动作重新应用于初始图像作为每幅校正器图像,但这次使用权重来修改初始图像中的每个图像点的校正动作。个体权重或权重因子由控制器140计算作为参数。在初始图像的图像平面中的每个待校正图像点周围,由控制器定义补片或邻域。对于这样定义的每个初始图像邻域,由控制器定义在校正器图像中的对应邻域。对于这样定义的邻域对中的每个,控制器140计算出数值权重。每个权重测量初始邻域中的图像信息能由校正器图像邻域中的图像信息补偿的程度。
之后控制器140在校正器图像中查找先前由MAR190估计的将要应用于任何给定点的校正动作。之后所估计的校正动作由各自的权重加权,并且之后这样重新加权的校正动作由控制器140应用到初始图像中的各自的图像点。通过针对初始图像中的每个图像点或选择的图像点以这种方式继续进行,建立最终经校正的图像。权重或参数由控制器140根据不同的原理来计算,并且每个实施例将在下面标题为“操作”的部分更详细地进行解释。
之后控制器140将最终经校正的图像传递到输出接口器件150。之后最终经校正的图像能够通过通信网络被调度到数据库180或者能够被转发到图像绘制器并被绘制以用于在屏幕上观看。
在另一个“装置外”实施方案中,控制器140指示MAR单元重新校正初始图像,但这次在MAR190处的校正操作是由伪影校正控制器140使用如上所述的校正器图像来控制的。在本实施例中,控制器140在运行时间期间与MAR190的算法通过接口连接,并指示将在MAR190的外来算法中使用的权重。在本实施例中,控制器140配备有适当编程的API来实现交互,并且在MAR190处的图像校正算法可能必须被适当地调整以实现控制器140在MAR190处的控制功能。
在装置外实施例中,一旦在控制器140的控制下在MAR190处产生最终经校正的图像,则之后最终经校正的图像被反馈给控制器140,控制器140进而将装置外经校正的图像传送到输出接口器件150。在又一实施例中,MAR190可以向数据库180或其他输出最终经校正的图像并且将其释放以进行进一步的后处理。
操作
初始图像μij被形成为包括i行和j列的矩阵。在第i行和第j列的矩阵条目是表示像素/体素的灰度值水平的数值。初始图像μij中的每个伪影被形成为由那些像素和体素的区域限定的独特图像特征。以下将仅参考像素进行描述,但应当理解,下面的内容同样适用于体素。
在MAR190处产生并且由校正采样器120请求的样本校正图像σij由与该图像μij数量相同的行和列形成,但因为在MAR190处经历了校正动作而一般具有不同的像素值条目。
比较器130被配置为基于初始图像μij与样本图像σij来生成校正器图像cij。在一个实施例中,比较器130形成初始图像μij与样本图像σij之间的逐像素的差,这导致校正器图像cij具有与初始图像μij或样本图像σij相同数量的行和列。这种差值则表示在该图像点i,j处的校正动作。之后校正器图像cij连同初始图像μij被转发到控制器140以用于处理。
控制器140被配置为使用校正器图像cij并且将重复的校正自适应地应用于初始图像μij。控制器140被配置为将校正图像cij仅局部应用于初始图像μij,亦即,应用被限制到校正器图像cij的校正动作将补偿伪影的区域,但在初始图像μij的将创建新伪影的其他区域中避免或充分衰减校正动作的应用。
由控制器140实施的算法允许区分将在下面更详细解释的这两种情况。
图2提供控制器140的操作的图解说明。
在图2的左侧,初始图像补片Ωmn包括示意性示为暗椭圆形的伪影。对应的校正器图像补片Ωij包含如图中以亮虚线椭圆示出的匹配或相似的图像结构或信息,其与初始图像补片Ωmn中的结构在形状上匹配,但以相反的灰度值示出。因此针对暗椭圆形伪影,MAR190正确地估计了校正动作,因为校正器图像cij在被应用于补片Ωmn中的图像点时将完全补偿暗椭圆形伪影。补片Ωmn中的图像信息因此具有高补偿程度。因此,控制器140计算大于单位值的参数A,以至少维持甚至放大由校正器图像cij针对该补片Ωmn中的图像点所记录的校正动作。
补片Ωkl示出了相反的情况。在这里,MAR190错误地估计了校正动作,因为它在应用时将会在完全没有伪影的地方(如图3中的左下方的“干净”补片Ωkl中所示)引入黑暗的矩形伪影(在图3中的右下方示出)。因此,控制器140针对补片Ωkl计算出非常低的补偿程度,因为两个补片中的结构不匹配。因此,控制参数A以小于单位值并接近于零的数值衰减,以便消灭或消除MAR190对补片Ωkl中的图像点提出的校正动作。
根据一个实施例,控制器140读取初始图像μij以为每个点i,j定义该点i,j周围的邻域Ωij(也被称为“补片”)。该点可以被称为该补片的中心点,每个补片具有至少一个这样的中心或“种子”点。邻域Ωij定义图像平面中的子集,并且可以由具有n×n像素的高度和宽度的矩形区域给出。在一个实施例中,11×11像素的正方形被选为针对邻域Ωij的默认尺寸。
在其他实施例中,邻域Ωij是每个图像点i,j周围的具有特定半径r的圆形。
之后控制器140将该邻域Ωij映射到校正器图像cij中的对应邻域Ωij。根据一个实施例,这是通过使用与针对校正器图像cij中的对应邻域所定义的像素坐标相同的像素坐标来实现的。
在一个实施例中,为了节省CPU时间,没有针对初始图像μij中的每一个像素图像像素点定义邻域,而是局限于很可能包括伪影的区域。图像平面中的该感兴趣区域能够例如通过使用适当编程的分割器来建立。作为各自的邻域Ωij的中心点的图像点i,j被选择为充分远离图像μij的边界,以确保该邻域被良好限定并且不会延伸超出当前边框。
在一个实施例中,控制器140将权重计算为初始图像μij的补片Ωij中的像素值与校正器图像cij的补片Ωij中的像素值之间的统计相关系数。根据下面的公式计算统计相关系数:
其中,和分别表示初始图像补片和校正器图像补片中的平均像素值。
之后根据以下公式计算最终经校正的图像fij的像素值:
fij=μij-tijcij
在一个实施例中,tij的值的范围被定窗宽并转变至可由削波函数A选择的范围。
在一个实施例中,部分正弦的函数被定义:
在该公式中,A(t)表示在正相关(tij=1)的0和完美反相关(tij=-1)(即负相关)的1之间变化的控制参数,并且定义了这两个机值之间的平滑正弦转变。
截止参数t0在-1和0之间变化,并定义一个截止点,在该截止点处相关系数tij被认为是“足够负的”,以保证由校正器图像cij在该点i,j处所应用的校正动作是不变的,即未衰减的。
之后使用削波函数A(t),将最终经校正的图像像素计算为:
fij=μij+A(t)cij
在一个实施例中,通过利用分割器跟踪表示骨骼或其他高密度物质的区域来使补片之间的相关性变得更鲁棒。之后在相对于该补片计算相关系数之前将这些区域从各自的补片排除。
在另一实施例中,控制器140被配置成根据相对于每对补片的熵来计算权重。
在基于熵的实施例中,针对每个补片对的权重由下面的公式确定:
换言之,控制器140被配置为求解:
其中,H表示熵函数
并且,μΩij和cΩij表示初始图像和校正图像,每个均被限制到所选择的补片Ωij。
熵函数H能够通过如上述公式(1)所示地针对每个补片建立像素水平强度直方图来计算。针对每个补片定义像素值的间隔(“箱”)并且记录数值落入任何一个箱中的像素的数目。
由具有在各自的箱范围中的值的像素的归一化分数针对每个箱近似P的数值pi,得到期望的直方图。
之后该直方图能够用作随机变量X的概率密度函数P=(pi)的近似,补片中的每个像素值i,j被视为独立同分布的像素值随机变量X的样本。
根据一个实施例,a=-0.5并且b=1.5,但是应理解,也可以根据即时伪影校正问题的需求选择另一种合理的解空间。允许大于1的值考虑了伪影强度被低估的情况。如果伪影结构的形状已被正确估计但具相反的强度,则Aij的负值可以出现。
在校正参数的公式(1)中,在矩阵的矢量空间中形成受限图像的线性组合。标量Aij测量应该在补片Ωij中添加多少校正器图像cij的校正动作,由此加入尽可能少的图像信息,即将线性组合的图像的熵最小化。
在一个实施例中,以10亨斯菲尔德单位(HU)的增量选择直方图的箱,并且默认邻域尺寸为11×11像素。
所述装置允许用户配置箱水平和默认邻域尺寸。
在任何给定图像点i,j处求解A之后,各自的邻域中心点i,j处的校正动作之后根据由作为图像结构度量的其邻近对的熵所测量的程度来由Aij加权:
fij=μij+Aijcij
使用两个补片的组合图像的熵作为权重确保得到的最终经校正的图像具有最小量的图像特征。由此,确保了校正器的应用不会增加新的伪影,亦即,初始图像中先前没有的新图像信息或特征。
在一个实施例中,权重Aij被限制于具有a=0和b=1的单位间隔。如果校正器图像仅包含已经存在于初始图像μij中的伪影结构,则被计算为单位间隔中的最小熵的参数将假定为值1或接近1。
如果校正器图像cij只包含不存在于初始图像μij中的新图像信息(新伪影),则权重Aij将被假设为值为0或接近0。
在另一个实施例中,在贝叶斯估计器框架中对所计算的权重进行合理性检查。不可能的值被丢弃,并且例如通过根据下式确定Aij而被复位到最可能的值:
其中,p是针对加权因子的先验概率函数(也被称为“先验”),例如预期值为1的高斯分布函数。
在任一上述实施例中,校正器140还能够被配置为当选择尺寸例如选择补片的边缘宽度或直径时使用所计算的熵。在这种方法中,补片是基于校正图像cij的该补片中的信息内容而被选择的。
根据由补片包围的像素区域的熵来选择补片将导致节省CPU时间并且改善最终经校正的图像的质量。使用熵允许区分信息丰富的图像区域和几乎没有结构的平坦区域。校正图像cij中的平坦区域具有为尖峰的像素值直方图,因此该区域具有高的熵,换言之,使用熵函数H作为结构度量将导致针对该区域选择比平均值更大的邻域或补片。另一方面,如果校正图像cij中的区域含有更多的图像结构,例如条纹、阴影或奇点,则该区域的像素值直方图很可能显示出几个峰,因此将具有更高的熵。在此高熵的情况下,控制器140被配置为在该区域中选择小于平均值的补片。根据一个实施例,因此自适应地根据由补片所包围的区域的熵来选择补片。
根据一个实施例,默认补片尺寸被选择并且根据在该补片中计算的熵而被缩放。控制器140使用可配置的平均熵值来影响相对于该平均值的比例。
在其他实施例中,控制器140实施不同于熵的其他图像结构度量,例如图像值的方差或与平均值的平均绝对差。
参考图3,示出针对图像伪影校正图像的方法的流程图。
在步骤S305中,接收初始图像。
在步骤S310中,初始图像被转发到现有伪影校正器以从那里接收经校正的图像的样本。该样本校正图像是由伪影校正器实施的伪影校正器算法应用于初始图像的估计校正动作的结果。
在步骤S320中,将初始图像与校正样本图像进行比较。之后该比较的结果被记录为表示伪影校正器对初始图像中的每个初始图像点i,j的校正动作的校正器图像。
在一个任选的步骤中,默认的邻域尺寸被调整(S330),直到在先前设置的默认尺寸的校正器图像邻域中的图像信息的熵超过预定阈值。阈值和默认尺寸都是可配置的,并且调整手段一般放大邻域以确保足够量的熵。熵的合理量可以是例如0.05或0.1。
之后,在步骤S340中,各自的校正动作被以加权方式自适应地重新应用于初始图像中的多个图像点中任何一个。每个影像点用于校正动作的加权的权重涉及该图像点周围的邻域中的图像信息能由校正器图像中的对应邻域中的对应图像信息补偿的程度。所述多个图像点可以包括初始图像中的所有图像点(除了沿边界或图像框架的图像点的带)或者包括已经确立为最有可能包括图像伪影的图像点的选集。分割步骤可以用于将图像平面适当地修剪到感兴趣区域。任选地,以类似的方式,图像平面可以被修整以排除图像背景区域。
在一个实施例中,权重被计算以针对每个邻域将来自各自的初始图象邻域和对应校正器图象邻域两者的图像信息的组合的熵最小化。
在一个备选实施例中,每个邻域的权重被计算为各自的初始图象邻域中的图像信息与对应校正器图象邻域中的对应图像信息之间的统计相关系数。
在本发明的另一个示范性实施例中,提供一种计算机程序或计算机程序单元,其特征在于,适于在适当的系统上执行根据前述实施例之一的方法的方法步骤。
计算机程序单元因此可以被存储在计算机单元上,所述计算机单元也可以是本发明的实施例的一部分。该计算单元可以适于执行或诱导执行上述方法的步骤。此外,其可以适于操作上述装置的部件。该计算单元能够适于自动操作和/或执行用户的命令。计算机程序可以被加载到数据处理器的工作存储器中。数据处理器因此可以被装配以执行本发明的方法。
本发明的示范性实施例涵盖以下两者,从一开始就使用本发明的计算机程序,以及借助更新将现有程序转化成使用本发明的程序的计算机程序。
进一步地,该计算机程序单元可以能够提供所有必要的步骤来完成如上所述的方法的示范性实施例的流程。
根据本发明的另一示范性实施例,提出一种计算机可读介质,例如CD-ROM,其中,所述计算机可读介质具有存储在其上的计算机程序单元,该计算机程序单元是由前一节所描述的。
计算机程序可以被存储和/或分布在适当的介质上,所述介质例如是与其他硬件一起供应或作为其他硬件的一部分供应的光学存储介质或固态介质,但是计算机程序也可以以其他形式分布,例如经由因特网或者其他有线或无线的远程通信系统。
然而,计算机程序也可以被提供在诸如万维网的网络上,并能够从这样的网络被下载到数据处理器的工作存储器中。根据本发明的另一个示范性实施例,提供一种用于使计算机程序单元可用于下载的介质,该计算机程序单元被布置为执行根据本发明的前述实施例之一的方法。
应注意,本发明的实施例是参考不同的主题进行描述的。具体而言,一些实施例是参考方法类型权利要求描述的,而其他实施例是参考设备类型权利要求描述的。然而,本领域技术人员将会从上文以及下文的描述中总结出,除非另外指出,除属于一种类型的主题的特征的任意组合外,涉及不同主题的特征之间的任何组合也被认为在本申请中公开。例如,实施例是参考CT扫描器进行描述的,但是其也能够被容易地应用于3D X射线成像。此外,所有的特征能够被组合以提供大于这些特征的简单相加的协同效应。
尽管已经在附图和前面的描述中详细说明并描述了本发明,但这样的说明和描述应当被认为是说明性或示范性的,而不是限制性的。本发明不限于所公开的实施例。通过研究附图、说明书以及从属权利要求,本领域技术人员在实施要求保护的发明时可以理解和实现所公开实施例的其他变型。
在权利要求中,词语“包括”不排除其他元件或步骤,并且量词“一”或“一个”不排除多个。单个处理器或其他单元可以实现权利要求中记载的几个项目的功能。在互不相同的从属权利要求中记载特定措施并不指示不能有利地使用这些措施的组合。权利要求书中的任何附图标记不得被解释为对范围的限制。
Claims (13)
1.一种用于针对伪影来校正图像的装置,所述装置包括:
输入单元(105),其用于接收初始图像;
处理单元(110),其用于处理所述图像以进行校正;
输出单元(150),其用于输出所述初始图像的经校正的版本;
所述处理单元包括:
校正采样器(120),其被配置为将所述初始图像转发到伪影校正器,并从那里接收经校正的图像的样本,样本图像是由所述伪影校正器对所述初始图像应用的校正动作的结果;
比较器(130),其被配置为比较所述初始图像和经校正的样本图像以便建立表示所述伪影校正器对所述初始图像的所述校正动作的校正器图像;
伪影校正控制器(140),其被配置为将各自的校正动作以加权方式自适应地重新应用于所述初始图像中的多个图像点中的任意一个,针对每个图像点的用于对所述校正动作进行加权的权重与该图像点周围的邻域中的图像信息能由所述校正器图像中的对应邻域中的对应图像信息补偿的程度相关,
所述输出单元被配置为输出被如此自适应地重新校正的初始图像作为所述经校正的图像。
2.根据权利要求1所述的装置,其中,所述权重是依赖于邻域的并且指示补偿程度,所述权重在任何给定图像点处衰减或放大所述校正动作,所述衰减与所述补偿程度相反地变化,并且所述放大直接随所述补偿程度变化。
3.根据权利要求1或2所述的装置,其中,所述控制器(140)被配置为计算所述权重以针对每个邻域将来自各自的初始图像邻域和对应的校正器图像邻域两者的图像信息的组合的熵最小化。
4.根据权利要求1或2所述的装置,其中,所述控制器(140)被配置为根据各自的初始图像邻域中的图像信息与对应的校正器图像邻域中的对应图像信息之间的统计相关系数来针对每个邻域计算所述权重。
5.根据权利要求1-4中任一项所述的装置,其中,所述对应校正器图像邻域的尺寸被所述控制器设置为默认尺寸,所述控制器被配置为调整所述尺寸直到默认的对应校正器图像邻域中的图像信息的熵超过预定阈值。
6.一种针对伪影校正图像的方法,所述方法包括:
接收(S305)初始图像;
将所述初始图像转发(S310)到伪影校正器,并从那里接收经校正的图像的样本,样本图像是由所述伪影校正器应用到所述初始图像的校正动作的结果;
比较(S320)所述初始图像和经校正的样本图像,以便建立表示所述伪影校正器对所述初始图像的所述校正动作的校正器图像;
将各自的校正动作以加权方式自适应地重新应用于(S340)所述初始图像中的多个图像点中的任意一个,针对每个图像点的用于对所述校正动作进行加权的权重与该图像点周围的邻域中的图像信息能由所述校正器图像中的对应邻域中的对应图像信息补偿的程度相关。
7.根据权利要求6所述的方法,其中,所述权重是依赖于邻域的并且指示补偿程度,所述权重在任何给定图像点处衰减或放大所述校正动作,所述衰减与所述补偿程度相反地变化,并且所述放大直接随所述补偿程度变化。
8.根据权利要求6或7所述的方法,其中,自适应地重新应用(S340)所述校正动作的步骤包括计算所述权重以针对每个邻域将来自各自的初始图像邻域和对应的校正器图像邻域两者的图像信息的组合的熵最小化。
9.根据权利要求6或7所述的方法,其中,自适应地重新应用(S340)所述校正动作的步骤包括根据各自的初始图像邻域中的图像信息与对应的校正器图像邻域中的对应图像信息之间的统计相关系数来针对每个邻域计算所述权重。
10.根据权利要求6-9中任一项所述的方法,其包括:
在自适应地重新应用(S340)所述校正动作之前,调整(S330)尺寸直到先前设置的默认尺寸的对应校正器图像邻域中的图像信息的熵超过预定阈值。
11.一种用于针对图像伪影校正图像的医学图像处理系统,所述系统包括:
根据权利要求1-5中任一项所述的装置;
伪影校正器(140);
保存所述初始图像的数据库系统(180)。
12.一种用于控制根据权利要求1-6中任一项所述的装置的计算机程序单元,所述计算机程序单元当由处理单元执行时,适于执行根据权利要求7-10所述的方法的步骤。
13.一种计算机可读介质,其上存储有根据权利要求12所述的程序单元。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201161556859P | 2011-11-08 | 2011-11-08 | |
US61/556,859 | 2011-11-08 | ||
PCT/IB2012/056041 WO2013068887A1 (en) | 2011-11-08 | 2012-10-31 | Adaptive application of metal artifact correction algorithms |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103918004A true CN103918004A (zh) | 2014-07-09 |
CN103918004B CN103918004B (zh) | 2016-12-14 |
Family
ID=47297339
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201280054611.3A Active CN103918004B (zh) | 2011-11-08 | 2012-10-31 | 金属伪影校正算法的自适应应用 |
Country Status (8)
Country | Link |
---|---|
US (1) | US10089720B2 (zh) |
EP (1) | EP2754122B1 (zh) |
JP (1) | JP6140716B2 (zh) |
CN (1) | CN103918004B (zh) |
BR (1) | BR112014010843A8 (zh) |
IN (1) | IN2014CN03500A (zh) |
RU (1) | RU2014123282A (zh) |
WO (1) | WO2013068887A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105225208A (zh) * | 2015-10-14 | 2016-01-06 | 上海联影医疗科技有限公司 | 一种计算机断层成像金属伪影校正方法及装置 |
CN106651986A (zh) * | 2016-01-21 | 2017-05-10 | 上海联影医疗科技有限公司 | 计算机断层成像伪影校正方法 |
CN107260195A (zh) * | 2016-02-29 | 2017-10-20 | 通用电气公司 | 用于计算机断层扫描成像的伪影去除的系统和方法 |
CN108937995A (zh) * | 2017-05-03 | 2018-12-07 | 西门子医疗有限公司 | 用于生成减少伪影的ct图像数据的自适应方法 |
CN113936068A (zh) * | 2020-07-14 | 2022-01-14 | 上海联影医疗科技股份有限公司 | 伪影校正方法、装置和存储介质 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9619889B2 (en) * | 2014-10-21 | 2017-04-11 | General Electric Company | Methods and systems for normalizing contrast across multiple acquisitions |
EP3362987B1 (en) | 2015-10-14 | 2021-09-22 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image correction |
CN106296615B (zh) * | 2016-08-16 | 2017-09-29 | 广州华端科技有限公司 | Ct图像中校正金属伪影的方法和系统 |
US11151703B2 (en) | 2019-09-12 | 2021-10-19 | International Business Machines Corporation | Artifact removal in medical imaging |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080089606A1 (en) * | 2006-10-16 | 2008-04-17 | Teradyne, Inc. | Adaptive background propagation method and device therefor |
CN101777177A (zh) * | 2009-12-29 | 2010-07-14 | 上海维宏电子科技有限公司 | 基于衰减滤波的ct图像去金属伪影混合重建法 |
CN102067170A (zh) * | 2007-08-31 | 2011-05-18 | 皇家飞利浦电子股份有限公司 | 在存在金属伪影的情况下用于分割的不确定性图 |
US7991243B2 (en) * | 2005-02-03 | 2011-08-02 | Koninklijke Philips Electronics N.V. | Radial adaptive filter for metal artifact correction |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5243664A (en) * | 1991-09-16 | 1993-09-07 | Picker International, Inc. | Post-processing technique for reducing metallic clip artifacts in CT images |
US5933471A (en) | 1997-02-26 | 1999-08-03 | International Business Machines Corporation | System and method for reducing reconstruction artifacts in computed tomography images |
US7001243B1 (en) | 2003-06-27 | 2006-02-21 | Lam Research Corporation | Neural network control of chemical mechanical planarization |
JP5405048B2 (ja) * | 2007-07-10 | 2014-02-05 | 株式会社東芝 | X線コンピュータ断層撮影装置、再構成処理装置及び画像処理装置 |
EP2014232B1 (en) * | 2007-07-10 | 2021-06-30 | Toshiba Medical Systems Corporation | X-ray computed tomography apparatus, reconstruction processing apparatus and image processing apparatus |
WO2010016425A1 (ja) * | 2008-08-07 | 2010-02-11 | 株式会社 日立メディコ | X線ct画像形成方法及びそれを用いたx線ct装置 |
US8244016B2 (en) | 2009-07-20 | 2012-08-14 | Wisconsin Alumni Research Foundation | Method for suppressing streak artifacts in images produced with an x-ray imaging system |
JP2011115404A (ja) * | 2009-12-03 | 2011-06-16 | Canon Inc | X線画像合成装置、およびx線画像合成方法 |
US8768027B2 (en) | 2010-02-23 | 2014-07-01 | Carestream Health, Inc. | Method and system for cone beam computed tomography high density object artifact reduction |
-
2012
- 2012-10-31 JP JP2014539451A patent/JP6140716B2/ja active Active
- 2012-10-31 US US14/354,685 patent/US10089720B2/en active Active
- 2012-10-31 WO PCT/IB2012/056041 patent/WO2013068887A1/en active Application Filing
- 2012-10-31 RU RU2014123282/08A patent/RU2014123282A/ru not_active Application Discontinuation
- 2012-10-31 CN CN201280054611.3A patent/CN103918004B/zh active Active
- 2012-10-31 BR BR112014010843A patent/BR112014010843A8/pt not_active IP Right Cessation
- 2012-10-31 IN IN3500CHN2014 patent/IN2014CN03500A/en unknown
- 2012-10-31 EP EP12797985.4A patent/EP2754122B1/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7991243B2 (en) * | 2005-02-03 | 2011-08-02 | Koninklijke Philips Electronics N.V. | Radial adaptive filter for metal artifact correction |
US20080089606A1 (en) * | 2006-10-16 | 2008-04-17 | Teradyne, Inc. | Adaptive background propagation method and device therefor |
CN102067170A (zh) * | 2007-08-31 | 2011-05-18 | 皇家飞利浦电子股份有限公司 | 在存在金属伪影的情况下用于分割的不确定性图 |
CN101777177A (zh) * | 2009-12-29 | 2010-07-14 | 上海维宏电子科技有限公司 | 基于衰减滤波的ct图像去金属伪影混合重建法 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105225208A (zh) * | 2015-10-14 | 2016-01-06 | 上海联影医疗科技有限公司 | 一种计算机断层成像金属伪影校正方法及装置 |
CN105225208B (zh) * | 2015-10-14 | 2017-10-03 | 上海联影医疗科技有限公司 | 一种计算机断层成像金属伪影校正方法及装置 |
CN106651986A (zh) * | 2016-01-21 | 2017-05-10 | 上海联影医疗科技有限公司 | 计算机断层成像伪影校正方法 |
CN106651984A (zh) * | 2016-01-21 | 2017-05-10 | 上海联影医疗科技有限公司 | 计算机断层成像伪影校正方法及装置 |
CN106683143B (zh) * | 2016-01-21 | 2021-03-16 | 上海联影医疗科技股份有限公司 | 图像金属伪影校正方法 |
CN106651984B (zh) * | 2016-01-21 | 2021-05-07 | 上海联影医疗科技股份有限公司 | 计算机断层成像伪影校正方法及装置 |
CN107260195A (zh) * | 2016-02-29 | 2017-10-20 | 通用电气公司 | 用于计算机断层扫描成像的伪影去除的系统和方法 |
CN108937995A (zh) * | 2017-05-03 | 2018-12-07 | 西门子医疗有限公司 | 用于生成减少伪影的ct图像数据的自适应方法 |
CN113936068A (zh) * | 2020-07-14 | 2022-01-14 | 上海联影医疗科技股份有限公司 | 伪影校正方法、装置和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
JP2014532506A (ja) | 2014-12-08 |
BR112014010843A2 (pt) | 2017-06-13 |
US20140286558A1 (en) | 2014-09-25 |
RU2014123282A (ru) | 2015-12-20 |
BR112014010843A8 (pt) | 2017-06-20 |
EP2754122A1 (en) | 2014-07-16 |
WO2013068887A1 (en) | 2013-05-16 |
US10089720B2 (en) | 2018-10-02 |
JP6140716B2 (ja) | 2017-05-31 |
EP2754122B1 (en) | 2019-01-23 |
CN103918004B (zh) | 2016-12-14 |
IN2014CN03500A (zh) | 2015-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103918004B (zh) | 金属伪影校正算法的自适应应用 | |
US10292672B2 (en) | Radiographic image processing device, method, and recording medium | |
JP6700273B2 (ja) | 雑音モデルに基づくマルチスケール雑音低減による線量増加のシミュレーション | |
US8233586B1 (en) | Iterative reduction of artifacts in computed tomography images using forward projection and an edge-preserving blur filter | |
EP2504811B1 (en) | Enhanced image data/dose reduction | |
US10395397B2 (en) | Metal artifacts reduction for cone beam CT | |
JP2007014759A (ja) | Pet/ct撮像における画像ベースのアーティファクト低減方法 | |
US20120148156A1 (en) | Enhanced contrast for scatter compensation in x-ray imaging | |
CN110533738B (zh) | 重建数据处理方法、装置、医学成像系统及存储介质 | |
US11010960B2 (en) | Method for enhanced display of image slices from 3-D volume image | |
US20210027430A1 (en) | Image processing apparatus, image processing method, and x-ray ct apparatus | |
CN111353958A (zh) | 图像处理方法、装置及系统 | |
US7680352B2 (en) | Processing method, image processing system and computer program | |
KR102414094B1 (ko) | Ct영상의 금속 아티팩트 보정 방법 및 장치 | |
KR102297972B1 (ko) | 전역 변형 저감화 기법을 이용한 저선량 콘빔 전산화 단층 촬영 시스템 | |
CN113077391A (zh) | 校正扫描图像的方法和装置以及图像扫描系统 | |
CN116563166A (zh) | 一种图像增强方法、装置、存储介质和设备 | |
KR102473009B1 (ko) | 비국부 전역 변형 저감화 기법을 이용한 저선량 콘빔 전산화 단층 촬영 시스템 | |
CN117252943A (zh) | 一种基于概率扩散模型(ddmp)的金属伪影去除方法 | |
Pan et al. | Combining recursive Retinex and eye adaptation for high dynamic range compression |
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 |