用于产生检查对象结果图像的方法

申请号 CN200410094232.4 申请日 2004-12-08 公开(公告)号 CN1666710A 公开(公告)日 2005-09-14
申请人 西门子公司; 发明人 马丁·坦克;
摘要 本 发明 涉及一种根据检查对象截面图像数据自动产生其结果图像的方法。根据诊断问题确定截面图像数据(D)中的目标结构(Z),相应于该目标结构选择解剖学标准模型(M),其几何形状根据模型参数(tx,ty,tz,rx,ry,rz,d0)可变。该标准模型自动与目标结构匹配。将截面图像数据根据匹配的标准模型分段,根据诊断问题将检查对象的相关解剖学结构分开,其中,选择截面图像数据中的所有那些图像点,它们位于相应于相关解剖学结构匹配的标准模型和/或至少一个模型部分(T1,T2,T3,T4,T5)的轮廓内或与其最大偏差特定的差值。最后将相关的解剖学结构分开显示和/或为过后的显示存储。此外还涉及相应的 图像处理 系统。
权利要求

1.一种方法,用于根据检查对象的截面图像数据自动产生该检查对象的结 果图像,其中:
-首先根据诊断问题确定截面图像数据(D)中的目标结构(Z),
-相应于该目标结构(Z)选择一个解剖学标准模型(M),其几何形状根 据模型参数(tx,ty,tz,rx,ry,rz,d0)可变,
-将标准模型(M)自动地与截面图像数据(D)中的目标结构(Z)进行 匹配,
-然后将截面图像数据基于匹配的标准模型(M)进行分段,其中,关于 诊断问题将检查对象的相关解剖学结构分开,其中,选择在截面图像数据中的 所有那些图像点,它们位于相应于相关解剖学结构的匹配的标准模型(M)和/ 或至少一个模型部分(T1,T2,T3,T4,T5)的轮廓内,或者与其最大偏差一 特定的差值,
-最后将相关的解剖学结构分开进行显示和/或为以后的可视化进行存储。
2.根据权利要求1的方法,其特征在于,在匹配期间分别根据特定的偏差 函数确定在修改的标准模型(M)和目标结构(Z)之间的当前偏差值。
3.根据权利要求2的方法,其特征在于,在自动匹配方法中的这样改变所 述模型参数(tx,ty,tz,rx,ry,rz,d0),使得偏差值最小。
4.根据权利要求2或3的方法,其特征在于,在分段前自动检验标准模型 (M)与目标结构(Z)的匹配是否达到最小值,该最小值小于预先规定的阈值, 否则中断该方法以手工处理截面图像数据。
5.根据权利要求1到4之一的方法,其特征在于,自动检验检查对象的至 少一个分开的解剖学结构的标准偏差。
6.根据权利要求5的方法,其特征在于,利用所属的分开的解剖学结构图 形显示所确定的标准偏差和/或用声音信号通知操作者。
7.根据权利要求5或6的方法,其特征在于,根据确定的标准偏差对检查 对象进行自动分类。
8.根据权利要求1到7中之一的方法,其特征在于,在多个迭代步骤中借 助模型参数(tx,ty,tz,rx,ry,rz,d0)将所述标准模型(M)与截面图像数据 (D)中的目标结构(Z)进行匹配,所述模型参数根据它们对模型(M)的整 体解剖学几何形状的影响被分成层次,并且随着迭代步骤数的增加可调模型参 数(tx,ty,tz,rx,ry,rz,d0)的数目也相应于它的层次级而升高。
9.根据权利要求8的方法,其特征在于,所述模型参数(tx,ty,tz,rx, ry,rz,d0)分别分配到一个层次级。
10.根据权利要求9的方法,其特征在于,所述将模型参数(tx,ty,tz,rx, ry,rz,d0)分配到层次级是根据在模型几何中的偏差进行的,该偏差当所涉及 的模型参数(tx,ty,tz,rx,ry,rz,d0)改变特定的值时产生。
11.根据权利要求10的方法,其特征在于,对不同的层次级分配特定值域 的偏差。
12.根据权利要求1到11中之一的方法,其特征在于,作为标准模型(M) 使用在三基础上产生的表面模型。
13.根据权利要求1到12中之一的方法,其特征在于,这样将所述模型参 数分别与至少一个解剖学标记的位置相关联,使得该模型(M)对每个参数数 据组都具有有解剖学意义的几何形状。
14.根据权利要求1到13中之一的方法,其特征在于,借助轮廓分析法至 少部分地自动确定所述截面图像数据(D)中的目标结构(Z)。
15.一种计算机程序产品,其可直接加载到可编程图像处理系统的存储器 中,具有程序代码工具,当该程序产品在所述图像处理系统中被执行时,其执 行根据权利要求1到14中之一的方法的所有步骤。
16.一种图像处理系统(1),用于根据检查对象的截面图像数据自动产生 该检查对象的结果图像,具有:-用于接收测量的截面图像数据的接口(13),
-目标结构确定单元(17),用于根据诊断问题确定截面图像数据(D)中 的目标结构(Z),
-存储设备(12),其具有多个用于截面图像数据(D)中不同目标结构(Z) 的解剖学标准模型(M),其几何形状可分别根据特定的模型参数(tx,ty,tz, rx,ry,rz,d0)改变,
-选择单元(14),用于相应于所确定的目标结构(Z)选择解剖学标准模 型(M),
-适配单元(15),用于把所选择的标准模型(M)与截面图像数据(D) 中的目标结构(Z)进行匹配,
-分段单元(16),用于将截面图像数据根据匹配的标准模型分段,并且根 据诊断问题把检查对象的相关解剖学结构分开,其中,选择在截面图像数据中 的所有那些图像点,它们位于相应于相关解剖学结构的匹配的标准模型(M) 和/或至少一个模型部分的轮廓内,或者最大与其偏差特定的差值,
-显示单元(18),用于把相关的解剖学结构自动分开,以显示和/或为过 后的显示进行存储。
17.一种用于测量检查对象的截面图像数据的模态机,包括根据权利要求 16所述的图像处理系统(1)。

说明书全文

技术领域

发明涉及根据所涉及的检查对象的截面图像数据自动产生检查对象的 结果图像的方法。此外,本发明涉及用以执行这样的方法的图像处理系统。

背景技术

通常借助如计算机X射线断层摄影、磁共振断层摄影和超声设备等模态机 的截面图像产生的检查结果产生所涉及的检查对象的多截面图像系列。为了计 划其它检查和/或设置诊断,在许多情况下这些截面图像必须在检查期间已经或 者在检查后立即被继续处理。在此这样的检查的运行通常通过诊断问题来确定。 在此在大多数情况下在产生概略图像后立即较精确地检查特定的器官或者器官 系统。为此的一个例子是患者临床有关膝部病变的检查。在产生膝部较少的截 面图像系列后首先确定关于膝部内部结构可能存在的病理学的中间诊断,并且 然后在这一基础上进行膝部所涉及的区域的继续检查。通常,为建立这一中间 诊断例如放射设备的操作者或者MTRA(医学技术放射设备助理)必须分析单 个的概略图像,然后决定采取的措施。这样的中间诊断的确定需要不可低估的 时间花费,因此使整个检查流程变坏。另一个问题在于,在截面图像数据中识 别特定的内部结构的病理学,特别是在非常复杂的解剖学结构中,可能特别困 难,而需要一些经验。因此可能容易导致错误的中间诊断。这可能导致截面图 像检查质量的变坏。
也就是说,已经知道不同的方法为在截面图像中的特定感兴趣的结构产生 单个模型,以及这些模型用于诊断支持和用于干涉计划。例如在WO 99/55233 中说明了一种用于心脏超声图像的基于模型的分析方法,其中半自动地(通过 把一个模型与三个手工检测的解剖学标记匹配)产生和分析被检查人员的单独 的心脏模型。此外在DE 10311319A1中说明了一种方法,其中根据CT图像同 样在使用三个手工确定的解剖学标记的情况下产生用于计划心脏干涉方法的心 脏的单个3D模型。此外,在US 2003/0097219中说明一种方法,其中根据解剖 学标记半自动地产生左心室的模型。最后在WO 00/32106中说明了一种方法, 借助呼吸道或者消化道的个性化模型执行虚拟的内窥镜检查。然而在所有这些 方法中总是只输出一个模型,并且在此基础上的诊断或者干涉计划与此相应非 常依赖于所产生的模型的质量。

发明内容

因此本发明的任务在于,提供一种根据已经建立的截面图像数据自动产生 检查对象的结果图像的方法和图像处理系统,使用所述方法和系统能够相当简 单地、迅速地、和可靠地建立诊断、特别是用于继续进行检查的中间诊断。
根据本发明的方法,首先根据诊断问题自动确定截面图像数据中的感兴趣 的目标结构。然后相应于这一目标结构选择解剖学标准模型,其几何形状根据 模型参数可变。在这件事情上可以把不同的解剖学模型在一个数据库中进行管 理,其中对每一个要检查的器官至少有一个包括该器官的相应的解剖学标准模 型与之对应。然后把该标准模型自动与截面图像数据中的目标结构匹配,亦即 相应于该目标结构个性化。接着根据匹配的标准模型把截面图像数据分段,其 中根据诊断问题把检查对象的感兴趣的有关解剖学结构分开,其中,选择在截 面图像数据中的所有那些图像点,它们位于相应于相关解剖学结构的匹配的模 型和/或至少一个模型部分的轮廓内,或者与其最大偏差特定的差值。在此这一 选择可以采取这种形式,即去除所涉及的图像点或者去除所涉及的模型或者模 型部分的所有其余的图像点,亦即把所涉及的图像点剪去。这里“模型部分”理 解为标准模型的一个部分,例如一个头盖骨模型的颅底骨。在此可以准确地把 该模型部分与真正要检查的器官(部分)相对应。接着把相关的解剖学结构分 开显示和/或为以后的显示进行存储。
在此这一显示可以以2维或者3维方式例如在所涉及的模态机的操作台的 屏幕上或者在通过网络连接的工作站的屏幕上进行。同样该结果图像可以向打 印机、电影制作站等输出。相关解剖学结构的分开显示可以以这种形式进行, 即例如按照爆破图方式在结果图像中把所涉及的器官的所有单个部分彼此分开 显示。此外单个结构也可以在单个结果图像上显示,它们可以由建立该诊断的 人交替地、连续地或者并行地在不同的打印件、屏幕窗口等上观看。在3维显 示中这一点优选这样进行,即操作者交互地在一个相应的用户接口上在空间虚 拟旋转该结构或者单个的结构,使得可以从所有侧面观察。此外除像上面已经 提到的简单地表示结构表面的所谓的SSD(Surface Shaded Display(表面影线 显示))表示方式外,对于单个的相关结构每次使用适当的表示方式,像VRT (Volume Rendering Technique(体积着色技术))、MPR(Multiplanar Reconstruction(多平面再现))、MIP(Maximum Intensity Projection(最大强度 投影))等。
通过建议的方法可以根据标准模型把截面图像数据分段,亦即拆成所有与 诊断相关的部分。通过随后分开显示结果图像中的不同的解剖学结构,特别是 对于较少经验的人员也能非常容易地得到正确的中间诊断。因此本方法在截面 图像检查期间导致快速建立并保证中间诊断,由此可以减少整个检查时间同时 改善检查结果的质量。此外该方法可以帮助原来由医生进行的诊断在检查之后 优化。与开始时提到的、迄今已知的方法不同,在此显示的是感兴趣的结构的 实际测量和分段的体积数据,而不是该结构的模型。
在此,与如在US 6,556,696 B1中所说明的经典的阈值方法或者区域增长方 法相反,根据个性化了的模型的分段具有该方法还可以应用在下述场合的优点, 在这种场合中要分开的结构不能通过截面图像数据中的强对比度跳跃而识别。
根据本发明的图像处理系统为此首先需要一个用于接收测量的截面图像 数据的接口;一个目标结构确定单元,用于根据诊断问题确定截面图像数据中 的目标结构;一个存储设备,具有多个用于截面图像数据中的不同目标结构的 解剖学标准模型,优选以数据库的形式,其几何形状可分别根据特定的模型参 数改变;一个选择单元,用于选择相应于确定的目标结构的一个解剖学标准模 型。此外该图像处理系统需要一个适配单元,用于把选择的标准模型与截面图 像数据中的目标结构进行匹配;一个分段单元,用于将截面图像数据根据匹配 的标准模型分段,其中根据诊断问题把检查对象的相关解剖学结构分开,在此 选择在截面图像数据中的所有那些图像点,它们位于相应于相关解剖学结构的 匹配的标准模型和/或至少一个模型部分的轮廓内,或者与其最大偏差特定的差 值。接着需要一个显示设备,用于把相关的解剖学结构自动分开显示和/或为后 来的显示存储。这里“显示设备”理解为一个这样的设备,它这样准备分段的截 面图像数据,使得能够把相关的结构例如在屏幕上或者在其它连接在图像处理 系统上的输出单元上分开一个接一个地显示和单个地观察。
在一个优选的变形中在将标准模型与目标结构进行匹配期间分别根据特 定的偏差函数确定被修改的标准模型和目标结构的几何形状之间的当前偏差 值。由此可以全自动执行这一匹配,以此简单地使偏差值最小。
在此自动匹配可完全在后台进行,使得操作者可以转向其它的工作,特别 可以在图像处理系统中进行希望的结果图像显示的操作台上并行处理其它的图 像数据和/或控制其它的测量。但是也可以在自动执行方法期间在屏幕上或者在 该屏幕的一个部分区域显示过程参数,使得使用者能够控制匹配过程的进展。
优选向操作者显示偏差函数的当前值。特别在用户界面的其余部分为操作 者的其它工作空闲时也可以把偏差值在屏幕上例如在一个任务栏或者类似栏目 中长久显示。
优选为操作者提供下面的可能性,即根据需要干涉自动匹配过程和手工建 立各个模型参数。在此以有利的方式为操作者显示当前的偏差值,使得他在改 变所涉及的模型参数时立即看到是否能够和以何种措施通过他的行动使几何偏 差最小。在此特别还可以为每一模型参数规定单个的偏差值以及代替总偏差值 或者在其之外将该单个的偏差值显示。为此的一个典型的例子是在终端的图形 用户界面上显示目标结构和/或要匹配的标准模型或者该对象的至少部分,其中 用户例如借助键盘或者借助于指示设备如鼠标或者类似设备能够匹配特定的模 型参数、例如在模型上的点之间的距离。然后借助连续条或者以相似的良好的 光学识别方式向用户显示,通过他的行动该偏差减小到何种程度,在此特别是 一方面显示该模型的总偏差另一方面显示有关具体的当前模型参数的匹配的偏 差,例如对于模型中两点之间距离的情况,将两点的差显示为目标结构中的所 涉及点之间的距离。
在一个特别优选的实施例中在分段前自动检验,在将标准模型与目标结构 匹配时是否达到低于预先规定的阈值的一个最小偏差值。亦即检验该模型对数 据组中的目标结构的偏差是否足够小。仅当在这种场合,才根据该模型进行所 测量的数据组的自动分段。否则为另外手工处理截面图像数据中断该方法。以 这种方式可靠地防止在模型与测量的数据组偏差很大时执行错误的自动分段, 根据该自动分段可能导致错误的诊断和显示的解剖学结构。
此外除简单地分开显示相关的解剖学结构外还特别优选能够根据标准偏 差检验解剖学结构。亦即自动确定所涉及的解剖学结构与一个个性化模型或者 模型部分的偏差。
为此优选使用仅以特定方式个性化了的标准模型或者模型部分。在个性化 为这样识别标准偏差而应该使用的比较标准模型时,必须保证只执行这样的变 换,即该比较标准模型或者所涉及的标准模型部分的几何形状自身不具有任何 病变。所确定的偏差然后可以与解剖学结构一起用图形显示。例如它可以以显 示的数据组在屏幕上为用户标记。此外还可以通过声音信号不致错误理解地向 用户显示这样的偏差。这样可以以简单的方式自动确定所检查的解剖学结构的 病变并向用户对此进行提示。
在本方法的一个改进中还可以根据确定的标准偏差将检查对象自动分类。 例如可以自动确定是否需要继续检查,是的话,执行什么样的检查。在此还给 操作者把该分类仅作为建议说明,以便使用户然后赞成该建议,并无需很大的 开销执行另外的检查,或者操作者可以简单地拒绝该建议,以便以常规方式主 动决定是否执行和执行哪些详细检查。
解剖学标准模型的个性化,亦即对目标结构的匹配原理上可以使用任意适 宜的个性化方法执行。解剖学模型个性化的思想一般可以简化地这样表达,即 在相应于3维模型的3维变换中寻找一个几何学变换,它使模型在具体的计算 机X射线断层摄影数据组、核磁共振断层摄影数据组、或者超声数据组中优化 地匹配。能够分配给该模型的几何形状的全部信息在此同样被个性化。在医学 图像处理中用于特定优化变换参数的这样的方法也称为登记或者匹配方法。其 中,根据使用哪种几何变换,区分为所谓的刚性的、仿射的(affinen)、透视的 和塑性的方法。这些登记方法例如迄今为下面的目的使用,即为使两个或者多 个图像结合在一个总的图像中,或者使解剖学图集与图像数据匹配。各种这样 的方法在WO 01/45047 A1、DE 69332042 T2、WO 01/43070 A1和DE 19953308 A1中说明。
为数学处理个性化问题,如已经说明的通常使用描写任意变换后的模型与 目标结构的偏差的偏差函数。在此偏差函数的类型取决于所使用的解剖学标准 模型的各类型。
可应用的数字解剖学标准模型原则上可以以不同方式构造。一种可能例如 是解剖学结构在体元基础上的模型化,其中为编辑这样的体积数据需要专软件,这些软件一般很贵并且较少推广。另一种可能是使用所谓的“有限元”进 行模型化,其中通常从四面体构造模型。但是对于这种模型也需要专门且昂贵 的软件。较多推广的是通过三测量的解剖学边界平面的简单模型化。相应的 数据结构通过多个标准程序由计算机图形学支持。根据这一原理构造的模型被 称为所谓面向表面的解剖学模型。这里涉及解剖学结构模型化的一个最小的公 分母,因为不仅从刚才提到的体积模型通过体元的测量而且也通过有限元方法 的四面体的转换可导出相应于三角形的表面模型。
因此提供了使用在三角形基础上构造的、面向表面的模型作为标准模型的 情况。一方面,使用这种方法能最简单和费用最适宜地产生模型,另一方面, 可以通过相应的变换接受以另外的形式产生的模型,特别是所述体积模型,使 得不用新建立相应的模型。
为新建立这样的表面模型,例如可以用相应的开销使用经典手工方法分段 截面图像数据。从如此获得的关于单个结构例如单个器官的信息,最后可以生 成模型。为得到人的骨骼模型,例如也可以借助激光扫描器测量人的骨骼,或 者使用计算机X射线断层摄影进行扫描和分段以及三角测量。
在这样的模型中例如可以根据最小平方法定义偏差函数,其中使用该函数 由相对于目标结构变换后的模型三角形的位置中计算出用于偏差的尺度。
在本发明的一个特别优选的实施例中使用一种弹性登记方法。为尽可能快 地找到偏差函数的最小值,在此优选使用一种多级方法。例如可以在一个三级 方法中首先借助于一个适当的定位,即平移、旋转和缩放,粗略匹配模型。接 着可以在第二步骤中执行体积变换,以便得到较好的校准。之后在第三步骤中 执行精调,以便使模型与结构局部最优匹配。
特别优选地,为了个性化使用层次参数化的标准模型,其中根据模型参数 对模型的总的解剖学几何形状的影响而对其进行分层次排序。然后以多个迭代 步骤进行标准模型的个性化,其中,随着迭代步骤数量的增加的同时,在各迭 代步骤中可以调整的模型参数的数目(由此在模型变换中自由度的数目)对应 于该参数分层的级提高。通过这一方法保证,在个性化时首先调整对模型总的 解剖学几何形状有最大影响的模型参数。然后,逐步调整仅对总几何形状的一 部分有影响的下级模型参数。由此在模型匹配时保证有效的并因此是节省时间 的优选方式,而不依赖该匹配是全自动进行还是操作者手工干预该匹配方法。 在一个(部分)手工的方法中,这点例如可以通过下述实现,即在每个迭代步 骤中仅允许操作者例如借助于图形化的用户界面、根据参数的分层级改变各模 型参数。
优选给每一层次级分配模型参数。这意味着,不同的模型参数在必要时也 可以分配在同一层次级,因为它们对该模型的总的解剖学几何形状具有大致相 同的影响。然后可以在一个规定的迭代步骤中把所有的模型参数为进行调整另 外新纳入一个规定的层次级。然后在下一个迭代步骤中把这些模型参数另外纳 入位于其下的层次级。
将模型参数分配到层次级可以根据在模型几何形状中当所涉及的模型参 数改变特定的值时所出现的偏差进行。在此在一个特别优选的方法中,使不同 的层次级对应特定的偏差范围,例如数值偏差范围。也就是说例如为将参数分 配到一个层次级,改变该参数并计算输出状态下的该几何形状改变了的模型所 产生的偏差。在此该偏差的度量取决于所使用的标准模型的类型。起决定作用 的仅仅是,确定一个精确定义的偏差度量,该偏差度量尽可能精确地对模型在 有关模型参数变化之前和之后的几何变化进行量化,以便保证将不同模型参数 对模型几何形状的影响进行实际的比较。因此,为了能够对几何形状影响进行 直接的比较,优选地对于每个产生类型(例如,对于其中模型的两点之间的距 离变换的距离参数,或者对于其中模型的三点之间的角度变化的角度参数)使 用统一的步长。然后,简单地通过对该偏差度量给出数字间隔来将这些参数划 分到分层级。在使用在三角形基础上产生的表面模型的情况下,对于在未改变 的标准模型和一个参数变动后改变了的标准模型之间的偏差,优选地在对在不 同状态下的模型的对应三角形的几何距离的和的基础上进行计算。
优选地,在最上面的分层级中(其模型参数在第一迭代步骤中就可以立刻 设置)至少配置这样的模型参数,即这些参数改变时规范模型全程发生变动。 此外,例如整个模型转动的所谓新参数是围绕三个模型轴的,平移是沿三个模 型轴的,而整个模型的缩放是沿三个模型轴的。
各个模型参数的层次安排基本可以在截面图像数据的分段期间进行。那时 例如在每个迭代步骤中首先检验哪些另外的参数对几何形状有最大的影响,然 后增加这些参数。然而因为由此与极大的计算开销相联系,因此特别优选事先 执行模型参数到层次级别中的分类或者安排,例如在产生该标准模型时已经进 行,不过至少在通常为以后选择把标准模型向模型数据库中存储之前进行。
也就是说优选事先在一个为产生标准模型的独立方法中把有关其对模型 的总的解剖学几何形状的影响的模型参数进行层次安排,所述标准模型为用于 上述产生结果图像的方法可供使用。在此同样可以相应于层次级分配模型参数, 其中一个参数对一个层次级的分配再次根据在模型几何形状中当所涉及的模型 参数改变特定的值时所出现的偏差进行。在用于产生标准模型的单独的方法中 模型参数的层次安排的这一方式具有下面的优点,即对每一标准模型只需执行 一次模型参数层次级别的计算,因此在分段期间可以有价值地节省计算时间。 层次级别可以以相对简单的方式与标准模型一起存储,例如通过把参数安排在 层次级别中、或者使用相应的标记等连接而在一个文件头或者在还包含所涉及 的标准模型的其它数据的文件的另一标准化位置存储。
在一个特别优选的实施例中模型参数分别这样与该模型的至少一个解剖 学标记的一个位置相关联,使得该模型对每一参数数据组具有一个有解剖学意 义的几何形状。为此的一个典型的例子一方面是全局参数如总模型的旋转或者 平移,此时全部模型参数相应适当地彼此在位置上改变。另外的模型参数例如 是在两个解剖学标记之间的距离或者在三个解剖学标记之间的角度,例如用于 确定膝盖位置。
模型参数对选择的医学上有意义的标记的这种耦合具有下面的优点,即在 个性化后总可能有一个诊断结论。而且在解剖学专业文献中这种解剖学标记的 位置被精确说明。通过这样的先行处理,分段的执行因此变得容易,因为一个 受过医学培训的用户,例如医生或者MTA熟悉这些解剖学标记,并基本上可 以用其确定解剖结构。
为自动确定在分层图像数据中的被分开的部分对象的目标几何形状,存在 不同的可能。一种选择在于使用所谓的“阈值方法”。这种方法以这种方式作用, 即把单个体元即单个3D图像点的强度值(在计算机X射线断层摄影中称为 “Houns域值”)与一个固定调准的阈值比较。如果该体元的值超过该阈值,则 把该体元认为是一个规定的结构。然而这种方法在磁共振方法中主要用在对患 者的造影剂检查或者识别患者的皮肤表面。在采用计算机X射线断层摄影时这 种方法另外还为识别特定的骨骼结构使用。该方法一般不适合于识别其它组织 结构。因此在一个优选的方法中至少部分地通过轮廓分析方法确定目标几何形 状。这样的轮廓分析方法根据在相邻图像点之间的梯度工作。专业人员熟悉不 同的轮廓分析方法。这种轮廓分析方法的优点在于不仅在计算机X射线断层摄 影截面图像数据而且在磁共振截面图像数据以及在超声截面图像数据中可稳定 使用。
图像处理系统的目标结构确定单元、选择单元、适配单元和分段单元以及 显示单元可以特别优选以软件的形式在图像计算机的一个相应适宜的处理器上 实现。这一图像计算机应该具有一个用于接收图像数据的相应的接口和一个用 于解剖学标准模型的适宜的存储设备。在此该存储设备不一定必须是该图像计 算机的一个集成部分,而只要该图像计算机能够访问一个适宜的外部存储设备 就已足够。为完整起见这里要提醒,不一定强制要求不同的部件在一个处理器 或者在一个图像计算机中存在,而是这些不同的部件也可以分布在多个处理器 或者彼此联网的计算机上。
以软件形式实现本发明的方法具有下述优点,即通过适宜的更新能够相对 简单地对已有图像处理系统进行相应的后装备。在本发明的图像处理系统中, 也可以特别涉及一个用于记载截面图像数据自身的模态机的控制单元,其具有 用于按照本发明处理截面图像数据所需要的部件。
附图说明
下面根据实施例参考附图详细说明本发明。其中,
图1是根据本发明的图像处理系统的一个实施例的原理表示,该图像处理 系统通过数据总线与一个模态机和一个图像数据存储器连接,
图2是根据本发明方法的一种可能的运行的流程图
图3是用于模型个性化的优选方法的详细流程图,
图4示出在计算机X射线断层摄影的截面图像中的人头盖骨的可能的目标 结构,
图5示出人头盖骨的表面模型,
图6A示出按照图4的目标结构,其具有按照图5尚未匹配的表面标准模 型(无下颌骨),
图6B示出按照图6A的目标结构和标准模型,但具有在目标结构上部分匹 配的标准模型,
图6C示出按照图6B的目标结构和标准模型,但具有与目标结构进一步匹 配的标准模型,
图7A示出按照图5的头盖骨标准模型,其以爆破图(Sprengzeichnung) 的形式分成多个部分模型分开显示,
图7B示出从另一观察方向按照图7A的头盖骨标准模型的一部分,
图8示出在按照图5的头盖骨标准模型上的解剖学标记。
图9示出在三角形基础上建立的人骨盆的表面模型。

具体实施方式

在图1中示出的根据本发明的图像处理系统1的实施例主要由图像计算机 10和与其连接的操作台5组成,后者例如具有图像显示器6、键盘7和指示装 置8,这里是一个鼠标。通过操作台5或者另外的用户接口,例如可以由操作 者输入诊断问题,或者用预先给定的诊断问题从数据库中选择。
图像计算机10可以是以通常方式构造的计算机,例如工作站,它还可以 用于其它图像分析和/或控制图像拍摄设备(模态机),如计算机X射线断层摄 影仪、磁共振断层摄影仪、超声设备等。图像计算机10中的主要部件通常有处 理器11和接口13,以便接收患者P的截面图像数据D,这些数据由模态机2、 这里是磁共振断层摄影仪2测得。
在图1示出的实施例中模态机2与控制设备3连接,后者又与总线4连接, 图像处理系统1也连接于该总线4上。此外,在总线4上连接有海量存储器9, 用于中间存储或者长期存储由模态机2所记录的图像和/或由图像处理系统1进 一步处理的图像数据D。不言而喻,在总线4上在构造一个较大的网络的条件 下还可以连接其它在通常的放射学信息系统(RIS)中存在的部件,例如另外 的模态机、海量存储器、工作站、输出设备(像打印机)、电影制作站 (Filming-Station)等。同样也可以与外部网络或者与其它的RIS连接。在此, 优选将所有数据按所谓的DICOM标准(Digital Imaging and Communication in Medicine,医学中的数字成像和通信)格式化,以便在各个部件之间进行通信。
对模态机2的控制以通常方式在控制设备3的控制下实现,控制设备3还 从模态机2获取数据。控制设备3为现场操作可以具有自己的操作台,不过它 在这里没有示出。但是还可以使操作例如通过总线借助一个位于模态机附近的 单独的工作站实现。
图2表示根据本发明方法产生检查对象结果图像的典型流程。
首先,在第一方法步骤I中根据预先给定的诊断问题在截面图像数据D中 确定目标结构。这一点优选全自动实现,例如借助已经提到的轮廓分析。对于 特定的结构和特定的拍摄方法,也可以使用阈值方法,像已经在再前面提到的。 截面图像数据D例如可以直接从模态机2或者其控制设备3经由总线4传送到 图像计算机10。但是也可以是在一定时间前已经采集并存储在海量存储器9中 存储的截面图像数据D。
然后在步骤II选择一个相应于目标结构Z的标准模型M。该步骤也可以并 行于或者先于确定目标结构的方法步骤I进行,因为要确定的目标结构Z通过 现有技术的诊断问题已知。在此图像计算机10具有一个存储器12,其具有用 于不同的可能的解剖学结构的各种标准模型。在此所涉及的通常是由多个模型 部分组成的模型。
为此,一个典型的例子可以根据膝盖检查说明,其中诊断问题针对查明膝 盖中的特定结构。于是首先在拍摄的截面图像数据中确定膝盖的目标结构,例 如该膝盖的外骨表面。为此适用的膝盖模型例如由模型部分“大腿骨”、“胫骨”、 “膝盖骨”和单个的半月板组成。然而在涉及患者头部例如对怀疑头盖骨破裂进 行检查的诊断问题中,可以从截面图像数据中将头盖骨的骨表面结构确定为目 标结构。图4示出了从患者的计算机X射线断层摄影数据得到的一个这样的目 标结构。图5示出一合适的头盖骨标准模型,其包括作为(在该图中可识别) 模型部分的前额骨T1、右头顶骨T2、左头顶骨T3、颜面骨T4和下颌骨T7。该 模型因为利用连贯的表面显示所以其辨识性更好。实际上该模型构建在三角形 基础上。图9示出骨盆的相应表面模型。
通过选择单元14选择合适的模型M,而确定结构通过结构确定单元17进 行,它们在这里以软件的形式在图像计算机10的处理器11中实现。其在图1 中原理示出。
接着在方法步骤III中通过所谓的“弹性记录方法”进行该模型的个性化。但 是原则上也可以使用其它个性化方法。标准模型M与目标结构的这种匹配在适 配单元15中进行,如在图1中原理示出的那样,它同样以软件模的形式在图 像计算机10的处理器11中实现。
个性化过程的优选实施例在图3中以流程图的形式较精确地原理性示出。 在该匹配过程中,在多个迭代步骤S中改变各个模型参数,直到最后所有的参 数都被个性化,或者个性化被足够地进行,亦即在标准模型M和目标结构Z 之间的偏差最小,或者低于预定的阈值。在此每一迭代步骤S包括多个处理步 骤IIIa、IIIb、IIIc、IIId,它们以循环方式运行。
该循环或者说第一迭代步骤以方法步骤IIIa开始,在其中,首先为平移、 旋转和缩放确定优化参数。这些是最上层(下面称“0”级)层次级的参数,因为 这些参数对总的几何形状起作用。平移的三个参数tx、ty、tz和围绕三个模型轴 旋转的三个参数rx、ry、rz在图5中原理性示出。
如果这一匹配能够尽可能地执行,则在后一步骤IIIb中通过已经确定的参 数来估计尚未确定的模型参数。这意味着从确立的上层参数估计下层参数的开 始值。对此的一个例子是从用于身高的缩放参数的设置来估计膝盖宽度。该值 为后面所涉及的参数的设置作为输出值输出。以这种方式可以极大地加速这一 方法。
然后在方法步骤IIIc中优化地设置所涉及的参数。
在所示实施例中根据参数对模型的总的解剖学几何形状的影响而对其进 行排序。一个参数对几何的影响越大,则其在层次结构中就处于越上层。随着 迭代步骤S的数目的增加,相应于该层次级别的可调整的模型参数的数目也增 大。
亦即在第一迭代步骤S或者在该循环的第一次运行中在步骤IIIc为设置该 模型仅使用第0层次级以下的第1层次级的参数。然后在第二次运行时能够首 先使该模型在方法步骤IIIa中进行新的平移、旋转和缩放。接着在方法步骤IIIb 中通过已经确定的参数估计第2层次级尚未确定的参数,这些参数然后在步骤 IIIc中被引入用于调整。这一方法然后被n次重复,其中在第n次迭代步骤中 第n级的所有参数都被优化,并且在迭代步骤S的最后步骤IIId再次声明,是 否还有迄今尚未被优化的可供使用的其它参数。接着又开始一次新的第(n+1) 次迭代步骤,其中首先重新把该模型进行移位、旋转和缩放,接着可以依次再 次调整所有参数,在此现在第(n+1)级参数也可供使用。接着在方法步骤IIId重 新检验是否所有参数都被个性化,亦即是否还存在尚未优化的参数,或者是否 已经达到所希望的匹配。
图6A到图6C示出这样的匹配过程的非常简单的情况。在这些图中,为了 更加明了再次以连续的表面重新示出模型M。图6A示出具有相对移位的模型 M的目标结构。通过一次简单的平移、旋转和缩放然后可以得到在图6B中表 示的图像,在该图像中模型M已经较好地与目标图像Z匹配。通过调整另外的、 下层的参数,最终得到图6C中示出的匹配。
通过上述迭代方法可以保证实现尽可能节省时间的和有效的匹配。其中在 匹配期间可以在操作台5的图像屏幕上随时显示目标结构Z和所属模型M以及 当前计算的偏差值或者当前计算的一个偏差函数的值。此外,该偏差也可以像 在图6A到6C示出的那样被显示。另外,偏差的显示还可以通过给出相应的颜 色实现。
下层的层次级产生于对几何形状影响的定量分析。为此要改变每一个参数 并对输出状态计算几何形状改变了的模型的结果偏差。当使用如图9所示基于 三角形的表面模型时,这些偏差例如可以通过对相应模型三角形的几何距离的 和来量化。通过对偏差的数值区间的预设,然后可以把参数分为层次级。在此 完全有可能不同的参数落入同一层次级。除别的事项外这点依赖于偏差的数值 区间的宽度。在同一层次级中的这些参数如上所述在一个特定的迭代步骤S内 为了改变而同时被首次提供,或者在自动匹配方法中相应地自动改变。
如已经说明的,优选在该方法中使用与该模型的特定的解剖学标记的一个 或者多个位置直接相关的模型参数。这一方面具有只执行该模型的有医学意义 的变换的优点。另一方面,它具有受过医学培训的用户通常了解这些解剖学标 记、并因此可以很好处理这些参数的优点。对于这些参数例如有图8中在头盖 骨模型上标记的解剖学标记L、L1、L2的位置或者各个标记之间的距离,例如 在眼眶(眼窝)中点的解剖学标记L1、L2之间的距离d0。为在自动匹配过程中 由操作者手工操作调整眼窝的这一距离d0,操作者例如可以通过鼠标选择一个 解剖学标记L1、L2并交互地改变其位置。于是该模型的几何形状一起自动地适 当改变。
在包含标准模型M的两个解剖学标记之间的距离的模型参数的变形中,优 选将该标准模型的几何形状在沿这些解剖学标记之间的直线的区域内正比于距 离的变化而改变。在包含将第一解剖学标记的位置相对于相邻标记变化的模型 参数的变形中,优选将标准模型M的几何形状在围绕所涉及的第一解剖学标记 的环境中在所涉及的相邻标记的方向上一起适当地改变。在此该变形有利地随 与所涉及的第一解剖学标记的距离的增加而减少。亦即,该变形在围绕该标记 的近区域内要比在远离它的区域中强,以便达到在图中所示的效果。但是也可 以设想另外的变换规定,只要这些规定能执行有解剖学意义的变换。通常这取 决于各选择的模型。
根据头盖骨模型上的解剖学标记L、L1、L2,还可以说明一个典型的例子, 其中设置两个位于不同层次级的标记之间的距离。于是在图8中所示的头盖骨 模型不仅通过两个眼窝的距离d0确定,也可通过两个茎突(processi styloidei) 之间的距离参数化,这里茎突涉及的是颅底内骨上的小骨突起(在图8的视图 中不可见)。这里给出眼窝距离的第一参数的几何学效果大于给出茎突之间距离 的第二参数的几何学效果。这允许在参数改变时通过模型的几何形状变化可检 查到1毫米。因为茎突涉及相对小的结构,因此模型的几何学改变限制在这些 骨突起的小区域内。与此相对眼窝则大得多。在眼眶距离改变时该模型的多个 部分将改变它们的几何形状从而导致偏差增大。因此眼眶距离的参数比茎突距 离的改变位于高许多的层次级,因为原则上具有对参数层次较大的几何作用距 离的参数比更多是具有局部影响的参数处于更高的位置。
如果最后所有可调整的参数都被个性化或者偏差函数达到其最小值,则在 方法步骤IV检验,该个性化的标准模型的偏差对数据组亦即对于目标结构是否 充分小。这里例如可以检验当前达到的偏差值是否低于一个边界值。如果不是, 则该自动处理被中断,而以常规方式进行其它的处理,如在此作为方法步骤V 示意性示出的那样。即由操作者手工分析图像数据,并进行手工的中间诊断。 在这种中断的情况下以有意义的方式向操作者输出一个相应的信号,以便操作 者立即知道他必须手工继续处理运行的过程。
反之如果实现了标准模型M与目标结构Z的充分匹配,则在方法步骤VI 执行分段。这一点在分离单元16中进行,如在图1中所示,它同样作为在处理 器11中的软件模块实现。在此选择截面图像数据中所有那些图像点,其相应于 按照诊断问题相关的解剖学结构位于该模型或者模型的特定部分的轮廓内。在 此例如删除所有多余的数据,仅留下所希望的图像点。
在方法步骤VII,全自动准备全部分段的数据,使得可以采取希望的结果 图像形式分开显示诊断相关的解剖学结构。这一点借助于图形用户界面实现。 为此可以使用一个用于表示3维对象的商业程序,这样例如通过相应于一个这 样的程序的接口的显示单元来准备分开的、相关(部分)结构的数据。
在图7A和7B中示出例如在头盖骨检查中相关结构可以什么形式进行显 示。所显示的分别是按照图5的头盖骨标准模型。图7A示出按照爆破图显示 的该模型M,其中主要的模型部分T1、T2、T3、T4、T5、T6、T7在一个结果图 像上被彼此分开地示出。这里各部分具体涉及额骨T1、右顶骨T2、左顶骨T3、 面骨T4、后枕骨T5、包括后枕骨T5的一部分的颅底骨T6、和下颌骨T7。在图 7A中面骨T4和颅底骨T6(包括后枕骨T5)尚在一起作为一个部分。所有部分 结构或者模型部分T1、T2、T3、T4、T5、T6、T7可以由用户在图形用户界面上 分开标记,例如通过鼠标“点击”和通过虚拟转动和缩放在空间中从所有方面分 开检查。图7B中示出由面骨T4和颅底骨T6(包括后枕骨T5)组成的、连接在 一起的头盖骨部分从上面看去的情形。由图7A和7B与图5的比较很快看出, 由于相关结构的分开可视显示(亦即内结构)可以简单地确定复杂的结构内的 病理。于是在所示例子中对头盖骨的检查自身可以由没有经验的医学人员或者 甚至外行在按照图7B的显示上毫不困难地确定颅底骨的裂缝。而在对截面图 像数据的经典分析时与此相反只有有经验的医学人员才可能。
在图2所示实施例中像在大多数情况下一样立即实现显示。只要处理过程 在后台运行,则例如就可执行声学的和/或光学的指示,表明该过程仍在进行, 可进行可视化。另外还可选择的方案或者此外附带地,首先中间存储如此产生 的、将诊断上相关的解剖学结构彼此分开显示的结果图像,或作为这些图像基 础的所准备的数据,以便可以在后来任何时间调用它们。结果图像优选还可以 在打印机、电影制作站等类似设备输出,或者通过网络向其它地方发送,以在 那里在屏幕上或类似设备上显示。
此外在图2所示实施例中,在结果图像上还标记各所属标准模型或者模型 部分的不同的分开的结构的标准偏差,以便简化由操作者进行的诊断。这优选 结合声学信号进行,用该声学信号通知操作者在特定位置存在相应的标准偏差。
然后在方法步骤IX执行另外的检查步骤。这一点可以自动根据确定的标 准偏差或者也通过操作者手工进行。在一个特别优选的变形中,根据标准偏差 自动向操作者建议另外的检查步骤,操作者可以接受或者拒绝或者可以补充或 者改变这些步骤。
因此所建议的图像处理系统不仅用于像通常的图像处理系统为检查准备 要观察的图像,而且也作为基于模型的专家系统,该专家系统在运行的截面图 像检查中引导快速建立和保证中间诊断。因此本发明的方法或者图像处理系统 能够有助于使总检查时间极大减少和改善检查结果的质量。特别是也可以使实 际的医师的诊断在通过所述应用的检查之后得到优化,因为通过具有分开的相 关解剖学结构的可供使用的结果图像,必要时同时具有对标准偏差已进行的标 记,极大地简便了医生对可能的病理的识别。
在此尚需特别指出,在附图中示出的系统结构和过程仅是实施例,本领域 技术人员能够毫不困难地改变细节。特别是控制设备3,只要它例如装备一个 相应的操作台,就也可以具有图像计算机10的所有相应的部件,以便在那里直 接执行根据本发明方法的图像处理。在这种情况下,控制设备3本身成为本发 明的图像处理系统,不需要另外的工作站或者分开的图像计算机。
此外,还可以为在其内已实现已知的后处理程序的现有图像处理系统后装 备本发明的过程控制单元,以便按照上述本发明的方法使用该设备。在许多情 况下在必要时用适宜的控制软件模块对控制软件更新一次已经足够。
QQ群二维码
意见反馈