CN110613476A - 超声信号的处理方法和装置 - Google Patents
超声信号的处理方法和装置 Download PDFInfo
- Publication number
- CN110613476A CN110613476A CN201810632043.XA CN201810632043A CN110613476A CN 110613476 A CN110613476 A CN 110613476A CN 201810632043 A CN201810632043 A CN 201810632043A CN 110613476 A CN110613476 A CN 110613476A
- Authority
- CN
- China
- Prior art keywords
- variable
- formula
- array element
- initialized
- iteration
- 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.)
- Granted
Links
Classifications
-
- 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
Abstract
本发明实施例提供一种超声信号的处理方法和装置。本发明的超声信号的处理方法，包括：根据所述初始化的第一迭代变量Am(0)和所述初始化的第一中间变量Bm(0)之间的大小关系，对第一迭代变量Am(n)和第二迭代变量am(n)进行迭代，确定阵元m发射的超声信号的回波信号对不同的探测距离的延时km(n)，对第一迭代变量Am(n)进行插值，分别根据插值变量和与第一中间变量Bm(n)之间的大小关系，对延时km(n+1)进行修正。本发明实施例可以实现降低波束合成过程中的延时计算所需的运算量，并且可以提升延时计算的精度值，提升波束合成质量。
Description
技术领域
本发明实施例涉及信号处理技术，尤其涉及一种超声信号的处理方法和装置。
背景技术
随着科学技术的不断发展，超声系统在医学上得到越来越广泛的应用。其中最重要的部分是波束合成，而波束合成的关键是对不同阵元不同探测距离上的延时计算。
超声探头阵列的多个阵元发出探测声波，由于不同的阵元相对于同一组织点的位置不同，该次声波达到组织的路径长短不同，不同的组织对声波的反射也不同，最终声波被反射回来，到达聚焦点位置的时长并不相同，不同的声波束叠加形成一个聚焦点，聚焦图像是由一个个聚焦点形成的。这些用于形成一个聚焦点的声波由于并不同时返回，因此需要计算各自的延时，从而找到属于本次声波的反馈信号的声波以形成聚焦点。波束合成就是对声波进行延时，变迹处理，形成聚焦图像。
对于波束合成要对每一个探测距离进行聚焦(波束合成)，计算延时值，会导致运算量很大，且运算精度不高。
发明内容
本发明实施例提供一种超声信号的处理方法和装置，以实现降低波束合成过程中的延时计算所需的运算量，并且可以提升延时计算的精度值。
第一方面，本发明实施例提供一种超声信号的处理方法，包括：
初始化阵元m的延时变量km(n)，获取阵元m的km(0)；
根据所述km(0)确定初始化的第一迭代变量Am(0)、初始化的第二迭代变量am(0)、初始化的第一中间变量Bm(0)和初始化的第二中间变量bm(0)；
根据所述初始化的第一迭代变量Am(0)和所述初始化的第一中间变量Bm(0)之间的大小关系，对第一迭代变量Am(n)和第二迭代变量am(n)进行迭代，确定阵元m发射的超声信号的回波信号对不同的探测距离的延时km(n)；
对第一迭代变量Am(n)进行插值，分别根据插值变量和与第一中间变量Bm(n)之间的大小关系，对延时km(n+1)进行修正；
其中，m的不同取值分别对应不同的阵元，n分别取0至X。
第二方面，本发明实施例提供一种超声信号的处理装置，包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序，以实现如第一方面所述的方法。
第三方面，本发明实施例提供一种计算机存储介质，包括：所述计算机存储介质用于存储计算机程序，所述计算机程序执行时用于实现如第一方面所述的方法。
本发明实施例的超声信号的处理方法和装置，通过初始化阵元m的延时变量km(n)，获取阵元m的km(0)，根据所述km(0)确定初始化的第一迭代变量Am(0)、初始化的第二迭代变量am(0)、初始化的第一中间变量Bm(0)和初始化的第二中间变量bm(0)，根据所述初始化的第一迭代变量Am(0)和所述初始化的第一中间变量Bm(0)之间的大小关系，对第一迭代变量Am(n)和第二迭代变量am(n)进行迭代，确定阵元m发射的超声信号的回波信号对不同的探测距离的延时km(n)，对第一迭代变量Am(n)进行插值，分别根据插值变量和与第一中间变量Bm(n)之间的大小关系，对延时km(n+1)进行修正，以修正后的延时值对超声信号的回波信号进行延时，以便后续进行空间复合成像，实现降低波束合成过程中的延时计算所需的运算量，并且可以提升延时计算的精度值，提升波束合成质量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的超声信号的处理方法的应用场景示意图；
图2为本发明的超声信号的处理方法实施例一的流程图；
图3为本发明的超声信号的处理方法实施例二的流程图；
图4为探头的类型信息为线性类型时的延时计算示意图；
图5为探头的类型信息为凸阵类型时的延时计算示意图；
图6为本发明的超声信号的处理方法实施例三的流程图；
图7为本发明的超声信号的处理装置实施例一的结构示意图；
图8为本发明的超声信号的处理装置实施例二的结构示意图；
图9为本发明的超声信号的处理装置实施例三的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为本发明的超声信号的处理方法的应用场景示意图，如图1所示，该应用场景包括探头和超声信号的处理装置。该超声信号的处理装置和探头连接。
该超声信号的处理装置可以执行本发明的超声信号的处理方法，以降低波束合成过程中的延时计算所需的运算量，并且可以提升延时计算的精度值，提升波束合成质量。
该探头用于发射和接收超声信号，即可以接收超声信号的回波信号。该超声信号的处理装置可以对超声信号的回波信号进行延时计算，以便后续进行空间复合成像，并输出图像。
需要说明的是，该超声信号的处理装置可以是芯片，也可以是芯片中的一个功能模块，该超声信号的处理装置可以与显示器分开设置，也可以设置在一起。
下面采用几个具体的实施例对本发明的超声信号的处理方法进行具体解释说明。
图2为本发明的超声信号的处理方法实施例一的流程图，如图2所示，本实施例的方法可以包括：
步骤101、初始化阵元m的延时变量km(n)，获取阵元m的km(0)。
其中，m的不同取值分别对应不同的阵元，n分别取0至X，n可以指示不同的探测距离。以线性类型的探头的阵元m为例，该xm为阵元m的坐标。以凸阵类型的探头的阵元m为例，该av(m)为阵元m的曲率半径。其中，表示向下取整。X为任意整数。
步骤102、根据所述km(0)确定初始化的第一迭代变量Am(0)、初始化的第二迭代变量am(0)、初始化的第一中间变量Bm(0)和初始化的第二中间变量bm(0)。
其中，Am(0)为初始化的第一迭代变量，Am(n)为第一迭代变量，Am(0)即为Am(n)在n＝0的取值，Am(n)的具体表达是为下述公式(1)。am(0)为初始化的第二迭代变量，am(n)为第二迭代变量，am(0)即为am(n)在n＝0的取值，am(n)的具体表达是为下述公式(2)。Bm(0)为初始化的第一中间变量，Bm(n)为第一中间变量，Bm(0)即为Bm(n)在n＝0的取值，Bm(n)的具体表达是为下述公式(3)。bm(0)为初始化的第二中间变量，bm(n)为第二中间变量，bm(0)即为bm(n)在n＝0的取值，bm(n)的具体表达是为下述公式(4)。
Am(n)＝n2+2k(n)n+k2(n) (1)
Bm(n)＝n2+αmn+βm (2)
am(n)＝2n+2k(n)+1 (3)
bm(n)＝2n+1+αm (4)
其中，以线性类型的探头的阵元m为例，以凸阵类型的探头的阵元m为例，
可以根据上述公式(1)和km(0)确定初始化的第一迭代变量Am(0)，根据公式(3)和km(0)确定初始化的第二迭代变量am(0)，根据公式(2)和km(0)确定初始化的第一中间变量Bm(0)，根据公式(4)和km(0)确定初始化的第二中间变量bm(0)。
步骤103、根据所述初始化的第一迭代变量Am(0)和所述初始化的第一中间变量Bm(0)之间的大小关系，对第一迭代变量Am(n)和第二迭代变量am(n)进行迭代，确定阵元m发射的超声信号的回波信号对不同的探测距离的延时km(n)。
具体的，本发明实施例采用迭代的方式对上述第一迭代变量Am(n)和第二迭代变量am(n)进行迭代，即使用Am(0)和am(0)确定Am(1)和am(1)，使用Am(1)和am(1)确定Am(2)和am(2)，依次类推，直至确定Am(X)和am(X)。并在迭代过程中可以根据公式(3)确定n在不同取值时的km(n)。
步骤104、对第一迭代变量Am(n)进行插值，分别根据插值变量和与第一中间变量Bm(n)之间的大小关系，对延时km(n+1)进行修正。
具体的，本实施例还可以对第一迭代变量Am(n)进行插值，确定插值变量和在通过比较Am(n)和Bm(n)、和Bm(n)、和Bm(n)之间的大小关系，对延时km(n+1)进行修正，以修正后的延时值对超声信号的回波信号进行延时，以便后续进行空间复合成像。
本实施例，通过初始化阵元m的延时变量km(n)，获取阵元m的km(0)，根据所述km(0)确定初始化的第一迭代变量Am(0)、初始化的第二迭代变量am(0)、初始化的第一中间变量Bm(0)和初始化的第二中间变量bm(0)，根据所述初始化的第一迭代变量Am(0)和所述初始化的第一中间变量Bm(0)之间的大小关系，对第一迭代变量Am(n)和第二迭代变量am(n)进行迭代，确定阵元m发射的超声信号的回波信号对不同的探测距离的延时km(n)，对第一迭代变量Am(n)进行插值，分别根据插值变量和与第一中间变量Bm(n)之间的大小关系，对延时km(n+1)进行修正，以修正后的延时值对超声信号的回波信号进行延时，以便后续进行空间复合成像，实现降低波束合成过程中的延时计算所需的运算量，并且可以提升延时计算的精度值，提升波束合成质量。
下面采用几个具体的实施例，对图2所示方法实施例的技术方案进行详细说明。
图3为本发明的超声信号的处理方法实施例二的流程图，如图3所示，本实施例的方法可以包括：
步骤201、初始化阵元m的延时变量km(n)，获取阵元m的km(0)。
步骤202、根据所述km(0)确定初始化的第一迭代变量Am(0)、初始化的第二迭代变量am(0)、初始化的第一中间变量Bm(0)和初始化的第二中间变量bm(0)。
步骤2031、根据公式δm(n)＝Bm(n)-Am(n)计算δm(n)。
步骤2032、如果δm(n)<0，则使用公式(5)确定am(n+1)，根据公式(6)确定Am(n+1)。
am(n+1)＝am(n) (5)
Am(n+1)＝Am(n) (6)
步骤2033、如果δm(n)≥0，则使用公式(7)确定am(n+1)，根据公式(8)确定Am(n+1)。
am(n+1)＝am(n)+2 (7)
Am(n+1)＝Am(n)+am(n) (8)
步骤2034、根据公式(9)计算km(n)。
km(n)＝(am(n)-1-2n)/2 (9)
步骤2041、根据公式(10)计算根据公式(11)计算
步骤2042、当δm(n)<0，且时，根据公式(12)对km(n+1)进行修正。
k′(n+1)＝k(n+1)+corr (12)
步骤2043、当δm(n)<0，且时，根据公式(13)对km(n+1)进行修正。
k′(n+1)＝k(n+1) (13)
步骤2044、当δm(n)≥0，且时，根据公式(12)对km(n+1)进行修正。
步骤2045、当δm(n)≥0，且时，根据公式(13)对km(n+1)进行修正。
其中，corr＝0.5。k′(n+1)即为k(n+1)的修正后的值。
需要说明的是，执行一次步骤2041至步骤2045，可以提升2倍延时计算精度，具体的，并修正之前的延时计算精度为Tu为所述超声信号的采样周期。通过执行一次步骤2041至步骤2045，可以将延时计算精度提升至
可选的，在执行一次步骤2041至步骤2045之后，还可以根据公式(当δ(n)<0时)或者公式 (当δ(n)≥0时)，对δm(n)进行更新。根据公式对am(n)进行更新。根据公式对corr进行更新。可以使用更新后的δm(n)、am(n)以及corr，通过如上述步骤2041至步骤2044的方式对km(n+1)进行再修正。将延时计算精度提升至N为执行步骤2041的次数。
本实施例，通过初始化阵元m的延时变量km(n)，获取阵元m的km(0)，根据所述km(0)确定初始化的第一迭代变量Am(0)、初始化的第二迭代变量am(0)、初始化的第一中间变量Bm(0)和初始化的第二中间变量bm(0)，根据所述初始化的第一迭代变量Am(0)和所述初始化的第一中间变量Bm(0)之间的大小关系，对第一迭代变量Am(n)和第二迭代变量am(n)进行迭代，确定阵元m发射的超声信号的回波信号对不同的探测距离的延时km(n)，对第一迭代变量Am(n)进行插值，分别根据插值变量和与第一中间变量Bm(n)之间的大小关系，对延时km(n+1)进行修正，以修正后的延时值对超声信号的回波信号进行延时，以便后续进行空间复合成像，实现降低波束合成过程中的延时计算所需的运算量，并且可以提升延时计算的精度值，提升波束合成质量。
图4为探头的类型信息为线性类型时的延时计算示意图，如图4所示，该探头包括多个阵元，每个阵元用一个小方格表示，该多个阵元沿直线排布。以0点为接收中心线位置，θ为接收中心线偏角，两边的阵元同时接收超声信号的回波信号，在示意图中，对于x+的振元来说，偏转角(即扫描波束与振元夹角)大于90°；对于x-的振元来说，偏转角小于90°。
Tu为所述超声信号的采样周期，c0为所述超声信号的传播速度、θ为接收中心线偏角(如图所示)。
如图4所涉及的参数(如r、l等)都使用d来量化(d代表在一个采样周期内往返的声波经历的一半的距离长度，为仅往或返的距离长度)，其中n＝r(nTu)/d，为整数；k(n)＝l(nTu)/d，使用d来量化。
当探头的类型信息为如图4所示的线性类型时，上述实施例中的
并且根据公式Am(0)＝km 2(0)对阵元m的延时变量Am(n)进行初始化，根据公式am(0)＝2km(0)+1对阵元m的延时变量am(n)进行初始化，根据公式Bm(0)＝xm 2/d2对阵元m的延时变量Bm(n)进行初始化，根据公式bm(0)＝2xmsinθ/d+1对阵元m的延时变量bm(n)进行初始化。
需要说明的是，另一种可实现方式，根据预设的数据表和聚焦点的偏转角度，获取与所述偏转角度对应的阵元m的km(0)；其中，所述预设的数据表中包括不同的阵元位置和偏转角度对应的初始化的值。即可以通过查表的方式确定km(0)。
图5为探头的类型信息为凸阵类型时的延时计算示意图，如图5所示，探头的多个阵元不在同一个水平线上，存在一个角度差γ，所以，上述线性类型中的θ就变为了θ+γ，x变为了av，其中av表示为探头的属性数值。对于不同的阵元m(m∈[-K,K])，存在不同的θ+γ(m)和av(m)。m>0对应于线阵中x+的部分，m<0对应于线阵中x-的部分。
当探头的类型信息为如图5所示的凸阵类型时，上述实施例中的
根据公式对阵元m的延时变量Am(n)进行初始化，根据公式am(0)＝2km(0)+1对阵元m的延时变量am(n)进行初始化，根据公式对阵元m的延时变量Bm(n)进行初始化，根据公式bm(0)＝2av(m)sin(θ+γ(m))/d+1对阵元m的延时变量bm(n)进行初始化。
需要说明的是，另一种可实现方式，根据预设的数据表和聚焦点的偏转角度，获取与所述偏转角度对应的阵元m的km(0)；其中，所述预设的数据表中包括不同的阵元位置和偏转角度对应的初始化的值。即可以通过查表的方式确定km(0)。
图6为本发明的超声信号的处理方法实施例三的流程图，如图6所示，本实施例在上述任一实施例的基础上，本实施例的方法还可以包括：
步骤301、根据中心线与法线之间的夹角范围确定有效阵元的范围。
具体的，当探头的类型信息为线性类型时，根据公式计算根据公式计算所述有效阵元的范围为
当探头的类型信息为凸阵类型时，根据公式计算根据公式计算所述有效阵元的范围为
其中，θ为所述中心线与法线之间的夹角，也可称为接收中心线偏角。
步骤302、根据所述有效阵元的范围确定m的取值范围。
具体的，上述实施例中的xm需要在该范围内。上述实施例中的av(m)需要在该范围内。
本实施例，通过根据中心线与法线之间的夹角范围确定有效阵元的范围，根据所述有效阵元的范围确定m的取值范围，从而可以在不影响后续成像的基础上，去除无效阵元的计算，可以进一步提升运算效率。
图7为本发明的超声信号的处理装置实施例一的结构示意图，如图7所示，本实施例的装置可以包括：初始化模块11、迭代计算模块12和修正模块13，其中，初始化模块11用于初始化阵元m的延时变量km(n)，获取阵元m的km(0)；初始化模块11还用于根据所述km(0)确定初始化的第一迭代变量Am(0)、初始化的第二迭代变量am(0)、初始化的第一中间变量Bm(0)和初始化的第二中间变量bm(0)；迭代计算模块12用于根据所述初始化的第一迭代变量Am(0)和所述初始化的第一中间变量Bm(0)之间的大小关系，对第一迭代变量Am(n)和第二迭代变量am(n)进行迭代，确定阵元m发射的超声信号的回波信号对不同的探测距离的延时km(n)；修正模块13用于对第一迭代变量Am(n)进行插值，分别根据插值变量和与第一中间变量Bm(n)之间的大小关系，对延时km(n+1)进行修正；其中，m的不同取值分别对应不同的阵元，n分别取0至X。
该迭代计算模块12用于根据公式δm(n)＝Bm(n)-Am(n)计算δm(n)；
如果δm(n)<0，则am(n+1)＝am(n),Am(n+1)＝Am(n)；
如果δm(n)≥0，则am(n+1)＝am(n)+2,Am(n+1)＝Am(n)+am(n)；
km(n)＝(am(n)-1-2n)/2。
该修正模块13用于对第一迭代变量Am(n)进行插值，包括：根据公式计算根据公式 计算
该修正模块13用于分别根据插值变量和与第一中间变量Bm(n)之间的大小关系，对延时km(n+1)进行修正，包括：
当δm(n)<0，且时，根据公式k′(n+1)＝k(n+1)+corr对km(n+1)进行修正；
当δm(n)<0，且时，根据公式k′(n+1)＝k(n+1)对km(n+1)进行修正；
当δm(n)≥0，且时，根据公式k′(n+1)＝k(n+1)+corr对km(n+1)进行修正；
当δm(n)≥0，且时，根据公式k′(n+1)＝k(n+1)对km(n+1)进行修正；
其中，corr＝0.5。
该初始化模块11用于初始化阵元m的延时变量km(n)，获取阵元m的km(0)，包括：获取探头的类型信息；当所述探头的类型信息为线性类型时，则根据公式对阵元m的延时变量km(n)进行初始化；当所述探头的类型信息为凸阵类型时，则根据公式对阵元m的延时变量km(n)进行初始化；其中，表示向下取整。
当所述探头的类型信息为线性类型时，该初始化模块11用于根据所述km(0)确定初始化的第一迭代变量Am(0)、初始化的第二迭代变量am(0)、初始化的第一中间变量Bm(0)和初始化的第二中间变量bm(0)，包括：
根据公式Am(0)＝km 2(0)对阵元m的延时变量Am(n)进行初始化；
根据公式am(0)＝2km(0)+1对阵元m的延时变量am(n)进行初始化；
根据公式Bm(0)＝xm 2/d2对阵元m的延时变量Bm(n)进行初始化；
根据公式bm(0)＝2xmsinθ/d+1对阵元m的延时变量bm(n)进行初始化。
当所述探头的类型信息为凸阵类型时，该初始化模块11用于根据所述km(0)确定初始化的第一迭代变量Am(0)、初始化的第二迭代变量am(0)、初始化的第一中间变量Bm(0)和初始化的第二中间变量bm(0)，包括：
根据公式对阵元m的延时变量Am(n)进行初始化；
根据公式am(0)＝2km(0)+1对阵元m的延时变量am(n)进行初始化；
根据公式对阵元m的延时变量Bm(n)进行初始化；
根据公式bm(0)＝2av(m)sin(θ+γ(m))/d+1对阵元m的延时变量bm(n)进行初始化。
所述初始化阵元m的延时变量km(n)，获取阵元m的km(0)，包括：
根据预设的数据表和聚焦点的偏转角度，获取与所述偏转角度对应的阵元m的km(0)；
其中，所述预设的数据表中包括不同的阵元位置和偏转角度对应的初始化的值。
本实施例的装置，可以用于执行上述方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
图8为本发明的超声信号的处理装置实施例二的结构示意图，如图8所示，本实施例的装置在图7所示装置结构的基础上，进一步地，还可以包括：有效阵元确定模块14，该有效阵元确定模块14用于根据中心线与法线之间的夹角范围确定有效阵元的范围；根据所述有效阵元的范围确定m的取值范围。
其中，根据中心线与法线之间的夹角范围确定有效阵元的范围，包括：
当探头的类型信息为线性类型时，根据公式计算根据公式计算所述有效阵元的范围为
当探头的类型信息为凸阵类型时，根据公式计算根据公式计算所述有效阵元的范围为
其中，θ为所述中心线与法线之间的夹角。
本实施例的装置，可以用于执行上述方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
图9为本发明的超声信号的处理装置实施例三的结构示意图，如图9所示，本实施例的装置可以包括：存储器21和处理器22，其中，存储器21用于存储计算机程序，处理器22用于执行所述计算机程序，以实现上述各个实施例所述的方法。
本实施例的装置，可以用于执行上述方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
本发明实施例还提供一种计算机存储介质，所述计算机存储介质用于存储计算机程序，所述计算机程序执行时用于实现如上述任一实施例所述的方法。其实现原理和技术效果类似，此处不再赘述。
本领域普通技术人员可以理解：实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时，执行包括上述各方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (12)
1.一种超声信号的处理方法，其特征在于，包括：
初始化阵元m的延时变量km(n)，获取阵元m的km(0)；
根据所述km(0)确定初始化的第一迭代变量Am(0)、初始化的第二迭代变量am(0)、初始化的第一中间变量Bm(0)和初始化的第二中间变量bm(0)；
根据所述初始化的第一迭代变量Am(0)和所述初始化的第一中间变量Bm(0)之间的大小关系，对第一迭代变量Am(n)和第二迭代变量am(n)进行迭代，确定阵元m发射的超声信号的回波信号对不同的探测距离的延时km(n)；
对第一迭代变量Am(n)进行插值，分别根据插值变量和与第一中间变量Bm(n)之间的大小关系，对延时km(n+1)进行修正；
其中，m的不同取值分别对应不同的阵元，n分别取0至X。
2.根据权利要求1所述的方法，其特征在于，所述根据所述初始化的第一迭代变量Am(0)和所述初始化的第一中间变量Bm(0)之间的大小关系，对第一迭代变量Am(n)和第二迭代变量am(n)进行迭代，确定阵元m发射的超声信号的回波信号对不同的探测距离的延时km(n)，包括：
根据公式δm(n)＝Bm(n)-Am(n)计算δm(n)；
如果δm(n)＜0，则am(n+1)＝am(n)，Am(n+1)＝Am(n)；
如果δm(n)≥0，则am(n+1)＝am(n)+2，Am(n+1)＝Am(n)+am(n)；
km(n)＝(am(n)-1-2n)/2。
3.根据权利要求2所述的方法，其特征在于，所述对第一迭代变量Am(n)进行插值，包括：
根据公式计算
根据公式计算
4.根据权利要求3所述的方法，其特征在于，所述分别根据插值变量和与第一中间变量Bm(n)之间的大小关系，对延时km(n+1)进行修正，包括：
当δm(n)＜0，且时，根据公式k′(n+1)＝k(n+1)+corr对km(n+1)进行修正；
当δm(n)＜0，且时，根据公式k′(n+1)＝k(n+1)对km(n+1)进行修正；
当δm(n)≥0，且时，根据公式k′(n+1)＝k(n+1)+corr对km(n+1)进行修正；
当δm(n)≥0，且时，根据公式k′(n+1)＝k(n+1)对km(n+1)进行修正；
其中，corr＝0.5。
5.根据权利要求1至4任一项所述的方法，其特征在于，所述初始化阵元m的延时变量km(n)，获取阵元m的km(0)，包括：
获取探头的类型信息；
当所述探头的类型信息为线性类型时，则根据公式对阵元m的延时变量km(n)进行初始化；
当所述探头的类型信息为凸阵类型时，则根据公式对阵元m的延时变量km(n)进行初始化；
其中，表示向下取整。
6.根据权利要求5所述的方法，其特征在于，当所述探头的类型信息为线性类型时，所述根据所述km(0)确定初始化的第一迭代变量Am(0)、初始化的第二迭代变量am(0)、初始化的第一中间变量Bm(0)和初始化的第二中间变量bm(0)，包括：
根据公式Am(0)＝km 2(0)对阵元m的延时变量Am(n)进行初始化；
根据公式am(0)＝2km(0)+1对阵元m的延时变量am(n)进行初始化；
根据公式Bm(0)＝xm 2/d2对阵元m的延时变量Bm(n)进行初始化；
根据公式bm(0)＝2xmsinθ/d+1对阵元m的延时变量bm(n)进行初始化。
7.根据权利要求5所述的方法，其特征在于，当所述探头的类型信息为凸阵类型时，所述根据所述km(0)确定初始化的第一迭代变量Am(0)、初始化的第二迭代变量am(0)、初始化的第一中间变量Bm(0)和初始化的第二中间变量bm(0)，包括：
根据公式对阵元m的延时变量Am(n)进行初始化；
根据公式am(0)＝2km(0)+1对阵元m的延时变量am(n)进行初始化；
根据公式对阵元m的延时变量Bm(n)进行初始化；
根据公式bm(0)＝2av(m)sin(θ+γ(m))/d+1对阵元m的延时变量bm(n)进行初始化。
8.根据权利要求1至4任一项所述的方法，其特征在于，所述初始化阵元m的延时变量km(n)，获取阵元m的km(0)，包括：
根据预设的数据表和聚焦点的偏转角度，获取与所述偏转角度对应的阵元m的km(0)；
其中，所述预设的数据表中包括不同的阵元位置和偏转角度对应的初始化的值。
9.根据权利要求1至4任一项所述的方法，其特征在于，所述方法还包括：
根据中心线与法线之间的夹角范围确定有效阵元的范围；
根据所述有效阵元的范围确定m的取值范围。
10.根据权利要求9所述的方法，其特征在于，根据中心线与法线之间的夹角范围确定有效阵元的范围，包括：
当探头的类型信息为线性类型时，根据公式计算根据公式计算所述有效阵元的范围为
当探头的类型信息为凸阵类型时，根据公式计算根据公式计算所述有效阵元的范围为
其中，θ为所述中心线与法线之间的夹角。
11.一种超声信号的处理装置，其特征在于，包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序，以实现如权利要求1至10中任一项所述的方法。
12.一种计算机存储介质，其特征在于，包括：所述计算机存储介质用于存储计算机程序，所述计算机程序执行时用于实现如权利要求1至10中任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810632043.XA CN110613476B (zh) | 2018-06-19 | 2018-06-19 | 超声信号的处理方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810632043.XA CN110613476B (zh) | 2018-06-19 | 2018-06-19 | 超声信号的处理方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110613476A true CN110613476A (zh) | 2019-12-27 |
| CN110613476B CN110613476B (zh) | 2022-09-23 |
Family
ID=68920240
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810632043.XA Active CN110613476B (zh) | 2018-06-19 | 2018-06-19 | 超声信号的处理方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110613476B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113925528A (zh) * | 2020-06-29 | 2022-01-14 | 青岛海信医疗设备股份有限公司 | 多普勒成像方法和超声设备 |
Citations (32)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3048031A (en) * | 1956-02-17 | 1962-08-07 | Siderurgie Fse Inst Rech | Device for automatically correcting the variables in the ultrasonic testing of materials |
| KR19990039095A (ko) * | 1997-11-10 | 1999-06-05 | 이민화 | 초음파영상화시스템의 초음파신호집속방법 및장치 |
| US6141293A (en) * | 1997-10-30 | 2000-10-31 | Netmor Ltd. | Ultrasonic positioning and tracking system |
| CN1412551A (zh) * | 2001-10-17 | 2003-04-23 | 中国石油天然气管道科学研究院 | 相控阵超声波仪器及其检测方法 |
| CN1788685A (zh) * | 2004-12-15 | 2006-06-21 | 深圳迈瑞生物医疗电子股份有限公司 | 基于双波束及合成孔径的接收方法及其装置 |
| CN101097256A (zh) * | 2006-06-28 | 2008-01-02 | 深圳迈瑞生物医疗电子股份有限公司 | 超声波束合成聚焦参数实时计算方法及装置 |
| CN101264022A (zh) * | 2007-03-16 | 2008-09-17 | 通用电气公司 | 用于在超声成像中精确时间延迟估计的方法和系统 |
| CN101371789A (zh) * | 2007-08-22 | 2009-02-25 | 深圳迈瑞生物医疗电子股份有限公司 | 实时计算波束合成聚焦延迟参数的方法及装置 |
| KR20090077104A (ko) * | 2008-01-10 | 2009-07-15 | 주식회사 메디슨 | 초음파 시스템 |
| CN101612049A (zh) * | 2008-06-26 | 2009-12-30 | 株式会社东芝 | 超声波诊断装置 |
| CN101683274A (zh) * | 2008-09-25 | 2010-03-31 | 深圳迈瑞生物医疗电子股份有限公司 | 偏转凸阵数字扫描变换和凸阵空间复合的方法及成像系统 |
| CN101785684A (zh) * | 2010-01-27 | 2010-07-28 | 天津大学 | 超声内窥镜合成孔径成像系统及较大孔径的合成方法 |
| CN101940479A (zh) * | 2010-07-15 | 2011-01-12 | 重庆大学 | 基于分段动态变迹系数的超声成像波束合成方法及其装置 |
| US20110237950A1 (en) * | 2010-03-23 | 2011-09-29 | Shenzhen Landwind Industry Co., Ltd. | Method and apparatus for multi-beam beamformer based on real-time calculation of time delay and pipeline design |
| CN102657918A (zh) * | 2012-05-07 | 2012-09-12 | 中国科学院声学研究所 | 一种聚焦方法和装置 |
| CN102809610A (zh) * | 2012-06-04 | 2012-12-05 | 北京航空航天大学 | 一种基于改进的动态深度聚焦的相控阵超声检测方法 |
| US20130211254A1 (en) * | 2010-05-28 | 2013-08-15 | Norwegian University Of Science And Technology (Ntnu) | Ultrasound acquisition |
| CN103310472A (zh) * | 2013-06-21 | 2013-09-18 | 中国科学院自动化研究所 | 基于正则化迭代的有限角度光声成像重建方法及装置 |
| CN104306023A (zh) * | 2014-10-24 | 2015-01-28 | 西安电子科技大学 | 基于压缩感知的超声成像快速实现方法 |
| US20150032005A1 (en) * | 2011-10-28 | 2015-01-29 | Beijing East Whale Image Technology Co., Ltd. | Color ultrasound system and method and device thereof for obtaining beam-forming data |
| CN104546003A (zh) * | 2015-01-27 | 2015-04-29 | 哈尔滨工业大学 | 一种基于平面波发射的可变帧率彩色超高速超声成像方法 |
| CN104605888A (zh) * | 2015-01-22 | 2015-05-13 | 飞依诺科技(苏州)有限公司 | 超声信号的时延估计方法及系统 |
| CN104739449A (zh) * | 2015-04-15 | 2015-07-01 | 惠尔图像公司 | 超声波束合成方法 |
| CN105520752A (zh) * | 2014-10-21 | 2016-04-27 | 爱飞纽医疗机械贸易有限公司 | 波束成形设备、超声成像设备以及波束成形方法 |
| KR20160054303A (ko) * | 2014-11-06 | 2016-05-16 | 삼성메디슨 주식회사 | 초음파 진단 장치 및 초음파 진단 방법 |
| CN105722462A (zh) * | 2013-09-13 | 2016-06-29 | 决策科学国际公司 | 在合成孔径图像形成中的相干展频 |
| CN105832366A (zh) * | 2016-03-18 | 2016-08-10 | 深圳安盛生物医疗技术有限公司 | 一种用于波束合成过程中的延时实时计算方法 |
| KR20170088534A (ko) * | 2016-01-25 | 2017-08-02 | 삼성메디슨 주식회사 | 빔포밍 장치, 이를 포함하는 프로브, 초음파 진단 장치 및 그 제어방법 |
| WO2017202362A1 (zh) * | 2016-05-25 | 2017-11-30 | 深圳开立生物医疗科技股份有限公司 | 一种超声成像方法、系统及装置 |
| US20170364658A1 (en) * | 2013-01-15 | 2017-12-21 | Cathworks Ltd | Vascular flow assessment |
| CN107970042A (zh) * | 2018-01-03 | 2018-05-01 | 声泰特（成都）科技有限公司 | 一种超声非线性定量实时成像方法与系统 |
| CN107997784A (zh) * | 2017-11-29 | 2018-05-08 | 声泰特（成都）科技有限公司 | 一种基于声速自适应修正的超声波束合成方法和系统 |
-
2018
- 2018-06-19 CN CN201810632043.XA patent/CN110613476B/zh active Active
Patent Citations (32)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3048031A (en) * | 1956-02-17 | 1962-08-07 | Siderurgie Fse Inst Rech | Device for automatically correcting the variables in the ultrasonic testing of materials |
| US6141293A (en) * | 1997-10-30 | 2000-10-31 | Netmor Ltd. | Ultrasonic positioning and tracking system |
| KR19990039095A (ko) * | 1997-11-10 | 1999-06-05 | 이민화 | 초음파영상화시스템의 초음파신호집속방법 및장치 |
| CN1412551A (zh) * | 2001-10-17 | 2003-04-23 | 中国石油天然气管道科学研究院 | 相控阵超声波仪器及其检测方法 |
| CN1788685A (zh) * | 2004-12-15 | 2006-06-21 | 深圳迈瑞生物医疗电子股份有限公司 | 基于双波束及合成孔径的接收方法及其装置 |
| CN101097256A (zh) * | 2006-06-28 | 2008-01-02 | 深圳迈瑞生物医疗电子股份有限公司 | 超声波束合成聚焦参数实时计算方法及装置 |
| CN101264022A (zh) * | 2007-03-16 | 2008-09-17 | 通用电气公司 | 用于在超声成像中精确时间延迟估计的方法和系统 |
| CN101371789A (zh) * | 2007-08-22 | 2009-02-25 | 深圳迈瑞生物医疗电子股份有限公司 | 实时计算波束合成聚焦延迟参数的方法及装置 |
| KR20090077104A (ko) * | 2008-01-10 | 2009-07-15 | 주식회사 메디슨 | 초음파 시스템 |
| CN101612049A (zh) * | 2008-06-26 | 2009-12-30 | 株式会社东芝 | 超声波诊断装置 |
| CN101683274A (zh) * | 2008-09-25 | 2010-03-31 | 深圳迈瑞生物医疗电子股份有限公司 | 偏转凸阵数字扫描变换和凸阵空间复合的方法及成像系统 |
| CN101785684A (zh) * | 2010-01-27 | 2010-07-28 | 天津大学 | 超声内窥镜合成孔径成像系统及较大孔径的合成方法 |
| US20110237950A1 (en) * | 2010-03-23 | 2011-09-29 | Shenzhen Landwind Industry Co., Ltd. | Method and apparatus for multi-beam beamformer based on real-time calculation of time delay and pipeline design |
| US20130211254A1 (en) * | 2010-05-28 | 2013-08-15 | Norwegian University Of Science And Technology (Ntnu) | Ultrasound acquisition |
| CN101940479A (zh) * | 2010-07-15 | 2011-01-12 | 重庆大学 | 基于分段动态变迹系数的超声成像波束合成方法及其装置 |
| US20150032005A1 (en) * | 2011-10-28 | 2015-01-29 | Beijing East Whale Image Technology Co., Ltd. | Color ultrasound system and method and device thereof for obtaining beam-forming data |
| CN102657918A (zh) * | 2012-05-07 | 2012-09-12 | 中国科学院声学研究所 | 一种聚焦方法和装置 |
| CN102809610A (zh) * | 2012-06-04 | 2012-12-05 | 北京航空航天大学 | 一种基于改进的动态深度聚焦的相控阵超声检测方法 |
| US20170364658A1 (en) * | 2013-01-15 | 2017-12-21 | Cathworks Ltd | Vascular flow assessment |
| CN103310472A (zh) * | 2013-06-21 | 2013-09-18 | 中国科学院自动化研究所 | 基于正则化迭代的有限角度光声成像重建方法及装置 |
| CN105722462A (zh) * | 2013-09-13 | 2016-06-29 | 决策科学国际公司 | 在合成孔径图像形成中的相干展频 |
| CN105520752A (zh) * | 2014-10-21 | 2016-04-27 | 爱飞纽医疗机械贸易有限公司 | 波束成形设备、超声成像设备以及波束成形方法 |
| CN104306023A (zh) * | 2014-10-24 | 2015-01-28 | 西安电子科技大学 | 基于压缩感知的超声成像快速实现方法 |
| KR20160054303A (ko) * | 2014-11-06 | 2016-05-16 | 삼성메디슨 주식회사 | 초음파 진단 장치 및 초음파 진단 방법 |
| CN104605888A (zh) * | 2015-01-22 | 2015-05-13 | 飞依诺科技(苏州)有限公司 | 超声信号的时延估计方法及系统 |
| CN104546003A (zh) * | 2015-01-27 | 2015-04-29 | 哈尔滨工业大学 | 一种基于平面波发射的可变帧率彩色超高速超声成像方法 |
| CN104739449A (zh) * | 2015-04-15 | 2015-07-01 | 惠尔图像公司 | 超声波束合成方法 |
| KR20170088534A (ko) * | 2016-01-25 | 2017-08-02 | 삼성메디슨 주식회사 | 빔포밍 장치, 이를 포함하는 프로브, 초음파 진단 장치 및 그 제어방법 |
| CN105832366A (zh) * | 2016-03-18 | 2016-08-10 | 深圳安盛生物医疗技术有限公司 | 一种用于波束合成过程中的延时实时计算方法 |
| WO2017202362A1 (zh) * | 2016-05-25 | 2017-11-30 | 深圳开立生物医疗科技股份有限公司 | 一种超声成像方法、系统及装置 |
| CN107997784A (zh) * | 2017-11-29 | 2018-05-08 | 声泰特（成都）科技有限公司 | 一种基于声速自适应修正的超声波束合成方法和系统 |
| CN107970042A (zh) * | 2018-01-03 | 2018-05-01 | 声泰特（成都）科技有限公司 | 一种超声非线性定量实时成像方法与系统 |
Non-Patent Citations (3)
| Title |
|---|
| 张永宏等: "超声相控线阵探头参数对波束形成的影响研究", 《仪表技术与传感器》 * |
| 王平等: "超声成像中基于虚拟阵元的双聚焦波束合成方法", 《重庆大学学报》 * |
| 齐雁等: "基于FPGA的医学超声成像数字波束合成器设计", 《计算机测量与控制》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113925528A (zh) * | 2020-06-29 | 2022-01-14 | 青岛海信医疗设备股份有限公司 | 多普勒成像方法和超声设备 |
| CN113925528B (zh) * | 2020-06-29 | 2023-11-03 | 青岛海信医疗设备股份有限公司 | 多普勒成像方法和超声设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110613476B (zh) | 2022-09-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Asen et al. | Implementing capon beamforming on a GPU for real-time cardiac ultrasound imaging | |
| CN101116622B (zh) | 波束合成的接收变迹参数的实时计算方法及其装置 | |
| US4989143A (en) | Adaptive coherent energy beam formation using iterative phase conjugation | |
| US8968199B2 (en) | Spatial compound imaging in an ultrasound system | |
| Hansen et al. | An object-oriented multi-threaded software beamformation toolbox | |
| JP5497821B2 (ja) | 流体流速検出装置及びプログラム | |
| JP2009532089A (ja) | 超音波受信ビームフォーマ用の遅延コントローラ | |
| WO2019127621A1 (zh) | 超声成像方法、系统和设备 | |
| JP2012170826A (ja) | 超音波ビーム形成方法及びその装置と超音波画像システム並びにコンピュータで読み取り可能な記録媒体 | |
| US20230324273A1 (en) | Method and device for quantifying viscoelasticity of a medium | |
| JP2000166925A (ja) | 信号集束遅延方法 | |
| JP4107840B2 (ja) | 機敏型ビーム形成器を使用した屈折遅延誤差補正 | |
| CN110507355B (zh) | 一种超声成像系统、方法、设备及介质 | |
| Choo et al. | Compressive spherical beamforming for localization of incipient tip vortex cavitation | |
| Bae et al. | Ultrasonic sector imaging using plane wave synthetic focusing with a convex array transducer | |
| CN107205729B (zh) | 超声波成像装置以及超声波信号的处理方法 | |
| CN110613476B (zh) | 超声信号的处理方法和装置 | |
| CN110731795B (zh) | 空间复合成像的处理方法和装置 | |
| CN112998745A (zh) | 一种用于超声成像的发射波束形成方法、系统及诊断设备 | |
| JP2017000547A (ja) | 超音波診断装置 | |
| CN102871685A (zh) | 超声探头几何参数的校正方法和装置及系统 | |
| US20190046162A1 (en) | Ultrasonic signal processor, ultrasonic diagnostic device, and ultrasonic signal processing method | |
| Yan et al. | A submatrix spatial coherence approach to minimum variance beamforming combined with sign coherence factor for coherent plane wave compounding | |
| Zhong et al. | A modified convolution model for calculating the far field directivity of a parametric array loudspeaker | |
| Moshavegh et al. | Optimization of synthetic aperture image quality |
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 |