CN104188689B - Based on displacement of tissue evaluation method and the system of ultrasonic echo radiofrequency signal - Google Patents
Based on displacement of tissue evaluation method and the system of ultrasonic echo radiofrequency signal Download PDFInfo
- Publication number
- CN104188689B CN104188689B CN201410378475.4A CN201410378475A CN104188689B CN 104188689 B CN104188689 B CN 104188689B CN 201410378475 A CN201410378475 A CN 201410378475A CN 104188689 B CN104188689 B CN 104188689B
- Authority
- CN
- China
- Prior art keywords
- ultrasonic echo
- radiofrequency signal
- echo radiofrequency
- displacement
- frame
- 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.)
- Active
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The present invention relates to a kind of displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal and system. Described method, comprises the following steps: switch process, and the ultrasonic echo radiofrequency signal of whole frame is converted to analytic signal; Partiting step, carries out segmentation by the analytic signal of whole frame according to preset window length, obtains segmentation window data section; Sum step, to two adjacent frame signals, is undertaken corresponding segments window data section obtaining target plural number once without the sum computing offset; Phase calculation step, calculates the phase place of described target plural number; Relative displacement estimation steps, utilizes described phase calculation to obtain the relative displacement between two frame signals. The above-mentioned displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal, the ultrasonic echo radiofrequency signal of whole frame is once converted to analytic signal, greatly improve computing velocity, the relative displacement calculating two frame signals only can need to obtain through a sum computing, reduce calculated amount, improve counting yield, and complicacy is low.
Description
Technical field
The present invention relates to ultra sonic imaging process field, particularly relate to a kind of displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal and system.
Background technology
Acoustic radiation force pulse elastogram (AcousticRadiationForceImpulseImaging, ARFI) it is the acoustic radiation force utilizing the focused ultrasound beams in safe power range to produce, make biological tissue's regional area produce microdeformation and also form the shearing wave laterally propagated, miniature deformation is followed the trail of and is calculated the lateral propagation speed of shearing wave by recycling ultrasonic echo tracking technique, finally utilizes a kind of ultrasonic elastograph imaging technology of the Young's modulus that elasticity reconstruct algorithm organized by shear wave velocity quantitative estimation. This acoustic radiation force pulse elastogram technology is applied to the aspects such as breast cancer detection, the diagnosis of artery atherosclerotic patch, focused ultrasound therapy security monitoring.
Acoustic radiation force pulse elastogram is adopted to measure in displacement of tissue, there is the signal relative time-delay caused by displacement of tissue between two adjacent frame ultrasonic radio frequency echoed signals, usually can not just be divided exactly by the digital sample cycle of rf echo signal, if directly using time domain computing cross-correlation, the time lag corresponding to signal numeral sampling period integral multiple can only be obtained, and the time lag part being less than a digital sampling period will be ignored, displacement of tissue calculation result now will produce relatively big error. In order to solve the bigger problem of error, (1) Cespedes etc. proposed in nineteen ninety-five, after carrying out the time domain cross-correlation of signal, find the position of the maximum point of cross-correlation coefficient, maximum point and front and back two point (totally 3 points) thereof are carried out parabola interpolation, the position of its maximum point is picked up again within the scope of interpolation, the size of this position deviation starting position, it is the relative time-delay of signal, but the method needs to carry out fitting of parabola and interpolation, calculated amount is big and the Displacement Estimation precision lower than a sampling interval is not high; (2) propose in 1981 by Cabot, the real number numerary signal of the cross-correlation coefficient of gained is configured to the analytic signal of plural number, the phase place near the numerary signal maximum value of cross-correlation coefficient is adopted to cross zero point, corresponding to the maximum point of accurate cross-correlation coefficient continuous signal, the accurate relative time-delay of signal can be obtained by the size of the position of accurate maximum point deviation starting position, but the method needs to calculate three sum computings, calculated amount is big; (3) based on the Displacement Estimation algorithm of Doppler (doppler) velocity estimation, need original ultrasonic radio frequency echoed signal is carried out IQ decomposition, namely In-phase (real part of complex-envelope signal) and Quadraturephase (imaginary part of complex-envelope signal) two paths of signals is obtained, calculated amount is big, and complicacy height, counting yield are low.
Summary of the invention
Based on this, it is necessary to the counting yield for existing estimation displacement of tissue is low, the problem that complicacy is high, it is provided that a kind of counting yield height and the low displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal of complicacy.
In addition, a kind of counting yield height and the low displacement of tissue estimating system based on ultrasonic echo radiofrequency signal of complicacy are provided provide.
Based on a displacement of tissue evaluation method for ultrasonic echo radiofrequency signal, comprise the following steps:
Switch process, is converted to analytic signal by the ultrasonic echo radiofrequency signal of whole frame;
Partiting step, carries out segmentation by the analytic signal of whole frame according to preset window length, obtains segmentation window data section;
Sum step, to two adjacent frame signals, is undertaken corresponding segments window data section obtaining target plural number once without the sum computing offset;
Phase calculation step, calculates the phase place of described target plural number;
Relative displacement estimation steps, utilizes described phase calculation to obtain the relative displacement between two frame signals. Wherein in an embodiment, before described switch process, also comprise:
Piecemeal step, is divided into multiple data matrix by ultrasonic echo radiofrequency signal, and each data matrix comprises the ultrasonic echo radiofrequency signal of the whole frame of some frames, and each frame ultrasonic echo radiofrequency signal is corresponding to a ultrasonic scanning line of ultrasound probe axial direction due;
Described switch process comprises:
The ultrasonic echo radiofrequency signal of the whole frame of each data matrix is converted to analytic signal;
After described relative displacement estimation steps, also comprise:
Absolute displacement estimation steps, carries out the relative displacement of multiple data matrix adding up on time orientation respectively, obtains absolute displacement.
Wherein in an embodiment, described switch process comprises:
Fourier transformation and inverse transformation is adopted to obtain analytic signal the ultrasonic echo radiofrequency signal of whole frame.
Wherein in an embodiment, described switch process adopts the first kernel function realized by graphic process unit to process, the block number of described first kernel function is the line number of data matrix, described line number refers to the frame number of the ultrasonic echo radiofrequency signal of whole frame, the thread number of described first kernel function is the row number of data matrix, and described row number refers to the number of data points on ultrasound probe axial direction due sweep trace.
Wherein in an embodiment, described sum step, phase calculation step and relative displacement estimation steps adopt the 2nd kernel function realized by graphic process unit to calculate; The line number that the block number of described 2nd kernel function is data matrix subtracts 1, the thread number of each block is segmentation window number, described segmentation window number is by the difference of the row number of described data matrix and window width, divided by step-length, adding 1 again to obtain, described step-length refers to the number of data points that two adjacent segmentation windows comprise except the part overlapped each other.
Wherein in an embodiment, described relative displacement estimation steps also comprises: if the relative displacement between the segmentation window data section of two each correspondences of frame signal is less than 1/4th of ultrasonic echo radiofrequency signal wavelength, then without the need to carrying out extra computation; If the relative displacement between two frame signals is more than or equal to 1/4th of ultrasonic echo radiofrequency signal wavelength, then carry out time domain cross-correlation calculation and modified result calculating.
Based on a displacement of tissue estimating system for ultrasonic echo radiofrequency signal, comprising:
Conversion module, for being converted to analytic signal by the ultrasonic echo radiofrequency signal of whole frame;
Divide module, for the analytic signal of whole frame is carried out segmentation according to preset window length, obtain segmentation window data section;
Sum module, for two adjacent frame signals, being undertaken corresponding segments window data section obtaining target plural number once without the sum computing offset;
Phase calculation module, for calculating the phase place of described target plural number;
Relative displacement estimation block, for utilizing described phase calculation to obtain the relative displacement between two frame signals.
Wherein in an embodiment, described system also comprises piecemeal module and absolute displacement estimation block,
Described piecemeal module is for being divided into multiple data matrix by ultrasonic echo radiofrequency signal, and each data matrix comprises the ultrasonic echo radiofrequency signal of the whole frame of some frames, and each frame ultrasonic echo radiofrequency signal is corresponding to a ultrasonic scanning line of ultrasound probe axial direction due;
Described conversion module is also for being converted to analytic signal by the ultrasonic echo radiofrequency signal of the whole frame of each data matrix;
Described absolute displacement estimation block, for the relative displacement of multiple data matrix carries out adding up on time orientation respectively, obtains absolute displacement.
Wherein in an embodiment, described conversion module is also for adopting fourier transformation and inverse transformation to obtain analytic signal the ultrasonic echo radiofrequency signal of whole frame.
Wherein in an embodiment, described conversion module adopts the first kernel function realized by graphic process unit to process, the block number of described first kernel function is the line number of data matrix, described line number refers to the frame number of the ultrasonic echo radiofrequency signal of whole frame, the thread number of described first kernel function is the row number of data matrix, and described row number refers to the number of data points on ultrasound probe axial direction due sweep trace.
Wherein in an embodiment, described sum module, phase calculation module and relative displacement estimation block adopt the 2nd kernel function realized by graphic process unit to calculate; The line number that the block number of described 2nd kernel function is data matrix subtracts 1, the thread number of each block is segmentation window number, described segmentation window number is by the difference of the row number of described data matrix and window width, divided by step-length, adding 1 again to obtain, described step-length refers to the number of data points that two adjacent segmentation windows comprise except the part overlapped each other.
Wherein in an embodiment, if described relative displacement estimation block is also less than 1/4th of ultrasonic echo radiofrequency signal wavelength for the relative displacement between the segmentation window data section of two each correspondences of frame signal, then without the need to carrying out extra computation; And if the relative displacement between two frame signals is more than or equal to 1/4th of ultrasonic echo radiofrequency signal wavelength, then carry out time domain cross-correlation calculation and modified result calculating.
The above-mentioned displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal and system, the ultrasonic echo radiofrequency signal of whole frame being once converted to analytic signal, greatly improves computing velocity, the relative displacement calculating two frame signals only can need to obtain through a sum computing, do not need to carry out multiple sum computing, simplify calculation procedure, and do not need to carry out signal decomposition, simplify calculating, reduce calculated amount, substantially reduce computing time, it is to increase counting yield, and complicacy is low.
In addition, the ultrasonic echo radiofrequency signal of whole frame is once carried out fourier transformation and Fourier's inverse transformation, compared to segmentation rear hatch data section is carried out tens times to several hundred times fourier transformation and Fourier's inverse transformation, greatly improve computing velocity; When adopting graphic process unit to realize, it is possible to use multiple block number and thread number parallel processing, it is to increase processing efficiency.
Accompanying drawing explanation
Fig. 1 is the schema of displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal in an embodiment;
Fig. 2 is separately the schema of displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal in an embodiment;
Fig. 3 is the structure block diagram of displacement of tissue estimating system based on ultrasonic echo radiofrequency signal in an embodiment;
Fig. 4 is the structure block diagram of displacement of tissue estimating system based on ultrasonic echo radiofrequency signal in an embodiment.
Embodiment
In order to make the object of the present invention, technical scheme and advantage clearly understand, below in conjunction with drawings and Examples, the present invention is further elaborated. It is to be understood that specific embodiment described herein is only in order to explain the present invention, it is not intended to limit the present invention.
Fig. 1 is the schema of displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal in an embodiment. Based on the displacement of tissue evaluation method of ultrasonic echo radiofrequency signal, should comprise the following steps:
Step 102, switch process, is converted to analytic signal by the ultrasonic echo radiofrequency signal of whole frame.
In the present embodiment, adopt fourier transformation and inverse transformation to obtain analytic signal the ultrasonic echo radiofrequency signal of whole frame, be converted to complex signal by real number signal. Concrete switching process is: first transfer time-domain signal to frequency domain signal by fourier transformation, and then real part by frequency domain signal is multiplied by 2, and imaginary part sets to 0, and then gained signal is carried out Fourier's inverse transformation, obtains the analytic signal of corresponding original time-domain signal. The ultrasonic echo radiofrequency signal of whole frame is once carried out fourier transformation and Fourier's inverse transformation, compared to segmentation rear hatch data section is carried out tens times to several hundred times fourier transformation and Fourier's inverse transformation, greatly improves computing velocity.
Step 104, partiting step, carries out segmentation by the analytic signal of whole frame according to preset window length, obtains segmentation window data section.
Concrete, preset window length can set as required.
Step 106, sum step, to two adjacent frame signals, is undertaken corresponding segments window data section obtaining target plural number once without the sum computing offset.
Concrete, first the data section of correspondence is carried out complex multiplication (plural number of the first frame data point is multiplied by the plural conjugation of the 2nd frame corresponding data point), the more all products that will obtain are added, obtain a plural form and a+jb, i.e. target plural number. Wherein, a, b are real number, and j is imaginary part identifier. N-th segmentation window data section of such as the first frame comprises data point x1+jy1,��xn+jyn, the n-th segmentation window data section of the correspondence of the 2nd frame comprises data point z1+jc1,��zn+jcn, then a+jb=(x1+jy1)��(z1-jc1)+��+(xn+jyn)��(zn-jcn)��
Step 108, phase calculation step, calculates the phase place of this target plural number.
Concrete, the phase place of target plural numberBy the positive minus symbol of a, b, can by phase placeRange expansion to [-�� ,-��]. The scope of the relative displacement corresponding to this scope is [-L, L], and L represents 1/4th of ultrasonic echo radiofrequency signal wavelength.
Step 110, relative displacement estimation steps, utilizes this phase calculation to obtain the relative displacement between two frame signals.
Concrete, the calculation formula of relative displacement is:
Wherein, d is relative displacement, and c is velocity of ultrasonic sound, and f is ultrasound emission frequency,For phase place.
The above-mentioned displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal, the ultrasonic echo radiofrequency signal of whole frame being once converted to analytic signal, greatly improves computing velocity, the relative displacement calculating two frame signals only can need to obtain through a sum computing, do not need to carry out multiple sum computing, simplify calculation procedure, and do not need to carry out signal decomposition, simplify calculating, reduce calculated amount, substantially reduce computing time, it is to increase counting yield, and complicacy is low.
In an embodiment, above-mentioned relative displacement estimation steps also comprises: if the relative displacement between the segmentation window data section of two each correspondences of frame signal is less than 1/4th of ultrasonic echo radiofrequency signal wavelength, then without the need to carrying out extra computation; If the relative displacement between two frame signals is more than or equal to 1/4th of ultrasonic echo radiofrequency signal wavelength, then also need to carry out extra time domain cross-correlation calculation and modified result calculating. We represent 1/4th of ultrasonic echo radiofrequency signal wavelength with L, K represents that the length of a sampling interval (can be calculated by formula 0.5*c/fs, c is velocity of ultrasonic sound, fs is the sample frequency of ultrasonic echo radiofrequency signal), then relative displacement d can represent the multiple for L, i.e. d=(N+m) * L, and wherein N is unknown integral part, m is known mantissa, and m is namely corresponding to above-mentionedMeanwhile, can also representing to be the multiple of K by the known relative displacement d of the result of cross-correlation calculation, M*K < d < (M+1) * K, namely the span of d is between M*K and (M+1) * K. By the span with co-relation and m, it is possible to easily calculate N, and the relative displacement d after being revised.
As shown in Figure 2, in another embodiment based on the schema of the displacement of tissue evaluation method of ultrasonic echo radiofrequency signal. Wherein in an embodiment, before this switch process, also comprise:
Step 202, piecemeal step, is divided into multiple data matrix by ultrasonic echo radiofrequency signal, and each data matrix comprises the ultrasonic echo radiofrequency signal of the whole frame of some frames, and each frame ultrasonic echo radiofrequency signal is corresponding to a ultrasonic scanning line of ultrasound probe axial direction due.
Concrete, ultrasonic echo radiofrequency signal is divided into N number of data matrix by side direction collection position, each data matrix size is row �� col, represents and comprises row frame data (time orientation), and every frame has col element (depth direction).
Step 204, switch process, is converted to analytic signal by the ultrasonic echo radiofrequency signal of the whole frame of each data matrix.
Switch process adopts by graphic process unit (GraphicsProcessUnit, GPU) the first kernel function realized processes, the block number of the first kernel function is the line number of data matrix, the i.e. frame number of the ultrasonic echo radiofrequency signal of whole frame, the thread number of the first kernel function is the row number of data matrix, i.e. number of data points on ultrasound probe axial direction due sweep trace.
Concrete, adopting GPU optimization design, fourier transformation and Fourier's inverse transformation all use the CUFFT storehouse function cufftExecC2C of CUDA (CPU+GPU). Process for frequency domain signal adopts kernel function kernel_1 to realize. CPU is CentralProcessingUnit, central processing unit.
First kernel function call order: 1) cufftExecC2C (just converting); 2) kernel_1; 3) cufftExecC2C (inverse transformation).
The blocks quantity that in first kernel function, each step calculating adopts is identical with threads quantity, is blocks=dim3 (input-> rows); Threads=dim3 (input-> cols), input-> rows is the line number of data matrix, and input-> cols is the row number of data matrix.
Step 206, partiting step, carries out segmentation by the analytic signal of whole frame according to preset window length, obtains segmentation window data section.
Concrete, in the depth direction the analytic signal of whole frame is carried out segmentation according to preset window length. Wherein, preset window length can set as required.
Step 208, sum step, to two adjacent frame signals, is undertaken corresponding segments window data section obtaining target plural number once without the sum computing offset.
Concrete, to two frame signals adjacent on time orientation, undertaken segmentation window data section corresponding on depth direction obtaining target plural number once without the sum computing offset.
Step 210, phase calculation step, calculates the phase place of this target plural number.
Concrete, the phase place of target plural numberAnd can incite somebody to action according to the symbol of a, bSpan expand to [-��, ��]. The scope of the relative displacement corresponding to this scope is [-L, L], and L represents 1/4th of ultrasonic echo radiofrequency signal wavelength.
Step 212, relative displacement estimation steps, utilizes this phase calculation to obtain the relative displacement between two frame signals.
Concrete, the calculation formula of relative displacement is:
Wherein, d is relative displacement, and c is velocity of ultrasonic sound, and f is ultrasound emission frequency,For phase place.
Step 214, absolute displacement estimation steps, carries out the relative displacement of multiple data matrix adding up on time orientation respectively, obtains absolute displacement.
Concrete, each segmentation window data section in the depth direction can be obtained the relative displacement between adjacent two frame data according to step 204 to 212. So can obtaining relative displacement matrix and be of a size of (row-1) �� Win_Num, wherein, Win_Num is segmentation window number.
In an embodiment, sum step, phase calculation step and relative displacement estimation steps adopt the 2nd kernel function realized by graphic process unit (GraphicsProcessUnit, GPU) to calculate; The line number that the block number of the 2nd kernel function is data matrix subtracts 1, the thread number of each block is segmentation window number, this segmentation window number is by the difference of the row number of this data matrix and window width, divided by step-length, adding 1 again to obtain, this step-length refers to the number of data points that two adjacent segmentation windows comprise except the part overlapped each other.
Concrete, the 2nd kernel function adopts ComputeDisp. Segmentation window number WinNum=(InputWidth-WindowHW)/Step+1, wherein, WindowHW is window width, the row number of step-length to be Step, InputWidth be data matrix.
The thread (i.e. Thread Count) that 2nd kernel function ComputeDisp designs each block equals segmentation window number WinNum, the calculating of whole matrix is realized by (rows-1) individual block, and each block calculates the relative displacement between one group of two adjacent frame signal. Therefore kernel function Thread Count is as follows:
2nd kernel function Thread Count ComputeDisp<<<blocks,threads>>>(), wherein, blocks=dim3 (rows-1); Threads=dim3 (WinNum). Wherein,<<<��>>>it is CUDA grammer symbol, represents that GPU performs the Thread Count of kernel needs startup, it may also be useful to mode:<<<��ʾ�����block��������ʾ1��block������߳���>>>. Input is a structure type ,->represent the element getting structure. Adopt multiple block number and thread number parallel processing, it is to increase processing efficiency.
As shown in Figure 3, it is the structure block diagram of displacement of tissue estimating system based on ultrasonic echo radiofrequency signal in an embodiment. Based on the displacement of tissue estimating system of ultrasonic echo radiofrequency signal, conversion module 310 should be comprised, divides module 320, sum module 330, phase calculation module 340 and relative displacement estimation block 350. Wherein:
Conversion module 310, for being converted to analytic signal by the ultrasonic echo radiofrequency signal of whole frame.
Concrete, this conversion module 310 is also for adopting fourier transformation and inverse transformation to obtain analytic signal the ultrasonic echo radiofrequency signal of whole frame.
Divide module 320, for the analytic signal of whole frame is carried out segmentation according to preset window length, obtain segmentation window data section.
Concrete, preset window length can set as required.
Sum module 330, for two adjacent frame signals, being undertaken corresponding segments window data section obtaining target plural number once without the sum computing offset.
Concrete, first the data section of correspondence is carried out complex multiplication, the more all products that will obtain are added, obtain a plural form and a+jb, i.e. target plural number.
Phase calculation module 340, for calculating the phase place of this target plural number.
Concrete, the phase place of target plural numberAnd can incite somebody to action according to the symbol of a, bSpan expand to [-��, ��]. The scope of the relative displacement corresponding to this scope is [-L, L], and L represents 1/4th of ultrasonic echo radiofrequency signal wavelength.
Relative displacement estimation block 350, for utilizing this phase calculation to obtain the relative displacement between two frame signals.
Concrete, the calculation formula of relative displacement is:
Wherein, d is relative displacement, and c is velocity of ultrasonic sound, and f is ultrasound emission frequency,For phase place.
The above-mentioned displacement of tissue estimating system based on ultrasonic echo radiofrequency signal, the ultrasonic echo radiofrequency signal of whole frame being once converted to analytic signal, greatly improves computing velocity, the relative displacement calculating two frame signals only can need to obtain through a sum computing, do not need to carry out multiple sum computing, simplify calculation procedure, and do not need to carry out signal decomposition, simplify calculating, reduce calculated amount, substantially reduce computing time, it is to increase counting yield, and complicacy is low.
As shown in Figure 4, it is the structure block diagram of displacement of tissue estimating system based on ultrasonic echo radiofrequency signal in an embodiment. Should based on the displacement of tissue estimating system of ultrasonic echo radiofrequency signal, except comprising conversion module 310, divide module 320, sum module 330, phase calculation module 340 and relative displacement estimation block 350, also comprise piecemeal module 360 and absolute displacement estimation block 370.
This piecemeal module 360 is for being divided into multiple data matrix by ultrasonic echo radiofrequency signal. Concrete, ultrasonic echo radiofrequency signal is divided into N number of data matrix by side direction collection position, each data matrix size is row �� col, represents and comprises row frame data (time orientation), and every frame has col element (depth direction).
This conversion module 310 is also for being converted to analytic signal by the ultrasonic echo radiofrequency signal of the whole frame of each data matrix. This conversion module 310 adopts the first kernel function realized by graphic process unit to process, the block number of this first kernel function is the line number of data matrix, the i.e. frame number of the ultrasonic echo radiofrequency signal of whole frame, the thread number of this first kernel function is the row number of data matrix, i.e. number of data points on ultrasound probe axial direction due sweep trace.
This division module 320 also in the depth direction the analytic signal of whole frame being carried out segmentation according to preset window length, obtains segmentation window data section.
This sum module 330 is also for two frame signals adjacent on time orientation, being undertaken segmentation window data section corresponding on depth direction obtaining target plural number once without the sum computing offset.
Relative displacement estimation block 350 also obtains the relative displacement between adjacent two frame data for each segmentation window data section in the depth direction.
This absolute displacement estimation block 370, for the relative displacement of multiple data matrix carries out adding up on time orientation respectively, obtains absolute displacement. Concrete, absolute displacement estimation block 370 carries out adding up on time orientation respectively for the relative displacement of data matrix just multiple on time orientation, obtains absolute displacement.
In an embodiment, this sum module 330, phase calculation module 340 and relative displacement estimation block 350 adopt the 2nd kernel function realized by graphic process unit to calculate; The line number that the block number of the 2nd kernel function is data matrix subtracts 1, the thread number of each block is segmentation window number, this segmentation window number is by the difference of the row number of this data matrix and window width, divided by step-length, adding 1 again to obtain, this step-length refers to the number of data points that two adjacent segmentation windows comprise except the part overlapped each other.
The above-mentioned displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal and system, can be applicable to any need to utilize ultrasonic echo radiofrequency signal to carry out displacement of tissue estimation equipment or method in, such as the two-dimensional ultrasound elastogram based on acoustic radiation force, quasistatic ultrasonic elastograph imaging (quasi-staticelastography), other ultrasonic elastograph imaging methods such as transient state ultrasonic elastograph imaging (transientelastography), or ultrasonic temperature imaging (temperatureimaging), the methods such as ultrasonic thermal strain imaging (thermalstrainimaging).
Adopting Matlab simulation, calculate the time needed for one group of ARFI data (4 matrixes), the algorithm (comprising cross-correlation calculation) proposed with Cabot needs about 7s, does not carry out cross-correlation calculation, it is necessary to approximately 2s; With the Displacement Estimation algorithm based on Doppler velocity estimation, if decomposed with software simulating IQ, it is necessary to approximately 11s, if decomposed with hardware implementing IQ, then needs about 0.5s. Adopt the algorithm of the present invention, it is also desirable to approximately 0.5s, but relative to the Displacement Estimation algorithm based on Doppler velocity estimation, it is not necessary to Hardware I Q decomposition circuit, reduces cost and the complicacy of system. The displacement of tissue obtained by algorithms of different is all 10-6To 10-7The order of magnitude, and the difference between algorithms of different gained displacement of tissue is 10-8To 10-9The order of magnitude, therefore the accuracy that displacement of tissue is estimated is not caused remarkably influenced by this algorithm.
In addition, in an embodiment, if relative displacement estimation block 350 is also less than 1/4th of ultrasonic echo radiofrequency signal wavelength for the relative displacement between the segmentation window data section of two each correspondences of frame signal, then without the need to carrying out extra computation; If the relative displacement between two frame signals is more than or equal to 1/4th of ultrasonic echo radiofrequency signal wavelength, then also need after calculating the phase place of this target plural number, increase extra time domain cross-correlation calculation and modified result calculating. 1/4th of ultrasonic echo radiofrequency signal wavelength are represented with L, K represents that the length of a sampling interval (can be calculated by formula 0.5*c/fs, c is velocity of ultrasonic sound, fs is the sample frequency of ultrasonic echo radiofrequency signal), then relative displacement d can represent the multiple for L, i.e. d=(N+m) * L, and wherein N is unknown integral part, m is known mantissa, and m is namely corresponding to above-mentionedMeanwhile, can also representing to be the multiple of K by the known relative displacement d of the result of time domain cross-correlation calculation, M*K < d < (M+1) * K, namely the span of d is between M*K and (M+1) * K. By the span with co-relation and m, it is possible to easily calculate N, and the relative displacement d after being revised.
Described time domain cross-correlation calculation comprises normalization method correlation time-delay estimate method, normalization method covariance Time Delay Estimation Method, non-normalized correlation time-delay estimate method, mixed signs correlation time-delay estimate method and Fast Fourier Transform (FFT)/inverse transformation and calculates cross-correlation coefficient method etc. Wherein,
Normalization method correlation time-delay estimate method:
Normalization method covariance Time Delay Estimation Method:
Wherein:
Non-normalized correlation time-delay estimate method:
Mixed signs correlation time-delay estimate method:
Wherein:
Fast Fourier Transform (FFT)/inverse transformation calculates cross-correlation coefficient: first the time-domain signal that two to be carried out computing cross-correlation is carried out fourier transformation, transformed to frequency domain, calculate the cross-power spectrum of these two frequency domain signals again, cross-power spectrum signal is carried out Fourier's inverse transformation and switches back to time domain, what obtain is exactly the time domain cross-correlation coefficient signal of signal, and formula summary is as follows:
Sr=FFT (sr)
Sd=FFT (sd)(8)
Srd=Sr��Sd *
Rrd=IFFT (Srd)
Wherein, sr��sdFor time-domain signal, Sr, SdFor carrying out the corresponding frequency domain signal after fast Flourier conversion, SrdBeing both cross-power spectrums, * represents the conjugation of complex signal, RrdBeing the cross-correlation coefficient of time-domain signal, FFT is fourier transformation, and IFFT is Fourier's inverse transformation.
In addition, when the ultrasonic echo radiofrequency signal of whole frame is converted to analytic signal, it is also possible to adopt Hilbert transform to obtain. Hilbert transform formula is:
In formula (10), r^ (m, t) is the Hilbert transform of signal r (m, t), and m is the coordinate of sampling point, and t is the time, and * is convolution algorithm.
Then analytic signal corresponds to:
One of ordinary skill in the art will appreciate that all or part of flow process realizing in above-described embodiment method, it is can be completed by the hardware that computer program carrys out instruction relevant, described program can be stored in a computer read/write memory medium, this program, when performing, can comprise the flow process of the embodiment such as above-mentioned each side method. Wherein, described storage media can be magnetic disc, CD, read-only storage memory body (Read-OnlyMemory, ROM) or random storage and remembers body (RandomAccessMemory, RAM) etc.
The above embodiment only have expressed several enforcement modes of the present invention, and it describes comparatively concrete and detailed, but therefore can not be interpreted as the restriction to patent scope of the present invention. , it is also possible to make some distortion and improvement, it should be appreciated that for the person of ordinary skill of the art, without departing from the inventive concept of the premise these all belong to protection scope of the present invention. Therefore, the protection domain of patent of the present invention should be as the criterion with claims.
Claims (10)
1., based on a displacement of tissue evaluation method for ultrasonic echo radiofrequency signal, comprise the following steps:
Piecemeal step, is divided into multiple data matrix by ultrasonic echo radiofrequency signal, and each data matrix comprises the ultrasonic echo radiofrequency signal of the whole frame of some frames, and each frame ultrasonic echo radiofrequency signal is corresponding to a ultrasonic scanning line of ultrasound probe axial direction due;
Switch process, is converted to analytic signal by the ultrasonic echo radiofrequency signal of the whole frame of each data matrix;
Partiting step, carries out segmentation by the analytic signal of whole frame according to preset window length, obtains segmentation window data section;
Sum step, to two adjacent frame signals, is undertaken corresponding segments window data section obtaining target plural number once without the sum computing offset;
Phase calculation step, calculates the phase place of described target plural number;
Relative displacement estimation steps, utilizes described phase calculation to obtain the relative displacement between two frame signals;
Absolute displacement estimation steps, carries out the relative displacement of multiple data matrix adding up on time orientation respectively, obtains absolute displacement.
2. the displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal according to claim 1, it is characterised in that, described switch process comprises:
Fourier transformation and inverse transformation is adopted to obtain analytic signal the ultrasonic echo radiofrequency signal of whole frame.
3. the displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal according to claim 1, it is characterized in that, described switch process adopts the first kernel function realized by graphic process unit to process, the block number of described first kernel function is the line number of data matrix, described line number refers to the frame number of the ultrasonic echo radiofrequency signal of whole frame, the thread number of described first kernel function is the row number of data matrix, and described row number refers to the number of data points on ultrasound probe axial direction due sweep trace.
4. the displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal according to claim 1, it is characterised in that, described sum step, phase calculation step and relative displacement estimation steps adopt the 2nd kernel function realized by graphic process unit to calculate; The line number that the block number of described 2nd kernel function is data matrix subtracts 1, the thread number of each block is segmentation window number, described segmentation window number is by the difference of the row number of described data matrix and window width, divided by step-length, adding 1 again to obtain, described step-length refers to the number of data points that two adjacent segmentation windows comprise except the part overlapped each other.
5. the displacement of tissue evaluation method based on ultrasonic echo radiofrequency signal according to claim 1, it is characterized in that, described relative displacement estimation steps also comprises: if the relative displacement between the segmentation window data section of two each correspondences of frame signal is less than 1/4th of ultrasonic echo radiofrequency signal wavelength, then without the need to carrying out extra computation; If the relative displacement between two frame signals is more than or equal to 1/4th of ultrasonic echo radiofrequency signal wavelength, then carry out time domain cross-correlation calculation and modified result calculating.
6. the displacement of tissue estimating system based on ultrasonic echo radiofrequency signal, it is characterised in that, comprising:
Conversion module, for being converted to analytic signal by the ultrasonic echo radiofrequency signal of whole frame;
Divide module, for the analytic signal of whole frame is carried out segmentation according to preset window length, obtain segmentation window data section;
Sum module, for two adjacent frame signals, being undertaken corresponding segments window data section obtaining target plural number once without the sum computing offset;
Phase calculation module, for calculating the phase place of described target plural number;
Relative displacement estimation block, for utilizing described phase calculation to obtain the relative displacement between two frame signals;
Described system also comprises piecemeal module and absolute displacement estimation block,
Described piecemeal module is for being divided into multiple data matrix by ultrasonic echo radiofrequency signal, and each data matrix comprises the ultrasonic echo radiofrequency signal of the whole frame of some frames, and each frame ultrasonic echo radiofrequency signal is corresponding to a ultrasonic scanning line of ultrasound probe axial direction due;
Described conversion module is also for being converted to analytic signal by the ultrasonic echo radiofrequency signal of the whole frame of each data matrix;
Described absolute displacement estimation block, for the relative displacement of multiple data matrix carries out adding up on time orientation respectively, obtains absolute displacement.
7. the displacement of tissue estimating system based on ultrasonic echo radiofrequency signal according to claim 6, it is characterised in that, described conversion module is also for adopting fourier transformation and inverse transformation to obtain analytic signal the ultrasonic echo radiofrequency signal of whole frame.
8. the displacement of tissue estimating system based on ultrasonic echo radiofrequency signal according to claim 6, it is characterized in that, described conversion module adopts the first kernel function realized by graphic process unit to process, the block number of described first kernel function be data matrix line number described in line number refer to the frame number of ultrasonic echo radiofrequency signal of whole frame, the thread number of described first kernel function is the row number of data matrix, and described row number refers to the number of data points on ultrasound probe axial direction due sweep trace.
9. the displacement of tissue estimating system based on ultrasonic echo radiofrequency signal according to claim 6, it is characterised in that, described sum module, phase calculation module and relative displacement estimation block adopt the 2nd kernel function realized by graphic process unit to calculate; The line number that the block number of described 2nd kernel function is data matrix subtracts 1, the thread number of each block is segmentation window number, described segmentation window number is by the difference of the row number of described data matrix and window width, divided by step-length, adding 1 again to obtain, described step-length refers to the number of data points that two adjacent segmentation windows comprise except the part overlapped each other.
10. the displacement of tissue estimating system based on ultrasonic echo radiofrequency signal according to claim 6, it is characterized in that, if described relative displacement estimation block is also less than 1/4th of ultrasonic echo radiofrequency signal wavelength for the relative displacement between the segmentation window data section of two each correspondences of frame signal, then without the need to carrying out extra computation; And if the relative displacement between two frame signals is more than or equal to 1/4th of ultrasonic echo radiofrequency signal wavelength, then carry out time domain cross-correlation calculation and modified result calculating.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410378475.4A CN104188689B (en) | 2013-10-25 | 2014-08-01 | Based on displacement of tissue evaluation method and the system of ultrasonic echo radiofrequency signal |
Applications Claiming Priority (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310514494.0A CN103519848A (en) | 2013-10-25 | 2013-10-25 | Tissue displacement estimation method and system based on ultrasonic echo radio frequency signals |
| CN2013105144940 | 2013-10-25 | ||
| CN201310514494.0 | 2013-10-25 | ||
| CN201410378475.4A CN104188689B (en) | 2013-10-25 | 2014-08-01 | Based on displacement of tissue evaluation method and the system of ultrasonic echo radiofrequency signal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104188689A CN104188689A (en) | 2014-12-10 |
| CN104188689B true CN104188689B (en) | 2016-06-01 |
Family
ID=49922450
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310514494.0A Pending CN103519848A (en) | 2013-10-25 | 2013-10-25 | Tissue displacement estimation method and system based on ultrasonic echo radio frequency signals |
| CN201410378475.4A Active CN104188689B (en) | 2013-10-25 | 2014-08-01 | Based on displacement of tissue evaluation method and the system of ultrasonic echo radiofrequency signal |
Family Applications Before (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310514494.0A Pending CN103519848A (en) | 2013-10-25 | 2013-10-25 | Tissue displacement estimation method and system based on ultrasonic echo radio frequency signals |
Country Status (1)
| Country | Link |
|---|---|
| CN (2) | CN103519848A (en) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110507360B (en) * | 2014-08-28 | 2022-06-03 | 深圳迈瑞生物医疗电子股份有限公司 | Shear wave imaging method and system |
| CN104367346B (en) * | 2014-11-18 | 2016-07-06 | 声泰特（成都）科技有限公司 | Based on acoustic radiation force echo displacement detecting method and formation method |
| CN104306026B (en) * | 2014-11-18 | 2016-09-28 | 声泰特（成都）科技有限公司 | Based on acoustic radiation force echo displacement detection system and imaging system |
| CN105748100B (en) * | 2014-12-19 | 2019-04-16 | 深圳开立生物医疗科技股份有限公司 | Quasi-static ultrasonic elastograph imaging displacement calculates method and apparatus |
| CN107530051B (en) * | 2015-03-02 | 2021-03-16 | B-K医疗公司 | Non-invasive estimation of intravascular pressure changes using vector velocity Ultrasound (US) |
| US10582911B2 (en) * | 2015-08-11 | 2020-03-10 | Siemens Medical Solutions Usa, Inc. | Adaptive motion estimation in acoustic radiation force imaging |
| US20170135675A1 (en) * | 2015-11-12 | 2017-05-18 | Vanderbilt University | Adaptive clutter demodulation for ultrasound imaging |
| CN106037814B (en) * | 2016-05-13 | 2019-04-02 | 江西师范大学 | Ultrasonic elastography method based on fractal interpolation |
| CN106501775A (en) * | 2016-10-10 | 2017-03-15 | 南京威卡尔软件有限公司 | Continuous quick voice sound localization method for embedded platform |
| CN108553763B (en) * | 2018-01-19 | 2020-02-18 | 北京工业大学 | Microwave thermotherapy monitoring method based on ultrasonic echo decorrelation imaging technology |
| US11452503B2 (en) * | 2018-05-18 | 2022-09-27 | Siemens Medical Solutions Usa, Inc. | Shear wave imaging based on ultrasound with increased pulse repetition frequency |
| CN111399051B (en) * | 2019-01-02 | 2021-10-08 | 无锡海斯凯尔医学技术有限公司 | Data processing method, apparatus and medium |
| CN109615677B (en) * | 2019-02-13 | 2023-05-12 | 南京广慈医疗科技有限公司 | Method for calculating thermal strain distribution based on low sampling rate B ultrasonic image |
| CN110251159A (en) * | 2019-05-31 | 2019-09-20 | 深圳先进技术研究院 | A kind of tissue testing method, device and storage medium |
| CN110811689B (en) * | 2019-10-31 | 2020-11-27 | 汕头市超声仪器研究所股份有限公司 | First-order estimation method for shear wave motion speed |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1678243A (en) * | 2002-07-31 | 2005-10-05 | 株式会社日立医药 | Ultrasonic diagnosis system and distortion distribution display method |
| CN102423264A (en) * | 2011-09-01 | 2012-04-25 | 中国科学院深圳先进技术研究院 | Image-based biological tissue elasticity measuring method and device |
-
2013
- 2013-10-25 CN CN201310514494.0A patent/CN103519848A/en active Pending
-
2014
- 2014-08-01 CN CN201410378475.4A patent/CN104188689B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1678243A (en) * | 2002-07-31 | 2005-10-05 | 株式会社日立医药 | Ultrasonic diagnosis system and distortion distribution display method |
| CN102423264A (en) * | 2011-09-01 | 2012-04-25 | 中国科学院深圳先进技术研究院 | Image-based biological tissue elasticity measuring method and device |
Non-Patent Citations (1)
| Title |
|---|
| 声辐射力脉冲弹性成像技术及其算法研究;沈洋 等;《声学技术》;20111231;第30卷(第6期);149-150 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104188689A (en) | 2014-12-10 |
| CN103519848A (en) | 2014-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104188689B (en) | Based on displacement of tissue evaluation method and the system of ultrasonic echo radiofrequency signal | |
| CN102551791B (en) | A kind of ultrasonic imaging method and device | |
| CN103235041B (en) | Based on the Cavitation inciption threshold value distribution method for reconstructing of ultrasonic active cavitation imaging | |
| Smith et al. | Multistatic micro-Doppler radar signatures of personnel targets | |
| CN103439692B (en) | STAP method based on wide symmetrical characteristic of covariance matrix | |
| CN105411577B (en) | Fetal ECG signal separating method and system | |
| CN103116162B (en) | High-resolution sonar location method based on sparsity of objective space | |
| CN103431874B (en) | Method and system for estimating acoustic radiation force pulse imaging | |
| CN103300890A (en) | System and method for measuring tissue mechanical characteristics | |
| CN103675759A (en) | Modified FRFT (fractional Fourier transform) maneuvering weak target detection method | |
| CN104360251B (en) | A kind of ultrasonic signal delay time estimation method of partial discharge of transformer | |
| CN103519847A (en) | Doppler blood flow velocity estimation method and system based on ultrasonic echo radio frequency signals | |
| US20160018364A1 (en) | Methods, systems and computer program products for estimating shear wave speed using statistical inference | |
| CN110368031B (en) | Elasticity measurement detection method and system | |
| CN103245976A (en) | Human body target and surrounding structure compatibility detecting method based on UWB (Ultra Wideband) bio-radar | |
| CN104867163A (en) | Marginal distribution passing measurement-driven target tracking method and tracking system thereof | |
| CN102860842A (en) | Real-time quasi-static ultrasonic elastography method | |
| CN104138277B (en) | Based on the tissue microstructure detection method of empirical mode decomposition in quantitative ultrasound system | |
| CN103954962B (en) | A kind of ISAR imaging pulse algorithm for estimating based on compressed sensing | |
| EP2979645A1 (en) | Ultrasonic diagnostic device and ultrasonic diagnostic method | |
| CN110507361B (en) | Shear wave imaging method and system | |
| CN106373166A (en) | Multi-shape plane shear wave composite imaging method based on Mach cone effects | |
| CN103792527A (en) | Method for applying M sequence to phase encoding system imaging radar pulse compression | |
| CN104622507A (en) | Elasticity modulus measurement method and system | |
| CN103654847B (en) | Subject information acquisition device, display methods |
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 |