专利名称：组织恶性肿瘤检测方法、组织恶性肿瘤检测装置的制作方法以往以来，开发了检测组织的癌性肿瘤的许多技术(例如，专利文献1 7)。作为用来对组织的灌注进行摄影的造影剂的例子，可以举出在计算机断层摄影法(CT)中使用的碘造影剂、在磁共振影像法(MRI)中使用的钆螯合剂、在正电子发射断层摄影法(PET)中使用的放射性医药品(参照专利文献4、5)。CT和MRI典型地讲，利用起因于肿瘤性血管的透过性提高而造成的造影剂迅速积存到血管外的空间的性质来检测癌性组织。PET利用因肿瘤性细胞的代谢需要增加而造成的葡萄糖分子的吸收增加的性质来检测肿瘤性组织。超声波造影剂能够利用更大的典型地是填充了气体的泡来进行肿瘤性的微小循环的摄影。但是，在这些方法中，需要造影剂的侵袭性注入。进而，这些方法仅依存于癌性肿瘤部位的微小循环的形态构造，而没有利用微小血管的功能性特征、例如血流及博动。所以，作为不使用造影剂的侵袭性注入而根据癌性肿瘤自身的血液灌注来检测癌性肿瘤的别的方法，已知有利用超声波多普勒摄像的方法(另外，超声波多普勒摄像中的造影剂的使用也非常普遍)。在专利文献1中，公开了将组织密度与血流比较组合来检测初始的前列腺癌的系统。在专利文献2中，公开了使用能量多普勒摄像计算在肿瘤内新增殖的微小血管的程度、检测位置完全清楚的癌性肿瘤的恶性度和转移能力的方法。在专利文献3中，公开了基于向组织的血液灌注通过热而增加的前提、通过将组织内的血液灌注在安放热源前后进行比较来检测组织的异常的方法。专利文献1 美国专利申请公开第2003/0149358号说明书专利文献2 美国专利申请公开第2006/0M1463号说明书专利文献3 美国专利第5，233，994号说明书专利文献4 美国专利第4，276，885号、国际公开第80/02365号专利文献5 国际公开第01/28594号专利文献6 美国专利申请公开第2004/0199077号说明书、美国专利第6，984，211 号说明书专利文献7 美国专利申请公开第2005/0027188号说明书、美国专利第7，466，848 号说明书发明概要发明要解决的问题但是，这些方法都只是在血流速度的测量中利用多普勒摄像，而没有利用例如时间域中的组织的博动波形(以下也称作博动曲线)的特征。即，仅利用了与血管新生相关联地从组织得到的信息中的一部分。通过肿瘤的血管新生，产生出具有与正常的组织的缺陷构造不同的构造的新的血管。与癌性肿瘤关联的这些新的血管具有血管周围的分离、血管扩张、及不规则的形状的特征。肿瘤性的血管不具有存在于正常的组织的血管中的平滑肌细胞，此外，没有充分形成到对全组织给与氧的程度。进而，肿瘤性血管还与正常的组织的血管相比为多孔质，容易渗出ο但是，在以往技术中，尝试了仅根据通过基于多普勒摄像的血流速度的测量来检测癌性肿瘤的博动，关于使用能够充分捕捉癌性肿瘤的博动所具有的上述特征的特征量进行检测的方法没有公开。
所以，本发明的目的是提供一种通过更适当地检测癌性肿瘤的博动所具有的特征，从而更正确地判断恶性肿瘤的组织恶性肿瘤检测方法。解决问题的手段有关本发明的一技术方案的组织恶性肿瘤检测方法，通过用超声波扫描组织而得到的扫描信号检测包含在上述组织中的恶性肿瘤，包括块划分步骤，将对上述组织进行了扫描的区域划分为多个块；组织博动推测步骤，按上述多个块中的每个块，基于上述扫描信号推测组织博动，该组织博动是通过上述组织进行博动而发生的上述组织的变位的时间变化；博动关联特征提取步骤，按上述多个块中的每个块，从上述组织博动中提取作为关于上述组织的博动的参数的多个博动关联特征；分布特性计算步骤，按上述多个块中的每个块， 计算上述多个博动关联特征的分布特性；以及恶性分类步骤，基于上述分布特性，对上述多个块中的各个块是否是作为包含恶性肿瘤的块的恶性块进行分类。由此，能够按照包含在扫描区域中的每个块，根据表示组织的变位的时间变化的信息，计算与因恶性肿瘤造成的血管新生给组织的博动带来的多个影响分别对应的多个特征量。能够根据这样计算出的多个特征量将包含在扫描区域中的块分类为包含恶性肿瘤的准确度较高的恶性块、和除此以外的块。结果，通过更适当地检测癌性肿瘤的博动具有的特征，能够更正确地判断恶性肿瘤。另外，本发明不仅能够作为这样的组织恶性肿瘤检测方法实现，还能够作为组织恶性肿瘤检测装置实现、或使计算机执行组织恶性肿瘤检测方法中的特征性的步骤的程序实现。并且，当然可以使这样的程序经由CD_R0M(Compact Disc Read Only Memory)等的记录介质及因特网等的传送介质流通。进而，本发明可以作为实现这样的组织恶性肿瘤检测装置的功能的一部分或全部的半导体集成电路(LSI)实现、或作为包括这样的组织恶性肿瘤检测装置的组织恶性肿瘤检测系统实现。发明效果能够提供一种通过更适当地检测癌性肿瘤的博动具有的特征，从而更正确地判断恶性肿瘤的组织恶性肿瘤检测方法。图1是表示有关本发明的实施方式的组织恶性肿瘤检测装置的功能块的图。图2是表示有关本发明的实施方式的组织恶性肿瘤检测装置的处理的流程的流程图。图3是表示有关本发明的实施方式的组织博动推测部的结构的一例的图。图4是表示本发明的实施方式的超声波原RF数据的规格的一例的图。图5是表示有关本发明的实施方式的前处理部的结构的一例的图。图6是表示有关本发明的实施方式的主要成分提取部的结构的一例的图。图7是表示有关本发明的实施方式的极值识别部及博动调整部的处理结果的一例的图。图8是表示有关本发明的实施方式的振幅特征计算部的结构的一例的图。图9是表示有关本发明的实施方式的形状特征计算部的结构的一例的图。图10是表示有关本发明的实施方式的心周期识别部的结构的一例的图。图11是表示有关本发明的实施方式的临界点识别部进行的处理的一例的图。图12是表示本发明的实施方式的心博动的博动形状特征的一例的图。图13是表示肿瘤性组织的博动与正常的组织的博动的博动形状特征的比较的图。图14是表示有关本发明的实施方式的主要成分提取部的结构的一例的图。图15是表示有关本发明的实施方式的主要成分提取部的信号的分段、减缩及代表提取的各处理的结果的一例的图。图16是表示有关本发明的实施方式的空间分布计算部的结构的一例的图。图17是表示本发明的实施方式的灰度级同时发生概率矩阵的计算处理的一例的图。图18是表示本发明的实施方式的肿瘤位置确定部的结构的一例的图。图19是说明本发明的实施方式的肿瘤位置确定部进行的处理的一例的图。图20是表示本发明的实施方式的心周期延迟计算部的结构的一例的图。图21是表示本发明的实施方式的延迟计算部的结构的一例的图。图22是表示本发明的实施方式的自回归系数计算部的结构的一例的图。图23是表示在组织中没有肿瘤的情况下的有关本发明的实施方式的组织恶性肿瘤检测装置的输出结果的一例的图。图M是表示在组织中有肿瘤的情况下的有关本发明的实施方式的组织恶性肿瘤检测装置的输出结果的一例的图。图25是表示实现有关本发明的实施方式的组织恶性肿瘤检测装置的计算机系统的硬件结构的块图。
此外，也可以是，在上述组织博动推测步骤中，以包含在上述多个块中的各块的1 个或几个代表性博动、或上述代表性博动的组合，来推测上述组织博动。由此，能够将从多个扫描点得到的数据综合判断，来判断恶性肿瘤的有无。此外，也可以是，上述博动关联特征是作为关于上述组织博动的振幅的特征量的博动振幅特征、作为关于上述组织博动的波形的形状的特征量的博动形状特征、以及作为关于上述组织博动的波形的时间变化的特征量的博动时间特征中的至少1个。由此，能够根据博动波形的振幅、形状、时间变化等多个特征量进行更正确的判断。更具体地讲，也可以是，上述博动形状特征是以下中的至少1个L1/L2，该L1/L2 是作为心周期的收缩部分的期间的收缩期间(Li)与作为心周期的舒张部分的期间的舒张期间(L2)的比率；L3/L4，该L3/L4是收缩中间点及舒张中间点之间的期间(L3)与心博动的周期(L4)的比率，其中上述收缩中间点是在上述收缩期间中振幅相对于最大振幅等于预先定义的比率的收缩曲线上的点，上述舒张中间点是在上述舒张期间中振幅相对于上述最大振幅等于预先定义的比率的舒张曲线上的点；上述舒张曲线与连结收缩峰和舒张结束点的预先定义的曲线之间的偏差，其中上述收缩峰是上述收缩期间中的振幅的波峰，上述舒张结束点是上述舒张期间的结束点；偏度，表示心博动的非对称性；峰度，表示上述收缩峰的尖锐度；和上述舒张曲线中包含的极值的偏差。此外，也可以是，上述博动形状特征在形状特征计算步骤中计算，该形状特征计算步骤包括根据上述组织博动计算心周期的心博期间计算步骤；使用上述心周期识别临界点的临界点识别步骤；以及基于上述心周期和上述临界点提取上述博动形状特征的形状特征提取步骤。此外，也可以是，上述心周期的上述临界点包括收缩开始点，该收缩开始点是心周期的收缩部分的开始点；舒张结束点，该舒张结束点是心周期的舒张部分的结束点；收缩峰，该收缩峰是上述收缩期间中的振幅的波峰；收缩中间点，该收缩中间点是在心周期的收缩部分的期间中振幅相对于最大振幅等于预先定义的比率的收缩曲线上的点；和舒张中间点，在心周期的舒张部分的期间中振幅相对于最大振幅等于预先定义的比率的舒张曲线上的点。此外，也可以是，上述临界点识别步骤包括检索步骤，检索上述组织博动中的极小点和极大点；博动方向识别步骤，基于上述极小点和极大点，识别表示上述组织博动具有向上的收缩曲线、还是具有向下的收缩曲线的博动方向；以及临界点决定步骤，使用上述极大点、上述极小点、和上述博动方向，求出上述博动波形中的上述临界点。此外，也可以是，上述预先定义的曲线是直线。此外，也可以是，在博动关联特征提取步骤中，在上述博动具有向上的收缩曲线的情况下，计算上述舒张曲线与上述预先定义的曲线之间的正的差分和，来作为上述偏差，在上述博动具有向下的收缩曲线的情况下，计算上述舒张曲线与上述预先定义的曲线之间的负的差分和，来作为上述偏差。由此，通过根据博动处于舒张期还是处于收缩期计算不同的特征量，能够进行更正确的判断。此外，也可以是，上述博动时间特征是上述心周期的延迟、作为上述心博动的波形的心波形的差异、和上述组织博动的波形的自回归系数中的至少1个。具体而言，也可以是，上述心周期的延迟在心周期延迟计算步骤中计算，该心周期延迟计算步骤包括决定基准心周期的基准心周期决定步骤，该基准心周期是作为基准的心周期；以及计算对象心周期相对于上述基准心周期的延迟的延迟计算步骤。更具体地讲，也可以是，在上述基准心周期决定步骤中，选择扫描数据之中的振幅最大的心周期，作为上述基准心周期。此外，也可以是，上述基准心周期根据心电图检查ECG信号决定。此外，也可以是，上述心波形的差异在延迟计算步骤中计算，该延迟计算步骤包括计算基准心波形的基准心波形计算步骤，该基准心波形是作为基准的心波形；计算博动产生的多个心波形中的各个心波形与上述基准心波形之间的差异的个别差异计算步骤； 以及根据计算出的多个上述差异，计算作为代表上述多个心波形与上述基准心波形之间的差异的值的整体心波形差异值的整体差异计算步骤。此外，也可以是，上述整体心波形差异值是计算出的多个上述差异的标准偏差。此外，也可以是，上述自回归系数在自回归系数计算步骤中计算，该回归系数计算步骤包括将多个上述博动波形减缩以使上述心周期成为相同的期间的博动再采样步骤； 以及基于事前设定的自回归模型和减缩后的上述博动波形求出作为上述自回归模型具有的系数的自回归系数的自回归式运算步骤。此外，也可以是，上述分布特性是空间性分布参数和统计性分布参数中的至少1 个。由此，能够通过考虑到空间性的位置关系的分布、和考虑到统计性的特性的分布、 进行更正确的判断。具体而言，也可以是，上述空间性分布参数是上述博动关联特征的能量、上述博动关联特征的熵、上述博动关联特征的对比度、上述博动关联特征的均质性、和上述博动关联特征的相关性中的至少1个。此外，也可以是，上述统计性分布参数是上述博动关联特征的平均值、上述博动关联特征的中位数、上述博动关联特征的最大值、上述博动关联特征的最小值、上述博动关联特征的标准偏差、上述博动关联特征的峰度、和上述博动关联特征的偏度中的至少1个。此外，也可以是，上述博动关联特征及其分布特性是以下中的至少1个在上述多个块中的各块中包含的多个扫描点的博动振幅的中位数、熵、标准偏差、平均值、及最大值； 在上述多个块中的各块中包含的多个扫描点的心波形差分的中位数、熵、标准偏差、平均值、及最大值；在上述多个块中的各块中包含的多个扫描点的收缩期间与舒张期间的比率的中位数，其中上述收缩期间是心周期的收缩部分的期间，上述舒张期间是心周期的舒张部分的期间；上述多个块中的各块的代表性博动的上述收缩期间与上述舒张期间的比率； 在上述多个块中的各块中包含的多个扫描点的舒张曲线的偏差的最大值及标准偏差。具体而言，也可以是，上述肿瘤位置确定步骤还包括对象区域确定步骤，按照所扫描的上述组织的每个扫描点规定对象区域；肿瘤块划分步骤，确定上述多个块之中的、属于上述对象区域的块；以及癌概率计算步骤，根据属于上述对象区域的块的上述恶性分类步骤的分类结果，计算上述组织是癌的概率。此外，也可以是，上述肿瘤位置确定步骤还包括将所扫描的上述组织的扫描点的上述是癌的概率用2维或3维图像显示的图像化步骤。有关本发明的另一技术方案的组织恶性肿瘤检测装置，通过用超声波扫描组织而得到的扫描信号检测包含在上述组织中的恶性肿瘤，具备块划分部，将对上述组织进行了扫描的区域划分为多个块；组织博动推测部，按上述多个块中的每个块，基于上述扫描信号推测组织博动，该组织博动是通过上述组织进行博动而发生的上述组织的变位的时间变化；博动关联特征提取部，按上述多个块中的每个块，从上述组织博动中提取作为关于上述组织的博动的参数的多个博动关联特征；分布特性计算部，按上述多个块中的每个块，计算上述多个博动关联特征的分布特性；以及恶性分类部，基于上述分布特性，对上述多个块中的各个块是否是作为包含恶性肿瘤的块的恶性块进行分类。根据该结构，能够将癌性肿瘤的博动具有的特征根据多个特征量更适当地检测癌性肿瘤的博动具有的特征。结果，能够提供能够更正确地判断恶性肿瘤的组织恶性肿瘤检测装置。此外，也可以是，还具备根据在上述恶性分类部中被分类为是上述恶性块的块来确定癌性肿瘤的位置的肿瘤位置确定部。由此，能够简单且正确地确定恶性肿瘤的位置。作为以上所述的本发明的前提，首先，对由癌性肿瘤造成的血管新生、及由此发生的癌性肿瘤的博动详细地说明。癌性肿瘤需要形成用来向自身供给营养和氧、并且将代谢作用的结果生成的排泄物除去的、新的血管。将这样的形成新的血管的过程称作血管新生(angiogenesis)。周知的是，在癌性肿瘤的初期阶段中，当直径超过约2 3毫米而成长时需要血管形成。这些新的血管在肿瘤的周围及内部成长，此外，随着血管的增加，癌细胞达到体内的主要的血管， 结果通过这些新形成的血管转移的可能性也增加。通过检测血管新生，有可能有助于初期阶段的癌性肿瘤的检测，有可能有助于处置后的转移概率的预测。这里，癌性肿瘤(也称作恶性肿瘤)在许多方面与正常的组织或良性肿瘤不同。首先，在毛细血管床的级别发生血管新生。通过心脏将血液向身体的各种部分送出，产生组织内的博动，结果，通过血液的灌注，组织博动。随着向毛细血管的血液供给的增力口，可以预想到与周围的正常的组织相比，癌性组织的博动振幅会增加。同时，还可以预想到向癌性肿瘤的血液的出入也增加。除了振幅以外，与癌性肿瘤关联的这些新的血管还具有血管周围的分离、血管扩张、及不规则的形状的特征。肿瘤性的血管不具有存在于正常的组织的血管中的平滑肌细胞，此外，没有充分地形成到对全组织给与氧的程度。进而，肿瘤性血管还容易比正常的组织的血管更多孔质而容易渗出。通过这些全部的因素，癌性肿瘤的脉搏的到达时间、脉搏波形、因呼吸带来的脉搏振幅的变化、及因呼吸带来的基准值变化与正常的组织中的规则性的变化不同，有可能为不规则的。如果空间性的位置不同，则因为博动样式(pattern)不同等的理由，预测癌性肿瘤的博动波形具有比正常的组织的博动更加扩张的形状。微小血管系统中的上述差异也成为博动波形的出现定时的差异，表示时间变动的特性。此外，这也成为统计性且空间性表示的博动的分布特性的差异。由此，通过识别随着血液灌注而从组织检测的博动的差异，能够对癌性肿瘤组织和正常的组织进行分类。具体而言，在本发明的实施方式中，利用作为来自博动的主要的参数的(1)博动振幅、(2)博动形状、(3)博动时间特征、以及(4)它们的分布特性作为特征量，以便能够判断扫描的组织的恶性度。如果作为血管新生的结果而向毛细血管的血液供给和微小血管系统增加，则在癌性组织和正常的组织中有可能在博动振幅中发生差别。肿瘤的血管构造比正常的血管构造更多孔质而容易渗出，在癌性微小血管中没有平滑肌，所以它们成为血管扩张状态。由此，有可能发生通过从波形提取的特征——即Li/ L2、L3/L4、偏度、峰度、舒张曲线的偏差、及舒张曲线的极值的偏差测量出的关于博动形状的差异。这些在实施方式中进一步说明。另外，L1/L2是博动的收缩期间(Li)与舒张期间(L2)的比率，L1/L2及偏度表示收缩峰值以何种程度延迟。延迟越大，该博动来源于癌性肿瘤的可能性越高。L3/L4是例如一半的高度下的博动幅度(U)与全博动期间(L4)的比率，L3/L4及峰度表示博动波形以何种程度宽度扩大且扩张。越是扩张，则是癌性博动的可能性越高。并且，舒张曲线的偏差及舒张曲线的极值的偏差表示博动波形中的重复隆起的有无及舒张曲线的曲率。如果没有重复隆起，则表示有可能是癌性博动。博动时间特征是表示搏动怎样随着时间变化的关于博动波形的时间变化的特征量。具体而言，是心周期(也称作心博期间)的延迟、心波形的差异、及自回归系数。进而，特定的组织区域中的上述特征的空间性且统计性的特性对于因微小血管构造的差异而表现基本的博动分布是有效的，表示博动在区域内怎样分布。能够选择利用这些特征而分类恶性组织和良性组织。并且，可以考虑将分类结果组合，计算被扫描的组织的内部区域的恶性分数 (score)，通过这些分数将癌性肿瘤存在于扫描的组织中的概率表示给医生或患者等。对以上叙述的本发明的前提、及有关本发明的特征量的一例说明了概要。接着，参照附图对有关本发明的组织恶性肿瘤检测装置及组织恶性肿瘤检测方法的实施方式詳细地说明。另外，以下的实施方式仅是说明各种发明步骤的原理的。应理解的是，这里说明的具体例的各种各样的变形对于本领域的技术人员而言是显而易见的。在本实施方式中说明的方法有时使用2维数据表示实施方式的步骤及效果。但是，这些数据单纯以例示为目的，这种方法的3维数据下的实施方式对于本领域的技术人员而言是显而易见的。3维数据下的这样的实施方式也还是本发明的范围内。医疗用图像设备也可以具备以下的组件的一些或全部。信号收发机、数据取得、数据处理、及显示器。本发明指示在医疗用图像设备的数据处理组件中实施的方法。这样的数据处理组件的输入是从数据取得组件取得的扫描信号，将输出通过显示器组件显示给用户。应理解的是，这样的医疗用图像设备的他的组件适合于使用本发明的应用例。本发明的主要的实施例在图1中表示。图1是表示有关本发明的实施方式的组织恶性肿瘤检测装置90的功能块的图。这里，scarmedSig表示从在被扫描的组织中检测博动而使用的诊断图像设备接收到的扫描数据。作为这样的设备的例子，可以举出scarmedSig为高频(RF)数据或将该RF 数据解调而得到的表示正弦波的振幅及相位的变化的IQ数据的、医疗用超声波设备。超声波设备的动作以朝向被扫描的介质发送高频的超声波脉冲为基础，然后，脉冲一边反射及散射一边沿着进路与处于下层的构造相互影响，接收反射及散乱的信号。即，scarmedSig是表示朝向作为扫描对象的组织发送的超声波脉冲的反射波的信号。使用这些信号，推测处于下层的构造的形态学的特征及功能性特征。作为医疗用图像技术而使用超声波的优点之一是具有生成高帧速率图像、并且能检测组织内的较小的运动的能力。根据图1所示的主要的实施方式，组织恶性肿瘤检测装置90是通过用超声波将组织扫描而得到的扫描信号检测包含在组织中的恶性肿瘤的组织恶性肿瘤检测装置，更详细地讲，具备块划分部100、组织博动推测部101、前处理部102、博动关联特征提取部103、分布特性计算部104、和恶性分类部105。另外，前处理部102是可选的，出于例示的目的而包含。无前处理部102的实施方式也是本发明的范围内。块划分部100根据应用例，表示将通过超声波将组织扫描的区域——即扫描区域划分为特定的形状及尺寸的组织块(以后，也称作块)的处理。块的形状可以选择方形、长方形、多边形、或者圆，但并不限定于这些。尺寸可以选择1毫米，但并不限定于此。组织博动推测部101根据接收到的scarmedSig推测基础性的已扫描组织的博动 (以后，也称作组织博动)。即，组织博动推测部101对多个块中的每个块，分别基于扫描信号推测作为通过组织博动而发生的组织的变位的时间变化一即组织博动。这里，假设在各块中包含1个以上的扫描信号。组织博动推测部101的输出是作为表示未加工的组织博动的信号的rawPulse。为了表示组织内的博动而能够利用的参数包括组织变位和组织形变，但并不限定于这些。组织变位被按照零点时的原来的位置推测为组织的绝对运动。组织形变被推测为组织内的部分间的相对的变位。即，组织形变是组织变位的空间的梯度。可以选择特定的组织点(即， 从扫描对象的组织的某个地点取得的扫描信号)的、相对于时间的组织变位或相对于时间的组织形变当中的某个作为表示该点的组织博动的参数。即，所谓组织博动，是组织的空间位置的变位的时间变化。另外，组织博动也可以包括组织形变和组织变位的两者。前处理部102将推测的原博动中的干扰(所谓的噪声)排除。未加工的组织博动有包括各种各样的种类的组织的运动的情况。这样的运动之一，是作为按照心脏的心周期在血管中流动的血流的结果而产生的心博动。作为别的运动，有比心周期慢数倍的、通常振幅比心博动大幅增大的患者的呼吸周期。该运动在本发明中是不重要的。还有发生其他环境的干扰或电子的干扰的情况。前处理部102输出抑制了干扰的纯粹的博动、即 procPulse。另外，在包含在从组织博动推测部101输出的rawPulse中的干扰较少的情况下，即使不具备前处理部102，组织恶性肿瘤检测装置90也起到同样的效果。例如，在扫描对象者处于安静状态的情况下等，即使不具备前处理部102，博动关联特征提取部103也能够进行特征提取。但是，具备前处理部102能够更正确地检测恶性肿瘤。博动关联特征提取部103按多个块中的每个块，分别进行用来根据组织博动提取关于时间域中的博动波形的参数——即多种博动关联特征的计算处理。该块的输出是 pulsei^eature，为向分布特性计算部104及恶性分类部105的输入信号。分布特性计算部104计算关于所提取的特征的、作为各块中的每个块的代表特性的分布特性的参数。该分布特性计算部104的输出是distributior^eature，它也为向恶性分类部105的输入。
恶性分类部105计算为了将组织内的区域分类为良性或恶性而需要的信息，按多个块中的每个块，进行是否是作为包含恶性肿瘤的块的恶性块的分类。该恶性分类部105 的输出是maligScores。参照图2说明以上叙述的组织恶性肿瘤检测装置90的处理的流程。图2是表示组织恶性肿瘤检测装置90的处理的流程的流程图。首先，如果从外部设备取得组织的扫描信号，则块划分部100将由扫描信号表示的扫描区域划分为块(SlOO)。接着，组织博动推测部101按多个块中的每个块，分别推测组织博动(SlOl)。接着，博动关联特征提取部103根据推测的组织博动提取表示时间域的组织博动的特征的博动关联特征(S103)。接着，分布特性计算部104在扫描区域中，按多个块中的每个块，计算表示博动关联特征的值怎样分布的分布特性(S104)。最后，恶性分类部105判断位于与各块对应的组织内的区域中的肿瘤的恶性度， 将肿瘤分类为良性或恶性(S105)。接着，再次参照图1，表示基于图1所示的主要的实施方式的一些具体的实施方式。图1的块划分部100将通过超声波将组织扫描的区域划分为更小的区域——即组织块。进一步的处理全部对该组织块进行。组织块的形状及尺寸根据实际的应用例决定。 形状可以选择方形、长方形、或圆，但并不仅限定于这些。尺寸可以选择1毫米。根据发明者的实验，约1 2毫米尺寸的长方形组织块是有效的。块划分的输出是作为表示组织块的结构的信息的blocks。图1的组织博动推测部101由具备超声波换能器的图像装置将组织用超声波扫描，使用作为输出的扫描信号的超声波原RF数据，根据组织推测未处理(即，实施干扰的除去处理前)的博动。作为组织博动推测部101的一例，可以举出图3所示的结构，但并不限定于此。在图3中，说明推测由使用超声波的组织形变表示的组织博动的结构。如图3所示，有关本实施方式的组织博动推测部101具有组织变位推测部200和组织形变推测部201。这里，参照图4，如果设想超声波换能器的扫描方向朝向组织(垂直于身体的表面的方向)，则作为超声波换能器的输出的超声波原RF数据表示为scarmedSigW，1，f)。这里，d表示与扫描方向平行的深度方向，此外也称作快时间(fast-time)，1表示与扫描方向正交的线(line)方向，f表示时间或帧，此外也称作慢时间(low-time)。组织变位推测部200作为输入而取得scarmedSig (d，1，f)，计算表示组织变位的 disp(d, l，f)。在根据未加工的超声波原RF数据计算变位的各种方法中，包括自相关及互相关，但并不限定于这些。组织形变推测部201作为输入而取得表示组织变位的disp (d，1，f)，计算组织形变。作为计算组织形变的方法的一例，可以举出计算沿着深度方向的变位的空间的梯度、即 disp(d, l，f)相对于d的导数。在对时间标绘的特定的深度d和特定的线1中计算出的组织形变表示某个组织点的组织博动raWPulse(d，1，f)。各组织点由对应的深度及行记述。由此，它通过 rawPulse(d,l,f)的特定的博动波形记述。该博动波形也包含在组织中可能发生的任何种类的运动及干扰。图1的前处理部102通过对未处理的博动rawPulse实施进一步的处理，能够将干扰除去、使重要的博动的提取变得容易。前处理部102是可选的，不包含该块的实施方式也还是包含在本发明的范围中。前处理的输出是procPulse。在图5中表示这样的前处理部102的结构。如图5所示，前处理部102具有噪声滤波部400、和主要成分提取部401。另外，前处理部102的结构并不限定于这些。这是单纯以例示为目的的，这些结构部执行的步骤及其子步骤以任何顺序或组合执行都可以，这些步骤都能够不给其他步骤带来影响而可以省略，这些步骤都能够包含在图1的其他块的任何之中，这样的实施方式还是本发明的范围内。噪声滤波部400是有可能被实施的前处理部102具有的构成单元之一。心博动的频率是1 2Hz，呼吸周期的频率是不到0. 5Hz。由此，通过将噪声滤波部400构成为仅使心博动的频率1 2Hz通过的带通滤波器，能够将环境的高频的振动噪声、及低频的呼吸运动的两者高效率地除去而不会给心博动的特征带来影响。主要成分提取部401是另一有可能被实施的前处理部102具有的构成单元的例子，提取所选择的重要的区域中的主要博动成分(component)(即，包含在组织博动中的博动成分当中的主要的成分)。这在噪声滤波部400不能完全抑制干扰、特别是身体的运动那样的不规则的干扰的情况下是有效的。该块的输入也可以是未加工的博动、或已经实施了其他前处理步骤的博动。同样，该块的输出既可以为向其他前处理步骤的输入，也可以作为向博动关联特征提取部103、及分布特性计算部104的输入而原样使用。在1实施方式中，主要成分提取部401通过使用主成分分析(PCA;)或独立成分分析(ICA;)能够执行主要成分的提取，但并不限定于这些。这里，设想所谓主要成分，是心博动和呼吸运动，其他全部的成分是干扰。PCA或ICA取得预先定义的区域的博动，在对于其他成分(即干扰)不相关(PCA的情况)或独立(ICA的情况)的博动中取得最有力的成分。这伴随着削减空间分辨率的代价而执行。这样，能够选择重要的主要成分。图6是表示主要成分提取部401的另一例的块图。这里，假设重要的博动仅包含心博动。即，主要成分提取部401从组织博动中提取心博动。进而，在图6中，主要成分提取部401具有心功率计算部500、心博聚类部501、极值识别部502、和博动调整部503。心功率计算部500根据作为未处理的博动的rawPulse (d，1，f)，计算与心博动关联的信号功率、即cardiaCP0W(d，l)。另外，假设心博动的频率是已知的(用fc表示，在人的情况下通常是1 2Hz)。在实施方式的一例中，计算CardiaCP0W(d，1)作为与fc关联的未加工的博动的功率谱成分的值。在实施方式的另一例，作为CardiaCP0W(d，l)，作为与用全成分的功率的整体和或部分和进行归一化的fc关联的未加工的博动的功率谱成分的值计算。心博聚类(cluster)部501识别振幅或心博动的存在较高的组织内的区域、即 cardiacCluster(d, 1)。在实施方式的一例中，将cardiacCluster(d，1)作为具有包含在预先定义的范围中的大小的心功率的组织内的区域而被识别。在实施方式的另一例中，将cardiacCluster(d, 1)确定为在组织内包含心功率更高的区域，以使其相对于组织整体占预先定义的比例。极值识别部502决定由cardiacCluster (d，1)确定的组织内的区域中的博动波形的波峰及波谷(即极值点)，输出属于心博动的这样的波峰及波谷、即cardiacExtremaW， 1)。即，极值识别部502识别聚类后的组织内的区域中的心博动的振幅的极值。在图7中表示极值识别部502的处理。求出波峰及波谷的实施方式的一例是推测博动波形的相对于时间的一次导数和二次导数。表示为pulsePeaksW，1)的波峰被在一次导数是零、二次导数是负的时间实例中被识别，此外，表示为pulseTr0Ughs(d，1)的波谷被在一次导数是零、二次导数是正的时间实例中被识别。波峰的例子参照图7(A)，是点603、点604、点605、点606，波谷的例子是用点600、点607、点601、点602表示的点。设想属于心博动的波峰及波谷比干扰或其他运动的波峰及波谷更规则，所以作为执行该任务的方法，可以使用异常值丢弃，但并不限定于此。将异常值丢弃的情况下的心博的波峰和波谷的例子是点603、点607、点605、点606，干扰的波峰和干扰的波谷的例子是点 600、点604、点601、点602。即，极值识别部502也可以具有识别博动波形中的波峰和波谷的极值点识别部(未图示)、和为了将与其他波峰的偏差的大小比事前设定的阈值大的干扰波峰、以及与其他波谷的偏差的大小比事前设定的阈值大的干扰波谷删除而将博动波形的振幅的异常值丢弃的异常值丢弃部(未图示)。另外，所谓干扰的波峰(也称作干扰峰)，是指波峰的大小与其他波峰的大小相比，大得超过事前设定的阈值的波峰。此外，所谓干扰的波谷(也称作干扰谷)，是指波谷的大小与其他波谷的大小相比，大得超过事前设定的阈值的波谷。博动调整部503作为用来调整rawPulse (d，1，f)以使其仅具有心博动的信息而取得 cardiacCluster (d, 1)及 cardiacExtrema (d, 1)。该块的输出是 procPulse (d, 1)。博动调整部503的处理的效果在图7中表示。在博动调整部503的实施方式中，如图7(B)的曲线608所示，根据心脏的波峰及波谷识别仅对应于心博动的时间的长度。将对应于干扰的期间从博动中除外。在博动调整部503的另一实施方式中，参照图7(C)，还调整博动波形的振幅，以使心博动的波峰及波谷的振幅分别在时间方向上排列为一列(曲线611)。即，博动调整部503 调整博动波形的振幅的大小，以使博动波形的极小值及极大值的大小分别一致。此时，在使振幅一致时使用的峰振幅和谷振幅例如是用图7(B)的虚线609及虚线610表示的波峰及波谷的平均振幅，但并不限定于此。将不是波峰及波谷的其他全部的点通过线性映射的方法调整，但并不限定于此。在博动调整的另一实施方式中，如曲线608及曲线611那样，将片断化的心博时间组合。将结果得到的博动参照图7(D)表示为曲线612。通过这样的处理方法，保证片断化的心博期间的终点必须是与下个片断化的心博期间的起点相同的种类(即都是波峰或都是波谷)。在主要成分提取部401的其他实施方式中，使用分段、减缩、及平均化，取得组织块的代表性的博动。这在图14中表示。在图14所示的主要成分提取部401的实施方式的另一例中，主要成分提取部401具有心博动选择部1300、博动减缩部1301、和代表性博动提取部1302。心博动选择部1300以博动rawPulSe(d，1，f)为输入，基于一定的基准选择心博动。作为选择基准的一例，可以考虑心频率区域的信号能量比其他区域高的基准。作为选择基准的另一例，可以考虑心频率区域的信号波谱平坦度比阈值低的基准。在没有选择心博动、即输入的博动不满足选择基准的情况下，也可以将心博动设为零。心博动选择部1300 的输出是 cardiacPulse(d, 1，f)。博动减缩(日语歹一"丨J >夕’)部I3Ol以CardiaCRilSe(d，1，f)为输入，使各心周期逐渐变小，使这些心周期成为相同的期间以便能够正确地比较。该块的输出是 taperedCardiacPulse(d, 1，f)。代表性博动提取部1302以taperedCardiacPulse (d，1，f)为输入，根据博动波形的减缩后的心脉冲，计算该组织点的代表性的博动。作为这样的计算的一例，可以举出将减缩后的心脉冲遍及几个周期进行平均化的计算。作为这样的计算的另一例，可以举出使用 PCA或ICA的计算。将信号的分段、减缩、及代表提取的例子表示在图15中。首先，调整心脉冲以使波形振幅为全部是正或全部是负。并且，根据作为基准的博动(基准博动)的起点及终点，将不同的组织点的心脉冲划分为连续的心周期。将其结果表示为图15(A)的曲线1400。作为基准博动的一例，可以举出扫描面内的较强的心脉冲，作为基准博动的另一例，可以举出 ECG信号。在分段后，通过由窗函数使各心周期逐渐变小(减缩)而使其成为相同的期间。将其结果表示为图15(C)的曲线1402。作为该窗函数的一例，可以举出图15(B)所示的定制 Hann窗(customized hann window) 1401。另外，窗函数也可以使用汉明窗等其他任意的函数。最后，通过将这些减缩后的心脉冲遍及几个周期进行平均化，得到代表性的博动。 将得到的代表性的博动表示为图15(D)的曲线1403。将关于上述主要成分提取的3个实施方式分别使用或组合使用均可，这样的实施方式还是本发明的范围内。作为块划分部100、组织博动推测部101、及前处理部102的输出，已加工博动既可以是在包含在各组织块中的全扫描点(即，计测部位)计算的，也可以是各组织块的1个或几个代表性博动。此外，也可以是这样的代表性博动的组合。再次参照图1，图1的博动关联特征提取部103以pr0CPulSe(d，1)为输入，从那里取得作为关于时间域的博动波形的参数(即，特征量)的博动关联特征。这样的博动关联特征也可以是博动振幅特征、博动形状特征、及博动时间特征的全部、或它们的子集。将它们汇总称作pulseFeature。例如，博动关联特征提取部103也可以计算作为关于博动波形的振幅的特征量的博动振幅特征。癌性肿瘤需要形成用来供给营养和氧的新的血管。随着其在微小血管系统中的增加，通过血液的灌注，组织周边的博动振幅及其分布(统计性且空间性的分布)有可能比正常的组织高。在图8中表示计算博动的振幅特征的振幅特征计算部710的结构的一例。图8表示计算博动振幅特征的博动关联特征提取部103具备的振幅特征计算部710的功能块。如图8所示，振幅特征计算部710具有博动振幅计算部700和振幅直方图计算部701。在博动振幅计算部700中，计算通过d及1指定的各组织点处的心博动振幅、即 cardiacAmplitude(d, 1)。作为计算cardiacAmplitude (d，1)的方法，有使用博动的平均峰振幅与平均谷振幅的差的方法，但并不限定于此。在振幅直方图计算部701的一实施方式中，将在全部的组织点或组织块内计算出的cardiacAmplitude (d，l)分组为各种各样的预先定义的bin。振幅直方图计算部701在该博动振幅属于bin的组织点或组织块的情况下输出“1”、否则输出“0”。各bin的振幅范围可以分别设定，1个bin的范围不影响到其他bin的范围。可以使用cardiacAmp 1 itude (d， 1)的值和聚类结果的值的两者作为振幅特征计算部710的输出。这样的bin的一例是
，这在假定上与关于组织形变的肿瘤性博动振幅对应。另一例是
，这在假定上与关于组织形变的通常组织的博动振幅对应。也可以基于关于博动振幅的特定的假定来定义其他bin，这样的实施方式还是本发明的范围内。在振幅直方图计算部701的另一实施方式中，振幅直方图计算部701也可以将被扫描的组织整体划分为小区域，对各小区域计算振幅直方图，输出该直方图遵循预先定义的样式的小区域。组织的划分的一例是长方形的格子。作为振幅直方图的一例，可以举出遍及预先定义的bin的集合整体的博动振幅的概率分布。各直方图bin的范围可以分别设定，1个bin的范围不影响到其他bin的范围。作为预先定义的样式的一例，可以举出概率分布比
的范围偏斜。另外，作为振幅特征计算部710的输出pulseAmpFeatures，也可以选择博动的绝对值cardiacAmplitude (d，1)、或基于直方图的聚类结果ampHist (d，1)。此外，博动关联特征提取部103也可以计算作为关于组织博动的波形的形状的特征量的博动形状特征。在与癌性组织关联的新的血管中，有血管周围的分离、血管扩张、及不规则的形状。肿瘤性的血管不具有存在于正常的组织的血管中的平滑肌细胞，此外，没有充分形成为对整个癌性组织带来氧及营养的程度。进而，肿瘤性血管还比正常的组织的血管更多孔质而容易渗出。通过这些原因，肿瘤性的博动具有比正常的组织的博动更加扩张的形状。图9表示计算博动形状特征的博动关联特征提取部103具备的形状特征计算部 711的功能块。如图9所示，形状特征计算部711具有心周期识别部800和形状特征提取部 801。心周期识别部800根据procPulse (d，1，f)识别各心周期。能够从该块提取各心周期的临界点、即criticalPoints(d, 1)。在图10中表示心周期识别方法的1实施方式。 该实施方式的心周期识别部800具有心博期间计算部900和临界点识别部901。心博期间计算部900根据procPulse (d，1，f)计算输入博动的平均心博期间。另夕卜，心博期间计算部900也可以计算各输入博动的心周期。作为该计算的一例，可以举出将输入博动变换为频率域、检索基础性的频率成分、 选择它作为平均心博期间。另外，作为该运算的另一例，也可以在博动调整部503中计算两个邻接的波峰或两个邻接的波谷之间的平均距离、选择该平均距离作为平均心周期(即，平均心博期间)。心博期间计算部900的输出是cardiacPeriocKd，1)，这里，d表示深度，1表示线。临界点识别部901识别各测量点的心周期的临界点。在临界点中，可以考虑例如收缩开始点、舒张结束点、收缩峰、收缩中间点、及舒张中间点等，但并不仅限定于这些。该块的输入包括procPulse(d，1，f)及cardiacPeriod(d，1)，该块的输出是 criticalPoints (d, 1)。这里，d表示深度，1表示线。参照图11说明临界点识别部901进行的处理的例子。临界点识别部901进行的处理的第1步骤是检索procPulse (d，1)中的极大点及极小点。该步骤的一例是以下的次序。(1)从时间0秒起在心周期内检索第1极大点1000， (2)从对应于第1极大点1000的时间起在心周期内检索第1极小点1001，(3)从第1极小点1001起在心周期内检索第2极大点1002，接着，(4)反复处理上述(1) (3)，直到检索到希望数量的极大点和极小点。心博期间计算部900进行的处理的第2步骤是决定博动的方向。博动的方向可能有1个心周期整体的上方向(在朝上的收缩曲线后具有朝下的舒张曲线)、或下方向(在朝下的收缩曲线后具有朝上的舒张曲线)的两种。设想心博动的收缩部分的期间比舒张部分的期间短。另外，所谓收缩曲线，是心脏的收缩期间中的博动的波形，所谓舒张曲线，是心脏的舒张期间中的博动的波形。作为决定博动的方向的方法的一例，可以举出将图11的线ab与线be的梯度比较。点a及点c是两个连续的局部最大点，点b是点a与点c之间的局部最小点。线ab和线be的梯度也可以取预先定义的数量的心周期的平均。如果线ab的梯度比线be的梯度大，则博动的方向是朝下的，点a是该心周期的起点。如果小，则博动的方向是朝上的，点b 是1心周期的起点。作为决定博动的方向的方法的另一例，可以举出将Dl与D2的长度比较。Dl是从第1极大点1000到与第1极大点1000邻接的第1极小点1001的期间。D2是从上述第1 极小点1001到第2极大点1002的期间。Dl及D2也可以取预先定义的数量的心周期的平均。如果D2比Dl大，则博动的方向是朝下的，如果小，则博动的方向是朝上的。临界点识别部901进行的处理的最终步骤是识别各心周期中的临界点。所谓临界点，包括收缩开始点、舒张结束点、收缩峰、收缩中间点、及舒张中间点，但并不仅限定于这些。收缩中间点及舒张中间点是博动的振幅相对于最大博动振幅等于预先定义的比率时的收缩曲线及舒张曲线上的两点。作为该比率的例子，可以举出最高博动振幅的2分之1。在图11中，如果设想博动方向是朝下的，则收缩开始点是第1极大点1000。舒张结束点是下个第2极大点1002。收缩峰是第1极小点1001。收缩中间点1004及舒张中间点1005是与bf/bg为事前设定的值时的振幅f相对应的、收缩曲线及舒张曲线上的两点d及e。形状特征提取部801提取与博动的波形形状关联的有用的特征。参照图12说明该特征。在与博动的波形形状关联的特征中，有L1/L2、L3/L4、博动偏度、博动峰度、舒张曲线be相对于直线be的偏差、及舒张曲线be中的极值的偏差等，但并不仅限定于这些。另夕卜，直线be也可以是事前设定的曲线。L2(d，l)是心周期的舒张部分的期间。即，在朝上心周期的情况下，参照图12(A)，是从极大点1100到舒张结束点1101的期间，在朝下心周期的情况下，参照图12(B)，是从极小点1102到舒张结束点1103的期间。Ll(d，l)是心周期的收缩部分的期间。即，在朝上心周期的情况下，参照图12(A)， 是从收缩开始点1104到极大点1100的期间，在朝下心周期的情况下，参照图12(B)，是从收缩开始点1105到极小点1102的期间。通过癌性组织中的微小血管系统的不规则的构造，肿瘤性的博动的收缩峰延迟、 或收缩期间变长的可能性较高。如果L1/L2更高，则收缩期间变得更长。这意味着对应的博动为肿瘤性博动的概率更高。L3(d，1)被定义为博动的振幅是最高振幅的预先设定的比率时、即从收缩中间点到舒张中间点的期间的振幅。另外，所谓收缩中间点，是在收缩期间中、收缩曲线的振幅相对于最大振幅等于预先定义的比率时的收缩曲线上的点。此外，所谓舒张中间点，是在舒张期间中、舒张曲线的振幅相对于最大振幅等于预先定义的比率时的舒张曲线上的点。L4(d，1)是1脉冲的期间。比率L3/L4用于测量脉冲的尖锐度。肿瘤性的博动比正常的博动幅度要宽、即具有更扩张的形状的可能性较高。如果L3/L4更高，则意味着博动形状更加扩张。这意味着对应的博动是肿瘤性博动的概率更高。偏度是心博动的非对称性的尺度，表示为pulseSkewness (d，1)。作为计算 pulseSkewness(d, 1)的方法的一例，可以举出使用以下的数式1的方法。[数式1]pulseSkeMmess{d,l) = ！Y(ProcPulse、d，(数式丄)
NfTi σ这里，N是心周期中的帧数，μ及ο是某个心周期内的博动振幅的平均及标准偏差。与L1/L2同样，如果收缩期间更短、舒张期间更长，则该博动是正常的博动的可能性较高。如果正的偏度更大，则表示心周期的舒张期间比心周期的收缩期间长，表示是肿瘤性博动的概率更低。正的峰度是收缩峰相对于正态分布的尖锐度的尺度，表示为pulSeKurt0SiS(d， 1)。作为计算pulseKurtosisW，1)的方法的一例，可以举出使用以下的数式2的方法。[数式2]pulseKurtosis{dJ)=丄ff procPulse{d,lJ)-μ^ —3 (数式幻
L^/^Vσ)」与L3/L4同样，如果博动更窄、即更尖锐，则是正常的博动的可能性更高。如果正的峰度更大，则表示收缩峰周边的心博动的部分更尖锐，暗示是肿瘤性博动的可能性更低。参照图12，作为舒张曲线be相对于直线be的偏差，计算舒张曲线be的振幅间的差异、与直线be的振幅间的差异的和。这考虑称作pulseCurveDeviationW，1)的心博方向，通过以下的数式3计算(式中的“if cardiac direction is upward”意为若搏动方向为朝上，“if cardiac direction is downward”表示若搏动方向为朝下)。[数式3]
{procPidse(dJ, /) - StraightLinebc (d, /, /)) if cardiac direction is upwardpulseCm veDe iation{d,l) ^(d, I, /) - procFulse{d, /, /)) if cardiac direction is downward
,be
(数式3)将其结果用于测量曲线以何种程度弯曲。如果？111%011^劝扑1站1011((1，1)是正，则意味着博动曲线的幅度更宽，如果是负则意味着相反。根据该计算，如果博动曲线 1106具有正的pulseCurveDeviation(d，1)，则可知幅度更宽，如果博动1107具有负的 pulseCurveDeviation(d, 1)，则可知幅度不怎么宽。越是幅度更宽的博动，则是肿瘤性博动的概率越高。作为舒张曲线be中的极值的偏差，计算舒张曲线be内的极大值与极小值之间的博动振幅的差异，称作pulseExtremaDeviationW，1)。在舒张曲线be中没有极值的情况下，极值的偏差是零。一般而言，通过癌性组织中的微小血管系的血管扩张状态，在博动中不再存在重复隆起。因而，如果pulseExtremaDeviation (d，1)更高，则表示存在重复隆起的可能性更高、是肿瘤性博动的概率更低。在图13中表示肿瘤性博动的特征与正常组织的博动的特征的比较的例子。如图13(A)所示，更高的L1/L2、更高的L3/L4、更低的偏度、更低的峰度、正的博动曲线的偏差、以及舒张曲线中的更低的极值的偏差意味着波形的宽度较大(曲线1200)，表示是肿瘤性博动波形的概率更高。另一方面，如图13 (B)所示，更低的L1/L2、更低的L3/L4、 更高的偏度、更高的峰度、负的博动曲线的偏差、以及舒张曲线中的更高的极值的偏差意味着波形的宽度不宽(曲线1201)，是肿瘤性博动波形的概率更低。在形状特征提取部801的另一实施方式中，也可以将某个特征量(例如，Li/ L2、L3/L4、pulseSkewness (d, 1)、pulseKurtosis (d, 1)、pulseCurveDeviation (d, 1)、 pulseExtremaDeviation (d,l)中的至少1个)的值用预先定义的值(例如，从知道是正常组织的区域取得的值)归一化(normalization)，也可以归一化为预先定义的范围(例如、 W，l])。在形状特征提取部801的再另一实施方式中，也可以使得如果任意的特征量的值或将其归一化后的值，值更高时则是肿瘤性博动的概率更高，此外也可以相反使得如果值更低则是肿瘤性博动的概率更低。也可以更具体地，取倒数。这与pulseSkewness (d，1)、 pulseKurtosis (d, 1)、及 pulseExtremaDeviation (d，1)符合。将原始的博动关联特征的值(即，归一化之前的博动关联特征的值)、归一化后的博动关联特征的值、及与博动关联特征的倒数对应的值在以后的段落中也统称作“特征的值”。特征的值，用于关于肿瘤性博动的概率的高低将组织的扫描对象区域分段。该分段的实施方式是对各特征设定阈值、和基于这些阈值根据肿瘤性博动的概率的高低将区域分类，但并不限定于这些。作为阈值选择的一例，可以举出使用特征的值的平均值。作为形状特征提取部801的输出pulseSiapeFeatures，可以选择特征的值、或基于特征的值的聚类结果。再次参照图1，博动关联特征提取部103还可以计算作为关于组织博动的波形的时间轴的特征量的博动时间特征。如果博动的振幅及形状不同，则癌性肿瘤的博动波形整体与正常的组织不同。进而，如果因微小血管系构造而对于血流的阻力不同，则博动的到达时间对癌性组织和正常组织来说是不同的。该博动时间特征例如是心周期的延迟、心波形的差异、及心波形的自回归系数中的至少1个。将它们统称作pulseTemporalFeatures。图20表示计算心周期的延迟量的博动关联特征提取部103具备的心周期延迟计算部731的功能块的一例。如图20所示，心周期延迟计算部731具有基准心周期决定部 1900和延迟计算部1901。向心周期延迟计算部731的输入是博动、即pr0CPulSe(d，1，f)。并且，输出是相对于基准博动(或基准心周期)、即CarDelay(d，l)的(由心周期构成的)各博动的延迟。 所谓心周期的延迟，意味着对于血流的血管阻力。基准心周期决定部1900检测进行与全部博动的比较的基准的博动、或心周期的基准点。作为这样的基准的一例，可以举出心电图检查(ECG)信号。它表现心周期的正确的起点和收缩峰。在不能利用ECG信号的情况下，作为这样的基准的另一例，有扫描的组织中的最强的博动。基准心周期决定部1900的输出是refCarCycle(f)。延迟计算部1901计算相对于作为基准的refCarCycle(f)的、procPulse (d，1，f) 中的心周期的起点或收缩峰的延迟。另外，在计算前，也可以进行用来将心周期的长度标准化的再采样，以便能够与refCarCycle(f)比较。在图21中表示进行心波形的差异的计算的延迟计算部1901的功能块的一例。向延迟计算部1901的输入是博动procPulse (d, 1，f)和carDelay(d, 1)，输出是博动中的心周期的差值、即carDiff(d，1)。通过这样的值，决定心周期以何种程度相互类似或不同。如图21所示，延迟计算部1901具有个别差异计算部2001和整体差异计算部 2002。另外，基准心波形计算部2000计算与所取得的全部的心波形进行比较的基准心波形。在一实施方式中，将这样的基准心波形作为所取得的全部的心波形的平均来计算。 在另一实施方式中，将这样的基准心波形作为所取得的全部的心波形中的最强的心波形计算。在再另一实施方式中，也可以将这样的基准心波形不考虑所取得的心波形而定义为例如从文献引用的周知的心波形等。该块的输出是refCarWaveform(f)。按照各博动procPulse (d，1，f)的每个心波形，个别差异计算部2001计算心波形与基准心波形refCarWaveform(f)的差异。在进行该计算之前，也可以使这些心波形逐渐变小到相同的长度。作为由个别差异计算部2001计算的差异的一例，可以举出各博动 procPulse (d，1，f)与refCarWaveform(f)的差的绝对值的和(或积分值)。作为这样的差异的另一例，可以举出均方根的差。整体差异计算部2002根据由个别差异计算部2001计算出的全部的个别差异的值，按每个博动决定1个整体差异的值。作为这样的整体差异的值的一例，可以举出使用个别差异的值的标准偏差。图22示出用于计算表示心博动的波形的自回归系数的、博动关联特征提取部103 所具备的自回归系数计算部732的功能块的一例。如图22所示，自回归系数计算部732具有博动再采样部2100和自回归式运算部2101。一般而言，可以使用自回归系数将任何波形都记述为自回归过程。即，现在的样本是过去的样本的线性结合。因而，也可以使用这些自回归系数来记述博动。向自回归系数计算部732的输入是博动procPulse (d, 1，f)，输出是模型的系数、即arCoeffs (d, 1)。博动再采样部2100对博动进行再采样，以使各心周期包含标准数量的样本。例如，在每1秒40帧的帧速率(即，时间域的采样速率是40Hz)、每1秒1博动的心博数的情况下，在1个心周期中包含40个样本。但是，例如在超声波换能器的结构不同、并且生理条件不同的情况下，该数量有可能变动。在成为与40个样本不同的数量那样的情况下，对博动进行再采样，以使1个心周期中包含40样本。该块的输出是resampPulse(d，1，f)。自回归式运算部2101求出自回归式的解(S卩，自回归系数)。将其用到博动中。一般而言，自回归式如以下的数式4那样记述。
P[数式

1.一种组织恶性肿瘤检测方法，通过用超声波扫描组织而得到的扫描信号检测包含在上述组织中的恶性肿瘤，其特征在于，包括块划分步骤，将对上述组织进行了扫描的区域划分为多个块； 组织博动推测步骤，按上述多个块中的每个块，基于上述扫描信号推测组织博动，该组织博动是通过上述组织进行博动而发生的上述组织的变位的时间变化；博动关联特征提取步骤，按上述多个块中的每个块，从上述组织博动中提取作为关于上述组织的博动的参数的多个博动关联特征；分布特性计算步骤，按上述多个块中的每个块，计算上述多个博动关联特征的分布特性；以及恶性分类步骤，基于上述分布特性，对上述多个块中的各个块是否是作为包含恶性肿瘤的块的恶性块进行分类。
2.如权利要求1所述的组织恶性肿瘤检测方法，其特征在于，还包括根据在上述恶性分类步骤中被分类为是上述恶性块的块来确定癌性肿瘤的位置的肿瘤位置确定步骤。
3.如权利要求1所述的