CN110507355A - 一种超声成像系统、方法、设备及介质 - Google Patents
一种超声成像系统、方法、设备及介质 Download PDFInfo
- Publication number
- CN110507355A CN110507355A CN201910893370.5A CN201910893370A CN110507355A CN 110507355 A CN110507355 A CN 110507355A CN 201910893370 A CN201910893370 A CN 201910893370A CN 110507355 A CN110507355 A CN 110507355A
- Authority
- CN
- China
- Prior art keywords
- receiving
- line
- amplitude
- group
- array element
- 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
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000001427 coherent effect Effects 0.000 claims abstract description 45
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 26
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 17
- 239000000523 sample Substances 0.000 claims description 42
- 230000006870 function Effects 0.000 claims description 24
- 230000005540 biological transmission Effects 0.000 claims description 23
- 230000035945 sensitivity Effects 0.000 claims description 18
- 238000012285 ultrasound imaging Methods 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 14
- 230000015654 memory Effects 0.000 claims description 10
- 210000001503 joint Anatomy 0.000 claims description 7
- 238000012986 modification Methods 0.000 claims description 6
- 230000004048 modification Effects 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 16
- 230000008569 process Effects 0.000 description 9
- 238000003384 imaging method Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 238000002604 ultrasonography Methods 0.000 description 4
- 238000013178 mathematical model Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
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
- 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)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种超声成像系统、方法、设备及介质,根据每一条第一接收线对应的第一接收点同一轴向上,与第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;确定每个组对应的相干系数权重值,根据所述每个组对应的相干系数权重值与该组包含的所述接收点对所述第一接收点的相干能量和的乘积的和,对所述第一接收线的幅值进行修正;合并修正后的每条第一接收线,生成对应的图像数据,并发送给所述图像显示单元生成对应的图像,提高了并行波束合成的帧频与准确性,解决了现有技术中接收线数据在不可信的情况下对修正结果造成误差的问题。
Description
技术领域
本发明涉及数字图像处理领域,尤其涉及一种超声成像系统、方法、设备及介质。
背景技术
在超声系统中,为了提供帧频,也就是时间分辨率,往往需要进行并行波束合成。即发射一条线,接收多条线来进行波束合成。但是因为发射线和接收线无法在几何位置上重合,所以往往会在图像中呈现一些问题。具体表现在图像中会出现较为明显的栅格。这是由于接收线与发射线没有重合,导致能量发散,采集回来的数据包含了较多的其他邻近线的信息所造成的。因此需要对接收线进行修正,这是并行波束合成的核心部分。
现有技术中为了实现并行波束合成,许多超声设备公司提出了自己的并行波束合成的实现方法。其中,有些并行波束合成的实现方法是依据某条接收线周围多条接收线的信息,对该条接收线进行修正,但如果其中某条线不可靠,会对修正产生影响。
现有技术中还有的并行波束合成的实现方法是利用回顾式相加的方法,即先对相位进行修正,再对相应的接收线进行相加,去除误差。该实现方法是比较容易实现,但是其对接收线的处理还是比较简单。因为传统的波束合成方案中,每次发射和接收都是同一条线,这样往往造成帧频较低的问题,在某些情况下无法满足实时性的要求。此外,由于发射是一个从探头到探测组织的物理过程,故发射一条线的时间是固定的,而接收是在FPGA或者PC中对回波数据处理的,是一个计算过程,也是可以很快的完成。所以,为了节省发射时间,提高帧频。可以发射N条线,接收回来M条线,M>N。当M/N=2时,则称该并行波束合成为2倍并行波束合成;当M/N=4时,则称该并行波束合成为4倍并行波束合成。但是,如果在4倍并行波束合成时采用上述的回顾式相加的方法,则有一半接收线是重叠的,此时4倍并行波束合成的帧频和2倍并行波束合成的帧频一致,因此只能达到2倍并行波束合成的帧频,无法再提高,所以采用回顾式相加的方法还是有很大的局限性。
发明内容
本发明实施例提供了一种超声成像系统、方法、设备及介质,用以解决现有技术中接收线在不可信的情况下对修正结果造成误差的问题,提高并行波束合成的帧频与准确性。
本发明实施例提供了一种超声成像系统,所述系统包括探头、发射波束单元、多线接收波束单元以及图像显示单元,还包括:多线波束合成单元;
所述多线波束合成单元与所述图像显示单元连接,用于针对每条第一接收线,确定该第一接收线对应的第一接收点同一轴向距离上,与所述第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;确定每个组对应的相干系数权重值,根据所述每个组对应的相干系数权重值与该组包含的所述接收点对所述第一接收点的相干能量和的乘积的和,对所述第一接收线的幅值进行修正;合并修正后的每条第一接收线,生成对应的图像数据,并发送给所述图像显示单元生成对应的图像。
进一步地,所述多线波束合成单元,还用于在确定第二接收点之前,根据预先确定的相位差和/或幅度差,对在发射聚焦、接收聚焦得到的信号的相位和/或幅度进行修正。
进一步地,所述多线波束合成单元在同时对相位和幅度进行修正时,具体用于根据Scorr(xt,zt,xr,zr)=Aejφ=A·ΔAejφ+Δφ,对在发射聚焦、接收聚焦得到的信号的相位和幅度进行修正,其中,
Δφ(t,xt,zt,xr,zr)=φ(t,xt,zr,xr,zr)-φ(t,xt,zt,xr,zr),
其中A为所述发射聚焦、接收聚焦得到的信号的幅值,A(t,xt,zr,xr,zr)为幅度的参考点,φ(t,xt,zr,xr,zr)为相位的参考点,
且
A(t,xt,zt,xr,zr)=abs(IQ(H(t,xt,zt,xr,zr)),
φ(t,xt,zt,xr,zr)=angle(IQ(H(t,xt,zt,xr,zr))),
其中,A(t,xt,zt,xr,zr)是在时刻t得到的信号的幅度,φ(t,xt,zt,xr,zr)是在时刻t得到的信号的相位,IQ代表正交解调,abs代表复数取模,angle代表复数取角度,
其中,xt是发射聚焦点的横向坐标,zt是发射聚焦点的纵向坐标,xr是接收聚焦点的横向坐标,zr是接收聚焦点的纵向坐标,t是时间,是探头发射时第i个阵元对坐标(xt,zt)的指向性函数,代表探头发射时第j个阵元对坐标(xr,zr)的指向性函数,ξi是探头发射时第i个阵元的敏感系数,ηj是探头发射时第j个阵元的敏感系数,Qi代表发射时第i个阵元的变迹系数,Bj代表发射时第j个阵元的变迹系数,P代表发射信号,Yi(x,z)代表第i个阵元距离坐标(x,z)的距离,c是声速,μ是声传播时的衰减因子,τi(x,z)代表第i个阵元在坐标(x,z)上的延迟。
进一步地,所述多线波束合成单元,具体用于根据
确定每个组对应的相干系数权重值,其中,a为并行波束合成倍数的二分之一,为每个组对应的可信度,为接收线的相干系数。
进一步地,所述多线波束合成单元,具体用于根据以下公式,确定第一接收线本身的能量,
其中,为接收线Rxjn所包含各条线数的总能量,bi,jn是计算接收线Rxjn本身的能量时,对接收线Rxjn所包含各条线数的总能量参与合成时的比例系数。
本发明实施例还提供了一种超声成像方法,所述方法包括:
针对每条第一接收线,确定该第一接收线对应的第一接收点同一轴向距离上,与所述第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;
确定每个组对应的相干系数权重值,根据所述每个组对应的相干系数权重值与该组包含的所述接收点对所述第一接收点的相干能量和的乘积的和,对所述第一接收线的幅值进行修正;
合并修正后的每条第一接收线,生成对应的图像数据,并发送给所述图像显示单元生成对应的图像。
进一步地,在确定第二接收点之前,根据预先确定的相位差和/或幅度差,对在发射聚焦、接收聚焦得到的信号的相位和/或幅度进行修正。
进一步地,在同时对相位和幅度进行修正时,根据Scorr(xt,zt,xr,zr)=Aejφ=A·ΔAejφ+Δφ,对在发射聚焦、接收聚焦得到的信号的相位和幅度进行修正,其中,
Δφ(t,xt,zt,xr,zr)=φ(t,xt,zr,xr,zr)-φ(t,xt,zt,xr,zr),
其中A为所述发射聚焦、接收聚焦得到的信号的幅值,A(t,xt,zr,xr,zr)为幅度的参考点,φ(t,xt,zr,xr,zr)为相位的参考点,
且
A(t,xt,zt,xr,zr)=abs(IQ(H(t,xt,zt,xr,zr)),
φ(t,xt,zt,xr,zr)=angle(IQ(H(t,xt,zt,xr,zr))),
其中,A(t,xt,zt,xr,zr)是在时刻t得到的信号的幅度,φ(t,xt,zt,xr,zr)是在时刻t得到的信号的相位,IQ代表正交解调,abs代表复数取模,angle代表复数取角度,
其中,xt是发射聚焦点的横向坐标,zt是发射聚焦点的纵向坐标,xr是接收聚焦点的横向坐标,zr是接收聚焦点的纵向坐标,t是时间,是探头发射时第i个阵元对坐标(xt,zt)的指向性函数,代表探头发射时第j个阵元对坐标(xr,zr)的指向性函数,ξi是探头发射时第i个阵元的敏感系数,ηj是探头发射时第j个阵元的敏感系数,Qi代表发射时第i个阵元的变迹系数,Bj代表发射时第j个阵元的变迹系数,P代表发射信号,Yi(x,z)代表第i个阵元距离坐标(x,z)的距离,c是声速,μ是声传播时的衰减因子,τi(x,z)代表第i个阵元在坐标(x,z)上的延迟。
进一步地,所述确定每个组对应的相干系数权重值,包括:
根据
确定每个组对应的相干系数权重值,其中,a为并行波束合成倍数的二分之一,为每个组对应的可信度,为接收线的相干系数。
进一步地,所述确定每个组对应的相干系数权重值之后,还包括:
根据以下公式,确定第一接收线本身的能量,
其中,为接收线Rxjn所包含各条线数的总能量,bi,jn是计算接收线Rxjn本身的能量时,对接收线Rxjn所包含各条线数的总能量参与合成时的比例系数。
本发明实施例还提供了一种超声成像设备,所述超声成像设备包括处理器,所述处理器用于执行存储器中存储的计算机程序时实现如上述任一所述方法的步骤。
本发明实施例还提供了一种计算机可读存储介质,其存储有计算机程序,所述计算机程序被处理器执行时实现如上述任一所述方法的步骤。
由于在本发明实施例中,根据每一条第一接收线对应的第一接收点同一轴向上,与第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;根据每个组对应的相干系数权重值与该组包含的接收点对第一接收点的相干能量和的乘积的和,对第一接收线的幅值进行修正,提高了并行波束合成的准确性,解决了现有技术中接收线数据在不可信的情况下对修正结果造成误差的问题。
附图说明
图1为本发明实施例提供的一种超声成像系统的结构示意图;
图2为本发明实施例提供的2倍并行波束合成示意图;
图3为本发明实施例提供的一种超声成像的过程示意图;
图4为本发明实施例还提供一种超声成像设备结构示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
为了避免接收线数据在不可信的情况下对修正结果造成的误差,提高并行波束合成的帧频与准确性,本发明实施例提供了一种超声成像系统、方法、设备及介质。
实施例1:
图1为本发明实施例提供的一种超声成像系统的结构示意图,该系统包括:包括探头10、发射波束单元11、多线接收波束单元12以及图像显示单元13,还包括:多线波束合成单元14;
所述多线波束合成单元14与所述图像显示单元13连接,用于针对每条第一接收线,确定该第一接收线对应的第一接收点同一轴向距离上,与所述第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;确定每个组对应的相干系数权重值,根据所述每个组对应的相干系数权重值与该组包含的所述接收点对所述第一接收点的相干能量和的乘积的和,对所述第一接收线的幅值进行修正;合并修正后的每条第一接收线,生成对应的图像数据,并发送给所述图像显示单元13生成对应的图像。
该超声成像系统的探头10,其中包括有阵列换能器,用于将电信号转换成对应的超声波或者将接收到的超声波转换成对应的电信号。
发射波束单元11耦合在该阵列换能器上,用于发射波束。
多线接收波束单元12耦合在该阵列换能器上,用于根据需要响应的波束产生对应的多条接收线。
其中,探头10、发射波束单元11、多线接收波束单元12的工作方式与现有技术相同,这里不再赘述。
第一接收线是接收线中的任意一条接收线。因为发射线和接收线在几何位置上不重合,且两者本身就会有扩散,所以发射线和接收线并非是一条直线,即发射发射线时,发射线的能量会扩散到其他线的位置,而得到接收线时,接收线会将其他线的能量接收过来。因此,每条接收线实际上包含了许多其他线的能量,并将这些其他线的能量理解成相干能量,在得到某一接收线时,就需要把这些相干能量删除。而根据预先确定的相位差和/或幅度差,对在发射聚焦、接收聚焦得到的信号的相位和/或幅度进行修正后得到的该接收线,是不能保证该接收线所包含的总能量都是可靠的,其中有一些接收线受到一些相干而不可信,这时就需要根据各个相干能量的相干系数确定该接收线可信的程度,对该接收线进一步的修正。
其中,一般对第一接收线上需要修正的接收点相干比较大的接收点,为该待修正的接收点相邻的同一轴向距离上的接收点,即第一接收线对应的接收点同一轴向距离上,与该第一接收点相邻的第二接收点,所以在对该第一接收线对应的第一接收点进行修正时,只需考虑与该第一接收点相邻的,且与该第一接收点在同一轴向距离的第二接收点对该第一接收点的影响即可。因此,在得到第一接收线之后,根据该第一接收线对应的第一接收点的横、纵坐标,确定与该第一接收线相邻的,且与该第一接收点在同一轴向距离上设定数量的第二接收点。
例如,当前第一接收线对应的第一接收点的坐标为(xi,z),设定数量为4,则确定该第一接收点同一轴向距离上的相邻的4个点作为第二接收点,即该第一接收点附近的四个第二接收点为(xi-2,z)、(xi-1,z)、(xi+1,z)、(xi+2,z)。
确定了设定数量的第二接收点所在的第二接收线之后,将该第一接收点与
第二接收点随机划分成多个组,每个组中都包含了第二设定数量的接收点,且每个组内包含的接收点均不相同。其中,第二设定数量小于等于第一设定数量。
假设进行4倍并行波束合成,第二设定数量为4,第一接收点为(xi,z),根据公式确定每个组对应的相干系数权重值。
其中, 为每个组对应的可信度,为接收点所在的接收线的相干系数。
而相干系数可以采用下面这种相干系数的算法:
其中,代表第m帧的接收线的回波信号强度,代表第m帧的图像中该点的具体位置,m代表是第m帧图像,M为包含有接收线的图像的帧数。
需要说明的是,相干系数的计算不限于上面的形式,其他的相干系数计算方式还有GCF、SCF、GSCF、STF等。
其中,
j1=i-2,i-1,......i+2
j2=i-2,i-1,......i+2
j3=i-2,i-1,......i+2
j4=i-2,i-1,......i+2
j1≠j2≠j3≠j4
j1、j2、j3、j4代表该组内的4个接收点,其值为该组内对应位置处接收点的位置标识值。
进一步的,该组包含的接收点对第一接收点的相干能量和是指每个组内确定的对应的每个接收点所包含的总能量以及该接收点所在的接收线本身的能量与第一接收线总能量的比例系数的乘积的和。
具体的,确定第一接收线修正的幅值的公式为:
其中,为接收线Rxjn所包含各条线数的总能量,是计算接收线本身的能量时,对接收线所包含各条线数的总能量所采取的一个比例系数,即不同的在参与合成时的比例系数。
需要说明的是,bi,jn与纵向距离r有关,需要在实际操作中通过调试确定,具体的相信本领域技术人员能够基于本发明实施例的描述确定具体的bi,jn确定过程,在此不再赘述。
按照上述的方式,对该第一接收线进行幅值上的修正之后,合并每条修正后的第一接收线,生成对应的图像数据,并发送给所述图像显示单元生成对应的图像。
图像显示单元接收到多线波束合成单元发送的图像数据后,生成对应的图像。
由于在本发明实施例中,根据每一条第一接收线对应的第一接收点同一轴向上,与第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;根据每个组对应的相干系数权重值与该组包含的接收点对第一接收点的相干能量和的乘积的和,对第一接收线的幅值进行修正,提高了并行波束合成的准确性,解决了现有技术中接收线数据在不可信的情况下对修正结果造成误差的问题。
实施例2:
为了提高并行波束合成的准确性,在上述实施例的基础上,在本发明实施例中,所述多线波束合成单元,还用于在确定第二接收点之前,根据预先确定的相位差和/或幅度差,对在发射聚焦、接收聚焦得到的信号的相位和/或幅度进行修正。
进行修正时,可以只对相位进行修正,也可以只对幅度进行修正,较佳地,是同时对相位和幅度进行修正。
图2为本发明实施例提供的一种2倍并行波束合成示意图。在并行波束合成中,因为发射线和接收线在几何上并不重合,所以接收的接收线和发射线往往会有相位差和幅度差,如图2中弧线所示。因此在确定第二接收点之前,必须要将这个相位差和/或幅度差去除。
多线波束合成单元在同时对相位和幅度进行修正时,具体用于根据Scorr(xt,zt,xr,zr)=Aejφ=A·ΔAejφ+Δφ,对在发射聚焦、接收聚焦得到的信号的相位和幅度进行修正,其中,
Δφ(t,xt,zt,xr,zr)=φ(t,xt,zr,xr,zr)-φ(t,xt,zt,xr,zr),
其中A为所述发射聚焦、接收聚焦得到的信号的幅值,A(t,xt,zt,xr,zr)为幅度的参考点,φ(t,xt,zr,xr,zr)为相位的参考点,
且
A(t,xt,zt,xr,zr)=abs(IQ(H(t,xt,zt,xr,zr)),
φ(t,xt,zt,xr,zr)=angle(IQ(H(t,xt,zt,xr,zr))),
其中,A(t,xt,zt,xr,zr)是在时刻t得到的信号的幅度,φ(t,xt,zt,xr,zr)是在时刻t得到的信号的相位,IQ代表正交解调,abs代表复数取模,angle代表复数取角度,
其中,xt是发射聚焦点的横向坐标,zt是发射聚焦点的纵向坐标,xr是接收聚焦点的横向坐标,zr是接收聚焦点的纵向坐标,t是时间,是探头发射时第i个阵元对坐标(xt,zt)的指向性函数,代表探头发射时第j个阵元对坐标(xr,zr)的指向性函数,ξi是探头发射时第i个阵元的敏感系数,ηj是探头发射时第j个阵元的敏感系数,Qi代表发射时第i个阵元的变迹系数,Bj代表发射时第j个阵元的变迹系数,P代表发射信号,Yi(x,z)代表第i个阵元距离坐标(x,z)的距离,c是声速,μ是声传播时的衰减因子,τi(x,z)代表第i个阵元在坐标(x,z)上的延迟。
在本发明实施例中,采用的声场仿真数学模型是:
根据该模型中所需要的参数在超声成像系统中采集对应的各项参数,通过仿真,确定H(t,xt,zt,xr,zr)后,对该H(t,xt,zt,xr,zr)进行正交解调,对正交解调的结果分别进行复数取模和和复数取角度,得到一个幅值A(t,xt,zt,xr,zr)和角度φ(t,xt,zt,xr,zr)。
进一步地,为了确定对在发射聚焦、接收聚焦得到的信号进行修正的相位差和幅度差,还需要选取一个点作为进行幅度和相位校正的参考点。在选取过程中,可以选取距离发射聚焦点在设定范围内的点,在本申请中选取(xt,zr)作为参考点,取A(t,xt,zr,xr,zr)作为幅度调整的参考点,φ(t,xt,zt,xr,zr)作为相位调整的参考点。
然后将上述得到的幅值A(t,xt,xt,xr,zr)、角度φ(t,xt,zr,xr,zr)、幅度的参考点A(t,xt,zr,xr,zr)、相位的参考点φ(t,xt,zr,xr,zr)代入下式,确定的相位差和幅度差,
Δφ(t,xt,zt,xr,zr)=φ(t,xt,zr,xr,zr)-φ(t,xt,zt,xr,zr),
由于考虑了发射和接收特性共同造成的相位差和/或幅度差,对发射聚焦、接收聚焦得到的信号进行相位和/或幅值上的修正,一定程度上提高了并行波束合成的准确性。
实施例3:
为了确定该第一接收线是否可靠,在上述各实施例的基础上,在本发明实施例中,所述多线波束合成单元,具体用于根据
确定每个组对应的相干系数权重值,其中,a为并行波束合成倍数的二分之一,为每个组对应的可信度,为接收线的相干系数。
这里为了便于说明该公式,采用具体的例子进行说明。
假设对第一接收点为(xi,z)进行修正,确定在其同一轴向距离上,与在该第一接收点附近的四个第二接收点为(xi-2,z)、(xi-1,z)、(xi+1,z)、(xi+2,z),加上第一接收点(xi,z)一共5个点。
进一步地,令
j1=i-2,i-1,......i+2
j2=i-2,i-1,......i+2
j3=i-2,i-1,......i+2
j4=i-2,i-1,......i+2
j1≠j2≠j3≠j4
这里j1、j2、j3、j4是在[xi-2,xi+2]内的任意一个整数,且互不相等,即当前i为3,则j1可能为[1,5]之间任意一个数,其值可能为1、2、3、4、5中的一个,若j1为1,则j2为剩下的四个数中的任意一个数,其值可能为2、3、4、5中的一个,同时以此类推j3、j4的值。按照上述方法确定j1、j2、j3、j4,则可以确定5种不同的j1、j2、j3、j4值的组合,即可以确定5种不同的4个接收点的组合。
具体的,根据当前j1、j2、j3、j4的值,确定对应的4个接收点,并可以通过相干系数的算法公式确定该4个接收点中每一个点所在的接收线的总能量与该4个接收点所在的接收线的总能量和之间的相干性,即确定和
例如,可以采用CF算法计算相干系数:
其中,代表第m帧的接收线的回波信号强度,代表第m帧的图像中该点的具体位置,m代表是第m帧图像,M为包含有接收线的图像的帧数。
本发明实施例中相干系数是一种衡量接收线的回波信号强度相干性的指标,其定义为接收线中相干能量与接收线中包含有相干能量和非相干能量的总能量的比值。相干系数代表了第一接收线的一种可信程度,相干系数越高,说明第一接收线越可信。
需要说明的是,相干系数的计算不限于上面的形式,其他的相干系数计算方式还有GCF、SCF、GSCF、STF等。
得到上述4个接收点对应的相干系数后进行相乘,确定当前j1、j2、j3、j4对应的即
由于有5种不同的4个接收点之间的相干系数组合,则可以确定五个对应不同的4个接收点之间的相干系数组合的
其中,当越小,这4个接收点组合的信息的可信度越低,其对的影响力也越低,甚至完全没有影响力;越大,这4个接收点组合的信息的可信度越高,其对的影响力也越高,所以需要对得到的五个进行加权,即
是对应当前的进行加权得到的。
需要说明的是,也可以是另外的加权形式,并不限于上面的形式。
由于将接收点划分为多个组,其中每个组中包含第二设定数量的接收点,每个组中包含的所述接收点完全不同,并确定每个组对应的相干系数权重值,帮助更好的确定第一接收线是否可靠,从而保证后续对该第一接收线修正的准确性。
实施例4:
为了进一步提高并行波束合成的帧频与准确性,在上述各实施例的基础上,在本发明实施例中,所述多线波束合成单元,具体用于根据以下公式,确定第一接收线本身的能量,
其中,为接收线Rxjn所包含各条线数的总能量,是计算接收线本身的能量时,对接收线所包含各条线数的总能量所采取的一个比例系数,即不同的在参与合成时的比例系数。
根据确定该接收线所包含各条接收线的总能量
需要说明的是,需要在实际情况中确定,其实际与为纵向距离r有关的函数,需要在实际操作中通过调试确定,具体的相信本领域技术人员能够基于本发明实施例的描述确定具体的确定过程,在此不再赘述。
具体的,假设第二设定数量为4,有5种不同的包含有4个接收点的组合,确定每个组对应的相干系数权重值与该组包含的4个接收点对所述第一接收点的相干能量和的乘积的和,确定该第一接收线修正后的幅值,即
其中,就是第一接收点(xi,z)修正后的幅值。
下面通过具体的实施例,对本发明实施例提供的超声成像系统进行详细说明。
为了方便说明,以4倍并行波束合成为例,并以第一接收线Rx3进行具体说明。
当接收到第一接收线Rx3后,根据本发明实施例采用的声场仿真模型:
其中,xt是发射聚焦点的横向坐标,zt是发射聚焦点的纵向坐标,xr是接收聚焦点的横向坐标,zr是接收聚焦点的纵向坐标,t是时间,是探头发射时第i个阵元对坐标(xt,zt)的指向性函数,代表探头发射时第j个阵元对坐标(xr,zr)的指向性函数,ξi是探头发射时第i个阵元的敏感系数,ηj是探头发射时第j个阵元的敏感系数,Qi代表发射时第i个阵元的变迹系数,Bj代表发射时第j个阵元的变迹系数,P代表发射信号,Yi(x,z)代表第i个阵元距离坐标(x,z)的距离,c是声速,μ是声传播时的衰减因子,τi(x,z)代表第i个阵元在坐标(x,z)上的延迟。
并将上述采集的参数的值代入声场仿真数学模型中。
根据得到的H(t,xt,zt,xr,zr)后,取
A(t,xt,zt,xr,zr)=abs(IQ(H(t,xt,zt,xr,zr))),
φ(t,xt,zt,xr,zr)=angle(IQ(H(t,xt,zt,xr,zr)),
其中,A(t,xt,zt,xr,zr)是在时刻t得到的信号的幅度,φ(t,xt,zt,xr,zr)是在时刻t得到的信号的相位,IQ代表正交解调,abs代表复数取模,angle代表复数取角度。
选取参考点为(xt,zr),则A(t,xt,zr,xr,zr)作为幅度的参考点,φ(t,xt,zr,xr,zr)作为相位的参考点。
则
Δφ(t,xt,zt,xr,zr)=φ(t,xt,zr,xr,zr)-φ(t,xt,zt,xr,zr),
所以,假设在发射聚焦(xt,zt)、接收聚焦(xr,zr)处得到的信号是:
S(xt,zt,xr,zr)=Aejφ
修正后的信号为:
Scorr(xt,zt,xr,zr)=Aejφ=A·ΔAejφ+Δφ
修正后得到第一接收线Rx3。
假设对第一接收线Rx3上的第一接收点(xi,z)进行修正,第一设定数量为4,第二设定数量也为4,确定在该第一接收点同一轴向距离上,且在该第一接收点相邻的4个接收点作为第二接收点,即该第一接收点的4个第二接收点为(xi-2,z)、(xi-1,z)、(xi+1,z)、(xi+2,z)。该第一接收点(xi,z)加上上述4个第二接收点(xi-2,z)、(xi-1,z)、(xi+1,z)、(xi+2,z),一共5个点,利用其中每4个接收点之间的相干系数进行修正。
令
j1=i-2,i-1,......i+2
j2=i-2,i-1,......i+2
j3=i-2,i-1,......i+2
j4=i-2,i-1,......i+2
j1≠j2≠j3≠j4
这里j1、j2、j3、j4是在[xi-2,xi+2]内的任意一个整数,且互不相等,即当前i为3,则j1可能为[1,5]之间任意一个数,其值可能为1、2、3、4、5中的一个,若j1为1,则j2为剩下的四个数中的任意一个数,其值可能为2、3、4、5中的一个,同时以此类推j3、j4的值。按照上述方法确定j1、j2、j3、j4,则可以确定5种不同的j1、j2、j3、j4值的组合,即可以确定5种不同的4个接收点的组合。
具体的,根据当前j1、j2、j3、j4的值,确定对应的4个接收点,并可以通过相干系数的算法公式确定该4个接收点中每一个点所在的接收线的总能量与该4个接收点所在的接收线的总能量和之间的相干性,即确定和
例如,采用CF算法计算相干系数:
其中,代表第m帧的接收线的回波信号强度,代表第m帧的图像中该点的具体位置,m代表是第m帧图像,M为包含有接收线的图像的帧数。
根据上式得到上述4个接收点对应的相干系数后进行相乘,确定当前j1、j2、j3、j4对应的即
由于有5种不同的四个接收点之间的相干系数组合,则可以确定五个对应不同的4个接收点之间的相干系数组合的
其中,当越小,这4个接收点组合的信息的可信度越低,其对的影响力也越低,甚至完全没有影响力;越大,这4个接收点组合的信息的可信度越高,其对的影响力也越高,所以需要对得到的五个进行加权,即
是对应当前的进行加权得到的。
需要说明的是,也可以是另外的加权形式,并不限于上面的形式。
最后,根据上述5种不同的包含有4个接收点的组合,确定每个组对应的相干系数权重值与该组包含的4个接收点对所述第一接收点的相干能量和的乘积的和,确定该第一接收线修正后的幅值,即
其中,
而就是第一接收点(xi,z)修正后的幅值。
按照上述的方式,确定每一条第一接收线本身的能量,对相应的第一接收线进行幅值上的处理之后,合并每条第一接收线本身的能量,生成对应的图像数据,并发送给所述图像显示单元生成对应的图像。
图像显示单元接收到多线波束合成单元发送的图像数据后,生成对应的图像。
实施例5:
图3为本发明实施例提供的一种超声成像的过程示意图,该过程包括以下步骤:
S301:针对每条第一接收线,确定该第一接收线对应的第一接收点同一轴向距离上,与所述第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;
S302:确定每个组对应的相干系数权重值,根据所述每个组对应的相干系数权重值与该组包含的所述接收点对所述第一接收点的相干能量和的乘积的和,对所述第一接收线的幅值进行修正;
S303:合并修正后的每条第一接收线,生成对应的图像数据,并发送给所述图像显示单元生成对应的图像。
其中,在确定第二接收点之前,根据预先确定的相位差和/或幅度差,对在发射聚焦、接收聚焦得到的信号的相位和/或幅度进行修正。
具体地,在同时对相位和幅度进行修正时,具体用于根据Scorr(xt,zt,xr,zr)=Aejφ=A·ΔAejφ+Δφ,对在发射聚焦、接收聚焦得到的信号的相位和幅度进行修正,其中,
Δφ(t,xt,zt,xr,zr)=φ(t,xt,zr,xr,zr)-φ(t,xt,zt,xr,zr),
其中A为所述发射聚焦、接收聚焦得到的信号的幅值,A(t,xt,zr,xr,zr)为幅度的参考点,φ(t,xt,zr,xr,zr)为相位的参考点,
且
A(t,xt,zt,xr,zr)=abs(IQ(H(t,xt,zt,xr,zr)),
φ(t,xt,zt,xr,zr)=angle(IQ(H(t,xt,zt,xr,zr))),
其中,A(t,xt,zt,xr,zr)是在时刻t得到的信号的幅度,φ(t,xt,zt,xr,zr)是在时刻t得到的信号的相位,IQ代表正交解调,abs代表复数取模,angle代表复数取角度,
其中,xt是发射聚焦点的横向坐标,zt是发射聚焦点的纵向坐标,xr是接收聚焦点的横向坐标,zr是接收聚焦点的纵向坐标,t是时间,是探头发射时第i个阵元对坐标(xt,zt)的指向性函数,代表探头发射时第j个阵元对坐标(xr,zr)的指向性函数,ξi是探头发射时第i个阵元的敏感系数,ηj是探头发射时第j个阵元的敏感系数,Qi代表发射时第i个阵元的变迹系数,Bj代表发射时第j个阵元的变迹系数,P代表发射信号,Yi(x,z)代表第i个阵元距离坐标(x,z)的距离,c是声速,μ是声传播时的衰减因子,τi(x,z)代表第i个阵元在坐标(x,y)上的延迟。
其中,所述确定每个组对应的相干系数权重值,包括:
根据
确定每个组对应的相干系数权重值,其中,a为并行波束合成倍数的二分之一,为每个组对应的可信度,为接收线的相干系数。
进一步地,所述确定每个组对应的相干系数权重值之后,还包括:
根据以下公式,确定第一接收线本身的能量,
其中,为接收线Rxjn所包含各条线数的总能量,bi,jn是计算接收线Rxjn本身的能量时,对接收线Rxjn所包含各条线数的总能量参与合成时的比例系数。
由于在本发明实施例中,根据每一条第一接收线对应的第一接收点同一轴向上,与第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;根据每个组对应的相干系数权重值与该组包含的接收点对第一接收点的相干能量和的乘积的和,对第一接收线的幅值进行修正,提高了并行波束合成的准确性,解决了现有技术中接收线数据在不可信的情况下对修正结果造成误差的问题。
实施例6:
如图4为本发明实施例还提供一种超声成像设备结构示意图,在上述各实施例的基础上,本发明实施例还提供了一种超声成像设备,所述超声成像设备包括:处理器41和存储器42;
处理器41用于执行存储器42中存储的计算机程序时实现上述跟踪聚焦的方法的步骤。
可选的,处理器41可以是CPU(中央处埋器)、ASIC(Application SpecificIntegrated Circuit,专用集成电路)、FPGA(Field-Programmable Gate Array,现场可编程门阵列)或CPLD(Complex Programmable Logic Device,复杂可编程逻辑器件)。
处理器41,用于按照存储器42中存储的计算机程序时执行如下步骤:
针对每条第一接收线,确定该第一接收线对应的第一接收点同一轴向距离上,与所述第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;
确定每个组对应的相干系数权重值,根据所述每个组对应的相干系数权重值与该组包含的所述接收点对所述第一接收点的相干能量和的乘积的和,对所述第一接收线的幅值进行修正;
合并修正后的每条第一接收线,生成对应的图像数据,并发送给所述图像显示单元生成对应的图像。
其中,在确定第二接收点之前,根据预先确定的相位差和/或幅度差,对在发射聚焦、接收聚焦得到的信号的相位和/或幅度进行修正。
具体地,在同时对相位和幅度进行修正时,具体用于根据Scorr(xt,zt,xr,zr)=Aejφ=A·ΔAejφ+Δφ,对在发射聚焦、接收聚焦得到的信号的相位和幅度进行修正,其中,
Δφ(t,xt,zt,xr,zr)=φ(t,xt,zr,xr,zr)-φ(t,xt,zt,xr,zr),
其中A为所述发射聚焦、接收聚焦得到的信号的幅值,A(t,xt,zr,xr,zr)为幅度的参考点,φ(t,xt,zr,xr,zr)为相位的参考点,
且
A(t,xt,zt,xr,zr)=abs(IQ(H(t,xt,zt,xr,zr)),
φ(t,xt,zt,xr,zr)=angle(IQ(H(t,xt,zt,xr,zr))),
其中,A(t,xt,zt,xr,zr)是在时刻t得到的信号的幅度,φ(t,xt,zt,xr,zr)是在时刻t得到的信号的相位,IQ代表正交解调,abs代表复数取模,angle代表复数取角度,
其中,xt是发射聚焦点的横向坐标,zt是发射聚焦点的纵向坐标,xr是接收聚焦点的横向坐标,zr是接收聚焦点的纵向坐标,t是时间,是探头发射时第i个阵元对坐标(xt,zt)的指向性函数,代表探头发射时第j个阵元对坐标(xr,zr)的指向性函数,ξi是探头发射时第i个阵元的敏感系数,ηj是探头发射时第j个阵元的敏感系数,Qi代表发射时第i个阵元的变迹系数,Bj代表发射时第j个阵元的变迹系数,P代表发射信号,Yi(x,z)代表第i个阵元距离坐标(x,z)的距离,c是声速,μ是声传播时的衰减因子,τi(x,z)代表第i个阵元在坐标(x,z)上的延迟。
其中,所述确定每个组对应的相干系数权重值,包括:
根据
确定每个组对应的相干系数权重值,其中,a为并行波束合成倍数的二分之一,为每个组对应的可信度,为接收线的相干系数。
进一步地,所述确定每个组对应的相干系数权重值之后,还包括:
根据以下公式,确定第一接收线本身的能量,
其中,为接收线Rxjn所包含各条线数的总能量,bi,jn是计算接收线Rxjn本身的能量时,对接收线Rxjn所包含各条线数的总能量参与合成时的比例系数。
由于在本发明实施例中,根据每一条第一接收线对应的第一接收点同一轴向上,与第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;根据每个组对应的相干系数权重值与该组包含的接收点对第一接收点的相干能量和的乘积的和,对第一接收线的幅值进行修正,提高了并行波束合成的准确性,解决了现有技术中接收线数据在不可信的情况下对修正结果造成误差的问题。
实施例7:
在上述各实施例的基础上,本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质内存储有可由电子设备执行的计算机程序,当所述程序在所述电子设备上运行时,使得所述电子设备执行时实现如下步骤:
针对每条第一接收线,确定该第一接收线对应的第一接收点同一轴向距离上,与所述第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;
确定每个组对应的相干系数权重值,根据所述每个组对应的相干系数权重值与该组包含的所述接收点对所述第一接收点的相干能量和的乘积的和,对所述第一接收线的幅值进行修正;
合并修正后的每条第一接收线,生成对应的图像数据,并发送给所述图像显示单元生成对应的图像。
其中,在确定第二接收点之前,根据预先确定的相位差和/或幅度差,对在发射聚焦、接收聚焦得到的信号的相位和/或幅度进行修正。
具体地,在同时对相位和幅度进行修正时,具体用于根据Scorr(xt,zt,xr,zr)=Aejφ=A·ΔAejφ+Δφ,对在发射聚焦、接收聚焦得到的信号的相位和幅度进行修正,其中,
Δφ(t,xt,zt,xr,zr)=φ(t,xt,zr,xr,zr)-φ(t,xt,zt,xr,zr),
其中A为所述发射聚焦、接收聚焦得到的信号的幅值,A(t,xt,zr,xr,zr)为幅度的参考点,φ(t,xt,zr,xr,zr)为相位的参考点,
且
A(t,xt,zt,xr,zr)=abs(IQ(H(t,xt,zt,xr,zr)),
φ(t,xt,zt,xr,zr)=angle(IQ(H(t,xt,zt,xr,zr))),
其中,A(t,xt,zt,xr,zr)是在时刻t得到的信号的幅度,φ(t,xt,zt,xr,zr)是在时刻t得到的信号的相位,IQ代表正交解调,abs代表复数取模,angle代表复数取角度,
其中,xt是发射聚焦点的横向坐标,zt是发射聚焦点的纵向坐标,xr是接收聚焦点的横向坐标,zr是接收聚焦点的纵向坐标,t是时间,是探头发射时第i个阵元对坐标(xt,zt)的指向性函数,代表探头发射时第j个阵元对坐标(xr,zr)的指向性函数,ξi是探头发射时第i个阵元的敏感系数,ηj是探头发射时第j个阵元的敏感系数,Qi代表发射时第i个阵元的变迹系数,Bj代表发射时第j个阵元的变迹系数,P代表发射信号,Yi(x,z)代表第i个阵元距离坐标(x,z)的距离,c是声速,μ是声传播时的衰减因子,τi(x,z)代表第i个阵元在坐标(x,z)上的延迟。
其中,所述确定每个组对应的相干系数权重值,包括:
根据
确定每个组对应的相干系数权重值,其中,a为并行波束合成倍数的二分之一,
为每个组对应的可信度,为接收线的相干系数。
进一步地,所述确定每个组对应的相干系数权重值之后,还包括:
根据以下公式,确定第一接收线本身的能量,
其中,为接收线Rxjn所包含各条线数的总能量,bi,jn是计算接收线Rxjn本身的能量时,对接收线Rxjn所包含各条线数的总能量参与合成时的比例系数。
上述计算机可读存储介质可以是电子设备中的处理器能够存取的任何可用介质或数据存储设备,包括但不限于磁性存储器如软盘、硬盘、磁带、磁光盘(MO)等、光学存储器如CD、DVD、BD、HVD等、以及半导体存储器如ROM、EPROM、EEPROM、非易失性存储器(NANDFLASH)、固态硬盘(SSD)等。
由于在本发明实施例中,根据每一条第一接收线对应的第一接收点同一轴向上,与第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;根据每个组对应的相干系数权重值与该组包含的接收点对第一接收点的相干能量和的乘积的和,对第一接收线的幅值进行修正,提高了并行波束合成的准确性,解决了现有技术中接收线数据在不可信的情况下对修正结果造成误差的问题。
对于系统/装置实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (12)
1.一种超声成像系统,所述系统包括探头、发射波束单元、多线接收波束单元以及图像显示单元,其特征在于,还包括:多线波束合成单元;
所述多线波束合成单元与所述图像显示单元连接,用于针对每条第一接收线,确定该第一接收线对应的第一接收点同一轴向距离上,与所述第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;确定每个组对应的相干系数权重值,根据所述每个组对应的相干系数权重值与该组包含的所述接收点对所述第一接收点的相干能量和的乘积的和,对所述第一接收线的幅值进行修正;合并修正后的每条第一接收线,生成对应的图像数据,并发送给所述图像显示单元生成对应的图像。
2.根据权利要求1所述的系统,其特征在于,所述多线波束合成单元,还用于在确定第二接收点之前,根据预先确定的相位差和/或幅度差,对在发射聚焦、接收聚焦得到的信号的相位和/或幅度进行修正。
3.根据权利要求2所述的系统,其特征在于,所述多线波束合成单元在同时对相位和幅度进行修正时,具体用于根据Scorr(xt,zt,xr,zr)=Aejφ=A·ΔAejφ+Δφ,对在发射聚焦、接收聚焦得到的信号的相位和幅度进行修正,其中,
Δφ(t,xt,zt,xr,zr)=φ(t,xt,zr,xr,zr)-φ(t,xt,zt,xr,zr),
其中A为所述发射聚焦、接收聚焦得到的信号的幅值,A(t,xt,zr,xr,zr)为幅度的参考点,φ(t,xt,zr,xr,zr)为相位的参考点,
且
A(t,xt,zt,xr,zr)=abs(IQ(H(t,xt,zt,xr,zr)),
φ(t,xt,zt,xr,zr)=angle(IQ(H(t,xt,zt,xr,zr))),
其中,A(t,xt,zt,xr,zr)是在时刻t得到的信号的幅度,φ(t,xt,zt,xr,zr)是在时刻t得到的信号的相位,IQ代表正交解调,abs代表复数取模,angle代表复数取角度,
其中,xt是发射聚焦点的横向坐标,zt是发射聚焦点的纵向坐标,xr是接收聚焦点的横向坐标,zr是接收聚焦点的纵向坐标,t是时间,是探头发射时第i个阵元对坐标(xt,zt)的指向性函数,代表探头发射时第j个阵元对坐标(xr,zr)的指向性函数,ξi是探头发射时第i个阵元的敏感系数,ηj是探头发射时第j个阵元的敏感系数,Qi代表发射时第i个阵元的变迹系数,Bj代表发射时第j个阵元的变迹系数,P代表发射信号,Yi(x,z)代表第i个阵元距离坐标(x,z)的距离,c是声速,μ是声传播时的衰减因子,τi(x,z)代表第i个阵元在坐标(x,z)上的延迟。
4.根据权利要求1所述的系统,其特征在于,所述多线波束合成单元,具体用于根据
确定每个组对应的相干系数权重值,其中,a为并行波束合成倍数的二分之一,为每个组对应的可信度,为接收线的相干系数。
5.根据权利要求4所述的系统,其特征在于,所述多线波束合成单元,具体用于根据以下公式,确定第一接收线本身的能量,
其中,为接收线Rxjn所包含各条线数的总能量,bi,jn是计算接收线Rxjn本身的能量时,对接收线Rxjn所包含各条线数的总能量参与合成时的比例系数。
6.一种超声成像方法,其特征在于,所述方法包括:
针对每条第一接收线,确定该第一接收线对应的第一接收点同一轴向距离上,与所述第一接收点相邻的第一设定数量的第二接收点,将接收点划分为多个组,其中每个组中包含第二设定数量的所述接收点,每个组中包含的所述接收点完全不同;
确定每个组对应的相干系数权重值,根据所述每个组对应的相干系数权重值与该组包含的所述接收点对所述第一接收点的相干能量和的乘积的和,对所述第一接收线的幅值进行修正;
合并修正后的每条第一接收线,生成对应的图像数据,并发送给所述图像显示单元生成对应的图像。
7.根据权利要求6所述的方法,其特征在于,在确定第二接收点之前,根据预先确定的相位差和/或幅度差,对在发射聚焦、接收聚焦得到的信号的相位和/或幅度进行修正。
8.根据权利要求7所述的方法,其特征在于,在同时对相位和幅度进行修正时,根据Scorr(xt,zt,xr,zr)=Aejφ=A·ΔAejφ+Δφ,对在发射聚焦、接收聚焦得到的信号的相位和幅度进行修正,其中,
Δφ(t,xt,zt,xr,zr)=φ(t,xt,zr,xr,zr)-φ(t,xt,zt,xr,zr),
其中A为所述发射聚焦、接收聚焦得到的信号的幅值,A(t,xt,zr,xr,zr)为幅度的参考点,φ(t,xt,zr,xr,zr)为相位的参考点,
且
A(t,xt,zt,xr,zr)=abs(IQ(H(t,xt,zt,xr,zr)),
φ(t,xt,zt,xr,zr)=angle(IQ(H(t,xt,zt,xr,zr))),
其中,A(t,xt,zt,xr,zr)是在时刻t得到的信号的幅度,φ(t,xt,zt,xr,zr)是在时刻t得到的信号的相位,IQ代表正交解调,abs代表复数取模,angle代表复数取角度,
其中,xt是发射聚焦点的横向坐标,zt是发射聚焦点的纵向坐标,xr是接收聚焦点的横向坐标,zr是接收聚焦点的纵向坐标,t是时间,是探头发射时第i个阵元对坐标(xt,zt)的指向性函数,代表探头发射时第j个阵元对坐标(xr,zr)的指向性函数,ξi是探头发射时第i个阵元的敏感系数,ηj是探头发射时第j个阵元的敏感系数,Qi代表发射时第i个阵元的变迹系数,Bj代表发射时第j个阵元的变迹系数,P代表发射信号,Yi(x,z)代表第i个阵元距离坐标(x,z)的距离,c是声速,μ是声传播时的衰减因子,τi(x,z)代表第i个阵元在坐标(x,z)上的延迟。
9.根据权利要求6所述的方法,其特征在于,所述确定每个组对应的相干系数权重值,包括:
根据
确定每个组对应的相干系数权重值,其中,a为并行波束合成倍数的二分之一,为每个组对应的可信度,为接收线的相干系数。
10.根据权利要求9所述的方法,其特征在于,所述确定每个组对应的相干系数权重值之后,还包括:
根据以下公式,确定第一接收线本身的能量,
其中,为接收线Rxjn所包含各条线数的总能量,bi,jn是计算接收线Rxjn本身的能量时,对接收线Rxjn所包含各条线数的总能量参与合成时的比例系数。
11.一种超声成像设备,其特征在于,所述超声成像设备包括处理器,所述处理器用于执行存储器中存储的计算机程序时实现如权利要求6-10中任一所述方法的步骤。
12.一种计算机可读存储介质,其特征在于,其存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求6-10中任一所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910893370.5A CN110507355B (zh) | 2019-09-20 | 2019-09-20 | 一种超声成像系统、方法、设备及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910893370.5A CN110507355B (zh) | 2019-09-20 | 2019-09-20 | 一种超声成像系统、方法、设备及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110507355A true CN110507355A (zh) | 2019-11-29 |
CN110507355B CN110507355B (zh) | 2022-04-15 |
Family
ID=68631646
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910893370.5A Active CN110507355B (zh) | 2019-09-20 | 2019-09-20 | 一种超声成像系统、方法、设备及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110507355B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021114071A1 (zh) * | 2019-12-10 | 2021-06-17 | 深圳大学 | 一种经颅超声成像方法、装置及计算机可读存储介质 |
CN113030983A (zh) * | 2021-03-17 | 2021-06-25 | 中国科学院声学研究所 | 一种基于测深侧扫声纳的近场逐点聚焦doa方法 |
CN115153622A (zh) * | 2022-06-08 | 2022-10-11 | 东北大学 | 一种基于虚源的基带延时乘累加超声波束形成方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6544184B1 (en) * | 2001-08-28 | 2003-04-08 | Acuson Corporation | Imaging with reduced artifacts for medical diagnostic ultrasound |
US20060173313A1 (en) * | 2005-01-27 | 2006-08-03 | Siemens Medical Solutions Usa, Inc. | Coherence factor adaptive ultrasound imaging |
CN102435992A (zh) * | 2011-09-26 | 2012-05-02 | 重庆博恩富克医疗设备有限公司 | 基于广义相干系数的合成聚焦成像方法 |
JP2013079949A (ja) * | 2011-09-30 | 2013-05-02 | Ge Inspection Technologies Ltd | 撮像システムおよび撮像方法 |
CN107260217A (zh) * | 2017-07-17 | 2017-10-20 | 西安交通大学 | 用于脑部聚焦超声空化实时监控的三维无源成像方法及系统 |
CN109363714A (zh) * | 2018-09-17 | 2019-02-22 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像设备及其超声成像方法 |
-
2019
- 2019-09-20 CN CN201910893370.5A patent/CN110507355B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6544184B1 (en) * | 2001-08-28 | 2003-04-08 | Acuson Corporation | Imaging with reduced artifacts for medical diagnostic ultrasound |
US20060173313A1 (en) * | 2005-01-27 | 2006-08-03 | Siemens Medical Solutions Usa, Inc. | Coherence factor adaptive ultrasound imaging |
CN102435992A (zh) * | 2011-09-26 | 2012-05-02 | 重庆博恩富克医疗设备有限公司 | 基于广义相干系数的合成聚焦成像方法 |
JP2013079949A (ja) * | 2011-09-30 | 2013-05-02 | Ge Inspection Technologies Ltd | 撮像システムおよび撮像方法 |
CN107260217A (zh) * | 2017-07-17 | 2017-10-20 | 西安交通大学 | 用于脑部聚焦超声空化实时监控的三维无源成像方法及系统 |
CN109363714A (zh) * | 2018-09-17 | 2019-02-22 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像设备及其超声成像方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021114071A1 (zh) * | 2019-12-10 | 2021-06-17 | 深圳大学 | 一种经颅超声成像方法、装置及计算机可读存储介质 |
CN113030983A (zh) * | 2021-03-17 | 2021-06-25 | 中国科学院声学研究所 | 一种基于测深侧扫声纳的近场逐点聚焦doa方法 |
CN115153622A (zh) * | 2022-06-08 | 2022-10-11 | 东北大学 | 一种基于虚源的基带延时乘累加超声波束形成方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110507355B (zh) | 2022-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110507355B (zh) | 一种超声成像系统、方法、设备及介质 | |
US10976422B2 (en) | Ultrasound imaging methods and systems | |
CN107789008B (zh) | 一种基于通道数据的自适应超声波束合成方法和系统 | |
CN106500820B (zh) | 一种用于二维波达方向估计的声速测量方法及装置 | |
CN109975814A (zh) | 超声成像方法、系统和设备 | |
CN107942284A (zh) | 基于二维正交非均匀线阵的水下波达方向估计方法与装置 | |
CN104605890A (zh) | 一种剪切波峰值波形校正的方法、装置、系统及其应用 | |
Rindal et al. | Double adaptive plane-wave imaging | |
JPWO2015166869A1 (ja) | 超音波撮像装置 | |
JP5767002B2 (ja) | 超音波送受信装置、および魚量検出方法 | |
Handegard et al. | Tracking individual fish from a moving platform using a split-beam transducer | |
CN106373103A (zh) | 一种超声数据的复合方法及装置 | |
JP4761999B2 (ja) | 超音波診断装置およびその画像処理方法、その画像処理プログラム | |
JP2022552781A (ja) | 組織弾性検出方法及び装置 | |
CN110554393B (zh) | 基于深度学习的高对比度最小方差成像方法 | |
CN103841897A (zh) | 信号处理装置及方法、记录介质和程序 | |
JP2017000547A (ja) | 超音波診断装置 | |
CN112120734B (zh) | 血流方向的多普勒频谱生成方法、装置及相关设备 | |
CN110693524B (zh) | 一种超声医学成像聚焦校正方法和装置 | |
CN110613476B (zh) | 超声信号的处理方法和装置 | |
Filipik et al. | Time-of-flight based calibration of an ultrasonic computed tomography system | |
JP5950534B2 (ja) | 超音波距離画像生成装置 | |
KR101510678B1 (ko) | 고조파 영상 형성 방법과 그를 위한 초음파 의료 장치 | |
CN113925528B (zh) | 多普勒成像方法和超声设备 | |
JP6311230B2 (ja) | 目標物検出装置、目標物検出方法、プログラム及び記録媒体 |
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 |