专利名称：用于在根据成像数据的参数重构中数据驱动的自适应网格生成时的用户交互的系统和方法
技术领域：
本发明整体上涉及一种用于在利用自适应网栅来根据图像数据进行直 接迭代重构中的用户交互的系统和方法。
背景技术：
可以利用从多种(分子)成像模式中聚集的数据来重构用于描述在人体 内的示踪化学剂的浓度(例如，浓度的动态特性)的模型参数，所述成像模式
例如为正电子发射断层造影(PET)和单光子发射计算机断层造影(SPECT)扫 描器。在"逐个体素"的基础上描述这种参数，其中，体素是在研究对象 上所施加的三维(3D)栅格内的小体积元素。栅格内的体素的大小决定了可以 用来估计模型参数分布的空间精确度或分辨率。几种重构方法(例如，根据 列表模式数据的直接重构和后验估计最大化方法)允许不规则的体素栅格 (grid),例如，包含有各种形状和大小的体素的栅格。对这种体素栅格布局 的最佳选择仍然是悬而未决的问题。
重构图像数据的技术现状包括在具有分辨率(例如，体素大小)局部变化 的"不规则的"体素栅格上进行重构，而且包括在感兴趣区域中具有较高 分辨率的固态栅格，这是在重构前、在例如，初步重构或以另一个模态(例 如，CT扫描)为基础的重构上用户指示的。技术现状还包括用于图像重构的 内容自适应网格生成，其中，在较高空间变化的区域内自动提高分辨率。
对图像的网格建模包括将图像域划分成许多不重叠的(一般为多边形的) 碎片，其称为网格元素，(在此，如图1所示，使用了三角形)；然后通过根 据元素的网格节点进行内插，来在每一个元素上确定图像函数。 一个节点 对于图像的贡献限于附着到该节点的那些元素的程度。利用网格模型，可 以有策略地在包括重要特征的区域内最为密集地放置网格节点，结果生成 相比于体素表现更紧凑的图像表现。通过非常细的体素栅格来实现的高分辨率需要过长的计算时间。由较 粗的体素栅格来实现的低分辨率会导致空间信息的损失和精确度较低的系 统输出(例如，参数图)。此外，包括例如感兴趣区域、空间变化、足量统计 数据的可用性和计算时间的可用性或需求在内的各个方面都影响了在速度 与高分辨率之间的最佳折衷。在考虑中的感兴趣区域影响了对较高分辨率 的需求。研究对象的特定区域也许比其它区域更为重要，由此需要较高的 分辨率。而且，在"不那么"感兴趣区域(例如，背景)中的较高分辨率不产 生附加信息值，但却仍然会使重构处理慢了下来。关于空间变化，模型参 数会以一个区域内的逐个体素间的强空间变化为特征，然而在其它区域中 变化较为缓慢。在模型参数变化相对缓慢的区域中，只需要有限的分辨率， 反之，通过细化的网栅来对强模型参数变化区域进行最佳建模。每一个体素的模型参数估计都依赖于大量充足的与该特定体素有关的 事件(例如，检测器测量数据)。如果由于体素太小而存在太少事件，例如，具有较差的信噪比(SNR)，则由此会产生较差的估计值。由此可见，充足的 统计数据的可用性影响了在速度与高分辨率之间的最佳折衷。当然，在可以使用无限的时间时，如果在较小的体素大小上检测到大 量充足的事件，则整个栅格中的较小的体素大小将总会是优选的。然而， 医师的可利用时间是有限的，而且并不喜欢等待结果。因此速度与精确度 的相对重要性就直接影响了对分辨率的选择。因此，将要求在重构处理中进行控制以最小化计算时间，从而避免不 必要的等待时间并在临床上可接受的重构时间的范围内确保感兴趣区域的 最大分辨率。发明内容本发明涉及一种方法，用于在根据成像数据进行模型参数重构中、数 据驱动的自适应网格生成时的具有用户交互的迭代重构。在一个示范性实施例中，所述方法包括读取来自用户的(先验的(110)和即时的(115))输入，并且在每一次迭代之后检查所重构参数的收敛性(130)。在每一次迭代之后， 基于当前网栅和所期望的迭代次数，来估计所需要的计算时间(BO)，然后 更新所述网栅(140)。在所述重构的过程中，显示对所述所重构参数和调节后的网栅(mesh grid)的即时呈现(170)，并使所述重构的下一次迭代以所述调 节后的网栅为基础(145)。
在另一示范性实施例中，公开了一种系统，用于在根据成像数据进行 模型参数重构中的数据驱动的自适应网格生成时的具有用户交互的迭代重 构。所述系统包括重构器，配置为在每一次迭代之后检査所重构参数的 收敛性，并且在每一次迭代之后，基于当前网栅和所期望的迭代次数来估 计所需要的计算时间。用户界面被配置为接受用户输入以供所述重构器读 取；显示装置17在所述重构更新所述网栅14的过程中，显示对所述所重 构参数和调节后的网栅的即时呈现，其中，所述重构的下一次迭代是以所 述调节后的网栅为基础的。
在再另一示范性实施例中，公开了计算机软件产品，用于在根据成像 数据进行模型参数重构中的数据驱动的自适应网格生成时的具有用户交互 的迭代重构。所述产品包括计算机可读介质，其中存储了程序指令，所述 指令在由计算机读取时使所述计算机读取来自用户的(先验的(110)和即时 的(115))输入，并在每一次迭代之后检査所重构参数的收敛性(130)。在每一 次迭代之后，在每一次迭代之后，基于当前网栅和所期望的迭代次数来估 计所需要的计算时间(130)，然后更新所述网栅(140)。然后，所述计算机指 示在重构的过程中在显示装置上显示对所述所重构参数和调节后的网栅的 即时呈现(170)，并使所述重构的下一次迭代以所述调节后的网栅为基础 (145)。
与公开的系统和方法有关的其他特征、功能和优势将从以下的详细说 明中显而易见，尤其是在结合其附图进行评述时。
为了帮助本领域普通技术人员制造和使用该公开的系统和方法，参考
了附图，其中
图1描绘了根据本发明示范性实施例所述的图形用户界面的平面图， 用户利用该图形用户界面来选择感兴趣区域、最大可允许计算时间和其它 重构选项；以及
图2是表示根据本发明示范性实施例所述的利用图1中的图形用户界面的重构处理的流程图。
具体实施方式
正如在此所描述的，本发明有利地提供了一种使用自适应网栅的直接 迭代重构方法。由对感兴趣区域的先验指示和该重构处理的状态来确定栅 格布局。在较早的迭代中，参数估计仍然是粗糙的，并且以低分辨率栅格 为特征。分辨率随着每一次迭代而提高，并且在参数估计开始收敛时达到 其峰值。还根据每个体素的可用数据来确定栅格布局。在不活跃的区域中， 对体素进行合并(例如，汇聚)，以形成在其中对于每个体素都具有更佳的信 噪比的粗糙栅格。还利用所重构参数的空间变化来确定栅格布局。用更细 的体素栅格来覆盖高度变化的区域。还通过对最大可允许计算时间的选择 来确定栅格布局。在重构开始之前，用户定义一个最大计算时间。在每一 次迭代之后，估计剩余的计算时间，并根据是否会(例如，容易地)超过或达 到该可允许计算时间，来调节栅格分辨率。正如以下更充分讨论的，还利 用其它用户交互来确定栅格布局。参照图1，先验的和即时的用户交互通过图形用户界面(GUI)10来进行。 在重构开始之前，提示用户在先前进行的重构数据中指定感兴趣区域，该 重构数据可以例如是利用不同的模式所获得的，例如计算机断层噪声(CT)。 用户可以利用整体显示在图形用户界面10右下方的导航窗12，来指定感兴 趣区域，如图1所示。可选地，可以使用鼠标和/或具有快捷键的键盘(未显 示)。用户还可以设置最大可允许计算时间和其它重构选项。在(例如，即时的)重构处理过程中，用户观看当前使用的(3D)栅格14 和在16上的逐个体素进行亮度编码的所重构模型的参数值，并定义一个重 构参数的图15。用户在显示了当前估计的所重构参数的图15和当前使用的 网栅14的显示器17上査看它们。利用箭头按钮18来引导3D光标19穿过 栅格14，并利用縮放按钮20来调整栅格14的大小，其中用户可以选择要 提高或降低分辨率的区域。还可以利用相应的按钮22来围绕着三个轴旋转 整个图像。整体显示在GUI 10左侧的多个按钮用于整体上由24表示的全 局操作，日志消息窗口 26显示了与用户动作有关的重构进程和反馈。另外， 日志消息窗口 26提供了与所估计参数的收敛性、所估计的剩余时间和当前分辨率有关的信息，这也是基于用户动作的。用户还可以分别利用按钮28 和30来选择提高或降低整体分辨率以及提高或降低重构参数处理的速度。 现在参照图2，将参照以100总体表示的流程图来描述该参数重构处理。 在块110中，用户输入列表模式数据、感兴趣区域的定义/初始分割、最大 重构时间周期和初始模式参数。在块120中，将这些用户输入传送到重构 器，用以进行最初的迭代。在重构器的每一次迭代之后，在块130中，检 查所重构参数的收敛性，估计所需的计算时间，并在块130中读取用户在 块115中的任何即时输入。接下来，在块140中更新网栅。在块140中， 以所重构参数(en)的局部变化、仍然需要的计算时间与剩余的可允许计算时 间之间的比率(ETA/Tmax)、以及来自用户的命令(用户输入)为基础，来更新 网栅。重构器的下一次迭代是以指向块120的直线145所表示的经调节的 网格为基础的。在块170中，用户通过显示器17来接收与虚线160所表示的当前参数 估计(6 )和虚线150所表示的网栅14相关的信息。如上所述，用户可以在 块IIO和115中主动改变网栅14。在完成重构或重构循环已经收敛时，如 直线175所示，在块180中输出重构图像。本技术领域技术人员将清楚， 尽管显示器17显示为GUI 10的一部分，但是该显示器17可以是一个独立 显示器，并且与如图1所示的、位于显示器17的左下方的用户输入按钮彼 此分离。对所需计算时间的每一次估计都取决于当前网栅和所期望的迭代次 数。后者易于通过所谓的一次通过(one-pass)算法或者迭代次数固定的其它算法来计算，其中，在所述一次通过算法中，仅观看一次所有的数据。明 确取决于对所重构参数的估计的收敛性的算法必须基于收敛性统计来估计 剩余的迭代次数。在本领域中已知的重构算法在每一个事件后、在整组事件中的一个事 件子集之后、或者在一次包括了所有事件的迭代之后，产生用于模型参数 的更新的en。为了确保正确的用户交互，必须选择每次迭代所使用的事件 数量(例如，用于每一次参数更新)，使其足够小以便为用户提供以合理时间 间隔进行交互的机会。应当注意，用户交互不仅发生在重构过程中，还会发生在该重构结束之后。在重构循环已经收敛时，在计算机存储器内存储该系统的状态(例如， 当前所使用的体素栅格和所估计的模型参数值)。图形用户界面10还允许用 户以该重构图像为基础来提高在感兴趣区域内的空间分辨率，随后重构循 环可以继续进行。利用了该特征的实例将会进行最初的"快速重构"，提高 在感兴趣区域(可能是患者特有的)内的分辨率，然后允许系统继续进行"主 要的重构"。有利地，本发明的各个实施例使该系统、方法和计算机软件产品的用 户能够在重构过程中可视地检査重构参数和网栅的即时显示。此外，本发 明的系统、方法和计算机软件产品通过手动调节局部和全局网栅分辨率来 促进对重构处理的用户用户交互，并使用所估计的剩余计算时间作为网格 调节中的判定因子。其它优势包括对重构数据的先验可用性的较小依赖。例如，可以在重 构数据变为可用时，对感兴趣区域的粗糙的第一指定进行即时细化。本发 明的系统、方法和计算机软件产品还在重构处理中提供更多的控制。例如， 自动的数据驱动的网格分割可以不同于人类专家的选择。虽然关于用重构 算法选择应具有高分辨率的区域的选项仍然存在，但是用户界面还是添加 了用来使人类专家了解判定处理的主动部分的选项。 一旦感兴趣的特征开 始显露在重构的数据中，就可以"更精确地"(例如，在较高分辨率下)检查 在所重构参数的图中非预期出现的感兴趣的特征。因此，就不必完成整个 重建、添加或改变感兴趣区域、以及重新运行该重构。本发明上述系统、 方法和计算机软件产品所提供的另一个优势包括用来设置最大计算时间的 选项，用于避免不必要的等待以及在可允许时间范围内确保最大分辨率。尽管已经参照本发明的示范性实施例描述了本发明的方法、系统和软 件产品，但是本发明并不局限于这些示范性实施例。而是，该公开的方法、 系统和软件产品在不脱离其精神或范围的情况下可以进行各种修改、增强 和/或变化。因此，本发明采用和包含了在所附权利要求范围内的这些修改、 增强和/或变化。
权利要求
1、一种方法，用于在根据成像数据进行模型参数重构中的数据驱动的自适应网格生成时的具有用户交互的迭代重构，所述方法包括读取来自用户的输入(110，115)；在每一次迭代之后，检查所重构参数的收敛性(130)；在每一次迭代之后，基于当前网栅和所期望的迭代次数来估计所需要的计算时间(130)；更新所述网栅(140)；在所述重构的过程中，显示对所述所重构参数和调节后的网栅的即时呈现(170)；以及使所述重构的下一次迭代以所述调节后的网栅为基础(145)。
2、 如权利要求l所述的方法，其中，更新所述网栅的步骤是基于以下 至少之一的所述所重构参数(en)的局部变化、仍然需要的计算时间与剩余 的可允许计算时间之间的比率、以及所述用户输入。
3、 如权利要求l所述的方法，其中，读取用户输入(110， 115)的步骤在以下情况中的至少一种情况下进行在重构过程中和在所述重构完成之后。
4、 如权利要求3所述的方法，其中，所述用户输入包括以下之一 感兴趣区域；所述感兴趣区域的期望分辨率； 用于完成重构的最大可允许计算时间；以及 所述最大可允许计算时间的变化。
5、 如权利要求1所述的方法，还包括在所述所重构参数已经收敛时， 在一存储器装置中存储当前使用的体素栅格和所估计的模型参数值。
6、 如权利要求1所述的方法，还包括形成初始重构图像；基于初始重构图像，提高在感兴趣区域内的空间分辨率；以及 继续进行主要的重构循环。
7、 如权利要求1所述的方法，还包括使用所估计的剩余计算时间作 为对所述调节后的网栅的网格调节中的判定因子。
8、 如权利要求l所述的方法，还包括通过手动调节局部和全局网栅 分辨率中至少一个，来促进对所述重构的即时用户交互。
9、 如权利要求l所述的方法，其中，所述所重构参数是根据多种成像 形式中聚集的数据推导的。
10、 如权利要求1所述的方法，其中，在每一次迭代之后估计所需要 的计算时间的步骤基于来自所述用户输入的对于最大可允许计算时间的输 入，来估计剩余的计算时间，并且根据预计是否能够在所述最大可允许计 算时间内完成重构，来调节所述网栅。
11、 如权利要求1所述的方法，其中，读取来自用户的输入的步骤包 括读取通过图形用户输入而输入的输入内容。
12、 一种方法，用于在根据成像数据进行模型参数重构中的数据驱动 的自适应网格生成时的具有用户交互的迭代重构，所述方法包括利用图形用户界面，在先前重构的图像数据中指示感兴趣区域(110); 通过所述图形用户界面，输入重构器完成重构的最大计算时间(110); 通过所述图形用户界面，输入初始模式参数(110); 初始化一迭代(120); 检查所重构参数的收敛性(130);针对所述最大计算时间来估计所需的剩余计算时间(130); 利用对所述图形用户界面的任何用户输入(115)以及所需的计算时间与剩余的计算时间之间的比率，来更新所述网栅(140);在所述图形用户界面上显示当前参数估计值(170);以及 当在以下情况之一时输出所述重构(180):所述重构完成以及所述重构参数已经收敛。
13、 一种系统，用于在根据成像数据进行模型参数重构中的数据驱动 的自适应网格生成时的具有用户交互的迭代重构，所述系统包括重构器，配置为在每一次迭代之后检查所重构参数的收敛性，并且在 每一次迭代之后，基于当前网栅和所期望的迭代次数来估计所需要的计算 时间；用户界面，配置为接受用户输入以供所述重构器读取； 显示装置17，配置为在所述重构更新所述网栅14的过程中，显示对所 述所重构参数和调节后的网栅的即时呈现；以及其中，所述重构的下一次迭代是以所述调节后的网栅为基础的。
14、 如权利要求13所述的系统，其中，所述用户界面和所述显示装置 17是图形用户界面10。
15、 如权利要求13所述的系统，其中，所述更新后的网栅是基于以下 至少之一的所述所重构参数(6n)的局部变化、仍然需要的计算时间与剩余 的可允许计算时间之间的比率、以及所述用户输入。
16、 如权利要求13所述的系统，其中，在以下至少一种情况中读取用 户输入在重构过程中和在所述重构完成之后。
17、 如权利要求16所述的系统，其中，所述用户输入包括以下之一 感兴趣区域；所述感兴趣区域的期望分辨率； 用于完成重构的最大可允许计算时间；以及 所述最大可允许计算时间的变化。
18、 如权利要求13所述的系统，还包括存储装置，用于在所述所重构 参数已经收敛时，存储当前使用的体素栅格和所估计的模型参数值。
19、 如权利要求13所述的系统，其中，所述重构器通过基于来自所述 用户输入的对于最大可允许计算时间的输入而估计剩余的计算时间，来在 每一次迭代之后估计所需要的计算时间，并且所述网栅是根据预计是否能 够在所述最大可允许计算时间内完成重构而进行调节的。
20、 一种计算机软件产品，用于在根据成像数据进行模型参数重构中 的数据驱动的自适应网格生成时的具有用户交互的迭代重构，所述产品包 括计算机可读介质，其中存储了程序指令，所述指令在由计算机读取时使 所述计算机能够执行以下操作读取来自用户输入设备的输入(IIO, 115);在每一次迭代之后，检查所重构参数的收敛性(130);在每一次迭代之后，基于当前网栅和所期望的迭代次数来估计所需要的计算时间(130);更新所述网栅(140);在所述重构的过程中，显示对所述所重构参数和调节后的网栅的即时呈现(170);以及使所述重构的下一次迭代以所述调节后的网栅为基础(145)。
全文摘要
公开了一种系统与方法，用于在根据成像数据进行模型参数重构中的数据驱动的自适应网格生成时的具有用户交互的迭代重构。所述方法包括读取来自用户的输入(110，115)，并且在每一次迭代之后，检查所重构参数的收敛性(130)。在每一次迭代之后，基于当前网栅和所期望的迭代次数来估计所需要的计算时间(130)，并更新所述网栅(140)。在所述重构的过程中，显示对所述所重构参数和调节后的网栅的即时呈现(170)，并且使所述重构的下一次迭代以所述调节后的网栅为基础(145)。
文档编号G06T17/40GK101322157SQ200680045286
公开日2008年12月10日 申请日期2006年11月15日 优先权日2005年12月2日
发明者A·韦伯, B·巴克, M·纳拉亚南 申请人:皇家飞利浦电子股份有限公司