CN111882515B - 超声信号处理方法、超声信号处理设备及存储介质 - Google Patents
超声信号处理方法、超声信号处理设备及存储介质 Download PDFInfo
- Publication number
- CN111882515B CN111882515B CN202011037549.XA CN202011037549A CN111882515B CN 111882515 B CN111882515 B CN 111882515B CN 202011037549 A CN202011037549 A CN 202011037549A CN 111882515 B CN111882515 B CN 111882515B
- Authority
- CN
- China
- Prior art keywords
- ultrasonic image
- image signal
- signal
- ultrasonic
- confidence
- 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
Links
- 238000012545 processing Methods 0.000 title claims abstract description 35
- 238000003672 processing method Methods 0.000 title claims abstract description 25
- 238000002604 ultrasonography Methods 0.000 claims abstract description 135
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 20
- 238000005070 sampling Methods 0.000 claims description 41
- 238000000034 method Methods 0.000 claims description 27
- 230000006870 function Effects 0.000 claims description 15
- 238000012937 correction Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 abstract description 7
- 239000000523 sample Substances 0.000 description 10
- 230000008569 process Effects 0.000 description 7
- 230000000295 complement effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000013507 mapping Methods 0.000 description 4
- 230000001629 suppression Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000012285 ultrasound imaging Methods 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 230000001427 coherent effect Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000009210 therapy by ultrasound Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Quality & Reliability (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种超声信号处理方法，包括以下步骤：获取各个阵元通道接收到的射频信号，以及至少两组预存的变迹权重值；基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号；获取每一所述超声图像信号对应的置信度估计结果，以及每一所述超声图像信号相对于其它超声图像信号的相关性估计结果，其中，所述其它超声图像信号为所述至少两个超声图像信号中，除自身之外的其它超声图像信号；基于所述置信度估计结果、所述相关性估计结果以及所述超声图像信号合成目标超声图像信号。本发明还公开了一种超声信号处理设备及计算机可读存储介质，达成了提升超声图像的分辨率的效果。
Description
技术领域
本发明涉及超声技术领域，尤其涉及超声信号处理方法、超声信号处理设备及计算机可读存储介质。
背景技术
超声成像被广泛应用于多科室的医学实践中，用于辅助医生进行诊断以及治疗。现有的基于多变迹技术的杂波抑制方法一般通过采用2个互补的变迹来得到一对互补的阵元全通道信号，进而通过这一对信号的互相关来得到一个对应的信号权重，并最后将权重与原本正常的DAS（Delay And Sum，延时求和）结果相乘，以得到杂波抑制后的图像信号。
但是，上述的信号处理方法，只会将图像杂波压制掉, 效果与手动降低图像增益区别不大，无法提升超声图像的分辨率，导致最终得出的超声图像的分辨率较低。
上述内容仅用于辅助理解本发明的技术方案，并不代表承认上述内容是现有技术。
发明内容
本发明的主要目的在于提供一种超声信号处理方法、超声信号处理设备及计算机可读存储介质，旨在实现提升超声图像的分辨率的效果。
为实现上述目的，本发明提供一种超声信号处理方法，所述超声信号处理方法包括以下步骤：
获取各个阵元通道接收到的射频信号，以及至少两组预存的变迹权重值；
基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号；
获取每一所述超声图像信号对应的置信度估计结果，以及每一所述超声图像信号相对于其它超声图像信号的相关性估计结果，其中，所述其它超声图像信号为所述至少两个超声图像信号中，除自身之外的其它超声图像信号；
基于所述置信度估计结果、所述相关性估计结果以及所述超声图像信号合成目标超声图像信号。
可选地，所述基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号的步骤包括：
基于第一组所述预存的变迹权重值对所述射频信号进行变迹处理，得到第一变迹射频信号；
基于所述第一变迹射频信号合成第一超声图像信号；
基于第二组所述预存的变迹权重值对所述射频信号进行变迹处理，得到第二变迹射频信号；
基于所述第二变迹射频信号合成第二超声图像信号。
可选地，所述获取每一所述超声图像信号对应的置信度估计结果，以及每一所述超声图像信号相对于其它超声图像信号的相关性估计结果的步骤包括：
根据每一所述超声图像信号中的各个采样点的置信度，确定所述超声图像信号对应的置信度估计结果；
获取每一所述超声图像信号相对于所述其它超声图像信号的相关性值，作为所述超声图像信号对应的相关性估计结果。
可选地，所述根据每一所述超声图像信号中的各个采样点的置信度，确定所述超声图像信号对应的置信度估计结果的步骤包括：
获取第一超声图像信号中的各个采样点的置信度，作为所述第一超声图像信号对应的置信度估计结果；
获取第二超声图像信号中的各个所述采样点的置信度，作为所述第二超声图像信号对应的置信度估计结果。
可选地，所述获取第一超声图像信号中的各个采样点的置信度，作为所述第一超声图像信号对应的置信度估计结果的步骤包括：
获取所述第一超声图像信号中的各个所述采样点对应的局域信号的统计参数；
计算各个所述采样点对应的局域信号的标准统计参数；
根据每一所述采样点对应的所述统计参数和所述标准统计参数，确定每一所述采样点对应的置信度，并将各个所述采样点的置信度，作为所述第一超声图像信号对应的置信度估计结果。
可选地，所述基于所述置信度估计结果、所述相关性估计结果以及所述超声图像信号合成目标超声图像信号的步骤包括：
根据每一所述超声图像信号对应的所述置信度估计结果和所述相关性估计结果，确定每一所述超声图像信号对应的权重函数；
根据所述权重函数及所述超声图像信号合成所述目标超声图像信号。
可选地，所述基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号的步骤之前，还包括：
对所述射频信号进行延时校正，以基于校正后的所述射频信号以及所述预存的变迹权重值生成所述超声图像信号。
可选地，所述基于所述置信度估计、所述相关性估计以及所述超声图像信号合成目标超声图像信号的步骤之后，还包括：
基于所述目标超声图像信号生成超声图像。
此外，为实现上述目的，本发明还提供一种超声信号处理设备，所述超声信号处理设备包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的超声信号处理程序，所述超声信号处理程序被所述处理器执行时实现如上所述的超声信号处理方法的步骤。
此外，为实现上述目的，本发明还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有超声信号处理程序，所述超声信号处理程序被处理器执行时实现如上所述的超声信号处理方法的步骤。
本发明实施例提出的一种超声信号处理方法、超声信号处理设备及计算机可读存储介质，先获取各个阵元通道接收到的射频信号，以及至少两组预存的变迹权重值，然后基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号，进而获取每一所述超声图像信号对应的置信度估计结果，以及每一所述超声图像信号相对于其它超声图像信号的相关性估计结果，其中，所述其它超声图像信号为所述至少两个超声图像信号中，除自身之外的其它超声图像信号，并基于所述置信度估计结果、所述相关性估计结果以及所述超声图像信号合成目标超声图像信号。由于超声图像信号综合利用了每组变迹处理后的射频信号的信息，因而能够有效地抽取同一个采样点在不同情况下的特征,经过相关性估计结果和置信度估计结果加权后能够更好地抑制杂波, 而保留原本杂波区域中的信号, 因此在提高对比度的同时, 能够一定程度上提升图像分辨率, 并且由于加权的鲁棒性,使得图像中不会产生肉眼可见的伪影。
附图说明
图1是本发明实施例方案涉及的硬件运行环境的终端结构示意图；
图2为本发明超声信号处理方法的一实施例的流程示意图；
图3为本发明超声信号处理方法的另一实施例的流程示意图；
图4为本发明超声信号处理方法的又一实施例的流程示意图。
本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明实施例提出的超声信号的处理方法的主要解决方案是：
获取各个阵元通道接收到的射频信号，以及至少两组预存的变迹权重值；
基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号；
获取每一所述超声图像信号对应的置信度估计结果，以及每一所述超声图像信号相对于其它超声图像信号的相关性估计结果，其中，所述其它超声图像信号为所述至少两个超声图像信号中，除自身之外的其它超声图像信号；
基于所述置信度估计结果、所述相关性估计结果以及所述超声图像信号合成目标超声图像信号。
由于超声图像信号综合利用了每组变迹处理后的射频信号的信息，因而能够有效地抽取同一个采样点在不同情况下的特征,经过相关性估计结果和置信度估计结果加权后能够更好地抑制杂波, 而保留原本杂波区域中的信号, 因此在提高对比度的同时, 能够一定程度上提升图像分辨率, 并且由于加权的鲁棒性,使得图像中不会产生肉眼可见的伪影。
如图1所示，图1是本发明实施例方案涉及的硬件运行环境的终端结构示意图。
本发明实施例终端可以是超声信号处理设备。
如图1所示，该终端可以包括：处理器1001，例如CPU，用户接口1003，存储器1004，通信总线1002。其中，通信总线1002用于实现这些组件之间的连接通信。用户接口1003可以包括显示屏（Display）等，可选用户接口1003还可以包括标准的有线接口、无线接口。存储器1004可以是高速RAM存储器，也可以是稳定的存储器（non-volatile memory），例如磁盘存储器。存储器1004可选的还可以是独立于前述处理器1001的存储装置。
本领域技术人员可以理解，图1中示出的终端结构并不构成对终端的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
如图1所示，作为一种计算机存储介质的存储器1004中可以包括操作系统、用户接口模块以及超声信号处理程序。
在图1所示的终端中，处理器1001可以用于调用存储器1004中存储的超声信号处理程序，并执行以下操作：
获取各个阵元通道接收到的射频信号，以及至少两组预存的变迹权重值；
基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号；
获取每一所述超声图像信号对应的置信度估计结果，以及每一所述超声图像信号相对于其它超声图像信号的相关性估计结果，其中，所述其它超声图像信号为所述至少两个超声图像信号中，除自身之外的其它超声图像信号；
基于所述置信度估计结果、所述相关性估计结果以及所述超声图像信号合成目标超声图像信号。
进一步地，处理器1001可以调用存储器1004中存储的超声信号处理程序，还执行以下操作：
基于第一组所述预存的变迹权重值对所述射频信号进行变迹处理，得到第一变迹射频信号；
基于所述第一变迹射频信号合成第一超声图像信号；
基于第二组所述预存的变迹权重值对所述射频信号进行变迹处理，得到第二变迹射频信号；
基于所述第二变迹射频信号合成第二超声图像信号。
进一步地，处理器1001可以调用存储器1004中存储的超声信号处理程序，还执行以下操作：
根据每一所述超声图像信号中的各个采样点的置信度，确定所述超声图像信号对应的置信度估计结果；
获取每一所述超声图像信号相对于所述其它超声图像信号的相关性值，作为所述超声图像信号对应的相关性估计结果。
进一步地，处理器1001可以调用存储器1004中存储的超声信号处理程序，还执行以下操作：
获取第一超声图像信号中的各个采样点的置信度，作为所述第一超声图像信号对应的置信度估计结果；
获取第二超声图像信号中的各个所述采样点的置信度，作为所述第二超声图像信号对应的置信度估计结果。
进一步地，处理器1001可以调用存储器1004中存储的超声信号处理程序，还执行以下操作：
获取所述第一超声图像信号中的各个所述采样点对应的局域信号的统计参数；
计算各个所述采样点对应的局域信号的标准统计参数；
根据每一所述采样点对应的所述统计参数和所述标准统计参数，确定每一所述采样点对应的置信度，并将各个所述采样点的置信度，作为所述第一超声图像信号对应的置信度估计结果。
进一步地，处理器1001可以调用存储器1004中存储的超声信号处理程序，还执行以下操作：
根据每一所述超声图像信号对应的所述置信度估计结果和所述相关性估计结果，确定每一所述超声图像信号对应的权重函数；
根据所述权重函数及所述超声图像信号合成所述目标超声图像信号。
进一步地，处理器1001可以调用存储器1004中存储的超声信号处理程序，还执行以下操作：
对所述射频信号进行延时校正，以基于校正后的所述射频信号以及所述预存的变迹权重值生成所述超声图像信号。
进一步地，处理器1001可以调用存储器1004中存储的超声信号处理程序，还执行以下操作：
基于所述目标超声图像信号生成超声图像。
超声成像技术广泛应用于多科室的医学实践中，用于辅助医生进行诊断以及治疗。一般来说，现代医学超声成像系统的基础均为DAS的波束合成成像法。其基本原理在于，先通过超声探头向人体组织发射超声波，然后超声探头的各阵元接收不同时间的散射回波，并对各阵元的回波信号作延时校正处理后，再进行相干叠加，以获得高分辨率的超声图像。
在上述过程中，超声探头的阵元数量对于最终图像的效果影响比较大，因为每个阵元都能获得组织的信息。因此，阵元越多获得的组织信息越丰富，图像质量越好。另外由于声波的干涉特性，在进行相干叠加时，如果各阵元的延时不准确, 则导致图像的信噪比显著降低。而延时不准确的原因主要是系统设计和真实情况的不一致。在超声系统设计时，往往是基于一系列基本假设（如组织均匀和只考虑主瓣等）进行设计的，而在实际超声成像时，真实组织显然是非均匀的。并且，真实探头发射超声时的旁瓣不可忽略。正是由于这些设计与真实使用状态下的差异导致图像信号中含有大量不可忽略的杂波，从而导致图像劣化。
为了提高图像质量，现有的基于多变迹技术的杂波抑制方法一般通过采用2个互补的变迹来得到一对互补的阵元全通道信号，进而通过这一对信号的互相关来得到一个对应的信号权重，并最后将权重与原本正常的DAS（Delay And Sum，延时求和）结果相乘，以得到杂波抑制后的图像信号。
但是，上述的信号处理方法，只会将图像杂波压制掉, 效果与手动降低图像增益区别不大，无法提升超声图像的分辨率，导致最终得出的超声图像的分辨率较低。并且，由于上述方案中，通过单纯的压制杂波占比大的信号，还会将原本的杂波区域的真实组织也压制掉许多，从而导致最终输出的超声图像出现artifacts（伪影）。
为解决现有滤波方案导致超声图像分辨率较低，且会导致最终输出的超声图像出现artifacts的缺陷，本发明提出一种超声信号处理方法。以下通过结合附图，通过具体的实施例对本发明提出的超声信号处理方法进行解释说明，以使本领域技术人员理解本发明权利要求要求保护的范围。需要强调的是，以下实施例并不对本发明的保护范围进行限定，只用于解释对本发明提出的超声信号处理方法的具体内容进行解释说明。
参照图2，在本发明超声信号处理方法的一实施例中，所述超声信号处理方法包括以下步骤：
步骤S10、获取各个阵元通道接收到的射频信号，以及至少两组预存的变迹权重值；
步骤S20、基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号；
步骤S30、获取每一所述超声图像信号对应的置信度估计，以及每一所述超声图像信号相对于其它超声图像信号的相关性估计，其中，所述其它超声图像信号为所述至少两个超声图像信号中，除自身之外的其它超声图像信号；
步骤S40、基于所述置信度估计、所述相关性估计以及所述超声图像信号合成目标超声图像信号。
在本实施例中，上述超声信号处理方法的各个步骤执行于超声信号处理设备中。
作为一种实现方案，上述超声处理设备可以包括探头和显示组件，上述探头用于发射和接收超声波。其中，上述压电换能器将施加在其上面的电压脉冲激励转换成机械振动从而对外发射超声波。上述显示组件用于显示超声图像。可以理解的是，探头一般由若干大小相同的长条形压电换能器（单个的压电换能器称为阵元）间隔排列组成。但是本实施例并不对上述阵元的排列方式进行限定，例如，上述阵元也可以排列成二维矩阵形状。
需要说明的是，上述超声处理设备也可以仅作为一个数据处理单元，而不包括上述探头和显示组件。其仅用于接收由其它具备探头的终端采集的RF(Radio Frequency，射频)信号，并进行处理后，输出一目标超声图像信号至用于显示超声图像的终端。
上述超声信号处理设备可以先获取各个阵元通道接收到的射频信号，其中，射频信号的个数由探头设置的阵元数确定，当探头设置有n个阵元时，对应的射频频信号也为n个。
当获取到多个射频信号后，还可以获取预先保存于超声信号处理设备的存储介质中的，预存的变迹权重值。可以理解的是，在本实施例中，超声信号处理设备中，预先保存有多组变迹权重值。使得超声信号处理设备可以基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号。
具体地，每一组变迹权重值中，均设置有与各个射频信号一一对应的变迹权重值。因此，可以基于一组变迹权重值中与各个射频信号一一对应的变迹权重值对各个射频信号进行变迹处理后，通过对变迹处理后的各个射频信号进行通道求和，得到一超声图像信号。其中，由于预先保存有多组不同的变迹权重值，因此，可以对应得到多个不同模式的超声图像信号。
需要说明的是，上述变迹处理优选为相位加权变迹。但是，在一些实施方案中，也可以采用幅度加权变迹或者孔径加权变迹。在此本实施例不作具体限定，本领域技术人员可以根据实际需求，选择性修改。另外，当上述变迹处理过程设置为相位变迹时，上述与射频信号一一对应的权重值为一复数。此外，上述变迹权重值的组数可以一般可以设置为2-10组，其具体组数设置，可以由本领域技术人员基于用于实施本方案的设备的硬件参数确定。
示例性地，可以用以下式子表示第j个超声图像信号S j ：
其中，n为射频信号的个数（即阵元个数），apo(j,i)为第j个超声图像中，第i个射频信号对应的权重值，ch(i)为第i个射频信号。
进一步地，在得到多个不同模式的超声图像信号后，可以获取每一所述超声图像信号对应的置信度估计结果，以及每一所述超声图像信号相对于其它超声图像信号的相关性估计，其中，所述其它超声图像信号为所述至少两个超声图像信号中，除自身之外的其它超声图像信号。即可以获取第一超声图像信号中的各个采样点的置信度，作为所述第一超声图像信号对应的置信度估计结果，获取第二超声图像信号中的各个所述采样点的置信度，作为所述第二超声图像信号对应的置信度估计结果。
具体地，可以获取所述第一超声图像信号中的各个所述采样点对应的局域信号的统计参数，然后计算各个所述采样点对应的局域信号的标准统计参数，以根据每一所述采样点对应的所述统计参数和所述标准统计参数，确定每一所述采样点对应的置信度，并将各个所述采样点的置信度，作为所述第一超声图像信号对应的置信度估计结果。
示例性地，可以根据预设的统计量对超声图像信号标准统计分布的距离来确定，置信度。即可以基于一下公式计算第j个超声图像信号S j 的第x个采样点的置信度P x ：
其中，m为统计量的数量（m的具体数值不作限定，可以根据需求设置多个统计量），W k 为第k个统计量的权重值，为预设的经验值，Q xk 为第j个超声图像信号S j 的第x个采样点的第k个统计量的数值（例如，为第x个采样点对应的局域均值），Q standard 为代表通过标准分布或者模拟计算给出的标准统计量。
需要说明的是，超声图像信号标准统计分布可以通过理想的瑞利分布或者Nakagami分布来近似计算, 也可以通过蒙特卡罗模拟做大量实验来确定。
进一步地，还可以获取每一所述超声图像信号相对于所述其它超声图像信号的相关性值，作为所述超声图像信号对应的相关性估计结果。即可以将全部超声图像信号集合在一起计算所述相关性估计结果, 用以表征某个超声图像信号与其它超声图像信号的相似性, 从而表征其信号有多接近真实信号。可以理解的是, 这里的相关性估计的输入是所有的超声图像信号, 输出也会针对性地输出每一个超声图像信号相对于其它超声图像信号相关性估计结果corr。
需要说明的是，本实施例并不限定具体的相关性估计方式，例如，可以通过欧氏距离、KL散度估计等方式确定相关性估计结果。但是，在计算相关性时, 要采用局域ROI(region of interest，感兴趣区域)的信号综合判断。例如，在有n个超声图像信号时，n个超声图像信号中的第i个超声图像信号对应的相关性估计结果corr（i；1,2，…i-1，i+1…n）可以表示为如下式所示：
进一步地，在确定每一超声图像信号对应的相关性估计结果和置信度估计结果后，可以基于所述置信度估计结果、所述相关性估计结果以及所述超声图像信号合成目标超声图像信号。即根据每一所述超声图像信号对应的所述置信度估计结果和所述相关性估计结果，确定每一所述超声图像信号对应的权重函数，并根据所述权重函数及所述超声图像信号合成所述目标超声图像信号。
示例性地，对每个超声图像信号，可以结合其对应的置信度估计结果与相关性估计结果，得到其对应的权重函数。例如，所述权重函数Fi可以表示为：
其中，A、B和G代表置信度估计结果和相关性估计结果均作一次函数映射后, 再相乘, 并将结果再作一次函数映射得到最终的权重函数。这里的函数映射A，B和G的选择是多样化的, 最简单的是直接令F i =corr/P，也可以有一定效果。或者也可以根据需要采取指数映射得到F i =exp(-corr/P)，这样有归一化的效果, 在实现中可以选择和设计的函数方式非常多，因此，本实施例并不限定一定要采用哪一种。
进一步地，得到各个超声图像信号对应的权重函数后，可以基于所述权重函数及所述超声图像信号合成所述目标超声图像信号。其合成过程可以表示为下式：
其中，output为所述目标超声图像信号。
在本实施例公开的技术方案中，先获取各个阵元通道接收到的射频信号，以及至少两组预存的变迹权重值，然后基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号，进而获取每一所述超声图像信号对应的置信度估计结果，以及每一所述超声图像信号相对于其它超声图像信号的相关性估计结果，其中，所述其它超声图像信号为所述至少两个超声图像信号中，除自身之外的其它超声图像信号，并基于所述置信度估计结果、所述相关性估计结果以及所述超声图像信号合成目标超声图像信号。由于超声图像信号综合利用了每组变迹处理后的射频信号的信息，因而能够有效地抽取同一个采样点在不同情况下的特征,经过相关性估计结果和置信度估计结果加权后能够更好地抑制杂波, 而保留原本杂波区域中的信号, 因此在提高对比度的同时, 能够一定程度上提升图像分辨率, 并且由于加权的鲁棒性,使得图像中不会产生肉眼可见的伪影。
参照图3，基于上述实施例，在另一实施例中，所述步骤S20之前，还包括：
步骤S50、对所述射频信号进行延时校正，以基于校正后的所述射频信号以及所述预存的变迹权重值生成所述超声图像信号。
在本实施例中，当超声信号处理设备获取到各个阵元对应的射频信号后，可以先对射频信号进行延时校正来各对应采样点的正确信号值，以提高信号准确度。
可选地，在一些实现方案中，由于各通道中的模拟信号会经过低噪放及时间增益补偿等常规模拟信号放大处理过程, 以方便后续的信号处理，然后通过模数转换，将之前的模拟信号转换为数字信号, 得到各阵元通道的数字RF信号。
参照图4，基于上述任一实施例，在又一实施例中，所述步骤S40之后，还包括：
步骤S60、基于所述目标超声图像信号生成超声图像。
在本实施例中，超声信号处理设备还可以基于所述目标超声图像信号生成超声图像，并将该超声图像显示于显示组件中。
可选地，在一些实施方案中，还可以对上述超声图像进行后处理，其中，所述后处理包括空间复合，增益，动态范围调整和/或图像增强等过程。
此外，本发明实施例还提出一种超声信号处理设备，所述超声信号处理设备包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的超声信号处理程序，所述超声信号处理程序被所述处理器执行时实现如上各个实施例所述的超声信号处理方法的步骤。
此外，本发明实施例还提出一种计算机可读存储介质，所述计算机可读存储介质上存储有超声信号处理程序，所述超声信号处理程序被处理器执行时实现如上各个实施例所述的超声信号处理方法的步骤。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者系统不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者系统所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者系统中还存在另外的相同要素。
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在如上所述的一个存储介质(如ROM/RAM、磁碟、光盘)中，包括若干指令用以使得一台终端设备(可以是超声信号处理设备等)执行本发明各个实施例所述的方法。
以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。
Claims (10)
1.一种超声信号处理方法，其特征在于，所述超声信号处理方法包括以下步骤：
获取各个阵元通道接收到的射频信号，以及至少两组预存的变迹权重值；
基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号；
获取每一所述超声图像信号对应的置信度估计结果，以及每一所述超声图像信号相对于其它超声图像信号的相关性估计结果，其中，所述其它超声图像信号为所述至少两个超声图像信号中，除自身之外的其它超声图像信号；
基于所述置信度估计结果、所述相关性估计结果，以及基于至少两组所述预存的变迹权重值及所述射频信号生成的至少两个超声图像信号，合成目标超声图像信号。
2.如权利要求1所述的超声信号处理方法，其特征在于，所述基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号的步骤包括：
基于第一组所述预存的变迹权重值对所述射频信号进行变迹处理，得到第一变迹射频信号；
基于所述第一变迹射频信号合成第一超声图像信号；
基于第二组所述预存的变迹权重值对所述射频信号进行变迹处理，得到第二变迹射频信号；
基于所述第二变迹射频信号合成第二超声图像信号。
3.如权利要求1所述的超声信号处理方法，其特征在于，所述获取每一所述超声图像信号对应的置信度估计结果，以及每一所述超声图像信号相对于其它超声图像信号的相关性估计结果的步骤包括：
根据每一所述超声图像信号中的各个采样点的置信度，确定每一所述超声图像信号对应的置信度估计结果；
获取每一所述超声图像信号相对于所述其它超声图像信号的相关性值，作为所述超声图像信号对应的相关性估计结果。
4.如权利要求3所述的超声信号处理方法，其特征在于，所述根据每一所述超声图像信号中的各个采样点的置信度，确定每一所述超声图像信号对应的置信度估计结果的步骤包括：
获取第一超声图像信号中的各个采样点的置信度，作为所述第一超声图像信号对应的置信度估计结果；
获取第二超声图像信号中的各个所述采样点的置信度，作为所述第二超声图像信号对应的置信度估计结果。
5.如权利要求4所述的超声信号处理方法，其特征在于，所述获取第一超声图像信号中的各个采样点的置信度，作为所述第一超声图像信号对应的置信度估计结果的步骤包括：
获取所述第一超声图像信号中的各个所述采样点对应的局域信号的统计参数；
计算各个所述采样点对应的局域信号的标准统计参数；
根据每一所述采样点对应的所述统计参数和所述标准统计参数，确定每一所述采样点对应的置信度，并将各个所述采样点的置信度，作为所述第一超声图像信号对应的置信度估计结果。
6.如权利要求1所述的超声信号处理方法，其特征在于，所述基于所述置信度估计结果、所述相关性估计结果，以及基于至少两组所述预存的变迹权重值及所述射频信号生成的至少两个超声图像信号，合成目标超声图像信号的步骤包括：
根据每一所述超声图像信号对应的所述置信度估计结果和所述相关性估计结果，确定每一所述超声图像信号对应的权重函数；
根据所述权重函数及基于至少两组所述预存的变迹权重值及所述射频信号生成的至少两个超声图像信号，合成所述目标超声图像信号。
7.如权利要求1所述的超声信号处理方法，其特征在于，所述基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号的步骤之前，还包括：
对所述射频信号进行延时校正，以基于校正后的所述射频信号执行所述基于至少两组所述预存的变迹权重值及所述射频信号生成至少两个超声图像信号的步骤。
8.如权利要求1所述的超声信号处理方法，其特征在于，所述基于所述置信度估计、所述相关性估计，以及基于至少两组所述预存的变迹权重值及所述射频信号生成的至少两个超声图像信号，合成目标超声图像信号的步骤之后，还包括：
基于所述目标超声图像信号生成超声图像。
9.一种超声信号处理设备，其特征在于，所述超声信号处理设备包括：存储器、处理器及存储在所述存储器上并可在所述处理器上运行的超声信号处理程序，所述超声信号处理程序被所述处理器执行时实现如权利要求1至8中任一项所述的超声信号处理方法的步骤。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有超声信号处理程序，所述超声信号处理程序被处理器执行时实现如权利要求1至8中任一项所述的超声信号处理方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011037549.XA CN111882515B (zh) | 2020-09-28 | 2020-09-28 | 超声信号处理方法、超声信号处理设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011037549.XA CN111882515B (zh) | 2020-09-28 | 2020-09-28 | 超声信号处理方法、超声信号处理设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111882515A CN111882515A (zh) | 2020-11-03 |
| CN111882515B true CN111882515B (zh) | 2020-12-29 |
Family
ID=73199201
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011037549.XA Active CN111882515B (zh) | 2020-09-28 | 2020-09-28 | 超声信号处理方法、超声信号处理设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111882515B (zh) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109069119A (zh) * | 2016-04-26 | 2018-12-21 | 皇家飞利浦有限公司 | 用于超声胎儿成像的3d图像合成 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7771356B2 (en) * | 2005-11-14 | 2010-08-10 | Spentech, Inc. | Ultrasound tympanoscope |
| US8235900B2 (en) * | 2009-03-23 | 2012-08-07 | Imsonic Medical, Inc. | Method and apparatus for an automatic ultrasound imaging system |
| CN101940479B (zh) * | 2010-07-15 | 2012-10-03 | 重庆大学 | 基于分段动态变迹系数的超声成像波束合成方法及其装置 |
| CN105193453B (zh) * | 2015-10-15 | 2018-01-30 | 武汉超信电子工程有限公司 | 一种超声成像的区域合成方法 |
| US20200187910A1 (en) * | 2017-05-09 | 2020-06-18 | The University Of North Carolina At Chapel Hill | Adaptive multifocus beamforming ultrasound methods and systems for improved penetration and target sensitivity at high frame-rates |
| CN107997783B (zh) * | 2017-11-29 | 2021-03-30 | 声泰特（成都）科技有限公司 | 一种基于超声波方向性的自适应超声波束合成方法和系统 |
| CN108113703B (zh) * | 2017-12-21 | 2021-02-09 | 飞依诺科技(苏州)有限公司 | 一种用于超声扫描的生成变迹值的方法和装置 |
| CN109754407B (zh) * | 2019-01-10 | 2021-06-01 | 青岛海信医疗设备股份有限公司 | 超声图像处理方法、装置及设备 |
-
2020
- 2020-09-28 CN CN202011037549.XA patent/CN111882515B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109069119A (zh) * | 2016-04-26 | 2018-12-21 | 皇家飞利浦有限公司 | 用于超声胎儿成像的3d图像合成 |
Non-Patent Citations (2)
| Title |
|---|
| 超声医学图像的去噪及增强研究;王绍波;《中国优秀硕士学位论文全文数据库 医药卫生科技辑》;20120515;第E080-48页 * |
| 超声相控阵成像关键算法研究;赵大丹;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20170315;第B022-9页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111882515A (zh) | 2020-11-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Xu et al. | Spatio-temporally smoothed coherence factor for ultrasound imaging [correspondence] | |
| US20060241456A1 (en) | Ultrasonic imaging apparatus and ultrasonic imaging method | |
| EP2673657B1 (en) | An imaging device with image acquisition rate optimization | |
| JP2006346459A (ja) | 超音波空間複合イメージングのリアルタイム運動補正のための方法及び装置 | |
| US20060079780A1 (en) | Ultrasonic imaging apparatus | |
| EP1777543A2 (en) | Aberration correction beam patterns | |
| WO2014115377A1 (ja) | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム | |
| JP2017093913A (ja) | 超音波診断装置、信号処理装置及び解析プログラム | |
| Wang et al. | Short-lag spatial coherence combined with eigenspace-based minimum variance beamformer for synthetic aperture ultrasound imaging | |
| KR101552427B1 (ko) | 초음파 영상의 스페클 저감장치 | |
| JP7199972B2 (ja) | 情報処理装置、情報処理方法、プログラム | |
| KR20160046669A (ko) | 빔포밍 장치, 초음파 이미징 장치 및 빔포밍 방법 | |
| EP2269091A2 (en) | Image forming method using ultrasound and aberration correction method | |
| CN113631953A (zh) | 使用基波和谐波信号的定量超声 | |
| JP2022023982A (ja) | 超音波画像を処理するための方法及びシステム | |
| JP6444518B2 (ja) | 超音波撮像装置 | |
| CN107997784B (zh) | 一种基于声速自适应修正的超声波束合成方法和系统 | |
| CN111882515B (zh) | 超声信号处理方法、超声信号处理设备及存储介质 | |
| JP7171625B2 (ja) | 超音波画像を処理するための方法及びシステム | |
| US20220022848A1 (en) | Ultrasound imaging system using coherence estimation of a beamformed signal | |
| CN111278363B (zh) | 超声成像设备、系统及其超声造影成像的图像增强方法 | |
| US11986356B2 (en) | Reduction of reverberation artifacts in ultrasound images and associated devices, systems, and methods | |
| CN112370073A (zh) | 超声声速校正方法、超声信号处理设备及存储介质 | |
| US11982741B2 (en) | Methods and systems for processing an ultrasound image | |
| KR20160046777A (ko) | 빔포밍 장치, 초음파 이미징 장치 및 빔포밍 방법 |
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 |