诊断成像应用中的基于随着时间的统计分析的分割的制作方法
诊断成像应用中的基于随着时间的统计分析的分割的制作方法
【技术领域】
[0001] 根据本公开的方案设及诊断成像领域。更具体地说，本方案设及诊断成像应用中 的分割。
【背景技术】
[0002] 分割技术广泛用于诊断成像应用（或简称为成像应用）中；一般来说，任何分割技 术旨在将所分析的患者的身体部分的表示（例如，其数字图像）分割为两个或更多个部分， 每个部分包括该身体部分的该表示的具有基本上同质的特性的一部分。
[0003] 例如，分割技术可W在动态对比度增强超声值CE-U巧成像应用中用于评估血液 灌注。在该种情况下，充当高效超声反射体的超声造影剂扣CA)或简称为造影剂被施用于 患者。造影剂在患者中W与红血球相同的速率流动，W使得其通过记录响应于超声波返回 的回声信号的检测提供关于相应的血液灌注的信息。具体地说，对于身体部分的其中回声 信号已经被记录的每个位置，（例如，从随着时间拟合回声信号的参数化函数）计算表征它 的灌注参数（诸如峰值增强）的值。通过将相应的灌注参数值分配给表示身体部分的位置 的每个像素，构建参数化图像。该参数化图像可W用于计算身体部分的包括病变的分析区 域中的灌注参数值的合并值（例如，它们的均值）；通常，分析区域的合并值用相对于身体 部分的包括健康问题的控制区域中的灌注参数值的合并值的相对项来表示。该相对的合并 值提供病变的量化表征，该量化表征可W用于例如其治疗随访中。
[0004] 然而，分析区域和控制区域的合并值本质上受一定程度的不精确度的影响，因为 它们是在身体部分的不完美地同质的区域中整体计算的；例如，分析区域可能包括低灌注 的坏死组织，控制区域可能包括高灌注的大血管。分析区域和控制区域的该些低灌注部分 和高灌注部分分别在它们的合并值中引起在分析区域的相对合并值中反映的误差。该缺点 在治疗随访中尤其严重，其中，病变在其处理期间的灌注中即使是微小的变化都可能指示 处理的有效性。
[0005] 可替代地，W0-A-2011/026866(其全部公开内容通过引用并入本文）公开了通过 下述方式确定关于分析区域（限定感兴趣区域）的每个像素的分析函数，即，用时间模型函 数（例如，对数正态分布函数）拟合其回声-功率信号，并且通过用相同的时间模型函数拟 合其回声-功率信号的平均值来确定关于参考区域（包括健康实质）的参考函数。对于分 析区域的每个像素，通过从其分析函数减去参考函数来计算差函数；然后根据像素的差函 数将该像素分配给不同的类（即，正单极类、负单极类、正到负双极类和负到正双极类）。为 了该个目的，可W计算差函数的正能量和负能量；当正能量和负能量的和高于显著性阔值 时，计算相对正能量和相对负能量。如果相对正能量或相对负能量超过鉴别阔值，则将差函 数分别分配给正单极类或负单极类；相反，根据差函数的主要极性变化的数量级将该差函 数分配到正到负双极类或负到正双极类。
[0006] 分割技术还可W用于超声分子成像扣MI)应用中。在该种情况下，使用目标特定 的造影剂；该目标特定的造影剂适于到达特定生物目标（例如，病变），然后通过特定的相 互作用在其上保持固定。目标特定的造影剂的在身体部分中固定的任何颗粒的检测使得可W识别否则将难W发现的相应病变；而且，该些固定颗粒的量化使得可W确定病变的状态， 该状态可W用于例如该病变的治疗随访中。
[0007] 然而，目标特定的造影剂的固定颗粒的检测和量化受到如下事实的阻碍，即，其总 量中只有一小部分实际上到达目标，并且在其上保持固定；大部分目标特定的造影剂而是 继续循环相当长的时间（诸如，几分钟）一一例如，直到它被肺脏滤出和/或在患者的肝脏 中被滤出为止。因此，身体部分中的目标特定的造影剂的颗粒包括固定的颗粒的贡献和仍 然循环的颗粒的贡献两者；因此，区别该两个贡献是相当困难的（尤其是在施用目标特定 的造影剂之后的早期）。
【发明内容】
[0008] 本公开的精简的总结在本文中是为了提供其基本理解而呈现的；然而，本总结的 唯一目的是作为其后面的更详细描述的序言W精简的方式介绍本公开的一些概念，并且本 总结不应被解释为标识其关键元素，也不应被解释为其范围的划定。
[0009] 一般地说，本公开基于使分割基于随着时间的统计分析的构思。
[0010] 具体地说，一方面提供一种用于诊断成像应用中的数据处理分割方法，其中，被灌 注造影剂的身体部分的表示的选定位置在每个选定位置的分析时间段期间根据与对询问 信号的响应相关的至少一个统计参数的值而分割。
[0011] 另一方面提供一种相应的计算机程序。
[0012] 另一方面提供一种相应的计算机程序产品。
[0013] 另一方面提供一种相应的系统。
[0014] 另一方面提供一种相应的诊断成像方法。
[0015] 更具体地说，本公开的一个或多个方面在独立权利要求中阐述，其有利的特征在 从属权利要求中阐述，所有权利要求的措辞在本文中通过引用逐字地并入（参照任何特定 方面提供的任何有利的特征适用于每一个其它方面）。
【附图说明】
[0016] 根据本公开的一个或多个实施例的解决方案及其进一步的特征和优点将参照W 下详细描述被最好地理解，W下详细描述仅仅作为非限制性指示而给出，应结合附图进行 阅读（其中，为了简单起见，相应的元件用相同的或相似的标号表示，并且它们的说明不被 重复，并且每个实体的名称一般用于表示其类型和其属性两者一-诸如值、内容和表示）。 具体地说：
[0017] 图1示出可W用于实施根据本公开的实施例的解决方案的系统的图形表示，
[001引图2示出根据本公开的实施例的解决方案的应用的体外情况的例子，
[0019] 图3A-图3B示出根据本公开的实施例的解决方案的应用的例子，
[0020] 图4A-图4D示出根据本公开的另一实施例的解决方案的应用的例子，
[0021] 图5A-图5C不出根据本公开的另一实施例的解决方案的应用的例子，
[0022] 图6示出表示可W用于实现根据本公开的实施例的解决方案的主要组件的作用 的协作图，和
[0023] 图7、图8和图9示出根据本公开的实施例的解决方案的体内应用的不同例子。
【具体实施方式】
[0024] 具体参照图1，示出了可W用于实施根据本公开的实施例的解决方案的系统100。 具体地说，系统100是超声扫描仪，其包括中央单元105W及连接到中央单元105的、例如 阵列类型的手持发射-接收成像探头（或简称为成像探头）11〇。成像探头110发射包括超 声脉冲序列的超声波（例如，具有1和50MHz之间的中屯、频率），并且接收由该些超声脉冲 在选定的扫描平面中的反射而导致的射频（RF)回声信号；为了该个目的，成像探头110设 有发射/接收复用器，其允许在上述脉冲-回声模式下使用成像探头110。
[00巧]中央单元105容纳主板115,在主板115上，安装控制超声扫描仪100的操作的电 子电路（例如，微处理器、工作存储器和硬盘驱动器）。而且，一个或多个子板（整个用标号 120表示）被插入到主板115中；子板120提供用于驱动成像探头110并且用于处理回声信 号的电子电路。中央单元105还配备有用于可移动盘130 (诸如CD或DVD)的驱动器125。 用于显示与进行中的分析处理相关的图像的监视器135连接到中央单元105。超声扫描仪 100的操作通过键盘140 (W常规方式连接到中央单元105)控制；优选地，键盘140设有用 于在监视器135上操纵指针（该图中未示出）的位置的轨迹球145。
[0026] 超声扫描仪100用于分析患者155的身体部分150;为了该个目的，（超声）造影 剂被施用于患者155。
[0027] 造影剂包括充当良好的超声反射体的颗粒。例如，造影剂是充气气泡在液态载体 中的悬浮液；通常，充气气泡具有大约0. 1-5ym的直径，W便使得它们可W保留在患者155 的血管系统内，但是同时使得它们可W通过他/她的毛细血管。充气气泡通常通过将气体 或者其前体夹带或封装到各种系统中来使其稳定，所述各种系统包括磯脂、乳化剂、油、增 稠剂、糖、蛋白质或聚合物；稳定的充气气泡一般被称为微泡。具体地说，分散在水介质中 并且在气体/液体界面处被包含表面活性剂（即，两亲性材料）的非常薄的包膜束缚的气 泡也被称为微气泡。可替代地，被由脂类或（天然或合成）聚合物形成的固态材料包膜包 围的微泡也被称为微气球或微胶囊。另一种类的造影剂包括聚合物或其它固体的多孔微粒 的悬浮液，该悬浮液携带包埋在微粒的孔内的气体的气泡。EP-A-0458745、W0-A-91/15244、 EP-A-0554213、W0-A-94/09829和W0-A-95/16467(其全部公开内容通过引用并入本文）中 描述了微泡（具体地说，微气泡和微气球）的合适的水性悬浮液及其制备的例子。商业的 包含微泡的造影剂的例子是化accoInternationalBV(商标）的Sono化e。
[0028] 造影剂还可W是目标特定的类型。目标特定的造影剂基本上自由地在患者155内 循环；然而，目标特定的造影剂还能够固定在选定的（生物）目标上，W便保留在基本固定 的位置上一定时间段。为了该个目的，目标特定的造影剂W通过与期望目标的特定相互作 用选择性地绑定到该期望目标的该样的方式配制。例如，该行为可W通过合并能够选择性 地绑定到期望组织或受体（例如，通过生物化学亲合性和/或静电作用）的目标特定的配 体来实现。目标特定的配体（其可W插入到微气泡的膜中）的例子是单克隆抗体、肤或多 糖。术语组织包括（在其意义内）单个的细胞W及细胞的聚合（诸如细胞膜或器官）。该 术语指的是正常的（健康的）或异常的（病态的）细胞或细胞聚合。组织的例子是屯、肌 组织（包括屯、室肌细胞和屯、肌细胞）、膜组织（诸如内皮和上皮）W及结缔组织；病态组 织的例子是梗塞屯、脏组织、血凝块、动脉粥样硬化斑块、炎性组织和肿瘤组织。受体包括位 于组织上（例如，在细胞内或者在它们的表面上）的、能够选择性地绑定到特定物质的任 何分子结构。示例性受体是糖蛋白GPnbllla或纤维蛋白（例如，位于血凝块或血栓上）、 P选择素（例如，位于发炎组织的激活内皮上）或KDR(例如，位于肿瘤组织）上。W下 文献中描述了合适的目标特定的造影剂和目标特定的配体的例子；"化chon等人，BR55:A lipopeptide-basedVEGFR2-targetedultrasoundcontrastagentformolecular imagingofangiogenesis.InvesRadiol2010;45:89-95"和"Tardy等人，叫trasound molecularimagingofVEGFR2inaratprostatetumormodelusingBR55,Inves Radiol2010 ;45:573-578" ;W及WO-A-2006018433 (该些文献的全部公开内容都通过引用 并入本文）。
[0029] 可W通过应用超声波并且记录响应于该些超声波返回的回声信号来检测造影剂； 因为造影剂在患者155中W与红血球相同的速率流动，所W其在身体部分150中的检测和 跟踪提供关于相应的血液灌注的信息。
[0030] 例如，造影剂被W大剂量一一即，在（大约2-20秒的）短时间段内用手通过注射 器提供的一次剂量一一静脉注射到患者 155。造影剂在患者155的血管系统内循环，W便灌 注身体部分150。同时，成像探头110被放置为在身体部分150的区域中与患者155的皮 肤接触。超声脉冲的成像帖序列施加于身体部分150 ;在分析（时间）段期间随着时间的 过去在连续的采集时刻（例如，在2-25秒期间W4-20HZ的帖率）响应于不同成像帖的超 声脉冲记录的回声信号提供身体部分150在该分析时间段期间的相应位置的表示（在其在 选定的扫描平面中的切片中）。回声信号由造影剂（如果存在）所产生的不同贡献和周围 组织的叠加得到。优选地，超声扫描仪100W对比度特定的成像模式操作，W便相对于造影 剂的（非线性）贡献大幅减小回声信号中的组织的主导（线性）贡献一一例如，基于谐波 成像化I)、脉冲反向（PI)、功率调制（PM)、或如"Rafter等人，Imagingtechnologiesand techniques,CardiologyClinics22 (2004),pp. 181-197"(该文献的全部公开内容通过引 用并入本文）中所述的其它对比脉冲序列（CP巧技术。然后产生相应的（对比度特定的） 视频图像的序列；对于每个采集时刻，视频图像包括关于表示身体部分的相应位置的每个 可视化元素（即，像素）的数字值；像素值由在该采集时刻关于该位置记录的回声信号的强 度值定义。该样，视频图像序列表示身体部分150在分析时间段期间的灌注的演变。
[0031] 图2中示出了关于根据本公开的实施例的解决方案的应用的体外情况的例子。
[0032] 具体地说，使（磯脂稳定的）微气泡WImm/s的平均流速循环通过仿组织流动体 模内所包含的直径为3mm的管。用SiemensMedicalSystems(商标）的Sequoia超声扫 描仪在(PS模式下在25秒内W4Hz的帖率（即，100个成像帖）采集相应的视频图像序列。 该图中分别用标号205a、20化和205c指示与成像帖No. 25、No. 43和No. 66相应的视频图 像。可W看出，管的上壁处的微气泡因为层流条件，基本上是静止的一一例如，对于用标号 210s表示的区域可见；相反，管的中屯、（与其上壁间隔开）处的微气泡W平均流速移动一一 例如，对于用标号210m表示的区域可见。
[0033] 在根据本公开的实施例的解决方案中，如W下详细描述的，分割（在该种情况下， 在静止微气泡的像素和移动微气泡的像素之间）基于像素值随着时间的统计分析。
[0034] 具体地说，图3A-图3B中示出了根据本公开的实施例的解决方案的应用的例子。
[003引从图3A开始，在纵轴上绘制像素值（任意单元，a.u.)对比在横轴上绘制帖数 (adimensional)的常见示图中，关于静止微气泡的区域和移动微气泡的区域中的代表性像 素的（线性化）像素值分别用曲线305s(实线）和曲线305m(点线）表示。
[0036] 可W看出，关于静止微气泡的像素（曲线305s)和移动微气泡的像素（曲线305m) 的像素值随着时间的动态大不相同；具体地说，关于静止微气泡的曲线305s几乎是恒定 的，而关于移动微气泡的曲线305m表现出高变化（由于微气泡的出现和消失）。因此，每个 像素的像素值随着时间的统计分布的统计参数可W用于区分它们。
[0037]例如，对于每个像素，计算其统计分布的（算术）均值、中间值和标准差。
[0038] 具体地说，具有离散值Xi(其中i=1…脚的统计变量的均值y(即，在该种情况 下是随时间的像素值）等于：
[0039]
[0040] 均值X是表示当统计变量具有对称形状的统计分布时该统计变量的期望值的集 中趋势统计参数（即，指示统计变量值围绕其聚集的值）。
[0041] 同一统计变量Xi的中间值y(从其值的最小一个Xi到最大一个Xw排序）等于：
[0044] 中间值y是表示当统计变量具有歪斜（即，不对称）形状的统计分布时该统计变 量的期望值的另一个集中趋势统计参数。
[0045] 同一统计变量Xi的标准差0等于：
[0046]
[0047] 标准差0是表示统计变量值如何从它们的期望值离散的离差统计参数（即，指示 统计变量值的变化性或伸展性）。
[004引现在参照图3B，该些统计参数可W用于分割图2的视频图像，W便鉴别静止微气 泡和移动微气泡。为了该个目的，计算=个不同的参数化图像305y、305d和305 0。具体 地说，参数化图像305y、305d和305 0具有与视频图像相同的大小；参数化图像305U、 305d和305 0的每个像素具有参数值，该参数值分别等于相应像素值随着时间的均值、中 间值和标准差的值。
[0049] 可W看出，静止微气泡的像素和移动微气泡的像素之间的对比度（或差异）（例 如，对于分别用标号310s和310m表示的区域可见）相当大地改进。具体地说，在关于均值 的参数化图像305U中W及在关于中间值的参数化图像305d中，参数值在静止微气泡（区 域310s)的像素中较高，而在移动微气泡（区域310m)的像素中较低；该是由于如下的事 实而导致的，即，像素值的期望值（用均值和中间值表达）在静止微气泡的像素中高于移动 微气泡的像素中。对比度在关于中间值的参数化图像305d中高于关于均值的参数化图像 305y中，因为在移动微气泡的像素中，它们的像素值的离群值（将均值拉向其值）不影响 中间值。相反，在关于标准差的参数化图像305 0中，参数值在静止微气泡（区域310s)的 像素中较低，而在移动微气泡（区域310m)的像素中较高；该是由于如下的事实而导致的， 即，像素值在静止微气泡的像素中随着时间的离差小于移动微气泡的像素中。例如，可W通 过计算相应像素的参数值的平均值之间的比率来测量静止微气泡和移动微气泡之间的对 比度。具体地说，该比率在关于均值的参数化图像305y中等于1. 69,在关于中间值的参数 化图像305d中等于2. 72,在关于标准差的参数化图像305 0中等于0. 84 ;该确认提供最好 结果的统计参数是中间值（之后是均值和标准差）。
[0050] W上提及的参数值（尤其是中间值）可W用于用阔值化技术分割原始图像。例 如，用于参数值的分割阔值可W被设置为等于静止微气泡的像素的参数值的平均值和移动 微气泡的像素的参数值的平均值之间的均值（手动选择）；然后，当每个像素的参数值高于 分割阔值时，将该像素分配给静止微气泡的部分，或者否则，将该像素分配给移动微气泡的 部分。
[0051]一般地，上述解决方案方便分割，有益效果是所获得的结果的准确度。
[0052]具体地说，该解决方案可W用于改进DCE-US成像应用中的量化表征的性能。例 如，该使得可W排除（或至少大幅减小）将被分析的包括病变的分析区域中的低灌注部分 (例如，具有坏死组织）的贡献和/或包括将被与病变进行比较的健康组织的控制区域中的 高灌注部分（例如，具有大血管）的贡献。结果，在分析区域的相对合并值中可能反映的相 应误差大幅减小（使得即使当分析/控制区域由操作者手动选择时，分割的可再现性也更 高，因为其质量现在不太取决于个人技能和经验）。
[0053]而且，相同的解决方案可W在UMI应用中用于改进目标特定的造影剂的固定颗粒 的检测和量化。例如，该使得可W移除（或至少大幅减小）目标特定的造影剂的循环颗粒 的贡献，尤其是在其施用之后的早期。
[0054]W上提及的优点在治疗随访的上下文中尤其感受得到，因为它们方便识别病变在 其处理期间的灌注中的即使是微小的变化，该些微小的变化可能指示处理的有效性。
[0055] 作为进一步的改进，分割可W基于与每个像素随着时间的像素值相关的两个或更 多个统计参数。
[0056]具体地说，图4A-图4D中示出了根据本公开的另一实施例的解决方案的应用的例 子。
[0057] 从图4A开始，在横轴上绘制均值并且在纵轴上绘制标准差值（adimensional)的 常见的散布图中，图3B的静止微气泡的区域（310s)和移动微气泡的区域（310m)中的像素 的均值和标准差值分别用点405AS(圆形）和点405Am(十字形）表示；具体地说，用标号 405AS(黑色实屯、圆）和405Am(黑色实屯、星）表示的点分别对应于图2的静止微气泡的区 域（210s)中的代表性像素和移动微气泡的区域（210m)中的代表性像素。
[0058] 计算用于均值的分割阔值THy= 97和用于标准差值的分割阔值THO= 57 (例 如，通过求取点405AS和405Am的相应参数值的平均值）。分割阔值THy和THo在散布图 中定义四个分割象限。具体地说，右上分割象限QAur由分割标准（均值〉THy和标准差值 〉TH0 )定义，左上分割象限QAul由分割标准（均值《THy和标准差值〉TH0 )定义，左下 分割象限QA11由分割标准（均值《THy和标准差值《TH0 )定义，右下分割象限QAlr由 分割标准（均值〉THy和标准差值《THo)定义。
[0059]同样地，在图4B中，在横轴上绘制均值并且在纵轴上绘制标准差值 (adimensional)的另一常见的散布图中，图3B的静止微气泡的区域（310s)和移动微气泡 的区域（310m)中的像素的均值和标准差值分别用点405BS(圆形）和点405Bm(十字形） 表不；具体地说，用t不号405Bs(黑色实屯、圆）和405Bm(黑色实屯、星）表不的点分别对应于 图2的静止微气泡的区域（210s)中的代表性像素和移动微气泡的区域（210m)中的代表性 像素。
[0060] 除了相同的用于标准差值的分割阔值THO= 57之外，还如上那样计算用于中间 值的分割阔值THd= 79。分割阔值THd和THO在散布图中定义其它四个分割象限。具体 地说，右上分割象限地ur由分割标准（均值〉THd和标准差值〉TH0 )定义，左上分割象限 地ul由分割标准（均值《THd和标准差值〉THO)定义，左下分割象限地11由分割标准（均 值《THd和标准差值《TH0 )定义，右下分割象限地Ir由分割标准（均值〉THd和标准差 值《THo)定义。
[0061] 通过比较图4A-图4B的散布图可W看出，静止微气泡的像素和移动微气泡的像素 之间的对比度在关于中间值/标准差的图4B的散布图中高于关于均值/标准差的图4A的 散布图中（例如，对于图4B中的点405Bm、405Bs相对于图4A中的点405Am、405As显而易 见）；该确认与标准差组合提供更好结果的统计参数再次是中间值。
[0062] W上参数值可W用于根据相应的分割标准来分割原始图像。
[0063] 例如，如图4C中所示，四个分割掩模MAur、MAul、Mll和Mir使用关于均值/标准 差的图4A的分别与分割象限QAur、QAul、QA11和QAlr相应的分割标准而构建。每个分割 掩模MAur、MAul、Mll和Mir具有与视频图像相同的大小；具有满足相应的分割标准的参 数值的像素是白色，而其它像素是黑色。而且，在所有的分割掩模MAur、MAul、MAll和MAlr 中，区域410表示用于采集视频图像的流动体模的3mm管（由表示其壁的两个虚线定界）。
[0064] 可W看出，分割掩模MAur和MAlr中的白色像素检测管410的上壁处的静止微气 泡的像素的一部分（即使管410的中屯、处的一些白色像素存在于分割掩模MAur中）；另 一方面，分割掩模MAul和Mil中的白色像素 检测管410的中屯、处的移动微气泡的像素的 一部分（分割掩模Mil中的管410外部的白色像素可能是由于噪声和体模组织信号导致 的）。
[006引 同样地，在图4D中，其它四个分割掩模MBur、MBul、MBll和MBlr如上使用关于中 间值/标准差的图4B的分别与分割象限地ur、地ul、地11和地Ir相应的分割标准而构建 (管的同一区域再次用标号410标识）。如上，分割掩模MBur和MBlr中的白色像素检测管 410的上壁处的静止微气泡的像素的一部分（即使管410的中屯、处的一些白色像素存在于 (但是程度较低）分割掩模MAur中）；另一方面，分割掩模MBul和MB11中的白色像素检测 管410的中屯、处的移动微气泡的像素的一部分（分割掩模MB11中的管410外部的白色像 素再次可能是由于噪声和体模组织信号导致的）。
[0066] 通过比较图4C-图4D的分割掩模可W看出，静止微气泡的像素和移动微气泡的像 素之间的分割在关于中间值/标准差的图4D的分割掩模中好于关于均值/标准差的图4C 的分割掩模中--例如，对于图4D中的分割掩模MBur相对于图4C中的分割掩模MAur是 显而易见的；该再次确认与标准差组合提供更好结果的统计参数是中间值。
[0067] 作为更进一步的改进，分割可W基于W上提及的与每个像素随着时间的像素值相 关的统计参数，该些统计参数被变换W提供具有歪斜形状的、（变换的）像素值的（变换 的）统计分布。
[006引具体地说，图5A-图5C中示出了根据本公开的另一实施例的解决方案的应用的例 子。
[0069] 从图5A开始，假定（变换的）统计分布是对数正态统计分布。在对数正态统计分 布中，统计变量的自然对数具有正态统计分布。具体地说，在连续的统计变量X的情况下， 表示（连续的）正态统计分布的统计密度函数是正态统计密度函数（或简称为正态函数）：
[0070]
[0071] 其中，相同的符号义和0分别用于指示统计变量X的均值和标准差；因此，表示 (连续的）对数正态统计分布的统计密度函数是对数正态统计密度函数（或简称为对数正 态函数）；
[0072]
[0073] 其中，m和S现在分别是统计变量X的自然对数（即，InOO)的均值和标准差。
[0074] 在该种情况下，通过应用W下变换公式，可W从（正态）均值支和（正态）中间值 X直接计算（对数正态）均值m和（对数正态）标准差S:
[0078] 分割现在基于由对数正态均值和对数正态标准差定义的参数值。为了该个目的， 通过将W上提及的变换公式应用于每个像素的随着时间的像素值的实际的均值和中间值 来计算该像素的对数正态均值和对数正态标准差值。在横轴上绘制对数正态均值并且在 纵轴上绘制对数正态标准差值的常见的散布图（adimensional)中，静止微气泡的区域 (310s)和移动微气泡的区域（310m)中的像素的对数正态均值和对数正态标准差值分别 用点505s(圆形）和点505m(十字形）表示；具体地说，用标号505s'(黑色实屯、圆）和 505m'(黑色实屯、星）表示的点分别对应于图2的静止微气泡的区域（210s)中的代表性像 素和移动微气泡的区域（210m)中的代表性像素。
[0079] 如上计算用于对数正态均值的分割阔值T血=4. 4和用于对数正态标准差值的分 割阔值THs= 0. 54。分割阔值T血和ras在散布图中定义四个分割象限。具体地说，右上 分割象限QLur由分割标准（对数正态均值〉T血和对数正态标准差值〉THs)定义，左上分 割象限QLul由分割标准（对数正态均值《T血和对数正态标准差值〉THs)定义，左下象限 QL11由分割标准（对数正态均值《T血和对数正态标准差值《THs)定义，右下象限化Ir 由分割标准（对数正态均值〉T血和对数正态标准差值《THs)定义。
[0080] 可W看出，静止微气泡的像素和移动微气泡的像素之间的对比度进一步改进；实 际上，在该种情况下，静止微气泡和移动微气泡的点聚集在散布图的不同区域中，重叠减 小。该进一步方便分割，因为分割阔值T血、ras现在可W独立地定义，而不会损害分割的性 能。
[0081] 移至图5B，四个分割掩模MLur、MLul、ML11和MLlr如上使用图5A的分别与分割 象限化ur、QLul、化11和化Ir相应的分割标准而构建（管的同一区域再次用标号410标 识）。
[0082] 在该种情况下，分割掩模MLlr中的白色像素检测管410的上壁处的静止微气泡的 大多数像素；另一方面，分割掩模MLul中的白色像素检测管410的中屯、处的移动微气泡的 大多数像素（掩模MLul中的管410外部的白色像素可能是由于CPS的声学伪像与管410 中的微气泡的高浓度组合而导致的，所述声波伪像可W通过降低微气泡的浓度来减小）。相 反，分割掩模MLur和ML11的大多数像素是黑色的（分割掩模MB11中的管410外部的白色 像素再次可能是由于噪声和体模组织信号而导致的）。该确认像素值的对数变换改进了分 割性能。
[0083] 如图5C中所示，W上定义的分割掩模可W应用于图3B的参数化图像305y(表示 每个像素的均值）。具体地说，四个分割图像SIur、SIul、SIll和Sllr通过分别将图5B的 分割掩模MLur、MLul、ML11和MLlr应用于该参数化图像上而构建；分割图像Slur、SIul、 Sill和Sllr的每个像素被分配其参数值，即，如果它在相应的分割掩模中为白色，则被分 配均值（为了显示的目的，用压缩因子40地进行对数压缩），否则，它保持黑色（W便仅显 示分割像素的均值）。结果，分割图像Sllr主要示出管410的上壁处的静止微气泡的像素， 而分割图像SIul主要示出管410的中屯、处的移动微气泡的像素（W便方便区分它们）。
[0084] 图6中示出了表示可W用于实现根据本公开的实施例的解决方案的主要组件的 作用的协作图。
[0085] 具体地说，所有的组件作为整体用标号600表示；该图描述了组件600的静态结构 及其动态行为（通过一系列交换的消息，每个交换的消息表示相应的动作，该动作用前面 为符号"A"的序号表示）。
[0086] 信息（程序和数据）通常被存储在硬盘中，并且当程序连同操作系统和其它应用 程序（该图中未示出）一起运行时，被（至少部分）加载到超声扫描仪的工作存储器中。程 序最初例如从移动盘安装到硬盘上。在该方面，组件600中的每个可W表示包括用于实现 指定的一个（或多个）逻辑函数的一个或多个可执行指令的模块、代码段或代码部分。
[0087] 更具体地说，采集器603用于采集视频图像。为了该个目的，采集器603包括控制 成像探头的驱动器。例如，成像探头驱动器设有发射波束形成器和脉冲器，它们用于在每个 采集时刻产生将施加于所分析的身体部分的超声脉冲的成像帖；超声脉冲具有低声能（诸 如具有机械指数MI=0.01-0. 3)，W便引入造影剂的可忽略的破坏（诸如小于其在连续的 成像帖之间的局部浓度的5%，优选地，小于1 % )。成像探头然后在每个采集时刻接收身体 部分的不同位置（在选定的扫描平面中）所反射的（模拟RF)回声信号。该些模拟RF回 声信号被供给接收处理器，该接收处理器预先放大模拟RF回声信号，并且施加初步时域补 偿（TGC);模拟RF回声信号然后被模数转换器（ADC)转换为数字值，并且被接收波束形成 器组合为聚焦的波束信号。如此获得的数字RF回声信号优选地通过进一步的数字算法和 其它的线性或非线性信号调理器（例如，后波束形成TGC)进行处理；具体地说，接收处理器 应用对比度特定的算法（诸如基于(PS技术）来减小组织的贡献。数字RF回声信号然后 被解调、对数压缩（为了获得具有很好平衡的对比度的图像）、扫描转换为视频格式并且根 据给定的色表调色盘擅染。该样，产生每个采集时刻的视频图像。该视频图像由单元矩阵 定义，每个单元存储相应的像素值。像素值定义像素的亮度；例如，在灰阶视频图像中，像素 值可W在8个比特上编码，作为相应的回声信号的强度的函数，从0 (黑色）增大到255 (白 色）（表示身体部分的相应位置在该采集时刻的声学响应）。
[0088] 在分析处理开始时，超声扫描仪的操作者启动成像探头，并且围绕将被分析的身 体部分移动它（在施用任何造影剂之前）。相应的视频图像一被采集，就被提供给显示器 606 ;显示器606包括驱动器，该驱动器控制超声扫描仪的监视器，W便使得实时地显示该 些视频图像（动作"A1.显示"）。
[0089] 操作者选择表示身体部分的特定切片的扫描平面（例如，包括将被分析的病变）， 并且将成像探头保持固定在该位置上一段预定义的分析时间段。造影剂现在被施用于患 者，并且相应的视频图像序列（表示身体部分在分析时间段期间在所选扫描平面中的灌 注）被如上那样采集。该些视频图像（除了被显示器606实时显示之外，在该图中为了清 楚起见未示出显示器606的连接）还通过将它们保存到贬存器609而被登记（动作"A2. 登记"）。
[0090] 选择器612被操作者用于在贬存器609的一个任意选择的视频图像中选择身体部 分的将被分析的期望的感兴趣区域（R0I)。感兴趣区域用R0I掩模表示，R0I掩模包括具有 与视频图像相同大小的单元的矩阵；R0I掩模的每个单元存储R0I标志（即，二进制值），当 相应像素在感兴趣区域的内部（即，它已经被选择）时，该R0I标志被断言（例如，为逻辑 值1)，否则，该R0I标志被解除断言（例如，为逻辑值0)。R0I掩模被保存到贬存器615中 (动作"A3.选择"）。
[0091] 预处理器618访问贬存器609中的视频图像和贬存器615中的R0I掩模。对于每 个视频图像，预处理器618对其在R0I掩模中的R0I标志被断言（即，在感兴趣区域的内部） 的每个单元的像素值进行处理，W便使得它与造影剂的相应的局部浓度成正比。例如，该结 果可W通过下述方式实现，即，应用逆序调色板-照度表和逆序对数压缩（W使采集器603 应用它们的效果逆转），然后求取如此获得的值的平方值（如W0-A-2004/110279中所述那 样，该申请的全部公开内容通过引用并入本文）。预处理器618然后创建像素值矢量图，该 像素值矢量图包括具有与视频图像相同大小的单元的矩阵；像素值矢量图的其在R0I掩模 中的R0I标志被断言的每个单元存储像素值矢量，该像素值矢量包括沿着视频图像的相应 的线性化像素值的序列（其长度对应于分析时间段的持续时间）。该像素值矢量图被保存 到贬存器621中（动作"A4.预处理"）。
[0092] 贬存器621中的像素值矢量图被供给计算器624,计算器624也访问贬存器 615(在该图中为了清楚起见未示出其连接）中的R0I掩模。对于像素值矢量图的其在R0I 掩模中的R0I标志被断言（即，在感兴趣区域的内部）的每个单元，计算器624计算相应的 (线性化的）像素值的序列的均值和中间值的值。计算器624然后创建关于均值的参数化 图像和关于中间值的参数化图像，每个图像包括具有与视频图像相同大小的单元的矩阵； 关于均值的参数化图像和关于中间值的参数化图像的其在ROI掩模中的ROI标志被断言的 每个单元分别存储相应的均值和中间值。该些参数化图像被保存到贬存器627中（动作 "A5.计算"）。
[0093] 贬存器627中的参数化图像被供给变换器630,变换器630也访问贬存器615(在 该图中为了清楚起见未示出其连接） 中的R0I掩模。对于参数化图像的其在R0I掩模中的 R0I标志被断言（即，在感兴趣区域的内部）的每个像素，变换器630通过应用相应的变换 公式从相应的均值和中间值计算对数正态均值和对数正态标准差的值。变换器630然后创 建统计图，该统计图包括与视频图像相同大小的单元的矩阵；该统计图的其在R0I掩模中 的R0I标志被断言的每个单元存储相应的统计参数（即，对数正态均值和对数正态标准差 值）。该统计图被保存到贬存器633中（动作"A6.变换"）。
[0094] 贬存器633中的统计图被供给分割器636,分割器636也访问贬存器615(在该图 中为了清楚起见未示出其连接）中的R0I掩模。分割器636从贬存器639提取选定的一个 (或多个）分割标准，贬存器639存储若干个基于用于对数正态均值和对数正态标准差值的 分割阔值的预定分割标准；例如，该些分割阔值可能已经在学习阶段期间通过分析具有公 知特性的样品身体部分而被设置。对于统计图的其在R0I掩模中的R0I标志被断言（即， 在感兴趣区域的内部）的每个像素，分割器636验证相应的统计参数是否满足选定的分割 标准（意指它是否属于表示身体部分的具有相应的同质特性的部分的相应部分）。分割器 636然后创建用于选定的分割标准的分割掩模，该分割掩模包括具有与视频图像相同大小 的单元的矩阵；分割掩模的每个单元存储分割标志（即，二进制值），当满足相应的分割标 准时，该分割标志被断言（目P，相应的像素属于该部分），否则，该分割标志被解除断言。如 此创建的每个分割掩模被保存到贬存器642中（动作"A7.分割"）。
[0095] 施加器645从贬存器642提取选定的一个（或多个）分割掩模；施加器645还从 贬存器627提取选定的一个（或多个）参数化图像，和/或从贬存器609提取选定的一个 (或多个）视频图像。施加器645通过将它们逐个单元地相乘来将（选定的）分割掩模施 加于（选定的）参数化/视频图像。该操作创建分割图像，该分割图像包括与视频/参数 化图像相同大小的单元的矩阵；分割图像的属于由分割掩模限定的部分（目P，其分割标志 被断言）的每个单元存储相应的参数化/像素值，而其它单元则被重置为0。如此创建的每 个分割图像被保存到贬存器648中（动作"A8.施加"）。
[0096] 贬存器648中的每个分割图像现在可W被传递给显示器606W供其显示W及任何 下一次使用（动作"A9.显示")。例如，区分目标特定的造影剂的固定颗粒与其余的循环颗 粒的分割图像可W用于检测固定颗粒，因此识别身体部分中的相应病变。另外地或者可替 代地，合并器651可W从贬存器642提取选定的一个（或多个）分割标准，和/或从贬存器 648提取选定的一个（或多个）分割图像。合并器651计算身体部分的由（选定的）分割 掩模/分割图像限定的部分所表示的部分的一个（或多个）合并值；例如，合并值可W被设 置为等于在分割掩模中断言的分割标志的数量（测量身体部分的所述部分），或者合并值 可W被设置为等于分割图像的像素/参数值的和或平均值（表征身体部分的所述部分）。 如此计算的每个合并值被保存到贬存器654中W供下一次使用（动作"A10.合并"）。例 如，其中低灌注部分被负向区分的分割图像中的与包括病变的分析区域相关的R0I的参数 值的平均值自动地不受它们的分布的影响，并且其中高灌注部分被负向区分的分割图像中 的与包括将被与病变进行比较的健康组织的控制区域相关的ROI的参数值的平均值自动 地不受它们的分布的影响；因此，该些合并值可W组合为用相对于控制区域的合并值的相 对项表达的、分析区域的合并值（例如，等于它们的比率）。而且，在分割掩模中被断言的、 区分目标特定的造影剂的固定颗粒的分割标志的数量W及分割图像中的再次区分目标特 定的造影剂的固定颗粒的像素/参数值的和，该两者都量化该些固定颗粒。
[0097] 图7、图8和图9中示出了根据本公开的实施例的解决方案的体内应用的不同例 子。
[0098] 从图7开始，肝肿瘤的视频图像705是在注入造影剂之后（在峰值造影期间），在 对比度特定的成像模式下20秒采集的；两个（低灌注）区域710a和710b(较暗）在视频 图像705上被手动勾画。
[0099] 分割掩模715是W分割标准（对数正态均值乂. 0和对数正态标准差值〉〇)，从在 视频图像705的采集时刻左右的12. 5秒采集的肝肿瘤的100个视频图像的序列创建的。可 W看出，分割掩模715排除了两个区域720a和72化（黑色），如被排除的区域720a、72化和 在分割掩模715上用灰色虚线报告的低灌注区域710a、71化的轮廓之间的良好重叠所确认 的，区域720a和72化分别基本上对应于低灌注区域710a和710b。
[0100] 移至图8,将用作参考区域的肝脏的健康组织的视频图像805被如上那样采集；与 大血管相应的（高灌注）区域810(较亮）在视频图像805上被手动勾画。
[0101] 分割掩模815是W分割标准（对数正态均值巧.9和对数正态标准差值〉0)，从视 频图像805的采集时刻左右的3秒采集的肝脏的25个视频图像的序列创建的。可W看出， 分割掩模815排除了区域820 (黑色），如被排除区域820和在分割掩模815上用灰色虚线 报告的高灌注区域810的轮廓之间的良好重叠所确认的，区域820基本上对应于高灌注区 域 810。
[0102] 最后参照图9,具有（腺癌）肿瘤的模型鼠的前列腺的视频图像905是在注入针对 血管细胞粘附分子1 (VCAM-1)的目标特定的造影剂之后，在对比度特定的成像模式下2分 钟采集的；肿瘤区域910 (较亮）在视频图像905上被手动勾画。
[0103]分割掩模是W分割标准（对数正态均值M. 4和对数正态标准差值〉0. 54)，从在视 频图像905的采集时刻之后的15秒采集的、前列腺的60个视频图像的序列创建的；然后通 过将该分割掩模施加于从与上相同的视频图像序列创建的、关于均值的参数化图像来创建 分割图像915。可W看出，分割图像915基本上示出了仅肿瘤区域910中的目标特定的造 影剂的固定颗粒，目标特定的造影剂的其余的循环颗粒在视频图像905中存在于肿瘤区域 910的外部，肿瘤区域910在分割图像915中基本上不存在。而且，同样地在肿瘤区域910 中，目标特定的造影剂的其余的循环颗粒基本上被移除（通过将它与视频图像905中的肿 瘤区域910进行比较，显而易见）。
[0104] 作为控制，在相同的条件下，但是现在是使用（非针对性的）造影剂来采集另一视 频图像925,并且如上创建相应的分割图像935。
[0105] 在该种情况下，造影剂的（循环/固定）颗粒在分割图像935中基本上都是不可 见的（两者都在肿瘤区域910的内部和外部）。
[0106] 修改
[0107] 自然地，为了满足局部要求和特定要求，本领域技术人员可W将许多逻辑和/或 物理修改和改变应用于上述解决方案。更具体地说，尽管已经参照该解决方案的一个或多 个实施例一定程度地具体地描述了该解决方案，但是应理解，形式和细节上的各种省略、替 换和改变是可能的。具体地说，本公开的不同实施例甚至可W在没有在前面的描述中所阐 述的、提供其更透彻的理解的特定细节（诸如数值）的情况下实施；相反，公知特征可能已 经被省略或简化W便不使不必要的细节模糊本描述。而且，明确的意图是，与所公开的解决 方案的任何实施例结合描述的特定的元件和/或方法步骤可W作为一般设计选择合并在 任何其它实施例中。在任何情况下，顺序限定词或其它限定词仅用作区分具有相同名称的 元件的标签，而不独自暗示任何优先权、优先级或次序。而且，术语包括、具有、包含和设及 (及其任何形式）的意图应是开放式的、非穷举的意义（即，不限于所记载的项），术语基 于、取决于、根据、…的功能（及其任何形式）的意图应是作为非排他性的关系（即，设及 可能的另外的变量），术语一个的意图应是作为一个或多个项（除非另有明确指示）。
[010引例如，实施例提供一种用于诊断成像应用中的数据处理分割方法。该方法包括W下步骤。提供身体部分（被灌注造影剂）在非零分析时间段期间的表示；对于身体部分的 一组位置中的每个位置，该表示包括该位置在分析时间段期间对询问信号的响应的指示。 对于一组选定位置中的每个选定位置，计算该选定位置在分析时间段期间的响应的统计分 布的至少一个统计参数的值；所述一组选定位置包括所有的位置或者其一部分。根据关于 选定位置的所述至少一个统计参数的值与至少一个分割阔值之间的比较，分割选定位置。
[0109] 然而，身体部分的表示可任何方式（参见下面）提供，W便提供其在具有任 何非零长度（长得足W在统计上显著）的分析时间段期间的表示；该表示可W与任何类型 的灌注（例如，基于任何类型的造影剂（非针对性或针对性类型）的任何类型的灌注（例 如，基于通过破坏飞瓣的连续注射））相关。所述方法可W应用于任何类型的位置的级别上 (例如，像素、体素或者其群组）；具体地说，还可W对所述表示在空间上进行二次采样，确 定二次采样的表示上的分割，然后通过插值技术恢复分割的原尺寸W使其应用于原始表示 上。而且，所述方法可W应用于任何选定位置（例如，在W任何方式定义的一个或多个感兴 趣区域中）、或者甚至身体部分的整个表示。统计参数可W是任何类型和任何数量（参见下 面）；然而，可W根据统计参数的值与任何类型和数量的分割阔值之间的任何比较将选定 位置分割为任何数量的部分（两个或更多个）（参见下面）。
[0110] 在实施例中，所述提供身体部分在非零分析时间段期间的表示的步骤包括：提供 在所述分析时间段期间登记的身体部分的数字图像序列；每个数字图像包括关于每个位置 的强度值，该强度值基于该位置在分析时间段中的相应采集时刻的响应。所述计算至少一 个统计参数的值的步骤包括从沿着数字图像的选定位置的强度值计算所述至少一个统计 参数的值。
[0111] 然而，数字图像可WW任何方式（例如，通过对它们进行预处理W补偿运动伪像 或减去背景图像）登记。数字图像可W是任何类型（例如，单独地或组合地、3D类型、基于 任何对比度特定的成像模式、或者甚至基于基板B模式）、任何大小和任何非零数量（W便 在统计上显著）；强度值可W是任何类型（例如，对数压缩值）并且W任何方式编码（例如， 对任何数量的比特进行颜色编码）。在任何情况下，不排除将同一方法应用于身体部分的任 何表示（甚至不是图像形式）的可能性。例如，身体部分在非零分析时间段期间的表示可 W由原始（未处理RF)回声信号定义。在该种情况下，对于每个选定位置，（响应的统计分 布的）所述至少一个统计参数的值通过下述方式获得，即，计算选定位置的原始回声信号 的HiAed变换，然后在其幅值或相位的级别上计算（复数）相关系数或者其方差（然后 如上那样，根据该些统计参数的值和至少一个分割阔值之间的比较来分割选定位置）。
[0112] 在实施例中，所述至少一个统计参数是多个统计参数。
[0113] 然而，统计参数可W是任何数量（下至一个，或者相反，=个或更多个）。
[0114] 在实施例中，统计参数是集中趋势统计参数和离差统计参数。
[0115] 然而，可W使用任何集中趋势统计参数（例如，除了（算术）均值和中间值、模式、 几何均值、谐波均值、截断均值、中间范围均值等之外）和/或任何离差统计参数（例如，除 了标准差、均值差值、平均偏差、方差等之外）。在任何情况下，可W使用不同的、附加的或替 代的类型的统 计参数（例如，位置或形状统计参数）。
[0116] 在实施例中，所述计算至少一个统计参数的值的步骤包括计算与响应相关的所述 至少一个统计参数的值，该响应被变换W提供变换的响应在分析时间段期间的、具有歪斜 形状的变换的统计分布。
[0117] 然而，该计算可W根据任何变换的统计分布W任何方式实现（参见下面）。在任何 情况下，不排除仅从响应的（原始）统计分布计算统计参数的可能性。
[0118] 在实施例中，变换的统计分布是对数正态统计分布。
[0119] 然而，可W使用任何其它的具有歪斜形状的变换的统计分布（例如，伽马或局部 密度随机游走分布）。
[0120] 在实施例中，所述计算至少一个统计参数的值的步骤包括；从选定位置在分析时 间段期间的响应的统计分布计算至少一个另一统计参数的值，并且通过将根据变换的统计 分布的相应变换应用于所述至少一个另一统计参数的值来计算每个统计参数的值。
[0121] 然而，变换可W是任何类型的（参见下面）。在任何情况下，没有什么能阻止变换 响应、然后从变换的响应直接计算期望的统计参数；例如，可W如下计算对数正态均值和对 数正态标准差：
[0125] 如果Xi<l，则s= 0
[0126] 可替代地，还可W通过参数化函数（直接地或者从其直方图）拟合每个响应、然后 从该参数化函数计算期望的统计参数，该参数化函数是表示变换的统计分布的分布函数的 参数化模型的实例。
[0127] 在实施例中，所述通过应用相应变换计算每个统计参数的值的步骤包括将对数变 换应用于所述至少一个另一统计参数的值。
[012引然而，变换可W是任何类型的（例如，指数类型）。
[0129] 在实施例中，所述至少一个统计参数是对数正态均值和对数正态标准差，并且所 述至少一个另一统计参数是正态均值和正态中间值。
[0130] 然而，所述统计参数和/或另一统计参数可W是任何类型和任何数量的（例如，如 w上所指出的）。
[0131] 在实施例中，所述通过应用相应变换来计算每个统计参数的值的步骤包括：设置 对数正态均值的值等于正态中间值的值的自然对数，并且如果正态均值的值和正态中间值 的值之间的比率严格地高于1，则设置对数正态标准差的值等于该比率的自然对数的两倍 的平方根，否则，等于0。
[0132] 然而，可W根据变换的统计分布和/或统计参数的类型使用其它变换。
[0133] 在实施例中，所述分割选定位置的步骤包括根据每个统计参数的值和所述至少一 个分割阔值中的相应的一个分割阔值之间的比较来分割选定位置。
[0134] 然而，分割阔值可任何方式定义。例如，可W应用主成分分析（PCA)技术来解 除统计参数的值的相关性；而且，分割阔值可W从相应的统计参数的值的重屯、计算。可W根 据每个统计参数和相应的分割阔值（例如，就绝对值而言或者根据预定义百分比或者任何 组合，较高或较低）之间的任何比较来分割选定位置。在任何情况下，没有什么能阻止使用 与统计参数之一不同的若干个分割阔值（例如，当分割基于统计参数的任何（线性或非线 性）组合时）。
[0135] 在实施例中，所述分割选定位置的步骤包括；当集中趋势统计参数的值高于相应 的分割阔值并且离差统计参数的值高于相应的分割阔值时，将身体部分的分析区域（包括 将被分析的病变）的每个选定位置分配给第一部分（与分析区域的低灌注相关），否则，分 配给第二部分；所述方法还包括计算分析区域的不包括在第一部分中的选定位置的响应的 合并值。
[0136] 然而，该分割可W基于其它分割标准（例如，原始回声信号的化化6的变换的高相 关系数和低方差）；而且，它可W用于计算分析区域的任何合并值（例如，它们到达峰值的 平均时间、歪斜等）。
[0137] 在实施例中，所述分割选定位置的步骤包括；当集中趋势统计参数的值低于相应 的分割阔值并且离差统计参数的值高于相应的分割阔值时，将身体部分的控制区域（包括 将被与分析区域进行比较的健康组织）的每个选定位置分配给另一第一部分（与控制区域 的高灌注部分相关），否则，分配给另一第二部分；所述方法还包括计算控制区域的不包括 在所述另一第一部分中的选定位置的响应的另一合并值。
[013引然而，该分割可W基于其它分割标准（例如，原始回声信号的化化6的变换的低相 关系数和高方差）；然而，它可W用于如上那样计算控制区域的任何合并值。
[0139] 在实施例中，所述方法还包括根据所述合并值和另一合并值的组合来计算分析区 域相对于控制区域的相对合并值。
[0140] 然而，相对合并值可任何方式计算（例如，按照所述合并值和另一合并值之 间的差值计算）；在任何情况下，可任何其它方式使用所述另一合并值（例如，用于规 范化整个参数化图像的参数值）。
[0141] 在实施例中，造影剂是能够在所述患者内循环并且能够基本固定在生物目标上的 针对性的造影剂；所述分割选定位置的步骤包括；当集中趋势统计参数的值高于相应的分 割阔值并且离差统计参数的值低于相应的分割阔值时，将每个选定位置分配给第一部分 (与固定的针对性的造影剂相关），否则，分配给第二部分。
[0142] 然而，针对性的造影剂可W是任何类型（例如，基于特定的或非特定的相互作 用）；而且，该分割可W基于其它分割标准（例如，原始回声信号的HUbed变换的高相关 系数和低方差）。
[0143] 在任何情况下，相同的技术可W用于任何其它诊断成像应用中；例如，相同的过程 可W对选定的部分重复迭代（一次或多次），或者它可W重复地W不同的分割标准应用，其 中，所得的分割掩模通过逻辑运算符组合，诸如W获得相应部分的交集（逻辑AND)、并集 (逻辑OR)或异并狂OR)。
[0144] 一般地，如果相同的解决方案用等同的方法（通过使用其多个步骤或部分的具有 相同功能的类似步骤、移除不必要的一些步骤、或者添加另外的操作步骤）实现，则类似的 考虑因素适用；而且，该些步骤可W按不同次序、同时地或者（至少部分）W交错的方式执 行。
[0145] 在任何情况下，强调的是，上述方法是可W独立于与患者（W及特别是与可W在 执行该方法之前预先施用于患者的造影剂）的任何相互作用实现的数据处理（或计算）方 法。而且，造影剂可非侵入性的方式（例如，口服W用于使胃肠道成像或者经由雾化 器进入气道中）施用于患者，或者在对其没有任何实质的物理干预的任何情况下施用于患 者，所述物理干预将需要专业的医疗经验或者承担患者的健康风险（例如，肌肉注射）。在 任何情况下，尽管所提出的方法可W方便医生的任务，但是它一般仅提供中间结果，该些中 间结果可W帮助他/她分析身体部分一一例如，出于诊断的目的（即使用于治愈目的的诊 断在严格意义上总是由医生自己做出）。
[0146]另一实施例提供一种计算机程序，该计算机程序被配置为当该计算机程序在计算 系统上执行时使该计算系统执行W上提及的方法。
[0147] 另一实施例提供一种计算机程序产品，该计算机程序产品包括包含计算机程序的 非暂时性计算机可读介质；该计算机程序可加载到计算系统的工作存储器中，从而将该计 算系统配置为执行相同的方法。
[0148] 然而，相同的解决方案可W实现为独立的模块、用于超声扫描仪的控制程序的插 件、或者甚至直接在控制程序本身中实现；将显而易见的是，还可W将相同的解决方案部署 为通过网络访问（诸如在互联网中）的服务。在任何情况下，如果软件程序（其可W用于 实现本公开的每个实施例）W不同的方式构造，或者如果提供附加模块或功能，则类似的 考虑因素适用；同样地，存储器结构可W是其它类型，或者可W用等同的实体（不一定由物 理存储介质构成）取代。程序可W采取适合于被任何计算（或数据处理）系统使用的或者 与其结合使用的形式（例如，在虚拟机内），从而将系统配置为执行期望的操作；具体地说， 程序可W是外部或驻留软件、固件或微代码（用对象代码或源代码一一例如，将被编译或解 释）的形式。而且，可W在任何计算机可用介质上提供程序（具体地说，作为非暂时性介质 上的制造品）；介质可W是适合于包含、存储、传送、传播或传输程序的任何元件。例如，介 质可W是电子、磁、光、电磁、红外或半导体类型；该样的介质的例子是固定盘（其中程序可 W被预先加载）、可移动盘、带、卡、线、光纤、无线连接、网络、广播波等。在任何情况下，根据 本公开的实施例的解决方案适合于甚至用硬件结构实现（例如，集成在半导体材料的巧片 中），或者用被适当地编程的或者W其它方式配置的软件和硬件的组合实现。
[0149] 另一实施例提供一种包括被构造为执行W上提及的方法的步骤的装置的系统。
[0150] 然而，相同的解决方案可W应用于包括超声扫描仪和不同的计算机的系统（或任 何等同的系统）中；在该种情况下，所记录的信息从超声扫描仪传送到计算机w供对其进 行处理（例如，通过数字、模拟或网络连接）。在任何情况下，不排除将所提出的解决方案 应用于任何其它的诊断成像系统一-例如，基于磁共振成像（MRI)或X射线计算断层扫描 (CT)--的可能性。
[0151] 一般地，如果系统具有不同的结构或者包括等同的组件，或者它具有其它的操作 特性，类似的考虑因素适用。在任何情况下，其每一个组件可W划分为多个元件，或者两个 或更多个组件可W-起组合为单个元件；而且，每个组件可W被复制W支持相应操作的并 行执行。而且，除非另有指定，否则不同组件之间的任何相互作用一般不需要是连续的，并 且它可W是直接的或者通过一个或多个中间物间接的。
[0152] 另一实施例提供一种包括W下步骤的诊断成像方法。对患者施用造影剂W使造影 剂灌注患者的身体部分。在非零分析时间段期间将询问信号施加于身体部分。采集身体部 分在分析时间段期间的表示；对于身体部分的一组位置中的每个位置该表示包括该位置在 分析时间段期间对询问信号的响应的指示（该些表示根据与W上方法相同的方法进行处 理W分割一组选定位置，该一组选定位置包括所有的位置或者其一部分）。根据所述分割评 估身体部分的状况。
[0153] 然而，所述方法可W在任何种类的诊断应用中找到应用（从该术语的最广泛的意 义上来讲一一例如，旨在发现新的病变或监视已知的病变），并且用于分析任何（人类或动 物）患者的任何种类的身体部分（例如，器官（诸如肝脏、前列腺或屯、脏）、区域或组织）。 具体地说，分割可W用于识别身体部分的无论什么生理状况的区域一一例如，区分病变与 健康组织或者身体部分的具有不同的灌注参数的部分。
【主权项】
1. 一种用于诊断成像应用中的数据处理分割方法（Al-AlO)，所述方法包括： 提供（A1-A4)被灌注造影剂的身体部分在非零分析时间段期间的表示，对于所述身体 部分的一组位置中的每个位置所述表示包括所述位置在所述分析时间段期间对询问信号 的响应的指示； 对于一组选定位置中的每个选定位置，计算（A5-A6)所述选定位置在所述分析时间段 期间的响应的统计分布的至少一个统计参数的值，所述一组选定位置包括所有的所述位置 或者其一部分；和 根据关于所述选定位置的所述至少一个统计参数的值与至少一个分割阈值之间的比 较，分割（A7-A9)所述选定位置。2. 根据权利要求1所述的方法（Al-AlO)，其中，所述提供（A1-A4)身体部分在非零分 析时间段期间的表示包括： 提供（A2)在所述分析时间段期间登记的所述身体部分的 数字图像序列，每个数字图 像包括关于每个位置的强度值，所述强度值基于所述位置在所述分析时间段中的相应采集 时刻的响应； 并且其中，所述计算（A5-A6)至少一个统计参数的值包括： 从沿着所述数字图像的选定位置的强度值计算（A5)所述至少一个统计参数的值。3. 根据权利要求1或2所述的方法（Al-AlO)，其中，所述至少一个统计参数是多个统 计参数。4. 根据权利要求3所述的方法（Al-AlO)，其中，所述统计参数是集中趋势统计参数和 离差统计参数。5. 根据权利要求1至4中的任何一个所述的方法（Al-AlO)，其中，所述计算（A5-A6) 至少一个统计参数的值包括： 计算与所述响应相关的所述至少一个统计参数的值，所述响应被变换以提供所述变换 的响应在所述分析时间段期间的、具有歪斜形状的变换的统计分布。6. 根据权利要求5所述的方法（Al-AlO)，其中，所述变换的统计分布是对数正态统计 分布。7. 根据权利要求5或6所述的方法（A1-A10)，其中，所述计算（A5-A6)至少一个统计 参数的值包括： 从所述选定位置在所述分析时间段期间的响应的统计分布计算（A5)至少一个另一统 计参数的值；和 通过将根据所述变换的统计分布的相应变换应用于所述至少一个另一统计参数的值 来计算（A6)每个统计参数的值。8. 根据权利要求7所述的方法（A1-A10)，其中，所述通过应用相应变换来计算（A6)每 个统计参数的值包括： 将对数变换应用（A6)于所述至少一个另一统计参数的值。9. 根据权利要求7或8所述的方法（Al-AlO)，其中，所述至少一个统计参数是对数正 态均值和对数正态标准差，并且所述至少一个另一统计参数是正态均值和正态中间值。10. 根据权利要求9所述的方法（A1-A10)，其中，所述通过应用相应变换来计算（A6) 每个统计参数的值包括： 设置（A6)所述对数正态均值的值等于所述正态中间值的值的自然对数，并且如果所 述正态均值的值和所述正态中间值的值之间的比率严格地高于1，则设置所述对数正态标 准差的值等于所述比率的自然对数的两倍的平方根，否则，等于0。11. 根据权利要求1至10中的任何一个所述的方法（Al-AlO)，其中，所述分割（A7-A9) 选定位置包括： 根据每个统计参数的值和所述至少一个分割阈值中的相应的一个分割阈值之间的比 较，分割（A7)所述选定位置。12. 根据权利要求11所述的方法（Al-AlO)，其中，所述分割（A7-A9)选定位置包括： 当集中趋势统计参数的值高于相应的分割阈值并且离差统计参数的值高于相应的 分割阈值时，将所述身体部分的包括将被分析的病变的分析区域的每个选定位置分配 (A7-A8)给所述分析区域的与低灌注部分相关的第一部分，否则，分配给第二部分； 所述方法还包括： 计算（AlO)所述分析区域的不包括在所述第一部分中的选定位置的响应的合并值。13. 根据权利要求12所述的方法（Al-AlO)，其中，所述分割（A7-A9)选定位置包括： 当所述集中趋势统计参数的值低于相应的分割阈值并且所述离差统计参数的值高于 相应的分割阈值时，将所述身体部分的将被与所述分析区域进行比较的、包括健康组织的 控制区域的每个选定位置分配（A7-A8)给所述控制区域的与高灌注部分相关的另一第一 部分，否则，分配给另一第二部分； 所述方法还包括： 计算（AlO)所述控制区域的不包括在所述另一第一部分中的选定位置的响应的另一 合并值。14. 根据权利要求13所述的方法（Al-AlO)，还包括： 根据所述合并值和所述另一合并值的组合，计算（AlO)所述分析区域相对于所述控制 区域的相对合并值。15. 根据权利要求11所述的方法（Al-AlO)，其中，所述造影剂是能够在所述患者内循 环并且能够基本固定在生物目标上的针对性的造影剂，所述分割（A7-A9)选定位置包括： 当集中趋势统计参数的值高于相应的分割阈值并且离差统计参数的值低于相应的分 割阈值时，将每个选定位置分配（A7-A8)给与固定的针对性的造影剂相关的第一部分，否 贝1J，分配给第二部分。16. -种计算机程序（600)，被配置为当所述计算机程序在计算系统（100)上执行时使 所述计算系统执行根据权利要求1至15中的任何一个所述的方法（Al-AlO)。17. -种计算机程序产品（130)，所述计算机程序产品包括包含计算机程序（600)的非 暂时性计算机可读介质，所述计算机程序可加载到计算系统（100)的工作存储器中，从而 将所述计算系统配置为执行根据权利要求1至15中的任何一个所述的方法（Al-AlO)。18. -种用于诊断成像应用中的系统（100)，所述系统包括： 用于提供被灌注造影剂的身体部分在非零分析时间段期间的表示的装置（603-621)， 对于所述身体部分的一组位置中的每个位置所述表示包括所述位置在所述分析时间段期 间对询问信号的响应的指示； 用于对于一组选定位置中的每个选定位置计算所述选定位置在所述分析时间段期间 的响应的统计分布的至少一个统计参数的值的装置（624-633)，所述一组选定位置包括所 有的所述位置或者其一部分；和 用于根据关于所述选定位置的所述至少一个统计参数的值与至少一个分割阈值之间 的比较分割所述选定位置的装置（636-648)。19. 根据权利要求18所述的系统（100)，其中，所述用于提供身体部分在非零分析时间 段期间的表示的装置（603-621)包括： 用于提供在所述分析时间段期间登记的所述身体部分的数字图像序列的装置（603)， 每个数字图像包括关于每个位置的强度值，所述强度值基于所述位置在所述分析时间段中 的相应采集时刻的响应； 并且其中，所述用于计算至少一个统计参数的值的装置（624-633)包括： 用于从沿着所述数字图像的选定位置的强度值计算所述至少一个统计参数的值的装 置（624)。20. 根据权利要求18或19所述的系统（100)，其中，所述至少一个统计参数是多个统 计参数。21. 根据权利要求20所述的系统（100)，其中，所述统计参数是集中趋势统计参数和离 差统计参数。22. 根据权利要求18至21中的任何一个所述的系统（100)，其中，所述用于计算至少 一个统计参数的值的装置（624-633)包括： 用于计算与所述响应相关的所述至少一个统计参数的值的装置（630)，所述响应被变 换以提供所述变换的响应在所述分析时间段期间的、具有歪斜形状的变换的统计分布。23. 根据权利要求22所述的系统（100)，其中，所述变换的统计分布是对数正态统计分24. 根据权利要求22或23所述的系统（100)，其中，所述用于计算至少一个统计参数 的值的装置（624-633)包括： 用于从所述选定位置在所述分析时间段期间的响应的统计分布计算至少一个另一统 计参数的值的装置（624);和 用于通过将根据所述变换的统计分布的相应变换应用于所述至少一个另一统计参数 的值来计算每个统计参数的值的装置（630)。25. 根据权利要求24所述的系统（100)，其中，所述用于通过应用相应变换来计算每个 统计参数的值的装置（630)包括： 用于将对数变换应用于所述至少一个另一统计参数的值的装置（630)。26. 根据权利要求24或25所述的系统（100)，其中，所述至少一个统计参数是对数正 态均值和对数正态标准差，并且所述至少一个另一统计参数是正态均值和正态中间值。27. 根据权利要求26所述的系统（100)，其中，所述用于通过应用相应变换来计算每个 统计参数的值的装置（630)包括： 用于以下操作的装置（630):设置所述对数正态均值的值等于所述正态中间值的值的 自然对数，并且如果所述正态均值的值和所述正态中间值的值之间的比率严格地高于1，则 设置所述对数正态标准差的值等于所述比率的自然对数的两倍的平方根，否则，等于0。28. 根据权利要求18至27中的任何一个所述的系统（100)，其中，所述用于分割选定 位置的装置(636-648)包括： 用于根据每个统计参数的值和所述至少一个分割阈值中的相应的一个分割阈值之间 的比较分割所述选定位置的装置（636、645)。29. 根据权利要求28所述的系统（100)，其中，所述用于分割选定位置的装置（636、 645)包括： 用于以下操作的装置（645):当集中趋势统计参数的值高于相应的分割阈值并且离差 统计参数的值高于相应的分割阈值时，将所述身体部分的包括将被分析的病变的分析区域 的每个选定位置分配给所述分析区域的与低灌注部分相关的第一部分，否则，分配给第二 部分； 所述系统还包括： 用于计算所述分析区域的不包括在所述第一部分中的选定位置的响应的合并值的装 置（651)。30. 根据权利要求29所述的系统（100)，其中，所述用于分割选定位置的装置（636、 645)包括： 用于以下操作的装置（645):当所述集中趋势统计参数的值低于相应的分割阈值并且 所述离差统计参数的值高于相应的分割阈值时，将所述身体部分的将被与所述分析区域进 行比较的、包括健康组织的控制区域的每个选定位置分配给所述控制区域的与高灌注部分 相关的另一第一部分，否则，分配给另一第二部分； 所述系统还包括： 用于计算所述控制区域的不包括在所述另一第一部分中的选定位置的响应的另一合 并值的装置（651)。31. 根据权利要求30所述的系统（100)，还包括： 用于根据所述合并值和所述另一合并值的组合计算所述分析区域相对于所述控制区 域的相对合并值的装置（651)。32. 根据权利要求28所述的系统（100)，其中，所述造影剂是能够在所述患者内循环 并且能够基本固定在生物目标上的针对性的造影剂，所述用于分割选定位置的装置（630、 645)包括： 用于以下操作的装置（645):当集中趋势统计参数的值高于相应的分割阈值并且离差 统计参数的值低于相应的分割阈值时，将每个选定位置分配给与固定的针对性的造影剂相 关的第一部分，否则，分配给第二部分。33. -种诊断方法，包括： 对患者施用造影剂以使所述造影剂灌注所述患者的身体部分； 在非零分析时间段期间将询问信号应用于所述身体部分； 采集所述身体部分在所述分析时间段期间的表示，其中，对于所述身体部分的一组位 置中的每个位置，所述表示包括所述位置在所述分析时间段期间对所述询问信号的响应的 指示，所述表示根据权利要求1至15中的任何一个所述的方法进行处理以分割一组选定位 置，所述一组选定位置包括所有的所述位置或者其一部分；和 根据所述分割评估所述身体部分的状况。
【专利摘要】提出了一种用于诊断成像应用中的分割方案。相应的数据处理分割方法(A1-A10)包括：提供(A1-A4)被灌注造影剂的身体部分在非零分析时间段期间的表示，对于所述身体部分的一组位置中的每个位置所述表示包括所述位置在所述分析时间段期间对询问信号的响应的指示；对于一组选定位置中的每个选定位置，计算(A5-A6)所述选定位置在所述分析时间段期间的响应的统计分布的至少一个统计参数的值，所述一组选定位置包括所有的所述位置或者其一部分；并且根据关于所述选定位置的所述至少一个统计参数的值与至少一个分割阈值之间的比较，分割(A7-A9)所述选定位置。
【IPC分类】G06T7/00
【公开号】CN104885121
【申请号】CN201380066958
【发明人】P·弗林京
【申请人】博莱科瑞士股份有限公司
【公开日】2015年9月2日
【申请日】2013年12月18日
【公告号】EP2936433A1, US20150348277, WO2014096041A1
【技术领域】
[0001] 根据本公开的方案设及诊断成像领域。更具体地说，本方案设及诊断成像应用中 的分割。
【背景技术】
[0002] 分割技术广泛用于诊断成像应用（或简称为成像应用）中；一般来说，任何分割技 术旨在将所分析的患者的身体部分的表示（例如，其数字图像）分割为两个或更多个部分， 每个部分包括该身体部分的该表示的具有基本上同质的特性的一部分。
[0003] 例如，分割技术可W在动态对比度增强超声值CE-U巧成像应用中用于评估血液 灌注。在该种情况下，充当高效超声反射体的超声造影剂扣CA)或简称为造影剂被施用于 患者。造影剂在患者中W与红血球相同的速率流动，W使得其通过记录响应于超声波返回 的回声信号的检测提供关于相应的血液灌注的信息。具体地说，对于身体部分的其中回声 信号已经被记录的每个位置，（例如，从随着时间拟合回声信号的参数化函数）计算表征它 的灌注参数（诸如峰值增强）的值。通过将相应的灌注参数值分配给表示身体部分的位置 的每个像素，构建参数化图像。该参数化图像可W用于计算身体部分的包括病变的分析区 域中的灌注参数值的合并值（例如，它们的均值）；通常，分析区域的合并值用相对于身体 部分的包括健康问题的控制区域中的灌注参数值的合并值的相对项来表示。该相对的合并 值提供病变的量化表征，该量化表征可W用于例如其治疗随访中。
[0004] 然而，分析区域和控制区域的合并值本质上受一定程度的不精确度的影响，因为 它们是在身体部分的不完美地同质的区域中整体计算的；例如，分析区域可能包括低灌注 的坏死组织，控制区域可能包括高灌注的大血管。分析区域和控制区域的该些低灌注部分 和高灌注部分分别在它们的合并值中引起在分析区域的相对合并值中反映的误差。该缺点 在治疗随访中尤其严重，其中，病变在其处理期间的灌注中即使是微小的变化都可能指示 处理的有效性。
[0005] 可替代地，W0-A-2011/026866(其全部公开内容通过引用并入本文）公开了通过 下述方式确定关于分析区域（限定感兴趣区域）的每个像素的分析函数，即，用时间模型函 数（例如，对数正态分布函数）拟合其回声-功率信号，并且通过用相同的时间模型函数拟 合其回声-功率信号的平均值来确定关于参考区域（包括健康实质）的参考函数。对于分 析区域的每个像素，通过从其分析函数减去参考函数来计算差函数；然后根据像素的差函 数将该像素分配给不同的类（即，正单极类、负单极类、正到负双极类和负到正双极类）。为 了该个目的，可W计算差函数的正能量和负能量；当正能量和负能量的和高于显著性阔值 时，计算相对正能量和相对负能量。如果相对正能量或相对负能量超过鉴别阔值，则将差函 数分别分配给正单极类或负单极类；相反，根据差函数的主要极性变化的数量级将该差函 数分配到正到负双极类或负到正双极类。
[0006] 分割技术还可W用于超声分子成像扣MI)应用中。在该种情况下，使用目标特定 的造影剂；该目标特定的造影剂适于到达特定生物目标（例如，病变），然后通过特定的相 互作用在其上保持固定。目标特定的造影剂的在身体部分中固定的任何颗粒的检测使得可W识别否则将难W发现的相应病变；而且，该些固定颗粒的量化使得可W确定病变的状态， 该状态可W用于例如该病变的治疗随访中。
[0007] 然而，目标特定的造影剂的固定颗粒的检测和量化受到如下事实的阻碍，即，其总 量中只有一小部分实际上到达目标，并且在其上保持固定；大部分目标特定的造影剂而是 继续循环相当长的时间（诸如，几分钟）一一例如，直到它被肺脏滤出和/或在患者的肝脏 中被滤出为止。因此，身体部分中的目标特定的造影剂的颗粒包括固定的颗粒的贡献和仍 然循环的颗粒的贡献两者；因此，区别该两个贡献是相当困难的（尤其是在施用目标特定 的造影剂之后的早期）。
【发明内容】
[0008] 本公开的精简的总结在本文中是为了提供其基本理解而呈现的；然而，本总结的 唯一目的是作为其后面的更详细描述的序言W精简的方式介绍本公开的一些概念，并且本 总结不应被解释为标识其关键元素，也不应被解释为其范围的划定。
[0009] 一般地说，本公开基于使分割基于随着时间的统计分析的构思。
[0010] 具体地说，一方面提供一种用于诊断成像应用中的数据处理分割方法，其中，被灌 注造影剂的身体部分的表示的选定位置在每个选定位置的分析时间段期间根据与对询问 信号的响应相关的至少一个统计参数的值而分割。
[0011] 另一方面提供一种相应的计算机程序。
[0012] 另一方面提供一种相应的计算机程序产品。
[0013] 另一方面提供一种相应的系统。
[0014] 另一方面提供一种相应的诊断成像方法。
[0015] 更具体地说，本公开的一个或多个方面在独立权利要求中阐述，其有利的特征在 从属权利要求中阐述，所有权利要求的措辞在本文中通过引用逐字地并入（参照任何特定 方面提供的任何有利的特征适用于每一个其它方面）。
【附图说明】
[0016] 根据本公开的一个或多个实施例的解决方案及其进一步的特征和优点将参照W 下详细描述被最好地理解，W下详细描述仅仅作为非限制性指示而给出，应结合附图进行 阅读（其中，为了简单起见，相应的元件用相同的或相似的标号表示，并且它们的说明不被 重复，并且每个实体的名称一般用于表示其类型和其属性两者一-诸如值、内容和表示）。 具体地说：
[0017] 图1示出可W用于实施根据本公开的实施例的解决方案的系统的图形表示，
[001引图2示出根据本公开的实施例的解决方案的应用的体外情况的例子，
[0019] 图3A-图3B示出根据本公开的实施例的解决方案的应用的例子，
[0020] 图4A-图4D示出根据本公开的另一实施例的解决方案的应用的例子，
[0021] 图5A-图5C不出根据本公开的另一实施例的解决方案的应用的例子，
[0022] 图6示出表示可W用于实现根据本公开的实施例的解决方案的主要组件的作用 的协作图，和
[0023] 图7、图8和图9示出根据本公开的实施例的解决方案的体内应用的不同例子。
【具体实施方式】
[0024] 具体参照图1，示出了可W用于实施根据本公开的实施例的解决方案的系统100。 具体地说，系统100是超声扫描仪，其包括中央单元105W及连接到中央单元105的、例如 阵列类型的手持发射-接收成像探头（或简称为成像探头）11〇。成像探头110发射包括超 声脉冲序列的超声波（例如，具有1和50MHz之间的中屯、频率），并且接收由该些超声脉冲 在选定的扫描平面中的反射而导致的射频（RF)回声信号；为了该个目的，成像探头110设 有发射/接收复用器，其允许在上述脉冲-回声模式下使用成像探头110。
[00巧]中央单元105容纳主板115,在主板115上，安装控制超声扫描仪100的操作的电 子电路（例如，微处理器、工作存储器和硬盘驱动器）。而且，一个或多个子板（整个用标号 120表示）被插入到主板115中；子板120提供用于驱动成像探头110并且用于处理回声信 号的电子电路。中央单元105还配备有用于可移动盘130 (诸如CD或DVD)的驱动器125。 用于显示与进行中的分析处理相关的图像的监视器135连接到中央单元105。超声扫描仪 100的操作通过键盘140 (W常规方式连接到中央单元105)控制；优选地，键盘140设有用 于在监视器135上操纵指针（该图中未示出）的位置的轨迹球145。
[0026] 超声扫描仪100用于分析患者155的身体部分150;为了该个目的，（超声）造影 剂被施用于患者155。
[0027] 造影剂包括充当良好的超声反射体的颗粒。例如，造影剂是充气气泡在液态载体 中的悬浮液；通常，充气气泡具有大约0. 1-5ym的直径，W便使得它们可W保留在患者155 的血管系统内，但是同时使得它们可W通过他/她的毛细血管。充气气泡通常通过将气体 或者其前体夹带或封装到各种系统中来使其稳定，所述各种系统包括磯脂、乳化剂、油、增 稠剂、糖、蛋白质或聚合物；稳定的充气气泡一般被称为微泡。具体地说，分散在水介质中 并且在气体/液体界面处被包含表面活性剂（即，两亲性材料）的非常薄的包膜束缚的气 泡也被称为微气泡。可替代地，被由脂类或（天然或合成）聚合物形成的固态材料包膜包 围的微泡也被称为微气球或微胶囊。另一种类的造影剂包括聚合物或其它固体的多孔微粒 的悬浮液，该悬浮液携带包埋在微粒的孔内的气体的气泡。EP-A-0458745、W0-A-91/15244、 EP-A-0554213、W0-A-94/09829和W0-A-95/16467(其全部公开内容通过引用并入本文）中 描述了微泡（具体地说，微气泡和微气球）的合适的水性悬浮液及其制备的例子。商业的 包含微泡的造影剂的例子是化accoInternationalBV(商标）的Sono化e。
[0028] 造影剂还可W是目标特定的类型。目标特定的造影剂基本上自由地在患者155内 循环；然而，目标特定的造影剂还能够固定在选定的（生物）目标上，W便保留在基本固定 的位置上一定时间段。为了该个目的，目标特定的造影剂W通过与期望目标的特定相互作 用选择性地绑定到该期望目标的该样的方式配制。例如，该行为可W通过合并能够选择性 地绑定到期望组织或受体（例如，通过生物化学亲合性和/或静电作用）的目标特定的配 体来实现。目标特定的配体（其可W插入到微气泡的膜中）的例子是单克隆抗体、肤或多 糖。术语组织包括（在其意义内）单个的细胞W及细胞的聚合（诸如细胞膜或器官）。该 术语指的是正常的（健康的）或异常的（病态的）细胞或细胞聚合。组织的例子是屯、肌 组织（包括屯、室肌细胞和屯、肌细胞）、膜组织（诸如内皮和上皮）W及结缔组织；病态组 织的例子是梗塞屯、脏组织、血凝块、动脉粥样硬化斑块、炎性组织和肿瘤组织。受体包括位 于组织上（例如，在细胞内或者在它们的表面上）的、能够选择性地绑定到特定物质的任 何分子结构。示例性受体是糖蛋白GPnbllla或纤维蛋白（例如，位于血凝块或血栓上）、 P选择素（例如，位于发炎组织的激活内皮上）或KDR(例如，位于肿瘤组织）上。W下 文献中描述了合适的目标特定的造影剂和目标特定的配体的例子；"化chon等人，BR55:A lipopeptide-basedVEGFR2-targetedultrasoundcontrastagentformolecular imagingofangiogenesis.InvesRadiol2010;45:89-95"和"Tardy等人，叫trasound molecularimagingofVEGFR2inaratprostatetumormodelusingBR55,Inves Radiol2010 ;45:573-578" ;W及WO-A-2006018433 (该些文献的全部公开内容都通过引用 并入本文）。
[0029] 可W通过应用超声波并且记录响应于该些超声波返回的回声信号来检测造影剂； 因为造影剂在患者155中W与红血球相同的速率流动，所W其在身体部分150中的检测和 跟踪提供关于相应的血液灌注的信息。
[0030] 例如，造影剂被W大剂量一一即，在（大约2-20秒的）短时间段内用手通过注射 器提供的一次剂量一一静脉注射到患者 155。造影剂在患者155的血管系统内循环，W便灌 注身体部分150。同时，成像探头110被放置为在身体部分150的区域中与患者155的皮 肤接触。超声脉冲的成像帖序列施加于身体部分150 ;在分析（时间）段期间随着时间的 过去在连续的采集时刻（例如，在2-25秒期间W4-20HZ的帖率）响应于不同成像帖的超 声脉冲记录的回声信号提供身体部分150在该分析时间段期间的相应位置的表示（在其在 选定的扫描平面中的切片中）。回声信号由造影剂（如果存在）所产生的不同贡献和周围 组织的叠加得到。优选地，超声扫描仪100W对比度特定的成像模式操作，W便相对于造影 剂的（非线性）贡献大幅减小回声信号中的组织的主导（线性）贡献一一例如，基于谐波 成像化I)、脉冲反向（PI)、功率调制（PM)、或如"Rafter等人，Imagingtechnologiesand techniques,CardiologyClinics22 (2004),pp. 181-197"(该文献的全部公开内容通过引 用并入本文）中所述的其它对比脉冲序列（CP巧技术。然后产生相应的（对比度特定的） 视频图像的序列；对于每个采集时刻，视频图像包括关于表示身体部分的相应位置的每个 可视化元素（即，像素）的数字值；像素值由在该采集时刻关于该位置记录的回声信号的强 度值定义。该样，视频图像序列表示身体部分150在分析时间段期间的灌注的演变。
[0031] 图2中示出了关于根据本公开的实施例的解决方案的应用的体外情况的例子。
[0032] 具体地说，使（磯脂稳定的）微气泡WImm/s的平均流速循环通过仿组织流动体 模内所包含的直径为3mm的管。用SiemensMedicalSystems(商标）的Sequoia超声扫 描仪在(PS模式下在25秒内W4Hz的帖率（即，100个成像帖）采集相应的视频图像序列。 该图中分别用标号205a、20化和205c指示与成像帖No. 25、No. 43和No. 66相应的视频图 像。可W看出，管的上壁处的微气泡因为层流条件，基本上是静止的一一例如，对于用标号 210s表示的区域可见；相反，管的中屯、（与其上壁间隔开）处的微气泡W平均流速移动一一 例如，对于用标号210m表示的区域可见。
[0033] 在根据本公开的实施例的解决方案中，如W下详细描述的，分割（在该种情况下， 在静止微气泡的像素和移动微气泡的像素之间）基于像素值随着时间的统计分析。
[0034] 具体地说，图3A-图3B中示出了根据本公开的实施例的解决方案的应用的例子。
[003引从图3A开始，在纵轴上绘制像素值（任意单元，a.u.)对比在横轴上绘制帖数 (adimensional)的常见示图中，关于静止微气泡的区域和移动微气泡的区域中的代表性像 素的（线性化）像素值分别用曲线305s(实线）和曲线305m(点线）表示。
[0036] 可W看出，关于静止微气泡的像素（曲线305s)和移动微气泡的像素（曲线305m) 的像素值随着时间的动态大不相同；具体地说，关于静止微气泡的曲线305s几乎是恒定 的，而关于移动微气泡的曲线305m表现出高变化（由于微气泡的出现和消失）。因此，每个 像素的像素值随着时间的统计分布的统计参数可W用于区分它们。
[0037]例如，对于每个像素，计算其统计分布的（算术）均值、中间值和标准差。
[0038] 具体地说，具有离散值Xi(其中i=1…脚的统计变量的均值y(即，在该种情况 下是随时间的像素值）等于：
[0039]
[0040] 均值X是表示当统计变量具有对称形状的统计分布时该统计变量的期望值的集 中趋势统计参数（即，指示统计变量值围绕其聚集的值）。
[0041] 同一统计变量Xi的中间值y(从其值的最小一个Xi到最大一个Xw排序）等于：
[0044] 中间值y是表示当统计变量具有歪斜（即，不对称）形状的统计分布时该统计变 量的期望值的另一个集中趋势统计参数。
[0045] 同一统计变量Xi的标准差0等于：
[0046]
[0047] 标准差0是表示统计变量值如何从它们的期望值离散的离差统计参数（即，指示 统计变量值的变化性或伸展性）。
[004引现在参照图3B，该些统计参数可W用于分割图2的视频图像，W便鉴别静止微气 泡和移动微气泡。为了该个目的，计算=个不同的参数化图像305y、305d和305 0。具体 地说，参数化图像305y、305d和305 0具有与视频图像相同的大小；参数化图像305U、 305d和305 0的每个像素具有参数值，该参数值分别等于相应像素值随着时间的均值、中 间值和标准差的值。
[0049] 可W看出，静止微气泡的像素和移动微气泡的像素之间的对比度（或差异）（例 如，对于分别用标号310s和310m表示的区域可见）相当大地改进。具体地说，在关于均值 的参数化图像305U中W及在关于中间值的参数化图像305d中，参数值在静止微气泡（区 域310s)的像素中较高，而在移动微气泡（区域310m)的像素中较低；该是由于如下的事 实而导致的，即，像素值的期望值（用均值和中间值表达）在静止微气泡的像素中高于移动 微气泡的像素中。对比度在关于中间值的参数化图像305d中高于关于均值的参数化图像 305y中，因为在移动微气泡的像素中，它们的像素值的离群值（将均值拉向其值）不影响 中间值。相反，在关于标准差的参数化图像305 0中，参数值在静止微气泡（区域310s)的 像素中较低，而在移动微气泡（区域310m)的像素中较高；该是由于如下的事实而导致的， 即，像素值在静止微气泡的像素中随着时间的离差小于移动微气泡的像素中。例如，可W通 过计算相应像素的参数值的平均值之间的比率来测量静止微气泡和移动微气泡之间的对 比度。具体地说，该比率在关于均值的参数化图像305y中等于1. 69,在关于中间值的参数 化图像305d中等于2. 72,在关于标准差的参数化图像305 0中等于0. 84 ;该确认提供最好 结果的统计参数是中间值（之后是均值和标准差）。
[0050] W上提及的参数值（尤其是中间值）可W用于用阔值化技术分割原始图像。例 如，用于参数值的分割阔值可W被设置为等于静止微气泡的像素的参数值的平均值和移动 微气泡的像素的参数值的平均值之间的均值（手动选择）；然后，当每个像素的参数值高于 分割阔值时，将该像素分配给静止微气泡的部分，或者否则，将该像素分配给移动微气泡的 部分。
[0051]一般地，上述解决方案方便分割，有益效果是所获得的结果的准确度。
[0052]具体地说，该解决方案可W用于改进DCE-US成像应用中的量化表征的性能。例 如，该使得可W排除（或至少大幅减小）将被分析的包括病变的分析区域中的低灌注部分 (例如，具有坏死组织）的贡献和/或包括将被与病变进行比较的健康组织的控制区域中的 高灌注部分（例如，具有大血管）的贡献。结果，在分析区域的相对合并值中可能反映的相 应误差大幅减小（使得即使当分析/控制区域由操作者手动选择时，分割的可再现性也更 高，因为其质量现在不太取决于个人技能和经验）。
[0053]而且，相同的解决方案可W在UMI应用中用于改进目标特定的造影剂的固定颗粒 的检测和量化。例如，该使得可W移除（或至少大幅减小）目标特定的造影剂的循环颗粒 的贡献，尤其是在其施用之后的早期。
[0054]W上提及的优点在治疗随访的上下文中尤其感受得到，因为它们方便识别病变在 其处理期间的灌注中的即使是微小的变化，该些微小的变化可能指示处理的有效性。
[0055] 作为进一步的改进，分割可W基于与每个像素随着时间的像素值相关的两个或更 多个统计参数。
[0056]具体地说，图4A-图4D中示出了根据本公开的另一实施例的解决方案的应用的例 子。
[0057] 从图4A开始，在横轴上绘制均值并且在纵轴上绘制标准差值（adimensional)的 常见的散布图中，图3B的静止微气泡的区域（310s)和移动微气泡的区域（310m)中的像素 的均值和标准差值分别用点405AS(圆形）和点405Am(十字形）表示；具体地说，用标号 405AS(黑色实屯、圆）和405Am(黑色实屯、星）表示的点分别对应于图2的静止微气泡的区 域（210s)中的代表性像素和移动微气泡的区域（210m)中的代表性像素。
[0058] 计算用于均值的分割阔值THy= 97和用于标准差值的分割阔值THO= 57 (例 如，通过求取点405AS和405Am的相应参数值的平均值）。分割阔值THy和THo在散布图 中定义四个分割象限。具体地说，右上分割象限QAur由分割标准（均值〉THy和标准差值 〉TH0 )定义，左上分割象限QAul由分割标准（均值《THy和标准差值〉TH0 )定义，左下 分割象限QA11由分割标准（均值《THy和标准差值《TH0 )定义，右下分割象限QAlr由 分割标准（均值〉THy和标准差值《THo)定义。
[0059]同样地，在图4B中，在横轴上绘制均值并且在纵轴上绘制标准差值 (adimensional)的另一常见的散布图中，图3B的静止微气泡的区域（310s)和移动微气泡 的区域（310m)中的像素的均值和标准差值分别用点405BS(圆形）和点405Bm(十字形） 表不；具体地说，用t不号405Bs(黑色实屯、圆）和405Bm(黑色实屯、星）表不的点分别对应于 图2的静止微气泡的区域（210s)中的代表性像素和移动微气泡的区域（210m)中的代表性 像素。
[0060] 除了相同的用于标准差值的分割阔值THO= 57之外，还如上那样计算用于中间 值的分割阔值THd= 79。分割阔值THd和THO在散布图中定义其它四个分割象限。具体 地说，右上分割象限地ur由分割标准（均值〉THd和标准差值〉TH0 )定义，左上分割象限 地ul由分割标准（均值《THd和标准差值〉THO)定义，左下分割象限地11由分割标准（均 值《THd和标准差值《TH0 )定义，右下分割象限地Ir由分割标准（均值〉THd和标准差 值《THo)定义。
[0061] 通过比较图4A-图4B的散布图可W看出，静止微气泡的像素和移动微气泡的像素 之间的对比度在关于中间值/标准差的图4B的散布图中高于关于均值/标准差的图4A的 散布图中（例如，对于图4B中的点405Bm、405Bs相对于图4A中的点405Am、405As显而易 见）；该确认与标准差组合提供更好结果的统计参数再次是中间值。
[0062] W上参数值可W用于根据相应的分割标准来分割原始图像。
[0063] 例如，如图4C中所示，四个分割掩模MAur、MAul、Mll和Mir使用关于均值/标准 差的图4A的分别与分割象限QAur、QAul、QA11和QAlr相应的分割标准而构建。每个分割 掩模MAur、MAul、Mll和Mir具有与视频图像相同的大小；具有满足相应的分割标准的参 数值的像素是白色，而其它像素是黑色。而且，在所有的分割掩模MAur、MAul、MAll和MAlr 中，区域410表示用于采集视频图像的流动体模的3mm管（由表示其壁的两个虚线定界）。
[0064] 可W看出，分割掩模MAur和MAlr中的白色像素检测管410的上壁处的静止微气 泡的像素的一部分（即使管410的中屯、处的一些白色像素存在于分割掩模MAur中）；另 一方面，分割掩模MAul和Mil中的白色像素 检测管410的中屯、处的移动微气泡的像素的 一部分（分割掩模Mil中的管410外部的白色像素可能是由于噪声和体模组织信号导致 的）。
[006引 同样地，在图4D中，其它四个分割掩模MBur、MBul、MBll和MBlr如上使用关于中 间值/标准差的图4B的分别与分割象限地ur、地ul、地11和地Ir相应的分割标准而构建 (管的同一区域再次用标号410标识）。如上，分割掩模MBur和MBlr中的白色像素检测管 410的上壁处的静止微气泡的像素的一部分（即使管410的中屯、处的一些白色像素存在于 (但是程度较低）分割掩模MAur中）；另一方面，分割掩模MBul和MB11中的白色像素检测 管410的中屯、处的移动微气泡的像素的一部分（分割掩模MB11中的管410外部的白色像 素再次可能是由于噪声和体模组织信号导致的）。
[0066] 通过比较图4C-图4D的分割掩模可W看出，静止微气泡的像素和移动微气泡的像 素之间的分割在关于中间值/标准差的图4D的分割掩模中好于关于均值/标准差的图4C 的分割掩模中--例如，对于图4D中的分割掩模MBur相对于图4C中的分割掩模MAur是 显而易见的；该再次确认与标准差组合提供更好结果的统计参数是中间值。
[0067] 作为更进一步的改进，分割可W基于W上提及的与每个像素随着时间的像素值相 关的统计参数，该些统计参数被变换W提供具有歪斜形状的、（变换的）像素值的（变换 的）统计分布。
[006引具体地说，图5A-图5C中示出了根据本公开的另一实施例的解决方案的应用的例 子。
[0069] 从图5A开始，假定（变换的）统计分布是对数正态统计分布。在对数正态统计分 布中，统计变量的自然对数具有正态统计分布。具体地说，在连续的统计变量X的情况下， 表示（连续的）正态统计分布的统计密度函数是正态统计密度函数（或简称为正态函数）：
[0070]
[0071] 其中，相同的符号义和0分别用于指示统计变量X的均值和标准差；因此，表示 (连续的）对数正态统计分布的统计密度函数是对数正态统计密度函数（或简称为对数正 态函数）；
[0072]
[0073] 其中，m和S现在分别是统计变量X的自然对数（即，InOO)的均值和标准差。
[0074] 在该种情况下，通过应用W下变换公式，可W从（正态）均值支和（正态）中间值 X直接计算（对数正态）均值m和（对数正态）标准差S:
[0078] 分割现在基于由对数正态均值和对数正态标准差定义的参数值。为了该个目的， 通过将W上提及的变换公式应用于每个像素的随着时间的像素值的实际的均值和中间值 来计算该像素的对数正态均值和对数正态标准差值。在横轴上绘制对数正态均值并且在 纵轴上绘制对数正态标准差值的常见的散布图（adimensional)中，静止微气泡的区域 (310s)和移动微气泡的区域（310m)中的像素的对数正态均值和对数正态标准差值分别 用点505s(圆形）和点505m(十字形）表示；具体地说，用标号505s'(黑色实屯、圆）和 505m'(黑色实屯、星）表示的点分别对应于图2的静止微气泡的区域（210s)中的代表性像 素和移动微气泡的区域（210m)中的代表性像素。
[0079] 如上计算用于对数正态均值的分割阔值T血=4. 4和用于对数正态标准差值的分 割阔值THs= 0. 54。分割阔值T血和ras在散布图中定义四个分割象限。具体地说，右上 分割象限QLur由分割标准（对数正态均值〉T血和对数正态标准差值〉THs)定义，左上分 割象限QLul由分割标准（对数正态均值《T血和对数正态标准差值〉THs)定义，左下象限 QL11由分割标准（对数正态均值《T血和对数正态标准差值《THs)定义，右下象限化Ir 由分割标准（对数正态均值〉T血和对数正态标准差值《THs)定义。
[0080] 可W看出，静止微气泡的像素和移动微气泡的像素之间的对比度进一步改进；实 际上，在该种情况下，静止微气泡和移动微气泡的点聚集在散布图的不同区域中，重叠减 小。该进一步方便分割，因为分割阔值T血、ras现在可W独立地定义，而不会损害分割的性 能。
[0081] 移至图5B，四个分割掩模MLur、MLul、ML11和MLlr如上使用图5A的分别与分割 象限化ur、QLul、化11和化Ir相应的分割标准而构建（管的同一区域再次用标号410标 识）。
[0082] 在该种情况下，分割掩模MLlr中的白色像素检测管410的上壁处的静止微气泡的 大多数像素；另一方面，分割掩模MLul中的白色像素检测管410的中屯、处的移动微气泡的 大多数像素（掩模MLul中的管410外部的白色像素可能是由于CPS的声学伪像与管410 中的微气泡的高浓度组合而导致的，所述声波伪像可W通过降低微气泡的浓度来减小）。相 反，分割掩模MLur和ML11的大多数像素是黑色的（分割掩模MB11中的管410外部的白色 像素再次可能是由于噪声和体模组织信号而导致的）。该确认像素值的对数变换改进了分 割性能。
[0083] 如图5C中所示，W上定义的分割掩模可W应用于图3B的参数化图像305y(表示 每个像素的均值）。具体地说，四个分割图像SIur、SIul、SIll和Sllr通过分别将图5B的 分割掩模MLur、MLul、ML11和MLlr应用于该参数化图像上而构建；分割图像Slur、SIul、 Sill和Sllr的每个像素被分配其参数值，即，如果它在相应的分割掩模中为白色，则被分 配均值（为了显示的目的，用压缩因子40地进行对数压缩），否则，它保持黑色（W便仅显 示分割像素的均值）。结果，分割图像Sllr主要示出管410的上壁处的静止微气泡的像素， 而分割图像SIul主要示出管410的中屯、处的移动微气泡的像素（W便方便区分它们）。
[0084] 图6中示出了表示可W用于实现根据本公开的实施例的解决方案的主要组件的 作用的协作图。
[0085] 具体地说，所有的组件作为整体用标号600表示；该图描述了组件600的静态结构 及其动态行为（通过一系列交换的消息，每个交换的消息表示相应的动作，该动作用前面 为符号"A"的序号表示）。
[0086] 信息（程序和数据）通常被存储在硬盘中，并且当程序连同操作系统和其它应用 程序（该图中未示出）一起运行时，被（至少部分）加载到超声扫描仪的工作存储器中。程 序最初例如从移动盘安装到硬盘上。在该方面，组件600中的每个可W表示包括用于实现 指定的一个（或多个）逻辑函数的一个或多个可执行指令的模块、代码段或代码部分。
[0087] 更具体地说，采集器603用于采集视频图像。为了该个目的，采集器603包括控制 成像探头的驱动器。例如，成像探头驱动器设有发射波束形成器和脉冲器，它们用于在每个 采集时刻产生将施加于所分析的身体部分的超声脉冲的成像帖；超声脉冲具有低声能（诸 如具有机械指数MI=0.01-0. 3)，W便引入造影剂的可忽略的破坏（诸如小于其在连续的 成像帖之间的局部浓度的5%，优选地，小于1 % )。成像探头然后在每个采集时刻接收身体 部分的不同位置（在选定的扫描平面中）所反射的（模拟RF)回声信号。该些模拟RF回 声信号被供给接收处理器，该接收处理器预先放大模拟RF回声信号，并且施加初步时域补 偿（TGC);模拟RF回声信号然后被模数转换器（ADC)转换为数字值，并且被接收波束形成 器组合为聚焦的波束信号。如此获得的数字RF回声信号优选地通过进一步的数字算法和 其它的线性或非线性信号调理器（例如，后波束形成TGC)进行处理；具体地说，接收处理器 应用对比度特定的算法（诸如基于(PS技术）来减小组织的贡献。数字RF回声信号然后 被解调、对数压缩（为了获得具有很好平衡的对比度的图像）、扫描转换为视频格式并且根 据给定的色表调色盘擅染。该样，产生每个采集时刻的视频图像。该视频图像由单元矩阵 定义，每个单元存储相应的像素值。像素值定义像素的亮度；例如，在灰阶视频图像中，像素 值可W在8个比特上编码，作为相应的回声信号的强度的函数，从0 (黑色）增大到255 (白 色）（表示身体部分的相应位置在该采集时刻的声学响应）。
[0088] 在分析处理开始时，超声扫描仪的操作者启动成像探头，并且围绕将被分析的身 体部分移动它（在施用任何造影剂之前）。相应的视频图像一被采集，就被提供给显示器 606 ;显示器606包括驱动器，该驱动器控制超声扫描仪的监视器，W便使得实时地显示该 些视频图像（动作"A1.显示"）。
[0089] 操作者选择表示身体部分的特定切片的扫描平面（例如，包括将被分析的病变）， 并且将成像探头保持固定在该位置上一段预定义的分析时间段。造影剂现在被施用于患 者，并且相应的视频图像序列（表示身体部分在分析时间段期间在所选扫描平面中的灌 注）被如上那样采集。该些视频图像（除了被显示器606实时显示之外，在该图中为了清 楚起见未示出显示器606的连接）还通过将它们保存到贬存器609而被登记（动作"A2. 登记"）。
[0090] 选择器612被操作者用于在贬存器609的一个任意选择的视频图像中选择身体部 分的将被分析的期望的感兴趣区域（R0I)。感兴趣区域用R0I掩模表示，R0I掩模包括具有 与视频图像相同大小的单元的矩阵；R0I掩模的每个单元存储R0I标志（即，二进制值），当 相应像素在感兴趣区域的内部（即，它已经被选择）时，该R0I标志被断言（例如，为逻辑 值1)，否则，该R0I标志被解除断言（例如，为逻辑值0)。R0I掩模被保存到贬存器615中 (动作"A3.选择"）。
[0091] 预处理器618访问贬存器609中的视频图像和贬存器615中的R0I掩模。对于每 个视频图像，预处理器618对其在R0I掩模中的R0I标志被断言（即，在感兴趣区域的内部） 的每个单元的像素值进行处理，W便使得它与造影剂的相应的局部浓度成正比。例如，该结 果可W通过下述方式实现，即，应用逆序调色板-照度表和逆序对数压缩（W使采集器603 应用它们的效果逆转），然后求取如此获得的值的平方值（如W0-A-2004/110279中所述那 样，该申请的全部公开内容通过引用并入本文）。预处理器618然后创建像素值矢量图，该 像素值矢量图包括具有与视频图像相同大小的单元的矩阵；像素值矢量图的其在R0I掩模 中的R0I标志被断言的每个单元存储像素值矢量，该像素值矢量包括沿着视频图像的相应 的线性化像素值的序列（其长度对应于分析时间段的持续时间）。该像素值矢量图被保存 到贬存器621中（动作"A4.预处理"）。
[0092] 贬存器621中的像素值矢量图被供给计算器624,计算器624也访问贬存器 615(在该图中为了清楚起见未示出其连接）中的R0I掩模。对于像素值矢量图的其在R0I 掩模中的R0I标志被断言（即，在感兴趣区域的内部）的每个单元，计算器624计算相应的 (线性化的）像素值的序列的均值和中间值的值。计算器624然后创建关于均值的参数化 图像和关于中间值的参数化图像，每个图像包括具有与视频图像相同大小的单元的矩阵； 关于均值的参数化图像和关于中间值的参数化图像的其在ROI掩模中的ROI标志被断言的 每个单元分别存储相应的均值和中间值。该些参数化图像被保存到贬存器627中（动作 "A5.计算"）。
[0093] 贬存器627中的参数化图像被供给变换器630,变换器630也访问贬存器615(在 该图中为了清楚起见未示出其连接） 中的R0I掩模。对于参数化图像的其在R0I掩模中的 R0I标志被断言（即，在感兴趣区域的内部）的每个像素，变换器630通过应用相应的变换 公式从相应的均值和中间值计算对数正态均值和对数正态标准差的值。变换器630然后创 建统计图，该统计图包括与视频图像相同大小的单元的矩阵；该统计图的其在R0I掩模中 的R0I标志被断言的每个单元存储相应的统计参数（即，对数正态均值和对数正态标准差 值）。该统计图被保存到贬存器633中（动作"A6.变换"）。
[0094] 贬存器633中的统计图被供给分割器636,分割器636也访问贬存器615(在该图 中为了清楚起见未示出其连接）中的R0I掩模。分割器636从贬存器639提取选定的一个 (或多个）分割标准，贬存器639存储若干个基于用于对数正态均值和对数正态标准差值的 分割阔值的预定分割标准；例如，该些分割阔值可能已经在学习阶段期间通过分析具有公 知特性的样品身体部分而被设置。对于统计图的其在R0I掩模中的R0I标志被断言（即， 在感兴趣区域的内部）的每个像素，分割器636验证相应的统计参数是否满足选定的分割 标准（意指它是否属于表示身体部分的具有相应的同质特性的部分的相应部分）。分割器 636然后创建用于选定的分割标准的分割掩模，该分割掩模包括具有与视频图像相同大小 的单元的矩阵；分割掩模的每个单元存储分割标志（即，二进制值），当满足相应的分割标 准时，该分割标志被断言（目P，相应的像素属于该部分），否则，该分割标志被解除断言。如 此创建的每个分割掩模被保存到贬存器642中（动作"A7.分割"）。
[0095] 施加器645从贬存器642提取选定的一个（或多个）分割掩模；施加器645还从 贬存器627提取选定的一个（或多个）参数化图像，和/或从贬存器609提取选定的一个 (或多个）视频图像。施加器645通过将它们逐个单元地相乘来将（选定的）分割掩模施 加于（选定的）参数化/视频图像。该操作创建分割图像，该分割图像包括与视频/参数 化图像相同大小的单元的矩阵；分割图像的属于由分割掩模限定的部分（目P，其分割标志 被断言）的每个单元存储相应的参数化/像素值，而其它单元则被重置为0。如此创建的每 个分割图像被保存到贬存器648中（动作"A8.施加"）。
[0096] 贬存器648中的每个分割图像现在可W被传递给显示器606W供其显示W及任何 下一次使用（动作"A9.显示")。例如，区分目标特定的造影剂的固定颗粒与其余的循环颗 粒的分割图像可W用于检测固定颗粒，因此识别身体部分中的相应病变。另外地或者可替 代地，合并器651可W从贬存器642提取选定的一个（或多个）分割标准，和/或从贬存器 648提取选定的一个（或多个）分割图像。合并器651计算身体部分的由（选定的）分割 掩模/分割图像限定的部分所表示的部分的一个（或多个）合并值；例如，合并值可W被设 置为等于在分割掩模中断言的分割标志的数量（测量身体部分的所述部分），或者合并值 可W被设置为等于分割图像的像素/参数值的和或平均值（表征身体部分的所述部分）。 如此计算的每个合并值被保存到贬存器654中W供下一次使用（动作"A10.合并"）。例 如，其中低灌注部分被负向区分的分割图像中的与包括病变的分析区域相关的R0I的参数 值的平均值自动地不受它们的分布的影响，并且其中高灌注部分被负向区分的分割图像中 的与包括将被与病变进行比较的健康组织的控制区域相关的ROI的参数值的平均值自动 地不受它们的分布的影响；因此，该些合并值可W组合为用相对于控制区域的合并值的相 对项表达的、分析区域的合并值（例如，等于它们的比率）。而且，在分割掩模中被断言的、 区分目标特定的造影剂的固定颗粒的分割标志的数量W及分割图像中的再次区分目标特 定的造影剂的固定颗粒的像素/参数值的和，该两者都量化该些固定颗粒。
[0097] 图7、图8和图9中示出了根据本公开的实施例的解决方案的体内应用的不同例 子。
[0098] 从图7开始，肝肿瘤的视频图像705是在注入造影剂之后（在峰值造影期间），在 对比度特定的成像模式下20秒采集的；两个（低灌注）区域710a和710b(较暗）在视频 图像705上被手动勾画。
[0099] 分割掩模715是W分割标准（对数正态均值乂. 0和对数正态标准差值〉〇)，从在 视频图像705的采集时刻左右的12. 5秒采集的肝肿瘤的100个视频图像的序列创建的。可 W看出，分割掩模715排除了两个区域720a和72化（黑色），如被排除的区域720a、72化和 在分割掩模715上用灰色虚线报告的低灌注区域710a、71化的轮廓之间的良好重叠所确认 的，区域720a和72化分别基本上对应于低灌注区域710a和710b。
[0100] 移至图8,将用作参考区域的肝脏的健康组织的视频图像805被如上那样采集；与 大血管相应的（高灌注）区域810(较亮）在视频图像805上被手动勾画。
[0101] 分割掩模815是W分割标准（对数正态均值巧.9和对数正态标准差值〉0)，从视 频图像805的采集时刻左右的3秒采集的肝脏的25个视频图像的序列创建的。可W看出， 分割掩模815排除了区域820 (黑色），如被排除区域820和在分割掩模815上用灰色虚线 报告的高灌注区域810的轮廓之间的良好重叠所确认的，区域820基本上对应于高灌注区 域 810。
[0102] 最后参照图9,具有（腺癌）肿瘤的模型鼠的前列腺的视频图像905是在注入针对 血管细胞粘附分子1 (VCAM-1)的目标特定的造影剂之后，在对比度特定的成像模式下2分 钟采集的；肿瘤区域910 (较亮）在视频图像905上被手动勾画。
[0103]分割掩模是W分割标准（对数正态均值M. 4和对数正态标准差值〉0. 54)，从在视 频图像905的采集时刻之后的15秒采集的、前列腺的60个视频图像的序列创建的；然后通 过将该分割掩模施加于从与上相同的视频图像序列创建的、关于均值的参数化图像来创建 分割图像915。可W看出，分割图像915基本上示出了仅肿瘤区域910中的目标特定的造 影剂的固定颗粒，目标特定的造影剂的其余的循环颗粒在视频图像905中存在于肿瘤区域 910的外部，肿瘤区域910在分割图像915中基本上不存在。而且，同样地在肿瘤区域910 中，目标特定的造影剂的其余的循环颗粒基本上被移除（通过将它与视频图像905中的肿 瘤区域910进行比较，显而易见）。
[0104] 作为控制，在相同的条件下，但是现在是使用（非针对性的）造影剂来采集另一视 频图像925,并且如上创建相应的分割图像935。
[0105] 在该种情况下，造影剂的（循环/固定）颗粒在分割图像935中基本上都是不可 见的（两者都在肿瘤区域910的内部和外部）。
[0106] 修改
[0107] 自然地，为了满足局部要求和特定要求，本领域技术人员可W将许多逻辑和/或 物理修改和改变应用于上述解决方案。更具体地说，尽管已经参照该解决方案的一个或多 个实施例一定程度地具体地描述了该解决方案，但是应理解，形式和细节上的各种省略、替 换和改变是可能的。具体地说，本公开的不同实施例甚至可W在没有在前面的描述中所阐 述的、提供其更透彻的理解的特定细节（诸如数值）的情况下实施；相反，公知特征可能已 经被省略或简化W便不使不必要的细节模糊本描述。而且，明确的意图是，与所公开的解决 方案的任何实施例结合描述的特定的元件和/或方法步骤可W作为一般设计选择合并在 任何其它实施例中。在任何情况下，顺序限定词或其它限定词仅用作区分具有相同名称的 元件的标签，而不独自暗示任何优先权、优先级或次序。而且，术语包括、具有、包含和设及 (及其任何形式）的意图应是开放式的、非穷举的意义（即，不限于所记载的项），术语基 于、取决于、根据、…的功能（及其任何形式）的意图应是作为非排他性的关系（即，设及 可能的另外的变量），术语一个的意图应是作为一个或多个项（除非另有明确指示）。
[010引例如，实施例提供一种用于诊断成像应用中的数据处理分割方法。该方法包括W下步骤。提供身体部分（被灌注造影剂）在非零分析时间段期间的表示；对于身体部分的 一组位置中的每个位置，该表示包括该位置在分析时间段期间对询问信号的响应的指示。 对于一组选定位置中的每个选定位置，计算该选定位置在分析时间段期间的响应的统计分 布的至少一个统计参数的值；所述一组选定位置包括所有的位置或者其一部分。根据关于 选定位置的所述至少一个统计参数的值与至少一个分割阔值之间的比较，分割选定位置。
[0109] 然而，身体部分的表示可任何方式（参见下面）提供，W便提供其在具有任 何非零长度（长得足W在统计上显著）的分析时间段期间的表示；该表示可W与任何类型 的灌注（例如，基于任何类型的造影剂（非针对性或针对性类型）的任何类型的灌注（例 如，基于通过破坏飞瓣的连续注射））相关。所述方法可W应用于任何类型的位置的级别上 (例如，像素、体素或者其群组）；具体地说，还可W对所述表示在空间上进行二次采样，确 定二次采样的表示上的分割，然后通过插值技术恢复分割的原尺寸W使其应用于原始表示 上。而且，所述方法可W应用于任何选定位置（例如，在W任何方式定义的一个或多个感兴 趣区域中）、或者甚至身体部分的整个表示。统计参数可W是任何类型和任何数量（参见下 面）；然而，可W根据统计参数的值与任何类型和数量的分割阔值之间的任何比较将选定 位置分割为任何数量的部分（两个或更多个）（参见下面）。
[0110] 在实施例中，所述提供身体部分在非零分析时间段期间的表示的步骤包括：提供 在所述分析时间段期间登记的身体部分的数字图像序列；每个数字图像包括关于每个位置 的强度值，该强度值基于该位置在分析时间段中的相应采集时刻的响应。所述计算至少一 个统计参数的值的步骤包括从沿着数字图像的选定位置的强度值计算所述至少一个统计 参数的值。
[0111] 然而，数字图像可WW任何方式（例如，通过对它们进行预处理W补偿运动伪像 或减去背景图像）登记。数字图像可W是任何类型（例如，单独地或组合地、3D类型、基于 任何对比度特定的成像模式、或者甚至基于基板B模式）、任何大小和任何非零数量（W便 在统计上显著）；强度值可W是任何类型（例如，对数压缩值）并且W任何方式编码（例如， 对任何数量的比特进行颜色编码）。在任何情况下，不排除将同一方法应用于身体部分的任 何表示（甚至不是图像形式）的可能性。例如，身体部分在非零分析时间段期间的表示可 W由原始（未处理RF)回声信号定义。在该种情况下，对于每个选定位置，（响应的统计分 布的）所述至少一个统计参数的值通过下述方式获得，即，计算选定位置的原始回声信号 的HiAed变换，然后在其幅值或相位的级别上计算（复数）相关系数或者其方差（然后 如上那样，根据该些统计参数的值和至少一个分割阔值之间的比较来分割选定位置）。
[0112] 在实施例中，所述至少一个统计参数是多个统计参数。
[0113] 然而，统计参数可W是任何数量（下至一个，或者相反，=个或更多个）。
[0114] 在实施例中，统计参数是集中趋势统计参数和离差统计参数。
[0115] 然而，可W使用任何集中趋势统计参数（例如，除了（算术）均值和中间值、模式、 几何均值、谐波均值、截断均值、中间范围均值等之外）和/或任何离差统计参数（例如，除 了标准差、均值差值、平均偏差、方差等之外）。在任何情况下，可W使用不同的、附加的或替 代的类型的统 计参数（例如，位置或形状统计参数）。
[0116] 在实施例中，所述计算至少一个统计参数的值的步骤包括计算与响应相关的所述 至少一个统计参数的值，该响应被变换W提供变换的响应在分析时间段期间的、具有歪斜 形状的变换的统计分布。
[0117] 然而，该计算可W根据任何变换的统计分布W任何方式实现（参见下面）。在任何 情况下，不排除仅从响应的（原始）统计分布计算统计参数的可能性。
[0118] 在实施例中，变换的统计分布是对数正态统计分布。
[0119] 然而，可W使用任何其它的具有歪斜形状的变换的统计分布（例如，伽马或局部 密度随机游走分布）。
[0120] 在实施例中，所述计算至少一个统计参数的值的步骤包括；从选定位置在分析时 间段期间的响应的统计分布计算至少一个另一统计参数的值，并且通过将根据变换的统计 分布的相应变换应用于所述至少一个另一统计参数的值来计算每个统计参数的值。
[0121] 然而，变换可W是任何类型的（参见下面）。在任何情况下，没有什么能阻止变换 响应、然后从变换的响应直接计算期望的统计参数；例如，可W如下计算对数正态均值和对 数正态标准差：
[0125] 如果Xi<l，则s= 0
[0126] 可替代地，还可W通过参数化函数（直接地或者从其直方图）拟合每个响应、然后 从该参数化函数计算期望的统计参数，该参数化函数是表示变换的统计分布的分布函数的 参数化模型的实例。
[0127] 在实施例中，所述通过应用相应变换计算每个统计参数的值的步骤包括将对数变 换应用于所述至少一个另一统计参数的值。
[012引然而，变换可W是任何类型的（例如，指数类型）。
[0129] 在实施例中，所述至少一个统计参数是对数正态均值和对数正态标准差，并且所 述至少一个另一统计参数是正态均值和正态中间值。
[0130] 然而，所述统计参数和/或另一统计参数可W是任何类型和任何数量的（例如，如 w上所指出的）。
[0131] 在实施例中，所述通过应用相应变换来计算每个统计参数的值的步骤包括：设置 对数正态均值的值等于正态中间值的值的自然对数，并且如果正态均值的值和正态中间值 的值之间的比率严格地高于1，则设置对数正态标准差的值等于该比率的自然对数的两倍 的平方根，否则，等于0。
[0132] 然而，可W根据变换的统计分布和/或统计参数的类型使用其它变换。
[0133] 在实施例中，所述分割选定位置的步骤包括根据每个统计参数的值和所述至少一 个分割阔值中的相应的一个分割阔值之间的比较来分割选定位置。
[0134] 然而，分割阔值可任何方式定义。例如，可W应用主成分分析（PCA)技术来解 除统计参数的值的相关性；而且，分割阔值可W从相应的统计参数的值的重屯、计算。可W根 据每个统计参数和相应的分割阔值（例如，就绝对值而言或者根据预定义百分比或者任何 组合，较高或较低）之间的任何比较来分割选定位置。在任何情况下，没有什么能阻止使用 与统计参数之一不同的若干个分割阔值（例如，当分割基于统计参数的任何（线性或非线 性）组合时）。
[0135] 在实施例中，所述分割选定位置的步骤包括；当集中趋势统计参数的值高于相应 的分割阔值并且离差统计参数的值高于相应的分割阔值时，将身体部分的分析区域（包括 将被分析的病变）的每个选定位置分配给第一部分（与分析区域的低灌注相关），否则，分 配给第二部分；所述方法还包括计算分析区域的不包括在第一部分中的选定位置的响应的 合并值。
[0136] 然而，该分割可W基于其它分割标准（例如，原始回声信号的化化6的变换的高相 关系数和低方差）；而且，它可W用于计算分析区域的任何合并值（例如，它们到达峰值的 平均时间、歪斜等）。
[0137] 在实施例中，所述分割选定位置的步骤包括；当集中趋势统计参数的值低于相应 的分割阔值并且离差统计参数的值高于相应的分割阔值时，将身体部分的控制区域（包括 将被与分析区域进行比较的健康组织）的每个选定位置分配给另一第一部分（与控制区域 的高灌注部分相关），否则，分配给另一第二部分；所述方法还包括计算控制区域的不包括 在所述另一第一部分中的选定位置的响应的另一合并值。
[013引然而，该分割可W基于其它分割标准（例如，原始回声信号的化化6的变换的低相 关系数和高方差）；然而，它可W用于如上那样计算控制区域的任何合并值。
[0139] 在实施例中，所述方法还包括根据所述合并值和另一合并值的组合来计算分析区 域相对于控制区域的相对合并值。
[0140] 然而，相对合并值可任何方式计算（例如，按照所述合并值和另一合并值之 间的差值计算）；在任何情况下，可任何其它方式使用所述另一合并值（例如，用于规 范化整个参数化图像的参数值）。
[0141] 在实施例中，造影剂是能够在所述患者内循环并且能够基本固定在生物目标上的 针对性的造影剂；所述分割选定位置的步骤包括；当集中趋势统计参数的值高于相应的分 割阔值并且离差统计参数的值低于相应的分割阔值时，将每个选定位置分配给第一部分 (与固定的针对性的造影剂相关），否则，分配给第二部分。
[0142] 然而，针对性的造影剂可W是任何类型（例如，基于特定的或非特定的相互作 用）；而且，该分割可W基于其它分割标准（例如，原始回声信号的HUbed变换的高相关 系数和低方差）。
[0143] 在任何情况下，相同的技术可W用于任何其它诊断成像应用中；例如，相同的过程 可W对选定的部分重复迭代（一次或多次），或者它可W重复地W不同的分割标准应用，其 中，所得的分割掩模通过逻辑运算符组合，诸如W获得相应部分的交集（逻辑AND)、并集 (逻辑OR)或异并狂OR)。
[0144] 一般地，如果相同的解决方案用等同的方法（通过使用其多个步骤或部分的具有 相同功能的类似步骤、移除不必要的一些步骤、或者添加另外的操作步骤）实现，则类似的 考虑因素适用；而且，该些步骤可W按不同次序、同时地或者（至少部分）W交错的方式执 行。
[0145] 在任何情况下，强调的是，上述方法是可W独立于与患者（W及特别是与可W在 执行该方法之前预先施用于患者的造影剂）的任何相互作用实现的数据处理（或计算）方 法。而且，造影剂可非侵入性的方式（例如，口服W用于使胃肠道成像或者经由雾化 器进入气道中）施用于患者，或者在对其没有任何实质的物理干预的任何情况下施用于患 者，所述物理干预将需要专业的医疗经验或者承担患者的健康风险（例如，肌肉注射）。在 任何情况下，尽管所提出的方法可W方便医生的任务，但是它一般仅提供中间结果，该些中 间结果可W帮助他/她分析身体部分一一例如，出于诊断的目的（即使用于治愈目的的诊 断在严格意义上总是由医生自己做出）。
[0146]另一实施例提供一种计算机程序，该计算机程序被配置为当该计算机程序在计算 系统上执行时使该计算系统执行W上提及的方法。
[0147] 另一实施例提供一种计算机程序产品，该计算机程序产品包括包含计算机程序的 非暂时性计算机可读介质；该计算机程序可加载到计算系统的工作存储器中，从而将该计 算系统配置为执行相同的方法。
[0148] 然而，相同的解决方案可W实现为独立的模块、用于超声扫描仪的控制程序的插 件、或者甚至直接在控制程序本身中实现；将显而易见的是，还可W将相同的解决方案部署 为通过网络访问（诸如在互联网中）的服务。在任何情况下，如果软件程序（其可W用于 实现本公开的每个实施例）W不同的方式构造，或者如果提供附加模块或功能，则类似的 考虑因素适用；同样地，存储器结构可W是其它类型，或者可W用等同的实体（不一定由物 理存储介质构成）取代。程序可W采取适合于被任何计算（或数据处理）系统使用的或者 与其结合使用的形式（例如，在虚拟机内），从而将系统配置为执行期望的操作；具体地说， 程序可W是外部或驻留软件、固件或微代码（用对象代码或源代码一一例如，将被编译或解 释）的形式。而且，可W在任何计算机可用介质上提供程序（具体地说，作为非暂时性介质 上的制造品）；介质可W是适合于包含、存储、传送、传播或传输程序的任何元件。例如，介 质可W是电子、磁、光、电磁、红外或半导体类型；该样的介质的例子是固定盘（其中程序可 W被预先加载）、可移动盘、带、卡、线、光纤、无线连接、网络、广播波等。在任何情况下，根据 本公开的实施例的解决方案适合于甚至用硬件结构实现（例如，集成在半导体材料的巧片 中），或者用被适当地编程的或者W其它方式配置的软件和硬件的组合实现。
[0149] 另一实施例提供一种包括被构造为执行W上提及的方法的步骤的装置的系统。
[0150] 然而，相同的解决方案可W应用于包括超声扫描仪和不同的计算机的系统（或任 何等同的系统）中；在该种情况下，所记录的信息从超声扫描仪传送到计算机w供对其进 行处理（例如，通过数字、模拟或网络连接）。在任何情况下，不排除将所提出的解决方案 应用于任何其它的诊断成像系统一-例如，基于磁共振成像（MRI)或X射线计算断层扫描 (CT)--的可能性。
[0151] 一般地，如果系统具有不同的结构或者包括等同的组件，或者它具有其它的操作 特性，类似的考虑因素适用。在任何情况下，其每一个组件可W划分为多个元件，或者两个 或更多个组件可W-起组合为单个元件；而且，每个组件可W被复制W支持相应操作的并 行执行。而且，除非另有指定，否则不同组件之间的任何相互作用一般不需要是连续的，并 且它可W是直接的或者通过一个或多个中间物间接的。
[0152] 另一实施例提供一种包括W下步骤的诊断成像方法。对患者施用造影剂W使造影 剂灌注患者的身体部分。在非零分析时间段期间将询问信号施加于身体部分。采集身体部 分在分析时间段期间的表示；对于身体部分的一组位置中的每个位置该表示包括该位置在 分析时间段期间对询问信号的响应的指示（该些表示根据与W上方法相同的方法进行处 理W分割一组选定位置，该一组选定位置包括所有的位置或者其一部分）。根据所述分割评 估身体部分的状况。
[0153] 然而，所述方法可W在任何种类的诊断应用中找到应用（从该术语的最广泛的意 义上来讲一一例如，旨在发现新的病变或监视已知的病变），并且用于分析任何（人类或动 物）患者的任何种类的身体部分（例如，器官（诸如肝脏、前列腺或屯、脏）、区域或组织）。 具体地说，分割可W用于识别身体部分的无论什么生理状况的区域一一例如，区分病变与 健康组织或者身体部分的具有不同的灌注参数的部分。
【主权项】
1. 一种用于诊断成像应用中的数据处理分割方法（Al-AlO)，所述方法包括： 提供（A1-A4)被灌注造影剂的身体部分在非零分析时间段期间的表示，对于所述身体 部分的一组位置中的每个位置所述表示包括所述位置在所述分析时间段期间对询问信号 的响应的指示； 对于一组选定位置中的每个选定位置，计算（A5-A6)所述选定位置在所述分析时间段 期间的响应的统计分布的至少一个统计参数的值，所述一组选定位置包括所有的所述位置 或者其一部分；和 根据关于所述选定位置的所述至少一个统计参数的值与至少一个分割阈值之间的比 较，分割（A7-A9)所述选定位置。2. 根据权利要求1所述的方法（Al-AlO)，其中，所述提供（A1-A4)身体部分在非零分 析时间段期间的表示包括： 提供（A2)在所述分析时间段期间登记的所述身体部分的 数字图像序列，每个数字图 像包括关于每个位置的强度值，所述强度值基于所述位置在所述分析时间段中的相应采集 时刻的响应； 并且其中，所述计算（A5-A6)至少一个统计参数的值包括： 从沿着所述数字图像的选定位置的强度值计算（A5)所述至少一个统计参数的值。3. 根据权利要求1或2所述的方法（Al-AlO)，其中，所述至少一个统计参数是多个统 计参数。4. 根据权利要求3所述的方法（Al-AlO)，其中，所述统计参数是集中趋势统计参数和 离差统计参数。5. 根据权利要求1至4中的任何一个所述的方法（Al-AlO)，其中，所述计算（A5-A6) 至少一个统计参数的值包括： 计算与所述响应相关的所述至少一个统计参数的值，所述响应被变换以提供所述变换 的响应在所述分析时间段期间的、具有歪斜形状的变换的统计分布。6. 根据权利要求5所述的方法（Al-AlO)，其中，所述变换的统计分布是对数正态统计 分布。7. 根据权利要求5或6所述的方法（A1-A10)，其中，所述计算（A5-A6)至少一个统计 参数的值包括： 从所述选定位置在所述分析时间段期间的响应的统计分布计算（A5)至少一个另一统 计参数的值；和 通过将根据所述变换的统计分布的相应变换应用于所述至少一个另一统计参数的值 来计算（A6)每个统计参数的值。8. 根据权利要求7所述的方法（A1-A10)，其中，所述通过应用相应变换来计算（A6)每 个统计参数的值包括： 将对数变换应用（A6)于所述至少一个另一统计参数的值。9. 根据权利要求7或8所述的方法（Al-AlO)，其中，所述至少一个统计参数是对数正 态均值和对数正态标准差，并且所述至少一个另一统计参数是正态均值和正态中间值。10. 根据权利要求9所述的方法（A1-A10)，其中，所述通过应用相应变换来计算（A6) 每个统计参数的值包括： 设置（A6)所述对数正态均值的值等于所述正态中间值的值的自然对数，并且如果所 述正态均值的值和所述正态中间值的值之间的比率严格地高于1，则设置所述对数正态标 准差的值等于所述比率的自然对数的两倍的平方根，否则，等于0。11. 根据权利要求1至10中的任何一个所述的方法（Al-AlO)，其中，所述分割（A7-A9) 选定位置包括： 根据每个统计参数的值和所述至少一个分割阈值中的相应的一个分割阈值之间的比 较，分割（A7)所述选定位置。12. 根据权利要求11所述的方法（Al-AlO)，其中，所述分割（A7-A9)选定位置包括： 当集中趋势统计参数的值高于相应的分割阈值并且离差统计参数的值高于相应的 分割阈值时，将所述身体部分的包括将被分析的病变的分析区域的每个选定位置分配 (A7-A8)给所述分析区域的与低灌注部分相关的第一部分，否则，分配给第二部分； 所述方法还包括： 计算（AlO)所述分析区域的不包括在所述第一部分中的选定位置的响应的合并值。13. 根据权利要求12所述的方法（Al-AlO)，其中，所述分割（A7-A9)选定位置包括： 当所述集中趋势统计参数的值低于相应的分割阈值并且所述离差统计参数的值高于 相应的分割阈值时，将所述身体部分的将被与所述分析区域进行比较的、包括健康组织的 控制区域的每个选定位置分配（A7-A8)给所述控制区域的与高灌注部分相关的另一第一 部分，否则，分配给另一第二部分； 所述方法还包括： 计算（AlO)所述控制区域的不包括在所述另一第一部分中的选定位置的响应的另一 合并值。14. 根据权利要求13所述的方法（Al-AlO)，还包括： 根据所述合并值和所述另一合并值的组合，计算（AlO)所述分析区域相对于所述控制 区域的相对合并值。15. 根据权利要求11所述的方法（Al-AlO)，其中，所述造影剂是能够在所述患者内循 环并且能够基本固定在生物目标上的针对性的造影剂，所述分割（A7-A9)选定位置包括： 当集中趋势统计参数的值高于相应的分割阈值并且离差统计参数的值低于相应的分 割阈值时，将每个选定位置分配（A7-A8)给与固定的针对性的造影剂相关的第一部分，否 贝1J，分配给第二部分。16. -种计算机程序（600)，被配置为当所述计算机程序在计算系统（100)上执行时使 所述计算系统执行根据权利要求1至15中的任何一个所述的方法（Al-AlO)。17. -种计算机程序产品（130)，所述计算机程序产品包括包含计算机程序（600)的非 暂时性计算机可读介质，所述计算机程序可加载到计算系统（100)的工作存储器中，从而 将所述计算系统配置为执行根据权利要求1至15中的任何一个所述的方法（Al-AlO)。18. -种用于诊断成像应用中的系统（100)，所述系统包括： 用于提供被灌注造影剂的身体部分在非零分析时间段期间的表示的装置（603-621)， 对于所述身体部分的一组位置中的每个位置所述表示包括所述位置在所述分析时间段期 间对询问信号的响应的指示； 用于对于一组选定位置中的每个选定位置计算所述选定位置在所述分析时间段期间 的响应的统计分布的至少一个统计参数的值的装置（624-633)，所述一组选定位置包括所 有的所述位置或者其一部分；和 用于根据关于所述选定位置的所述至少一个统计参数的值与至少一个分割阈值之间 的比较分割所述选定位置的装置（636-648)。19. 根据权利要求18所述的系统（100)，其中，所述用于提供身体部分在非零分析时间 段期间的表示的装置（603-621)包括： 用于提供在所述分析时间段期间登记的所述身体部分的数字图像序列的装置（603)， 每个数字图像包括关于每个位置的强度值，所述强度值基于所述位置在所述分析时间段中 的相应采集时刻的响应； 并且其中，所述用于计算至少一个统计参数的值的装置（624-633)包括： 用于从沿着所述数字图像的选定位置的强度值计算所述至少一个统计参数的值的装 置（624)。20. 根据权利要求18或19所述的系统（100)，其中，所述至少一个统计参数是多个统 计参数。21. 根据权利要求20所述的系统（100)，其中，所述统计参数是集中趋势统计参数和离 差统计参数。22. 根据权利要求18至21中的任何一个所述的系统（100)，其中，所述用于计算至少 一个统计参数的值的装置（624-633)包括： 用于计算与所述响应相关的所述至少一个统计参数的值的装置（630)，所述响应被变 换以提供所述变换的响应在所述分析时间段期间的、具有歪斜形状的变换的统计分布。23. 根据权利要求22所述的系统（100)，其中，所述变换的统计分布是对数正态统计分24. 根据权利要求22或23所述的系统（100)，其中，所述用于计算至少一个统计参数 的值的装置（624-633)包括： 用于从所述选定位置在所述分析时间段期间的响应的统计分布计算至少一个另一统 计参数的值的装置（624);和 用于通过将根据所述变换的统计分布的相应变换应用于所述至少一个另一统计参数 的值来计算每个统计参数的值的装置（630)。25. 根据权利要求24所述的系统（100)，其中，所述用于通过应用相应变换来计算每个 统计参数的值的装置（630)包括： 用于将对数变换应用于所述至少一个另一统计参数的值的装置（630)。26. 根据权利要求24或25所述的系统（100)，其中，所述至少一个统计参数是对数正 态均值和对数正态标准差，并且所述至少一个另一统计参数是正态均值和正态中间值。27. 根据权利要求26所述的系统（100)，其中，所述用于通过应用相应变换来计算每个 统计参数的值的装置（630)包括： 用于以下操作的装置（630):设置所述对数正态均值的值等于所述正态中间值的值的 自然对数，并且如果所述正态均值的值和所述正态中间值的值之间的比率严格地高于1，则 设置所述对数正态标准差的值等于所述比率的自然对数的两倍的平方根，否则，等于0。28. 根据权利要求18至27中的任何一个所述的系统（100)，其中，所述用于分割选定 位置的装置(636-648)包括： 用于根据每个统计参数的值和所述至少一个分割阈值中的相应的一个分割阈值之间 的比较分割所述选定位置的装置（636、645)。29. 根据权利要求28所述的系统（100)，其中，所述用于分割选定位置的装置（636、 645)包括： 用于以下操作的装置（645):当集中趋势统计参数的值高于相应的分割阈值并且离差 统计参数的值高于相应的分割阈值时，将所述身体部分的包括将被分析的病变的分析区域 的每个选定位置分配给所述分析区域的与低灌注部分相关的第一部分，否则，分配给第二 部分； 所述系统还包括： 用于计算所述分析区域的不包括在所述第一部分中的选定位置的响应的合并值的装 置（651)。30. 根据权利要求29所述的系统（100)，其中，所述用于分割选定位置的装置（636、 645)包括： 用于以下操作的装置（645):当所述集中趋势统计参数的值低于相应的分割阈值并且 所述离差统计参数的值高于相应的分割阈值时，将所述身体部分的将被与所述分析区域进 行比较的、包括健康组织的控制区域的每个选定位置分配给所述控制区域的与高灌注部分 相关的另一第一部分，否则，分配给另一第二部分； 所述系统还包括： 用于计算所述控制区域的不包括在所述另一第一部分中的选定位置的响应的另一合 并值的装置（651)。31. 根据权利要求30所述的系统（100)，还包括： 用于根据所述合并值和所述另一合并值的组合计算所述分析区域相对于所述控制区 域的相对合并值的装置（651)。32. 根据权利要求28所述的系统（100)，其中，所述造影剂是能够在所述患者内循环 并且能够基本固定在生物目标上的针对性的造影剂，所述用于分割选定位置的装置（630、 645)包括： 用于以下操作的装置（645):当集中趋势统计参数的值高于相应的分割阈值并且离差 统计参数的值低于相应的分割阈值时，将每个选定位置分配给与固定的针对性的造影剂相 关的第一部分，否则，分配给第二部分。33. -种诊断方法，包括： 对患者施用造影剂以使所述造影剂灌注所述患者的身体部分； 在非零分析时间段期间将询问信号应用于所述身体部分； 采集所述身体部分在所述分析时间段期间的表示，其中，对于所述身体部分的一组位 置中的每个位置，所述表示包括所述位置在所述分析时间段期间对所述询问信号的响应的 指示，所述表示根据权利要求1至15中的任何一个所述的方法进行处理以分割一组选定位 置，所述一组选定位置包括所有的所述位置或者其一部分；和 根据所述分割评估所述身体部分的状况。
【专利摘要】提出了一种用于诊断成像应用中的分割方案。相应的数据处理分割方法(A1-A10)包括：提供(A1-A4)被灌注造影剂的身体部分在非零分析时间段期间的表示，对于所述身体部分的一组位置中的每个位置所述表示包括所述位置在所述分析时间段期间对询问信号的响应的指示；对于一组选定位置中的每个选定位置，计算(A5-A6)所述选定位置在所述分析时间段期间的响应的统计分布的至少一个统计参数的值，所述一组选定位置包括所有的所述位置或者其一部分；并且根据关于所述选定位置的所述至少一个统计参数的值与至少一个分割阈值之间的比较，分割(A7-A9)所述选定位置。
【IPC分类】G06T7/00
【公开号】CN104885121
【申请号】CN201380066958
【发明人】P·弗林京
【申请人】博莱科瑞士股份有限公司
【公开日】2015年9月2日
【申请日】2013年12月18日
【公告号】EP2936433A1, US20150348277, WO2014096041A1
最新回复(0)