确定迭代影像重建中的像元映射的方法

xiaoxiao2020-7-22  31

确定迭代影像重建中的像元映射的方法
【专利摘要】一种使用像元方法迭代重建包含噪声的信号的方法从用以在迭代中更新影像的变量中确定像元映射。更新变量是基于一个最优评价函数的,并且在该迭代期间将该更新变量平滑化。在该迭代结束时还可以使用在该迭代期间确定的该像元映射任选地进一步将所更新的影像平滑化。
【专利说明】确定迭代影像重建中的像元映射的方法
[0001]相关申请
[0002]本申请要求在2011年4月15日提交的美国临时申请号61/476,012的优先权,该申请的披露内容通过引用以其全部结合于此。
发明领域
[0003]本发明涉及一种用于信号重建和增强的方法,并且更确切地说,涉及一种自适应地减少输入对象中的噪声的方法。
_4] 发明背景
[0005]对包含在一个信号中的数据的最优提取需要消除多个缺陷例如噪声以及多个仪器限制。试图进行最优化提取和重建数据的主要领域是影像增强领域。即使在可以将多个仪器制造得实质上无噪声时,与有限的空间、频谱、或时间分辨率相关的多个仪器签名仍然存在。在此点上,需要影像重建来消除这些仪器签名。影像增强的多个应用和多个噪声源以及会负面地影响数据提取的多个其他因素覆盖一个很大的范围,包括:天文观察和行星探测,其中多个源可以是微弱的并且大气干扰引入噪声和失真;军事和安全监视,其中光可以是多个目标的较低且快速的移动,从而导致对比度差且模糊;医学成像,其经常清晰度不够;以及视频影像,其中传输和仪器限制以及对实时响应的需要会负面地影响影像清晰度和细节。
[0006]已经开发了数字影像处理来提供对由多种传感器收集的模糊且有噪声的数据的高质量、鲁棒性重建。该领域存在是因为创建产生未被测量噪声污染的清晰度一致的多个图像的多个成像仪器是不可能的。然而,从自现实世界的仪器获得的不理想的数据数学上重建底图像是可能的,这样使得可以将现存但隐藏在数据中的信息提取得具有较少模糊和噪声。许多此类方法利用如下过程:构建一个预测模型且将其与这些数据比较来评估该模型的数据拟合的准确度。
[0007]贯穿此书面描述,“数据”是指任何所测量的量值,通过影像重建过程从该所测量的量值估计出一个未知的“影像”。术语“影像”表示要么估计的解要么产生观察到的数据的真实的底图像。该讨论通常清楚地表明适用的上下文;在可能的歧义的情况下,“影像模型”用以表示估计的解。注意,数据和影像不需要是相似的,并且甚至可以具有不同的维度,例如,断层重建试图从投射的2D数据确定一个3D影像。影像的一个替代术语是“对象”,其传达了模型可以比影像更概括的思想。这两个术语此后用作同义词。
[0008]统计学家长期以来试图限制用以拟合数据的多个参数,以便提高该拟合的解释和预测准确度。标准的技术是:子集选择,其中这些参数中的一些被发现是不重要的并且被从该拟合中排除(例如,米勒(Mi 11 er ) 2002 );以及岭回归,其中这些参数的值通过向用于该拟合中的评价函数添加一个正则项而受到约束(例如,吉洪诺夫(Tikhonoν ), 1963 )。蒂施莱尼(Tibshirani) (1996)将这两种方法组合在一种技术中,该技术称为最小绝对缩减和选择算子(LASSO)。
[0009]限制一个拟合中的参数的数目的需要对欠定的或不充分地确定的问题而言是必要的,其中参数的数目大于数据点的数目或与其相当。如已经由最小二乘(高斯(Gauss),1809)和最大似然(费希尔(Fisher),1912,1922)的统计方法的创始人强调的,这些方法仅仅在渐近极限内有效,其中数据点的数目大大地超过拟合的参数的数目。远离此渐近极限,噪声被拟合为信号,并且该拟合失去其解释和预测能力。
[0010]佩特(Puetter)等人(2005)回顾了当前使用中的许多重建算法,包括迭代影像重建方法,其将多个影像模型迭代地拟合成数据。即使在其较早终止来避免过度拟合时,被设计来收敛到最大似然解的许多现有技术迭代方案缓慢地收敛。收敛可以通过使用评价函数相对于变量的二阶偏导数的海赛矩阵而更快速地达成(海赛(Hesse),1876)。不幸地,在处理在影像重建过程中经常遇到的大范围问题时,例如,在该影像包含具有显著发射的大量的像素的情况下,即使使用此方法,海赛矩阵也太大而不能被计算。在此类情况下,矩阵元素可以达到数万亿。这样大小的矩阵完全不能被现代计算机处理,或者甚至存储在存储器中。
[0011]限制欠定的或不充分地确定的问题中的参数的数目(如现今所谓的,稀疏性)的希望是基于更抽象的最小复杂性的概念(索洛莫诺夫(Solomonoff), 1964 ;柯尔莫戈洛夫(Kolmogorov), 1965 ;蔡廷(Chaitin), 1966),其追溯到提倡假说简约性的奥卡姆的威廉的中世纪工作。简单来说,在其他事物相等的情况下,一个较简单的解释比一个较复杂的解释更好。
[0012]例如通过向评价函数添加一个Itl范数正则项来找到一个稀疏的解是一个N-P困难问题,其中计算工作量比许多参数中的任何多项式增加得更迅速。这已经导致用一个I1范数来替换该Itl范数(陈、多纳霍(Donoho)和桑德斯(Saunders),1999),这样使得该拟合相对于这些参数的最优化变成一个可解的凸问题。坎德丝(Cand6s)、龙伯格(Romberg)和陶(2004)更进一步并且示出了如何在一个不连贯的条件下随机地减少该拟合所需的数据量,一种称为压缩感测的技术。多纳霍(2006)示出了在一个类似的不连贯的条件下在用于参数化的基函数当中最小I1范数解也是最稀疏的解。
[0013]I1范数方法的缺点有两部分。首先,许多重要的问题完全不满足不连贯条件并且不适用于I1范数方法。其次,即使在其不满足该不连贯条件时,大范围问题也需要过度的计算工作量。因此,尽管这些问题是凸问题且原则上可解,其实际上不能应用于具有数百万个或更多参数的当今问题,这是一个同样使传统的统计方法苦恼的问题。多纳霍等人(2006)讨论了如何将随机性更有效地应用于这些大范围问题而不使用一个I1范数。
[0014]像元方法是一种获得基于像素的数据的多个复杂性最小的解的有效技术,包括大范围问题,无需严格的稀疏性并且没有利用I1范数的方法所需的不连贯条件(参见例如碧娜(Pina)和佩特,1993 ;佩特和亚希勒(Yahil), 1999 ;佩特、格斯内尔(Gosnell)和亚希勒,2005 ;以及美国专利 5,912,993,6, 353,688,6, 490,374,6, 895,125,6, 993,204、7,863,574、7,928,727、8,014,580、8,026,846、8,058,601,8, 058,625、8,086,011、8,090, 179,8, 094, 898,8, 103, 487,这些专利的披露内容通过引用结合于此)。因此,像元方法对于大范围欠定的或不充分地确定的反问题例如影像重建或频谱分析而言是有用的。最小复杂性是通过在每个像素位置处通过一个内核库当中的最宽的内核自适应地平滑化而达成的,如此使得用此以及所有较窄的内核的平滑化提供对考虑中的像素的数据足迹的充分拟合。指明在每个像素处使用哪个内核的映射叫做像元映射。[0015]在其当前形式中,一个像元重建由三个步骤组成。首先,其重建一个没有任何像元约束的“伪影像”。其次,此伪影像用以确定像元映射。第三,通过由该像元映射引导的一个受约束的重建获得最终影像。可以将步骤二和三重复许多次,但实际上这通常是不必要的,条件是在第一步骤中获得了一个合理的伪影像。参见佩特和亚希勒(1999)、佩特等人(2005)以及美国专利号 5,912,993,6, 353,688,6, 490,374,6, 895,125,6, 993,204、7,863,574,7, 928,727,8, 014,580,8, 058,601,8, 058,625,8, 086,011,8, 090,179、8,094,898以及8,103,487来获得对像元方法以及其应用的更完整的讨论。
[0016]图2图示了具有一个成像检测器210和一个像元重建单元220的一个通用成像系统200。重建是基于一种像元方法的,该方法使用一个像元映射P,其与一个像元重建算法230互动。该像元方法是指一种通过向每个对象点指派一个形状或容量作为像元平滑化的基础来将对象空间中的每个点(此后,一个“对象点”)平滑化的方法。该对象空间是定义影像重建结果的空间并且对应于使用成像系统200所成像的一个域。(应注意,“影像空间”是与“对象空间”同义的一个术语,并且此后可交换地使用这两个术语。)通过用成像检测器210测量的多个数据点来给出一个对应的数据空间。
[0017]该像元方法提供对来自数据空间中的一个所测量的数据集d的对象空间中的一个影像对象I的高质量重建。作为一种空间上自适应的重建方法,该像元方法将一个数据激发的平滑化操作应用于每个对象点。在这样做时,该像元方法在向每个对象点指派一个像元内核函数时使用最小复杂性的原则,这是平滑化操作的基础。在像元重建单元220内,像元映射P定义将这些像元内核函数中的哪一个指派给这些对象点中的每一个。
[0018]在成像系统200中,成像检测器210检测所测量的数据集d并且将其传递给像元重建单元220。像元重建单元220使用多个特别地适配的像元重建算法230来将所获得的数据集d重建成一个影像对象I。在这样做时,像元重建算法230使用一个系统矩阵H来描述成像系统200的多个性质,并且通过调整数据模型估计一个迭代地改进的影像对象,该数据模型是影像对象I的基础。使用熟知的渲染技术将影像对象I例如显示在一个显示器240 上。
[0019]针对每个对象点,像元映射P提供一个基于最小复杂性方法确定的像元内核函数。此像元内核函数用于一个应用于对象空间中的像元平滑化操作中。
[0020]该像元方法还可以允许以超分辨率形式的压缩感测,从而利用影像的非负性以及最小复杂性来重建一个影像,该影像具有比用来获得数据的那些像素更精细的像素。因为非负性和最小复杂性的附加条件(例如,佩特等人,2005),这并没有因为奈奎斯特(Nyquist) (1928)和香农(Shanon) (1949)而违反取样定理。在数据中被截短的超过衍射极限的空间频率可以在该影像中类似地重建。
[0021]图3图示了像元方法的一个示例性流程。像元平滑化顺序地应用于一个标准的重
建算法。
[0022]使用一个标准的重建算法将输入影像拟合成一个所测量的数据集d (步骤300)。根据像元内核算子K的上文讨论的使用,该影像的所产生的估计叫做一个伪影像。使用该伪影像以及所测量的数据集d确定像元映射P (步骤310)。该伪影像还是像元平滑化操作(步骤320)的最初对象。在该像元平滑化操作(步骤320)期间,该伪影像的每个对象点通过一个像元内核函数平滑化。(在现有的像元方法的一些变体中,像元映射还可以在每个迭代中通过从更新的影像计算来更新。)
[0023]迭代的影像重建方法将多个影像模型迭代地拟合成所测量的数据,并且因此最小化噪声对最终影像的影响。一个重建算法的结果是根据该算法的多个规则拟合成所测量的数据集d的一个近似的影像。
[0024]在像元方法中,一个近似的影像可以用作像元平滑化的一个输入对象,以用于像元重建以及用于确定像元映射。
[0025]像元方法包括在对象空间中的每个点处对多个最广的可能像元内核函数的搜寻,这些像元内核函数一起支持一个对象向所测量的数据集d的一个充分拟合。具体地,该像元映射向每个对象点指派一个特定的像元内核函数。
[0026]首先计算一个伪影像,接着可以从该伪影像确定像元映射的步骤可能存在多个弊端。此过程需要更多的计算并且承受着将多个伪像引入该伪影像中的风险,这可以使该像元映射以及因此最终的影像重建产生偏差。此外,在从对象空间向数据空间的变换是非局部的情况下(巴特纳格尔(Bhatnagar)和康威尔(Cornwell), 2004),像元映射的确定效果较差。例如,在干涉法以及磁共振成像过程中,数据是影像(加上噪声)的多个傅立叶变换,其中每个傅立叶波(该影像的基函数)在整个影像之上传播。另一个示例是在断层成像中的部分非局部的变换。这些数据是一个3D影像(加上噪声)的多个2D投影,一个横向于投影方向是局部的但沿着投影方向是非局部的变换。
[0027]鉴于前述内容,需要一种改进的方法来确定该像元方法内的像元映射。
[0028]发明概沭
[0029]根据此中的发明,通过在迭代期间计算以及更新像元映射避免三阶段像元方法的上述多个困难。改进的方法从一个变量确定该像元映射,该变量用以在该迭代中更新影像,即,一个“更新变量”,并且在该迭代期间将此更新变量平滑化。在该迭代结束时还经常使用在该迭代期间确定的该像元映射进一步将更新的影像平滑化。相比之下,现有的像元方法从已经更新之后的影像确定像元映射并且继续用彼像元映射将该影像平滑化。
[0030]根据本发明,更新变量取决于重建方法,但典型地是一个评价函数的梯度或者一个乘法更新因子(例如,佩特等人,2005)。此更新变量通过多个像元内核平滑化,并且在每个影像位置处选择的内核是最宽的内核,其中该更新变量中因为平滑化的改变的平方与该更新变量的方差之间的比率小于该内核以及所有较窄的内核的一个预定的阈值。相同的内核接着通常还用以在继续进行到下一个迭代之前在彼位置处将更新的影像平滑化。一个进一步的细化允许“多个插入的内核”。
[0031 ] 在本发明的一个方面中,一种从自物理过程获得的数据集重建对象模型的方法,其中该数据集包含噪声,该方法包括:接收在一个数据空间中定义的该数据集;在一个对象空间中构建一个对象模型,其中,该对象模型包括多个对象点;设计该对象模型从对象空间向数据空间的一个变换,从而产生一个数据模型,其中,该变换对应于获得该数据集的物理过程;选择用于确定一个数据模型向该数据集的一个拟合的一个评价函数;基于该评价函数确定对象空间中的该对象模型的一个更新变量;将该更新变量平滑化以通过以下方法确定一个平滑化的更新变量:使该更新变量与多个像元内核中的每一个卷积;以及为该输入对象的每个对象点选择一个像元内核,该像元内核具有符合一个预定的最小标准的一个最大大小;通过在对应于所选的像元内核的该输入对象的每个对象点处指派多个索引来产生一个像元映射;以及基于该像元映射内的这些索引产生一个包括一个实质上降噪的对象模型的输出。在一个实施例中,使用一种共轭梯度法来确定该评价函数并且该更新变量是一个负梯度。
[0032]在本发明的另一个方面中,一种非瞬态计算机可读介质,具有附带的软件,该软件包括用于从自一个物理过程获得的一个数据集重建一个对象模型的多个指令,其中,该数据集包含噪声,并且这些指令包括:接收在一个数据空间中定义的该数据集;在一个对象空间中构建一个对象模型,其中,该对象模型包括多个对象点;设计该对象模型从对象空间向数据空间的一个变换,从而产生一个数据模型,其中,该变换对应于获得该数据集的物理过程;选择用于确定一个数据模型向该数据集的一个拟合的一个评价函数;基于该评价函数确定对象空间中的该对象模型的一个更新变量;将该更新变量平滑化以通过以下方法确定一个平滑化的更新变量:使该更新变量与多个像元内核中的每一个卷积;以及为该输入对象的每个对象点选择一个像元内核,该像元内核具有符合一个预定的最小标准的一个最大大小;以及通过在对应于所选的像元内核的该输入对象的每个对象点处指派多个索引来产生一个像元映射。在一个示例性实施例中,使用一种共轭梯度法来确定该评价函数并且该更新变量是一个负梯度。
[0033]附图简要说明
[0034]图1是可以实践本发明的一般计算环境的方块图。
[0035]图2是基于像元方法的一个现有技术成像系统的概述图。
[0036]图3是图示将现有技术方法用于像元重建的影像重建的一个示例的流程图。
[0037]图4是图示使用根据本发明的像元方法的影像重建的一个示例的流程图。
[0038]图5是图示使用像元方法的一个替代影像重建过程的一个示例的流程图。
[0039]图6示出了一个无噪声的、归一化的、质子密度(PD) SBD “真”的脑影像。
[0040]图7示出了图6的影像向数据空间的傅立叶变换的绝对值的对数(基数10)。
[0041]图8示出了具有随机(白)高斯噪声的图6的影像,该随机(白)高斯噪声添加到该影像,具有最亮的影像体素的3%的标准差。
[0042]图9示出了图8的影像向数据空间的傅立叶变换的绝对值的对数(基数10)。
[0043]图10示出了与图9相同的数据,其中高频率被截短,对应于限制到可用频率的三分之一的频率覆盖范围。图10中的数据包括向图11至图14中所示的影像的多个重建的输入。
[0044]图11示出了图10的数据的一个简单的逆傅立叶变换,揭示了两种类型的伪像:关联的(非白)噪声以及在多个锐边处的振影。
[0045]图12示出了图10的数据的一个像元重建的结果,其中像元因子为0.3。
[0046]图13示出了图10的数据的一个像元重建的结果,其中像元因子为0.5。
[0047]图14示出了图10的数据的多个像元重建的结果,其中像元因子为1.0。
[0048]发明详沭
[0049]在描述本发明的多个方面之前,提供对可以实现本发明的一个合适的计算系统环境100 (图1)的一个简要描述可能是有用的。计算系统环境100仅仅是一个合适的计算环境的一个示例,并且不旨在建议关于本发明的使用或功能性的范围的任何限制。也不应将计算环境100解释为具有与在示例性运行环境100中图示的多个部件中的任一者或组合有关的任何从属性或要求。
[0050]本发明是用许多其他一般目的或特殊目的计算系统环境或配置可运行的。可能适合于与本发明一起使用的多个熟知的计算系统、环境、和/或配置的多个示例包括(但不限于)个人计算机、服务器计算机、手持式或膝上型装置、多处理器系统、基于微处理器的系统、机顶盒、可编程的消费电子器件、网络PC、小型计算机、大型计算机、电话系统、包括上述系统或装置中的任一者的分布式计算环境,等等。
[0051]可以在多个计算机可执行的指令例如由一个计算机执行的多个程序模块的一般上下文中描述本发明。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、分量、数据结构等。本领域技术人员可以将此中的描述和/或多个图实现为多个计算机可执行的指令,这些指令可以在下文讨论的任何形式的计算机可读介质上体现。
[0052]本发明还可以在多个分布式计算环境中实践,其中通过多个远程处理装置执行多个任务,这些远程处理装置通过一个通信网络链接。在一个分布式计算环境中,多个程序模块可以定位在本地和远程计算机存储介质两者中,这些计算机存储介质包括存储器存储装置。
[0053]参见图1,一种用于实现本发明的示例性系统包括呈一个计算机110形式的一个通用计算装置。计算机110的多个部件可以包括(但不限于)一个处理单元120、一个系统存储器130、以及一个系统总线121,该系统总线将包括该系统存储器的不同系统部件耦合至处理单元120。系统总线121可以是若干类型的总线结构中的任何一种,包括存储器总线或存储器控制器、外围总线、以及使用多种总线架构中的任一种的局部总线。通过举例,并非限制,此类架构包括工业标准架构(ISA)总线、微通道架构(MCA)总线、增强型ISA(EISA)总线、视频电子标准协会(VESA)内部总线、以及外围部件互连(PCI)总线,也称为夹层(Mezzanine)总线。
[0054]计算机110典型地包括多种计算机可读介质。计算机可读介质可以是可以由计算机110访问的任何可用的介质并且包括易失性和非易失性介质、可移除和不可移除介质两者。通过举例,并非限制,计算机可读介质可以包括计算机存储介质和通信介质。计算机存储介质包括易失性和非易失性、可移除和不可移除介质两者,这些介质以用于存储信息例如计算机可读指令、数据结构、程序模块或其他数据的任何方法或技术实现。计算机存储介质包括(但不限于)RAM、ROM、EEPR0M、闪存或其他存储器技术、CD-ROM、数字多用磁盘(DVD)或其他光盘存储器、磁盒、磁带、磁盘存储器或其他磁性存储装置、或可用以存储期望的信息并且可以由计算机110访问的任何其他介质。通信介质典型地体现计算机可读指令、数据结构、程序模块或一个已调数据信号例如载波或其他传送机构中的其他数据,并且包括任何信号输送介质。术语“已调制数据信号”意味着具有其特征集中的一个或多个或被改变的方式为将该信号中的信息编码的一个信号。通过举例,并非限制,通信介质包括有线介质,例如有线网络或直接有线连接,以及无线介质,例如声学、RF、红外线以及其他无线介质。上述中的任一者的组合也应该包括在计算机可读介质的范围内。
[0055]系统存储器130包括呈易失性和/或非易失性存储器形式的计算机存储介质,例如只读存储器(R0M)131以及随机访问存储器(RAM)132。基本输入/输出系统133 (BIOS)典型地存储在R0M131中,该基本输入/输出系统包含多个基本例程,这些基本例程帮助例如在启动期间在计算机Iio内的多个元素之间传递信息。RAM132典型地包含处理单元120立即可访问和/或目前由该处理单元对其运行的数据和/或程序模块。通过举例,并非限制,图1图示了一个操作系统134、多个应用程序135、多个其他程序模块136、以及程序数据137。
[0056]计算机110还可以包括其他可移除/不可移除易失性/非易失性计算机存储介质。仅通过举例,图1图示了从不可移除、非易失性磁性介质读取或向该介质写入的一个硬磁盘驱动器141、从一个可移除、非易失性磁盘152读取或向该磁盘写入的一个磁盘驱动器151、以及从一个可移除、非易失性光盘156例如一个CD ROM或其他光学介质读取或向该光盘写入的一个光盘驱动器155。可用于示例性运行环境中的其他可移除/不可移除、易失性/非易失性计算机存储介质包括(但不限于)盒式磁带、闪存卡、数字多用磁盘、数字视频磁带、固态RAM、固态R0M,等等。硬磁盘驱动器141典型地通过一个不可移除的存储器接口例如接口 140连接至系统总线121,并且磁盘驱动器151和光盘驱动器155典型地通过一个可移除的存储器接口例如接口 150连接至系统总线121。
[0057]上文讨论并且在图1中图示的这些驱动器及其相关联的计算机存储介质为计算机110提供对计算机可读指令、数据结构、程序模块以及其他数据的存储。在图1中,例如,硬磁盘驱动器141图示为存储操作系统144、多个应用程序145、多个其他程序模块146以及程序数据147。注意,这些部件可以与操作系统134、多个应用程序135、多个其他程序模块136以及程序数据137要么相同要么不同。操作系统144、多个应用程序145、多个其他程序模块146、以及程序数据147在这里被给予不同的数字来说明在最低限度其为不同的副本。
[0058]使用者可以通过多个输入装置向计算机110输入命令和信息,这些输入装置例如键盘162、麦克风163 (其还可以表示通过电话提供的输入)、以及指点装置161,例如鼠标、轨迹球或触摸板。其他输入装置(未示出)可以包括操纵杆、游戏键盘、圆盘式卫星天线、扫描仪,等等。这些和其他输入装置经常通过耦合至该系统总线的使用者输入接口 160连接至处理单元120,但是可以通过其他接口和总线结构例如并行端口、游戏杆端口或通用串行总线(USB)连接。监视器191或其他类型的显示装置也经由一个接口例如视频接口 190连接至系统总线121。除该监视器之外,多个计算机还可以包括其他外围输出装置,例如扬声器197和打印机196,这些装置可以通过一个输出外围接口 195连接。
[0059]计算机110可以在一个联网环境中使用至一个或多个远程计算机例如一个远程计算机180的逻辑连接来运行。远程计算机180可以是个人计算机、手持式装置、服务器、路由器、网络PC、对等装置或其他公用网络节点,并且典型地包括上文相对于计算机110描述的多个元素中的许多或全部。图1中所描绘的多个逻辑连接包括一个局域网(LAN) 171以及一个广域网(WAN) 173,但是还可以包括其他网络。此类联网环境在办公室、企业范围的计算机网络、企业内部网以及互联网中是司空见惯的。
[0060]当用于一个LAN联网环境中时,计算机110通过一个网络接口或适配器170连接至LAN171。当用于一个WAN联网环境中时,计算机110典型地包括一个调制解调器172或其他装置来在WAN173例如互联网之上建立通信。可以是内部的或外部的调制解调器172可以经由使用者输入接口 160或其他适当的机构连接至系统总线121。在一个联网的环境中,相对于计算机110或其多个部分描绘的多个程序模块可以存储在远程存储器存储装置中。通过举例,并非限制,图1将多个远程应用程序185图示为驻存在远程计算机180上。应了解,所示的多个网络连接是示例性的并且可以使用在多个计算机之间建立一个通信链路的其他装置。
[0061]在此申请内,“像元”用以指示一个术语、一种方法、一个对象等是指该像元方法,SP,是指多个对象点特定的形状在将如在美国专利5,912,993,6, 895,125以及其他相关专利中所描述的一个影像对象平滑化时的使用。例如,多个所指派的形状由多个像元内核函数定义,并且像元映射P存储关于这些像元内核函数中的哪一个/哪些被指派给每个对象点的信息。
[0062]数据的拟合是一个噪声逆问题的解。这些数据被模型化成
[0063]d=H(l)+ ε , (I)
[0064]其中d是包含所测量的数据的一个阵列,I是包含多个拟合参数的一个阵列,H是从参数空间到数据空间的变换函数,并且ε是包括向数据的噪声贡献的一个阵列。出于图示的目的,在Eq.(I)中使用适合于影像重建的一个符号,即,d为像素数据并且I为重建的影像。然而,此示例不旨在是限制性的,因为此中的发明具有向可以使用该像元方法解的任何噪声逆问题的广泛的应用。
[0065]针对许多应用,变换H是线性的,并且Eq.⑴可以写成
[0066]d=HI+ ε , (2)
[0067]其中d和I被当作向量,并且H被当作一个矩阵,称为系统矩阵。以下讨论局限于Eq.(2)的线性情况,但是其可以通过线性化推广到非线性问题,即,考虑H(I)是近于线性的一个有限的影像域。
[0068]注意,事实上,d和I可以是多维的并且甚至不需要具有相同的维度(例如,在断层成像中)。但是多维的像素总是可以安排成长向量,并且人们可以接着将H当作一个矩阵。
[0069]数据模型
[0070]m=HI=d- ε (3)
[0071]是数据的无噪声的信号部分。其还是以下数据的期望值:
[0072]m=E (d), (4)
[0073]因为在没有通用性损失的情况下可以为以下所有数据点将噪声的期望值同样地设置成零:
[0074]Ε( ε ) ξ 0ο (5)
[0075]噪声的协方差矩阵一般地是非零的并且可以是位置(像素位置)和/或数据模型值的一个函数
[0076]V ( ε ) =V (m)。(6)
[0077]正常地,多个数据点是独立的,所以V是一个对角矩阵,但是还可以存在关联的数据,在此情况下,V具有多个非零的、非对角的元素。
[0078]给定数据d,系统矩阵H以及噪声ε的统计模型,影像重建是针对影像I解Eq.(2)——或者,更一般地,Eq.(I)——的逆问题。给定以下数据模型,影像重建是通过最优化数据的一个评价函数而求解的:
[0079]L ε L (d I m) =L (d IHI)。 (7)
[0080]评价函数的选择取决于噪声统计。评价函数通常被看作是对数似然函数(LLF),其与最大似然(例如,佩特等人,2005)相关。
[0081]针对高斯噪声,LLF减少至X2
[0082]
【权利要求】
1.一种从自物理过程获得的数据集重建对象模型的方法,其中,该数据集包含噪声,该方法包括: 接收在一个数据空间中定义的该数据集; 在一个对象空间中构建一个对象模型,其中,该对象模型包括多个对象点; 设计该对象模型从对象空间向数据空间的一个变换,从而产生一个数据模型,其中,该变换对应于获得该数据集的物理过程; 选择用于确定一个数据模型向该数据集的一个拟合的一个评价函数; 基于该评价函数确定对象空间中的该对象模型的一个更新变量; 将该更新变量平滑化以通过以下方法确定一个平滑化的更新变量: 使该更新变量与多个像元内核中的每一个卷积;以及 为该输入对象的每个对象点选择一个像元内核,该像元内核具有符合一个预定的最小标准的一个最大大小; 通过在对应于所选的像元内核的每个对象点处指派多个索引来产生一个像元映射;以及 基于该像元映射内的这些 索引产生一个包括一个实质上降噪的对象模型的输出。
2.如权利要求1所述的方法,进一步包括根据该像元映射通过使每个对象点处的对象模型与该所选的像元内核卷积来将该对象模型平滑化。
3.如权利要求1所述的方法,其中,使用一种共轭梯度法来确定该评价函数。
4.如权利要求3所述的方法,其中,该更新变量是一个负梯度。
5.如权利要求3所述的方法,进一步包括添加一个预处理算子。
6.如权利要求1所述的方法,进一步包括: 迭代地将该更新变量平滑化并且基于该平滑化的更新变量更新该像元映射以及该对象模型,直到将该评价函数最优化为止。
7.如权利要求1所述的方法,进一步包括在对该更新变量进行卷积的步骤之后,减去该更新变量来确定一个更新变量改变。
8.如权利要求7所述的方法,其中,该预定的最小标准是基于该更新变量改变的平方与该更新变量的一个方差之间的比率的。
9.如权利要求1所述的方法,其中,该数据集包括结合一个分析过程产生的干涉数据,该分析过程选自由以下各项组成的群组:孔径合成、干涉合成孔径雷达(IfSAR或InSAR)、合成孔径雷达(SAR)和逆合成孔径雷达(ISAR)、合成孔径声纳、波束成形、以及合成孔径磁场定位法。
10.如权利要求1所述的方法,其中,该数据集包括磁共振成像数据。
11.如权利要求1所述的方法,其中,该数据集包括计算机或发射断层成像数据。
12.如权利要求1所述的方法,其中,该数据集包括频谱分析数据。
13.如权利要求1所述的方法,其中,该数据集包括来自不同仪器的多个输入信号,并且该降噪的对象模型包括对应于这些输入信号的一个混合物的一个信号输出。
14.一种非瞬态计算机可读介质,其上具有附带的软件,该软件包括用于从自物理过程获得的数据集重建对象模型的多个指令,其中,该数据集包含噪声,该方法包括: 接收在一个数据空间中定义的该数据集;在一个对象空间中构建一个对象模型,其中,该对象模型包括多个对象点; 设计该对象模型从对象空间向数据空间的一个变换,从而产生一个数据模型,其中,该变换对应于获得该数据集的物理过程; 选择用于确定一个数据模型向该数据集的一个拟合的一个评价函数; 基于该评价函数确定对象空间中的该对象模型的一个更新变量; 将该更新变量平滑化以通过以下方法确定一个平滑化的更新变量: 使该更新变量与多个像元内核中的每一个卷积;以及 为该输入对象的每个对象点选择一个像元内核,该像元内核具有符合一个预定的最小标准的一个最大大小;以及 通过在对应于所选的像元内核的该输入对象的每个对象点处指派多个索引来产生一个像元映射。
15.如权利要求14所述的非瞬态计算机可读介质,进一步包括根据该像元映射通过使每个对象点处的对象模型与该所选的像元内核卷积来将该对象模型平滑化。
16.如权利要求14所述的非瞬态计算机可读介质,其中,使用一种共轭梯度法来将该评价函数最优化。
17.如权利要求16所述的非瞬态计算机可读介质,其中,该更新变量是一个负梯度。
18.如权利要求17所述的非瞬态计算机可读介质,进一步包括添加一个预处理算子。
19.如权利要求14所述的非瞬态计算机可读介质,进一步包括: 迭代地将该更新变量平滑化并且基于该迭代地平滑化的更新变量更新该像元映射以及该对象模型,直到将该评价函数最优化为止。
20.如权利要求19所述的非瞬态计算机可读介质,进一步包括在对该更新变量进行卷积的步骤之后,减去该更新变量来确定一个更新变量改变。
21.如权利要求20所述的非瞬态计算机可读介质,其中,该预定的最小标准是基于该更新变量改变的平方与该更新变量的一个方差之间的比率的。
【文档编号】G06K9/40GK103608825SQ201280025138
【公开日】2014年2月26日 申请日期:2012年4月16日 优先权日:2011年4月15日
【发明者】埃莫斯·亚希尔 申请人:伊麦格莱肯有限责任公司

最新回复(0)