CN102764142B - 基于双层相关函数法的正交弱相关双极性编码激励方法 - Google Patents

基于双层相关函数法的正交弱相关双极性编码激励方法 Download PDF

Info

Publication number
CN102764142B
CN102764142B CN201210268006.8A CN201210268006A CN102764142B CN 102764142 B CN102764142 B CN 102764142B CN 201210268006 A CN201210268006 A CN 201210268006A CN 102764142 B CN102764142 B CN 102764142B
Authority
CN
China
Prior art keywords
coded
signal
correlation function
amplitude
coding
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.)
Expired - Fee Related
Application number
CN201210268006.8A
Other languages
English (en)
Other versions
CN102764142A (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.)
Harbin Hangkong Bochuang Science & Technology Co Ltd
Original Assignee
Harbin Hangkong Bochuang Science & Technology 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 Harbin Hangkong Bochuang Science & Technology Co Ltd filed Critical Harbin Hangkong Bochuang Science & Technology Co Ltd
Priority to CN201210268006.8A priority Critical patent/CN102764142B/zh
Publication of CN102764142A publication Critical patent/CN102764142A/zh
Application granted granted Critical
Publication of CN102764142B publication Critical patent/CN102764142B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8959Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes

Abstract

基于双层相关函数法的正交弱相关双极性编码激励方法,涉及一种透射式超声断层成像的编码激励方法。为了解决由于各阵元发出的超声信号相互影响引入的超声参数误差的问题,本发明的编码激励方法包括如下步骤:根据所需要的同时发射阵元数生成编码激励所需要的正交弱相关双极性编码集;使用生成的编码集对断层成像的发射过程进行调制,并采集穿过被测物体的波形信号;使用双层相关函数法从接收到的波形信号中提取渡越时间和幅值。本发明与目前不使用编码激励的方法相比,可以大幅度减小测量时间;与传统的基于匹配滤波的编码激励方法相比,该方法可以最大限度减小各阵元发射波形之间的影响,编码激励过程不影响参数提取的精度。

Description

基于双层相关函数法的正交弱相关双极性编码激励方法
技术领域
[0001] 本发明涉及一种透射式超声断层成像的编码激励方法,具体涉及一种对超声信号发射过程进行调制和从接收到的信号中提取渡越时间和超声幅值的方法。
背景技术
[0002] 发生癌变的乳腺肿瘤组织与正常的乳腺组织在声速、声吸收系数等声学参数方面存在差异,基于声学参数的乳腺癌诊断研究一直是超声医学的一个重要研究方向。癌变组织的密度、声速和声衰减系数均高于正常组织,此外,具有较高乳腺密度的妇女患乳腺癌的可能性较高,而乳腺组织的密度和声速有很强的相关性,因此能够测量乳腺组织声速、衰减系数的透射式超声断层成像在乳腺癌的诊断和预防方面有巨大的应用前景。
[0003] 作为一种医学成像设备,其测量过程需要患者的参与,测量时间作为影响患者舒适度的重要指标是该技术能否从实验室走向临床应用的重要因素。此外,由于乳腺处于胸腔的外侧,心跳和呼吸都会造成乳腺位置和形状发生变化,影响测量精度和诊断效果,因此透射式超声断层成像设备进行数据采集时要求患者在整个测量时间内屏住呼吸并且每一个切片的测量周期要小于心动周期,通过实验给出的整体测量时间和每一个切片的测量时间分别不超过I分钟和0.1秒,可见测量时间不仅关系到患者的舒适程度,也是决定透射式超声CT设备能否进入临床使用的一个关键指标。
[0004] 编码激励是减小透射式超声断层成像的测量时间的有效方法,但是传统的编码激励方法主要用于雷达、B超等反射式成像设备中,其原理是不同的能量发射装置发射不同的波形,然后使用匹配滤波 器对接收到的信号进行滤波以减小不同发射装置发射信号的相互影响。但由于不存在有限长的、互相关函数恒为零的有限支撑非零函数集,直接使用这些编码激励方法并不能完全消除各发射信号之间的相互影响。
发明内容
[0005] 本发明的目的在于提出一种基于双层相关函数法的正交弱相关双极性编码激励方法,它解决了由于各阵元发出的超声信号相互影响引入的超声参数误差的问题,从而在不影响参数提取精度的同时减小测量时间。
[0006] 本发明的基于双层相关函数法的正交弱相关双极性编码激励方法包括如下步骤:
[0007] 步骤一:根据所需要的同时发射阵元数生成编码激励所需要的正交弱相关双极性编码集;
[0008] 步骤二:使用生成的编码集对断层成像的发射过程进行调制,并采集穿过被测物体的波形信号;
[0009] 步骤三:使用双层相关函数法从接收到的波形信号中提取渡越时间和幅值。
[0010] 本发明提供了一种适用于透射式超声断层成像的编码激励方法,与目前不使用编码激励的方法相比,可以大幅度减小测量时间;与传统的基于匹配滤波的编码激励方法相t匕,该方法可以最大限度减小各阵元发射波形之间的影响,编码激励过程不影响参数提取的精度。
附图说明
[0011] 图1为基于双层相关函数法的正交弱相关双极性编码激励方法的流程图;
[0012] 图2为双层相关函数法流程图;
[0013]图3为原子信号波形/CO ;
[0014] 图4为编码信号波形F4 © ;
[0015]图5为编码信号波形&⑶;
[0016] 图6为到达时间差小于一个原子波形长度时的接收信号波形KO ;
[0017] 图7为到达时间差小于一个原子波形长度时的接收信号的互相关函数,其中为实线为馬&(»,虚线为&_^(» ;
[0018] 图8为到达时间差大于一个原子波形长度时的接收信号波形KO ;
[0019] 图9为到达时间差大于于一个原子波形长度时的接收信号与编码信号的互相关
函数,其中为实线为馬4®,虚线为;
[0020] 图10为到达时间差大于一个原子波形长度时原子信号与接收信号的相关函数
KA£y。
具体实施方式
[0021] 具体实施方式一:本实施方式的基于双层相关函数法的正交弱相关双极性编码激励方法是通过以下技术方案实现的:首先使用正交弱相关双极性编码对透射式超声断层成像的超声发射过程进行调制,在接收到超声信号后使用双层相关函数法进行参数的提取。
[0022] 如图1所示,具体步骤如下:
[0023] 步骤一:生成正交弱相关双极性编码集。正交弱相关双极性编码集的具体描述如下。
[0024] M个长度为I:的序列组成的集合
[0025]
Figure CN102764142BD00041
[0026] 定义两个序列的互相关为
[0027]
Figure CN102764142BD00042
[0028] 如果满足以下条件就被称为正交编码集:
[0029] 1、码元值只能取+1或者-1两个值,且规定每个序列的第一个码元必须是+1。
[0030] 2、集合中任意两序列是正交的,SP
[0031]
Figure CN102764142BD00043
(3);
[0032]其中5(-)为 Duric 算子。[0033] 3、各编码序列满足,对于任意整数a和非O整数》,有
Figure CN102764142BD00051
[0035] 编码集的生成方法为:在确定所需要的同时发射通道数以后,按照码元长度逐步递增的顺序列出每种长度下所有的编码,然后采用穷举法判断能否找到满足上述条件的编码集,不断增加编码的长度直到找到满足要求的编码集为止。
[0036] 步骤二:使用生成的编码集对断层成像的发射过程进行调制,调制方法与普通的双极性编码相同,编码信号与原子信号的数学关系如公式(5)所示,并采集穿过被测物体的波形信号,接收信号的数学描述见公式(6)。
[0037] 设不使用编码激励方法时发射的波形为原子波形/CO,其不为零的长度为T,则各发射阵元按照时间间隔力T根据编码的正负号发射正负原子波形/(«,单个阵元发出的波形为
Figure CN102764142BD00052
[0039] 其中Ιϊ(-)为延时算子,表示信号延后2*T ;4为公式I中给出编码。
[0040] 假设每次有两个阵元同时发射,均从时刻O开始发射,所使用的编码分别为4和A2,则任意位置处接收的信号可表示
Figure CN102764142BD00053
[0042] 其中#4©为使用编码4生成的编码信号,巧和%为对应的延时和波形幅值,
为使用编码為生成的编码信号,巧和皂为对应的延时和波形幅值。
[0043] 步骤三:使用双层相关函数法从接收到的超声信号提取渡越时间和幅值参数。双层相关函数法的具体步骤如下:首先假设接收到的超声信号中两个编码信号的到达时间差小于一个原子信号的长度,计算两个编码信号的渡越时间和幅值,然后对上述假设进行验证,如不符,使用原子信号与接收到的超声信号的相关函数确定渡越时间较小的编码信号的渡越时间和幅值,然后确定另一个编码信号的参数。双层相关函数法的流程图如图2所示,具体描述如下:
[0044] I)计算接收到的信号与编码信号的互相关函数,计算公式为:
Figure CN102764142BD00054
[0046] 令公式(7)取得最大值的时刻V即为相应的渡越时间,对应的幅值可表示为:
Figure CN102764142BD00055
[0048] 分别使用两个编码信号进行计算,可得到两个渡越时间和幅值组和(r2'Λ);
[0049] 2)从接收信号中减去由幅值较大的参数组确定的编码信号,重新计算另一组的渡越时间和幅值。假设则从接收信号中减去参数组(¾:¼)代表的编码信号并重新
计算另一个编码信号的参数得到
[0050] 3)如果计算出的传播时间差Ilr2.-1^l超过一个原子信号长度y,转第5步;
[0051] 4)从接收信号中减去计算出的编码信号,并求取残差的能量与接收信号能量的比值,如果不超过设定的阈值,则计算过程结束,(ϊ.,ίΟ和即为最终结果;
[0052] 5)计算原子信号与接收信号的相关函数,计算公式为:
[0053]
Figure CN102764142BD00061
[0054] 搜索第一个超过指定阈值的点,在该点前后一个原子波形长度范围内搜索最大值点,对应的时间和函数值即为先到达编码信号的渡越时间和幅值(5» ;
[0055] 6)分别假设第一个和第二个为先到达的编码信号,从接收信号中减去该编码信号,并使用残差与另一个编码信号的互相关函数求取另一个编码信号的到达时间和幅值,并计算残差,选择残差能量小的假设作为最终结果。
[0056] 在编码激励方法使用前,首先应使用常规方法进行一次数据采集,得到计算需要使用的原子信号的波形,然后就可以使用设计出的编码对发射过程进行调制,最后使用双层相关函数法提取超声参数。使用过程中需要注意,如果三个编码信号发生混叠,且三者的到达时间差有两个小于原子信号的持续时间时,双层相关函数法将会失效,因此使用时应该控制编码的数量以保证不会有三个信号同时到达。
[0057] 具体实施方式二:本实施方式以使用环形阵列的透射式超声断层成像的数据采集为例,对本发明的基于双层相关函数法的正交弱相关双极性编码激励方法进一步解释说明,该方法可以应用在任意的透射式断层成像中。
[0058] 假设超声阵元每次发出的信号可近似为高斯包络信号:
Figure CN102764142BD00062
[0060] 其中a为发射波形幅值,;T为发射时间,为衰减因子,石为探头中心频率。在|ί-Ι.|>Κ时,由于指数衰减的原因,信号非常小,可以近似为0,本文中取*:=;!】。参数I1、!,0.和/0由发射电路和换能器特性决定,由于原子波形只要满足紧支撑这个条件即可,因此参数值对算法没有影响,本文中使用的参数值为、. = 0二Miz,此时原子波形的长度为:
[0061]
Figure CN102764142BD00063
[0062] 即3.5us,原子信号波形如图3所示。但是,应当记住,本发明中对原子信号的波形没有限制,只要要满足紧支撑这一条件即可。[0063] 使用的编码为:
[0064]
Figure CN102764142BD00071
[0065] 可以验证上述编码满足前述编码集的条件,使用这两个编码得到的编码后信号波形分别如图4和5所示。
[0066] 当两个编码信号的到达时间相差不到一个原子信号长度时,如τ1:O.1ms ,r = 0_102ms,此时Δτ =0.002ms= 2us <T = 3 .5us,此时接收信号以及接收信号与编码信号
的互相关函数分别如图6和7所示,可以看出ryag(t)和rfag(t)分别在t=τ1和t=τ2处取得最大值6。此时使用双层相关函数法,第一步中计算出的幅值和渡越时间分别为b1' =1、τ1'=0.1ms和b2''=1,τ2''=0.102ms,由于幅值相同;第二步任选一个作为幅值较大的编码
信号,重新计算结果不变;第三步中可以判断出到达时间差小于一个原子信号的长度;第四步判断出残差为0,计算过程结束。
[0067] 当两个编码信号的到达时间相差超过一个原子信号长度时,如=
τ 1= 0.1ms ,τ2=0.10678ma,此时 Δτ=0.00678ms = 6.78us>T=3.5us, 此时接收信号以及接收信
号与编码信号的互相关函数分别如图8和9所示,计算出的参数值为t1' =0.10003ms 、
t2'=0.10675ms 和
Figure CN102764142BD00072
,我们认为b1'>b2'.重新计算后得到t2= O.10675ms,因
此渡越时间差超过一个原子信号的长度T,直接使用相关函数法求解会存在误差。原子信号与接收信号的相关函数如图10所示,t1 = 0.1ms恰好是它的一个极大值点,且是区间[τi'-Tτi+Tl内的极大值点。因此只要通过阈值等方法找到τi'的大致位置且误差在一个
原子信号长度内,就可以找到编码信号的精确时间和幅值,通过原子信号与接收信号的相关函数得到精确的传播时间和幅值是可行。 [0068] 上述接收信号中两个编码信号的幅值是相同的,若两者幅值差别较大不可忽略时,如果时间差小于一个原子长度,则使用相关函数法得到的两个幅值相差较大,此时可认为幅值较大的为准确值,从接收信号减去其对应的编码信号后再使用相关函数法计算幅值较小的编码信号的参数。如果时间差大于一个原子的长度,一般情况下传播时间较短的信号衰减较小,不会对结果造成影响。在极特殊情况下,先到的编码信号远小于后到的编码信号时,也可以选择使用原子信号匹配后到编码信号的最后一个编码,此时要特别注意最后一个编码的系数是+1还是_1。
[0069] 使用此实施方案时,每次可以有两个阵元同时向外发射超声信号,因此所需要的发射次数为原来的一半,由于发射的间隔主要取决于超声在被测物体中的传播时间,因此理论上这种方案可以把测量时间降低为原来的一半。考虑到使用编码激励方法时,编码波形的长度略长于原子波形,实际中的节省时间略小于一半,具体比例与超声传播时间和原子信号长度有关。
[0070] 本实施方案中只使用了两个编码,因此实际节省测量时间比例约为二分之一,如果使用包含更多编码的编码集,会 进一步提高节约时间的比例。

Claims (2)

1.基于双层相关函数法的正交弱相关双极性编码激励方法,其特征在于所述编码激励方法适用于透射式超声断层成像,包括如下步骤: 步骤一:根据所需要的同时发射阵元数生成编码激励所需要的正交弱相关双极性编码集,其中设编码集力Af个长度为I:的序列组成的集合:
Figure CN102764142BC00021
定义两个序列的互相关为:
Figure CN102764142BC00022
正交弱相关双极性编码集满足如下要求: (1)码元值只能取+1或者-1两个值,且规定每个序列的第一个码元必须是+1; (2)集合中任意两序列是正交的,即(O)
Figure CN102764142BC00023
其中尺/ —J.)为 Duric 算子; (3)对于任意整数W和非O整数BI,有
Figure CN102764142BC00024
步骤二:使用生成的编码集对断层成像的发射过程进行调制,并采集穿过被测物体的波形信号:
Figure CN102764142BC00025
,其中(0为使用编码 生成的编码信号,T1和&为对应的延时和波形幅值,Fa、(f)为使用编码J2生成的编码信号,T2 为对应的延时和波形幅值; 步骤三:使用双层相关函数法从接收到的波形信号中提取渡越时间和幅值,如果接收到的波形中两个编码信号的时间差小于一个原子信号的长度,使用编码信号和接收信号的相关函数计算编码信号的幅值和渡越时间;如果接收到的波形中两个编码信号的时间差大于一个原子信号的长度,使用原子和接收信号的相关函数计算首先到达的编码信号的幅值和渡越时间,然后计算另一个编码信号的幅值和渡越时间。
2.根据权利要求1所述的基于双层相关函数法的正交弱相关双极性编码激励方法,其特征在于所述正交弱相关双极性编码集的生成方法为:在确定所需要的同时发射阵元数以后,按照码元长度逐步递增的顺序列出每种长度下所有的编码,然后采用穷举法判断能否找到满足(I) - (3)条件的编码集,不断增加编码的长度直到找到满足要求的编码集为止。
CN201210268006.8A 2012-07-31 2012-07-31 基于双层相关函数法的正交弱相关双极性编码激励方法 Expired - Fee Related CN102764142B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210268006.8A CN102764142B (zh) 2012-07-31 2012-07-31 基于双层相关函数法的正交弱相关双极性编码激励方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210268006.8A CN102764142B (zh) 2012-07-31 2012-07-31 基于双层相关函数法的正交弱相关双极性编码激励方法

Publications (2)

Publication Number Publication Date
CN102764142A CN102764142A (zh) 2012-11-07
CN102764142B true CN102764142B (zh) 2014-02-12

Family

ID=47091885

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210268006.8A Expired - Fee Related CN102764142B (zh) 2012-07-31 2012-07-31 基于双层相关函数法的正交弱相关双极性编码激励方法

Country Status (1)

Country Link
CN (1) CN102764142B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10092192B2 (en) 2014-12-24 2018-10-09 Bahman LASHKARI Methods for generating multiple mismatched coded excitation signals
CN105973752B (zh) * 2016-04-22 2019-10-01 上海理工大学 集成透射式和反射式超声过程层析成像方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101271154A (zh) * 2008-03-28 2008-09-24 天津大学 用于消除机器人超声测距系统串扰的超声激励方法
CN101637395A (zh) * 2008-03-13 2010-02-03 超声成像 一种用于超声波合成成像的方法和装置
CN101869485A (zh) * 2010-06-23 2010-10-27 深圳大学 超声成像方法及装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004108826A (ja) * 2002-09-13 2004-04-08 Mie Tlo Co Ltd 超音波を用いた距離測定方法および距離測定装置
US20070239002A1 (en) * 2005-12-28 2007-10-11 Alam Sheikh K Superfast, High-Resolution Ultrasonic Imaging Using Coded Excitation

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101637395A (zh) * 2008-03-13 2010-02-03 超声成像 一种用于超声波合成成像的方法和装置
CN101271154A (zh) * 2008-03-28 2008-09-24 天津大学 用于消除机器人超声测距系统串扰的超声激励方法
CN101869485A (zh) * 2010-06-23 2010-10-27 深圳大学 超声成像方法及装置

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
JP特开2004-108826A 2004.04.08
Lizhen Wu et al.Novel Ultrasonic Ranging Approach Based on Spread Spectrum Modulation.《IEEE 2006 International Conference on Communication Technology》.2006,
Novel Ultrasonic Ranging Approach Based on Spread Spectrum Modulation;Lizhen Wu et al;《IEEE 2006 International Conference on Communication Technology》;20061130;全文 *
刘凯,高上凯.编码激励超声成像系统中二进制最优编码序列的研究.《中国生物医学工程学报》.2007,第26卷(第1期), *
姚振静.超声扩频测距脉冲编码激励与渡越时间估计.《中国博士学位论文全文数据库工程科技Ⅱ辑》.2011,(第7期),
超声扩频测距脉冲编码激励与渡越时间估计;姚振静;《中国博士学位论文全文数据库工程科技Ⅱ辑》;20110524(第7期);全文 *

Also Published As

Publication number Publication date
CN102764142A (zh) 2012-11-07

Similar Documents

Publication Publication Date Title
CN102551791B (zh) 一种超声成像方法和装置
US4249539A (en) Ultrasound needle tip localization system
CN102469980B (zh) 空间上精细的剪切波分散超声振动测定采样
CA2249030A1 (en) Catheter tracking system
CN104739448A (zh) 一种超声成像方法及装置
CN102764142B (zh) 基于双层相关函数法的正交弱相关双极性编码激励方法
KR20130076031A (ko) 결정 데이터에 기초하여 벡터 도플러 영상을 제공하는 초음파 시스템 및 방법
CN104407340A (zh) 拖曳线列阵阵形标定装置及方法
CN103529436A (zh) 基于hht的无接触生命探测中呼吸和心跳信号的分离及时频分析方法
CN104622507B (zh) 弹性模量测量方法和系统
CN109330597A (zh) 一种基于信道状态信息的人体呼吸追踪方法
CN104406642B (zh) 一种时差法超声波流量计精确测量方法
US8574158B2 (en) Time of flight estimation method using beamforming for acoustic tomography
CN105249994A (zh) 超声骨密度检测设备
US8162838B2 (en) Doppler signal spectrum calculating method and apparatus
Wang et al. Reconstruction of microwave absorption properties in heterogeneous tissue for microwave-induced thermo-acoustic tomography
CN105073021A (zh) 超声波诊断装置和超声波诊断方法
CN104146729B (zh) 一种骨骼超声系统中编码增强的聚焦超声骨组织微观结构检测方法
US9535041B2 (en) Propagation rate measurement device, propagation rate measurement program, and propagation rate measurement method
Solberg et al. Experimental investigation into radar-based central blood pressure estimation
RU2559420C1 (ru) Способ просветного акустического зондирования легких
Burov et al. Reconstruction of the blood flow velocity vector as an additional mode of ultrasonic tomograph
CN105919624A (zh) 一种高脉冲重复扫查频率信号的编码、成像方法和装置
KR20150130093A (ko) 횡파 속도를 추정하는 방법 및 시스템과 이를 이용한 조직 내 병변 진단 방법 및 시스템
US4817618A (en) Range-continuity anti-aliasing

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140212

Termination date: 20180731

CF01 Termination of patent right due to non-payment of annual fee