A METHOD AND AN APPARATUS FOR INVESTIGATING MATERIAL PROPERTIES OF BONE USING ULTRASOUND

Abstract: A method for investigating the mechanical properties of bone inside a live animal or human being, said method comprises launching an ultrasound pulse wave into the body of said being and establishing a trace related to the magnitude of a reflected ultrasound wave versus the time lapsed since the launching of the pulse wave. On the trace an interval in which the magnitude of the reflected ultrasound wave exhibits a steady decline versus time, is identified, and the attenuation of the ultrasound wave based on the readings of the trace within this interval is computed.

Description: Title of the invention A method and an apparatus for investigating material properties of bone using ultrasound. Background of the invention 1. Field of the invention The present invention relates to a method and an apparatus for investigating the mechanical properties of tissue inside a live animal or human being. In particular, the present invention relates to the use of ultrasound for assessment of the mechanical properties of tissue. The invention more particularly relates to methods of investigating the mechanical properties of a material which may not be easily accessible and which may not be adequately described by simple models. The invention further relates to a method of investigating the mechanical properties of bone in vivo. The investigation of the mechanical properties of bone in vivo is of great interest in view of the occurrence of osteoporosis. Osteoporosis may be defined as a long term metabolic deficiency causing an unbalance in the natural process of bone resorbtion and bone rebuilding with the result of a loss of mechanical strength and increased risk of fractures. The degradation of the mechanical strength of the bones may proceed to a stage where even minimal trauma result in bone fractures. It has been estimated that osteoporosis affects from 10-20% or maybe even more of the female population and about 10% of the male population above the age of 50. Treatments do exist which may delay or reverse the progression of osteoporosis. However, effective methods permitting an accurate diagnosis to be established at an early stage and permitting accurate monitoring of patient response to medication would be of great value. Nevertheless, the accurate assessment of osteoporosis is difficult. The bones in the skeleton are by nature non homogeneous and different parts of the skeleton may not be affected to the same degree. The material strength of the bones naturally changes over time, reaching a maximum about the age from 20-30 years and gradually declining later on. Individual differences may be substantial. Furthermore, the mechanisms involved and the factors controlling this process may not be fully understood. The prior art US reissue patent No. Re32782 discloses a method for determining in vivo strength of bone in a live being. According to this method, a first transducer launches an acoustic pulse though a bone and surrounding soft tissue which pulse is received by a second transducer. The distance between the transducers is measured and so is the transit time of the pulse between the transducers. These measurements form the basis for determining the effective velocity of the pulse through the bone and surrounding soft tissue which velocity provides one piece of information used for evaluating the strength of the bone . The measurements of the distance and of the transit time in respect of an acoustic pulse launched through the bone with surrounding soft tissue provides no positive information about the path followed by the acoustic pulse. This may disturb the conclusions as the path may be non-linear due to refraction effects and as the acoustic transducers may exhibit strong sensitivity variations on varying orientations, both of which factors introduce potential sources of errors in case of accidental misalignments between the transducers. As the pulse is likely to have travelled through sections of soft tissue and through sections of bone tissue, any effect due to variations in bone structure may be buried in variations caused by other factors. US patent No. 5 038 787 discloses a method for analyzing material properties using reflected ultrasound. According to this method, a bone is placed in a separating medium which includes water and soft tissue and then subjected to ultrasound pulse waves from varying angles of incidence while the reflected wave is received and analyzed. The varying angles of incidence are obtained by shifting the transmitting transducer around. Critical angles of reflection are established and used in the calculation of a matrix of mechanical parameters such as elasticity. The prior art also comprises a suggestion for providing a pair of transducers coaxially in a water tank, submerging a bone sample between the transducers launching ultrasound pulses in various frequencies ranging from 0.2 .MHz to 0.8 MHz and measuring the attenuation as a function of frequency. In this method, any loss of signal will for the sake of the calculation be presumed as due to attenuation. This may amount to a misinterpretation, e.g. in case the loss was due to reflection, likely to arise whenever the pulse wave crosses some interface, e.g. water to soft tissue. In addition, this method does not separate contributions caused by effects in the soft tissue from contributions caused by effects in the bone tissue. US patent No. 4 408 492 which relates to ultrasonic echoscopy, refers in this context to the principle of time gain compensation understood as the concept of increasing the gain of the receiver in proportion to the echo time delay in order to compensate for the progressive attenuation sustained by echoes from deeper reflecting surfaces within the object. The publication mentions the definition of a first and a second contour within the object subjected to the examination, determining the attenuation exhibited between the first and second contours and using the time gain compensation facility as necessary to compensate for the attenuation as determined. The purpose of this method is to overcome the problem of shadowing in a situation where a local highly absorbent area obscures deeper lying information. US patent No. 4 414 850 pertains to a measurement method and a system for measuring characteristics of attenuation of domains in an object, by which method ultrasonic waves are transmitted into the object and reflected ultrasonic waves are received and analyzed. Attenuation of the reflected waves is presumed to be a linear function of frequency. Attenuation coefficients are determined for various domains within the body. US patent No. 4 941 474 relates to a multi-variable analysis of bones for the purpose of detecting abnormal bone conditions. One method explained comprises launching into the patient an ultrasonic signal having components in the spectrum from about 100 kHz to about 600 kHz, receiving reflected signals, and evaluating the magnitude of the received signals. Summary of the invention The invention in one aspect provides a method as recited in claim 1. This method produces information useful for characterizing bone inside a live animal or human being which information may provide a valuable part of the background for diagnosing e.g. osteoporosis. The method according to the invention permits retrieval of information related to trabecular bone, singling out these data from signals related to other kinds of tissue. The method according to the invention makes it feasible to use a high range of frequencies, by which it is possible to focus investigation on minute zones of interest . Trabecular bone, sometimes also referred to as cancellous bone forms only part of the bone structure. Trabecular bones may be explained as a porous structure comprising a network of mineralized collagen fibrils and marrow (fat) . Trabecular bones are generally sheathed within or enclosed by a layer of cortical bone which comprises a more dense structure. In long bones, the central part of the bone has a tubular shape mainly comprised of a wall of cortical bone. Towards the ends of the long bones, the cortical walls become thinner and increase in diameter as they form the bone end portions where most of the volume enclosed by the cortical walls is filled with trabecular bone consisting of plates, spicules or small septums of cancellous bone. The pattern of age related bone loss differs between cancellous and cortical bone. Commonly, trabecular bone decreases in number rather than thickness. Cortical bone usually becomes more porous and the cortical walls also become thinner. These changes in structure give rise to corresponding changes in acoustical properties, the increasing porosity of cancellous bone giving rise to a decrease in acoustic attenuation. Some data of acoustic properties of different kinds of tissue are listed in the table below. Table of physical properties of tissue Density Speed of Attenuation Sound g/cm 3 m/μs dB/cm Cortical 2.00 3.50 6.90 tissue Trabecular 1.06 1.89 1.9 to 15.7 tissue Fat 0.93 1.48 1.00 In the progress of osteoporosis, structural bone tissue is replaced with fat which exhibits a lower value of acoustical damping. Although the changes take place in cortical bone as well as in trabecular bone, trabecular bone due to its occurrence in more bulky regions of a homogeneous structure, seems to offer the prospect of a more accurate characterization by acoustic methods and thus the chance of detecting any change at an early stage. Acoustic characterization of trabecular bone has, however, been difficult in methods according to the prior art since the cortical bone reflects most of the acoustic energy, making it difficult to retrieve any useful data from the layer below. However, these problems are overcome by the method according to the invention. According to this method, a trace of the logaritm of the magnitude of an ultrasound pulse, plotted versus ime of flight is relied on for retrieving information about a bone tissue inside a live animal or human being. The method is non- invasive and does not require direct access to the tissue under investigation. Providing a plot of the logaritm of the amplitude permits overviewing a wide dynamic range within a single plot and permits direct reading of attenutation factors. As known in the art, the traces are likely to exhibit sections of steady decay but possibly also peaks or sections of level off-sets. Such sections of peaks or level off-sets are likely to arise in cases where the ultrasound pulse traverses an interface between regions of different acoustic impedance. Sections of steady decay are likely to be related to regions of tissue which are homogeneous on a macroscopic level, although likely to be non-homogeneous on a microscopic level. The microscopic inhomogenities will give rise to scattering and reflection of the ultrasound pulse. The transmitted pulse as well as any reflections are subjected to attenuation on passing through the tissue. In sections of steady decay of the trace, it may be presumed that the corresponding regions of tissue comprise distributed but mutually similar, microscopic non-homogenities in order that the ratio of reflected energy to transmitted energy is the same in any sub- region thereof. Based on this presumption, the attenuation in the corresponding region of tissue may be computed by comparing readings along the trace which in effect amounts to comparing levels of reflections from different minute non-homogenities to each other and relating them to the distance from the transducer. In an object comprising a boundary with a rather steep change of acoustic impedance, this transition will give rise to reflection of a substantial proportion of the acoustic energy. As the acoustic transducer itself in transducers according to the state of the art generally provides a reflecting surface, any pulse reflected into the transducer may be reflected from the transducer to re- enter the test object where it may be subjected to a renewed reflection. The result is that one boundary in the object may manifest itself by a trace which includes several spikes, effectively repeating or mirroring the image of the first spike and confusing the display. The step of identifying an interval of the trace exhibiting a steady decline is important with the view of avoiding falsifications of the data due to such mirroring. According to a preferred embodiment, the trace of the reflected wave magnitude is established by detecting peak levels of a high frequency signal received from the transducer. Evaluation of the peak detected signals permits a simplification of the data processing. According to preferred embodiments, the launching of the ultrasound pulse wave is carried out by placing an ultrasound transducer in contact with the skin. The transducer coupling to the skin may be enhanced by applying a layer of gel to the skin prior to the application of the transducer. This provides a simple procedure which combines accuracy in the measurement with convenience to the person who is the subject for the analysis . According to the preferred embodiments, the same ultrasound transducer is used for launching the ultrasound pulse wave and for monitoring the reflected ultrasound wave. The use of a single transducer simplifies the handling and avoids potential errors related to malalignment or mismatching. According to the preferred embodiments, the ultrasound pulse wave is launched by using a transducer which is adapted for delaying the transmitted pulse wave as well as the received pulse wave. The delay feature spaces the advent of multiple reflections on the time scale from the advent of the initial reflection, usually preferred for purposes of interpreting the results. The delay feature is important in applications where the task object comprises a transition zone of strong reflection close to the transducer surface and where the information wanted relates to structures beyond such transition zone. The delay should preferably by adapted to space any mirror images of pronounced spikes away from the region to be analyzed. This simplifies data interpretation and enhances accuracy in the measurements . According to further preferred embodiments, the transducer may be adapted to produce focused waves in order that the zone of investigation is narrowed down for enhanced resolution. A focused transducer is adapted to emit a focused beam of ultrasonic waves. Due to the reciprocal nature of the transducer, the sensitivity may be concentrated in a similar area. Focusing the transducer narrows down the zone of sensitivity and suppresses sensitivity to signals from other directions. Focusing an ultrasonic transducer may e.g. be achieved by adapting the transducer front surface to be curved concavely towards the focal point or by designing the active surface of the transducer in a pattern of planar, annular, concentric rings, which are controlled electrically so as to be fired sequentially in order to provide a focusing effect. Due to well known geometric aspects the width of the active area covered by the beam is linked to the spacing from the transducer, i.e. the effectiveness in the focusing is dependent on the distance. The focused transducer may to advantage be combined with a delay line, provided the focal length is adapted to compensate for the delay line. The delay line spaces the transducer from the test object and makes it possible to achieve a proper focusing by means of a transducer which is adapted to operate in more like a far field condition rather than in a near field condition. A focused transducer used in a far field situation provides the advantage of the effective beam width being uniform over a greater range of distances. According to the preferred embodiments, the ultrasound wave used comprises a pulse enveloping a signal oscillating at a frequency in the range between 0.5 and 5.0 MHz. The lower boundary of this interval is restricted in view of considerations related to resolution and of the lack of focusing possible at lower frequencies. The upper boundary in this interval is restricted in view of considerations concerning attenuation in relation to the required range. In another aspect, the invention provides a method as recited in claim 11. By this method, the relation between attenuation and frequency is established in respect of the investigated structure. This together with the data of attenuarion at the respective frequencies provides a comprehensive characterization of the tissue under investigation, widening the base of data available for drawing any conclusions. According to a preferred embodiment the attenuation vs frequency is presumed to be proportionate to a function of the type f b , i.e. f lifted to the power of b, where the base f represents the frequency, the exponent b is a constant related to a physical property of the test object, and wherein values obtained by measurements at a first frequency and at a second frequency are used to calculate the exponent b. The inventor has discovered that the expression f D seems to provide a more useful mathematical approximation of the attenuation vs frequency than the presumption whcih tends to prevail in the prior art, namely that attenuation is simply proportionate to frequency. The inventor has discovered that the exponent b seems to provide a most useful piece of information for estimating the bone status such as a state of osteoporosis. The invention makes it possible to obtain a reliable estimate of the exponent b and thereby offers the prospect of an early indication of a situation of importance. The invention, in a still further aspect, provides an apparatus according to claim 14. This provides an apparatus for implementing the method according to the preceding claims. Further features and advantages of the invention will appear from the following detailed description of preferred embodiments which is given with reference to the drawings . Brief description of the drawings Fig. 1 is a cross section through an object contacted by an acoustic transducer, the section including the axis of the transducer, Fig. 2 is a block diagram illustrating basic components of the instrumentation used, Fig. 3 is a trace produced on the display of the instrumentation, using a single frequency of excitation, Fig. 4 is a double trace produced on the display of the instrumentation in respect of two different frequencies of excitation, Fig. 5 is a symbolic figure illustrating in its upper part a cross section through an object contacted by a strongly focused acoustic transducer, the section including the axis of the transducer, while the figure in its lower part illustrates a time plot of the reflected signal received by the same transducer, and Fig. 6 illustrates a figure similar to Fig. 5 but modified by including a delay line in contact with a mildly focused transducer. Description of the preferred embodiments All figures are schematic and not to scale and illustrate only the details necessary for enabling those skilled in the art to practice the invention while other details have been omitted for the sake of clarity. In all figures, the same references are used to designate identical or corresponding items. Figure 1 illustrates in section a transducer Tr coupled by means of an intermediate layer of ultrasonic gel G to an object, e.g. a human extremity (calcaneum, patella, ulna, phalanx or other) . For the purposes of this explanation, it is sufficient to say that the object 0, explained as occurring in the direction from the transducer, comprises a layer of skin tissue S, a layer of cortical bone C, a layer of trabecular bone T, another layer of cortical bone C and another layer of skin tissue S. The transducer emits an acoustic pulse beam generally focused to propagate with a maximum of intensity along the axis A. The transducer may be focused and it may include a delay line feature as known in the art. Generally speaking, the acoustic pulse beam will be smeared out transversely in some region centered along the axis A. It lies within the capabilities of those skilled in the art to estimate the width of the beam at respective depths of penetration. Figure 1 indicates two locations, Xi and x 2 , spaced apart and located on the beam axis A. These points will be referred to later in the explanation. Figure 2 illustrates a block diagram of the instrumentation used. The instrumentation basically comprises a transducer, a transmitter, a receiver, a storage means, a display and a controller. The transducer serves the purpose of transmission as well as the purpose of reception of the acoustical wave. In the preferred embodiment, the transducer is specified as being resonant at 2 MHz, with a bandwidth of 30 to 40 % and provided with a damping serving the purpose of enhancing impulse response. The transducer is connected to a transmitter which provides the drive signal for causing the transducer to launch an ultrasonic pulse. The transmitter is adapted for outputting short pulse bursts, e.g. a pulse enveloping a few cycles of oscillation of the transducer. The transmitter is triggered by a signal emitted from the controller. The transducer is also connected to a receiver which processes the signal received by the transducer in order to provide on its output a rectified peak signal which may be averaged over time in order to suppress random error signals. The receiver also receives a synchronization signal from the controller enabling it to suppress the powerful signal caused by the transmitter drive pulse. A display is connected to receive the receiver output signal and adapted to present a trace of peak detected amplitude versus time. The output signal is also fed to a storage means where the signal is recorded for later processing. Reference is now made to figure 3 for a depiction of the trace on the display and for an explanation of the processing of the data. In figure 3, the x-axis is calibrated in microseconds. Although actually representing time, the x-axis is for many purposes interpreted as a scale of distance from the transducer. In view of the different acoustic velocities of the tissues involved, this may not be strictly true, however, it may for many purposes represent an acceptable approximation . The y-axis used for the amplitude is calibrated in decibels, signifying that it represents the logarithm of the relative amplitude. The trace in figure 3 shows an initial steep rise and then a more gradual decay superimposed with a number of localized peaks or deflections. The first peak L is regarded to signify the spill-over effect from the launching of the pulse wave, a second peak S/C is attributed to a reflection at the skin/cortical bone interface, a third peak C/T is attributed to a reflection at the cortical bone/trabecular bone interface and a fourth peak T/C is attributed to be caused by the trabecular bone/cortical bone interface. After the fourth peak the trace settles around a fairly constant level, which is presumed to signify a back-ground of noise. Between the third and the fourth peak, it is possible to identify an interval of steady decline in the trace. The operator in charge of carrying out the measurement will evaluate the trace and will refer to his knowledge concerning the anatomy of the site investigated in order to correlate various aspects of the trace with the topology and various properties concerning the site investigated. Thus the person in charge of the analysis will use his skill to verify that the interval of steady decline indeed relates to a section of trabecular bone, and that the trace within this interval represents valid data, not biased by a background of noise and not affected by reflections related to other causes. Two points ti and t 2 within this interval are selected. From the trace, corresponding amplitude values Ui and U 2 are read. The pulse peak amplitude as a function of distance within a homogeneous region is presumed to obey the equation U (x) = U (xo) * exp (-α * x) where U is the amplitude, X is the distance and α is the spatial attenuation factor. Writing into this formula measurement data from the readings at two points Xi and x 2 corresponding to the time instance of ti and t 2 , it may be shown that: \nU( X] ) - \nU(x 2 ) l ~ X 2 The distance Xi - x 2 may be computed from the time delay multiplied with the speed of sound of the respective tissue . Thus from the example in the figure 3 where ti = 32 μs, t 2 = 59 μs, Ui = 32 dB and U 2 = 4 dB, we have α = -(32- )dB / (32μs - 59μs) = 1.037 dB/μs. Assuming a speed of sound of 0.189 cm/μs, this converts into α = (1.037 dB/μs) / (0.189 dB/μs) = 5.49 dB/cm Reference is now made to figure 4 for the explanation of a procedure involving measurements taken at two different frequencies. Figure 4 illustrates a trace fl identical to the trace of figure 1 and recorded using a transducer pulse frequency of 5 MHz. In addition, figure 4 includes a second trace f2 recorded under identical circumstances except that the recording has been taken using a transducer pulse frequency of 1 MHz. The trace produced at 1 MHz exhibits a shape basically similar to that of the first trace but with a more flat slope attributed to the lower attenuation of the tissue as measured at the lower frequency. The frequency dependency of the attenuation may be accounted for in the above formula by taking α to be instead of a constant value, a function of the frequency as α(f) where a and b are two parameters specific to the object whereas f is the frequency. Once two values of alpha, i.e. cti and α, have been established in respect of two different frequencies, f and f 2 , the physical parameters a and b may be computed according to the following: α- from which b = ln(αα/α 2 ) / ln(fι/f 2 ) The second last of these equations is equivalent to b = log (αι/0-2) / log(fι/f 2 ) In respect of figure 4 and using for the first frequency (fi = 5MHz) data referred to in relation to figure 3, we may complement these data with the amplitude readings 43 dB and 33 dB in respect of the second frequency (f ; = 1MHz) from which we have α : = -(43-33)dB / (32-59)μs = 0.370 dB/μs α 2 = 0.370 dB / 0.189 cm = 1.95 dB/cm b = log . 5.49/1.95)/ log(5/l) b = 0.4495/0.694 = 0.643 a = 5.49 / 5 0 - 643 = 1.95 dB cm "1 MHz " Thus the value of the exponent b has been calculated to 0.643. Inserting these values into the expression for α as a function of the frequency, we have α(f) = 1.95 * f 0 - 6 " dB cm "1 MHz "0 - 643 Reference is now made to Figs. 5 and 6 for an explanation regarding aspects of combining a transducer with a delay line, Fig. 5 depicting a situation with a sharply forcused transducer in contact with a test object, and Fig. 6 depicting a situation with a mildly focused transducer and with a delay line interposed between the transducer and the test object. Fig. 5 is effectively in two parts. The upper part depicts a cross section through the transducer and through the object contacted by the transducer, the section taken by a plane which includes the transducer axis. The transducer is presumed to launch a pulse, which propagates through a layer of soft tissue into a layer of hard tissue, which includes the area of interest. Waves reflected from any boundaries on the way are picked up by the same transducer and plotted versus the time lapsed. The lower part of Fig. 5 is the time plot depicting the reflected wave received by the transducer. In this plot the scaling on the time axis has been adjusted m order to make the time plot immediately comparable to the section above. Thus the time scale has been multiplied by a factor which reflects the propagation velocity and which takes into account the fact that the ultrasonic wave travels forth and back. The time plot reveals four spikes referred to as al, a2, a3, and a4, respectively. The spikes have decaying peak amplitudes with their spacing from the starting point to the left. The first spike, al, is easily correlated with the interface between the transducer and the soft tissue contacted by the transducer. The second spike, a2, is also easily correlated with the boundary between the soft tissue and the hard tissue. The remaining spikes, a3 and a4, seem to suggest the presence of further boundaries within the specimen, however, they more likely represent repeated reflections, i.e. spikes produced by the launched pulse travelling from the transducer to the most imminent surface of the hard tissue, and back more than once. Reference is now made to Fig. 6 which illustrates a situation silimar to the one illustrated m Fig. 5 except for the facts that a delay line has now been interposed between the transducer front surface and the surface of the soft tissue contact and that the type of transducer is different. The delay line is a body of matter capable of transmitting ultrasound without causing scattering or reflections and adapted for spacing the advent of reflections m time. A preferred delay time is 20-30 μs . The practical embodiments may comprise a body of water enclosed within a container, of a polymere or a slab of metal or ceramics with excellent sound propagation capability. A body of water of a length of 50 mm is capable of delaying the signals by 50 mm/1.89 mm/μs = 26 μs . A preferred embodiment comprises a container holding a body of water and with a front which provides an acoustic window for interfacing with the object to be investigated by the ultrasound method. The acoustic window comprises a stratum or a disk of a matter selected for providing an acoustic impedance which matches the acoustic impedances of the water and of the object. In case of a human body, the acoustic impedance, i.e. sound velocity, multiplied by density, is estimated at 1.53 kg/ ( 2 • s). In case of water, the acoustic impedance is 150 kg/ (m 2 • s) . A disk of a TBX polymere, found to have an acoustic impedance of 1.51 kg/ (m 2 • s), has been found to work well. Generally a window within an acoustic impedance intermediate the values of the acoustic impedances of adjacent layers is transparent to the ultrasonic waves and thus transmits a maximum of energy while causing a minimum of reflections. The time plot in Fig. 6 illustrates spikes bl, b2, b3, and b4. Bl correlates with the interface between the transducer and the delay line front surface. B2 correlates with the interface between the delay line back surface and the front surface of the soft tissue. B3 correlates with the interface between soft tissue and hard tissue, and b4 illustrates a mirror image similar to a3 in Fig. 5. However, Fig. 6 clearly shows how the mirror image b4 has been spaced in time by an offset sufficient to avoid confusing this spike with any signals received from the area of interest. This simplifies interpretation of the plot and provides for enhanced accuracy of any readings. The Figs. 5 and 6 also show computed curves delimiting the zone covered by respective focused transducers in order to show the difference of type. In either case the criterion for selecting the transducer type has been to provide a focused transducer adapted for covering an area of a predetermined beam width within the area of interest. However, due to reasons of wave propagation geometry, the beam width varies with the distance, thus the required width only prevails at a particular distance . In the case of Fig. 5, the transducer is rated at covering a beam width at 4.3 mm within distances ranging from 23 mm to 36 mm. However, the short distance requires a transducer with a rather pronounced focusing feature as the transducer must in this case operate within or adjacent to the near field. This has the consequence that the true effective beam width varies substantially within the area of interest as shown in Fig. 5, depending on the distance from the transducer. In the case of Fig. 6, the transducer is rated for covering a beam width of 11 mm within distances ranging from 43 mm to 103 mm. Due to the longer distance this kind of beam width pattern has been obtained by means of a transducer with a moderate degree of focusing. The result is that the true beam width shows less variation within the area of interest, as will appear from Fig. 6. Thus the sensitivity tends to be more uniform over the area of interest and more easily predictable. Although specific methods have been described above in specific contexts, such methods are not excluded from use in other contexts, from combination in other ways and for being used for other purposes. The preceding description is offered with the sole purpose of illustrating the invention and is not intended to limit the scope thereof which is exclusively defined by the appended claims.

Claims: CLAIMS 1. A method for investigating the mechanical properties of bone inside a live animal or human being, said method comprising - launching an ultrasound pulse wave into the body of said being and into said bone, - monitoring the ultrasound wave reflected from said body and establishing a trace related to the magnitude of the reflected ultrasound wave versus the time lapsed since the launching of the pulse wave, identifying a pair of points on the time scale delimiting an interval in which the logaritm of the magnitude of the reflected ultrasound wave exhibits a steady decay versus time, and which interval corresponds with a section of bone estimated to predominantly comprise trabecular bone, - selecting at least two points on the time scale within the interval defined by the pair of points and reading from the trace corresponding values of reflected wave magnitude, - reading the time delay between the selected points, - computing the distance between those points in the bone which correspond with the selected points, using the sound velocity estimated in respect of wave propagation through trabecular bone, and - computing the attenuation of the ultrasound wave based on the readings corresponding with the selected points. 2. The method according to claim 1, wherein the trace of reflected wave magnitude is established by detecting peak levels of a high-frequency signal received from the transducer. 3. The method according to claim 1, wherein the step of launching of the ultrasound pulse wave is carried out by placing an ultrasound transducer in contact with the skin and operating the transducer to transmit the ultrasound pulse wave. 4. The method according to claim 3, wherein a thin layer of gel is applied to the skin prior to the step of placing the transducer in order to enhance the coupling of acoustic energy from the transducer into the skin. 5. The method according to claim 1, wherein the monitoring of the reflected ultrasound wave is carried out by using the same ultrasound transducer for launching the ultrasound pulse wave and for picking up the reflected ultrasound wave. 6. The method according to claim 1, wherein the step of launching of the ultrasound pulse wave is carried out by using a transducer which is adapted for delaying the launching and the picking up of pulse waves. 7. The method according to claim 1, wherein the step of launching of the ultrasound pulse wave is carried out by using a transducer which is adapted for transmitting a focused beam of ultrasound waves. 8. The method according to claim 1, wherein the ultrasound wave comprises a burst pulse, comprising a signal oscillating at a frequency in the range between 0.5 and 5 MHz. 9. The method according to claim 1, wherein the trace of reflected ultrasound wave magnitude versus time is recorded and stored in a recording medium. 10. The method according to claim 1, wherein the attenuation is computed as an average based on a number of points selected within the interval defined by said pair of points. 11. The method according to claim 1, wherein the steps are carried out using ultrasound waves at a first frequency and at a second frequency and wherein the respective computed values are combined to establish an expression relating the attenuation to frequency. 12. The method according to claim 11, wherein attenuation versus frequency is presumed to be proportional to a function of the type f lifted to the power of b, where the base f represents the frequency and the exponent b a constant, and wherein the computed values are used to calculate the exponent. 13. The method according to any of the preceding claims, wherein the transducer comprises an acoustic window formed by a disk for interfacing between the object to be investigated and the remainder of the transducer, which disk comprises a material selected for providing an acoustic impedance intermediate the respective values of acoustic impedance pertaining to said object and to the remainder of the transducer. 14. An apparatus comprising means for implementing the method according to any of the claims 1 through 12.

Inventors: CHRISTIANSEN KAARE, PETERSEN JOHN FINNICH, VAMMEN KLAUS, CHRISTIANSEN, KAARE, PETERSEN, JOHN, FINNICH, VAMMEN, KLAUS
Applicants: CHRISTIANSEN KAARE, PETERSEN JOHN FINNICH, VAMMEN KLAUS, CHRISTIANSEN, KAARE, PETERSEN, JOHN, FINNICH, VAMMEN, KLAUS
Application year: 1998
Main classification: A61B8
IPC codes: A61B8