ULTRASONIC METHOD AND APPARATUS FOR ASSESSMENT OF BONE

Abstract

Methods and apparatus are disclosed for non-invasive bone evaluation based on a broadband ultrasonic transducer emitting a train of ultrasonic wave packets of multiple carrier frequencies ranging from about 50 kHz to about 2 MHz. Receiving broadband ultrasonic transducer accepts broadband ultrasonic signal propagated through the bone. Computer processing means provide for data analysis and feature extraction allowing diagnostic evaluation of the bone, including comparing the features of the received signal to a database of known bone conditions.

Description

BRIEF DESCRIPTION OF DRAWINGS

The invention is described in detail thereafter in connection with the schematic drawings in which:

FIG. 1 is a flow chart of the method for assessment of bone conditions based on calculating multiple parameters of received waveforms according to the present invention;

FIG. 2 is a flow chart of a method for assessment of bone conditions based on the use of plurality of transfer functions for all transmitted and corresponding received waveforms according to the present invention;

FIG. 3A shows a train of transmitted ultrasound pulses of multiple frequencies and corresponding train of acquired pulses during one cycle of measurement;

FIG. 3B shows a sequence of received signal waveforms at multiple frequencies received on tibial diaphysis of a healthy male;

FIG. 4 shows gradual changes of ultrasonic signal at different frequencies in 0.25-1 MHz band in 4 mm thick double-layered phantoms with increasing porosity;

FIG. 5 shows exemplary emitted and acquired low-frequency (100 kHz) and high-frequency (1.0 MHz) signals in the part of the tibia containing predominantly compact bone;

FIG. 6 shows exemplary emitted and acquired low- and high-frequency signals in the part of the tibia containing predominantly spongeous bone;

FIG. 7 shows the steps in the evaluation of a phase profile by processing the 2-D map of acquired waveforms, where the vertical coordinate is distance along the scanned bone and the horizontal coordinate is time;

FIG. 8 is an illustration of using the wavelet technique to detect the time position of an acquired wave packet at a predetermined site of examined bone so as to capture a reference phase in the acquired waveform to further calculate the time position and profile of said reference phase along the length examined bone;

FIG. 9 is an illustration of an algorithm for detection of the received signal true maximum;

FIG. 10 is an illustration of employing the wavelet technique to detect the time position of the acquired wave packet and a center frequency of this wave packet corresponding to the maximum value of wavelet modulus;

FIG. 11 illustrates the calculation of relative strength for signals related to low- and high-frequency waveforms used for determination of signals strength ratio;

FIG. 12 presents a set of waveform parameters displayed in the form of profile graphs;

FIG. 13 is an illustration of processing of 2-D map of the received waveforms;

FIG. 14 is a chart showing the determination of proximity index;

FIG. 15 is an illustration of assessment of bone osteopenia and osteoporosis using the data on axial profile of acoustic wave fronts;

FIG. 16 is an illustration of assessment of signal fronts for evaluating the stages of bone osteopenia and osteoporosis using the data on signal intensities in predetermined informative areas of the axial profile;

FIG. 17 is a perspective view of one embodiment of the apparatus for assessment of bone according to the present invention;

FIG. 18 is a perspective view of another embodiment of the apparatus for assessment of bone with positioning means made in the form of a tape with a series of holes defining the points of successive measurements along the chosen trajectory over the tested bone;

FIG. 19 is a block-diagram of yet another embodiment of the apparatus for assessment of bone according to the present invention;

FIG. 20 shows two cross-sectional views of a broadband ultrasonic transducer in the accordance to the present invention; and finally

FIG. 21 shows two cross-sectional views of another broadband ultrasonic transducer in accordance to the present invention.

DETAILED DESCRIPTION OF THE INVENTION

A more complete appreciation of the invention and many of the advantages thereof will be accomplished by reference to the following description when considered in connection with the accompanying drawings. Wherever possible, the same reference numerals will be used throughout the drawings and descriptions to refer to the same or like parts.

FIG. 1 is a flow chart illustrating the main steps and procedures of the method for non-invasive and quantitative assessment of bone conditions according to the present invention.

As illustrated in FIG. 1, this method refers to the following specific exemplary steps taken for each evaluated bone:

1) 11—emitting a train of multiple wave packets in the ultrasonic frequency range from 50 kHz to 2 MHz by the emitting transducer;

2) 12—acquiring a train of ultrasonic pulses of multiple frequencies passed along the examined bone by the receiving transducer and recording the received waveforms;

3) 13—repeating the measurement at multiple points along the predetermined trajectory of scanning;

4) 14—calculating the waveform parameters for ultrasonic signals related to different modes of ultrasonic waves identified as longitudinal (bulk), surface, and guided waves;

5) 15—plotting axial profiles of waveform parameters;

6) 16—comparing the axial profiles with statistical reference data; and

7) 17—evaluating the level of abnormality or similarity to known bone pathology.

Steps 11-17 may be carried out in available zones of several types of long bones like the tibia, radius, ulna etc. to provide more comprehensive diagnostics of bone condition.

FIG. 2 is a flow chart illustrating another method for non-invasive and quantitative assessment of bone conditions according to the present invention, which has the same initial steps and procedures coinciding with those of the method illustrated in FIG. 1 but instead of steps 14-17 it adds steps 18 and 19. In the step 18, the signals acquired by the receiving transducer are related to corresponding transmitted signals and the acoustic transfer function for every point of measurement along the scanning trajectory is calculated. The acoustic transfer function contains comprehensive information on plurality of parameters of various modes of acoustic waves propagating in the tested site of the bone. In the step 19, the spatial profile of the acoustic transfer functions is used for assessment of the tested bone by correlating it with corresponding transfer functions from the database developed using systematic statistical data obtained in clinical studies.

FIG. 3A shows an example of a train of transmitted ultrasound pulses of multiple frequencies 209 and corresponding train of acquired pulses 210 during one cycle of measurement. The transmitted train 209 comprises ultrasonic pulses with carrier frequencies of 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.0 MHz. The acquired train of pulses 210 was obtained on tibial diaphysis of a young healthy male.

FIG. 3B shows close-ups of separate ultrasonic pulses 211-220 of the acquired train 210 shown in FIG. 3A. The sequence of received signal waveforms 211-220 at multiple frequencies in the range of 100 kHz-1 MHz was obtained on tibial diaphysis of a healthy male. Different types of acoustic waves become apparent at different frequencies and can be identified as bulk waves 223 (the fastest wave component at high frequencies); surface waves 222 (the slower wave component of greater amplitude at high frequencies) and guided waves 221 (the dominant wave component at low frequencies).

FIG. 4 shows typical changes in the received ultrasonic waveforms 227-230 in bone phantoms with varying cortical thickness and porosity level at different frequencies in 0.25-1 MHz band. Composite two-layer phantoms modeling advanced porosity in deep lying bone layers from which osteoporosis progresses were measured. The pores were mimicked by water soaked poppy seeds embedded in the epoxy. Thicknesses of the layers were 2.5 mm solid/1.5 mm porous (panel 224), 1.5 mm solid/2.5 mm porous (panel 225), 3.5 mm/0.5 mm porous (panel 226). Gradual substitution of the solid layer by the porous one results in redistribution of ultrasound energy between various acoustic wave modes in the received signal. Experimental data clearly shows dramatic changes in the spatio-temporal parameters of signals at every one of the tested frequencies upon variations of cortical thickness and porosity, quantitative evaluation of which is the key issue in osteoporosis assessment.

FIG. 5 shows low-frequency excitation signal 21 and high-frequency excitation signal 24 positioned on the same time scale with corresponding ultrasonic response signals obtained from bone in areas containing predominantly compact tissue. Typical signals in the tibia middle shaft are presented as an example. The following portions of ultrasonic signals passed through the bone are analyzed and used for wave profile calculation:

• • 1) for the low-frequency signal, a portion 22 (or the first wave packet) and a portion 23 characterized by maximum intensity and related to the guided wave,
  • 2) for the high-frequency signal, a portion of 26 related to the bulk wave and the slower portion 27 of greater amplitude related to the surface wave is analyzed. A predetermined time interval 25 preceding the arrival of first signal is used for calculating noise levels as control references for the electronic measurement circuit.

FIG. 6 shows low-frequency excitation signal 21 and high-frequency excitation signal 24 positioned on the same time scale with corresponding ultrasonic response signals obtained from bone in areas containing predominantly spongeous tissue. Typical signals in the tibia epiphisis are presented as an example. At both low and high frequencies, due to a thin cortex, maximum energy of the wave packets corresponds to the bulk wave or longitudinal sub-surface wave propagating in the spongeous bone. The corresponding portions of low-frequency and high-frequency signals 31 and 32, respectively, related to the bulk wave are used in this case for parameters calculation.

FIG. 7 illustrates the steps in composing the waveform parameter profile by processing a 2-D map of acquired ultrasonic signals 41 obtained by stepwise scanning the probe along a bone. At first, the 2-D signal map 42 is constructed from initial ultrasonic signals 41, where the vertical coordinate is distance along the scanned bone and the horizontal coordinate is time. Signal amplitudes are depicted by the gray scale of the image. Every horizontal line in 2-D signal map 42 corresponds to a single ultrasonic signal recorded at a specified distance along the bone from a starting reference point. Further, 2-D signal map 42 is processed by smoothing using rectangular and Gaussian filters to eliminate random deviations and to obtain smoothed 2-D signals map 43. In the next step, the image binarization is performed and 2-D map 44 is produced, where signal intensities below or above a predefined threshold allow separating a narrow border line 46 characterizing the wave profile. Additional filtering applied to map 44 provides smoothed binary map 45 and wave phase profile 47.

FIG. 8 is an illustration of using the wavelet technique to detect times-of-flight for the low-frequency guided wave. Initial signal 51 is processed by complex wavelet functions having waveforms similar to the excitation signal, and the maximum of calculated wavelet modulus 52 is detected. Estimating the temporal position of this maximum allows calculating the group velocity of the wave packet. This algorithm is also used to choose certain reference phase and track this reference phase for the same wave packet in every successive measurement point along the scanning trajectory. The temporal position of the zero intersection point 53 on the negative slope of the signal 51, nearest to the position of the wavelet modulus maximum 52, is calculated. The signal 54 recorded at the next measurement point along the scanning trajectory is analyzed and zero-intersection point 55 is detected by the same method as the point 53. A specified time interval 56 around the point 53 provides a range to search for the zero-intersection point 55. If point 55 is located outside the interval 56, it is rejected and the next measurement point along the bone is analyzed relative to the interval 56. If point 55 is located inside the interval 56, it is considered correct. This constraint helps to eliminate erroneous measurement data and makes the wave phase axial profile calculation more robust.

FIG. 9 is an illustration of signal strength calculations used as an algorithm for detection of a true signal maximum when the wave maximums fluctuate among signals in adjacent points along the bone due to measurement errors and noise of various origins. Signals 61, 62, and 63 are ultrasonic signals at adjacent measurement points along the bone with zero-crossing points A, B, C related to adjacent wave maxima. To identify zero-intersection point in each ultrasonic signal, the algorithm presented above in the description of FIG. 8 is used. The signal strength calculation algorithm includes the steps of determination of amplitude span around the zero-interception points A, B, and C. Nearest maximum and minimum around the corresponding point on each period of the waveform are detected, such as maximum 64 and minimum 65 around point C1. The sum of the amplitude spans for A1-A3, B1-B3, and C1-C3 are calculated and the largest sum represents the maximum signal strength. This procedure is applied to all measurement points along the bone.

FIG. 10 illustrates the use of a wavelet technique for detecting both the temporal position and the center frequency of acquired wave packet passed through a bone. Panel 72 shows a 3-D surface 74 representing a wavelet modulus function along time coordinate calculated for the frequencies ranging from about 100 kHz to about 200 kHz. The vertical axis is the normalized amplitude of wavelet modulus applied to the analyzed wave packet. The projection of the maximum amplitude 75 determines the temporal position 76 and the center frequency 77 for the wave packet. A frequency shift is defined as a difference between the center frequency of the acquired wave packet and the center frequency of the emitted pulse.

FIG. 11 illustrates the calculation of relative strength for signals related to low- and high-frequency waveforms used for determination of signals strength ratio. Low-frequency and high-frequency signals displayed in panels 111 and 112 accordingly, were recorded simultaneously during a single measurement. Thus, influence of variations of probe contact pressure and contact condition on relative power ratios is eliminated. At low frequency, signal strength W1 is calculated as the peak amplitude 113 of the wavelet envelope 114 corresponding to the guided wave in the compact bone or bulk wave in the spongy bone. At high frequency, signal strengths W2 and W3 are determined as peak amplitudes 116 and 117 of the wavelet envelopes 115 and 118 related accordingly for the bulk and surface waves.

FIG. 12 represents a set of calculated axial profiles of the waveform parameters along a bone including:

• • 1) 125—a profile of frequency shift for the high-frequency surface wave,
  • 2) 126—a profile of frequency shift for the high-frequency longitudinal wave,
  • 3) 127—a profile of frequency shift for the low-frequency guided wave,
  • 4) 128—a profile of amplitudes ratio W2/W1 (see FIG. 8), and
  • 5) 129—a profile of amplitudes ratio W3/W1 (see FIG. 9).

These profiles are also used for evaluation of the measurement data quality. It is assumed that a profile graph along a bone is a smooth line. Possible deviations of the measured parameters can be partly caused by signal noise, measurement error, poor contact of the ultrasound transducer with a skin surface or by improper probe handling. The standard deviation of the calculated parameters from the smooth curve is calculated for selected profile graphs and the obtained value serves as a “quality index” for the entire bone examination procedure. The quality index should be within certain predetermined limits to qualify the examination data as satisfactory. This index may also serve as a feedback for operator training.

FIG. 13 illustrates steps of processing of a 2-D map of the received waveforms. Initial map 131 is composed of received ultrasonic signals. The vertical coordinate in the map 131 is the distance along the scanned trajectory, the horizontal coordinate is time, and the signal amplitude is depicted by the gray scale of the image. Map 132 is obtained by smoothing and interpolation of map 131. The interpolation procedure eliminates sharp deviations in the pattern. Smoothing is made by spatial filtering of the pattern and normalization of the amplitude in every signal line. Signal rectification for every line of the map 132 results in the map 133 which is then used for selection of the region of interest 134. This region of interest 134 is defined as an area of the map corresponding to the time interval with the most pronounced changes of the pattern. The pattern 135 is obtained by enveloping procedure or low pass filtering of the pattern in the region of interest 134. The pattern 136 is obtained by binarization of the pattern 135, where faster front 137 and slower front 138 of enveloped ultrasonic signals are determined by selecting two amplitude thresholds. The pattern 139 clearly showing the signal front profile 140 is obtained from the pattern 136 by adjusting the brightness and contrast.

FIG. 14 is a chart illustrating the determination of proximity indexes for assessment of bone conditions. Waveform maps 141 and 142, which are similar to map 132 obtained by the procedure illustrated in FIG. 10, are statistically mean waveform maps stored in the database on classified bone conditions. For example, the map 141 corresponds to tibia of young, healthy individuals, and the map 142 corresponds to tibia for an osteoporotic group of subjects at elderly age. Waveform map 143 obtained on the tibia of a patient is compared with maps 141 and 142 using images matching algorithms and then certain similarity factors are calculated. Proximity index is defined as the ratio of these factors, which show compliance of the waveform map 143 to the maps 141 or 142.

FIG. 15 is an illustration of the use of received wave packet fronts for evaluating the stages of bone osteopenia and osteoporosis as progressing processes of bone degradation when it is developing along the bone length. Derivation of a wavefront profile (pattern 140) was illustrated above in FIG. 10. Wavefront profiles 151, 152 and 153 correspond to norm, osteopenia and osteoporosis, respectively. Zones 154, 155 and 156 are diagnostically informative zones of the wavefront profile. Signal arrival times 157, 158 and 159 shown by white arrows are quantitative characteristics of bone condition and pathologic processes.

FIG. 16 is an illustration of the use of signal strength for evaluating of stages of bone osteopenia and osteoporosis as progressing processes of bone degradation when it is developing along the bone length. Patterns 161, 162 and 163 correspond to norm, osteopenia and osteoporosis, respectively. Areas 164 and 165 are predetermined as diagnostically informative zones of the pattern. Average or integral intensities of areas 164 and 165 are quantitative characteristics of bone condition and pathologic processes.

FIG. 17 is a perspective view of an embodiment of the apparatus for quantitative and non-invasive assessment of bone in the accordance with one aspect of the present invention. The apparatus comprises a hand-held ultrasonic scanning probe 171, an electronic unit 172 for excitation and acquisition of ultrasonic and position signals and a computer 173 as a data processor and display. An operator places the probe 171 on a skin surface above a patient's bone at a predefined starting anatomical reference point, then presses the probe against the bone to obtain the first record, and moves the probe step-by-step, repeating the measurements along the specified trajectory.

FIG. 18 is a perspective view of another embodiment of the apparatus where the probe positioning means is realized in the form of a tape with a series of markings (such as holes) defining the points of successive measurements along the chosen trajectory over the tested bone.

FIG. 19 is a block-diagram of the apparatus for quantitative and non-invasive assessment of bone in the accordance with yet another aspect of the present invention. The ultrasonic scanning probe 171a being positioned on a skin surface above a bone 180 emits a series of broadband ultrasonic pulses through a thin layer of soft tissues by means of an ultrasonic transmitting transducer 182. An ultrasonic transducer 184 receives the ultrasonic signal passed through a bone 180. The acquired signal via the amplifier 187 is transmitted to an electronic unit 188. The probe 171a includes a distance metering means 185 realized in the form of a roller-based motion tracker. The motion tracker obviates the need for a tape as described above and allows monitoring and selecting proper successive locations for measurements based on a distance from a starting point. In a further improvement of this version of the device, a non-contact optical motion tracker similar to that used in the computer optical mouse can be employed.

An optional light emitting diode 186 is used to provide feedback to the operator about the amplitude of received ultrasonic signals to control the probe's applied force. The electronic unit 188 provides initial series of wave packets to transmitting transducer, acquires and processes the signals from the receiving transducer, acquires the data from the positioning system and communicates in real-time with the computer 189.

FIG. 20 shows two orthogonal cross-sections of a broadband ultrasonic transducer in accordance with the present invention. Frontal 190 and lateral 191 cross-sectional views are represented. The preferred design of a transducer consists of:

• • 1) a miniature piezoceramic transducer 192 shaped as a prism with inclined facets pointing towards the bone and mounted on a metal support rod 193, which is a hot electrode connected to the output of the electronic unit;
  • 2) a compound fill up 194 (for example, an epoxy compound);
  • 3) a metal shield housing 195;
  • 4) and an electrical insulation 196, for example, a Teflon cylinder.

A silver wire 199 connects the piezoceramic transducer 192 with the shield housing 195, which is used as the ground electrode. A skin contact surface 197 is flat with rounded edges to avoid local overpressure on the patient's body with a linear ridge 198 providing reliable acoustic contact with the tissue.

FIG. 21 shows two orthogonal cross-sections of another embodiment of a broadband ultrasonic transducer in accordance with the present invention. Frontal 200 and lateral 201 cross-sectional views are represented. The transducer consists of a miniature ellipsoidal prism-shaped piezoceramic element 202 with inclined facets pointing towards the bone. The piezoceramic element 202 is mounted on a metal support rod 203, which is a hot electrode connected to the output of the electronic unit. Inner transducer volume is filled up with a compound 204. A metal shield housing 205 provides mechanical protection of the transducer and is used as a ground electrode. A thin metal wire 206 connects the transducer 202 with the shield housing 205.

Although the invention herein has been described with respect to particular embodiments, it is understood that these embodiments are merely illustrative of the principles and applications of the present invention. It is therefore to be understood that numerous modifications may be made to the illustrative embodiments and that other arrangements may be devised without departing from the spirit and scope of the present invention as defined by the appended claims.

Claims

1. A method for quantitative non-invasive assessment of bone conditions comprising the steps of: (a) acoustically coupling of an ultrasonic probe with a body surface over the bone, said probe comprising a broadband emitting transducer and a broadband receiving transducer,(b) emitting a train of ultrasonic wave packets of multiple carrier frequencies and propagating said wave packets in axial propagation mode along a predetermined trajectory on said bone with said broadband emitting transducer,(c) acquiring a train of ultrasonic signals from said broadband receiving transducer,(d) repeating the steps (a), (b) and (c) at various locations along said trajectory,(e) defining said probe positions along said trajectory at the same time as acquiring said series of ultrasonic signals,(f) calculating spatial profiles of selected parameters of the received waveforms for at least a portion of said bone,(g) evaluating diagnostically relevant features of said spatial profiles, and(h) evaluating bone condition from said diagnostically relevant features of said spatial profiles.

2. The method as in claim 1, wherein the carrier frequencies for said ultrasonic wave packets range from about 50 kHz to about 2 MHz.

3. The method as in claim 1, wherein said trajectory covers spongeous, transition, and compact parts of said bone.

4. The method as in claim 1 including assessment of tibia of a human.

5. The method as in claim 1, wherein said step (f) further includes extracting data from a selected region of interest along said trajectory.

6. The method as in claim 1, wherein said step (h) includes relating said diagnostically relevant features of said spatial profiles to a database of statistical wave profiles for various bone conditions.

7. The method as in claim 1, wherein said step (f) further including analyzing said acquired ultrasonic signals to determine portions thereof related to different modes of ultrasonic waves, said different modes defined as longitudinal, surface, and guided waves.

8. The method as in claim 7, wherein said step (f) further including calculating waveform parameter profiles for said determined portions related to different modes of ultrasonic waves.

9. The method as in claim 8, wherein said waveform parameter profile is composed by processing a 2-D map of said acquired ultrasonic signals.

10. The method as in claim 1, wherein said step (f) further including using a wavelet technique for detecting both a temporal position and a center frequency of said acquired ultrasonic signals.

11. The method as in claim 8, wherein said step (f) further including calculating of a quality index defined as deviation of calculated parameters from a predetermined smooth curve.

12. The method as in claim 11, wherein said quality index is provided as a feedback for an operator, said quality index indicating a degree of operator error in handling said probe during examination of said bone.

13. The method as in claim 9, wherein said step (h) further including calculating proximity indexes of said 2-D map of said acquired train of ultrasonic signals to predetermined statistically mean 2-D maps of ultrasonic signals indicating various bone conditions.

14. The method as in claim 13, wherein said predetermined maps include maps for assessment of various stages of osteopenia and osteoporosis.

15. A method for quantitative non-invasive assessment of bone conditions comprising the steps of: (a) acoustically coupling of an ultrasonic probe with a body surface over the bone, said probe comprising a broadband emitting transducer and a broadband receiving transducer,(b) emitting a train of ultrasonic wave packets of multiple carrier frequencies and propagating said wave packets in axial propagation mode along a predetermined trajectory on said bone with said emitting transducer,(c) acquiring a train of ultrasonic signals from said receiving transducer,(d) repeating the steps (a), (b) and (c) at various locations along said trajectory,(e) defining said probe positions along said trajectory at the same time as acquiring said series of ultrasonic signals,(f) calculating an acoustic transfer function for every point of measurement along said trajectory by relating said acquired train of ultrasonic signals from the receiving transducer to the emitted train of ultrasonic signals from the transmitting transducer,(g) calculating spatial profiles of said acoustic transfer functions for at least a portion of said bone, and(h) evaluating said bone condition by correlating said spatial profiles with corresponding spatial profiles for various bone conditions obtained from a clinical database of spatial profiles.

16. An apparatus for quantitative non-invasive assessment of bone conditions comprising a hand-held probe including a broadband emitting transducer and a broadband receiving transducer, an electronic means for excitation and acquisition of ultrasonic signals and a data processing and display means, said broadband emitting and receiving transducers adapted to respectively emit and receive a train of ultrasonic wave packets of multiple carrier frequencies.

17. The apparatus as in claim 16, wherein the carrier frequencies for said ultrasonic wave packets range from about 50 kHz to about 2 MHz.

18. The apparatus as in claim 16, wherein said probe further including a positioning means adapted to relate measured ultrasonic data to position of said probe on a predetermined trajectory of bone evaluation.

19. The apparatus as in claim 18, wherein said positioning means is a tape with markings indicating evaluation locations for said bone assessment.

20. The apparatus as in claim 18, wherein said positioning means including a distance measuring means.