[go: up one dir, main page]

CN103735284B - Based on RF signal estimation method in the three-D ultrasonic elastogram of linear scanning - Google Patents

Based on RF signal estimation method in the three-D ultrasonic elastogram of linear scanning Download PDF

Info

Publication number
CN103735284B
CN103735284B CN201310737468.4A CN201310737468A CN103735284B CN 103735284 B CN103735284 B CN 103735284B CN 201310737468 A CN201310737468 A CN 201310737468A CN 103735284 B CN103735284 B CN 103735284B
Authority
CN
China
Prior art keywords
signal
estimated
phase
closest
amplitude
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
CN201310737468.4A
Other languages
Chinese (zh)
Other versions
CN103735284A (en
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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN201310737468.4A priority Critical patent/CN103735284B/en
Publication of CN103735284A publication Critical patent/CN103735284A/en
Application granted granted Critical
Publication of CN103735284B publication Critical patent/CN103735284B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了基于线性扫描的三维超声弹性成像中RF信号估计方法,包含以下顺序的步骤:首先由线性移动探头扫描得到的RF信号序列得到待估计RF信号的位置;然后取距离待估计RF信号最近的前后两条RF信号幅度的距离加权平均作为其幅度,取距离待估计RF信号最近的一条RF信号的相位作为其相位,利用其幅度和其相位计算得到待估计的RF信号。本发明的方法,可以解决因探头移动速率不同而导致采集的RF信号不均匀的问题,方法简单,实用性强。

The invention discloses a method for estimating RF signals in three-dimensional ultrasonic elastography based on linear scanning, which includes the following steps: first, the position of the RF signal to be estimated is obtained from the RF signal sequence scanned by a linear moving probe; and then the distance to the RF signal to be estimated is obtained The distance-weighted average of the amplitudes of the two closest RF signals before and after is taken as its amplitude, and the phase of the RF signal closest to the RF signal to be estimated is taken as its phase, and the RF signal to be estimated is calculated by using its amplitude and its phase. The method of the invention can solve the problem of inhomogeneous collected RF signals caused by different moving speeds of the probes, and the method is simple and practical.

Description

基于线性扫描的三维超声弹性成像中RF信号估计方法RF signal estimation method in 3D ultrasound elastography based on linear scan

技术领域technical field

本发明涉及涉及计算机辅助医学成像技术领域,特别涉及基于线性扫描的三维超声弹性成像中RF信号估计方法。The invention relates to the technical field of computer-aided medical imaging, in particular to a method for estimating RF signals in three-dimensional ultrasonic elastography based on linear scanning.

背景技术Background technique

超声弹性成像是一种新型超声诊断技术,能够对传统超声无法探测的肿瘤及扩散疾病进行硬度信息成像。在三维超声弹性成像中,一种简单的扫描方法是线性扫描(详见已公开的专利:一种用于三维超声弹性成像的扫查装置和方法,公开号:CN102908166A),这种扫描方法是通过将超声探头固定到带有位置传感器的直线型滑轨上,沿着滑轨移动探头扫描采集信号并记录探头的位置。Ultrasound elastography is a new ultrasonic diagnostic technology that can image the hardness information of tumors and diffuse diseases that cannot be detected by traditional ultrasound. In three-dimensional ultrasound elastography, a simple scanning method is linear scanning (see the published patent for details: a scanning device and method for three-dimensional ultrasound elastography, publication number: CN102908166A), this scanning method is By fixing the ultrasonic probe on a linear slide rail with a position sensor, moving the probe along the slide rail to scan and collect signals and record the position of the probe.

由于人手难以控制探头匀速运动,采集的信号在某些位置会比较密集而在另一些位置则很稀疏,而且通常还会出现,在某一位置处压缩前有信号而压缩后没采集到信号这种情况。在计算弹性图像时,如果仅仅根据压缩前和压缩后采集到的两组信号序列中位置最近邻的两条信号作为匹配信号,计算结果并不太准确。Because it is difficult for the human hand to control the uniform movement of the probe, the collected signal will be denser in some positions and sparse in other positions, and usually there will be a signal at a certain position before compression but no signal after compression. situation. When calculating the elastic image, if only the two nearest neighbor signals in the two sets of signal sequences collected before compression and after compression are used as matching signals, the calculation result is not very accurate.

因此,在线性扫描三维超声弹性成像中,若想得到更加准确的弹性图像,需要一种合理有效的RF信号估计方法。Therefore, in linear scanning 3D ultrasound elastography, a reasonable and effective RF signal estimation method is needed if more accurate elastic images are to be obtained.

发明内容Contents of the invention

本发明的目的在于克服现有技术的缺点与不足,提供基于线性扫描的三维超声弹性成像中RF信号估计方法,可用于解决压缩前、后信号的匹配问题,使成像效果更佳。The purpose of the present invention is to overcome the shortcomings and deficiencies of the prior art, and provide a method for estimating RF signals in three-dimensional ultrasonic elastography based on linear scanning, which can be used to solve the matching problem of signals before and after compression, and make the imaging effect better.

本发明的目的通过以下的技术方案实现:The purpose of the present invention is achieved through the following technical solutions:

基于线性扫描的三维超声弹性成像中RF信号估计方法,包含以下顺序的步骤:The RF signal estimation method in the three-dimensional ultrasound elastography based on linear scanning comprises the steps in the following sequence:

1)由线性移动探头扫描得到的RF信号序列得到待估计RF信号的位置;1) Obtain the position of the RF signal to be estimated from the RF signal sequence obtained by scanning the linear moving probe;

2)取距离待估计RF信号最近的前后两条RF信号幅度的距离加权平均作为其幅度,取距离待估计RF信号最近的一条RF信号的相位作为其相位,利用其幅度和其相位计算得到待估计的RF信号。2) Take the distance-weighted average of the two RF signal amplitudes closest to the RF signal to be estimated as its amplitude, take the phase of the RF signal closest to the RF signal to be estimated as its phase, and use its amplitude and phase to calculate the target estimated RF signal.

所述的步骤2)具体包含以下步骤:Described step 2) specifically comprises the following steps:

a、待估计RF信号的位置为x,取距离x最近的前后两条RF信号RF(a)和RF(b),其中a、b表示信号所在位置,a<x<b;a. The position of the RF signal to be estimated is x, and the two RF signals RF(a) and RF(b) closest to x are taken, where a and b indicate the position of the signal, a<x<b;

b、对RF(a)进行Hilbert变换,求得其幅度amp(a)和其相位ang(a),对RF(b)进行Hilbert变换,求得其幅度amp(b)和其相位ang(b);b. Perform Hilbert transform on RF(a), obtain its amplitude amp(a) and its phase ang(a), perform Hilbert transform on RF(b), obtain its amplitude amp(b) and its phase ang(b );

取二者幅度的距离加权平均为 A = b - x b - a . amp ( a ) + x - a b - a . amp ( b ) ; Taking the distance-weighted average of the magnitudes of the two is A = b - x b - a . amp ( a ) + x - a b - a . amp ( b ) ;

取二者中距x最近者的相位为φ=ang(i),i=a或b;Take the phase of the closest to x among the two as φ=ang(i), i=a or b;

采用公式RF(x)=A·cos(φ)计算出待估计的信号。The signal to be estimated is calculated using the formula RF(x)=A·cos(φ).

步骤1)中,所述的RF信号序列为压缩前的RF信号序列或压缩后的RF信号序列。In step 1), the RF signal sequence is a pre-compressed RF signal sequence or a compressed RF signal sequence.

本发明与现有技术相比,具有如下优点和有益效果:Compared with the prior art, the present invention has the following advantages and beneficial effects:

本发明所述方法简单适用,能够估计出相邻两条信号之间任意位置的信号,可以有效地解决线性扫描过程中因探头移动速率的变化导致采集到的信号不均匀的问题。另外,采用幅度距离加权平均可使待估计的信号更接近真实信息。The method of the invention is simple and applicable, can estimate the signal at any position between two adjacent signals, and can effectively solve the problem of inhomogeneous collected signals caused by changes in the moving speed of the probe during the linear scanning process. In addition, the weighted average of amplitude and distance can make the signal to be estimated closer to the real information.

附图说明Description of drawings

图1为本发明所述的基于线性扫描的三维超声弹性成像中RF信号估计方法的流程图;Fig. 1 is the flow chart of the RF signal estimation method in the three-dimensional ultrasound elastography based on linear scanning according to the present invention;

图2为图1所述方法的处理过程示意图;Fig. 2 is the processing schematic diagram of the method described in Fig. 1;

图3为图1所述方法的相邻位置的两条RF信号和位于这两条信号之间的待估计RF信号;Fig. 3 is two RF signals of adjacent positions of the method described in Fig. 1 and the RF signal to be estimated between these two signals;

图4为图1所述方法所得到的三维超声弹性图像。Fig. 4 is a three-dimensional ultrasonic elastic image obtained by the method described in Fig. 1 .

具体实施方式detailed description

下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。The present invention will be further described in detail below in conjunction with the embodiments and the accompanying drawings, but the embodiments of the present invention are not limited thereto.

如图1、2、3,基于线性扫描的三维超声弹性成像中RF信号估计方法,包含以下顺序的步骤:As shown in Figures 1, 2, and 3, the RF signal estimation method in 3D ultrasound elastography based on linear scanning includes the following sequential steps:

1)由线性移动探头扫描得到的RF信号序列得到待估计RF信号的位置x,其中RF信号序列为压缩前的RF信号序列或压缩后的RF信号序列;1) Obtain the position x of the RF signal to be estimated from the RF signal sequence scanned by the linear moving probe, wherein the RF signal sequence is the RF signal sequence before compression or the RF signal sequence after compression;

2)如图3,a、待估计RF信号的位置为x,取距离x最近的前后两条RF信号RF(a)和RF(b),其中a、b表示信号所在位置,a<x<b;2) As shown in Figure 3, a, the position of the RF signal to be estimated is x, take the two RF signals RF(a) and RF(b) closest to x, where a and b indicate the position of the signal, a<x< b;

b、对RF(a)进行Hilbert变换,求得其幅度amp(a)和其相位ang(a),对RF(b)进行Hilbert变换,求得其幅度amp(b)和其相位ang(b);b. Perform Hilbert transform on RF(a), obtain its amplitude amp(a) and its phase ang(a), perform Hilbert transform on RF(b), obtain its amplitude amp(b) and its phase ang(b );

取二者幅度的距离加权平均为 A = b - x b - a . amp ( a ) + x - a b - a . amp ( b ) ; Taking the distance-weighted average of the magnitudes of the two is A = b - x b - a . amp ( a ) + x - a b - a . amp ( b ) ;

取二者中距x最近者的相位为φ=ang(i),i=a或b;Take the phase of the closest to x among the two as φ=ang(i), i=a or b;

采用公式RF(x)=A·cos(φ)计算出待估计的信号。The signal to be estimated is calculated using the formula RF(x)=A·cos(φ).

由于每个信号帧包含许多条信号,为了估计出一帧信号,必须分别估计出每一条信号,每一条信号的估计方法都一样,重复步骤1)、2),直至估计出位于x处的完整一帧信号,最终得到最终的三维超声弹性图像,如图4,图像清晰可见,有效解决了因难以保持探头均匀移动导致的压缩前、后信号不匹配的问题。Since each signal frame contains many signals, in order to estimate a frame of signals, each signal must be estimated separately, and the estimation method of each signal is the same, repeating steps 1) and 2) until the complete signal at x is estimated One frame of signal, and finally get the final three-dimensional ultrasonic elastic image, as shown in Figure 4, the image is clearly visible, which effectively solves the problem of signal mismatch between before and after compression caused by the difficulty in keeping the probe moving uniformly.

上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。The above-mentioned embodiment is a preferred embodiment of the present invention, but the embodiment of the present invention is not limited by the above-mentioned embodiment, and any other changes, modifications, substitutions, combinations, Simplifications should be equivalent replacement methods, and all are included in the protection scope of the present invention.

Claims (3)

1.基于线性扫描的三维超声弹性成像中RF信号估计方法,其特征在于:包含以下顺序的步骤:1. RF signal estimation method in the three-dimensional ultrasound elastography based on linear scanning, it is characterized in that: comprise the step of following order: 1)由线性移动探头扫描得到的RF信号序列得到待估计RF信号的位置;1) Obtain the position of the RF signal to be estimated from the RF signal sequence obtained by scanning the linear moving probe; 2)取距离待估计RF信号最近的前后两条RF信号幅度的距离加权平均作为其幅度,取距离待估计RF信号最近的一条RF信号的相位作为其相位,利用其幅度和其相位计算得到待估计的RF信号。2) Take the distance-weighted average of the two RF signal amplitudes closest to the RF signal to be estimated as its amplitude, take the phase of the RF signal closest to the RF signal to be estimated as its phase, and use its amplitude and phase to calculate the target estimated RF signal. 2.根据权利要求1所述的基于线性扫描的三维超声弹性成像中RF信号估计方法,其特征在于:所述的步骤2)具体包含以下步骤:2. RF signal estimation method in the three-dimensional ultrasonic elastography based on linear scanning according to claim 1, is characterized in that: described step 2) specifically comprises the following steps: a、待估计RF信号的位置为x,取距离x最近的前后两条RF信号RF(a)和RF(b),其中a、b表示信号所在位置,a<x<b;a. The position of the RF signal to be estimated is x, and the two RF signals RF(a) and RF(b) closest to x are taken, where a and b indicate the position of the signal, a<x<b; b、对RF(a)进行Hilbert变换,求得其幅度amp(a)和其相位ang(a),对RF(b)进行Hilbert变换,求得其幅度amp(b)和其相位ang(b);b. Perform Hilbert transform on RF(a), obtain its amplitude amp(a) and its phase ang(a), perform Hilbert transform on RF(b), obtain its amplitude amp(b) and its phase ang(b ); 取二者幅度的距离加权平均为 Taking the distance-weighted average of the magnitudes of the two is 取二者中距x最近者的相位为φ=ang(i),i=a或b;Take the phase of the closest to x among the two as φ=ang(i), i=a or b; 采用公式RF(x)=A·cos(φ)计算出待估计的信号。The signal to be estimated is calculated using the formula RF(x)=A·cos(φ). 3.根据权利要求1所述的基于线性扫描的三维超声弹性成像中RF信号估计方法,其特征在于:步骤1)中,所述的RF信号序列为压缩前的RF信号序列或压缩后的RF信号序列。3. RF signal estimation method in the three-dimensional ultrasonic elastography based on linear scan according to claim 1, is characterized in that: in step 1), described RF signal sequence is the RF signal sequence before compression or RF after compression signal sequence.
CN201310737468.4A 2013-12-26 2013-12-26 Based on RF signal estimation method in the three-D ultrasonic elastogram of linear scanning Expired - Fee Related CN103735284B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310737468.4A CN103735284B (en) 2013-12-26 2013-12-26 Based on RF signal estimation method in the three-D ultrasonic elastogram of linear scanning

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310737468.4A CN103735284B (en) 2013-12-26 2013-12-26 Based on RF signal estimation method in the three-D ultrasonic elastogram of linear scanning

Publications (2)

Publication Number Publication Date
CN103735284A CN103735284A (en) 2014-04-23
CN103735284B true CN103735284B (en) 2015-10-28

Family

ID=50492410

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310737468.4A Expired - Fee Related CN103735284B (en) 2013-12-26 2013-12-26 Based on RF signal estimation method in the three-D ultrasonic elastogram of linear scanning

Country Status (1)

Country Link
CN (1) CN103735284B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113375597B (en) * 2021-06-28 2022-04-19 华中科技大学 3D laser scanning method and equipment for structural parts based on RF-NSGA-II

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1323570A (en) * 2000-05-15 2001-11-28 富士胶片株式会社 Method for producing image device, and ultrasonic diagnosis equipment using said method
CN1586408A (en) * 2004-08-20 2005-03-02 清华大学 Multile size biological tissue displacement evaluating method
CN1954235A (en) * 2003-10-08 2007-04-25 艾克提斯有源传感器有限公司 Improved method and device for local spectral analysis of an ultrasonic signal
CN102908166A (en) * 2012-08-28 2013-02-06 华南理工大学 Scanning device and method for three-dimensional ultrasound elasticity imaging
CN102920482A (en) * 2012-11-26 2013-02-13 重庆理工大学 Multi-frequency alternative-ejecting real-time ultrasonic elastography method
CN102940510A (en) * 2012-08-31 2013-02-27 华南理工大学 Automatic focusing method for ultrasonic elastography

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3734192B2 (en) * 1996-12-26 2006-01-11 株式会社日立メディコ Ultrasonic diagnostic equipment
US8254654B2 (en) * 2007-10-31 2012-08-28 University Of Southern California Sidelobe suppression in ultrasound imaging using dual apodization with cross-correlation

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1323570A (en) * 2000-05-15 2001-11-28 富士胶片株式会社 Method for producing image device, and ultrasonic diagnosis equipment using said method
CN1954235A (en) * 2003-10-08 2007-04-25 艾克提斯有源传感器有限公司 Improved method and device for local spectral analysis of an ultrasonic signal
CN1586408A (en) * 2004-08-20 2005-03-02 清华大学 Multile size biological tissue displacement evaluating method
CN102908166A (en) * 2012-08-28 2013-02-06 华南理工大学 Scanning device and method for three-dimensional ultrasound elasticity imaging
CN102940510A (en) * 2012-08-31 2013-02-27 华南理工大学 Automatic focusing method for ultrasonic elastography
CN102920482A (en) * 2012-11-26 2013-02-13 重庆理工大学 Multi-frequency alternative-ejecting real-time ultrasonic elastography method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A strain-based ultrasound elastography using phase shift with prior estimats and meshfree shape function;Yuan Jinwei 等;《IEEE International Symposium on Biomedical Imaging: From Nano to Macro》;20111231;532-535 *
基于先验估计和相关相位的实时超声弹性成像;崔少国 等;《中国生物医学工程学报》;20110820;第30卷(第4期);636-640 *
基于相移的实时超声弹性成像方法;孙瑞超 等;《中国生物医学工程学报》;20130630;第32卷(第3期);339-347 *

Also Published As

Publication number Publication date
CN103735284A (en) 2014-04-23

Similar Documents

Publication Publication Date Title
Song et al. Two-dimensional shear-wave elastography on conventional ultrasound scanners with time-aligned sequential tracking (TAST) and comb-push ultrasound shear elastography (CUSE)
CN102592137B (en) Multi-modality image registration method and operation navigation method based on multi-modality image registration
WO2014174305A3 (en) A method and apparatus for determining the location of a medical instrument with respect to ultrasound imaging, and a medical instrument to facilitate such determination
CN102314712B (en) Pelvis modeling and pelvis parameter measurement method based on three-dimensional imaging technology
KR102448397B1 (en) Bone age calculation system and method
CN103750860A (en) Wireless three-dimensional ultrasound imaging method and device
CN103295234B (en) Based on the medical image segmentation system and method for deformation surface model
CN103040488B (en) System and method for real-time ultrasonic elastography displacement estimation
CN104306021B (en) The ultrasonoscopy three-dimensional rebuilding method that global registration optimizes
CN102908166B (en) Scanning device and method for three-dimensional ultrasound elasticity imaging
CN105796101B (en) Displacement of tissue measurement method and system based on the imaging of magnetic resonance acoustic radiation force
CN102860842A (en) Real-time quasi-static ultrasonic elastography method
CN103006263B (en) Position calibrating method of medical ultrasonic three-dimensional imaging based on linear scanning
CN102608219B (en) Device for expanding ultrasonic detection region and increasing detection precision and method
CN103735284B (en) Based on RF signal estimation method in the three-D ultrasonic elastogram of linear scanning
CN104622507A (en) Elasticity modulus measurement method and system
CN103815932A (en) Ultrasonic quasi-static elastic imaging method based on optical flow and strain
CN202636976U (en) Elastic imaging equipment of shearing wave interference method
Yuan et al. A strain-based ultrasound elastography using phase shift with prior estimates and meshfree shape function
CN103156604B (en) Biological tissue&#39;s electrical impedance imaging method of supersonic synergic
KR20200058441A (en) Method and apparatus for measuring semi-elastic modulus and elasticity of human tissue
CN103767733B (en) The method of estimation of RF volume data in freedom-arm, three-D ultrasonic elastograph imaging
CN202512242U (en) Apparatus expanding ultrasonic detection area and improving ultrasonic detection precision
Ingle et al. Ultrasonic tracking of shear waves using a particle filter
US20130116560A1 (en) Ultrasound temperature mapping system and method

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20151028

Termination date: 20211226