CN1586409A - Biological tissue displacement evaluating method using two kinds of size - Google Patents

Abstract

The present invention belongs to the field of supersonic elastic imaging technology. The two size biological tissue displacement estimating method includes: taking the data d1 of the size Ta in the first scanning line in the 2D RF signal before tissue compression, finding the cross correlation function of the data d1 and the scanning data after tissue compression, taking data e1 in the size Tb at the center of data section d1, finding the cross correlation function of the data e1 and the scanning line data to obtain the position t''1 corresponding to maximum data, and finally finding out the fine regulated displacement t1 of the small data section d1 to obtain the tissue displacement evaluation corresponding to the scanning lines. In case of great tissue compression amount, the present invention will combine the displacement estimation methods with great size and small size and thus has raised tissue displacement estimation precision.

Description

Adopt the biological tissue displacement method of estimation of two kinds of yardsticks

Technical field

The invention belongs to the ultrasonic elastograph imaging technical field, particularly the biological tissue displacement method of estimation.

Background technology

The variation of biological tissue elasticity modulus is relevant with its pathological phenomenon usually.For example, virulent pathological lesion, for example breast inocarcinoma, carcinoma of prostate, thyroid carcinoma and hepatic metastases etc. are usually expressed as hard brief summary.The breast inocarcinoma is the most common form of breast carcinoma, accounts for 3/4ths of breast carcinoma sum greatly, shows as fine and close lump owing to its density of matrix increases.The breast carcinoma of other types such as intraductal carcinoma and papilloma then show as softish tissue, and benign fibrocystic disease of breast also seldom shows as lump.

The elastic modelling quantity information of biological tissue has important value for the diagnostic procedure of disease.Yet the traditional medicine image mode that comprises x-ray imaging, ultra sonic imaging, computer tomography (CT) and nuclear magnetic resonance (MRI) etc. all can not directly provide the information about the basic mechanical attribute of this tissue of elastic modelling quantity.1991, J.Ophir proposed the method for ultrasonic elastograph imaging (ultrasound elastography), the elastic modelling quantity of tissue is distributed quantitatively estimate, imaging.At present, ultrasonic elastic modulus has become one of medical ultrasound image research focus, is widely used in the detection and the assessment of the infringement (lesion) that breast, prostate, atherosclerotic plaque, myocardium kinetics and high intensity focused ultrasound and radio-frequency (RF) ablation cause.

The ultimate principle of ultrasonic elastograph imaging is: ultrasonic probe is embedded in the extruding flat board, along vertical compress tissue of probe, gathers the forward and backward radiofrequency signal of tissue compression respectively; Tissue will produce a strain along compression direction in the tissue when being compressed, if organization internal elastic modelling quantity skewness, in-house stress distribution is difference to some extent also; The zone that elastic modelling quantity is bigger, the strain ratio that causes is less; Otherwise, the zone that elastic modelling quantity is less, corresponding strain ratio is bigger.Estimate the displacement of organization internal diverse location by certain methods, thereby calculate the stress distribution situation of organization internal, be used for the elastic modelling quantity of intermediate description organization internal to distribute, thereby describe physiology, the pathological state of tissue.

For the two-dimensional ultrasound elastogram, the general Type B ultrasonic probe that adopts linear array, gather the radiofrequency signal of forward and backward each the bar scanning line of probe of tissue compression, carry out above-described Displacement Estimation respectively, thereby the uniaxial train that calculates the corresponding tissue of each bar scanning line distributes.At last the uniaxial train of all scanning line correspondences is distributed and press the scanning line order and form a two-dimentional stress distribution, represent, be used for the elastic modelling quantity distribution of intermediate description organization internal with the form of gray-scale map or pcolor.

General ultrasonic elastograph imaging method may further comprise the steps:

1. utilize commercial B-mode ultrasonic apparatus device (generally adopting linear array probe) to obtain biological tissue to be measured and (be generally tissue, also can be animal tissue, hereinafter to be referred as tissue) a digitized two-dimentional radiofrequency signal before the compression (can adopt analog radio-frequency signal output termination signal amplifier, connect high-speed data acquisition card again, obtain digitized two-dimentional radiofrequency signal; Also can on digitized B-mode ultrasonic apparatus device, directly obtain digitized two-dimentional radiofrequency signal);

2. the probe of hand-held this B-mode ultrasonic apparatus device or utilize motor or screw drives this probe, along probe vertically this tissue is applied a small extruding (it is 1% the order of magnitude that the decrement that tissue is applied generally is controlled at), obtain a digitized two-dimentional radiofrequency signal after the tissue compression;

3. the data of taking out article one scanning line respectively from the forward and backward two-dimentional radiofrequency signal of the tissue compression that obtains of step 1 and 2 are made as s ₁(n) and s ₂(n), n represents the data sequence number on these two scanning lines, 1≤n≤n _Max, the maximum n of n _MaxProbing depth, the spread speed of ultrasonic waves transmitted in tissue and the sample frequency decision of radiofrequency signal by this B-mode ultrasonic apparatus device;

4. from this scan-line data s ₁(n) get the data d that a bit of length is T in ₁, its data number is U, U=round (T * U ₁), wherein, the unit of T is mm, U ₁Represent the data number of the tissue correspondence of 1mm, by spread speed and the sample frequency of the radiofrequency signal decision of ultrasonic waves transmitted in tissue, round () representative rounds up rounds operation, these data d ₁Sequence number from n ₁To n ₁+ U-1, n ₁Can be at 1≤n ₁Select in the scope of≤U; At τ ₁To τ ₂Ask this little segment data and scan-line data s in the hunting zone of determining ₂(n) cross-correlation function R (τ), computing formula is as follows

$R\left(\mathrm{\τ}\right)=\frac{\underset{i=n_1}{\overset{n_1+U-1}{\mathrm{\Σ}}}{s}_1\left(i\right)s_2(i-\mathrm{\τ})}{\sqrt{\underset{i=n_1}{\overset{n_1+U-1}{\mathrm{\Σ}}}{s}_1^2\left(i\right)\·\underset{i=n_1}{\overset{n_1+U-1}{\mathrm{\Σ}}}{s}_2^2(i-\mathrm{\τ})}}({\mathrm{\τ}}_1\≤\mathrm{\τ}\≤{\mathrm{\τ}}_2)$

Wherein i is the cyclic variable of computational process invading the exterior registration according to sequence number, τ ₁Be 0, τ ₂Be the decrement that tissue is applied, represent with the number of sampled data; In order to improve the precision of Displacement Estimation, generally also need to carry out interpolation, as conventional interpolation methods such as parabolas to calculating cross-correlation function;

5. the position t of the maximum correspondence of the cross-correlation function R (τ) that obtains of determining step 4 ₁t ₁Be exactly data d ₁Displacement after tissue compression (is s ₁(n) sequence number in is from n ₁To n ₁The segment data d of+U-1 ₁After tissue compression, move to s ₂(n) sequence number in is from n ₁-t ₁To n ₁+ U-1-t ₁The position);

6. successively from scan-line data s ₁(n) getting a bit of length in is that T is that the data number is the data d of U ₂, d ₃..., d _N, the sequence number of every segment data V the sampled data that stagger successively, 1≤V≤U, V sampled data will exceed s up to staggering again ₁(n) scope, 4,5 identical methods obtain the displacement t of each segment data correspondence successively set by step ₂, t ₃..., t _N, wherein N is the sum of little segment data; Then displacement sequence t ₁, t ₂..., t _NBe article one scan-line data s ₁(n) Displacement Estimation of Dui Ying tissue;

7. utilize the method identical with step 3-6, obtain the 2nd, 3 successively ..., M bar scan-line data correspondence the Displacement Estimation of tissue, wherein M is the scanning line sum of expression probe, by the probe decision;

8. to article one scan-line data s ₁(n) the Displacement Estimation sequence t of Dui Ying tissue ₁, t ₂..., t _NAsk difference, obtain organizing article one scanning line s ₁(n) stress distribution of corresponding tissue, computing formula is as follows,

${\mathrm{\ϵ}}_1=\frac{t_2-t_1}{V},{\mathrm{\ϵ}}_2=\frac{t_3-t_2}{V},\·\·\·,{\mathrm{\ϵ}}_{N-1}=\frac{t_N-t_{N-1}}{V}$

Wherein, ε ₁, ε ₂..., ε _N-1Be respectively d ₁, d ₂..., d _N-1Corresponding organize strain;

9. utilize the method identical with step 8, obtain organizing the 2nd, 3 successively ..., M bar scan-line data correspondence the stress distribution of tissue;

10. the stress distribution of the M bar scan-line data correspondence that step 9 is obtained is synthesized a 2-D data according to the der group of scanning line, and is showed with the form of gray-scale map or pcolor, just obtains the two-dimentional diagram of strains of organizing.

In ultrasonic elastograph imaging, crucial problem is the Displacements Distribution of tissue is estimated, just the step 3-7 of above-described method.The value of cross-correlation function is big more, illustrate that the forward and backward little segment data of compression is identical well more, the maximum value position of cross-correlation function has been represented the position of the correspondence after compression of the little segment data before the compression, thereby can obtain the displacement of this little segment data, just the displacement of the tissue of this segment data correspondence.

In ultrasonic elastograph imaging, the little segment data in the compression front signal of selection is used for following the tracks of the displacement of corresponding segment tissue, is called the tracking wave band.Its length is called follows the tracks of wave band length, perhaps yardstick.It is extremely important in the Displacement Estimation of ultrasonic elastograph imaging to choose suitable yardstick.When the decrement that tissue is applied in the time of smaller (as less than 1%), yardstick is big more, and the information that comprises is many more, and the precision of Displacement Estimation is high more, and is also insensitive more to the interference of random noise.When the decrement that tissue is applied bigger (as 1-5%) time, yardstick is big more, because tissue is compressed, the forward and backward signal superposed part of tissue compression is more little, it is low more to utilize the cross-correlation function maximum value position to carry out the precision of Displacement Estimation, simultaneously, yardstick is big more, is not easy to be subjected to the interference of random noise again more; And when yardstick was less, the forward and backward signal superposed part of tissue compression was bigger, and the precision of Displacement Estimation is higher, yet the Displacement Estimation of this moment is responsive to the interference ratio of random noise again.

In order to increase the degree that the forward and backward signal of tissue compression overlaps, generally adopt approach based on linear interpolation, the signal that is ' compressed ' " stretching " (stretching) is become long together with original signal, thus the registration of waveform before and after the increase compression.In waveform stretching, drawing coefficient is relevant with decrement that tissue is applied or strain (representing with percentage ratio).When drawing coefficient near the tissue strain the time, the registration maximum of waveform, thereby cross-correlation function maximum maximum.But the strain of tissue is unknown before time delay is estimated, the selection of drawing coefficient can only be done roughly to judge according to the degree of depth of displacement that applies and tissue.And because the inhomogeneities of tissue, the strained distribution of organization internal is also inhomogeneous, and therefore specific drawing coefficient is in some part of tissue deficiency that may stretch, and it is excessive to stretch again in some part, and this introduces error conversely again.Therefore, the someone proposes a kind of the wave band of diverse location to be carried out the method for self adaptation drawing coefficient, makes the cross-correlation function maximum maximum of wave band before wave band and the compression after stretching, utilizes the drawing coefficient direct estimation to go out the strain of this position tissue then.But the tensile method of this self adaptation is longer computation time, may be restricted in actual applications.

From angle of practical application, wishing the decrement that tissue is applied, bigger good, can more can " see " details of the stress distribution difference that the organization internal elastic modulus difference causes like this clearly.The general repeatedly method of compression (multi-compression) that adopts is decomposed into little displacement sum with big displacement, to the tissue little displacement of exerting pressure, calculates corresponding stress distribution one by one.At last with the stack of several stress distribution, the stress distribution when obtaining big displacement.Yet this method has increased data acquisition time.

Summary of the invention

The present invention proposes the displacement of tissue method of estimation of two kinds of yardsticks of a kind of employing in order to solve the contradiction of the decrement that tissue is applied large scale and small scale bigger the time.When the tissue compression amount is bigger, two kinds of methods selecting large scale or small scale to carry out the displacement of tissue estimation are combined, can realize the pluses and minuses that these two kinds of displacement of tissue are estimated maximizing favourable factors and minimizing unfavourable ones, improve the precision that displacement of tissue is estimated.

The displacement of tissue method of estimation of two kinds of yardsticks of a kind of employing that the present invention proposes may further comprise the steps:

1. the data of taking out article one scanning line respectively from the forward and backward two-dimentional radiofrequency signal of tissue compression are made as s ₁(n) and s ₂(n), n represents the data sequence number on these two scanning lines, 1≤n≤n _Max, the maximum n of n _MaxProbing depth, the spread speed of ultrasonic waves transmitted in tissue and the sample frequency decision of radiofrequency signal by this B-mode ultrasonic apparatus device;

2. from this scan-line data s ₁(n) getting a bit of yardstick in is T _aData d ₁, its data number is U _a, U _a=round (T _a* U ₁), wherein, T _aUnit be mm, U ₁Represent the data number of the tissue correspondence of 1mm, by spread speed and the sample frequency of the radiofrequency signal decision of ultrasonic waves transmitted in tissue, round () representative rounds up rounds operation, these data d ₁Sequence number from n ₁To n ₁+ U _a-1, n ₁Can be at 1≤n ₁≤ U _aScope in select; At τ ₁To τ ₂The hunting zone in ask this little segment data and scan-line data s ₂(n) cross-correlation function R (τ), computing formula is as follows

$R^{\′}\left(\mathrm{\τ}\right)=\frac{\underset{i=n_1}{\overset{n_1+U-1}{\mathrm{\Σ}}}{s}_1\left(i\right)s_2(i-\mathrm{\τ})}{\sqrt{\underset{i=n_1}{\overset{n_1+U-1}{\mathrm{\Σ}}}{s}_1^2\left(i\right)\·\underset{i=n_1}{\overset{n_1+U-1}{\mathrm{\Σ}}}{s}_2^2(i-\mathrm{\τ})}}({\mathrm{\τ}}_1\≤\mathrm{\τ}\≤{\mathrm{\τ}}_2)$

Wherein i is the cyclic variable of computational process invading the exterior registration according to sequence number, τ ₁Be 0, τ ₂Be the decrement that tissue is applied, represent with the number of sampled data; In order to improve the precision of Displacement Estimation, preferably also need to carry out interpolation processing, as conventional interpolation methods such as parabola interpolations to calculating cross-correlation function;

3. the position t ' of the cross-correlation function R ' maximum correspondence (τ) that obtains of determining step 2 ₁, t ' ₁For utilizing yardstick T _aData d when carrying out Displacement Estimation ₁Displacement after tissue compression (is s ₁(n) sequence number in is from n ₁To n ₁+ U _a-1 segment data d ₁After tissue compression, move to s ₂(n) sequence number in is from n ₁-t ' ₁To n ₁+ U _a-1-t ' ₁The position);

4. from this scan-line data s ₁(n) the data segment d in ₁Middle position to get yardstick again be T _b, sequence number from $n_1+\mathrm{round}\left(\frac{U_a-U_b}{2}\right)$ Arrive $n_1+\mathrm{round}\left(\frac{U_a-U_b}{2}\right)+U_b-1$ A bit of data e ₁, this data number is U _b, U _b=round (T _b* U ₁), wherein, T _bUnit be mm, and satisfy T _b＜T _a, preferably get $\frac{1}{16}{T}_a\≤T_b\≤\frac{1}{2}{T}_a;$ At τ ₃To τ ₄Ask this little segment data and scan-line data s in the hunting zone of determining ₂(n) " (τ), computing formula is as follows for cross-correlation function R

$R^{\′\′}\left(\mathrm{\τ}\right)=\frac{\underset{i=m_1}{\overset{m_1+U_b-1}{\mathrm{\Σ}}}{s}_1\left(i\right)s_2(i-{t^{\′}}_1-\mathrm{\τ})}{\sqrt{\underset{i=m_1}{\overset{m_1+U_b-1}{\mathrm{\Σ}}}{s}_1^2\left(i\right)\·\underset{i=m_1}{\overset{m_1+U_b-1}{\mathrm{\Σ}}}{s}_2^2(i-{t^{\′}}_1-\mathrm{\τ})}}({\mathrm{\τ}}_3\≤\mathrm{\τ}\≤{\mathrm{\τ}}_4)$

τ wherein ₃And τ ₄Be the scope of the meticulous adjustment of displacement of tissue, satisfy following relation: τ ₄＜τ ₂, τ ₃=-τ ₄, preferably get $\frac{1}{16}{\mathrm{\τ}}_2{\≤\mathrm{\τ}}_4\≤\frac{1}{2}{\mathrm{\τ}}_2;$ In order to improve the precision of Displacement Estimation, preferably also carry out interpolation, as conventional interpolation methods such as parabola interpolations to calculating cross-correlation function;

5. the cross-correlation function R " the position t of maximum correspondence (τ) " that obtains of determining step 4 ₁, t " ₁For utilizing yardstick T _bData e when carrying out Displacement Estimation ₁Displacement with respect to t ' ₁Skew;

6. utilize yardstick T _bThe displacement bias t that estimation obtains " ₁To yardstick T _aThe displacement t ' that estimation obtains ₁Adjust the segment data d that obtains ₁Meticulous adjusted displacement t ₁, that is:

t ₁＝t′ ₁+t″ ₁

7. successively from scan-line data s ₁(n) getting a bit of yardstick in is T _aBe that the data number is U _aData d ₂, d ₃..., d _N, the sequence number of every segment data V the sampled data that stagger successively, 1≤V≤U _a, V sampled data will exceed s up to staggering again ₁(n) scope, the method that 2-6 is identical obtains the displacement t of each segment data correspondence successively set by step ₂, t ₃..., t _N, wherein N is the sum of little segment data; Then displacement sequence t ₁, t ₂..., t _NBe article one scan-line data s ₁(n) Displacement Estimation of Dui Ying tissue;

8. utilize the method identical with step 1-7, obtain the 2nd, 3 successively ..., M bar scan-line data correspondence the Displacement Estimation of tissue, wherein M is the scanning line sum of expression probe, by the probe decision.

Principle of the present invention:

The present invention adopts the tracking wave band of two kinds of yardsticks to carry out Displacement Estimation, uses the tracking wave band T of large scale earlier _aCarry out the first estimation of displacement, use the tracking wave band T of small scale then _bCarry out fine tuning.The mode that the thought of the method for two kinds of yardsticks of this employing has been imitated the eye-observation things: see slightly that with large scale the reuse small scale is scrutinized earlier, just the partial thought in whole earlier back.

The present invention utilizes big yardstick to eliminate the influence of random noise, thereby has determined the approximate range of displacement of tissue, utilizes little yardstick to come displacement of tissue is accurately estimated then.Because when utilizing the tracking wave band of small scale to carry out Displacement Estimation, be on large scale is followed the tracks of the basis of estimated result of wave band, to carry out, can in a little scope, calculate cross-correlation function, so the influence of random noise is smaller.This method has been broken the constant thinking of fixed size, has solved the contradiction of the decrement that tissue is applied large scale and small scale bigger the time well, has improved the precision that displacement of tissue is estimated, does not allow to be subject to the influence of random noise again.

From the angle of error analysis, the error of Displacement Estimation can be divided into phase error and pseudo-peak error.Phase error performance is near the disturbance actual displacement value of correspondence tissue of the maximum value position of the cross-correlation function between the forward and backward radio frequency reading line data of compression, and error is less.And the maximum value position that pseudo-peak error shows as this cross-correlation function actual displacement value of the corresponding tissue of substantial deviation, the displacement of tissue of this moment estimates not comprised enough quantity of information.When the decrement that tissue is applied was bigger, yardstick was big more, did not allow to be subject to the interference of random noise more, was not easy to occur pseudo-peak error more, but because the forward and backward signal superposed part of tissue compression is more little, the phase error of this moment is bigger; And yardstick hour, and phase error is reduced, but the pseudo-peak error of this moment is bigger.

The method of two kinds of yardsticks of employing that the present invention proposes is exactly to utilize large scale to reduce pseudo-peak error earlier, utilizes small scale to reduce phase error then, thereby combines the advantage of two kinds of yardsticks, has reduced the interference of error, has improved the precision that displacement of tissue is estimated.

Characteristics of the present invention:

1 adopts the tracking wave band of two kinds of different scales, utilizes the meticulous again adjustment of first rough estimate, elder generation's partial thinking in arrangement back to carry out Displacement Estimation;

2 carry out the rough estimate of displacement earlier with the tracking wave band of large scale, the tracking wave band with small scale carries out meticulous adjustment then;

3 results that utilize meticulous adjustment revise the result of rough estimate, estimate thereby obtain more accurate displacement of tissue.

Description of drawings

The computer artificial result of organizing stress distribution that Fig. 1 obtains for general displacement of tissue method of estimation;

The computer artificial result of organizing stress distribution that the embodiment of the displacement of tissue method of estimation of two kinds of yardsticks of employing that Fig. 2 proposes for the present invention obtains.

The specific embodiment

The displacement of tissue method of estimation of two kinds of yardsticks of employing that the present invention proposes reaches accompanying drawing in conjunction with specific embodiments and is described in detail as follows:

Embodiments of the invention utilize computer program and general ultrasonic scattering model emulation obtain one mimic be organized in the compression before and the compression after the one dimension radiofrequency signal.Tissue thickness is 60mm, and wherein the elastic modelling quantity at 0-30mm place is 3 times of elastic modelling quantity at 30-60mm place, adopt simple one-dimensional model, so the strain at 0-30mm place is strained 1/3 of 30-60mm place.The decrement that tissue is applied is 3.33%, i.e. 1.98mm; The center probe frequency is 3.5MHz, and-three dB bandwidth is 2.0MHz, because adopt simple one-dimensional model, so have only a scanning line; The sample frequency of radiofrequency signal is 20MHz, supposes that ultrasound wave is 1540m/s in in-house spread speed, so the corresponding 1mm/ (1540 * 10 of the tissue of 1mm ³Mm/s/ (20 * 10 ⁶Hz)/2 26 data of ≈, because tissue depth is 60mm, so the data of each bar scanning line are 60 * 26=1560, the decrement that tissue is applied is expressed as 52 sampled datas of 60 * 3.33% * 26 ≈ with the number of sampled data.

The large scale T that present embodiment adopts _aBe 4mm, the hunting zone is 0 to 52; The small scale T that adopts _bBe 1mm (also can adopt 2mm or any yardstick in 1mm to 2mm scope) that the hunting zone (also can be-20 to-20 or satisfied 10≤τ for-10 to 10 ₄≤ 26 and τ ₃=-τ ₄From τ ₃To τ ₄Any scope).

The concrete steps of present embodiment are as follows:

1. establish the forward and backward one dimension radiofrequency signal (being scan-line data) of tissue tissue compression and be respectively s ₁(n) and s ₂(n); N represents the data sequence number on these two scanning lines, 1≤n≤1560;

2. from this scan-line data s ₁(n) getting a bit of yardstick in is T _aData d ₁, T _a=4mm, its data number is U, U=104, the sequence number of these data from 32 to 135; In 0 to 52 hunting zone, ask this little segment data and scan-line data s ₂(n) cross-correlation function R (τ), computing formula is as follows

$R^{\′}\left(\mathrm{\τ}\right)=\frac{\underset{i=n_1}{\overset{n_1+U-1}{\mathrm{\Σ}}}{s}_1\left(i\right)s_2(i-\mathrm{\τ})}{\sqrt{\underset{i=n_1}{\overset{n_1+U-1}{\mathrm{\Σ}}}{s}_1^2\left(i\right)\·\underset{i=n_1}{\overset{n_1+U-1}{\mathrm{\Σ}}}{s}_2^2(i-\mathrm{\τ})}}(0\≤\mathrm{\τ}\≤52)$

3. determine the position t ' of this cross-correlation function R ' maximum correspondence (τ) ₁, t ' ₁Utilize yardstick T exactly _aData d when carrying out Displacement Estimation ₁Displacement after tissue compression (is s ₁The segment data d of the sequence number (n) from 32 to 135 ₁After tissue compression, move to s ₂(n) sequence number in is from 32-t ' ₁To 135-t ' ₁The position);

4. get s ₁The a bit of data e of the sequence number (n) from 71 to 96 ₁, these data e ₁Be positioned at data d ₁Middle position; These data e ₁Yardstick be T _b, T _b=1mm, the data number is U _b, U _b=26; In-10 to 10 hunting zone, ask this little segment data and scan-line data s ₂(n) " (τ), computing formula is as follows for cross-correlation function R

$R^{\′\′}\left(\mathrm{\τ}\right)=\frac{\underset{i=m_1}{\overset{m_1+U_b-1}{\mathrm{\Σ}}}{s}_1\left(i\right)s_2(i-{t^{\′}}_1-\mathrm{\τ})}{\sqrt{\underset{i=m_1}{\overset{m_1+U_b-1}{\mathrm{\Σ}}}{s}_1^2\left(i\right)\·\underset{i=m_1}{\overset{m_1+U_b-1}{\mathrm{\Σ}}}{s}_2^2(i-{t^{\′}}_1-\mathrm{\τ})}}(-10\≤\mathrm{\τ}\≤10)$

5. the cross-correlation function R " the position t of maximum correspondence (τ) " that obtains of determining step 4 ₁, t " ₁Utilize yardstick T exactly _bData e when carrying out Displacement Estimation ₁Displacement with respect to t ' ₁Skew;

6. utilize yardstick T _bThe displacement bias t that estimation obtains " ₁To yardstick T _aThe displacement t ' that estimation obtains ₁Adjust the segment data d that obtains ₁Meticulous adjusted displacement t ₁, that is:

t ₁＝t′ ₁+t″ ₁

7. successively from scan-line data s ₁(n) getting a bit of yardstick in is T _aBe that the data number is U _aData d ₂, d ₃..., d _N, the sequence number of every segment data 32 sampled datas that stagger successively, 32 sampled datas will exceed s up to staggering again ₁(n) scope, the method that 2-6 is identical obtains the displacement t of each segment data correspondence successively set by step ₂, t ₃..., t _N, wherein N is the sum of little segment data, N=45; Then displacement sequence t ₁, t ₂, t _NBe scan-line data s ₁(n) Displacement Estimation of Dui Ying tissue.

The Displacement Estimation effect of present embodiment and conventional method is compared as follows:

Fig. 1 and Fig. 2 are the result's of the result of present embodiment and conventional method comparison.Among Fig. 1 and Fig. 2, abscissa is represented the diverse location organized, and vertical coordinate is represented strained size; The actual strain of dotted line 11 among Fig. 1 and the 21 expression tissues of the dotted line among Fig. 2 distributes; Among Fig. 1 12 estimates the stress distribution of organizing that obtains for adopting conventional method, and 22 among Fig. 2 estimates the stress distribution of organizing that obtains for the method for two kinds of yardsticks of employing of adopting the present invention and proposing.

By Fig. 1 and Fig. 2 as seen, the method of two kinds of yardsticks of employing that the present invention proposes adopts the tracking wave band of two kinds of different scales to carry out Displacement Estimation, the result is better than adopting the constant tracking section of fixed size to carry out the result of conventional method of Displacement Estimation, the precision of Displacement Estimation is higher, is particularly organizing the bigger zone of strain ratio (30-60mm position) more obvious.

Claims (3)

1, the displacement of tissue method of estimation of two kinds of yardsticks of a kind of employing may further comprise the steps:

1) data of taking out article one scanning line respectively from the forward and backward two-dimentional radiofrequency signal of tissue compression are made as s ₁(n) and s ₂(n), n represents the data sequence number on these two scanning lines, 1≤n≤n _Max, n _MaxMaximum for n;

2) from this scan-line data s ₁(n) getting a bit of yardstick in is T _aData d ₁, its data number is U _a, U _a=round (T _a* U ₁), wherein, T _aUnit be mm, U ₁Represent the data number of the tissue correspondence of 1mm, round () representative rounds up rounds operation, these data d ₁Sequence number from n ₁To n ₁+ U _a-1,1≤n ₁≤ U _aAt τ ₁To τ ₂The hunting zone in ask this little segment data and scan-line data s ₂(n) cross-correlation function R (τ), computing formula is as follows

$R^{\′}\left(\mathrm{\τ}\right)=\frac{\underset{i=n_1}{\overset{n_1+U-1}{\mathrm{\Σ}}}{s}_1\left(i\right)s_2(i-\mathrm{\τ})}{\sqrt{\underset{i=n_1}{\overset{n_1+U-1}{\mathrm{\Σ}}}{s}_1^2\left(i\right)\·\underset{i=n_1}{\overset{n_1+U-1}{\mathrm{\Σ}}}{s}_2^2(i-\mathrm{\τ})}}({\mathrm{\τ}}_1\≤\mathrm{\τ}\≤{\mathrm{\τ}}_2)$

Wherein i is the cyclic variable of computational process invading the exterior registration according to sequence number, τ ₁Be 0, τ ₂Be the decrement that tissue is applied, represent with the number of sampled data;

3) determining step 2) the position t ' of the cross-correlation function R ' maximum correspondence (τ) that obtains ₁, t ' then ₁For utilizing yardstick T _aData d when carrying out Displacement Estimation ₁Displacement after tissue compression (is s ₁(n) sequence number in is from n ₁To n ₁+ U _a-1 segment data d ₁After tissue compression, move to s ₂(n) sequence number in is from n ₁-t ' ₁To n ₁+ U _a-1-t ' ₁The position);

4) again from this scan-line data s ₁(n) the data segment d in ₁Middle position to get yardstick be T _b, sequence number from

$n_1+\mathrm{round}\left(\frac{U_a-U_b}{2}\right)$ Arrive $n_1+\mathrm{round}\left(\frac{U_a-U_b}{2}\right)+U_b-1$ A bit of data e ₁, this data number is U _b, U _b=round (T _b* U ₁), wherein, T _bUnit be mm, and satisfy T _b＜T _aAt τ ₃To τ ₄The hunting zone in ask this little segment data and scan-line data s ₂(n) " (τ), computing formula is as follows for cross-correlation function R

$R^{\′\′}\left(\mathrm{\τ}\right)=\frac{\underset{i=m_1}{\overset{m_1+U_b-1}{\mathrm{\Σ}}}{s}_1\left(i\right)s_2(i-{t^{\′}}_1-\mathrm{\τ})}{\sqrt{\underset{i=m_1}{\overset{m_1+U_b-1}{\mathrm{\Σ}}}{s}_1^2\left(i\right)\·\underset{i=m_1}{\overset{m_1+U_b-1}{\mathrm{\Σ}}}{s}_2^2(i-{t^{\′}}_1-\mathrm{\τ})}}({\mathrm{\τ}}_3\≤\mathrm{\τ}\≤{\mathrm{\τ}}_4)$

τ wherein ₃And τ ₄Be the scope of the meticulous adjustment of displacement of tissue, satisfy following relation: τ ₄＜τ ₂

5) the cross-correlation function R " the position t of maximum correspondence (τ) " that obtains of determining step 4 ₁, t then " ₁For utilizing yardstick T _bData e when carrying out Displacement Estimation ₁Displacement with respect to t ' ₁Skew;

6) utilize yardstick T _bThe displacement bias t that estimation obtains " ₁To yardstick T _aThe displacement t ' that estimation obtains ₁Adjust the segment data d that obtains ₁Meticulous adjusted displacement t ₁, that is:

t ₁＝t′ ₁+t″ ₁

7) successively from scan-line data s ₁(n) getting a bit of yardstick in is T _aBe that the data number is U _aData d ₂, d ₃..., d _N, the sequence number of every segment data V the sampled data that stagger successively, 1≤V≤U _a, V sampled data will exceed s up to staggering again ₁(n) scope, the method that 2-6 is identical obtains the displacement t of each segment data correspondence successively set by step ₂, t ₃..., t _N, wherein N is the sum of little segment data; Then displacement sequence t ₁, t ₂..., t _NBe article one scan-line data s ₁(n) Displacement Estimation of Dui Ying tissue;

8) utilize the method identical with step 1-7, obtain the 2nd, 3 successively ..., M bar scan-line data correspondence the Displacement Estimation of tissue, wherein M is the scanning line sum of expression probe.

2, as claimed in claim 1, it is characterized in that, also further be included in the described the 2nd), 4) the described cross-correlation function that calculates is carried out interpolation processing in the step.

3, as claimed in claim 1, it is characterized in that the described the 4th) yardstick T in the step _bSpan be:

$\frac{1}{16}{T}_a\≤T_b\≤\frac{1}{2}{T}_a;$ The τ of described hunting zone ₄Value is $\frac{1}{16}{\mathrm{\τ}}_2\≤{\mathrm{\τ}}_4\≤\frac{1}{2}{\mathrm{\τ}}_2.$