Method and apparatus for elastographic measurement and imaging

• US 5,474,070 A
• Filed: 03/11/1994
• Issued: 12/12/1995
• Est. Priority Date: 11/17/1989
• Status: Expired due to Term

• Current Assignee: Board of Regents of the University of Texas System (University of Texas System)
• Original Assignee: Board of Regents of the University of Texas System (University of Texas System)
• Inventors: Cespedes, Ignacio, Ophir, Jonathan
• Primary Examiner(s): JAWORSKI, FRANCIS J
• Application Number:
  US08/212,734
• Time in Patent Office: 641 Days
• Field of Search: 128/660.01, 128/660.06, 128/660.07, 128/661.03, 128/660.08, 73/597
• US Class Current: 600/437
• CPC Class Codes: A61B 5/0048 Detecting, measuring or recording by applying mechanical forces or stimuli A61B5/021, A61B5/4884 and A61B9/00 take precedence
• CPC Class Codes: G01N 2203/0089 Biorheological properties
• CPC Class Codes: G01N 2291/02827 Elastic parameters, strength or force
• CPC Class Codes: G01N 2291/101 one transducer
• CPC Class Codes: G01N 29/07 by measuring propagation velocity or propagation time of acoustic waves
• CPC Class Codes: G01N 29/223 Supports, positioning or alignment in fixed situation mounting transducers per se G10K11/004
• CPC Class Codes: G01N 29/30 Arrangements for calibrating or comparing, e.g. with standard objects
• CPC Class Codes: G01N 29/38 by time filtering, e.g. using time gates
• CPC Class Codes: G01N 29/50 using auto-correlation techniques or cross-correlation techniques
• CPC Class Codes: G01S 15/8945 using transducers mounted for linear mechanical movement
• CPC Class Codes: G01S 7/52036 using analysis of echo signal for target characterisation
• CPC Class Codes: G01S 7/52042 determining elastic properties of the propagation medium or of the reflective target diagnostic techniques involving the measurement of strain A61B8/485

Abstract

An improved ultrasonic pulse-echo method and apparatus that has particular application in making precision measurements of compressibility in any backscattering material, in particular organic tissue, is disclosed. The method employs a standard transducer or transducer containing device which is translated transaxially, thereby compressing or displacing a proximal region of a target body in small known increments. At each increment, a pulse is emitted and an echo sequence (A-line) is detected from regions within the target along the sonic travel path or beam of the transducer. Resulting time shifts in echo segments corresponding to features in the target, corrected for regions of varying sonic speed along the sonic path, provide relative and quantitative information concerning the strain caused by the compressions. The stress imparted by the transducer and containing device is also determined, corrected for depth along the sonic path. The appropriate values for stress are divided into the respective values for strain along each path to yield an elastogram, or array of compressibility values, of the target.

Claims

1. A method for taking elastographic measurements of a heterogeneous target body comprising the steps of:

(a) sonically coupling an ultrasonic source to the surface of a target body;

(b) emitting a first pulse of ultrasonic energy from the source along a radiation axis in the target body;

(c) recording the arrival time of the first echo sequence having at least one echo segment arriving in response to said first pulse of ultrasonic energy;

(d) displacing the surface of the target body a known displacement, Δ

X, along a force axis so as to compress the target body;

(e) emitting a next pulse of ultrasonic energy from the source along said radiation axis in the target body following said movement;

(f) recording the arrival time of a next echo sequence arriving in response to said next pulse of ultrasonic energy, said echo sequence having at least one echo segment that is congruent with at least one said echo segment that arrived in response to the preceding pulse;

(g) measuring the differential displacement from at least one congruent echo segment of the two most recent pulses of ultrasonic energy;

(h) calculating the strain from the differential displacement data obtained in step (g);

(i) storing strain data that meets a predetermined quality standard in a retrievable medium;

(j) normalizing the calculated strain data stored in step (i) using the formula;

NS_n =CS_n /n where;

NS_n =normalized strain at displacement interval n

CS_n =calculated strain at displacement interval n; and

(k) repeating steps (d)-(j) until enough strain data has been stored to produce an image of the tissue body region under interrogation.

2. The method of claim 1, wherein said target body is a tissue body.

3. The method of claim 1, wherein the differential displacement is measured by a cross-correlation technique.

4. The method of claim 1, wherein the differential displacement is measured by a phase of the analytic signal.

5. The method of claim 1, wherein the differential displacement is measured by a phase of the cross-spectrum.

6. The method of claim 1, wherein said sonic coupling is achieved by placing an adjustable volume coupling device between the transducer and the surface of the target body.

7. The method of claim 6, wherein said coupling device is volume adjustable by inflation and deflation.

8. The method of claim 6, wherein the surface of the target body is displaced by changing the volume of the coupling device.

9. The method of claim 1, further comprising the steps of:

(a) displacing the surface of the target body a known displacement, AX, along the force axis so as to decompress the target body;

(b) emitting a next pulse of ultrasonic energy from the source along the radiation axis in the target body following said movement;

(c) recording the arrival time of a next echo sequence arriving in response to said next pulse of ultrasonic energy, said echo sequence having at least one echo segment that is congruent with at least one said echo segment that arrived in response to the preceding pulse;

(d) measuring the differential displacement from at least one congruent echo segment of the two most recent pulses of ultrasonic energy;

(e) calculating the strain from the differential displacement data obtained in step (d);

(f) storing strain data that meets a predetermined quality standard in a retrievable medium;

(g) normalizing the calculated strain data stored in step (f) using the formula;

NS_n =CS_n /n where;

NS_n =normalized strain at displacement interval n

CS_n =calculated strain at displacement interval n; and

(h) repeating steps (a)-(g) until the number of decompression intervals equal the number of compression intervals.

10. The method of claim 9, wherein each displacement of the surface of the target body results in an elastic wave traveling within the target body.

11. The method of claim 10, wherein the elapsed time between the emission of the first and last pulse of ultrasonic energy from the source is large enough to allow the elastic wave to travel to the desired depth of interrogation within the target body.

12. The method of claim 1 further comprising the step of calculating the displacement shift of the interrogated target body region from its initial position at each displacement interval, Δ

TB_n, using the following formula;

##EQU6## where;

Δ

TB_n is the interrogated target body region displacement at displacement interval n;

X_o is the initial displacement of the target body surface along the force axis;

X_n is the displacement of the target body surface along the force axis at displacement interval n;

NS_n is the normalized strain at displacement interval n; and

dx is the differential displacement of the target body surface along the force axis.

13. The method of claim 12, wherein the equation for calculating the displacement shift of the interrogated target body region is solved numerically.

14. The method of claim 1, wherein said ultrasound source comprises an array of transducers capable of emitting ultrasound energy in a preselected sequence at a preselected emission interval.

15. A method for taking elastographic measurements of a heterogeneous target body comprising the steps of:

(a) sonically coupling an ultrasonic source to a target body;

(b) emitting a first pulse of ultrasonic energy from the source along a radiation axis in the target body;

(c) recording the arrival time of the first echo sequence having at least one echo segment arriving in response to said first pulse of ultrasonic energy;

(d) moving the ultrasonic source a known displacement, Δ

X, along a force axis so as to compress the target body;

(e) emitting a next pulse of ultrasonic energy from the source along the radiation axis in the target body following said movement;

(f) recording the arrival time of a next echo sequence arriving in response to said next pulse of ultrasonic energy, said echo sequence having at least one echo segment that is congruent with at least one said echo segment that arrived in response to the preceding pulse;

(g) measuring the differential displacement from at least one congruent echo segment of the first pulse and most recent pulse of ultrasonic energy;

(h) calculating the strain from the differential displacement data obtained in step (g);

(i) storing strain data that meets a predetermined quality standard in a retrievable medium;

(j) normalizing the calculated strain data stored in step (i) using the formula;

NS_n =CS_n /n where;

NS_n =normalized strain at displacement interval n

CS_n =calculated strain at displacement interval n; and

(k) repeating steps (d)-(j) until the ultrasonic source has traveled a desired compression distance;

(l) moving the ultrasonic source a known displacement, Δ

X, along the force axis so as to decompress the target body;

(m) emitting a next pulse of ultrasonic energy from the source along the radiation axis in the target body following said movement;

(n) recording the arrival time of a next echo sequence arriving in response to said next pulse of ultrasonic energy, said echo sequence having at least one echo segment that is congruent with at least one said echo segment that arrived in response to the preceding pulse;

(o) measuring the strain from at least one congruent echo segment of the first pulse and most recent pulse of ultrasonic energy;

(p) calculating the strain from the differential displacement data obtained in step (o);

(q) storing strain data that meets a predetermined quality standard in a retrievable medium;

(r) normalizing the calculated strain data stored in step (q) using the formula;

NS_n =CS_n /n where;

NS_n =normalized strain at displacement interval n

CS_n =calculated strain at displacement interval n; and

(s) repeating steps (l)-(r) until the ultrasonic source has traveled a desired decompression distance.

16. The method of claim 15, wherein the displacement is measured by a cross-correlation technique.

17. The method of claim 15, wherein the strain is measured by a phase of the analytic signal.

18. The method of claim 15, wherein the strain is measured by a phase of the cross-spectrum.

19. The method of claim 15, further comprising the step of displaying an image of the stored strain data.

20. The method of claim 15, wherein the compression distance that the ultrasonic source is moved is equal to the decompression distance that the ultrasonic source was moved.