3D QUANTITATIVE-IMAGING ULTRASONIC METHOD FOR BONE INSPECTIONS AND DEVICE FOR ITS IMPLEMENTATION

• US 20100185089A1
• Filed: 04/07/2009
• Published: 07/22/2010
• Est. Priority Date: 10/24/2006
• Status: Abandoned Application

• Current Assignee: Alla Gourevitch
• Original Assignee: Alla Gourevitch
• Inventors: Gourevitch, Alla
• Application Number:
  US12/444,599
• Publication Number:
  US 20100185089A1
• Time in Patent Office: Days
• US Class Current: 600/443
• CPC Class Codes: A61B 5/417 the bone marrow
• CPC Class Codes: A61B 8/0875 for diagnosis of bone A61B5/4504 takes precedence
• CPC Class Codes: A61B 8/13 Tomography A61B8/10, A61B8/12 take precedence; tomography for radiation diagnosis A61B6/02
• CPC Class Codes: A61B 8/4477 using several separate ultrasound transducers or probes
• CPC Class Codes: A61B 8/4488 the transducer being a phased array
• CPC Class Codes: A61B 8/4494 characterised by the arrangement of the transducer elements devices for short-range imaging using particular transducer elements arrangements under G01S15/8909; arrangements of transducers in generation of mechanical vibrations of infrasonic, sonic or ultrasonic frequency B06B1/0607
• CPC Class Codes: A61B 8/483 involving the acquisition of a 3D volume of data
• CPC Class Codes: A61B 8/485 involving measuring strain or elastic properties
• CPC Class Codes: A61B 8/5223 for extracting a diagnostic or physiological parameter from medical diagnostic data for algorithms to analyse biomedical images G06T7/0012
• CPC Class Codes: G16H 50/30 for calculating health indices; for individual health risk assessment

Abstract

The Ultrasonic Tomographical method and system is provided using measurements of time of flight low frequency acoustic waves. Differences in first signal arrival times from plurality of known transmitters'"'"' locations to plurality of known receivers'"'"' location are used, wherein the transmitters and receivers are at an angle to the surface of the observed object. 3D mapping of the acoustic propagation speed is reconstructed, revealing anatomical details and physiological properties.

Claims

1. A method of 3D Quantitative-Imaging Ultrasonic Tomography for inspecting a heterogeneous object;

wherein said method comprises the step of;

a. providing a 3D quantitative imaging ultrasound tomography system, said system comprises at least;

i. a three dimensional ultrasonic unit characterized by;

(a) a grid array of evenly spaced ultrasonic transducers capable of transmitting an ultrasonic wave at an angle to said grid in response to excitation pulses, and capable of producing a signals in response to received ultrasonic waves at an angle to said grid;

said transducer grid is in acoustic contact with said inspected object; and

,

(b) a layered system comprising the inspected object, said layered system is characterized by acoustic impedance gradient causing the ultrasonic waves to propagate through said system and the object along a non-linear paths;

ii. a signal generator generating short excitation pulses;

iii. a scanning position controller adapted for consecutively emitting said generated excitation pulses and directing said pulses to a selected transmitting transducer in said transducer'"'"'s grid array and receiving signals created in other receiving transducers of grid'"'"'s array in response to said ultrasonic waves emitted by said transmitted transducer and propagating along the non-linearly paths passed through said inspected system according to a predetermined protocol;

iv. a measuring time unit receiving said signals from said receiving transducers and measuring a time of wave travel between corresponding pair of transmitting transducer and receiving transducer;

v. a processor adapted for acquiring a plurality of measured travel times corresponding to a plurality of paths between said transmitting and receiving transducers;

calculating according to the differential approach plurality of time values corresponding to plurality elementary volumes composing of said object and calculating length of elastic wave paths in elementary cells and calculating of said longitudinal wave velocity and porosity corresponding to a plurality of travel times further corresponding to each combination of a pair of adjacent transmitting transducers and a pair of adjacent receiving transducers for direct and reciprocal directions; and

evaluating longitudinal wave velocity in material matrix part of said heterogeneous object;

vi. an image formation unit;

vii. memory

(b) providing non-linear ultrasonic waves propagation through the layered system by;

i. establishing acoustical contact between said grid of transducers and a surface of said layered system in a unilateral way; and

ii. consecutively transmitting ultrasonic waves by each transducer and receiving ultrasonic waves by the rest of transducers of said grid'"'"'s array said ultrasonic waves being transmitted and received at an angle to the layered system;

(c) measuring travel times of said ultrasonic waves transmitted and received at said step of consecutively transmitting ultrasonic waves;

(d) interpreting said layered system containing inspected object as a plurality of elementary cells arranged in columns and rows;

(e) by means of differential approach calculating travel times corresponding to each elementary cell;

i. calculating changes in travel times Δ

t_m corresponding to each subsequent cell relatively to the previous cell along each column by combining the average values τ

_m, τ

<sub>m−

1</sub>, τ

_dir, τ

_rec from eight travel times of the direct and reciprocal directions according to equation;

Δ

τ

_m=τ

_m+τ

<sub>m−

1</sub>−

τ

_dir−

τ

_rec, where τ

_m is the average value between longitudinal wave travel time for direct and reciprocal directions for the transducers arrangement from points m to point (−

)m, were m is a number of a transducer location;

τ

<sub>m−

1</sub> is an average value between longitudinal wave travel time for direct and reciprocal directions for transducers arrangement from point (m−

1) to point (−

m+1);

τ

_dir is an average value between longitudinal wave travel time for direct and reciprocal directions for transducers arrangement from point (−

m) to point (m−

1) and τ

_rec is an average value between longitudinal wave travel time for direct and reciprocal directions for transducers arrangement from point m to point (−

m+1);

ii. calculating a sequence of travel time values corresponding to each of said elementary cell of said column by summing said values of said changes travel times in said elementary cells Δ

τ

_m (Δ

τ

₁, Δ

τ

₂, Δ

τ

₃, Δ

τ

₄ . . . Δ

τ

_m) and the travel times in said previous elementary cell in said column t_m (t₁, t₂, t₃, . . . t_m) according to equation;

τ

₂=t₁+Δ

τ

₁, τ

₃=t₂+Δ

τ

₂, τ

₃=t₂+Δ

τ

₂ . . . and τ

_m=t<sub>m−

1</sub>+Δ

τ

<sub>m−

1</sub>;

(f) calculating of longitudinal wave velocity values associated with ultrasonic waves propagating in each elementary cell, said calculating is carried out by dividing a length of a beam travel in the elementary cell by said travel time calculated at the step (e).ii;

wherein said travel length associated with a first cells of a columns is equal to 2b/sin α and

said travel length associated with the rest of cells of the same columns is equal to π

b, where b is a distance between adjacent transducers and α

is the incident angle on the layered system surface;

(g) statistical evaluating of longitudinal wave velocity value associated with ultrasonic waves propagating in the material matrix portion of said inspected object, said evaluating is carried out by means of histogramming of obtained longitudinal wave velocity values corresponding to said plurality of elementary cells;

(h) calculating of porosity values n for said plurality of elementary cells according to the formula;

2. The method of the claim 1, in which said object is a bone and said method further comprising estimating a bone mineral matrix portion and porous portions of human body parts and mapping distributions of said longitudinal wave velocity values and porosity in said body parts for diagnosing bone diseases.

3. The method of the claim 1, wherein said step of evaluating longitudinal wave velocity value in said material matrix portion of said inspected object further comprises histogramming said obtained longitudinal wave velocity values corresponding to said plurality of said elementary cells by maximizing thereof.

4. The method of the claim 3, comprising evaluating the mineral matrix part of said bone by a maximum value of said longitudinal wave velocity obtained by means of histogramming.

5. The method of the claim 1, wherein said step of pixel-by-pixel 2D or 3D image mapping distributions of said longitudinal wave velocity and porosity values further comprises diagnosing of bone and surrounding soft tissues according to inspected object material and structural properties, geometry, sizes, micro-architecture, bone mineral matrix condition, anatomical configurations of a patient'"'"'s organs.

6. The method of the claim 1, further comprising step of diagnosing bone diseases according to temporal changes of said material and structural properties, geometry, sizes, micro-architecture, bone matrix condition, anatomical configurations of a patient'"'"'s organs.

7. The method of the claim 1, wherein said step of diagnosing diseases is applicable to an object selected from the group consisting of a cortical bone, a trabecular bone, a bone marrow tissue;

a cartilage, and any combination thereof.

8. The method of claim 1, wherein said step of consecutively transmitting and receiving ultrasonic waves is performed at an angle to said inspected layered system selected between 0 up to 90°

.

9. The method of claim 1, used for inspecting of heterogeneous materials, wherein ultrasonic wave frequency is chosen according to condition that a length λ

_wave of said ultrasonic wave used for measurements is commensurable with the dimension d of said inspected object accommodated in the layered system such that inspected object in the layered system is considered as local heterogeneity.

10. The method of claim 1 further comprises the step of calculating travel times changes Δ

t_m for each subsequent cell relative to previous cell of the column between the corresponding four travel times for one direction τ

_m, τ

<sub>m−

1</sub>, τ

_dir, and τ

_rec, where τ

_m, is a value of longitudinal wave velocity travel between the m-th and (−

m)-th transducers, τ

<sub>m−

1</sub> a value of longitudinal wave velocity travel between the (m−

1)-th and (−

m+1)-th transducers, τ

_dir a value of longitudinal wave velocity travel between the m-th and (−

m+1)-th transducers, τ

_rec a value of longitudinal wave velocity travel between the (m−

1)-th and (−

m)-th transducers;

enumerating the transducers is performed relative to an axis of the column;

the time changes Δ

t_m are calculated for each combination a pair of adjacent transmitting transducers and a pair of adjacent receiving transducers in direct and reciprocal directions according to the equation;

Δ

t=τ

_m+τ

<sub>m−

1</sub>−

τ

_dir−

τ

_rec,

and travel time in an elementary cell;

τ

_m=t<sub>m−

1</sub>+Δ

τ

<sub>m−

1</sub>.

21. The method of the claim 1, further comprising step of analyzing the obtained maps thereby providing bone diagnostics according to bone and soft tissues quantitative material and structure properties estimations by anatomical pictures revealing geometry size, micro-architecture, bone mineral matrix condition, anatomical configurations of a patient organs.

22. The method of the claim 1, wherein additionally comprising step of providing natural layered system containing heterogeneous object and characterized by acoustic impedance gradient;

said system provides non-linear ultrasound beam travel in an inspected object.

23. The methods of the claim 1, further comprising step of integrating a number of 2D images performed along a third axis thereby providing said 3D image mapping distributions of longitudinal wave velocity and porosity values in said inspected heterogeneous object.

11. A method of 3D quantitative—

Imaging Ultrasonic Tomography for inspecting of a homogeneous object;

wherein said method comprising the step of;

(a) providing a 3D quantitative imaging ultrasound tomography system, said system comprises at least;

i. a three dimensional ultrasonic unit characterized by;

a. a grid array of evenly spaced ultrasonic transducers capable of transmitting an ultrasonic wave at an incident angle to said grid in response to excitation pulses, and capable of producing a signals in response to received ultrasonic waves at an angle to said grid;

said transducer grid array is in acoustic contact with said inspected homogeneous object; and

,

b. a layered system adapted for containing the inspected object;

the system is characterized by acoustic impedance gradient to provide non-linear beam paths in the heterogeneous object;

ii. a signal generator generating short excitation pulse;

iii. a scanning position controller adapted for consecutively emitting said generated excitation pulses and directing said pulses to a selected transmitting transducer in said transducer'"'"'s grid and receiving signals created in other receiving transducers in response to said ultrasonic wave emitted by said transmitted transducer and passed through said system according to a predetermined protocol;

iv. a measuring time unit capable of receiving said signals from said receiving transducers and measuring a time of wave travel between corresponding pair of transmitting transducer and receiving transducer;

v. a processor adapted for acquiring a plurality of measured travel times corresponding to a plurality of paths between said transmitting and receiving transducers;

calculating according to the differential approach plurality of time values corresponding to plurality elementary volumes composing of said object;

according to the times and length of elastic wave paths in elementary cells calculating of said longitudinal wave velocity and porosity corresponding to a plurality of travel times further corresponding to each combination of a pair of adjacent transmitting transducers and a pair of adjacent receiving transducers for direct and reciprocal directions; and

evaluating longitudinal wave velocity in material matrix part of object;

vi. an image formation unit;

vii. memory

(b) providing a refracted ultrasonic waves in said inspected object by;

i. acoustically contacting said grid of transducers to the surface of said system with inspected object; and

,

ii. consecutively transmitting ultrasonic waves by at least one of said transducers and receiving ultrasonic waves by other transducers of said grid;

said ultrasonic wave is refracted by said inspected homogeneous object;

said transmitting and receiving ultrasound waves are performed angularly to said surface of said object;

(c) measuring travel times of said ultrasonic waves transmitted and received at said step consecutively transmitting ultrasonic waves;

(d) dividing said homogeneous object into a plurality of elementary cells arranged in columns and rows;

(e) differentially calculating travel times corresponding to each elementary cell of said inspected object by means of the differential approach;

i. calculating changes in travel times Δ

t_m corresponding to each subsequent cell relatively to the previous cell along each column by combining the average values τ

_m, τ

<sub>m−

1</sub>, τ

_dir, τ

_rec from at least 8 travel times of the direct and reciprocal directions according to equation;

Δ

τ

_m=τ

_m+τ

<sub>m−

1</sub>−

τ

_dir−

τ

_rec, where τ

_m is the average value between longitudinal wave travel time for direct and reciprocal directions for the transducers arrangement from points m to point (−

m), were m is a number of a transducer location;

τ

<sub>m−

1</sub> is an average value between longitudinal wave travel time for direct and reciprocal directions for transducers arrangement from point (m−

1) to point (−

m+1);

τ

_dir is an average value between longitudinal wave travel time for direct and reciprocal directions for transducers arrangement from point (−

m) to point (m−

1) and τ

_rec is an average value between longitudinal wave travel time for direct and reciprocal directions for transducers arrangement from point m to point (−

m+1);

(ii) calculating a sequence of travel time values corresponding to each of said elementary cell of said column by summing said values of said changes travel times in said elementary cells Δ

τ

_m (Δ

τ

₁, Δ

τ

₂, Δ

τ

₃, Δ

τ

₄ . . . Δ

τ

_m) and the travel times in said previous elementary cell in said column t_m (t₁, t₂, t₃, . . . t_m) according to equation;

τ

₂=t₁+Δ

τ

₁, τ

₃=t₂+Δ

τ

₂, τ

₃=t₂+Δ

τ

₂ . . . and τ

_m=t<sub>m−

1</sub>+Δ

τ

<sub>m−

1</sub>;

f. calculating longitudinal wave velocity values corresponding to each elementary cell by dividing a length of beams travel in elementary cell by said travel time;

said length within said cells of columns is equal to 2b/sin α

where b is a distance between said transducers in said grid and α

is the incident angle and;

g. evaluating longitudinal wave velocity value in a material matrix portion of said inspected object by means of histogramming of said obtained longitudinal wave velocity values corresponding to said plurality of said elementary cells by maximizing thereof;

h. calculating porosity values n for said plurality of said cells according to the following formula;

12. The method of claim 11, further comprising the step of estimating risk of fracture in said maps of said inspected object in terms of relative changes of said longitudinal wave velocity values for said area divided by said longitudinal wave velocity in object material matrix.

13. An imaging ultrasound tomography system for inspecting of an object, wherein said system comprising:

(a) a three dimensional ultrasonic unit comprising

i. a grid of evenly spaced ultrasonic transducers capable of transmitting an ultrasonic wave at an incident angle to said grid in response to excitation pulses, and capable of producing signals in response to received ultrasonic waves by other transducers of said grid at an angle to said grid;

said transducer grid is in acoustic contact with said inspected object; and

ii. a layered system adapted for accommodating said inspected object;

said system is characterized by an acoustic impedance gradient causing the ultrasonic waves to propagate through the object along a non-linear beam paths in said heterogeneous object;

(b) a signal generator capable of generating short excitation pulses;

(c) a scanning position controller adapted for consecutively emitting said excitation pulses and directing said pulses to a selected transmitting transducer in said transducer'"'"'s grid and receiving signals created in other receiving transducers in response to said ultrasonic wave emitted by said transmitted transducer and non-linearly propagating through said inspected system according to a predetermined protocol;

(d) a measuring time unit capable of receiving said signals from said receiving transducers and measuring a time of wave travel between corresponding pair of transmitting transducer and receiving transducer;

(e) a processor adapted for (i) acquiring a plurality of measured travel times corresponding to a plurality of paths between said transmitting and receiving transducers;

(ii) calculating according to the differential approach plurality of time values corresponding to elementary cells composing an inspected object and lengths of ultrasonic wave in elementary cells;

(iii) calculating values of longitudinal wave velocity and porosity corresponding to a plurality of travel times corresponding to a plurality of elementary cells; and

, (iv) evaluating longitudinal wave velocity in material matrix of said heterogeneous object;

(f) an image formation unit adapted for;

(i) two- and three-dimensional mapping distributions of said longitudinal wave velocity and porosity in an inspected object volume;

(ii) contouring a physiologically distinct area of a human body;

(iii) contouring a defected area of risk of object fracture;

(iv) evaluating longitudinal wave velocity within said defected area; and

(v) estimating a fracture risk value for said inspected object.

14. The ultrasonic system of claim 13, wherein said layered artificial system is formed from a group consisting of at least one plate, pillow, sleeve, and any combination thereof.

15. The ultrasonic system of claim 13 wherein the image formation unit is further capable of mapping said longitudinal wave velocity and bone porosity.

16. The ultrasonic system of claim 13, wherein the image formation unit is further capable of evaluate risk of object fracture and to reveal its location and to detect a defected area dimension and location.

17. The ultrasonic system of claim 13 wherein said transducer'"'"'s grid is adapted to transmit and receive said ultrasonic wave beams at a constant incident and receiving angle selected in the region from 0 up to 90 degrees to a longitudinal axis thereof.

18. The ultrasonic system of claim 13, wherein said object is a bone and the image formation unit further is capable of evaluating risk of bone fracture, revealing its location and detecting a sore area dimension and location thereof.

19. The ultrasonic system of claim 13, in which said heterogeneous object accommodating in layered system is characterized by an acoustic impedance gradient providing the ultrasonic waves propagation through the heterogeneous object along a non-linear paths.

20. The ultrasonic system of claim 13 wherein the ultrasonic waves have a length λ

_wave of the ultrasonic wave used for measurements which is commensurable to a dimension d of the inspected object, such that inspected object in layered system is considered as a local heterogeneity.