CN113164160B - 超声造影成像方法、装置和存储介质 - Google Patents

超声造影成像方法、装置和存储介质 Download PDF

Info

Publication number
CN113164160B
CN113164160B CN202080006417.2A CN202080006417A CN113164160B CN 113164160 B CN113164160 B CN 113164160B CN 202080006417 A CN202080006417 A CN 202080006417A CN 113164160 B CN113164160 B CN 113164160B
Authority
CN
China
Prior art keywords
data
contrast
frame
difference
calculating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202080006417.2A
Other languages
English (en)
Other versions
CN113164160A (zh
Inventor
章希睿
桑茂栋
朱磊
何绪金
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Mindray Bio Medical Electronics Co Ltd
Original Assignee
Shenzhen Mindray Bio Medical Electronics Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen Mindray Bio Medical Electronics Co Ltd filed Critical Shenzhen Mindray Bio Medical Electronics Co Ltd
Publication of CN113164160A publication Critical patent/CN113164160A/zh
Application granted granted Critical
Publication of CN113164160B publication Critical patent/CN113164160B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/481—Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Hematology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

一种超声造影成像方法、装置和存储介质，该方法包括：获取至少两帧造影数据(S310)；计算相邻两帧造影数据之间的信息差异，并计算所述信息差异与所述相邻两帧造影数据中的一帧造影数据之间的关系(S320)；基于所述关系将所述一帧造影数据进行分类，并基于所述分类的结果对所述一帧造影数据的至少部分数据进行抑制处理(S330)。该方法从时域差分(帧间差分)的角度挖掘组织残留与造影剂信号的差异性，能够对造影数据中的组织残留部分与造影剂信号部分进行有效区分，从而能够有效抑制组织残留部分且尽可能不影响造影剂强度的效果。

Description

超声造影成像方法、装置和存储介质
技术领域
本申请涉及超声造影成像技术领域，更具体地涉及一种超声造影成像方法、装置和存储介质。
背景技术
在理想条件下，未注射造影剂的超声造影图像只存在背景噪声。但实际情况中，提取造影剂信号的同时难免引入组织成分，主要由系统发射电路的非对称性导致。以正反谐波方法为例说明，理论上两个幅度、频率一致，相位相反的波形叠加后可完全抵消，然而实际发射电路输出的正、负脉冲并非严格对称，二者叠加后会产生残留信号，体现在造影图像中就是组织残留。
此外，还有三个因素可导致组织残留。一是组织自身也会产生非线性成分，与造影剂的非线性回波一同被检测到；二是由于信号饱和，非线性成像方法中各个脉冲的幅度经过调制后，组织的线性成分无法被抵消，多出现在血管壁、脏器包膜和骨头等强反射面；三是扫查过程中的运动(探头位移、病人呼吸等)，会导致非线性成像方法中各脉冲回波的失相关，进而不满足组织信号相消的条件。
组织残留会干扰用户对造影剂信号的观察，进而影响诊断结果的准确性。对此，现有解决方案一方面从前端硬件入手，研究如何进一步完善发射波形的对称性，但改进难度极大；另一方面则通过新型成像方法(次谐波等)或信号处理手段(滤波等)尽可能抑制组织残留，但除了组织自身产生的非线性成分以外，其他原因导致的组织残留与造影剂信号在频域无法区分，抑制能力受限。
发明内容
本申请提供一种超声造影成像方案，能够有效对造影数据中的组织残留进行抑制。下面简要描述本申请提出的超声造影成像方案，更多细节将在后续结合附图在具体实施方式中加以描述。
本申请一方面，提供了一种超声造影成像方法，所述方法包括：控制超声探头向目标对象发射超声波，接收所述超声波的回波，并基于所述超声波的回波获取超声回波信号；基于所述超声回波信号获取至少两帧造影信号数据，其中，每帧造影信号数据用于生成一帧造影图像；计算相邻两帧造影信号数据之间的信息差异，并计算所述信息差异与所述相邻两帧造影信号数据中的一帧造影信号数据之间的关系；基于所述关系将所述一帧造影信号数据进行区域划分，并对不同区域的造影信号数据采用不同的抑制因子进行处理；以及基于经处理的所述一帧造影信号数据生成造影图像。
本申请另一方面，提供了一种超声造影成像方法，所述方法包括：控制超声探头向目标对象发射超声波，接收所述超声波的回波，并基于所述超声波的回波获取超声回波信号；基于所述超声回波信号获取造影信号数据，并基于所述造影信号数据生成造影图像；获取至少两帧造影图像，计算相邻两帧造影图像之间的信息差异，并计算所述信息差异与所述相邻两帧造影图像中的一帧造影图像之间的关系；以及基于所述关系将所述一帧造影图像进行区域划分，并对所述一帧造影图像的不同区域采用不同的抑制因子进行处理。
本申请再一方面，提供了一种超声造影成像装置，所述装置包括：超声探头、发射/接收序列控制器和处理器，其中，所述发射/接收序列控制器用于：激励所述超声探头向目标对象发射超声波，接收所述超声波的回波，基于所述超声波的回波获取超声回波信号；所述处理器用于执行上述超声造影成像方法。
本申请又一方面，提供了一种超声造影成像方法，所述方法包括：获取至少两帧造影数据；计算相邻两帧造影数据之间的信息差异，并计算所述信息差异与所述相邻两帧造影数据中的一帧造影数据之间的关系；基于所述关系将所述一帧造影数据进行分类，并基于所述分类的结果对所述一帧造影数据的至少部分数据进行抑制处理。
本申请又一方面，提供了一种超声造影成像方法，所述方法包括：控制超声探头向目标对象发射超声波，接收所述超声波的回波，并基于所述超声波的回波获取超声回波信号；基于所述超声回波信号获取至少两帧造影信号数据，其中，每帧造影信号数据用于生成一帧造影图像；计算相邻两帧造影信号数据之间的信息差异，并计算所述信息差异与所述相邻两帧造影信号数据中的一帧造影信号数据之间的关系；基于所述关系将所述一帧造影信号数据进行分类，并基于所述分类的结果对所述一帧造影信号数据的至少部分数据进行抑制处理；以及基于经处理的所述一帧造影信号数据生成造影图像。
本申请又一方面，提供了一种超声造影成像方法，所述方法包括：控制超声探头向目标对象发射超声波，接收所述超声波的回波，并基于所述超声波的回波获取超声回波信号；基于所述超声回波信号获取造影信号数据，并基于所述造影信号数据生成造影图像；计算相邻两帧造影图像之间的信息差异，并计算所述信息差异与所述相邻两帧造影图像中的一帧造影图像之间的关系；以及基于所述关系将所述一帧造影图像进行分类，并基于所述分类的结果对至少部分的所述一帧造影图像进行抑制处理。
本申请再一方面，提供了一种超声造影成像装置，所述装置包括超声探头、发射/接收序列控制器和处理器，其中：所述发射/接收序列控制器用于：激励所述超声探头向目标对象发射超声波，接收所述超声波的回波，基于所述超声波的回波获取超声回波信号；所述处理器用于：基于所述超声回波信号获取至少两帧造影数据；计算相邻两帧造影数据之间的信息差异，并计算所述信息差异与所述相邻两帧造影数据中的一帧造影数据之间的关系；基于所述关系将所述一帧造影数据进行分类，并基于所述分类的结果对所述一帧造影数据的至少部分数据进行抑制处理。
本申请再一方面，提供了一种超声造影成像装置，所述装置包括超声探头、发射/接收序列控制器和处理器，其中：所述发射/接收序列控制器用于：激励所述超声探头向目标对象发射超声波，接收所述超声波的回波，基于所述超声波的回波获取超声回波信号；所述处理器用于执行上述超声造影成像方法。
本申请又一方面，提供了一种存储介质，所述存储介质上存储有计算机程序，所述计算机程序在运行时执行上述超声造影成像方法。
根据本申请实施例的超声造影成像方法、装置和存储介质从时域差分(帧间差分)的角度挖掘组织残留与造影剂信号的差异性，能够对造影数据中的组织残留部分与造影剂信号部分进行有效区分，从而能够有效抑制组织残留部分且尽可能不影响造影剂强度的效果。
附图说明
图1示出造影图像中的组织残留的示意图。
图2示出根据本申请一个实施例的超声造影成像装置的示意性框图。
图3示出根据本申请一个实施例的超声造影成像方法的示意性流程图。
图4示出未注射造影剂时的人体甲状腺造影双实时图像。
图5A示出一帧造影数据中的第一类别数据点在帧间差分前后信号幅度的对比结果。
图5B示出一帧造影数据中的第二类别数据点在帧间差分前后信号幅度的对比结果。
图5C示出一帧造影数据中的第三类别数据点在帧间差分前后信号幅度的对比结果。
图6A示出未注射造影剂时且未采用根据本申请实施例的超声造影成像方法处理的人体甲状腺造影双实时图像的示例。
图6B示出未注射造影剂时且采用根据本申请实施例的超声造影成像方法处理后的人体甲状腺造影双实时图像的示例。
图7A示出注射造影剂后且未采用根据本申请实施例的超声造影成像方法处理的犬肝脏造影双实时图像的示例。
图7B示出注射造影剂后且采用根据本申请实施例的超声造影成像方法处理后的犬肝脏造影双实时图像的示例。
图8示出根据本申请另一个实施例的超声造影成像方法的示意性流程图。
图9示出根据本申请再一个实施例的超声造影成像方法的示意性流程图。
图10示出根据本申请又一个实施例的超声造影成像方法的示意性流程图。
图11示出根据本申请再一个实施例的超声造影成像方法的示意性流程图。
图12示出根据本申请另一个实施例的超声造影成像装置的示意性框图。
具体实施方式
为了使得本申请的目的、技术方案和优点更为明显，下面将参照附图详细描述根据本申请的示例实施例。显然，所描述的实施例仅仅是本申请的一部分实施例，而不是本申请的全部实施例，应理解，本申请不受这里描述的示例实施例的限制。基于本申请中描述的本申请实施例，本领域技术人员在没有付出创造性劳动的情况下所得到的所有其它实施例都应落入本申请的保护范围之内。
在下文的描述中，给出了大量具体的细节以便提供对本申请更为彻底的理解。然而，对于本领域技术人员而言显而易见的是，本申请可以无需一个或多个这些细节而得以实施。在其他的例子中，为了避免与本申请发生混淆，对于本领域公知的一些技术特征未进行描述。
应当理解的是，本申请能够以不同形式实施，而不应当解释为局限于这里提出的实施例。相反地，提供这些实施例将使公开彻底和完全，并且将本申请的范围完全地传递给本领域技术人员。
在此使用的术语的目的仅在于描述具体实施例并且不作为本申请的限制。在此使用时，单数形式的“一”、“一个”和“所述/该”也意图包括复数形式，除非上下文清楚指出另外的方式。还应明白术语“组成”和/或“包括”，当在该说明书中使用时，确定所述特征、整数、步骤、操作、元件和/或部件的存在，但不排除一个或更多其它的特征、整数、步骤、操作、元件、部件和/或组的存在或添加。在此使用时，术语“和/或”包括相关所列项目的任何及所有组合。
为了彻底理解本申请，将在下列的描述中提出详细的步骤以及详细的结构，以便阐释本申请提出的技术方案。本申请的较佳实施例详细描述如下，然而除了这些详细描述外，本申请还可以具有其他实施方式。
在理想条件下，未注射造影剂的超声造影图像只存在背景噪声。但实际情况中，提取造影剂信号的同时难免引入组织成分，体现在造影图像中就是组织残留，如图1所示的。组织残留会干扰用户对造影剂信号的观察，进而影响图像分析结果的准确性。
对此，本申请提供一种超声造影成像方案，其从时域差分(帧间差分)的角度挖掘组织残留与造影剂信号的差异性，能够有效对造影数据中的组织残留进行抑制。下面参照附图2到附图12对本申请的超声造影成像方案进行详细描述。
图2为本申请实施例中的超声造影成像装置的示意图。该超声造影成像装置100可以包括探头110、发射电路112、接收电路114、波束合成电路116、处理器118和显示器120。其中，发射电路112和接收电路114可以通过发射/接收选择开关122与超声探头110连接。发射电路112可以激励探头110向含造影剂的目标组织发射超声波；接收电路114可以通过探头110接收从目标组织返回的超声回波，从而获得超声回波信号/数据；该超声回波信号/数据经过波束合成电路116进行波束合成处理后，送入处理器118。处理器118对该超声回波信号/数据进行处理，以获得目标组织的超声图像以及造影图像。处理器118获得的超声图像和造影图像可以存储于存储器124中。这些超声图像和造影图像可以在显示器120上显示。更详细的描述可以参见本说明书的后续实施例。
应理解，图2所示的超声造影成像装置100所包括的部件只是示意性的，其可以包括更多或更少的部件。本发明对此不限定。图2所示的超声造影成像装置100可以用于执行下文中的根据本发明实施例的超声造影成像方法。
图3示出根据本申请一个实施例的超声造影成像方法300。如图3所示，超声造影成像方法300可以包括如下步骤：
在步骤S310，获取至少两帧造影数据。
在步骤S320，计算相邻两帧造影数据之间的信息差异，并计算信息差异与相邻两帧造影数据中的一帧造影数据之间的关系。
在步骤S330，基于关系将一帧造影数据进行分类，并基于分类的结果对一帧造影数据的至少部分数据进行抑制处理。
在本申请的实施例中，在步骤S310所获取的至少两帧造影数据可以是至少两帧造影信号数据，其中每帧造影信号数据用于生成一帧造影图像；此外，在步骤S310所获取的至少两帧造影数据也可以是至少两帧造影图像数据。此外，在步骤S310获取的至少两帧造影数据可以是实时产生的，也可以是来自任何源的。
在本申请的实施例中，基于相邻两帧造影数据之间的信息差异与相邻两帧造影数据中的一帧造影数据之间的关系，能够将一帧造影数据中组织残留强的数据点和组织残留弱(甚至无组织残留)的数据点明显地区分开来(下文中将通过示例进一步详细描述)，使得能够基于该关系对一帧造影数据进行分类(例如分类为组织残留强的数据点、组织残留弱的数据点、无组织残留的数据点等)，并基于分类的结果对一帧造影数据的至少部分数据(例如组织残留强的数据点和组织残留弱的数据点)进行抑制处理，从而达到能够有效抑制组织残留部分且尽可能不影响造影剂强度的效果。
在本申请的实施例中，相邻两帧造影数据包括前一帧造影数据(例如称为第n帧造影数据Dn，其中n为自然数)和后一帧造影数据(例如称为第n+1帧造影数据Dn+1。基于此，在一个实施例中，步骤S320中的计算相邻两帧造影数据之间的信息差异，可以包括：计算后一帧造影数据与前一帧造影数据的前向差分结果，即Dn+1-Dn的结果。相应地，步骤S320中的计算信息差异与相邻两帧造影数据中的一帧造影数据之间的关系，可以包括：计算前向差分结果(即Dn+1-Dn)与前一帧造影数据(即Dn)的数据关系，即计算(Dn+1-Dn)与Dn之间的数据关系。在另一个实施例中，步骤S320中的计算相邻两帧造影数据之间的信息差异，可以包括：计算前一帧造影数据与后一帧造影数据的后向差分结果，即Dn-Dn+1的结果。相应地，步骤S320中的计算信息差异与相邻两帧造影数据中的一帧造影数据之间的关系，可以包括：计算后向差分结果(即Dn-Dn+1)与后一帧造影数据(即Dn+1)的数据关系，即计算(Dn-Dn+1)与Dn+1之间的数据关系。因此，步骤S320中计算信息差异与相邻两帧造影数据中的一帧造影数据之间的关系，也可以描述为：计算信息差异与相邻两帧造影数据中被差分的一帧造影数据之间的关系。
在本申请的实施例中，当造影数据为诸如波束合成后的造影信号数据时，一帧造影数据(无论是前述的Dn还是Dn+1)可以为包含幅度和相位等信息的复解析信号，例如为I×J维复矩阵，其中I和J可以分别表示纵向点数和横向线数。基于此，前述的相邻两帧造影数据之间的信息差异可以包括幅度差异、相位差异或频率差异。当造影数据为诸如模数转换后的造影信号数据时，一帧造影数据(无论是前述的Dn还是Dn+1)可以为至少包含幅度信息的解析信号，例如为I×J维矩阵，其中I和J可以分别表示纵向点数和横向线数。基于此，前述的相邻两帧造影数据之间的信息差异可以包括幅度差异。
相应地，当相邻两帧造影数据之间的信息差异为幅度差异时，在一个实施例中，步骤S320中的计算相邻两帧造影数据之间的信息差异，并计算信息差异与相邻两帧造影数据中的一帧造影数据之间的关系，可以包括：计算相邻两帧造影数据之间的差分数据，并计算差分数据和一帧造影数据各自的幅度数据；基于差分数据的幅度数据和一帧造影数据的幅度数据计算差分数据与一帧造影数据的幅度比值，得到差分数据与一帧造影数据之间的关系结果。在该实施例中，先计算差分数据与被差分的一帧造影数据各自的幅度数据(例如通过包络检测的方式)，再计算这两者的比值，并将计算得到的幅度比值作为差分数据与被差分的一帧造影数据之间的关系结果。在本申请的进一步的实施例中，可以在对幅度数据进行对数压缩后再求比值，以减小计算量。例如，比值(关系结果)表示为R，则R可以表示为下式：
R＝20log10(abs(Dn+1-Dn))./20log10(abs(Dn))
其中，abs表示求幅值(绝对值)，“log10”表示以10为底的对数运算，符号“./”表示矩阵的点对点除法运算，R为I×J维实矩阵，I和J如前文的。应理解，该式仅是示例性的。
当相邻两帧造影数据之间的信息差异为幅度差异时，在另一个实施例中，步骤S320中的计算相邻两帧造影数据之间的信息差异，并计算信息差异与相邻两帧造影数据中的一帧造影数据之间的关系，可以包括：计算相邻两帧造影数据之间的差分数据，对差分数据与一帧造影信号数据求比值；对比值进行幅度求解得到差分数据与一帧造影数据的幅度比值，作为差分数据与一帧造影数据之间的关系结果。在该实施例中，先计算差分数据与被差分的一帧造影数据之间的比值，再对该比值进行幅度求解，最终也可以得到差分数据与被差分的一帧造影数据之间的幅度比值，以作为差分数据与被差分的一帧造影数据之间的关系结果。应该理解，除了对差分数据与被差分的一帧造影信号数据求比值以外，在另外的实施例中，也可以计算差分数据与被差分的一帧造影信号数据的差值，再基于差值与预设阈值的关系确定差分数据与被差分的一帧造影数据之间的关系结果。
在本申请的实施例中，当相邻两帧造影数据之间的信息差异为相位差异时，在一个实施例中，步骤S320中的计算相邻两帧造影数据之间的信息差异，并计算信息差异与相邻两帧造影数据中的一帧造影数据之间的关系，可以包括：计算相邻两帧造影数据之间的差分数据，并计算差分数据和一帧造影数据各自的相位数据；基于差分数据的相位数据和一帧造影数据的相位数据计算差分数据与一帧造影数据的相位差值，得到差分数据与一帧造影数据之间的关系结果。在该实施例中，先计算差分数据与被差分的一帧造影数据各自的相位数据，再计算这两者的差值，并将计算得到的相位差值作为差分数据与被差分的一帧造影数据之间的关系结果。
在本申请的实施例中，当相邻两帧造影数据之间的信息差异为相位差异时，在另一个实施例中，步骤S320中的计算相邻两帧造影数据之间的信息差异，并计算信息差异与相邻两帧造影数据中的一帧造影数据之间的关系，可以包括：计算相邻两帧造影数据之间的差分数据，对差分数据与一帧造影数据求比值；对比值进行相位求解，得到差分数据与一帧造影数据的相位差值，作为差分数据与一帧造影数据之间的关系结果。在该实施例中，先计算差分数据与被差分的一帧造影数据之间的比值，再对该比值进行相位求解，最终也可以得到差分数据与被差分的一帧造影数据之间的相位差值，以作为差分数据与被差分的一帧造影数据之间的关系结果。
在本申请的实施例中，当造影数据为造影图像数据时，一帧造影数据(无论是前述的Dn还是Dn+1)可以为包含灰度信息等的数据集合，例如为V×H维矩阵，其中V和H可以分别表示垂直像素点数和水平像素点数。基于此，在一个实施例中，步骤S320中的计算相邻两帧造影数据之间的信息差异，并计算信息差异与相邻两帧造影数据中的一帧造影数据之间的关系，可以包括：计算后一帧造影图像与前一帧造影图像的前向差分结果；计算前向差分结果与前一帧造影图像的灰度比值。在另一个实施例中，步骤S320中的计算相邻两帧造影数据之间的信息差异，并计算信息差异与相邻两帧造影数据中的一帧造影数据之间的关系，可以包括：计算前一帧造影图像与后一帧造影图像的后向差分结果；计算后向差分结果与后一帧造影图像的灰度比值。
以上示例性地描述了关于相邻两帧造影数据之间的信息差异与相邻两帧造影数据中的一帧造影数据之间的关系的计算方式的示例，如前所述的，基于相邻两帧造影数据之间的信息差异与相邻两帧造影数据中的被差分的一帧造影数据之间的关系，能够将一帧造影数据中组织残留强的数据点和组织残留弱(甚至无组织残留)的数据点明显地区分开来。下面结合图4和图5A到图5C来说明。
图4示出了未注射造影剂时的人体甲状腺造影双实时图像。如图4所示的，示出的是在未注射造影剂的条件下，共采集100帧人体甲状腺的造影波束合成后I/Q数据中的其中1帧对应的造影双实时图像，其中图4的左半部分是组织图像，右半部分是造影图像。如图4的右半部分所示的，组织残留主要集中在框410和框420中所包括的数据点处，其余数据点处只有噪声。
基于此，可以将一帧造影数据(造影信号数据或造影图像数据)进行分类，例如框410中所包括的数据点(例如称为第一类别数据点)处主要表示为易导致信号饱和的皮质层，其为组织残留最明显的数据点位置区域，一般为近场区，与探头的相对位移较小，且远离待扫查器官，不易受病人呼吸的影响；帧间差分后，该数据点位置区域的信号幅度会明显降低，如图5A所示的。框420中所包括的数据点(例如称为第二类别数据点)处主要包括血管壁、器官包膜等强反射面导致的组织残留，一般为中场区，位于待扫查器官区域，较易受到病人呼吸的影响；帧间差分后，该数据点位置区域信号幅度的下降不如框410中的数据点位置区域中的明显，如图5B所示的。框430中所包括的数据点(例如称为第三类别数据点)处主要包括帧间噪声一般为高斯-白噪声，一般为远场区；帧间差分后，整体会呈现方差变大的趋势，如图5C所示的。
由此可见，一帧造影数据中不同类别的数据点在帧间差分前后所显示的特性有明显的不同，因此本申请的实施例以此为依据，能够将一帧造影数据中组织残留强的数据点和组织残留弱(甚至无组织残留)的数据点明显地区分开来，使得能够基于该关系对一帧造影数据进行分类，并基于分类的结果对一帧造影数据的至少部分数据进行抑制处理，从而能够有效抑制组织残留部分且尽可能不影响造影剂强度的效果。
如前所述的，一帧造影数据可以表示为一个矩阵，因此相邻两帧造影数据的差分数据也可以表示为一个矩阵，为了区分，可以将相邻两帧造影数据的差分数据表示为第一矩阵，并将要与该第一矩阵比较关系的一帧造影数据表示为第二矩阵，则步骤S330中的基于关系将一帧造影数据进行分类，可以包括：确定第一矩阵中的数据点与第二矩阵中的相应位置处的数据点的关系结果所属的数值区间；基于关系结果所属的数值区间将第二矩阵中的相应位置处的数据点划分为属于相应的类别。在该实施例中，可以通过预设若干个数值区间的方式来划分第二矩阵(即被差分的一帧造影数据)中属于不同数值区间的不同类别的数据点。
在一个实施例中，数值区间可以是基于一个或更多个预设阈值来确定的。很明显，预设阈值的数目越多，所划分的数值区间越多，对一帧造影数据中数据点的分类越精细。为了描述简单，下面以基于一个预设阈值T(其中T取值范围例如为0到1)的数值区间为例来描述。在该示例中，预设的数值区间可以包括第一数值区间和第二数值区间，第一数值区间例如为(0,T]，第二数值区间例如为(T,1]。第一数值区间(0,T]对应第一类别数据点，第二数值区间(T,1]对应第二类别数据点。当确定第一矩阵(相邻两帧造影数据的差分数据)中的数据点与第二矩阵(相邻两帧造影数据中被差分的一帧造影数据)中的相应位置处的数据点的关系结果(例如前述的幅度比值R)在第一数值区间内，即R∈(0,T]时，将第二矩阵中的相应位置处的数据点划分为第一类别数据点；当确定第一矩阵中的数据点与第二矩阵中的相应位置处的数据点的关系结果(例如前述的幅度比值R)在第二数值区间内，即R∈(T,1]时，将第二矩阵中的相应位置处的数据点划分为第二类别数据点。
如前文结合图4和图5A的，帧间差分后，信号幅度与被差分的一帧造影数据的信号幅度相比会明显降低，即差分数据的信号幅度与被差分的一帧造影数据的信号幅度之间的比值相对较小，即R值相对较小；如前文结合图4和图5B的，帧间差分后，信号幅度与被差分的一帧造影数据的信号幅度相比略弱或接近，即差分数据的信号幅度与被差分的一帧造影数据的信号幅度之间的比值相对较大，即R值相对较大。因此，对应于R∈(0,T]的第一类别数据点与对应于R∈(T,1]的第二类别数据点相比，第一类别数据点处的组织残留强于(多于)第二类别数据点处的组织残留。
在本申请的进一步的实施例中，基于一个预设阈值T的数值区间还可以包括第三数值区间，即第一数值区间和第二数值区间以外的数值区间(1,+∞]∪(-∞,0)。当确定第一矩阵(相邻两帧造影数据的差分数据)中的数据点与第二矩阵(相邻两帧造影数据中被差分的一帧造影数据)中的相应位置处的数据点的关系结果(例如前述的幅度比值R)在第三数值区间内，即R∈(1,+∞]∪(-∞,0)时，将第二矩阵中的相应位置处的数据点划分为第三类别数据点。如前文结合图4和图5C的，帧间差分后，信号幅度与被差分的一帧造影数据的信号幅度相比更强，即R大于1，表示第三类别数据点是噪声数据点(R小于0认为是噪声随机性的体现)。
根据上述的分类结果，可以有针对性地进行抑制处理。在本申请的实施例中，步骤S330中的基于分类的结果对一帧造影数据中的至少部分数据进行抑制处理，可以包括：基于分类的结果确定一帧造影数据中需要进行抑制处理的数据；对于需要进行抑制处理的数据中不同类别的数据采用不同的抑制因子进行处理。
在本申请的实施例中，可以对不同类别的数据进行不同的处理，例如对于前文中描述的组织残留强的数据点(如图4中的框410中的数据点或前文被划分为第一类别数据点的数据点)，可以采用较强的抑制因子进行处理(重点抑制)，以更好地抑制造影数据中的强组织残留，其中抑制因子可以表示为F，示例性地，F＝R2，对于正小数R而言，其求平方后将更接近于0，可增强抑制效果。对于前文中描述的组织残留弱的数据点(如图4中的框420中的数据点或前文被划分为第二类别数据点的数据点)，可以采用较弱的抑制因子进行处理(适度抑制)，以在不影响造影剂强度的前提下适当地抑制造影数据中的弱组织残留，其中抑制因子可以表示为F，示例性地，
Figure BDA0003082909340000131
对于正小数R而言，其开根号后更趋向于1，以求较为保守的抑制效果。对于前文中描述的无组织残留的数据点(如图4中的框430中的数据点或前文被划分为第三类别数据点的数据点)，可以无需进行抑制处理，也可以令抑制因子F等于1，即不作处理。基于此，可以实现组织残留的有效抑制且不影响造影剂强度的效果。应理解，以上框410-430表示的不同类别的数据点仅为一种示例，并不用以限制本发明。
由于被差分的一帧造影数据可以表示为第二矩阵，第二矩阵中的不同类别的数据点可以采用不同的抑制因子F来处理，因此，F也可以表示为一个矩阵，将该矩阵与第二矩阵进行点对点相乘的方式，即可实现对造影数据的组织残留印制处理。如前所述的，F趋向于0，意味着抑制效果的逐渐加强；F趋向于1，则意味着抑制效果的弱化，甚至不作处理。
当然，前述以单阈值划分的数值区间、抑制因子的计算方式以及抑制因子的作用方式仅是示例性的，还可以根据本申请的教导采用多个阈值划分数值区间，抑制因子也可以是其他的计算方式(例如平方和开根号以外的其他计算方式)，抑制因子的作用方式也可以是其他的方式(例如不是点对点相乘，而是相加、相除等方式)。
基于上面的描述，根据本申请实施例的超声造影成像方法300从时域差分(帧间差分)的角度挖掘组织残留与造影剂信号的差异性，能够对造影数据中的组织残留部分与造影剂信号部分进行有效区分，从而能够有效抑制组织残留部分且尽可能不影响造影剂强度的效果。
可以参照图6A到图7B来查看根据本申请实施例的超声造影成像方法所获得的组织残留抑制效果。其中，考虑了两个场景：未注射造影剂时的残留抑制效果；注射造影剂后的残留抑制效果以及对造影剂强度的影响程度。图6A示出了未注射造影剂时且未采用根据本申请实施例的超声造影成像方法处理的人体甲状腺造影双实时图像的示例，图6B示出了未注射造影剂时且采用根据本申请实施例的超声造影成像方法处理后的人体甲状腺造影双实时图像的示例。如图6A和图6B所示的，近场残留被显著抑制，中场残留轻微弱化，背景噪声没有变化，达到了实施例方法所设计的效果。图7A示出了注射造影剂后且未采用根据本申请实施例的超声造影成像方法处理的犬肝脏造影双实时图像的示例。图7B示出了注射造影剂后且采用根据本申请实施例的超声造影成像方法处理后的犬肝脏造影双实时图像的示例。如图7A和图7B所示的，近场残留被显著抑制，造影剂信号轻微损失，背景噪声无变化，同样达到了期望的效果。关于如何进一步减少造影剂信号的损失，可以在本申请的教导下通过优化实施例方法来达成，比如：采用处理效果更为平滑的多阈值、改进抑制因子的计算方法等。
以上示例性地示出了根据本发明实施例的超声造影成像方法300。下面参考图8到图12描述根据本发明其他实施例的超声造影成像方法。
图8示出了根据本申请另一个实施例的超声造影成像方法800的示意性流程图。如图8所示，超声造影成像方法800可以包括如下步骤：
在步骤S810，控制超声探头向目标对象发射超声波，接收超声波的回波，并基于超声波的回波获取超声回波信号。
在步骤S820，基于超声回波信号获取至少两帧造影信号数据，其中，每帧造影信号数据用于生成一帧造影图像。
在步骤S830，计算相邻两帧造影信号数据之间的信息差异，并计算信息差异与相邻两帧造影信号数据中的一帧造影信号数据之间的关系。
在步骤S840，基于关系将一帧造影信号数据进行区域划分，并对不同区域的造影信号数据采用不同的抑制因子进行处理。
在步骤S850，基于经处理的一帧造影信号数据生成造影图像。
在根据本申请实施例的超声造影成像方法800中，通过实时的方式获取至少两帧造影信号数据，并基于相邻两帧造影信号数据之间的信息差异与相邻两帧造影信号数据中的一帧造影信号数据之间的关系，将一帧造影信号数据进行区域划分，不同区域的造影信号数据采用不同的抑制因子进行处理，最后基于经处理的造影信号数据生成造影图像。
因此，与前文的根据本申请实施例的超声造影成像方法300相类似的，在根据本申请实施例的超声造影成像方法800中，也是基于相邻两帧造影数据之间的信息差异与相邻两帧造影数据中的被差分的一帧造影数据之间的关系对被差分的一帧造影数据进行组织残留抑制处理。不同的是，在根据本申请实施例的超声造影成像方法800中，是基于相邻两帧造影信号数据之间的信息差异与被差分的一帧造影信号数据之间的关系将被差分的一帧造影信号数据进行区域划分，并基于不同区域经不同抑制因子处理后得到的造影信号数据生成造影图像。为了简洁，此处主要描述本申请实施例的超声造影成像方法800与本申请实施例的超声造影成像方法300的不同之处，其余相似之处仅简明概括。
在本申请的实施例中，基于相邻两帧造影信号数据之间的信息差异与相邻两帧造影信号数据中的一帧造影信号数据之间的关系，能够将一帧造影信号数据中组织残留强的区域(如前文图4中所示的框410内的区域对应的造影信号数据)和组织残留弱的区域(如前文图4中所示的框420内的区域对应的造影信号数据)和无组织残留的区域(如前文图4中所示的框430内的区域对应的造影信号数据)明显地区分开来，使得能够基于该关系对一帧造影信号数据进行分区(例如分为组织残留强的区域、组织残留弱的区域、无组织残留的区域等)，并基于区域划分的结果对一帧造影信号数据的不同区域采用不同的抑制因子进行抑制处理，从而能够使得生成的造影图像组织残留部分被有效抑制且不影响造影剂强度的效果。
在本申请的实施例中，相邻两帧造影信号数据包括前一帧造影信号数据和后一帧造影信号数据，步骤S830中的计算相邻两帧造影信号数据之间的信息差异，可以包括：计算后一帧造影信号数据与前一帧造影信号数据的前向差分结果；步骤S830中的计算信息差异与相邻两帧造影信号数据中的一帧造影信号数据之间的关系，包括：计算前向差分结果与前一帧造影信号数据的大小关系。
在本申请的实施例中，相邻两帧造影信号数据包括前一帧造影信号数据和后一帧造影信号数据，步骤S830中的计算相邻两帧造影信号数据之间的信息差异，包括：计算前一帧造影信号数据与后一帧造影信号数据的后向差分结果；步骤S830中的计算信息差异与相邻两帧造影信号数据中的一帧造影信号数据之间的关系，包括：计算后向差分结果与后一帧造影信号数据的大小关系。
在本申请的实施例中，步骤S820中的基于超声回波信号获取的造影信号数据可以为对超声回波信号至少进行波束合成处理后得到的造影信号数据，相邻两帧造影信号数据之间的信息差异包括幅度差异、相位差异或频率差异。
在本申请的实施例中，步骤S820中的基于超声回波信号获取的造影信号数据为对超声回波信号至少进行模数转换处理后得到的造影信号数据，相邻两帧造影信号数据之间的信息差异包括幅度差异。
在本申请的实施例中，步骤S830中的计算相邻两帧造影信号数据之间的信息差异，并计算信息差异与相邻两帧造影信号数据中的一帧造影信号数据之间的关系，可以包括：计算相邻两帧造影信号数据之间的差分数据，对差分数据和被差分的一帧造影信号数据分别进行包络检测，得到差分数据和被差分的一帧造影信号数据各自的幅度数据；基于差分数据的幅度数据和被差分的一帧造影信号数据的幅度数据计算差分数据与被差分的一帧造影信号数据的幅度比值，得到差分数据与被差分的一帧造影信号数据之间的关系结果。
在本申请的实施例中，还可以对包络检测之后得到的结果进行对数压缩，基于幅度比值将被差分的一帧造影信号数据进行区域划分时可以是基于经对数压缩的幅度比值。
在本申请的实施例中，步骤S830中的计算相邻两帧造影信号数据之间的信息差异，并计算信息差异与相邻两帧造影信号数据中被差分的的一帧造影信号数据之间的关系，可以包括：计算相邻两帧造影信号数据之间的差分数据，对差分数据与被差分的一帧造影信号数据求比值；对比值进行包络检测，得到差分数据与被差分的一帧造影信号数据的幅度比值，作为差分数据与被差分的一帧造影信号数据之间的关系结果。
在本申请的实施例中，步骤S830中的计算相邻两帧造影信号数据之间的信息差异，并计算信息差异与相邻两帧造影信号数据中被差分的一帧造影信号数据之间的关系，可以包括：计算相邻两帧造影信号数据之间的差分数据，并计算差分数据和被差分的一帧造影信号数据各自的相位数据；基于差分数据的相位数据和被差分的一帧造影信号数据的相位数据计算差分数据与被差分的一帧造影信号数据的相位差值，得到差分数据与被差分的一帧造影信号数据之间的关系结果。
在本申请的实施例中，步骤S830中的计算相邻两帧造影信号数据之间的信息差异，并计算信息差异与相邻两帧造影信号数据中被差分的的一帧造影信号数据之间的关系，可以包括：计算相邻两帧造影信号数据之间的差分数据，对差分数据与被差分的一帧造影信号数据求比值；对比值进行辐角检测，得到差分数据与被差分的一帧造影信号数据的相位差值，作为差分数据与被差分的一帧造影信号数据之间的关系结果。
在本申请的实施例中，相邻两帧造影信号数据的差分数据可以包括第一矩阵，被差分的一帧造影信号数据可以包括第二矩阵，步骤S840中的基于关系将一帧造影信号数据进行区域划分，可以包括：确定第一矩阵中的数据点与第二矩阵中的相应位置处的数据点的关系结果所属的数值区间；基于关系结果所属的数值区间将第二矩阵中的相应位置处的数据点划分为属于相应的区域。在该实施例中，可以通过预设若干个数值区间的方式来划分第二矩阵(即被差分的一帧造影信号数据)中属于不同数值区间的不同区域的数据点。
在一个实施例中，数值区间是基于一个或更多个预设阈值来确定的。很明显，预设阈值的数目越多，所划分的数值区间越多，对一帧造影信号数据的区域划分越精细。为了描述简单，下面以基于一个预设阈值T(其中T取值范围例如为0到1)的数值区间为例来描述。在该示例中，预设的数值区间可以包括第一数值区间和第二数值区间，第一数值区间例如为(0,T]，第二数值区间例如为(T,1]。
当确定第一矩阵(相邻两帧造影信号数据的差分数据)中的任一数据点与第二矩阵(相邻两帧造影数据中被差分的一帧造影信号数据)中的相应位置处的数据点的关系结果(例如前述的幅度比值R)属于第一数值区间，即R∈(0,T]时，将第二矩阵中的相应位置处的数据点划分为第一区域的数据点，例如图示的近场残留区域中的数据点。当第一矩阵中的任一数据点与第二矩阵中的相应位置处的数据点的关系结果(例如前述的幅度比值R)属于第二数值区间，即R∈(T,1]时，将第二矩阵中的相应位置处的数据点划分为第二区域的数据点，例如图示的中场残留区域中的数据点。其中第一区域中的组织残留多于第二区域中的组织残留。
在本申请的进一步的实施例中，基于一个预设阈值T的数值区间还可以包括第三数值区间，即第一数值区间和第二数值区间以外的数值区间(1,+∞]∪(-∞,0)。当确定第一矩阵(相邻两帧造影信号数据的差分数据)中的数据点与第二矩阵(相邻两帧造影信号数据中被差分的一帧造影信号数据)中的相应位置处的数据点的关系结果(例如前述的幅度比值R)在第三数值区间内，即R∈(1,+∞]∪(-∞,0)时，将第二矩阵中的相应位置处的数据点划分为噪声区域中的数据点。
根据上述的区域划分结果，可以有针对性地进行抑制处理。如前文结合图4到图5C所示的，近场残留区域的组织残留最明显，中场残留区域的组织残留较弱，噪声区域基本没有组织残留。因此，在本申请的实施例中，近场残留区域中的数据点采用第一抑制因子进行处理，中场残留区域中的数据点采用第二抑制因子进行处理，第一抑制因子的抑制作用大于第二抑制因子的抑制作用。例如，第一抑制因子可以表示为F1，示例性地，F1＝R2，对于正小数R而言，其求平方后将更接近于0，可增强抑制效果。第二抑制因子可以表示为F2，示例性地，
Figure BDA0003082909340000191
对于正小数R而言，其开根号后更趋向于1，以求较为保守的抑制效果。此外，噪声区域中的数据点采用第三抑制因子进行处理，第三抑制因子的抑制作用小于第二抑制因子的抑制作用，第三抑制因子可以表示为F3，示例性地，F3＝1。基于此，可以实现组织残留的有效抑制且不影响造影剂强度的效果。应理解，以上近场残留区域、中场残留区域和远场噪声区域的区域划分仅为一种示例，并不用以限制本发明。
由于被差分的一帧造影信号数据可以表示为第二矩阵，第二矩阵中处于不同区域的数据点可以采用不同的抑制因子来处理，因此，抑制因子也可以表示为一个矩阵，将该矩阵与第二矩阵进行点对点相乘的方式，即可实现对造影数据的组织残留印制处理。如前所述的，抑制因子趋向于0，意味着抑制效果的逐渐加强；抑制因子趋向于1，则意味着抑制效果的弱化，甚至不作处理。
当然，前述以单阈值划分的数值区间、抑制因子的计算方式以及抑制因子的作用方式仅是示例性的，还可以根据本申请的教导采用多个阈值划分数值区间，抑制因子也可以是其他的计算方式(例如平方和开根号以外的其他计算方式)，抑制因子的作用方式也可以是其他的方式(例如不是点对点相乘，而是相加、相除等方式)。
基于上面的描述，根据本申请实施例的超声造影成像方法700从时域差分(帧间差分)的角度挖掘组织残留与造影剂信号的差异性，能够对造影信号数据中的组织残留区域与造影剂信号区域进行有效区分，从而能够有效抑制组织残留区域且尽可能不影响造影剂强度的效果。
图9示出了根据本申请再一个实施例的超声造影成像方法900的示意性流程图。如图9所示，超声造影成像方法900可以包括如下步骤：
在步骤S910，控制超声探头向目标对象发射超声波，接收超声波的回波，并基于超声波的回波获取超声回波信号。
在步骤S920，基于超声回波信号获取造影信号数据，并基于造影信号数据生成造影图像。
在步骤S930，获取至少两帧造影图像，计算相邻两帧造影图像之间的信息差异，并计算信息差异与相邻两帧造影图像中的一帧造影图像之间的关系。
在步骤S940，基于关系将一帧造影图像进行区域划分，并对一帧造影图像的不同区域采用不同的抑制因子进行处理。
在根据本申请实施例的超声造影成像方法900中，通过实时的方式获取至少两帧造影图像数据，并基于相邻两帧造影图像数据之间的信息差异与相邻两帧造影图像数据中的一帧造影图像数据之间的关系，将一帧造影图像数据进行区域划分，不同区域的造影图像数据采用不同的抑制因子进行处理。
因此，与前文的根据本申请实施例的超声造影成像方法800相类似的，在根据本申请实施例的超声造影成像方法900中，也是基于相邻两帧造影数据之间的信息差异与相邻两帧造影数据中的一帧造影数据之间的关系对一帧造影数据进行组织残留抑制处理。不同的是，在根据本申请实施例的超声造影成像方法800中，是基于相邻两帧造影图像数据之间的信息差异与相邻两帧造影图像数据中的一帧造影图像数据之间的关系将一帧造影图像数据进行区域划分，而不是基于相邻两帧造影信号数据之间的信息差异与相邻两帧造影信号数据中的一帧造影信号数据之间的关系将一帧造影信号数据进行区域划分。为了简洁，此处主要描述本申请实施例的超声造影成像方法900与本申请实施例的超声造影成像方法800的不同之处，其余相似之处仅简明概括。
在本申请的实施例中，基于相邻两帧造影图像数据之间的信息差异与相邻两帧造影图像数据中的一帧造影图像数据之间的关系，能够将一帧造影图像数据中组织残留强的区域(如前文图4中所示的框410内的区域)和组织残留弱的区域(如前文图4中所示的框420内的区域)和基本无组织残留的区域(如前文图4中所示的框430内的区域)明显地区分开来，使得能够基于该关系对将一帧造影图像数据进行分区(例如分为组织残留强的区域、组织残留弱的区域、无组织残留的区域等)，并基于区域划分的结果对一帧造影图像数据的不同区域采用不同的抑制因子进行抑制处理，从而能够使得经处理的造影图像组织残留部分被有效抑制且不影响造影剂强度的效果。
在本申请的实施例中，相邻两帧造影图像包括前一帧造影图像和后一帧造影图像，步骤S930中的计算相邻两帧造影图像之间的信息差异，可以包括：计算后一帧造影图像与前一帧造影图像的前向差分结果；步骤S930中的计算信息差异与相邻两帧造影图像中的一帧造影图像之间的关系，可以包括：计算前向差分结果与前一帧造影图像的灰度比值。
在本申请的实施例中，相邻两帧造影图像包括前一帧造影图像和后一帧造影图像，步骤S930中的计算相邻两帧造影图像之间的信息差异，可以包括：计算前一帧造影图像与后一帧造影图像的后向差分结果；步骤S930中的计算信息差异与相邻两帧造影图像中的一帧造影图像之间的关系可以包括：计算后向差分结果与后一帧造影图像的灰度比值。
在本申请的实施例中，信息差异包括第一矩阵，被差分的一帧造影图像包括第二矩阵，基于关系将被差分的一帧造影图像进行区域划分，可以包括：确定第一矩阵中的数据点与第二矩阵中的相应位置处的数据点的灰度比值所属的数值区间；基于灰度比值所属的数值区间将第二矩阵中的相应位置处的数据点划分为属于相应的区域。
在本申请的实施例中，数值区间是基于一个或更多个预设阈值来确定的。
在本申请的实施例中，数值区间包括基于预设阈值确定的第一数值区间和第二数值区间；第一数值区间对应第一区域，第二数值区间对应第二区域，且第一区域中的组织残留多于第二区域中的组织残留；当第一矩阵中的任一数据点与第二矩阵中的相应位置处的数据点的灰度比值属于第一数值区间时，对第一区域中的数据点采用第一抑制因子进行处理；当第一矩阵中的任一数据点与第二矩阵中的相应位置处的数据点的灰度比值属于第二数值区间时，对第二区域中的数据点采用第二抑制因子进行处理，第一抑制因子的抑制作用大于第二抑制因子的抑制作用。
在本申请的实施例中，第一区域为近场残留区域，第二区域为中场残留区域。
在本申请的实施例中，数值区间还包括对应噪声区域的第三数值区间；当第一矩阵中的任一数据点与第二矩阵中的相应位置处的数据点的灰度比值属于第三数值区间时，对噪声区域中的数据点采用第三抑制因子进行处理，第三抑制因子的抑制作用小于第二抑制因子的抑制作用。
在本申请的实施例中，第三抑制因子等于1。
基于上面的描述，根据本申请实施例的超声造影成像方法900从时域差分(帧间差分)的角度挖掘组织残留与造影剂信号的差异性，能够对造影图像数据中的组织残留区域与造影剂图像区域进行有效区分，从而能够有效抑制组织残留区域且尽可能不影响造影剂强度的效果。
图10是根据本申请又一个实施例的超声造影成像方法1000的示意性流程图。如图10所示，超声造影成像方法100可以包括如下步骤：
在步骤S1010，控制超声探头向目标对象发射超声波，接收所超声波的回波，并基于超声波的回波获取超声回波信号。
在步骤S1020，基于超声回波信号获取至少两帧造影信号数据，其中，每帧造影信号数据用于生成一帧造影图像。
在步骤S1030，计算相邻两帧造影信号数据之间的信息差异，并计算信息差异与相邻两帧造影信号数据中的一帧造影信号数据之间的关系。
在步骤S1040，基于关系将该一帧造影信号数据进行分类，并基于分类的结果对该一帧造影信号数据的至少部分数据进行抑制处理。
在步骤S1050，基于经处理的一帧造影信号数据生成造影图像。
根据本申请实施例的超声造影成像方法1000与前文的根据本申请实施例的超声造影成像方法300整体上类似，区别仅在于根据本申请实施例的超声造影成像方法1000在步骤S1010-S1020为实时获取两帧造影数据，且获取的造影数据为用于生成造影图像的造影信号数据。在本申请的实施例中，在步骤S1010-S1020获取的造影信号数据可以为基带数据，相邻两帧造影信号数据之间的信息差异包括幅度差异、相位差异或频率差异。
基于上面的描述，根据本申请实施例的超声造影成像方法1000从时域差分(帧间差分)的角度挖掘组织残留与造影剂信号的差异性，能够对造影信号数据中的组织残留区域与造影剂信号区域进行有效区分，从而能够有效抑制组织残留区域且尽可能不影响造影剂强度的效果。
图11示出了根据本申请再一个实施例的超声造影成像方法1100的示意性流程图。如图10所示，超声造影成像方法1100可以包括如下步骤：
在步骤S1110，控制超声探头向目标对象发射超声波，接收超声波的回波，并基于超声波的回波获取超声回波信号。
在步骤S1120，基于超声回波信号获取造影信号数据，并基于造影信号数据生成造影图像。
在步骤S1130，计算相邻两帧造影图像之间的信息差异，并计算信息差异与相邻两帧造影图像中的一帧造影图像之间的关系。
在步骤S1140，基于关系将该一帧造影图像进行分类，并基于分类的结果对至少部分的一帧造影图像进行抑制处理。
根据本申请实施例的超声造影成像方法1100与前文的根据本申请实施例的超声造影成像方法300整体上类似，区别仅在于根据本申请实施例的超声造影成像方法1100在步骤S1110-S1120实时获取至少两帧造影图像。在本申请的实施例中，步骤S1130中的计算相邻两帧造影图像之间的信息差异，可以包括：计算相邻两帧造影图像之间的灰度差分结果；步骤S1130中的计算信息差异与相邻两帧造影图像中的一帧造影图像之间的关系，可以包括：计算灰度差分结果与相邻两帧造影图像中的一帧造影图像之间的灰度比值。
基于上面的描述，根据本申请实施例的超声造影成像方法1100从时域差分(帧间差分)的角度挖掘组织残留与造影剂信号的差异性，能够对造影图像数据中的组织残留区域与造影剂图像区域进行有效区分，从而能够有效抑制组织残留区域且尽可能不影响造影剂强度的效果。
重新参照图2，本发明实施例还提供了一种超声造影成像装置100，超声造影成像装置可以用于实现上述造影成像方法。超声造影成像装置100可以包括探头110、发射电路112、接收电路114、波束合成电路116、处理器118和显示器120、发射/接收选择开关122以及存储器124中的部分或全部部件，各个部件的相关描述可以参照上文。
其中，发射电路112用于激励超声探头110向目标对象发射超声波。接收电路114，用于控制超声探头110接收从目标对象返回的超声回波，获得超声回波信号。处理器118可用于执行上述超声造影成像方法300、800、900、1000和1100。以下仅对超声造影成像装置100的主要功能进行描述，而省略以上已经描述过的细节内容。
在一个实施例中，处理器118计算相邻两帧造影数据之间的信息差异，并计算所述信息差异与所述相邻两帧造影数据中的一帧造影数据之间的关系；基于所述关系将所述一帧造影数据进行分类，并基于所述分类的结果对所述一帧造影数据的至少部分数据进行抑制处理。
在一个实施例中，相邻两帧造影数据之间的信息差异包括幅度差异、相位差异或频率差异。
在一个实施例中，相邻两帧造影数据之间的信息差异为幅度差异时，处理器118计算所述相邻两帧造影数据之间的差分数据，并计算所述差分数据和所述一帧造影数据各自的幅度数据；基于所述差分数据的幅度数据和所述一帧造影数据的幅度数据计算所述差分数据与所述一帧造影数据的幅度比值，得到所述差分数据与所述一帧造影数据之间的关系结果。
在一个实施例中，相邻两帧造影数据之间的信息差异为幅度差异时，处理器118计算所述相邻两帧造影数据之间的差分数据，对所述差分数据与所述一帧造影数据求比值；对所述比值进行幅度求解得到所述差分数据与所述一帧造影数据的幅度比值，作为所述差分数据与所述一帧造影数据之间的关系结果。
在一个实施例中，相邻两帧造影数据之间的信息差异为相位差异时，处理器118计算所述相邻两帧造影数据之间的差分数据，并计算所述差分数据和所述一帧造影数据各自的相位数据；基于所述差分数据的相位数据和所述一帧造影数据的相位数据计算所述差分数据与所述一帧造影数据的相位差值，得到所述差分数据与所述一帧造影数据之间的关系结果。
在一个实施例中，相邻两帧造影数据之间的信息差异为相位差异时，处理器118计算所述相邻两帧造影数据之间的差分数据，对所述差分数据与所述一帧造影数据求比值；对所述比值进行相位求解，得到所述差分数据与所述一帧造影数据的相位差值，作为所述差分数据与所述一帧造影数据之间的关系结果。
在一个实施例中，所述差分数据包括第一矩阵，所述一帧造影数据包括第二矩阵，处理器118用于确定所述第一矩阵中的数据点与所述第二矩阵中的相应位置处的数据点的关系结果所属的数值区间，并基于所述关系结果所属的数值区间将所述第二矩阵中的所述相应位置处的数据点划分为属于相应的类别。
在一个实施例中，处理器118用于基于所述分类的结果确定所述一帧造影数据中需要进行抑制处理的数据；对于所述需要进行抑制处理的数据中不同类别的数据采用不同的抑制因子进行处理。
在一个实施例中，数值区间是基于一个或更多个预设阈值来确定的。
在一个实施例中，所述数值区间包括基于所述预设阈值确定的第一数值区间和第二数值区间；所述第一数值区间对应第一类别数据点，所述第二数值区间对应第二类别数据点，且所述第一类别数据点处的组织残留强于所述第二类别数据点处的组织残留；当所述第一矩阵中的任一数据点与所述第二矩阵中的相应位置处的数据点的关系结果属于所述第一数值区间时，将所述第二矩阵中的相应位置处的数据点划分为第一类别数据点；当所述第一矩阵中的任一数据点与所述第二矩阵中的相应位置处的数据点的关系结果属于所述第二数值区间时，将所述第二矩阵中的相应位置处的数据点划分为第二类别数据点；采用第一抑制因子对所述第一类别数据点进行处理，采用第二抑制因子对所述第二类别数据点进行处理，所述第一抑制因子的抑制作用大于所述第二抑制因子的抑制作用。
在一个实施例中，处理器118计算相邻两帧造影数据之间的信息差异，并计算所述信息差异与所述相邻两帧造影数据中的一帧造影数据之间的关系；基于所述关系将所述一帧造影数据进行区域划分，并对不同区域的造影数据采用不同的抑制因子进行处理。造影数据可以是用于生成一帧造影图像的造影信号数据，也可以是一帧造影图像数据。
下面结合图12描述根据本申请另一个实施例的超声造影成像装置，其可以用于实施上文中的根据本发明实施例的超声造影成像方法。
下面结合图12描述本申请另一个实施例的超声造影成像装置的示意性框图。图12示出了根据本申请实施例的超声造影成像装置1200的示意性框图。超声造影成像装置1200包括存储器1210以及处理器1220。
其中，存储器1210存储用于实现根据本申请实施例的超声造影成像方法300、800、900、1000和1100中的相应步骤的程序。处理器1220用于运行存储器1210中存储的程序，以执行根据本申请实施例的超声造影成像方法300、800、900、1000和1100的相应步骤。
此外，根据本申请实施例，还提供了一种存储介质，在所述存储介质上存储了程序指令，在所述程序指令被计算机或处理器运行时用于执行本申请实施例的超声造影成像方法的相应步骤。所述存储介质例如可以包括智能电话的存储卡、平板电脑的存储部件、个人计算机的硬盘、只读存储器(ROM)、可擦除可编程只读存储器(EPROM)、便携式紧致盘只读存储器(CD-ROM)、USB存储器、或者上述存储介质的任意组合。所述计算机可读存储介质可以是一个或多个计算机可读存储介质的任意组合。
此外，根据本申请实施例，还提供了一种计算机程序，该计算机程序可以存储在云端或本地的存储介质上。在该计算机程序被计算机或处理器运行时用于执行本申请实施例的超声造影成像方法的相应步骤。
基于上面的描述，根据本申请实施例的超声造影成像方法、装置和存储介质从时域差分(帧间差分)的角度挖掘组织残留与造影剂信号的差异性，能够对造影数据中的组织残留部分与造影剂信号部分进行有效区分，从而能够有效抑制组织残留部分且尽可能不影响造影剂强度的效果。
尽管这里已经参考附图描述了示例实施例，应理解上述示例实施例仅仅是示例性的，并且不意图将本申请的范围限制于此。本领域普通技术人员可以在其中进行各种改变和修改，而不偏离本申请的范围和精神。所有这些改变和修改意在被包括在所附权利要求所要求的本申请的范围之内。
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本申请的范围。
在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个装置，或一些特征可以忽略，或不执行。
在此处所提供的说明书中，说明了大量具体细节。然而，能够理解，本申请的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。
类似地，应当理解，为了精简本申请并帮助理解各个发明方面中的一个或多个，在对本申请的示例性实施例的描述中，本申请的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而，并不应将该本申请的方法解释成反映如下意图：即所要求保护的本申请要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说，如相应的权利要求书所反映的那样，其发明点在于可以用少于某个公开的单个实施例的所有特征的特征来解决相应的技术问题。因此，遵循具体实施方式的权利要求书由此明确地并入该具体实施方式，其中每个权利要求本身都作为本申请的单独实施例。
本领域的技术人员可以理解，除了特征之间相互排斥之外，可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者装置的所有过程或单元进行组合。除非另外明确陈述，本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。
此外，本领域的技术人员能够理解，尽管在此所述的一些实施例包括其它实施例中所包括的某些特征而不是其它特征，但是不同实施例的特征的组合意味着处于本申请的范围之内并且形成不同的实施例。例如，在权利要求书中，所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
本申请的各个部件实施例可以以硬件实现，或者以在一个或者多个处理器上运行的软件模块实现，或者以它们的组合实现。本领域的技术人员应当理解，可以在实践中使用微处理器或者数字信号处理器(DSP)来实现根据本申请实施例的一些模块的一些或者全部功能。本申请还可以实现为用于执行这里所描述的方法的一部分或者全部的装置程序(例如，计算机程序和计算机程序产品)。这样的实现本申请的程序可以存储在计算机可读介质上，或者可以具有一个或者多个信号的形式。这样的信号可以从因特网网站上下载得到，或者在载体信号上提供，或者以任何其他形式提供。
应该注意的是上述实施例对本申请进行说明而不是对本申请进行限制，并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中，不应将位于括号之间的任何参考符号构造成对权利要求的限制。本申请可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中，这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。
以上所述，仅为本申请的具体实施方式或对具体实施方式的说明，本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本申请的保护范围之内。本申请的保护范围应以权利要求的保护范围为准。

Claims (42)

1.一种超声造影成像方法，其特征在于，所述方法包括：
控制超声探头向目标对象发射超声波，接收所述超声波的回波，并基于所述超声波的回波获取超声回波信号；
基于所述超声回波信号获取至少两帧造影信号数据，其中，每帧造影信号数据用于生成一帧造影图像；
计算相邻两帧造影信号数据之间的信息差异，并计算所述信息差异与所述相邻两帧造影信号数据中的一帧造影信号数据之间的关系；
基于所述关系将所述一帧造影信号数据进行区域划分，并对不同区域的造影信号数据采用不同的抑制因子进行处理；以及
基于经处理的所述一帧造影信号数据生成造影图像。
2.根据权利要求1所述的方法，其特征在于，所述相邻两帧造影信号数据包括前一帧造影信号数据和后一帧造影信号数据，
所述计算相邻两帧造影信号数据之间的信息差异，包括：计算所述后一帧造影信号数据与所述前一帧造影信号数据的前向差分结果；
所述计算所述信息差异与所述相邻两帧造影信号数据中的一帧造影信号数据之间的关系，包括：计算所述前向差分结果与所述前一帧造影信号数据的大小关系。
3.根据权利要求1所述的方法，其特征在于，所述相邻两帧造影信号数据包括前一帧造影信号数据和后一帧造影信号数据，
所述计算相邻两帧造影信号数据之间的信息差异，包括：计算所述前一帧造影信号数据与所述后一帧造影信号数据的后向差分结果；
所述计算所述信息差异与所述相邻两帧造影信号数据中的一帧造影信号数据之间的关系，包括：计算所述后向差分结果与所述后一帧造影信号数据的大小关系。
4.根据权利要求1所述的方法，其特征在于，基于所述超声回波信号获取的造影信号数据为对所述超声回波信号至少进行波束合成处理后得到的造影信号数据，所述相邻两帧造影信号数据之间的信息差异包括幅度差异、相位差异或频率差异。
5.根据权利要求1所述的方法，其特征在于，基于所述超声回波信号获取的造影信号数据为对所述超声回波信号至少进行模数转换处理后得到的造影信号数据，所述相邻两帧造影信号数据之间的信息差异包括幅度差异。
6.根据权利要求4或5所述的方法，其特征在于，计算相邻两帧造影信号数据之间的幅度差异，并计算所述幅度差异与所述相邻两帧造影信号数据中的一帧造影信号数据之间的关系，包括：
计算所述相邻两帧造影信号数据之间的差分数据，对所述差分数据和所述相邻两帧造影信号数据中的一帧造影信号数据分别进行包络检测，得到所述差分数据和所述一帧造影信号数据各自的幅度数据；
基于所述差分数据的幅度数据和所述一帧造影信号数据的幅度数据计算所述差分数据与所述一帧造影信号数据的幅度比值，得到所述差分数据与所述一帧造影信号数据之间的关系结果。
7.根据权利要求6所述的方法，其特征在于，所述方法还包括：
对所述包络检测之后得到的结果进行对数压缩，基于所述幅度比值将所述一帧造影信号数据进行区域划分时是基于经对数压缩的幅度比值。
8.根据权利要求4或5所述的方法，其特征在于，计算相邻两帧造影信号数据之间的幅度差异，并计算所述幅度差异与所述相邻两帧造影信号数据中的一帧造影信号数据之间的关系，包括：
计算所述相邻两帧造影信号数据之间的差分数据，对所述差分数据与所述一帧造影信号数据求比值；
对所述比值进行包络检测，得到所述差分数据与所述一帧造影信号数据的幅度比值，作为所述差分数据与所述一帧造影信号数据之间的关系结果。
9.根据权利要求4所述的方法，其特征在于，计算相邻两帧造影信号数据之间的相位差异，并计算所述相位差异与所述相邻两帧造影信号数据中的一帧造影信号数据之间的关系，包括：
计算所述相邻两帧造影信号数据之间的差分数据，并计算所述差分数据和所述一帧造影信号数据各自的相位数据；
基于所述差分数据的相位数据和所述一帧造影信号数据的相位数据计算所述差分数据与所述一帧造影信号数据的相位差值，得到所述差分数据与所述一帧造影信号数据之间的关系结果。
10.根据权利要求4所述的方法，其特征在于，计算相邻两帧造影信号数据之间的相位差异，并计算所述相位差异与所述相邻两帧造影信号数据中的一帧造影信号数据之间的关系，包括：
计算所述相邻两帧造影信号数据之间的差分数据，对所述差分数据与所述一帧造影信号数据求比值；
对所述比值进行辐角检测，得到所述差分数据与所述一帧造影信号数据的相位差值，作为所述差分数据与所述一帧造影信号数据之间的关系结果。
11.根据权利要求9-10中的任一项所述的方法，其特征在于，所述差分数据包括第一矩阵，所述一帧造影信号数据包括第二矩阵，所述基于所述关系将所述一帧造影信号数据进行区域划分，包括：
确定所述第一矩阵中的数据点与所述第二矩阵中的相应位置处的数据点的关系结果所属的数值区间；
基于所述关系结果所属的数值区间将所述第二矩阵中的所述相应位置处的数据点划分为属于相应的区域。
12.根据权利要求11所述的方法，其特征在于，所述数值区间是基于一个或更多个预设阈值来确定的。
13.根据权利要求12所述的方法，其特征在于，所述数值区间是基于一个预设阈值T确定的，所述数值区间包括第一数值区间和第二数值区间，所述第一数值区间为(0,T]，所述第二数值区间为(T,1]；
当所述第一矩阵中的任一数据点与所述第二矩阵中的相应位置处的数据点的关系结果属于所述第一数值区间时，将所述第二矩阵中的相应位置处的数据点划分为近场残留区域中的数据点，所述近场残留区域中的数据点采用第一抑制因子进行处理；
当所述第一矩阵中的任一数据点与所述第二矩阵中的相应位置处的数据点的关系结果属于所述第二数值区间时，将所述第二矩阵中的相应位置处的数据点划分为中场残留区域中的数据点，所述中场残留区域中的数据点采用第二抑制因子进行处理，所述第一抑制因子的抑制作用大于所述第二抑制因子的抑制作用。
14.根据权利要求13所述的方法，其特征在于，所述数值区间还包括第三数值区间，所述第三数值区间为(-∞,0]∪(1,+∞]；
当所述第一矩阵中的任一数据点与所述第二矩阵中的相应位置处的数据点的关系结果属于所述第三数值区间时，将所述第二矩阵中的相应位置处的数据点划分为噪声区域中的数据点，所述噪声区域中的数据点采用第三抑制因子进行处理，所述第三抑制因子的抑制作用小于所述第二抑制因子的抑制作用，其中所述第三抑制因子等于1。
15.一种超声造影成像方法，其特征在于，所述方法包括：
控制超声探头向目标对象发射超声波，接收所述超声波的回波，并基于所述超声波的回波获取超声回波信号；
基于所述超声回波信号获取造影信号数据，并基于所述造影信号数据生成造影图像；
获取至少两帧造影图像，计算相邻两帧造影图像之间的信息差异，并计算所述信息差异与所述相邻两帧造影图像中的一帧造影图像之间的关系；以及
基于所述关系将所述一帧造影图像进行区域划分，并对所述一帧造影图像的不同区域采用不同的抑制因子进行处理。
16.根据权利要求15所述的方法，其特征在于，所述相邻两帧造影图像包括前一帧造影图像和后一帧造影图像，
所述计算相邻两帧造影图像之间的信息差异，包括：计算所述后一帧造影图像与所述前一帧造影图像的前向差分结果；
所述计算所述信息差异与所述相邻两帧造影图像中的一帧造影图像之间的关系，包括：计算所述前向差分结果与所述前一帧造影图像的灰度比值。
17.根据权利要求15所述的方法，其特征在于，所述相邻两帧造影图像包括前一帧造影图像和后一帧造影图像，
所述计算相邻两帧造影图像之间的信息差异，包括：计算所述前一帧造影图像与所述后一帧造影图像的后向差分结果；
所述计算所述信息差异与所述相邻两帧造影图像中的一帧造影图像之间的关系，包括：计算所述后向差分结果与所述后一帧造影图像的灰度比值。
18.根据权利要求15-17中的任一项所述的方法，其特征在于，所述信息差异包括第一矩阵，所述一帧造影图像包括第二矩阵，所述基于所述关系将所述一帧造影图像进行区域划分，包括：
确定所述第一矩阵中的数据点与所述第二矩阵中的相应位置处的数据点的灰度比值所属的数值区间；
基于所述灰度比值所属的数值区间将所述第二矩阵中的所述相应位置处的数据点划分为属于相应的区域。
19.根据权利要求18所述的方法，其特征在于，所述数值区间是基于一个或更多个预设阈值来确定的。
20.根据权利要求19所述的方法，其特征在于，所述数值区间包括基于所述预设阈值确定的第一数值区间和第二数值区间；所述第一数值区间对应第一区域，所述第二数值区间对应第二区域，且所述第一区域中的组织残留多于所述第二区域中的组织残留；
当所述第一矩阵中的任一数据点与所述第二矩阵中的相应位置处的数据点的灰度比值属于所述第一数值区间时，对所述第一区域中的数据点采用第一抑制因子进行处理；
当所述第一矩阵中的任一数据点与所述第二矩阵中的相应位置处的数据点的灰度比值属于所述第二数值区间时，对所述第二区域中的数据点采用第二抑制因子进行处理，所述第一抑制因子的抑制作用大于所述第二抑制因子的抑制作用。
21.根据权利要求20所述的方法，其特征在于，所述第一区域为近场残留区域，所述第二区域为中场残留区域。
22.根据权利要求20或21所述的方法，其特征在于，所述数值区间还包括对应噪声区域的第三数值区间；
当所述第一矩阵中的任一数据点与所述第二矩阵中的相应位置处的数据点的灰度比值属于所述第三数值区间时，对所述噪声区域中的数据点采用第三抑制因子进行处理，所述第三抑制因子的抑制作用小于所述第二抑制因子的抑制作用。
23.根据权利要求22所述的方法，其特征在于，所述第三抑制因子等于1。
24.一种超声造影成像装置，其特征在于，所述装置包括超声探头、发射/接收序列控制器和处理器，其中，
所述发射/接收序列控制器用于：激励所述超声探头向目标对象发射超声波，接收所述超声波的回波，基于所述超声波的回波获取超声回波信号；
所述处理器用于执行权利要求1-23中的任一项所述的超声造影成像方法。
25.一种超声造影成像方法，其特征在于，所述方法包括：
获取至少两帧造影数据；
计算相邻两帧造影数据之间的信息差异，并计算所述信息差异与所述相邻两帧造影数据中的一帧造影数据之间的关系；
基于所述关系将所述一帧造影数据进行分类，并基于所述分类的结果对所述一帧造影数据的至少部分数据进行抑制处理。
26.根据权利要求25所述的方法，其特征在于，所述相邻两帧造影数据包括前一帧造影数据和后一帧造影数据，
所述计算相邻两帧造影数据之间的信息差异，包括：计算所述后一帧造影数据与所述前一帧造影数据的前向差分结果；
所述计算所述信息差异与所述相邻两帧造影数据中的一帧造影数据之间的关系，包括：计算所述前向差分结果与所述前一帧造影数据的大小关系。
27.根据权利要求25所述的方法，其特征在于，所述相邻两帧造影数据包括前一帧造影数据和后一帧造影数据，
所述计算相邻两帧造影数据之间的信息差异，包括：计算所述前一帧造影数据与所述后一帧造影数据的后向差分结果；
所述计算所述信息差异与所述相邻两帧造影数据中的一帧造影数据之间的关系，包括：计算所述后向差分结果与所述后一帧造影数据的大小关系。
28.根据权利要求25所述的方法，其特征在于，所述相邻两帧造影数据之间的信息差异包括幅度差异、相位差异或频率差异。
29.根据权利要求28所述的方法，其特征在于，相邻两帧造影数据之间的信息差异为幅度差异时，所述计算相邻两帧造影数据之间的幅度差异，并计算所述幅度差异与所述相邻两帧造影数据中的一帧造影数据之间的关系，包括：
计算所述相邻两帧造影数据之间的差分数据，并计算所述差分数据和所述一帧造影数据各自的幅度数据；
基于所述差分数据的幅度数据和所述一帧造影数据的幅度数据计算所述差分数据与所述一帧造影数据的幅度比值，得到所述差分数据与所述一帧造影数据之间的关系结果。
30.根据权利要求28所述的方法，其特征在于，相邻两帧造影数据之间的信息差异为幅度差异时，所述计算相邻两帧造影数据之间的幅度差异，并计算所述幅度差异与所述相邻两帧造影数据中的一帧造影数据之间的关系，包括：
计算所述相邻两帧造影数据之间的差分数据，对所述差分数据与所述一帧造影数据求比值；
对所述比值进行幅度求解得到所述差分数据与所述一帧造影数据的幅度比值，作为所述差分数据与所述一帧造影数据之间的关系结果。
31.根据权利要求28所述的方法，其特征在于，相邻两帧造影数据之间的信息差异为相位差异时，所述计算相邻两帧造影数据之间的相位差异，并计算所述相位差异与所述相邻两帧造影数据中的一帧造影数据之间的关系，包括：
计算所述相邻两帧造影数据之间的差分数据，并计算所述差分数据和所述一帧造影数据各自的相位数据；
基于所述差分数据的相位数据和所述一帧造影数据的相位数据计算所述差分数据与所述一帧造影数据的相位差值，得到所述差分数据与所述一帧造影数据之间的关系结果。
32.根据权利要求28所述的方法，其特征在于，相邻两帧造影数据之间的信息差异为相位差异时，所述计算相邻两帧造影数据之间的相位差异，并计算所述相位差异与所述相邻两帧造影数据中的一帧造影数据之间的关系，包括：
计算所述相邻两帧造影数据之间的差分数据，对所述差分数据与所述一帧造影数据求比值；
对所述比值进行相位求解，得到所述差分数据与所述一帧造影数据的相位差值，作为所述差分数据与所述一帧造影数据之间的关系结果。
33.根据权利要求29-32中的任一项所述的方法，其特征在于，所述差分数据包括第一矩阵，所述一帧造影数据包括第二矩阵，所述基于所述关系将所述一帧造影数据进行分类，包括：
确定所述第一矩阵中的数据点与所述第二矩阵中的相应位置处的数据点的关系结果所属的数值区间；
基于所述关系结果所属的数值区间将所述第二矩阵中的所述相应位置处的数据点划分为属于相应的类别。
34.根据权利要求25-32中的任一项所述的方法，其特征在于，所述基于所述分类的结果对所述一帧造影数据中的至少部分数据进行抑制处理，包括：
基于所述分类的结果确定所述一帧造影数据中需要进行抑制处理的数据；
对于所述需要进行抑制处理的数据中不同类别的数据采用不同的抑制因子进行处理。
35.根据权利要求33所述的方法，其特征在于，所述数值区间是基于一个或更多个预设阈值来确定的。
36.根据权利要求35所述的方法，其特征在于，所述数值区间包括基于所述预设阈值确定的第一数值区间和第二数值区间；所述第一数值区间对应第一类别数据点，所述第二数值区间对应第二类别数据点，且所述第一类别数据点处的组织残留强于所述第二类别数据点处的组织残留；
当所述第一矩阵中的任一数据点与所述第二矩阵中的相应位置处的数据点的关系结果属于所述第一数值区间时，将所述第二矩阵中的相应位置处的数据点划分为第一类别数据点；
当所述第一矩阵中的任一数据点与所述第二矩阵中的相应位置处的数据点的关系结果属于所述第二数值区间时，将所述第二矩阵中的相应位置处的数据点划分为第二类别数据点；
采用第一抑制因子对所述第一类别数据点进行处理，采用第二抑制因子对所述第二类别数据点进行处理，所述第一抑制因子的抑制作用大于所述第二抑制因子的抑制作用。
37.一种超声造影成像方法，其特征在于，所述方法包括：
控制超声探头向目标对象发射超声波，接收所述超声波的回波，并基于所述超声波的回波获取超声回波信号；
基于所述超声回波信号获取至少两帧造影信号数据，其中，每帧造影信号数据用于生成一帧造影图像；
计算相邻两帧造影信号数据之间的信息差异，并计算所述信息差异与所述相邻两帧造影信号数据中的一帧造影信号数据之间的关系；
基于所述关系将所述一帧造影信号数据进行分类，并基于所述分类的结果对所述一帧造影信号数据的至少部分数据进行抑制处理；以及
基于经处理的所述一帧造影信号数据生成造影图像。
38.根据权利要求37所述的方法，其特征在于，所述获取的造影信号数据为基带数据，所述相邻两帧造影信号数据之间的信息差异包括幅度差异、相位差异或频率差异。
39.一种超声造影成像方法，其特征在于，所述方法包括：
控制超声探头向目标对象发射超声波，接收所述超声波的回波，并基于所述超声波的回波获取超声回波信号；
基于所述超声回波信号获取造影信号数据，并基于所述造影信号数据生成造影图像；
计算相邻两帧造影图像之间的信息差异，并计算所述信息差异与所述相邻两帧造影图像中的一帧造影图像之间的关系；以及
基于所述关系将所述一帧造影图像进行分类，并基于所述分类的结果对至少部分的所述一帧造影图像进行抑制处理。
40.根据权利要求39所述的方法，其特征在于，所述计算相邻两帧造影图像之间的信息差异，包括：计算相邻两帧造影图像之间的灰度差分结果；
所述计算所述信息差异与所述相邻两帧造影图像中的一帧造影图像之间的关系，包括：计算所述灰度差分结果与所述相邻两帧造影图像中的一帧造影图像之间的灰度比值。
41.一种超声造影成像装置，其特征在于，所述装置包括超声探头、发射/接收序列控制器和处理器，其中：
所述发射/接收序列控制器用于：激励所述超声探头向目标对象发射超声波，接收所述超声波的回波，基于所述超声波的回波获取超声回波信号；
所述处理器用于：基于所述超声回波信号获取至少两帧造影数据；计算相邻两帧造影数据之间的信息差异，并计算所述信息差异与所述相邻两帧造影数据中的一帧造影数据之间的关系；基于所述关系将所述一帧造影数据进行分类，并基于所述分类的结果对所述一帧造影数据的至少部分数据进行抑制处理。
42.一种存储介质，其特征在于，所述存储介质上存储有计算机程序，所述计算机程序在运行时执行如权利要求1-23以及25-40中的任一项所述的超声造影成像方法。
CN202080006417.2A 2020-05-18 2020-05-18 超声造影成像方法、装置和存储介质 Active CN113164160B (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2020/090855 WO2021232192A1 (zh) 2020-05-18 2020-05-18 超声造影成像方法、装置和存储介质

Publications (2)

Publication Number Publication Date
CN113164160A CN113164160A (zh) 2021-07-23
CN113164160B true CN113164160B (zh) 2022-09-16

Family

ID=76879265

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202080006417.2A Active CN113164160B (zh) 2020-05-18 2020-05-18 超声造影成像方法、装置和存储介质

Country Status (2)

Country Link
CN (1) CN113164160B (zh)
WO (1) WO2021232192A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20240114218A (ko) * 2023-01-16 2024-07-23 삼성메디슨 주식회사 초음파 진단 장치 및 진단 방법

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106028947A (zh) * 2015-05-07 2016-10-12 深圳迈瑞生物医疗电子股份有限公司 血流成像方法及系统
WO2016161574A1 (zh) * 2015-04-08 2016-10-13 深圳迈瑞生物医疗电子股份有限公司 一种超声造影成像方法及装置
CN107361793A (zh) * 2017-07-18 2017-11-21 深圳开立生物医疗科技股份有限公司 超声波成像方法、系统及超声成像设备
CN108135566A (zh) * 2016-06-30 2018-06-08 北京深迈瑞医疗电子技术研究院有限公司 一种增强超声造影图像的方法、系统及超声造影成像设备
CN110458836A (zh) * 2019-08-16 2019-11-15 深圳开立生物医疗科技股份有限公司 一种超声造影成像方法、装置和设备及可读存储介质
CN111012381A (zh) * 2019-12-30 2020-04-17 无锡祥生医疗科技股份有限公司 超声造影方法、装置及存储介质

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102525616B1 (ko) * 2015-10-08 2023-04-26 삼성메디슨 주식회사 조영제 초음파 진단 장치 및 방법
US20190154821A1 (en) * 2017-11-21 2019-05-23 The Charles Stark Draper Laboratory, Inc. System and method for imaging and localization of contrast-enhanced features in the presence of accumulating contrast agent in a body

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016161574A1 (zh) * 2015-04-08 2016-10-13 深圳迈瑞生物医疗电子股份有限公司 一种超声造影成像方法及装置
CN106028947A (zh) * 2015-05-07 2016-10-12 深圳迈瑞生物医疗电子股份有限公司 血流成像方法及系统
CN108135566A (zh) * 2016-06-30 2018-06-08 北京深迈瑞医疗电子技术研究院有限公司 一种增强超声造影图像的方法、系统及超声造影成像设备
CN107361793A (zh) * 2017-07-18 2017-11-21 深圳开立生物医疗科技股份有限公司 超声波成像方法、系统及超声成像设备
CN110458836A (zh) * 2019-08-16 2019-11-15 深圳开立生物医疗科技股份有限公司 一种超声造影成像方法、装置和设备及可读存储介质
CN111012381A (zh) * 2019-12-30 2020-04-17 无锡祥生医疗科技股份有限公司 超声造影方法、装置及存储介质

Also Published As

Publication number Publication date
CN113164160A (zh) 2021-07-23
WO2021232192A1 (zh) 2021-11-25

Similar Documents

Publication Publication Date Title
KR101407425B1 (ko) 최적 스펙클 추적을 위한 이중 경로 처리
JP5472914B2 (ja) 超音波診断装置
CN102429684B (zh) 一种多普勒彩色血流成像方法和装置
US20120154400A1 (en) Method of reducing noise in a volume-rendered image
US9261485B2 (en) Providing color doppler image based on qualification curve information in ultrasound system
US10540769B2 (en) Method and system for enhanced ultrasound image visualization by detecting and replacing acoustic shadow artifacts
CN110800019B (zh) 用于复合超声图像生成的方法和系统
US20110002518A1 (en) Method and system for processing ultrasound data
US20150379700A1 (en) Ultrasound image displaying apparatus and method for displaying ultrasound image
US8705802B2 (en) Providing a motion image in an ultrasound system
US20090024033A1 (en) Ultrasound diagnostic apparatus
US12004908B2 (en) Ultrasound imaging system for high resolution wideband harmonic imaging
CN111513754A (zh) 一种超声成像设备、超声图像的质量评估方法
CN106600550A (zh) 超声图像处理方法及系统
CN113164160B (zh) 超声造影成像方法、装置和存储介质
US6909796B2 (en) Imaging system and method
US20080097212A1 (en) Automatic identification of orientation in medical diagnostic ultrasound
US20110190632A1 (en) Ultrasonic diagnostic appratus and ultrasonic image processing method
CN114727760A (zh) 一种光声成像方法及光声成像系统
CN109363722B (zh) 彩色血流成像中运动伪像的抑制方法及设备
CN114246610A (zh) 用于生成增强图像的超声成像系统和方法
Hyun et al. A GPU-based real-time spatial coherence imaging system
JP2018187126A (ja) 超音波診断装置及び超音波画像生成方法
JP3977779B2 (ja) 超音波診断装置
CN115990036A (zh) 超声成像方法及设备

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant