CN103425894A - 基于误差约束的动态聚焦最优分段方法及装置 - Google Patents

基于误差约束的动态聚焦最优分段方法及装置 Download PDF

Info

Publication number
CN103425894A
CN103425894A CN2013103811920A CN201310381192A CN103425894A CN 103425894 A CN103425894 A CN 103425894A CN 2013103811920 A CN2013103811920 A CN 2013103811920A CN 201310381192 A CN201310381192 A CN 201310381192A CN 103425894 A CN103425894 A CN 103425894A
Authority
CN
China
Prior art keywords
focusing
segmentation
delay
dynamic
dynamic focusing
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
CN2013103811920A
Other languages
English (en)
Other versions
CN103425894B (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.)
Chongqing University
State Grid Eastern Inner Mongolia Power Co Ltd
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN201310381192.0A priority Critical patent/CN103425894B/zh
Publication of CN103425894A publication Critical patent/CN103425894A/zh
Application granted granted Critical
Publication of CN103425894B publication Critical patent/CN103425894B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于误差约束的动态聚焦最优分段方法及装置;该方法首先根据给定聚焦通道数和采样率生成聚焦延时表,然后结合聚焦误差δ,对聚焦延时表进行依次分段计算,在每个分段区间选取一个特定聚焦点,使得该分段聚焦区域内聚焦延时数据相对于该聚焦点延时参数误差控制在±δ/2范围内,依此方法对整个聚焦延时表进行动态分段;根据分段聚焦延时和聚焦误差δ,生成分段聚焦延时数据修正信息,聚焦延时参数发生装置在动态分段聚焦过程中,根据聚焦延时修正信息实时修正聚焦延时参数;本发明对全程探测范围内动态聚焦进行最优分段,并保证在全程动态聚焦时候,将聚焦延时参数误差限定在约束范围内,实现动态聚焦最优动态分段,最大限度降低了存储量。

Description

基于误差约束的动态聚焦最优分段方法及装置
技术领域
本发明属于超声成像技术领域,具体涉及一种基于超声成像系统的分段方法和聚焦延时参数的实时生成装置。
背景技术
动态聚焦的聚焦延时是随阵元与聚焦点间距离变换的函数。成功的动态聚焦效果能达到保持每条扫描线上的所有点都在焦点上,从而最大限度的提高超声图像的横向分辨率。动态聚焦要求控制系统能以回波相同的速度沿扫描线追踪目标,以形成一个滑动的焦点。在工程实际中,理想的动态聚焦难以实现。为了弥补这一缺陷,一般采用分段动态聚焦,即将成像空间分成多个小段,对每一小段接收采用定点聚焦。但是传统的分段聚焦是简单地将探测深度均匀分段,这样就会造成在近场区,聚焦点与分段边缘探测点间聚焦延时差大,而使得聚焦误差较大,成像质量较差;而在中远场区,聚焦点与段边缘探测点间聚焦延时差小,不能凸显分段聚焦简化计算的优点,并且对聚焦延时参数的存储量也比较大。
发明内容
有鉴于此,本发明的目的是提供一种基于误差约束的动态聚焦最优分段方法和装置,在保证成像质量的前提下,合理、动态分段,实现了动态聚焦的最优分段;同时,本发明在此方法的基础上,提供了一种最优分段聚焦延时参数的实时生成装置,简化了聚焦延时参数实时生成的难题。
本发明的目的之一是通过以下技术方案来实现:基于误差约束的动态聚焦最优分段方法,包括如下步骤:
1)根据探测深度L、采样率fs和聚焦通道数N,通过公式
Figure BDA0000373655100000011
计算得到动态聚焦延时参数表TMjn,其中
Figure BDA0000373655100000012
1≤n≤N,c为超声波的声速,j为采样数据变量,n表示通道变量,N表示聚焦通道数,rjn为第n号阵元到第j个聚焦点的距离,Mj为聚焦中心起始点到全动态聚焦时第j个聚焦点距离;
2)根据聚焦误差δ,计算出第k分段位置在动态聚焦延时参数表中的位置序号Pk和第k分段区域内聚焦点对应的聚焦延时参数Tkn,其中1≤k≤K,K表示分段总数,k表示为分段变量数;
3)根据聚焦延时参数Tkn计算得到分段动态聚焦的修正参数Rkn
4)微控制器MCU装载第一组聚焦延时参数T1n到各个通道对应的递减器中,然后将分段动态聚焦的修正参数Rkn依次写入双口RAM中;将分段动态聚焦的分段位置数据Pk依次写入地址发生器中;
5)当分段动态聚焦开始时,根据修正参数Rkn对各个聚焦通道所需的聚焦延时参数Tkn进行实时修正。
进一步,所述步骤1)具体包括以下子步骤:
11)首先根据探测深度L的起始位置设定为分段聚焦第一分段区域的起始位置,记为P1,计算逐点聚焦时最边沿第1号阵元的聚焦延迟时间TMj1式中rj1为第1号阵元到全动态聚焦时第j个聚焦点距离,rj1计算公式为:式
Figure BDA0000373655100000022
其中Mj为全动态聚焦时聚焦中心起始点到第j个聚焦点距离,x1表示第1号阵元到聚焦中心线的距离;
12)根据聚焦通道数n计算得到动态聚焦延时参数表TMjn
进一步,所述步骤2)具体包括以下子步骤:
21):初始化j=1,i=1,k=1,其中j表示第1号阵元的聚焦延迟时间索引下标,i是索引下标的偏移量,k为分段变量;
22):如果TMj1-TM(j+i)1≤δ/2,其中约束条件:
Figure BDA0000373655100000023
那么i'=i+1,i'表示下次计算时i的取值,i表示本次计算时的取值,重复步骤22);否则进入步骤23);
23):提取第k分段的聚焦延时参数:Tkn=TM(j+i-1)n,程序进入步骤24);
24):如果TMj1-TM(j+i)1≤δ,其中约束条件:
Figure BDA0000373655100000024
那么i'=i+1,i'表示下次计算时i的取值,i表示本次计算时i的取值,重复步骤24);否则进入步骤25);
25):记录第k分段的位置序号,Pk=j+i-1,k'=k+1;k'表示下次计算时的取值,k表示本次计算时k的取值;
26):如果
Figure BDA0000373655100000031
程序退出,否则j'=j+i,j'表示下次计算时j的取值,i,j表示本次计算时的取值,然后重新初始化i=1,程序回到步骤22);
上述步骤中,δ为允许聚焦延迟误差,Pk为第k次分段在动态聚焦延时参数表中的位置序号,TMjn为第j个采样数据,第n通道的聚焦延迟时间,Tkn为第k分段、第n通道的聚焦延时参数。
进一步,所述步骤3)具体包括以下子步骤:
31):初始化变量i=1,k=1;
32):若Tkn-T(k+i)n<δ,其中约束条件为,1≤i<K,则分段动态聚焦的修正参数R(k+i)n=0,i'=i+1,i'表示下次计算时i的取值,i表示本次计算时的取值,重复步骤32),否则进入步骤33);
33):若Tkn-T(k+i)n≥δ,其中约束条件为,1≤i<K,则分段动态聚焦的修正参数R(k+i)n=1,T(k+i)n=Tkn-δ,k'=k+i,k'表示下次计算时k的取值,k表示本次计算时的取值,然后重新初始化i=1,程序回到步骤32)。
上述步骤中,δ为允许聚焦延迟误差,Tkn为第k分段、第n通道的聚焦延时参数,Rkn为分段动态聚焦的修正参数。
进一步,所述步骤4)具体包括以下子步骤:
41)根据聚焦延时参数Tkn得到第一组聚焦延时参数T1n,将T1n装载到各个通道对应的递减器中,然后将步骤3)得到的分段动态聚焦的修正参数Rkn依次写入双口RAM中;
42)将步骤2)得到总共K段的分段位置序号Pk依次写入地址发生器中。
进一步,所述步骤5)具体包括以下子步骤:
51)当分段动态聚焦开始时,在聚焦时钟的作用下,地址发生器内部的计数器开始计数工作,当计数器等于Pk时,地址发生器输出为k,从双口RAM中读出第k段聚焦区间的分段动态聚焦修正参数Rkn
52)地址发生器产生一个时钟脉冲,在时钟脉冲的上升沿,聚焦延时参数Tkn按下式进行修正
Figure BDA0000373655100000032
从而实现了分段聚焦延时参数Tkn的实时生成。
本发明的目的之二是通过以下技术方案来实现的,基于误差约束的动态聚焦最优分段装置,包括微控制器、双口RAM、地址发生器和递减器,所述微控制器的数据端分别与双口RAM的第一数据端和地址发生器的数据端连接,微控制器的地址端分别与双口RAM的第一地址端和地址发生器的地址端连接,微控制器的第一片选端与双口RAM的第一片选端,微控制器的第二片选端与地址发生器的第一片选端连接,微控制器的写控制端分别与双口RAM的写控制端和地址发生器的写控制端连接,微控制器的时钟端分别与双口RAM的第一时钟端和地址发生器的第一时钟端连接,微控制器的IO口与地址发生器的使能端连接;所述双口RAM输出的数据位分别输入到每个递减器中,地址发生器产生的时钟脉冲C_CLK分别输入到每个递减器中,地址发生器的读控制端分别与双口RAM的读控制端和第二片选端连接,时钟信号分别输入到双口RAM的第二时钟端和地址发生器的第二时钟端,地址发生器输出的CNT端与双口RAM的第二地址端连接。工作时,微控制器装载第一组聚焦延时参数到各个通道对应的递减器中,然后将分段动态聚焦的修正参数依次写入双口RAM中;将分段动态聚焦的分段位置数据Pk依次写入地址发生器中。
进一步,动态聚焦最优分段方法包括以下步骤:
1)微控制器根据给定的探测深度、聚焦通道数和采样率生成聚焦延时表;
2)结合给定的聚焦误差δ对整个聚焦延时表进行动态分段,并计算出分段位置数据和分段聚焦延时数据;
3)根据步骤2)生成分段聚焦延时数据的修正信息;
4)微控制器装载第一组聚焦延时参数到各个通道对应的递减器中,然后将分段动态聚焦的修正参数依次写入双口RAM中;将分段动态聚焦的分段位置数据依次写入地址发生器中;
5)根据聚焦延时修正信息实时修正聚焦延时参数。
有益技术效果:
该分段方法首先根据给定的聚焦通道数和采样率生成聚焦延时表,然后结合给定的聚焦误差δ,对聚焦延时表进行依次分段计算,在每一个分段区间选取一个特定的聚焦点,使得该分段聚焦区域内的聚焦延时数据相对于该聚焦点的延时参数误差控制在±δ/2范围内,从而完成对整个聚焦延时表进行动态分段;根据分段位置数据和分段聚焦延时数据,生成分段聚焦延时数据的修正信息,聚焦延时参数发生装置在动态分段聚焦过程中,根据聚焦延时修正信息实时修正聚焦延时参数。本发明对全程探测范围内的动态聚焦进行最优分段,并保证在全程动态聚焦的时候,将聚焦误差限定在约束的误差范围内,从而实现了动态聚焦的最优动态分段,最大限度地降低了存储量。在保证成像质量的前提下,合理、动态分段,实现了动态聚焦的最优分段。同时,本发明在此方法的基础上,提供了一种最优分段聚焦延时参数的实时生成装置,简化了聚焦延时参数实时生成的难题。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步的详细描述,其中:
图1为本发明基于误差约束的动态聚焦最优分段方法流程图;
图2示出了本发明基于误差约束的动态聚焦最优分段方法的原理示意图;
图3示出了本发明实施例中,采用32阵元聚焦分段聚焦延时参数图;
图4示出了本发明32通道的基于误差约束的动态聚焦最优分段方法的分段聚焦延时参数实时生成装置的原理框图;
具体实施方式
以下将结合附图,对本发明的优选实施例进行详细的描述;应当理解,优选实施例仅为了说明本发明,而不是为了限制本发明的保护范围。
参见图2,其中,o为坐标原点,d为阵元间距,A为探头阵元,B为分段聚焦分的某一段,C为某一段内聚焦点,阵元个数为32个。图中所示的分段聚焦希望得到的分段效果是使分段达到一种最优,在聚焦延时允许的范围内最大深度间距的分段,得到最少的分段段数,从而大大简化计算。在本实施例中,聚焦目标处于探头阵列的中心聚焦线位置。由于探头关于其垂直方向中心轴对称,所以本实例取一半阵元作为研究对象。同时,由于对确定焦点聚焦时,离中心聚焦线起始位置越远,阵元的聚焦延时越大,所以作者选取边沿第1号阵元聚焦延时作为分段判断依据。
本实施例基于误差约束的动态聚焦最优分段方法,针对动态聚焦中的分段聚焦进行改进,包括如下步骤:
1)根据探测深度L、采样率fs和聚焦通道数N,通过公式
Figure BDA0000373655100000061
计算得到动态聚焦延时参数表TMjn,其中1≤n≤N,c为超声波的声速,j为采样数据变量,n表示通道变量,N表示聚焦通道数,rjn为第n号阵元到第j个聚焦点的距离,Mj为聚焦中心起始点到全动态聚焦时第j个聚焦点距离;
11)首先将探测深度L的起始位置设定为分段聚焦第一分段区域的起始位置,记为P1,计算逐点聚焦时最边沿第1号阵元的聚焦延迟时间TMj1式中rj1为第1号阵元到全动态聚焦时第j个聚焦点距离,rj1计算公式为:式
Figure BDA0000373655100000064
其中Mj为聚焦中心起始点到全动态聚焦时第j个聚焦点距离,x1表示第1号阵元到聚焦中心线的距离;
12)根据聚焦通道数n(1≤n≤N)计算得到动态聚焦延时参数表TMjn
2)根据聚焦误差δ,计算出第k分段位置在动态聚焦延时参数表中的位置序号Pk和第k分段区域内聚焦点对应的聚焦延时参数Tkn,其中1≤k≤K,K表示分段总数,k表示为分段变量数;
21):初始化j=1,i=1,k=1,其中j表示第1号阵元的聚焦延迟时间索引下标,i是索引下标的偏移量,k为分段变量;
22):如果TMj1-TM(j+i)1≤δ/2,其中约束条件:
Figure BDA0000373655100000065
那么i'=i+1,i'表示下次计算时i的取值,i表示本次计算时i的取值,重复步骤22);否则进入步骤23);
23):提取第k分段的聚焦延时参数:Tkn=TM(j+i-1)n(1≤n≤N),程序进入步骤24);
24):如果TMj1-TM(j+i)1≤δ,其中约束条件:
Figure BDA0000373655100000066
那么i'=i+1,i'表示下次计算时i的取值,i表示本次计算时i的取值,重复步骤24);否则进入步骤25);
25):记录第k分段的位置序号,Pk=j+i-1,k'=k+1;k'表示下次计算时k的取值,k表示本次计算时k的取值;
26):如果
Figure BDA0000373655100000067
程序退出,否则j'=j+i,j'表示下次计算时j的取值,i,j表示本次计算时i,j的取值,然后重新初始化i=1,程序回到步骤22);
上述步骤中,δ为允许聚焦延迟误差,Pk为第k次分段在动态聚焦延时参数表中的位置序号,TMjn为第j个采样数据,第n通道的聚焦延迟时间,Tkn为第k分段、第n通道的聚焦延时参数。
3)根据聚焦延时参数Tkn计算得到分段动态聚焦的修正参数Rkn
31):初始化变量i=1,k=1;
32):若Tkn-T(k+i)n<δ,其中约束条件为,1≤i<K,则分段动态聚焦的修正参数R(k+i)n=0,i'=i+1,i'表示下次计算时i的取值,i表示本次计算时i的取值,重复步骤32),否则进入步骤33);
33):若Tkn-T(k+i)n≥δ,其中约束条件为,1≤i<K,则分段动态聚焦的修正参数R(k+i)n=1,T(k+i)n=Tkn-δ,k'=k+i,k'表示下次计算时k的取值,k表示本次计算时k的取值,然后重新初始化i=1,程序回到步骤32)。
上述步骤中,δ为允许聚焦延迟误差,Tkn为第k分段、第n通道的聚焦延时参数,Rkn为分段动态聚焦的修正参数。
4)微控制器MCU装载第一组聚焦延时参数T1n到各个通道对应的递减器中,然后将分段动态聚焦的修正参数Rkn依次写入双口RAM中;将分段动态聚焦的分段位置数据Pk依次写入地址发生器中;
41)根据聚焦延时参数Tkn得到第一组聚焦延时参数T1n(1≤n≤N),将T1n 载到各个通道对应的递减器中,然后将步骤3)得到的分段动态聚焦的修正参数Rkn依次写入双口RAM中;
42)将步骤2)得到总共K段的分段位置序号Pk依次写入地址发生器中。
5)当分段动态聚焦开始时,根据修正参数Rkn对各个聚焦通道所需的聚焦延时参数Tkn进行实时修正。
51)当分段动态聚焦开始时,在聚焦时钟FOCUS_CLK(输入到地址发生器中)的作用下,地址发生器内部的计数器COUNTER开始计数工作,当计数器COUNTER=Pk时,地址发生器输出CNT[7..0]=k,从双口RAM中读出第k段聚焦区间的分段动态聚焦修正参数Rkn
52)地址发生器产生一个时钟脉冲C_CLK,在时钟脉冲C_CLK的上升沿,聚焦延时参数Tkn按下式进行修正
Figure BDA0000373655100000071
从而实现了分段聚焦延时参数Tkn的实时生成。
本发明还公布了基于误差约束的动态聚焦最优分段装置,包括微控制器、双口RAM、地址发生器和递减器,所述微控制器的数据端分别与双口RAM的第一数据端和地址发生器的数据端连接,微控制器的地址端分别与双口RAM的第一地址端和地址发生器的地址端连接,微控制器的第一片选端与双口RAM的第一片选端,微控制器的第二片选端与地址发生器的第一片选端连接,微控制器的写控制端分别与双口RAM的写控制端和地址发生器的写控制端连接,微控制器的时钟端分别与双口RAM的第一时钟端和地址发生器的第一时钟端连接,微控制器的IO口与地址发生器的使能端连接;所述双口RAM输出的数据位分别输入到每个递减器中,地址发生器产生的时钟脉冲C_CLK分别输入到每个递减器中,地址发生器的读控制端分别与双口RAM的读控制端和第二片选端连接,时钟信号分别输入到双口RAM的第二时钟端和地址发生器的第二时钟端,地址发生器输出的CNT端与双口RAM的第二地址端连接。
图2是基于误差约束的动态聚焦最优分段方法及装置选取32阵元超声探头分段的效果,横轴为分段段数,纵轴为对应段数时的聚焦延时。其中探头参数为,阵元间距d=0.48mm,超声声速c=1540m/s,脉冲中心频率fc=3.5MHz,全聚焦采样频率为fs=50MHz,聚焦范围为2mm~240mm,聚焦延时允许误差为40ns。
图4是基于误差约束的动态聚焦最优分段方法及装置选取32阵元超声探头分段聚焦延时参数实时生成装置的原理框图,给出了实时生成分段动态聚焦的修正参数Rkn,和用一系列递减信号逼近图3(横坐标表示分段次数k,纵坐标表示分段延时参数Tkn)中曲线的原理图(从上往下,曲线依次表示为第1号阵元、第2号阵元……超声探头生成的曲线)。
该分段方法首先根据给定的探测深度L、聚焦通道数和采样率生成聚焦延时表;然后结合给定的聚焦误差δ,对聚焦延时表进行依次分段计算,在每一个分段区间选取一个特定的聚焦点,使得该分段聚焦区域内的聚焦延时数据相对于该聚焦点的延时参数误差控制在±δ/2范围内,从而完成对整个聚焦延时表进行动态分段;根据分段位置数据和分段聚焦延时数据,生成分段聚焦延时数据的修正信息;工作时,微控制器装载第一组聚焦延时参数到各个通道对应的递减器中,然后将分段动态聚焦的修正参数依次写入双口RAM中,将分段动态聚焦的分段位置数据Lk依次写入地址发生器中;聚焦延时参数发生装置在动态分段聚焦过程中,根据聚焦延时修正信息实时修正聚焦延时参数。本发明对全程探测范围内的动态聚焦进行最优分段,并保证在全程动态聚焦的时候,将聚焦误差限定在约束的范围内,从而实现了动态聚焦的最优动态分段,最大限度地降低了存储量。在保证成像质量的前提下,合理、动态分段,实现了动态聚焦的最优分段;同时,本发明在此方法的基础上,提供了一种最优分段聚焦延时参数的实时生成装置,简化了聚焦延时参数实时生成的难题。
以上所述仅为本发明的优选实施例,并不用于限制本发明,显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (8)

1.基于误差约束的动态聚焦最优分段方法及装置,其特征在于,包括如下步骤:
1)根据探测深度L、采样率fs和聚焦通道数N,通过公式
Figure FDA0000373655090000011
计算得到动态聚焦延时参数表TMjn,其中
Figure FDA0000373655090000012
1≤n≤N,c为超声波的声速,j为采样数据变量,n表示通道变量,N表示聚焦通道数,rjn为第n号阵元到第j个聚焦点的距离,Mj为全动态聚焦时聚焦中心起始点到第j个聚焦点距离;
2)根据聚焦误差δ,计算出第k分段位置在动态聚焦延时参数表中的位置序号Pk和第k分段区域内聚焦点对应的聚焦延时参数Tkn,其中1≤k≤K,K表示分段总数,k表示分段变量;
3)根据聚焦延时参数Tkn计算得到分段动态聚焦的修正参数Rkn
4)微控制器MCU装载第一组聚焦延时参数T1n到各个通道对应的递减器中,然后将分段动态聚焦的修正参数Rkn依次写入双口RAM中;将分段动态聚焦的分段位置数据Pk依次写入地址发生器中;
5)当分段动态聚焦开始时,根据修正参数Rkn对各个聚焦通道所需的聚焦延时参数Tkn进行实时修正,得到各个通道所需聚焦延时参数。
2.根据权利要求1所述的基于误差约束的动态聚焦最优分段方法及装置,其特征在于:所述步骤1)具体包括以下子步骤:
11)首先将探测深度L的起始位置设定为分段聚焦第一分段区域的起始位置,记为P1,计算逐点聚焦时最边沿第1号阵元的聚焦延迟时间TMj1
Figure FDA0000373655090000013
式中rj1为第1号阵元到全动态聚焦时第j个聚焦点距离,rj1计算公式为:式
Figure FDA0000373655090000014
其中Mj为全动态聚焦时聚焦中心起始点到第j个聚焦点距离,x1表示第1号阵元到聚焦中心线的距离;
12)根据聚焦通道数n计算得到动态聚焦延时参数表TMjn
3.根据权利要求2所述的基于误差约束的动态聚焦最优分段方法及装置,其特征在于:所述步骤2)具体包括以下子步骤:
21):初始化j=1,i=1,k=1,其中j表示第1号阵元的聚焦延迟时间索引下标,i是索引下标的偏移量;
22):如果TMj1-TM(j+i)1≤δ/2,其中约束条件:
Figure FDA0000373655090000015
那么i'=i+1,i'表示下次计算时i的取值,i表示本次计算时的取值,重复步骤22);否则进入步骤23);
23):提取第k分段的聚焦延时参数:Tkn=TM(j+i-1)n,程序进入步骤24);
24):如果TMj1-TM(j+i)1≤δ,其中约束条件:
Figure FDA0000373655090000021
那么i'=i+1,i'表示下次计算时i的取值,i表示本次计算时i的取值,重复步骤24);否则进入步骤25);
25):记录第k分段的位置序号,Pk=j+i-1,k'=k+1;k'表示下次计算时k的取值,k表示本次计算时的取值;
26):如果程序退出,否则j'=j+i,j'表示下次计算时j的取值,i,j表示本次计算时的取值,然后重新初始化i=1,程序回到步骤22);
上述步骤中,δ为允许聚焦延迟误差,Pk为第k次分段在动态聚焦延时参数表中的位置序号,TMjn为第j个采样数据,第n通道的聚焦延迟时间,Tkn为第k分段、第n通道的聚焦延时参数。
4.根据权利要求1所述的基于误差约束的动态聚焦最优分段方法及装置,其特征在于:所述步骤3)具体包括以下子步骤:
31):初始化变量i=1,k=1;
32):若Tkn-T(k+i)n<δ,其中约束条件为,1≤i<K,则分段动态聚焦的修正参数R(k+i)n=0,i'=i+1,i'表示下次计算时i的取值,i表示本次计算时的取值,重复步骤32),否则进入步骤33);
33):若Tkn-T(k+i)n≥δ,其中约束条件为,1≤i<K,则分段动态聚焦的修正参数R(k+i)n=1,T(k+i)n=Tkn-δ,k'=k+i,k'表示下次计算时的取值,k表示本次计算时的取值,然后重新初始化i=1,程序回到步骤32)。
上述步骤中,δ为允许聚焦延迟误差,Tkn为第k分段、第n通道的聚焦延时参数,Rkn为第k分段、第n通道的分段动态聚焦的修正参数。
5.根据权利要求1所述的基于误差约束的动态聚焦最优分段方法及装置,其特征在于:所述步骤4)具体包括以下子步骤:
41)根据聚焦延时参数Tkn得到第一组聚焦延时参数T1n,将T1n装载到各个通道对应的递减器中,然后将步骤3)得到的分段动态聚焦的修正参数Rkn依次写入双口RAM中;
42)将步骤2)得到总共K段的分段位置序号Pk依次写入地址发生器中。
6.根据权利要求5所述的基于误差约束的动态聚焦最优分段方法及装置,其特征在于:所述步骤5)具体包括以下子步骤:
51)当分段动态聚焦开始时,在聚焦时钟的作用下,地址发生器内部的计数器开始计数工作,当计数器等于Pk时,地址发生器输出为k,从双口RAM中读出第k段聚焦区间的分段动态聚焦修正参数Rkn
52)地址发生器产生一个时钟脉冲,在时钟脉冲的上升沿,聚焦延时参数Tkn按下式进行修正
Figure FDA0000373655090000031
从而实现了分段聚焦延时参数Tkn的实时生成。
7.使用权利要求1-6任意一项所述的基于误差约束的动态聚焦最优分段方法的装置,其特征在于:包括微控制器、双口RAM、地址发生器和递减器,所述微控制器的数据端分别与双口RAM的第一数据端和地址发生器的数据端连接,微控制器的地址端分别与双口RAM的第一地址端和地址发生器的地址端连接,微控制器的第一片选端与双口RAM的第一片选端,微控制器的第二片选端与地址发生器的第一片选端连接,微控制器的写控制端分别与双口RAM的写控制端和地址发生器的写控制端连接,微控制器的时钟端分别与双口RAM的第一时钟端和地址发生器的第一时钟端连接,微控制器的IO口与地址发生器的使能端连接;所述双口RAM输出的数据位分别输入到每个递减器中,地址发生器产生的时钟脉冲分别输入到每个递减器中,地址发生器的读控制端分别与双口RAM的读控制端和第二片选端连接,时钟信号分别输入到双口RAM的第二时钟端和地址发生器的第二时钟端,地址发生器输出的CNT端与双口RAM的第二地址端连接。
8.根据权利要求7所述的装置,其特征在于:动态聚焦最优分段方法包括以下步骤:
1)微控制器根据给定的探测深度、聚焦通道数和采样率生成聚焦延时表;
2)结合给定的聚焦误差δ对整个聚焦延时表进行动态分段,并计算出分段位置数据和分段聚焦延时数据;
3)根据步骤2)生成分段聚焦延时数据的修正信息;
4)微控制器装载第一组聚焦延时参数到各个通道对应的递减器中,然后将分段动态聚焦的修正参数依次写入双口RAM中;将分段动态聚焦的分段位置数据依次写入地址发生器中;
5)根据聚焦延时修正信息实时修正聚焦延时参数。
CN201310381192.0A 2013-08-28 2013-08-28 基于误差约束的动态聚焦最优分段方法及装置 Active CN103425894B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310381192.0A CN103425894B (zh) 2013-08-28 2013-08-28 基于误差约束的动态聚焦最优分段方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310381192.0A CN103425894B (zh) 2013-08-28 2013-08-28 基于误差约束的动态聚焦最优分段方法及装置

Publications (2)

Publication Number Publication Date
CN103425894A true CN103425894A (zh) 2013-12-04
CN103425894B CN103425894B (zh) 2016-08-10

Family

ID=49650623

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310381192.0A Active CN103425894B (zh) 2013-08-28 2013-08-28 基于误差约束的动态聚焦最优分段方法及装置

Country Status (1)

Country Link
CN (1) CN103425894B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040189499A1 (en) * 2003-03-28 2004-09-30 Ho-San Han Sigma delta beamformer and method with reduced artifact
JP2004321647A (ja) * 2003-04-28 2004-11-18 Aloka Co Ltd 超音波診断装置
CN101907706A (zh) * 2009-06-05 2010-12-08 重庆博恩富克医疗设备有限公司 一种实现超声波逐点聚焦的方法及装置
US20120083695A1 (en) * 2005-07-22 2012-04-05 Zonare Medical Systems, Inc. Continuous transmit focusing method and apparatus for ultrasound imaging system
CN102727256A (zh) * 2012-07-23 2012-10-17 重庆博恩富克医疗设备有限公司 基于虚拟阵元的双聚焦波束合成方法及装置
CN102895000A (zh) * 2012-11-06 2013-01-30 重庆大学 一种基于自适应加权的双聚焦波束合成方法
US20130190594A1 (en) * 2012-01-23 2013-07-25 Alexander A. Oraevsky Scanning Optoacoustic Imaging System with High Resolution and Improved Signal Collection Efficiency

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040189499A1 (en) * 2003-03-28 2004-09-30 Ho-San Han Sigma delta beamformer and method with reduced artifact
JP2004321647A (ja) * 2003-04-28 2004-11-18 Aloka Co Ltd 超音波診断装置
US20120083695A1 (en) * 2005-07-22 2012-04-05 Zonare Medical Systems, Inc. Continuous transmit focusing method and apparatus for ultrasound imaging system
CN101907706A (zh) * 2009-06-05 2010-12-08 重庆博恩富克医疗设备有限公司 一种实现超声波逐点聚焦的方法及装置
US20130190594A1 (en) * 2012-01-23 2013-07-25 Alexander A. Oraevsky Scanning Optoacoustic Imaging System with High Resolution and Improved Signal Collection Efficiency
CN102727256A (zh) * 2012-07-23 2012-10-17 重庆博恩富克医疗设备有限公司 基于虚拟阵元的双聚焦波束合成方法及装置
CN102895000A (zh) * 2012-11-06 2013-01-30 重庆大学 一种基于自适应加权的双聚焦波束合成方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王平 等: "数字多波束逐点聚焦延时参数的压缩算法", 《重庆大学学报》, vol. 35, no. 9, 30 September 2012 (2012-09-30), pages 99 - 105 *
王平 等: "超声成像系统的动态孔径控制方法", 《重庆大学学报》, vol. 34, no. 3, 31 March 2011 (2011-03-31), pages 89 - 93 *

Also Published As

Publication number Publication date
CN103425894B (zh) 2016-08-10

Similar Documents

Publication Publication Date Title
CN103300890B (zh) 用于测量组织机械特性的系统及方法
US6231511B1 (en) Ultrasonic signal focusing method and apparatus for ultrasonic imaging system
CN103837608B (zh) 一种相控阵接收动态聚焦补偿方法及系统
CN105559828B (zh) 血流成像方法及系统
CN110927254B (zh) 一种基于fpga实现的高帧率超声全聚焦成像系统
CN105841645B (zh) 一种基于电磁超声的测厚方法
CN105832366A (zh) 一种用于波束合成过程中的延时实时计算方法
CN113504306B (zh) 一种基于超声相控阵低秩矩阵恢复的钢轨缺陷检测方法
CN107655533A (zh) 一种基于回波能量积分的气体超声波流量计信号处理方法和系统
WO2019187649A1 (ja) 超音波診断装置および超音波診断装置の制御方法
CN108186050B (zh) 一种基于超声通道数据的多普勒血流速度成像方法和系统
CN101907706B (zh) 一种实现超声波逐点聚焦的方法及装置
KR101652727B1 (ko) 초음파 도플러 영상 장치 및 그 제어 방법
Handegard et al. Tracking individual fish from a moving platform using a split-beam transducer
CN110353729A (zh) 一种基于双向长短期记忆网络的声波渡越时间检测方法
CN101961251B (zh) 一种医学超声诊断系统中实时计算变迹曲线的方法及装置
CN103425894A (zh) 基于误差约束的动态聚焦最优分段方法及装置
US20170322061A1 (en) Measurement apparatus, measurement method, and measurement program
CN106289124B (zh) 一种超声测厚闸门的实时跟随回波的方法
CN107643434A (zh) 一种基于分段切比雪夫距离的复杂波形触发方法
CN105842347A (zh) 一种时分复用的大功率相控阵超声信号发生装置
US9390510B2 (en) Method, apparatus, and system for outputting information of blood flow included in region of interest based on selection information
Drira et al. Image source detection for geoacoustic inversion by the Teager-Kaiser energy operator
CN110507361B (zh) 剪切波成像方法及系统
CN1439898A (zh) 基于线性插值的波束合成器及其合成方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C41 Transfer of patent application or patent right or utility model
CB03 Change of inventor or designer information

Inventor after: Pan Xiubao

Inventor after: Lin Nan

Inventor after: Luo Hanwu

Inventor after: Wang Ping

Inventor after: Sun Guang

Inventor after: Li Mengke

Inventor after: Cui Shigang

Inventor before: Wang Ping

Inventor before: Wang Siqi

Inventor before: Wang Linhong

Inventor before: Zou Qiangxin

Inventor before: Cheng Na

Inventor before: Pan Zhen

Inventor before: Gong Zhihui

COR Change of bibliographic data
TR01 Transfer of patent right

Effective date of registration: 20161114

Address after: The western part of Tongliao city the Inner Mongolia Autonomous Region 028000 Huolinhe Street No. 2080

Patentee after: STATE GRID EAST INNER MONGOLIA ELECTRIC POWER COMPANY LIMITED

Patentee after: Chongqing University

Address before: 400044 Shapingba street, Shapingba District, Chongqing, No. 174

Patentee before: Chongqing University