CN103549978A - 一种瞬时弹性成像快速时移估计方法 - Google Patents
一种瞬时弹性成像快速时移估计方法 Download PDFInfo
- Publication number
- CN103549978A CN103549978A CN201310553321.XA CN201310553321A CN103549978A CN 103549978 A CN103549978 A CN 103549978A CN 201310553321 A CN201310553321 A CN 201310553321A CN 103549978 A CN103549978 A CN 103549978A
- Authority
- CN
- China
- Prior art keywords
- signal
- time shift
- ultrasound echo
- sad
- echo signal
- 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
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种瞬时弹性成像快速时移估计方法,该方法包括以下步骤:接收不同时刻的N个超声回波信号s1,…,sN;根据预先设定的窗长和重叠率,将所述N个超声回波信号s1,…,sN分为若干个信号片段,所述N个超声回波信号s1,…,sN上相同深度的信号片段分别表示为s1k,…,sNk,计算所述信号片段s2k相对所述超声回波信号信号s1的时移d2k,并以所述时移d2k作为先验值计算所述信号片段s3k相对所述超声回波信号s1的时移d3k,以此类推,直到以d(n-1)k作为先验值计算所述信号片段snk和所述超声回波信号s1的时移dnk。本发明可以有效地缩小时移估计的范围,提高计算速度。
Description
技术领域
本发明涉及快速时移估计方法,尤其涉及一种瞬时弹性成像的快速时移估计方法。
背景技术
瞬时弹性成像是一种实用的间接测量组织弹性的方法,该方法的原理是根据剪切波在组织里面的传播速度与组织剪切模量的理论关系,利用测得的剪切波的传播速度来计算组织的弹性系数。其通常的做法是采用外振动器在组织表面产生50-200Hz的低频脉冲振动,低频振动在组织中会产生剪切波,在软组织里面剪切波的传播速度只有几到几十米每秒,而在组织成像里面通常所用的超声信号是一种压缩波,其波速约为1540米每秒。因此可以用超声波跟踪低频振动产生的剪切波在组织中的传播过程。通过对超声M模式信号的分析,来确定低频振动产生的剪切波在组织中的传播速度,从而计算组织的弹性系数。
在上述瞬时弹性成像方法中,需要从剪切波引起的组织应变图上得到剪切波的传播速度,因此分析不同时刻的超声回波信号,来确定剪切波传播造成的组织偏移是瞬时弹性成像信号处理过程中的关键步骤之一。对于每次低频脉冲振动,发射接收N次超声波,根据这N次超声回波信号进行位移/应变估计的步骤为:1)确定一定的窗长度和重叠率;2)按照确定的窗长和重叠率,确定每一对对应射频信号片段的时移,即相对位移;3)利用得到的位移计算组织的应变分布。
在第二步中,确定对应射频信号片段的时移时,通常是以第一根射频信号为基准,即确定射频信号片段与第一根射频信号上的对应信号片段间的时移。有两种实现方式,一种是直接确定射频信号片段与第一根射频信号上的对应信号片段间的时移;另一种是先确定相邻两根射频信号对应信号片段的时移,再通过累加得到相对于第一根射频信号对应信号片段的时移。由于组织的位移是连续的,所以第二种方式搜索时移的范围比较小,计算量小,但是由于受信号采样频率的限制,会造成误差累积,影响精度。虽然可以通过插值来解决这一问题,但是通过插值来提高采样频率会增加计算量。而第一种方式没有误差累积问题,但是时移的搜索范围较第二种方式大,计算量也随之增大。
发明内容
本发明的目的是,针对上述直接计算射频信号片段相对于第一根射频信号的对应信号片段的时移计算量大的问题,提供一种用于瞬时弹性成像快速时移估计方法,利用组织运动连续性,缩小搜索范围,降低计算量,避免误差累积问题或插值带来的计算量增大问题。
为了实现上述目的,本发明提供了一种瞬时弹性成像快速时移估计方法,该方法包括以下步骤:
接收不同时刻的N个超声回波信号s1,…,sN;根据预先设定的窗长和重叠率,将所述N个超声回波信号s1,…,sN分为若干个信号片段,所述N个超声回波信号s1,…,sN上相同深度的信号片段分别表示为s1k,…,sNk,所述snk=sn(k:M1+k-1),其中信号片段snk为超声回波信号sn的第k至第M1+k-1个采样点,n=1,…,N,M1为窗长;计算所述信号片段s2k相对所述超声回波信号s1的时移d2k,并以所述时移d2k作为先验值计算所述信号片段s3k相对所述超声回波信号s1的时移d3k,以此类推,直到以d(n-1)k作为先验值计算所述信号片段snk和所述超声回波信号s1的时移dnk,其中n=2,…,N。
进一步地,所述计算所述信号片段s2k相对所述超声回波信号s1的时移d2k步骤包括:
分别计算所述信号片段s2k与所述超声回波信号s1(k+i:M+k+i-1)的绝对差和(Sum ofAbsolute Differences,简称SAD)值sad(i),i=-L,…,L,其中L是一个正整数,为时移估计的搜索范围,即
进一步地,所述以d(n-1)k作为先验值计算所述信号片段snk和所述超声回波信号s1的时移dnk步骤包括:
分另计算snk与s1(k+i:M+k+i-1)的SAD值sad(i),i=d(n-1)k-L,…,d(n-1)k+L,即
本发明可以有效地缩小时移估计的范围,提高计算速度。另外,因为直接估计了每次超声回波信号与第一次超声回波信号的时移,避免了误差累积。无需对信号进行插值,避免了计算量的增加。
附图说明
图1是本发明实施例提供的一种瞬时弹性成像快速时移估计方法流程图。
具体实施方式
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
图1为本发明实施例提供的一种瞬时弹性成像快速时移估计方法流程图。如图1所示,本发明实施例将接收不同时刻的N个超声回波信号s1,…,sN,根据预先设定的窗长和重叠率将N个超声回波信号s1,…,sN划分为若干信号片段,设信号s1,…,sN上相同深度的信号片段分别表示为s1k,…,sNk,长度都为M1(即窗长为M1),且snk=sn(k:M1+k-1),即信号片段snk为信号sn的第k至第M1+k-1个采样点,n=1,…,N。
本发明实施例首先计算信号片段s2k相对超声回波信号s1的时移d2k,并以时移d2k作为先验值计算信号片段s3k相对超声回波信号s1的时移d3k,以此类推,直到以d(n-1)k作为先验值计算所述信号片段snk和所述超声回波信号s1的时移dnk,其中n=2,…,N。
优选地,本发明实施例采用绝对差和(Sum of Absolute Differences,简称SAD)方法进行时移计算。
以下通过SAD方法计算信号片段s2k相对超声回波信号s1的时移d2k,其实现步骤包括:
1)分别计算信号片段s2k与超声回波信号s1(k+i:M+k+i-1)的SAD值sad(i),i=-L,…,L,其中L是一个正整数,为时移估计的搜索范围,即
以下以d(n-1)k作为先验值,用SAD方法估计信号段snk和超声回波信号s1的时移dnk,其实现步骤包括:
1)分别计算snk与s1(k+i:M+k+i-1)的SAD值sad(i),i=d(n-1)k-L,…,d(n-1)k+L,即
本发明实施例利用组织位移的连续性特征,以前一次时移估计值为先验值,每次估计时移的需要计算2L+1次SAD值,而如果不利用这个先验值,则搜索时移时需要计算SAD值的次数一定要远大于2L+1才能得到正确的时移估计。本发明实施例可以有效地缩小时移估计的范围,提高计算速度。另外,因为直接估计了每次超声回波信号与第一次超声回波信号的时移,避免了误差累积。无需对信号进行插值,避免了计算量的增加。
显而易见,在不偏离本发明的真实精神和范围的前提下,在此描述的本发明可以有许多变化。因此,所有对于本领域技术人员来说显而易见的改变,都应包括在本权利要求书所涵盖的范围之内。本发明所要求保护的范围仅由所述的权利要求书进行限定。
Claims (3)
1.一种瞬时弹性成像快速时移估计方法,其特征在于,包括以下步骤:
接收不同时刻的N个超声回波信号s1,…,sN;
根据预先设定的窗长和重叠率,将所述N个超声回波信号s1,…,sN分为若干个信号片段,所述N个超声回波信号s1,…,sN上相同深度的信号片段分别表示为s1k,…,sNk,所述snk=sn(k:M1+k-1),其中信号片段snk为超声回波信号信号sn的第k至第M1+k-1个采样点,n=1,…,N,M1为窗长;
计算所述信号片段s2k相对所述超声回波信号s1的时移d2k,并以所述时移d2k作为先验值计算所述信号片段s3k相对所述超声回波信号s1的时移dsk,以此类推,直到以d(n-1)k作为先验值计算所述信号片段snk和所述超声回波信号s1的时移dnk,其中n=2,…,N。
3.根据权利要求1所述的方法,其特征在于:所述以d(n-1)k作为先验值计算所述信号片段snk和所述超声回波信号s1的时移dnk步骤包括:
分另计算snk与s1(k+i:M+k+i-1)的SAD值sad(i),i=d(n-1)k-L,…,d(n-1)k+L,即
计算其中i=d(n-1)k-L,…,d(n-1)k+L,即找出sad(i)的最小值出现的位置d,所述信号片段snk相对s1的时移dnk=d。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310553321.XA CN103549978B (zh) | 2013-11-08 | 2013-11-08 | 一种瞬时弹性成像快速时移估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310553321.XA CN103549978B (zh) | 2013-11-08 | 2013-11-08 | 一种瞬时弹性成像快速时移估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103549978A true CN103549978A (zh) | 2014-02-05 |
CN103549978B CN103549978B (zh) | 2016-01-13 |
Family
ID=50004473
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310553321.XA Active CN103549978B (zh) | 2013-11-08 | 2013-11-08 | 一种瞬时弹性成像快速时移估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103549978B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103829975A (zh) * | 2014-03-17 | 2014-06-04 | 飞依诺科技(苏州)有限公司 | 超声弹性成像位移值的计算方法 |
CN105844645A (zh) * | 2016-03-31 | 2016-08-10 | 深圳大学 | 一种瞬时弹性成像编码检测方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6176827B1 (en) * | 1998-05-19 | 2001-01-23 | U.S. Philips Corporation | Method of detecting elasticity variations and echographic apparatus for carrying out the method |
WO2005120358A1 (ja) * | 2004-06-09 | 2005-12-22 | Hitachi Medical Corporation | 弾性画像表示方法及び超音波診断装置 |
CN102860842A (zh) * | 2012-09-26 | 2013-01-09 | 浙江大学 | 一种实时准静态超声弹性成像方法 |
CN103040488A (zh) * | 2012-12-21 | 2013-04-17 | 深圳大学 | 一种实时超声弹性成像位移估计方法和系统 |
-
2013
- 2013-11-08 CN CN201310553321.XA patent/CN103549978B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6176827B1 (en) * | 1998-05-19 | 2001-01-23 | U.S. Philips Corporation | Method of detecting elasticity variations and echographic apparatus for carrying out the method |
WO2005120358A1 (ja) * | 2004-06-09 | 2005-12-22 | Hitachi Medical Corporation | 弾性画像表示方法及び超音波診断装置 |
CN102860842A (zh) * | 2012-09-26 | 2013-01-09 | 浙江大学 | 一种实时准静态超声弹性成像方法 |
CN103040488A (zh) * | 2012-12-21 | 2013-04-17 | 深圳大学 | 一种实时超声弹性成像位移估计方法和系统 |
Non-Patent Citations (5)
Title |
---|
ZAHIRI-AZAR, R.等: "Motion Estimation in Ultrasound Images Using Time Domain Cross Correlation With Prior Estimates", 《BIOMEDICAL ENGINEERING, IEEE TRANSACTIONS ON》 * |
刘宝亮: "基于超声瞬时弹性成像技术的生物组织弹性研究", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑(月刊)》 * |
孙瑞超 等: "基于相移的实时超声弹性成像方法", 《中国生物医学工程学报》 * |
孙瑞超 等: "实时超声弹性成像原理与方法", 《中国医疗器械信息》 * |
邱广: "基于超声的组织弹性成像算法研究", 《中国优秀硕士学位论文 医药卫生科技辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103829975A (zh) * | 2014-03-17 | 2014-06-04 | 飞依诺科技(苏州)有限公司 | 超声弹性成像位移值的计算方法 |
CN105844645A (zh) * | 2016-03-31 | 2016-08-10 | 深圳大学 | 一种瞬时弹性成像编码检测方法及系统 |
CN105844645B (zh) * | 2016-03-31 | 2019-11-12 | 深圳大学 | 一种瞬时弹性成像编码检测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103549978B (zh) | 2016-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106970371B (zh) | 一种基于Keystone和匹配滤波的目标检测方法 | |
EP2648379B1 (en) | Method and apparatus for channel delay and phase error estimation | |
US9370339B2 (en) | Systems, methods and computer programs for detection of tissue mechanical property | |
CN105589066B (zh) | 一种利用垂直矢量阵估计水下匀速运动航行器参数的方法 | |
CN104678435A (zh) | 一种提取Rayleigh面波频散曲线的方法 | |
Yang et al. | Frequency-varying group delay estimation using frequency domain polynomial chirplet transform | |
JP6024012B1 (ja) | 二点間の波動伝播時間の推定方法 | |
JP5599371B2 (ja) | 測位装置 | |
JP2016225981A (ja) | マルチパスの時間遅延の推定装置、推定方法及び受信機 | |
CN103323815A (zh) | 一种基于等效声速的水下声学定位方法 | |
CN104849708A (zh) | 基于频域多项式相位变换的高速机动目标参数估计方法 | |
CN108459087A (zh) | 一种应用于板结构损伤检测的多模态Lamb波模态分离方法 | |
CN108205080A (zh) | 相干平均法谐波信号功率谱估计方法 | |
CN103549978A (zh) | 一种瞬时弹性成像快速时移估计方法 | |
CN103809170A (zh) | 一种多普勒频偏修正方法及系统 | |
CN105877783A (zh) | 二维剪切波弹性成像方法和装置 | |
CN105425298A (zh) | 一种消除有限差分正演过程中数值频散的方法和装置 | |
CN104731762B (zh) | 基于循环移位的立方相位信号参数估计方法 | |
KR20160024034A (ko) | 이동 노드의 위치 추정 방법 및 그 장치 | |
US20110231097A1 (en) | Acoustic radial profiling via frequency domain processing | |
JP2014211330A (ja) | 目標追尾装置及び目標追尾方法 | |
CN102624657B (zh) | 频差估计方法 | |
CN104135767A (zh) | 一种测定信号直达波到达时差的分段互相关方法 | |
Wu et al. | Parameter estimation for SAR moving target detection using fractional Fourier transform | |
KR101689628B1 (ko) | 고정 신호원의 위치 추정 장치 및 방법 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |