CN107320129A - 基于医学的超声波信号的延时参数实时计算方法及装置 - Google Patents

基于医学的超声波信号的延时参数实时计算方法及装置 Download PDF

Info

Publication number
CN107320129A
CN107320129A CN201710422627.XA CN201710422627A CN107320129A CN 107320129 A CN107320129 A CN 107320129A CN 201710422627 A CN201710422627 A CN 201710422627A CN 107320129 A CN107320129 A CN 107320129A
Authority
CN
China
Prior art keywords
duration
array element
signal
focus
mrow
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
Application number
CN201710422627.XA
Other languages
English (en)
Other versions
CN107320129B (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.)
Qingdao Hisense Medical Equipment Co Ltd
Original Assignee
Qingdao Hisense Medical Equipment 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 Qingdao Hisense Medical Equipment Co Ltd filed Critical Qingdao Hisense Medical Equipment Co Ltd
Priority to CN201710422627.XA priority Critical patent/CN107320129B/zh
Publication of CN107320129A publication Critical patent/CN107320129A/zh
Application granted granted Critical
Publication of CN107320129B publication Critical patent/CN107320129B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves

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)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明实施例提供了一种基于医学的超声波信号的延时参数实时计算方法及装置。当需要得到信号分别从每一个阵元发射至该焦点所需耗费的时长时,可以按照本发明实施例的方法实时计算得到信号分别从每一个阵元发射至该焦点所需耗费的时长。在计算信号从某一阵元发射至该焦点所需耗费的时长时,需要利用到位置相邻的两个阵元之间的间隔距离、该阵元与中心线之间包括的间隔距离的数量、信号的传输速度以及阵元发射信号的发射频率。也即,存储与阵元的数量相同的间隔距离的数量,1个间隔距离,1个传输速度以及1个发射频率而已,存储的数据的数据量远远小于现有技术中需要存储的数据的数据量,进而可以节省存储资源。

Description

基于医学的超声波信号的延时参数实时计算方法及装置
技术领域
本发明实施例涉及医学超声技术领域,尤其涉及一种基于医学的超声波信号的延时参数实时计算方法及装置。
背景技术
当前,在医学超声诊断系统中,超声波信号从各个阵元发射并在焦点聚焦,由于超声波信号的聚焦特性,需要使得各个阵元发射的超声波信号同时达到焦点,但是,各个阵元与焦点之间的距离不同,使得超声波信号分别从每一个阵元发射后达到焦点所需耗费的时间不一样,因此,需要控制每一个阵元发射超声波信号的发射时刻,例如,距离焦点较远的阵元先发射超声波信号,距离焦点较近的阵元后发射超声波信号。
对于任意一个焦点,为了精确地控制每一个阵元发射信号的发射时刻,以使从每一个阵元分别发射的超声波信号能够同时到达该焦点,需要事先计算超声波信号分别从每一个阵元发射至该焦点所需耗费的时长,并存储在本地。对于其他每一个焦点,同样如此。
之后在需要向任意一个焦点发射超声波时,可以计算每相邻的两个阵元发射的信号至该焦点所需耗费的时长之间的时间差值,并根据得到的时间差值来确定每一个阵元发射超声波信号的发射时刻,进而每个阵元根据各自对应的发射时刻发射超声波。
然而,发明人发现,任意一个阵元与任意一个焦点都对应一个时长,如此,存储的时长的数量就等于阵元的数量与焦点的数量的乘积,可见存储的时长的数量较大,会占用过多的存储资源,造成存储资源浪费。
发明内容
为克服相关技术中存在的问题,本发明实施例提供一种基于医学的超声波信号的延时参数实时计算方法及装置。
根据本发明实施例的第一方面,提供一种基于医学的超声波信号的延时参数实时计算方法,所述方法包括:
获取预先存储的信号从第一阵元发射至焦点所需耗费的第一时长;
获取位置相邻的两个阵元之间的间隔距离、所述第一阵元与中心线之间包括的间隔距离的数量、所述信号的传输速度以及阵元发射信号的发射频率;
根据所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,获取信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值,所述第二阵元的位置与所述第一阵元的位置相邻;
根据所述第一时长与所述时间差值计算所述第二时长。
其中,所述根据所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,获取信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值,具体为:
若所述第二阵元与所述中心线之间的距离大于所述第一阵元与所述中心线之间的距离,利用所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,按照如下公式,计算信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值;
其中,在上述公式中,△t为所述时间差值,n为所述数量,s为所述间隔距离,c为所述传输速度,f为所述发射频率,t为所述第一时长,o为第一预设数值,p为第二预设数值,以及q为第三预设数值。
其中,所述根据所述第一时长与所述时间差值计算所述第二时长,具体为:
将所述第一时长与所述时间差值求和得到所述第二时长。
其中,所述根据所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,获取信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值,具体为:
若所述第二阵元与所述中心线之间的距离小于所述第一阵元与所述中心线之间的距离,利用所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,按照如下公式,计算信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值;
其中,在上述公式中,△t为所述时间差值,n为所述数量,s为所述间隔距离,c为所述传输速度,f为所述发射频率,t为所述第一时长,r为第四预设数值,w为第五预设数值,以及y为第六预设数值。
其中,所述根据所述第一时长与所述时间差值计算所述第二时长,具体为:
将所述第一时长与所述时间差值相减得到所述第二时长。
根据本发明实施例的第二方面,提供一种基于医学的超声波信号的延时参数实时计算装置,所述装置包括:
第一获取模块,用于获取预先存储的信号从第一阵元发射至焦点所需耗费的第一时长;
第二获取模块,用于获取位置相邻的两个阵元之间的间隔距离、所述第一阵元与中心线之间包括的间隔距离的数量、所述信号的传输速度以及阵元发射信号的发射频率;
第三获取模块,用于根据所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,获取信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值,所述第二阵元的位置与所述第一阵元的位置相邻;
计算模块,用于根据所述第一时长与所述时间差值计算所述第二时长。
其中,所述第三获取模块具体用于:若所述第二阵元与所述中心线之间的距离大于所述第一阵元与所述中心线之间的距离,利用所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,按照如下公式,计算信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值;
其中,在上述公式中,△t为所述时间差值,n为所述数量,s为所述间隔距离,c为所述传输速度,f为所述发射频率,t为所述第一时长,o为第一预设数值,p为第二预设数值,以及q为第三预设数值。
其中,所述计算模块具体用于:将所述第一时长与所述时间差值求和得到所述第二时长。
其中,所述第三获取模块具体用于:若所述第二阵元与所述中心线之间的距离小于所述第一阵元与所述中心线之间的距离,利用所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,按照如下公式,计算信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值;
其中,在上述公式中,△t为所述时间差值,n为所述数量,s为所述间隔距离,c为所述传输速度,f为所述发射频率,t为所述第一时长,r为第四预设数值,w为第五预设数值,以及y为第六预设数值。
其中,所述计算模块具体用于:将所述第一时长与所述时间差值相减得到所述第二时长。
本发明实施例提供的技术方案可以包括以下有益效果:
在本发明实施例中,对于任意一个焦点,无需事先计算信号分别从每一个阵元发射至该焦点所需耗费的时长,对于其他每一个焦点,同样如此。
对于任意一个焦点,当需要得到信号分别从每一个阵元发射至该焦点所需耗费的时长时,可以按照本发明实施例的方法实时计算得到信号分别从每一个阵元发射至该焦点所需耗费的时长。在计算信号从某一阵元发射至该焦点所需耗费的时长时,需要利用到位置相邻的两个阵元之间的间隔距离、该阵元与中心线之间包括的间隔距离的数量、信号的传输速度以及阵元发射信号的发射频率。即使有多个阵元,也只需事先存储每一个阵元与中心线之间包括的间隔距离的数量,也即,存储与阵元的数量相同的间隔距离的数量,1个间隔距离,1个传输速度以及1个发射频率而已,存储的数据的数据量远远小于现有技术中需要存储的数据的数据量。可见,相比于现有技术,本发明实施例可以节省存储资源。
例如,假设有10个阵元以及中心线上分布有10个焦点,现有技术中需要存储100个时长,一共100条数据。
本发明实施例只需存储10个间隔距离的数量、1个间隔距离,1个传输速度以及1个发射频率而已,当需要得到信号分别从每一个阵元发射至焦点所需耗费的时长时,一共需要使用13条数据,对于某一焦点,每计算出信号从一个阵元发射至该焦点所需耗费的时长就会将其存储,一共存储10个时长,当阵元发射完毕信号时,就会删除存储的10个时长。也即,本发明实施例中通常只会存储13条数据,最多会存储23条数据,远远小于现有技术中100条数据,从而相比于现有技术,本发明实施例可以节省存储资源。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本发明实施例。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本发明的实施例,并与说明书一起用于解释本发明实施例的原理。
图1是根据一示例性实施例示出的一种基于医学的超声波信号的延时参数实时计算方法的流程图;
图2是根据一示例性实施例示出的一种基于医学的超声波信号的延时参数实时计算装置的框图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明实施例相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本发明实施例的一些方面相一致的装置和方法的例子。
图1是根据一示例性实施例示出的一种基于医学的超声波信号的延时参数实时计算方法的流程图,如图1所示,该方法包括以下步骤。
在步骤S101中,获取预先存储的信号从第一阵元发射至焦点所需耗费的第一时长;
首先,需要说明的是,信号从第一阵元发射至焦点所需耗费的第一时长可以预先存储在本地的终端设备内,也可以存储在相应的云端或者服务器中,当本地终端设备需要时,从相应的云端或者服务器中,获取信号从第一阵元发射至焦点所需耗费的第一时长。
在本发明实施例中,中心线上分布有多个焦点。与中心线垂直的平面上分布有多个阵元,任意位置相邻的两个阵元之间的间隔距离都相同。
本发明实施例仅以多个阵元中的任意两个位置相邻的阵元进行举例,但不作为对本发明实施例保护范围的限制。在本发明实施例中,位置相邻的两个阵元分别为第一阵元和第二阵元。
在本发明实施例中,对于任意焦点,每当计算得到信号从某一阵元发射至该焦点所需耗费的时长,就会将该阵元的阵元标识与该时长组成一条记录,并存储在该焦点对应的、阵元标识与时长之间的对应关系中。进一步地,在每一个阵元都向该焦点发射信号之后,可以从上述对应关系中删除该记录。
因此,在本步骤中,当要计算信号从第二阵元发射至焦点所需耗费的第二时长时,可以在该焦点对应的、阵元标识与时长之间的对应关系中查找与第一阵元的阵元标识相对应的时长,并作为信号从第一阵元发射至焦点所需耗费的第一时长,然后执行步骤S102。
其中,阵元的阵元标识可以为阵元的编号或者名称等等,本发明实施例对此不加以限定。
在步骤S102中,获取位置相邻的两个阵元之间的间隔距离、第一阵元与中心线之间包括的间隔距离的数量、信号的传输速度以及阵元发射信号的发射频率;
在本发明实施例中,由于位置相邻的两个阵元之间的间隔距离都相同,且,每一个阵元与中心线之间包括的间隔距离的数量都是固定的。因此,技术人员事先可以存储相邻两个阵元之间的间隔距离,例如,在本地或者云端存储相邻两个阵元之间的间隔距离。
以及,对于任意阵元,技术人员可以将该阵元的阵元标识,以及,该阵元与中心线之间包括的间隔距离的数量组成一条记录,并存储在预先存储的阵元标识与数量之间的对应关系中,对于其他每一阵元,同样执行上述操作。
因此,在本步骤中,可以直接获取已存储的位置相邻的两个阵元之间的间隔距离,以及在预先存储的阵元标识与数量之间的对应关系中,查找与第一阵元的阵元标识相对应的数量,并作为第一阵元与中心线之间包括的间隔距离的数量。
在步骤S103中,根据第一时长、位置相邻的两个阵元之间的间隔距离、第一阵元与中心线之间包括的间隔距离的数量、信号的传输速度以及阵元发射信号的发射频率,获取第一时长与信号从第二阵元发射至焦点所需耗费的第二时长之间的时间差值,第二阵元的位置与第一阵元的位置相邻;
其中,此处暂不对获取第一时长与第二时长之间的距离差值的具体获取方法做解释说明,具体可参见下文所做介绍。
在步骤S104中,根据第一时长与该时间差值计算信号从第二阵元发射至焦点所需耗费的第二时长。
在本发明实施例中,如果第二阵元与中心线之间的距离大于第一阵元与中心线之间的距离,则可以将第一时长与该距离差值求和得到第二时长。如果第二阵元与中心线之间的距离小于第一阵元与中心线之间的距离,则可以将第一时长与该距离差值相减得到第二时长。
在本发明实施例中,对于任意一个焦点,无需事先计算信号分别从每一个阵元发射至该焦点所需耗费的时长,对于其他每一个焦点,同样如此。
对于任意一个焦点,当需要得到信号分别从每一个阵元发射至该焦点所需耗费的时长时,可以按照本发明实施例的方法实时计算得到信号分别从每一个阵元发射至该焦点所需耗费的时长。在计算信号从某一阵元发射至该焦点所需耗费的时长时,需要利用到位置相邻的两个阵元之间的间隔距离、该阵元与中心线之间包括的间隔距离的数量、信号的传输速度以及阵元发射信号的发射频率。即使有多个阵元,也只需事先存储每一个阵元与中心线之间包括的间隔距离的数量,也即,存储与阵元的数量相同的间隔距离的数量,1个间隔距离,1个传输速度以及1个发射频率而已,存储的数据的数据量远远小于现有技术中需要存储的数据的数据量。可见,相比于现有技术,本发明实施例可以节省存储资源。
例如,假设有10个阵元以及中心线上分布有10个焦点,现有技术中需要存储100个时长,一共100条数据。
本发明实施例只需存储10个间隔距离的数量、1个间隔距离,1个传输速度以及1个发射频率而已,当需要得到信号分别从每一个阵元发射至焦点所需耗费的时长时,一共需要使用13条数据,对于某一焦点,每计算出信号从一个阵元发射至该焦点所需耗费的时长就会将其存储,一共存储10个时长,当阵元发射完毕信号时,就会删除存储的10个时长。也即,本发明实施例中通常只会存储13条数据,最多会存储23条数据,远远小于现有技术中100条数据,从而相比于现有技术,本发明实施例可以节省存储资源。
在本发明实施例中,第二阵元与中心线之间的距离可能大于第一阵元与中心线之间的距离,也可能小于第一阵元与中心线之间的距离。
如果第二阵元与中心线之间的距离大于第一阵元与中心线之间的距离,则可以利用第一时长、位置相邻的两个阵元之间的间隔距离、第一阵元与中心线之间包括的间隔距离的数量、信号的传输速度以及阵元发射信号的发射频率,按照如下公式(1)计算第一时长与信号从第二阵元发射至焦点所需耗费的第二时长之间的时间差值;
其中,在上述公式(1)中,△t为第一时长与信号从第二阵元发射至焦点所需耗费的第二时长之间的时间差值,n为第一阵元与中心线之间包括的间隔距离的数量,s为位置相邻的两个阵元之间的间隔距离,c为信号的传输速度,f为阵元发射信号的发射频率,t为第一时长,o为第一预设数值,p为第二预设数值,以及q为第三预设数值。
其中,在本发明一个实施例中,o可以为2,p可以为1,q可以为2,当然,o、p以及q也可以按照需求取其他数值,本发明实施例对此不加以限定。
相应地,由于第二阵元与中心线之间的距离大于第一阵元与中心线之间的距离,因此,第二阵元与焦点之间的距离会大于第一阵元与焦点之间的距离,进而信号从第二阵元发射至焦点所需耗费的时长大于信号从第一阵元发射至焦点所需耗费的时间,所以,在根据第一时长与该时间差值计算信号从第二阵元发射至焦点所需耗费的第二时长时,可以将第一时长与该时间差值求和得到第二时长。
其中,在本发明实施例中,对于位置相邻的两个阵元:阵元A和阵元B,假设焦点与中心线所在的阵元切面之间的距离为d,位置相邻的的两个阵元之间的间隔距离为s,阵元A与中心线之间包括的间隔距离的数量为n,阵元B与中心线之间包括的间隔距离的数量为n+1,c为信号的传输速度,f为阵元发射信号的发射频率,在本发明实施例中的每一个阵元发射信号的发射频率都相同,则阵元A与中心线之间的距离为n*s,阵元B与中心线之间的距离为(n+1)*s。
可得:阵元A与焦点之间的距离La为La等于信号从阵元A发射至焦点所需耗费的时长ta与c的乘积,也即:
阵元B与焦点之间的距离Lb为
Lb与La之间的距离差值
其中,端点为阵元A和焦点的线段、端点为阵元B和焦点的线段以及端点为阵元A和阵元B的线段组成三角形,端点为阵元A和阵元B的线段的长度为一个间隔距离s。由于三角形的两条边的长度之间的差值小于第三边的长度的特性,因此,0<△x<s,也即,△x为s的一部分,将△x利用替代,得到等式1:其中,K=1、2、3……N-1。
将等式1转换为等式2:
将等式2中的等号两边平方,得到等式3:
由于小于1,则远远小于1,所以远远小于1,因此,可以将等式3中的舍弃,对等式3化简之后,得到等式4:
将等式4转换为等式5:
由于等式5中的为阵元A与焦点之间的距离La;且为Lb与La之间的差值△x;因此,可以将等式5转换为等式6:
等式6用于计算Lb与La之间的距离差值;
由于c为信号的传输速度,以及f为阵元发射信号的发射频率,进一步还可以根据距离差值、信号的传播速度以及阵元发射信号的发射频率,获取信号从阵元A发射至焦点所需耗费的时长与信号从阵元B发射至焦点所需耗费的时长之间的时间差值的计算公式:
由于时间差值的计算公式中的La等于ta与c的乘积,因此,可以利用ta*c替代La,得到:
由于时间差值的计算公式中包含有信号从阵元A传输至焦点所需耗费的时长ta,因此,需要在得到ta之后,才能利用时间差值的计算公式计算出信号从阵元A发射至焦点所需耗费的时长与信号从阵元B发射至焦点所需耗费的时长之间的时间差值。
在得到信号从阵元B发射至焦点所需耗费的时长与信号从阵元A发射至焦点所需耗费的时长之间的时间差值△t’之后,就可以根据信号从阵元A发射至焦点所需耗费的时长ta和时间差值△t’计算得到信号从阵元B发射至焦点所需耗费的时长tb,例如,由于阵元B与中心线之间的距离大于阵元A与中心线之间的距离,因此,需要将ta与△t’相加得到tb。
可见,在得到信号从某一阵元发射至焦点所需耗费的时长之后,将该时长带入时间差值的计算公式,可以得到:信号从某一阵元发射至焦点所需耗费的时长与信号从另一阵元发射至焦点所需耗费的时长之间的时间差值,另一阵元的位置与该某一阵元的位置相邻,由于另一阵元的与中心线之间的距离大于该某一阵元与中心线之间的距离,进而可以将该时长与该时间差值相加得到信号从另一阵元发射至焦点所需耗费的时长。
其中,在现有技术中,在计算信号从某一阵元发射至焦点所需耗费的时长时,需要计算该阵元与中心线之间的距离的平方,然后计算焦点与中心线所在的阵元切面之间的距离的平方,计算这两个平方之和,并对这两个平方之和开方得到该阵元与焦点之间的距离,再根据该距离、信号的传输速度以及阵元发射信号的发射频率计算信号从该阵元发射至焦点所需耗费的时长。
然而,在本发明实施例中,当需要得到信号从某一阵元发射至焦点所需耗费的时长时,无需通过计算该某一阵元与中心线之间的距离的平方,计算焦点与中心线所在的阵元切面之间的距离的平方,计算这两个平方之和并对这两个平方之和开方得到该某一阵元与焦点之间的距离,再根据该距离、信号的传输速度以及阵元发射信号的发射频率计算信号从某一阵元发射至焦点所需耗费的时长。而是,根据已经计算出的信号从另一阵元发射至焦点所需耗费的时长、位置相邻的两个阵元之间的间隔距离、另一阵元与中心线之间包括的间隔距离的数量、信号的传输速度以及阵元发射信号的发射频率,利用上述公式(1)进行乘法以及除法运算得到信号从该某一阵元发射至焦点所需耗费的时长,整个计算过程中无需进行开方运算,只需进行乘法和除法运算,使得整个计算过程在FPGA中容易实现,且节省了FPGA的时钟资源。
在本发明又一实施例中,如果第二阵元与中心线之间的距离小于第一阵元与中心线之间的距离,则可以利用第一时长、位置相邻的两个阵元之间的间隔距离、第一阵元与中心线之间包括的间隔距离的数量、信号的传输速度以及阵元发射信号的发射频率,按照如下公式(2)计算第一时长与信号从第二阵元发射至焦点所需耗费的第二时长之间的时间差值;
其中,在上述公式(2)中,△t为第一时长与信号从第二阵元发射至焦点所需耗费的第二时长之间的时间差值,n为第一阵元与中心线之间包括的间隔距离的数量,s为位置相邻的两个阵元之间的间隔距离,c为信号的传输速度,f为阵元发射信号的发射频率,t为第一时长,r为第四预设数值,w为第五预设数值,以及y为第六预设数值。
其中,在本发明一个实施例中,r可以为2,w可以为1,y可以为2,当然,r、w以及y也可以按照需求取其他数值,本发明实施例对此不加以限定。
相应地,由于第二阵元与中心线之间的距离小于第一阵元与中心线之间的距离,因此,第二阵元与焦点之间的距离会小于第一阵元与焦点之间的距离,进而信号从第二阵元发射至焦点所需耗费的时长小于信号从第一阵元发射至焦点所需耗费的时间,所以,在根据第一时长与该时间差值计算信号从第二阵元发射至焦点所需耗费的第二时长时,可以将第一时长与该时间差值相减得到第二时长。
其中,在本发明实施例中,对于位置相邻的两个阵元:阵元A和阵元B,假设焦点与中心线所在的阵元切面之间的距离为d,位置相邻的的两个阵元之间的间隔距离为s,阵元A与中心线之间包括的间隔距离的数量为n,阵元B与中心线之间包括的间隔距离的数量为n-1,c为信号的传输速度,f为阵元发射信号的发射频率,在本发明实施例中的每一个阵元发射信号的发射频率都相同,则阵元A与中心线之间的距离为n*s,阵元B与中心线之间的距离为(n-1)*s。
可得:阵元A与焦点之间的距离La为La等于信号从阵元A发射至焦点所需耗费的时长ta与c的乘积,也即:
阵元B与焦点之间的距离Lb为
Lb与La之间的距离差值
其中,端点为阵元A和焦点的线段、端点为阵元B和焦点的线段以及端点为阵元A和阵元B的线段组成三角形,端点为阵元A和阵元B的线段的长度为一个间隔距离s。由于三角形的两条边的长度之间的差值小于第三边的长度的特性,因此,0<△x<s,也即,△x为s的一部分,将△x利用替代,得到等式1:其中,K=1、2、3……N-1。
将等式1转换为等式2:
将等式2中的等号两边平方,得到等式3:
由于小于1,则远远小于1,所以远远小于1,因此,可以将等式3中的舍弃,对等式3化简之后,得到等式4:
将等式4转换为等式5:
由于等式5中的为阵元A与焦点之间的距离La;且为Lb与La之间的差值△x;因此,可以将等式5转换为等式6:
等式6用于计算Lb与La之间的距离差值,
由于c为信号的传输速度,以及f为阵元发射信号的发射频率,进一步还可以根据距离差值、信号的传播速度以及阵元发射信号的发射频率,获取信号从阵元A发射至焦点所需耗费的时长与信号从阵元B发射至焦点所需耗费的时长之间的时间差值的计算公式:
由于时间差值的计算公式中的La等于ta与c的乘积,因此,可以利用ta*c替代La,得到:
由于时间差值的计算公式中包含有信号从阵元A传输至焦点所需耗费的时长ta,因此,需要在得到ta之后,才能利用时间差值的计算公式计算出信号从阵元B发射至焦点所需耗费的时长与信号从阵元A发射至焦点所需耗费的时长之间的时间差值。
在得到信号从阵元A发射至焦点所需耗费的时长与信号从阵元B发射至焦点所需耗费的时长之间的时间差值△t’之后,就可以根据信号从阵元A发射至焦点所需耗费的时长ta和时间差值△t’计算得到信号从阵元B发射至焦点所需耗费的时长tb,例如,由于阵元B与中心线之间的距离小于阵元A与中心线之间的距离,因此,需要将ta与△t’相减得到tb。
可见,在得到信号从某一阵元发射至焦点所需耗费的时长之后,将该时长带入时间差值的计算公式,可以得到:信号从某一阵元发射至焦点所需耗费的时长与信号从另一阵元发射至焦点所需耗费的时长之间的时间差值,另一阵元的位置与该某一阵元的位置相邻,由于另一阵元的与中心线之间的距离小于该某一阵元与中心线之间的距离,进而可以将该时长与该时间差值相减得到信号从另一阵元发射至焦点所需耗费的时长。
其中,在现有技术中,在计算信号从某一阵元发射至焦点所需耗费的时长时,需要计算该阵元与中心线之间的距离的平方,然后计算焦点与中心线所在的阵元切面之间的距离的平方,计算这两个平方之和,并对这两个平方之和开方得到该阵元与焦点之间的距离,再根据该距离、信号的传输速度以及阵元发射信号的发射频率计算信号从该阵元发射至焦点所需耗费的时长。
然而,在本发明实施例中,当需要得到信号从某一阵元发射至焦点所需耗费的时长时,无需通过计算该某一阵元与中心线之间的距离的平方,计算焦点与中心线所在的阵元切面之间的距离的平方,计算这两个平方之和并对这两个平方之和开方得到该某一阵元与焦点之间的距离,再根据该距离、信号的传输速度以及阵元发射信号的发射频率计算信号从某一阵元发射至焦点所需耗费的时长。而是,根据已经计算出的信号从另一阵元发射至焦点所需耗费的时长、位置相邻的两个阵元之间的间隔距离、另一阵元与中心线之间包括的间隔距离的数量、信号的传输速度以及阵元发射信号的发射频率,利用上述公式(2)进行乘法以及除法运算得到信号从该某一阵元发射至焦点所需耗费的时长,整个计算过程中无需进行开方运算,只需进行乘法和除法运算,使得整个计算过程在FPGA中容易实现,且节省了FPGA的时钟资源。
图2是根据一示例性实施例示出的一种基于医学的超声波信号的延时参数实时计算装置的框图。参照图2,该装置包括:
第一获取模块11,用于获取预先存储的信号从第一阵元发射至焦点所需耗费的第一时长;
第二获取模块12,用于获取位置相邻的两个阵元之间的间隔距离、所述第一阵元与中心线之间包括的间隔距离的数量、所述信号的传输速度以及阵元发射信号的发射频率;
第三获取模块13,用于根据所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,获取信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值,所述第二阵元的位置与所述第一阵元的位置相邻;
计算模块14,用于根据所述第一时长与所述时间差值计算所述第二时长。
其中,所述第三获取模块13具体用于:若所述第二阵元与所述中心线之间的距离大于所述第一阵元与所述中心线之间的距离,利用所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,按照如下公式,计算信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值;
其中,在上述公式中,△t为所述时间差值,n为所述数量,s为所述间隔距离,c为所述传输速度,f为所述发射频率,t为所述第一时长,o为第一预设数值,p为第二预设数值,以及q为第三预设数值。
其中,所述计算模块14具体用于:将所述第一时长与所述时间差值求和得到所述第二时长。
其中,所述第三获取模块13具体用于:若所述第二阵元与所述中心线之间的距离小于所述第一阵元与所述中心线之间的距离,利用所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,按照如下公式,计算信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值;
其中,在上述公式中,△t为所述时间差值,n为所述数量,s为所述间隔距离,c为所述传输速度,f为所述发射频率,t为所述第一时长,r为第四预设数值,w为第五预设数值,以及y为第六预设数值。
其中,所述计算模块14具体用于:将所述第一时长与所述时间差值相减得到所述第二时长。
在本发明实施例中,对于任意一个焦点,无需事先计算信号分别从每一个阵元发射至该焦点所需耗费的时长,对于其他每一个焦点,同样如此。
对于任意一个焦点,当需要得到信号分别从每一个阵元发射至该焦点所需耗费的时长时,可以按照本发明实施例的方法实时计算得到信号分别从每一个阵元发射至该焦点所需耗费的时长。在计算信号从某一阵元发射至该焦点所需耗费的时长时,需要利用到位置相邻的两个阵元之间的间隔距离、该阵元与中心线之间包括的间隔距离的数量、信号的传输速度以及阵元发射信号的发射频率。即使有多个阵元,也只需事先存储每一个阵元与中心线之间包括的间隔距离的数量,也即,存储与阵元的数量相同的间隔距离的数量,1个间隔距离,1个传输速度以及1个发射频率而已,存储的数据的数据量远远小于现有技术中需要存储的数据的数据量。可见,相比于现有技术,本发明实施例可以节省存储资源。
例如,假设有10个阵元以及中心线上分布有10个焦点,现有技术中需要存储100个时长,一共100条数据。
本发明实施例只需存储10个间隔距离的数量、1个间隔距离,1个传输速度以及1个发射频率而已,当需要得到信号分别从每一个阵元发射至焦点所需耗费的时长时,一共需要使用13条数据,对于某一焦点,每计算出信号从一个阵元发射至该焦点所需耗费的时长就会将其存储,一共存储10个时长,当阵元发射完毕信号时,就会删除存储的10个时长。也即,本发明实施例中通常只会存储13条数据,最多会存储23条数据,远远小于现有技术中100条数据,从而相比于现有技术,本发明实施例可以节省存储资源。
关于上述实施例中的装置,其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述,此处将不做详细阐述说明。
本领域技术人员在考虑说明书及实践这里公开的发明后,将容易想到本发明的其它实施方案。本申请旨在涵盖本发明实施例的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本发明实施例的一般性原理并包括本发明实施例未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本发明实施例的真正范围和精神由所附的权利要求指出。
应当理解的是,本发明实施例并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。本发明实施例的范围仅由所附的权利要求来限制。

Claims (10)

1.一种基于医学的超声波信号的延时参数实时计算方法,其特征在于,所述方法包括:
获取预先存储的信号从第一阵元发射至焦点所需耗费的第一时长;
获取位置相邻的两个阵元之间的间隔距离、所述第一阵元与中心线之间包括的间隔距离的数量、所述信号的传输速度以及阵元发射信号的发射频率;
根据所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,获取信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值,所述第二阵元的位置与所述第一阵元的位置相邻;
根据所述第一时长与所述时间差值计算所述第二时长。
2.根据权利要求1所述的方法,其特征在于,所述根据所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,获取信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值,具体为:
若所述第二阵元与所述中心线之间的距离大于所述第一阵元与所述中心线之间的距离,利用所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,按照如下公式,计算信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值;
<mrow> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <mi>o</mi> <mo>*</mo> <mi>n</mi> <mo>+</mo> <mi>p</mi> <mo>)</mo> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>*</mo> <mi>f</mi> </mrow> <mrow> <mi>q</mi> <mo>*</mo> <mi>t</mi> <mo>*</mo> <msup> <mi>c</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>;</mo> </mrow>
其中,在上述公式中,△t为所述时间差值,n为所述数量,s为所述间隔距离,c为所述传输速度,f为所述发射频率,t为所述第一时长,o为第一预设数值,p为第二预设数值,以及q为第三预设数值。
3.根据权利要求2所述的方法,其特征在于,所述根据所述第一时长与所述时间差值计算所述第二时长,具体为:
将所述第一时长与所述时间差值求和得到所述第二时长。
4.根据权利要求1所述的方法,其特征在于,所述根据所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,获取信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值,具体为:
若所述第二阵元与所述中心线之间的距离小于所述第一阵元与所述中心线之间的距离,利用所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,按照如下公式,计算信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值;
<mrow> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <mi>r</mi> <mi>n</mi> <mo>-</mo> <mi>w</mi> <mo>)</mo> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>*</mo> <mi>f</mi> </mrow> <mrow> <mi>y</mi> <mo>*</mo> <mi>t</mi> <mo>*</mo> <msup> <mi>c</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>;</mo> </mrow>
其中,在上述公式中,△t为所述时间差值,n为所述数量,s为所述间隔距离,c为所述传输速度,f为所述发射频率,t为所述第一时长,r为第四预设数值,w为第五预设数值,以及y为第六预设数值。
5.根据权利要求4所述的方法,其特征在于,所述根据所述第一时长与所述时间差值计算所述第二时长,具体为:
将所述第一时长与所述时间差值相减得到所述第二时长。
6.一种基于医学的超声波信号的延时参数实时计算装置,其特征在于,所述装置包括:
第一获取模块,用于获取预先存储的信号从第一阵元发射至焦点所需耗费的第一时长;
第二获取模块,用于获取位置相邻的两个阵元之间的间隔距离、所述第一阵元与中心线之间包括的间隔距离的数量、所述信号的传输速度以及阵元发射信号的发射频率;
第三获取模块,用于根据所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,获取信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值,所述第二阵元的位置与所述第一阵元的位置相邻;
计算模块,用于根据所述第一时长与所述时间差值计算所述第二时长。
7.根据权利要求6所述的装置,其特征在于,所述第三获取模块具体用于:
若所述第二阵元与所述中心线之间的距离大于所述第一阵元与所述中心线之间的距离,利用所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,按照如下公式,计算信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值;
<mrow> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <mi>o</mi> <mo>*</mo> <mi>n</mi> <mo>+</mo> <mi>p</mi> <mo>)</mo> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>*</mo> <mi>f</mi> </mrow> <mrow> <mi>q</mi> <mo>*</mo> <mi>t</mi> <mo>*</mo> <msup> <mi>c</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>;</mo> </mrow>
其中,在上述公式中,△t为所述时间差值,n为所述数量,s为所述间隔距离,c为所述传输速度,f为所述发射频率,t为所述第一时长,o为第一预设数值,p为第二预设数值,以及q为第三预设数值。
8.根据权利要求7所述的装置,其特征在于,所述计算模块具体用于:将所述第一时长与所述时间差值求和得到所述第二时长。
9.根据权利要求6所述的装置,其特征在于,所述第三获取模块具体用于:
若所述第二阵元与所述中心线之间的距离小于所述第一阵元与所述中心线之间的距离,利用所述第一时长、所述间隔距离、所述数量、所述传输速度以及所述发射频率,按照如下公式,计算信号从第二阵元发射至焦点所需耗费的第二时长与所述第一时长之间的时间差值;
<mrow> <mi>&amp;Delta;</mi> <mi>t</mi> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <mi>r</mi> <mi>n</mi> <mo>-</mo> <mi>w</mi> <mo>)</mo> <msup> <mi>s</mi> <mn>2</mn> </msup> <mo>*</mo> <mi>f</mi> </mrow> <mrow> <mi>y</mi> <mo>*</mo> <mi>t</mi> <mo>*</mo> <msup> <mi>c</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>;</mo> </mrow>
其中,在上述公式中,△t为所述时间差值,n为所述数量,s为所述间隔距离,c为所述传输速度,f为所述发射频率,t为所述第一时长,r为第四预设数值,w为第五预设数值,以及y为第六预设数值。
10.根据权利要求9所述的装置,其特征在于,所述计算模块具体用于:将所述第一时长与所述时间差值相减得到所述第二时长。
CN201710422627.XA 2017-06-07 2017-06-07 基于医学的超声波信号的延时参数实时计算方法及装置 Active CN107320129B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710422627.XA CN107320129B (zh) 2017-06-07 2017-06-07 基于医学的超声波信号的延时参数实时计算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710422627.XA CN107320129B (zh) 2017-06-07 2017-06-07 基于医学的超声波信号的延时参数实时计算方法及装置

Publications (2)

Publication Number Publication Date
CN107320129A true CN107320129A (zh) 2017-11-07
CN107320129B CN107320129B (zh) 2020-06-16

Family

ID=60194131

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710422627.XA Active CN107320129B (zh) 2017-06-07 2017-06-07 基于医学的超声波信号的延时参数实时计算方法及装置

Country Status (1)

Country Link
CN (1) CN107320129B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110720947A (zh) * 2018-07-16 2020-01-24 青岛海信医疗设备股份有限公司 一种超声探头信号延时发射的控制方法及控制装置
CN110974304A (zh) * 2019-12-13 2020-04-10 山东大学齐鲁医院 基于穿戴式柔性超声换能器的超声波束合成系统及方法
CN114052786A (zh) * 2021-10-22 2022-02-18 武汉联影医疗科技有限公司 一种超声波发射方法和系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101097256A (zh) * 2006-06-28 2008-01-02 深圳迈瑞生物医疗电子股份有限公司 超声波束合成聚焦参数实时计算方法及装置
CN103913747A (zh) * 2012-12-31 2014-07-09 深圳市蓝韵实业有限公司 多通道扫描延时计算方法及装置
CN105832366A (zh) * 2016-03-18 2016-08-10 深圳安盛生物医疗技术有限公司 一种用于波束合成过程中的延时实时计算方法
CN106774901A (zh) * 2016-10-14 2017-05-31 厦门大学 基于超声定位的远距离pc体感输入方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101097256A (zh) * 2006-06-28 2008-01-02 深圳迈瑞生物医疗电子股份有限公司 超声波束合成聚焦参数实时计算方法及装置
CN103913747A (zh) * 2012-12-31 2014-07-09 深圳市蓝韵实业有限公司 多通道扫描延时计算方法及装置
CN105832366A (zh) * 2016-03-18 2016-08-10 深圳安盛生物医疗技术有限公司 一种用于波束合成过程中的延时实时计算方法
CN106774901A (zh) * 2016-10-14 2017-05-31 厦门大学 基于超声定位的远距离pc体感输入方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110720947A (zh) * 2018-07-16 2020-01-24 青岛海信医疗设备股份有限公司 一种超声探头信号延时发射的控制方法及控制装置
CN110974304A (zh) * 2019-12-13 2020-04-10 山东大学齐鲁医院 基于穿戴式柔性超声换能器的超声波束合成系统及方法
CN114052786A (zh) * 2021-10-22 2022-02-18 武汉联影医疗科技有限公司 一种超声波发射方法和系统

Also Published As

Publication number Publication date
CN107320129B (zh) 2020-06-16

Similar Documents

Publication Publication Date Title
CN107320129A (zh) 基于医学的超声波信号的延时参数实时计算方法及装置
CN105911551A (zh) 一种基于加权集合卡尔曼滤波算法的声速剖面反演方法
Abadi et al. Blind deconvolution for robust signal estimation and approximate source localization
CN106896403B (zh) 弹性高斯束偏移成像方法和系统
CN103197317A (zh) 基于fpga的sar成像方法
CN109270518B (zh) 一种水声垂直矢量阵列声呐的多目标模拟方法、装置及声呐
CN105832366B (zh) 一种用于波束合成过程中的延时实时计算方法
CN111950438B (zh) 基于深度学习的天宫二号成像高度计有效波高反演方法
CN104155030B (zh) 一种考虑声线弯曲的声学ct温度场重建方法
EP2202933A3 (en) Technique for generating an SC-FDMA signal
Diamantis et al. A comparison between temporal and subband minimum variance adaptive beamforming
CN106037795A (zh) 一种优化全波反演法重建声速图像的方法
JP2010127771A (ja) 合成開口ソーナー、合成開口ソーナーの位相誤差補正方法及びプログラム
JP2006524330A (ja) 多重ビームの形成方法および装置
CN111737823B (zh) 线性动力学系统运动状态预测方法及系统
CN108872902A (zh) 波形输出方法和装置
CN105116383B (zh) 产生fda雷达的发射波束的方法、装置及fda雷达
CN106814350A (zh) 基于压缩感知的外辐射源雷达参考信号信杂比估计方法
CN108513355B (zh) 基于单基站的网络定位方法、装置及设备
CN106680823A (zh) 利用抹香鲸叫声脉冲的目标距离和速度探测方法与装置
CN106501786A (zh) 一种基于矩阵相关的微动目标参数估计方法
CN103544134B (zh) 传递函数估计设备和方法
CN110489800A (zh) 一种基于矩阵正则化的结构动荷载稀疏识别方法
CN109951253A (zh) 一种数据报文首部校验和生成方法及装置
CN103674815A (zh) 一种基于传感器特征的气溶胶mtf监测方法

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