书签 分享 收藏 举报 版权申诉 / 14

针对利用正则化的迭代图像重建和/或图像去噪的正则因子的自动确定.pdf

  • 上传人:000****221
  • 文档编号:4615895
  • 上传时间:2018-10-22
  • 格式:PDF
  • 页数:14
  • 大小:3MB
  • 摘要
    申请专利号:

    CN201280061612.0

    申请日:

    2012.12.04

    公开号:

    CN104025156A

    公开日:

    2014.09.03

    当前法律状态:

    授权

    有效性:

    有权

    法律详情:

    授权|||实质审查的生效IPC(主分类):G06T 11/00申请日:20121204|||公开

    IPC分类号:

    G06T11/00

    主分类号:

    G06T11/00

    申请人:

    皇家飞利浦有限公司

    发明人:

    F·贝格纳; B·J·布伦德尔; T·克勒

    地址:

    荷兰艾恩德霍芬

    优先权:

    2011.12.13 US 61/569,835

    专利代理机构:

    永新专利商标代理有限公司 72002

    代理人:

    王英;刘炳胜

    PDF完整版下载: PDF下载
    内容摘要

    一种处理部件(122),其基于利用正则化的迭代重建算法和/或去噪算法处理图像。所述处理部件包括设定点确定器(224),所述设定点确定器基于指示感兴趣图像质量的质量变量(228),来确定在预定质量上限和质量下限(226)之间的质量设定点(216)。所述处理部件还包括比较器(214),所述比较器在每次处理迭代,将当前生成的图像的质量度量与所述质量设定点进行比较,并生成指示所述质量度量与所述质量设定点之间的差的差值。所述处理部件还包括正则因子更新器(220),所述正则因子更新器响应于所述差值指示所述质量度量在关于所述质量设定点的预定范围之外,而基于所述正则因子的当前值(222)和至少所述质量度量,来生成针对下一处理迭代的更新的正则因子。

    权利要求书

    权利要求书
    1.  一种处理部件(122),其基于利用正则化的迭代重建算法和/或去噪算法来处理图像,所述处理部件包括:
    设定点确定器(224),其基于指示感兴趣图像质量的质量变量(228),来确定介于预定质量上限和质量下限(226)之间的质量设定点(216);
    比较器(214),其在每次处理迭代,将当前生成的图像的质量度量与所述质量设定点进行比较,并生成指示所述质量度量与所述质量设定点之间的差的差值;以及
    正则因子更新器(220),其响应于所述差值指示所述质量度量处于关于所述质量设定点的预定范围之外,而基于正则因子的当前值(222)和至少所述质量度量,来生成针对下一处理迭代的更新的正则因子。

    2.  如权利要求1所述的处理部件,其中,所述质量变量限定介于所述上限和所述下限之间的相对目标水平。

    3.  如权利要求2所述的处理部件,其中,所述相对目标水平对应于感兴趣的目标图像噪声或伪影水平降低。

    4.  如权利要求1至3中的任一项所述的处理部件,其中,所述上限对应于在正则化所述图像之前的所述图像的质量度量,并且所述下限对应于在完全正则化所述图像之后的所述图像的质量度量。

    5.  如权利要求1至4中的任一项所述的处理部件,其中,所述下限被近似而无需完全正则化所述图像。

    6.  如权利要求1至5中的任一项所述的处理部件,其中,所述上限和所述下限通过所述上限被归一化。

    7.  如权利要求1至6中的任一项所述的处理部件,其中,所述正则因 子更新器通过以下中的一种或多种方式来生成所述更新的正则因子:用所述正则因子的所述当前值乘以是所述质量度量的函数的更新项;或通过向所述当前正则因子增加是所述质量度量的函数的更新项。

    8.  如权利要求1至7中的任一项所述的处理部件,其中,所述正则因子更新器基于比例-积分-微分控制来生成所述更新的正则因子。

    9.  如权利要求1至8中的任一项所述的处理部件,其中,迭代的所述质量度量指示该迭代所生成的图像的图像噪声水平。

    10.  如权利要求1至8中的任一项所述的处理部件,其中,所述迭代的所述质量度量指示该迭代所生成的图像的图像噪声水平的标准偏差。

    11.  如权利要求1至8中的任一项所述的处理部件,其中,迭代的所述质量度量指示该迭代所生成的图像的伪影水平。

    12.  如权利要求1至11中的任一项所述的处理部件,还包括:
    去噪器(302),其在每次迭代,基于针对相对应的迭代的所述更新的正则因子,来对所述图像进行去噪。

    13.  如权利要求1至11中的任一项所述的处理部件,还包括:
    图像更新器(202),其在每次迭代,基于针对相对应的迭代的所述更新的正则因子,来更新所述图像。

    14.  如权利要求1至13中的任一项所述的处理部件,其中,所述正则因子更新器基于所述质量度量的预定行为,来生成所述更新的正则因子。

    15.  一种方法,包括:
    在第一处理迭代期间,使用采用正则化和初始正则因子的算法,来处理初始图像并生成第一图像;
    确定所述第一图像的第一质量度量;
    将所述第一图像的所述第一质量度量与质量设定点进行比较,并确定所述第一质量度量与质量设定点之间的差值;并且
    响应于所述第一质量度量处于关于所述质量设定点的预定范围之外,而针对所述图像的后续处理迭代生成更新的正则因子,
    其中,所述更新的正则因子是基于所述初始正则因子和所述第一质量度量而生成的,并且所述更新的正则因子被用于在所述后续处理迭代期间处理所述第一图像并生成后续图像。

    16.  如权利要求15所述的方法,其中,所述质量设定点基于介于质量下限与质量上限之间的用户定义的目标百分数。

    17.  如权利要求16所述的方法,其中,所述下限对应于在正则化之前的所述图像的质量值,并且所述上限对应于在完全正则化之后的所述图像的质量值。

    18.  如权利要求16至17中的任一项所述的方法,其中,所述质量度量在每次迭代,指示每个各自处理迭代所生成的图像的图像噪声水平或伪影水平。

    19.  如权利要求15至18中的任一项所述的方法,还包括:
    在每次处理迭代,生成更新的正则因子,直到所述第一质量度量处于关于所述质量设定点的所述预定范围之内。

    20.  一种处理部件(122),包括:
    使用采用正则化的迭代算法处理图像并生成图像的部件(202、302),其中,所述迭代算法至少包括乘以正则因子的正则化项;
    分析器(208),其在每次处理迭代之后,针对各自生成的图像生成质量度量;以及
    正则因子确定器(212),其基于当前处理迭代生成的图像的质量度量 和针对所生成的图像的预定质量设定点(216),来自动生成针对所述部件的下一处理迭代的更新的正则因子。

    说明书

    说明书针对利用正则化的迭代图像重建和/或图像去噪的正则因子的自动确定
    技术领域
    以下大体涉及正则化的迭代图像重建和/或图像去噪,并且具体应用于计算机断层摄影(CT),并且也可用于其他成像模态。
    背景技术
    CT扫描器包括安装在可旋转机架上的X射线管,所述可旋转机架关于z轴绕检查区域旋转。探测器阵列对向一角度弧,跨检查区域与所述X射线管相对。所述X射线管发射穿过所述检查区域的辐射。所述探测器阵列探测穿过所述检查区域的辐射,并生成指示其的投影数据。重建器使用迭代或非迭代重建算法,重建所述投影数据,生成指示所述检查区域的图像。所述图像的噪声水平取决于各种成像参数,并且可以通过去噪算法得以降低。
    已通过具有数据保真项和正则项的代价函数,来将利用正则化的迭代图像重建和/或图像去噪算法公式化。用正则项乘以正则因子β,并然后将两个项相加,以确定总代价,所述总代价被最小化。通用代价函数可以被表达为:代价(f)=F(f)+β·R(f),其中f表示图像,F(f)表示数据保真项,R(f)表示正则项,并且β为自由参数并且表示控制所述正则化的强度的正则因子。对于正则化的迭代重建,所述数据保真项将正投影的图像与所测量的投影进行比较,任选地包括所述投影的噪声模型。为了去噪,所述数据保真项可以为与所述初始的、有噪声图像的差。类似地,可以包括噪声模型。
    一般地,手动确定所述正则因子β。然而,β倾向于依赖于图像和数据集。因此,必须针对每个不同的图像和/或数据集手动确定β。这已通过用不同的β(平行地或连续地)运行几次迭代,并选择满足预定标准的β,来实现。遗憾的是,该过程可能是计算和时间密集型的。例如,假定典型的迭代耗费二(2)或三(3)小时,并且可能在到达合乎期望的最终图像之 前要运行十次(10)或更多次迭代,确定针对特定图像的β可能耗费二十(20)或更多小时。
    发明内容
    本文描述的各方面解决上述问题以及其他问题。
    在一方面中,一种处理部件基于利用正则化的迭代重建算法和/或去噪算法,来处理图像。所述处理部件包括设定点确定器,设定点确定器基于指示感兴趣图像质量的质量变量,来确定在预定质量上限和质量下限之间的质量设定点。所述处理部件还包括比较器,所述比较器在每次处理迭代,将当前生成的图像的质量度量与所述质量设定点进行比较,并生成指示所述质量度量与所述质量设定点之间的差的差值。所述处理部件还包括正则因子更新器,正则因子更新器响应于所述差值指示所述质量度量在关于所述质量设定点的预定范围之外,而基于所述正则因子的当前值和至少所述质量度量,来生成针对下一处理迭代的更新的正则因子。
    在另一方面中,一种方法,包括:在第一处理迭代期间,使用采用正则化和初始正则因子的算法,来处理初始图像并生成第一图像。所述方法还包括:确定所述第一图像的第一质量度量。所述方法还包括:将所述第一图像的所述第一质量度量与质量设定点进行比较并确定它们之间的差值。所述方法还包括:响应于所述第一质量度量在关于所述质量设定点的预定范围之外,而针对所述图像的后续处理迭代,来生成更新的正则因子。所述更新的正则因子是基于所述初始正则因子和所述第一质量度量生成的,并且所述更新的正则因子被用于在所述后续处理迭代期间处理所述第一图像并生成后续图像。
    在另一方面中,一种处理部件,包括使用采用正则化的迭代算法处理图像并生成图像的部件,其中,所述迭代算法至少包括乘以正则因子的正则化项。所述处理部件还包括分析器,所述分析器在每次处理迭代之后,针对各自生成的图像生成质量度量。所述处理部件还包括正则因子确定器,所述正则因子确定器基于当前处理迭代生成的图像的质量度量和针对所生成的图像的预定质量设定点,自动生成针对所述部件的下一处理迭代的更新的正则因子。
    附图说明
    本发明可以采取各种部件与各部件的布置以及各个步骤与各步骤的安排的形式。附图仅出于图示优选的实施例的目的并且不旨在被解释为限制本发明。
    图1示意性地图示了示范性成像系统。
    图2示意性地图示了可以与图1的成像系统一起使用的示范性重建器。
    图3示意性地图示了可以与图1的成像系统一起使用的示范性图像去噪器。
    图4以图形方式图示了自动确定的正则因子到对应于预定质量度量的值的收敛。
    图5以图形方式图示了针对所述预定质量度量的所述收敛的期望行为。
    图6图示了根据图2或图3中的至少一个的示范性方法。
    具体实施方式
    以下大体涉及利用正则化的迭代重建和/或图像去噪。如在下文更详细描述,用本文中描述的方法,所述正则因子β被自动确定,并且即使所述正则因子β可以随图像和/或数据集的不同而变化,用户也仅需要提供(例如在0%(即没有正则化,或不正则化的)与100%(即,完全正则化)之间的)正则化的相对水平。在一个实例中,这减轻了由用户通过对每个图像和/或数据集试错来手动确定合适的正则因子,并因此减少了达到合适的正则因子β和最终图像的时间量。
    图1图示了示范性成像系统100,例如计算机断层摄影(CT)扫描器。成像系统100包括大体固定机架102和旋转机架104,旋转机架104由固定机架102可旋转地支撑并关于z轴绕检查区域106旋转。
    辐射源110,例如X射线管,由旋转机架104可旋转地支撑。辐射源110与旋转机架104一起旋转,并且发射穿过检查区域106的辐射。一维或二维辐射敏感探测器阵列112对向一角度弧,跨检查区域106与辐射源110相对。探测器阵列112探测穿过检查区域106的辐射,并生成指示其的投影数据。
    对象支撑体118,例如躺椅,在检查区域106中支撑对象或受试者,例如人或动物受试者。通用计算系统或计算机充当操作者控制台120。控制台120包括人类可读的输出设备(例如监视器)和输入设备(例如键盘、鼠标等)。驻留在控制台120上的软件允许操作者经由图形用户界面(GUI)或以其他方式与扫描器100交互和/或操作扫描器100。
    重建器114重建所述投影数据,并生成指示其的图像数据。如在下文更详细描述,在一个实例中,重建器114采用利用正则化的迭代重建算法,并基于针对当前迭代生成的图像的图像质量度量和预定的感兴趣正则水平,来确定针对每次迭代的正则因子β,所述预定的感兴趣正则水平可以为缺省和/或用户定义的水平。在一个实例中,自动确定所述正则因子。
    任选的图像去噪器116去噪由重建器114和/或其他重建器生成的图像,所述其他重建器包括采用常规的滤波反投影重建算法和/或其他基于正则化或非正则化的重建算法的重建器。如在下文更详细描述,在一个实例中,去噪器116,类似地,基于针对当前迭代生成的当前生成的图像的图像质量度量和预定的感兴趣正则水平,自动确定针对后续迭代的正则因子β,所述预定的感兴趣正则水平可以为缺省和/或用户定义的水平。
    重建器114和图像去噪器116可以为相同部件(例如在图1中所述的处理部件122)的部分,或为单独且不同的部件。要认识到,针对所述正则化的迭代重建算法或所述图像去噪中的任一个或两者,自动确定正则因子β,可以提升正则化的迭代重建算法和/或图像去噪的可用性,因为用户不需要知晓β的绝对值,而是仅提供正则化的相对水平。
    重建器114和/或图像去噪器116经由处理器(或微处理器、控制器等等)运行被编码或嵌入在计算机可读存储介质(例如物理存储器)上的一个或多个计算机可读指令来实施。这样的处理器可以为控制台120和/或其他计算设备(例如专用计算机和/或其他计算设备)的部分。所述处理器也可以运行由载波、信号或其他非计算机可读存储介质(例如暂时性介质)承载的至少一个计算机可读指令。
    图2图示了重建器114的范例。
    图像更新器202接收初始图像,并使用利用正则化的迭代重建算法204处理所述图像。在Fessler的“Statistical image reconstruction methods for  transmission tomography”(Handbook of Medical Imaging,第2卷:Medical Image Processing and Analysis,1.70页,SPIE,贝灵翰姆,2000年)中描述了这样的算法的范例。例如,该引文讨论了用于通过最大化以下形式的惩罚似然目标函数,来确定对真实衰减图的估计的方法:其中Φ为所述目标函数,并且R(μ)为粗糙度惩罚。
    在所图示的实施例中,针对第一重建迭代的正则因子被设定为初始值β206,其可以为预定的缺省值和/或用户定义的,并且可以为零或非零的值。在这样的处理之后,图像更新器202生成图像,并且在后续重建迭代期间类似地处理所生成的图像。在所图示的范例中,图像更新器202使用β在每次重建迭代的不同更新值,迭代地重建最终图像。
    分析器208在每次重建迭代之后接收所生成的图像,并基于一个或多个给定质量标准(例如噪声、伪影等等),来分析所接收的图像。分析器208生成指示其的质量度量。所述质量度量可以为单个值或多个值。质量度量的范例包括,但不限于,全局或局部图像噪声标准偏差或方差。可以任选地使用其他量化标准。分析器208在每次重建迭代之后输出所述质量度量。
    β确定器212在每次重建迭代之后接收每个生成的图像的质量度量,并确定针对下一重建迭代的更新的β。β确定器212包括比较器214,其将每次重建迭代的所接收的质量度量与存储在存储器218中的预定质量设定点216进行比较,确定每次重建迭代指示所接收的质量度量与预定质量设定点216之间的误差或差的值。
    如果所述质量度量满足质量设定点216(例如,所述质量度量在关于设定点216的预定范围或容差之内),则β确定器212向图像更新器202传送一信号,其终止所述迭代重建,并且图像更新器202输出最终的重建图像。其他停止标准可以包括迭代的次数、流逝了一时间持续等。否则,β更新器220更新β的当前值222,并将所更新的β传送到图像更新器202,用于下一重建迭代。
    可以以各种方式确定设定点216。例如,设定点确定器224基于质量上限和质量下限226和质量变量228,来确定设定点216,质量变量228指示 在质量上限和质量下限226之间的水平。在一个非限制性的实例中,上限被确定为初始的未正则化的图像的所述质量度量,并且下限被确定为完全(或100%)正则化的初始图像的所述质量度量,其可以通过完全正则化所述图像来确定,或无需完全正则化所述初始图像地近似。可以通过所述上限来归一化质量上限和质量下限226。本文中也预期用于确定质量上限和质量下限的其他方法。
    质量变量228在该实施例中表示上限和下限226之间的相对目标水平(例如,被表达为百分数或其他方式)(例如在上限和下限226之间的50%)。该相对目标水平可以由用户(例如放射科医师或其他执业医师)或计算机确定或设定。
    图3图示图像去噪器116的范例。
    去噪器302接收初始(例如,未去噪的)图像,并基于去噪算法304来对所述图像进行去噪。所述初始图像可以由成像系统100或其他系统生成。在2011年10月17日递交的、受让于皇家飞利浦电子股份有限公司(荷兰埃因霍温)的题为“Low Dose CT De-Noising”的国际(PCT)申请序列号PCT/IB2011/054588中,描述了合适的正则化的去噪算法的非限制性范例,在此通过引用将其整体并入。
    去噪器302使用β的初始值306,用去噪算法304来去噪所接收的初始图像。可以如本文针对β的初始值206所述,来确定β的初始值306。在这样的处理之后,图像去噪器302生成去噪图像。如果所述去噪图像的图像质量不满足感兴趣图像质量,则去噪器302也可以进一步去噪所生成的去噪图像,如下文所述。在所图示的范例中,去噪器302以迭代的方式,使用每次迭代β的不同更新值,进一步去噪所生成的去噪图像。
    结合图2描述分析器208和β确定器212。
    亦即,分析器208在每次去噪迭代之后,接收所生成的图像,并基于一个或多个给定质量标准(例如噪声、伪影等等),来分析接收到的图像。分析器208生成指示其的质量度量。所述质量度量可以为单个值或多个值。质量度量的范围包括,但不限于,全局或局部图像噪声标准偏差或方差。可以任选地使用其他量化标准。分析器208在每次图像去噪迭代之后,输出所述质量度量。
    β确定器212在每次去噪迭代之后,接收每个生成的图像的所述质量度量,并确定针对每个下一去噪迭代的更新的β。β确定器212包括比较器214,其将所接收的每次去噪迭代的质量度量与存储在存储器218中的预定质量设定点216进行比较,确定指示每次去噪迭代所接收的质量度量与预定质量设定点216之间的误差或差的值。
    类似地,如果所述质量度量满足质量设定点216(例如,所述质量度量在关于设定点216的预定范围或容差之内),则β确定器212将信号传送到图像去噪器302,其终止去噪,并且图像去噪器302输出最终的去噪图像。类似地,其他停止标准可以包括迭代的次数、逝去了一时间持续,等等。否则,β更新器220更新β的当前值222,并且将更新的β传送到图像去噪器302,用于所述下一图像去噪迭代。
    在该范例中,所述上限被确定为初始(非去噪)图像的质量度量,并且所述下限被确定为完全(100%)去噪的初始图像的质量度量,其可以通过完全去噪所述图像来确定,或无需完全去噪所述初始图像而被近似。类似地,可以通过所述上限归一化质量上限和质量下限226。本文中也预期用于确定质量上限和质量下限226的其他方法。
    此外,可以基于剂量设定质量变量228,这允许针对给定图像质量的所述图像的剂量减少。例如,当执业医师要将质量变量设定在50%时,由于可以确定剂量与噪声之间的关系,针对给定图像质量患者剂量被减少至四分之一。在另一实例中,可以针对给定患者剂量,改善最终图像质量。在又另一范例中,改善图像质量且减小剂量。本文中也预期用于确定设定点216的其他方法。
    针对在图2和图3中描述的实施例之一或两者,在等式1中示出了基于质量设定点216,使用比例-积分-微分(PID)控制,用于在每次迭代更新β的非限制性方法:
    等式1:
    β(k)=β(k-1)·2P(e(k)+IE(k)+DΔe(k)),
    其中,β(k)为针对所述下一迭代的更新的β,β(k-1)为前一更新的β,P、I和D分别为针对PID控制的比例、积分和微分常数,e(k)表示当前迭代的质量度量与当前迭代的设定点之间的差,E(k)为直到第k次迭代的e(k) 的加合,并且Δe(k)表示在第(k-1)到第k次迭代之间,所述质量度量与所述设定点之间的所述差的改变(即,e(k)-e(k-1))。在另一实施例中,可以使用P、PI和/或其他控制来更新β。
    图4以图形方式图示了使用具有Huber惩罚的最大似然重建,在多次迭代上更新β的非限制性范例,其中,所述质量度量为所述图像的全局噪声水平。
    第一y轴402表示在上限和下限226之间的质量度量范围,其被归一化到所述上限,第二y轴404表示β的值,并且x轴406表示重建或去噪迭代的次数。针对图4,质量设定点408被设定到在归一化的上限和下限之间的50%(或0.5)。质量度量曲线410示出,所述质量度量在大致500次迭代之后,收敛到所述50%点。Β曲线412示出,一旦所述质量度量达到设定点,β的值收敛到并在大致250次时保持相当恒定。该范例是出于解释目的提供的,并且非限制性的。
    要认识到,β更新器220可以任选地应用一套更新规则,例如以缓解在所述质量度量过低和/或过冲。例如,图5以图形方式图示了针对所述质量度量的范例目标行为500。在图5中,y轴502表示在上限和下限226之间的质量范围,其被归一化到所述上限,并且x轴504表示重建或去噪迭代的次数。如在图5中所示,目标行为500包括大致线性的区域506(即,具有固定斜率的线性衰减),直到到达设定点508,以缓解设定点216的可能的过冲。本文中也预期其他目标行为。该范例行为是出于解释目的提供的,并且非限制性的。
    图6图示了根据本文描述的实施例的范例方法。
    要认识到,本文描述的所述方法中的动作的顺序不是限制性的。这样,本文也预期其他排序。此外,可以省略一个或多个动作和/或可以包括一个或多个外的动作。
    在600,获得图像,用于以采用利用正则因子的算法处理。
    在602,将所述正则因子设定到初始值。
    在604,经由正则化的迭代重建算法或去噪算法,处理所述图像,例如如本文中所述,产生经处理的图像。
    在606,确定针对所处理的图像的质量度量。
    在608,将所述质量度量与质量设定点进行比较,并确定指示所述质量度量与所述质量设定点之间的差的差值。
    在610,确定所述质量度量的绝对值是否在关于所述质量设定点的预定范围(例如,±容差)之内。如本文中描述的,可以基于质量下限与质量上限226之间的目标水平,来确定所述质量设定点。
    如果所述质量度量在所述预定范围之外,则在612,如本文中所述,更新所述正则因子,并且用所更新的正则因子重复动作604-610。如果所述质量度量在所述预定范围之内,则在614,输出最终图像。
    已参考优选的实施例描述了本发明。他人在阅读和理解前文的详细描述时,可以想到多种修改和变型。本发明旨在被解释为包括所有这样的修改和变型,只要它们落入权利要求书或其等价方案的范围之内。

    关 键  词:
    针对 利用 正则 图像 建和 因子 自动 确定
      专利查询网所有文档均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    0条评论

    还可以输入200字符

    暂无评论,赶快抢占沙发吧。

    关于本文
    本文标题:针对利用正则化的迭代图像重建和/或图像去噪的正则因子的自动确定.pdf
    链接地址:https://www.zhuanlichaxun.net/p-4615895.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    copyright@ 2017-2018 zhuanlichaxun.net网站版权所有
    经营许可证编号:粤ICP备2021068784号-1