Next Article in Journal
Indoor Visible Light Positioning: Overcoming the Practical Limitations of the Quadrant Angular Diversity Aperture Receiver (QADA) by Using the Two-Stage QADA-Plus Receiver
Previous Article in Journal
Application of a Novel Long-Gauge Fiber Bragg Grating Sensor for Corrosion Detection via a Two-level Strategy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ultrasound Radiation Force for the Assessment of Bone Fracture Healing in Children: An In Vivo Pilot Study

1
Department of Radiology, Mayo Clinic College of Medicine & Science, Rochester, MN 55905, USA
2
Department of Physiology and Biomedical Engineering, Mayo Clinic College of Medicine & Science, Rochester, MN 55905, USA
3
Department of Orthopedic Surgery, Mayo Clinic College of Medicine & Science, Rochester, MN 55905, USA
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(4), 955; https://doi.org/10.3390/s19040955
Submission received: 25 January 2019 / Revised: 12 February 2019 / Accepted: 19 February 2019 / Published: 24 February 2019
(This article belongs to the Section Biosensors)

Abstract

:
Vibrational characteristics of bone are directly dependent on its physical properties. In this study, a vibrational method for bone evaluation is introduced. We propose a new type of quantitative vibro-acoustic method based on the acoustic radiation force of ultrasound for bone characterization in persons with fracture. Using this method, we excited the clavicle or ulna by an ultrasound radiation force pulse which induces vibrations in the bone, resulting in an acoustic wave that is measured by a hydrophone placed on the skin. The acoustic signals were used for wave velocity estimation based on a cross-correlation technique. To further separate different vibration characteristics, we adopted a variational mode decomposition technique to decompose the received signal into an ensemble of band-limited intrinsic mode functions, allowing analysis of the acoustic signals by their constitutive components. This prospective study included 15 patients: 12 with clavicle fractures and three with ulna fractures. Contralateral intact bones were used as controls. Statistical analysis demonstrated that fractured bones can be differentiated from intact ones with a detection probability of 80%. Additionally, we introduce a “healing factor” to quantify the bone healing progress which successfully tracked the progress of healing in 80% of the clavicle fractures in the study.

1. Introduction

Speed of fracture healing may vary from one individual to another and among different fracture locations based on anatomic location, fracture severity, patient physiological state, and treatment plan. Immobilization is the choice treatment for a majority of fractures. Determining the optimum duration of immobilization is important, as prolonged cast immobilization has multiple drawbacks, including muscle atrophy, joint stiffness, and skin breakdown. Furthermore, each cast change exposes the patient to the risk of cast saw burns or lacerations [1]. Paradoxically, excessive immobilization can delay the healing process, since bone formation needed for healing the fracture site is stimulated by mechanical stresses through mechanosensors within the osteocytes [2]. Significant debate exists amongst musculoskeletal care providers as to identifying when a fracture is healed sufficiently to discontinue immobilization and resume activities. Currently, the assessment of bone fracture healing is performed by serial radiographic and physical examinations [3]. Studies have shown the poor sensitivity of radiographic measurement in the prediction of healing in a rabbit model of bone fracture [4,5,6].
Conventional ultrasound in humans has been reported for detecting clavicle and forearm fractures with variable sensitivity [7]. A 1998 study [8] used conventional ultrasound to assess fracture healing in a preclinical canine model and found that a hyperechoic ultrasound signal correlated 100% of the time with the presence of hard fracture callus biopsy tissue. The development of quantitative tools for bone healing monitoring may assist physicians to identify the healing endpoint, thus preventing prolongation of the treatment course or incorrect timing of the removal of an external or internal fixation device. Several noninvasive quantitative techniques have been investigated to monitor bone fracture healing [9], including mechanical vibrational analysis [10,11,12] and acoustic emission [13,14]. Quantitative ultrasound (QUS) based on axial transmission has been proposed to assess normal and fractured bone [15,16].
The QUS technique involves placing a set of linearly arranged transducers (transmitters and receivers) along the bone axis. The receivers detect the ultrasound waves that have propagated along the cortical shell, and some of the generated parameters can be estimated for tissue characterization. Two-dimensional simulations on bone-mimicking plates, as well as in vivo studies, have investigated the potential of ultrasound axial transmission to evaluate fracture healing [15,16,17]. A model for fracture healing monitoring using QUS has been validated by experimental studies on ex vivo bovine femur and simulation studies [18]. Additionally, the effect of callus mineral composition has been reported by using simulation and experimental phantom studies [19]. Overall, both conventional ultrasound and QUS have variable results due to the presence of soft tissue around the bone. The impacts of overlying soft tissue on QUS and vibrational analysis in evaluating bone quality have been reported [20,21].
The utility of ultrasound radiation force (URF) for the characterization of soft tissues has been of interest for the past two decades. Vibro-acoustography, in particular, uses URF for imaging and tissue characterization [22,23]. The focused URF of vibro-acoustography induces vibrations in the object of interest. The frequency of those vibrations can be much smaller than that of the primary ultrasound used for URF. Object vibrations produce an acoustic field that is detected by a hydrophone sensitive to the frequency band of the induced vibrations. To construct an image of the object, the URF is scanned across it, usually in a raster fashion. The resulting image provides information on the object’s vibrational and mechanical characteristics. The temporal waveform of the URF depends on the particular application. Examples of URF waveforms include single frequency [22,24], multi-frequency [25], tone burst [26,27], and chirp signal [28]. Depending on the application, different types of transducers have been used to produce the URF, such as dual-element confocal [29], linear array [30], multi-row array [26], and sector array [31]. Vibro-acoustography has been used for imaging tissues and implants [32,33,34,35,36]. An overview of the application of this technique on various soft tissues has been previously reported [27]. Additionally, vibro-acoustography may be used for bone evaluation, as described in a recent study by Denis and colleagues [37]. The initial application of vibro-acoustography in the assessment of bone fracture and bone healing was presented in 2006 [38]. In that study, continuous URF was used to evaluate the resonance frequencies of ex-vivo animal bone. The premise behind their method was that bone fracture and the level of repair could be determined based on changes in the bone resonance frequency(ies).
In this paper, we propose a quantitative URF method for detecting bone fracture and monitoring bone healing in human subjects. We applied focused URF excitations at four equi-distant points on the bone of interest. Resulting acoustic signals were used for wave arrival time and velocity estimation. The arrival time of the received signal at the hydrophone from each excitation included the effect of overlaying soft tissue. However, by measuring the time difference between the arrival times of the signals originating from two excitation points at known positions, it is possible to calculate the speed in bone in the vicinity of the excitation points. The signal measured by the hydrophone can be considered as a superposition of different modes or waves. Variational mode decomposition (VMD) [39] is an efficient method for decomposing a signal into an ensemble of band-limited intrinsic mode functions, and we used it to decompose the received signal into its constitutive modes and estimate the arrival time and velocity from the decomposed modes. Our proposed method successfully distinguished between fractured and normal bone with an 80% probability and was able to monitor the healing progress in 80% of the clavicular fractures.

2. Materials and Methods

2.1. Human Subjects

This prospective study has been conducted under an approved Mayo Clinic IRB (institutional review board) Protocol and is HIPPA (Health Insurance Portability and Accountability Act) Compliant. The study included 15 children with an obtained IRB approved informed consent/assent form, recruited between August 2016 and October 2017. Of the total, 12 had a single clavicle fracture and three had a single ulna fracture. Median age was 12 ± 4.1 years. There were 10 males (median age 11 ± 3.9 years) and five females (median age 8 ± 4.6). Signed written consent was obtained from the parents/guardians of each subject. Participants who could attend multiple visits post-fracture were considered for assessing bone healing and included 13 subjects with a clavicle fracture and three with an ulna fracture. The median times between the fracture date and the dates of the first, second, and third visits were 8 ± 2.5 days, 30 ± 8.7 days, and 45 ± 16.3 days, respectively. Clinical reports, including X-ray results and other clinical indicators reported by the physician, were used for the definitive assessment of fracture healing.

2.2. Remote Acoustic Measurements Setup

A fully programmable ultrasound platform (Verasonics V-1, Kirkland, WA, USA) equipped with a linear array transducer (L7-4) was used to guide and generate the URF excitation. The URF signal consists of a five-cycle tone burst at a 5 MHz center frequency. The magnitude radiation force exerted by an acoustic plane wave in a 5 MHz center frequency transducer depends on the absorption coefficient and average acoustic intensity wave. The absorption coefficient directly depends on the frequency and the average acoustic intensity wave directly depends on the voltage of the transducer. By increasing the frequency, the absorption coefficient is increased, while the maximum tolerable voltage of each element in the transducer is decreased. We selected the frequency and voltage in order to achieve the maximum radiation force. Among the linear array transducers (Verasonics research system), L7-4 was the best candidate for our experimental setup. A Reson TC4034-1 hydrophone (Denmark) was used to detect the radiated acoustic pressure response. The safety measurements (i.e., temperature index and mechanical index) for the setup were performed; both were within the Food and Drug Administration limit. The hydrophone output was then filtered with an analogue eight-pole band-pass Bessel filter with a center frequency of 500 kHz, digitized at a sampling frequency of 100 Mega samples per second, and recorded on a computer. The transmitting linear array and the receiving hydrophone (RX) were placed in a fixture at a distance of 0.5 cm apart. Acoustic gel was used for coupling the hydrophone/probe to the skin. The experimental setup can be seen in Figure 1.

2.3. Problem Statement

We defined two hypothesis tests: one for distinguishing between intact and fractured bones and another for identifying the progression of the healing process in fractured bone. The first hypothesis test was defined based on the velocity of the wave received by the hydrophone. The second hypothesis test was defined based on the comparison of the wave velocity distribution in the fractured bone against that of the corresponding contralateral intact bone at different visits.
In the first test, the velocity of the wave received by the hydrophone was considered as a decision variable for separating fractured from intact bones, which was defined as:
{ H F : v < v T H R H I : v v T H R ,
where H F and H I denote a set of fractured and intact bones, respectively; v is the wave velocity; and v THR is the threshold of decision based on the wave velocity [9]. This hypothesis test was defined based on the fact that we would expect the wave velocity in intact bone to be higher than the wave velocity in fractured bone. To test the hypothesis in (1), we needed to estimate the wave velocity of the received signal. Hence, in the next Section 2.4.1 and Section 2.4.2, two methods are described for wave velocity estimation. For each subject, we measured the velocity of the fractured and the corresponding contralateral intact bone multiple times. Thus, for fracture detection, it was possible to compare the distributions of the velocity (or the corresponding time delay) in each pair of intact and fractured bones.
The second hypothesis test was also based on comparing wave velocity distributions between fractured bone and the corresponding contralateral intact bone in each patient. To detect the progress of the healing process of each fractured bone over time, we compared the probability distribution function (PDF) of the estimated velocity in fractured bone to that of the contralateral intact bone at each visit using a rank sum test. The derived p value from the comparison of the pair of PDFs of velocity (or time delay) in patient l at visit m is denoted by p l ( m ) . Thus, we defined a decision variable for detecting the progress of the healing process in each patient as follows:
H D l :   = { 1 m { 1 , ... , N l } : p l ( m + 1 ) > p l ( m ) 0 m { 1 , ... , N l } : p l ( m + 1 ) p l ( m ) ,
where N l is the number of visits for patient l . If p l ( m ) is a monotonically increasing function in terms of visit m , we detect the healing process; otherwise, we cannot consider that sufficient healing progress has occurred in patient l . Moreover, if the p value in the first visit is less than 0.05, i.e., p l ( 1 ) < 0.05 , and in the last visit, p l ( N l ) , is sufficiently large, the decision variable of H D l = 1 is considered as healing detection; otherwise, H D l = 0 . The percentage of successful healing detection was defined by the quality healing factor (QHF) as:
Q H F = l = 1 L H D l L × 100 ,
where L is the number of patients. The advantage of comparing the PDF of velocities (or time delay) on different measurements of intact bone and bone with fracture as opposed to having only one statistic, such as the mean or median of velocity (or time delay), is that complete information on the PDF of velocity (or time delay) is maintained for comparing intact and fractured bone.
Based on these two hypothesis tests and the use of the rank sum test, we examined the performance of our method for detecting fractures and fracture healing in human subjects.

2.4. Acoustic Parameters Estimation

Parameters were estimated from the temporal characteristics of the acoustic response in order to capture changes in the bone. The wave velocity and arrival time differences of the acoustic waves are of particular interest for bone mechanical assessment. For simplicity of notations, we omitted the volunteer index over all the defined parameters; however, in the results section, the distribution of each defined parameter was compared for fractured and intact bones over groups of volunteers. The detailed methods for parameter estimation used for evaluating the performance of the hypothesis tests are presented in the subsections that follow.

2.4.1. Wave Velocity Estimation

The velocity of the acoustic wave was estimated from the arrival time difference measurements along the bone at each visit. The probe/hydrophone assembly was positioned along the bone, such that an imaginary line between the probe and the hydrophone was parallel to the long axis of the bone. In the case of fractured bone, the fracture position was marked on the skin by the clinical staff. Using this mark as the guide, the probe/hydrophone assembly was placed over the fracture, such that the focus of the URF was on one side of the fracture while the hydrophone was on the other side of the fracture, as shown in Figure 2a. To make the measurements on the contralateral intact bone, the probe/hydrophone assembly was placed at an approximately similar position.
Multiple measurements were performed by focusing the URF at different positions on the bone. The distance between the URF focus (TX) and the receiving hydrophone (RX) is d i = d 0 + i Δ d   cm , where d 0 = 1.2   cm , i { 0 , 1 , 2 , 3 } , and Δ d = 1.2   mm (Figure 2b). For the measured signal at distance d i , the k th measurement of m th patient visit is denoted by x i ( k , m ) ( n T s ) , where T s = 1 / F s is the sampling interval and F s is the sampling frequency, as follows:
x i ( k , m ) : = [ x i ( k , m ) [ 1 ] , ... , x i ( k , m ) [ N ] ] ,
where x i ( k , m ) [ n ] = x i ( k , m ) ( n T s ) is the n th sample at distance d i and k th measurement of m th patient visit, 1 k K , 1 n N , 1 m M . The variables K , N , and M are the number of measurements at each patient visit, the number of samples in each measurement, and the number of patient visits, respectively. We measured at K different sites on the bone. At all sites, the fracture was between the excitation point of the probe and the hydrophone.
By cross-correlating measured signals of x i ( k , m ) and x j ( k , m ) , where i , j { 0 , 1 , 2 , 3 } and i j , the bone wave velocity of k th measurement and m t h patient visit was estimated as follows:
v i j ( k , m ) :   = ( j i ) Δ d Δ τ i j ( k , m ) ,
where Δ τ i j ( k , m ) is the time delay between the signals of x i ( k , m ) and x j ( k , m ) , which was derived from the following optimization problem:
Δ τ i j ( k , m ) :   = T s arg max R ^ i j ( k , m ) ( Δ n ) s . t . ( j i ) Δ d v max Δ τ i j ( k ) ( j i ) Δ d v min
where v max and v min are, respectively, the maximum and minimum values of the allowed velocity in the propagation medium (i.e., bone); and R ^ i j ( k ) ( Δ n ) is the estimated weighted cross correlation between x i ( k , m ) and x j ( k , m ) , which was obtained by:
R ^ i j ( k , m ) ( Δ n ) = n = 1 N C x i ( k , m ) [ n ] w [ n ] x j ( k , m ) [ n + Δ n ] w [ n + Δ n ] ,
where i , j { 0 , 1 , 2 , 3 } ,   j < i , and N C are the number of samples used for the estimation of the cross correlation. Also, w [ n ] : = w ( n T s ) denotes the Tukey window (tapered cosine) in the time domain and was defined as:
w ( t ) :   = { 1 2 ( 1 + cos ( 2 π r ( t r 2 ) ) ) , 0 t < r 2 1 , r 2 t < 1 r 2 1 2 ( 1 + cos ( 2 π r ( t 1 + r 2 ) ) ) , 1 r 2 t 1 ,
where 0 r 1 is the roll-off factor. To compare the bones with and without fracture, we estimated the sample mean, sample median, and sample standard deviation of the estimated time delay and bone wave velocity for each measurement and each volunteer. The mean, median, and standard deviation of Δ τ i j ( k , m ) at m th patient visit were derived as follows:
Δ τ ¯ ( m ) :   = 1 6 K k = 1 K j = i 3 i = 0 3 Δ τ i j ( k , m )
Δ τ ˜ ( m ) :   = m e d i a n ( Δ τ i j ( k , m ) )
σ Δ τ ( m ) :   = 1 6 K 1 k = 1 K j = i 3 i = 0 3 ( Δ τ i j ( k , m ) Δ τ ¯ ( m ) ) 2 .
In a similar way the mean, median, and standard deviation of v i j ( k , m ) at m th patient visit were obtained by:
v ¯ ( m ) :   = 1 6 K k = 1 K j = i 3 i = 0 3 v i j ( k , m )
v ˜ ( m ) :   = m e d i a n ( v i j ( k , m ) )
σ v ( m ) :   = 1 6 1 k = 1 K j = i 3 i = 0 3 ( v i j ( k , m ) v ¯ ) .

2.4.2. Variational Mode Decomposition

The signal received by the hydrophone may be considered as a superposition of different components. The measured signal was decomposed into a discrete number of sub-signals (modes), u p , i ( k , m ) , which had specific sparsity properties related to the bandwidth of the mode. These modes are reproducing the measured signal. The bandwidth of u p , i ( k , m ) , using a single side-band frequency spectrum, was derived for u p , i ( k , m ) from the corresponding analytic signal with a Hilbert transform. By multiplying an exponential with the estimated center frequency, the mode’s frequency spectrum was demodulated and transferred to the baseband. Bandwidth was estimated using a squared L 2 -norm of the gradient. Therefore, the resulting constrained variational problem for x i ( k , m ) ( t ) was obtained by minimizing the sum of the gradients of the Hilbert transform of the base band modes [39] as follows:
min { u p , i ( k , m ) } , { ω p , i ( k , m ) } { p = 1 P t [ ( δ ( t ) + j π t ) u p , i ( k , m ) ( t ) ] e j ω p t 2 2 } s . t . p = 1 P u p , i ( k , m ) = x i ( k , m ) ( t ) ,
where { u p , i ( k , m ) } = { u 1 , i ( k , m ) , ... , u P , i ( k , m ) } is the set of all modes; { ω p , i ( k , m ) } = { ω 1 , i ( k , m ) , ... , ω P , i ( k , m ) } is the set of center frequencies; P is the number of modes; δ ( t ) is the Dirac delta function; and and 2 2 denote the convolution operation and L 2 norm, respectively. The augmented Lagrangian function for the optimization of (10) was obtained by:
L ( { u p , i ( k , m ) } , { ω p , i ( k , m ) } , λ i ( k , m ) ) : = α p = 1 P t ( δ ( t ) + j π t ) u p , i ( k , m ) ( t ) e j ω p t 2 2 + x i ( k , m ) ( t ) p = 1 P u p , i ( k , m ) ( t ) 2 2 + < λ i ( k , m ) ( t ) , x i ( k , m ) ( t ) p = 1 P u p , i ( k , m ) ( t ) >
where α is the weight of the bandwidth constraint and λ i ( k , m ) ( t ) is the Lagrange multiplier. The solution of the original minimization problem was a saddle point of the augmented Lagrangian in a sequence of iterative sub-optimizations called the alternate direction method of multipliers [40] (see the Appendix for more details about the steps of optimization in VMD).
Therefore, by estimating modes of the measured signal using the optimization algorithm of VMD for x i ( k , m ) , we selected two modes for velocity estimation: the first mode with the minimum center frequency, u 1 , i ( k , m ) , and the mode with the maximum power, u max , i ( k , m ) = max p | u p , i ( k , m ) ( t ) | . The arrival times difference and velocity estimation were derived, respectively, from (5) and (4) by replacing x i ( k , m ) with u 1 , i ( k , m ) or u max , i ( k , m ) in (6).

2.4.3. Accuracy of Time Delay Estimation

To assess the accuracy of the time delay, Δ τ , estimation in intact bones as the control group, we estimated the coefficient of variation (CV) of Δ τ for intact bone over all measurement sites and participant visits with different method of time delay estimation, denoted by C V t o t a l ( Δ τ ) , which was defined as:
C V t o t a l :   = σ t o t a l μ t o t a l ,
where
μ t o t a l :   = 1 6 K M m = 1 M k = 1 K j = i 3 i = 0 3 Δ τ i j ( k , m )
σ t o t a l :   = 1 6 K M 1 m = 1 M k = 1 K j = i 3 i = 0 3 ( Δ τ i j ( k , m ) μ t o t a l ) 2 .
Also, we compared C V for m th patient visit as follows:
C V ( m ) :   = σ Δ τ ( m ) Δ τ ¯ ( m ) ,
where
Δ τ ¯ ( m ) :   = 1 6 K k = 1 K j = i 3 i = 0 3 Δ τ i j ( k , m )
σ Δ τ ( m ) :   = 1 6 K 1 k = 1 K j = i 3 i = 0 3 ( Δ τ i j ( k , m ) Δ τ ¯ ( m ) ) 2 .
To assess the accuracy of Δ τ ¯ ( m ) over different patient visits, we calculated the coefficient of variation of Δ τ ¯ ( m ) , denoted by C V m e a n ( Δ τ ) , which was defined as:
C V m e a n :   = σ m e a n μ m e a n ,
where
μ m e a n :   = 1 M m = 1 M Δ τ ¯ ( m )
σ m e a n :   = 1 M 1 m = 1 M ( Δ τ ¯ ( m ) μ m e a n ) 2 .

3. Results

The vibro-acoustic measurements for the fractured and intact bones demonstrated significant differences in velocity and arrival times. The number of measurements, K , was equal to three in each bone. The roll-off factor of the Tukey window in the estimation of velocity and time delay, r , was 0.25. The processing window duration, T W , was considered to be 25 µs to analyze only the first part of the received signal by the hydrophone, which mainly includes the bone component wave because the wave propagates faster in bone than in soft tissue. In this section, v max = 4080   m / s and v min = 1000   m / s are considered for solving the optimization problem in (5); v max is obtained from the maximum velocity of the wave in the bone, and v min is obtained from a lower bound for the minimum velocity of the wave in soft tissue. Typical signals received by the hydrophone from intact and fractured bones are shown in Figure 3.
Figure 4 depicts Pr { C V t o t a l T H R } , Pr { C V m e a n T H R } , and min m Pr { C V m T H R } for the first mode, which is the mode with the minimum center frequency, the mode with the maximum power, and the received signal for participants with clavicle fracture who completed at least three visits, respectively. Here, Pr{ } denotes the probability distribution function and THR is a threshold value. As shown by the example in the figure, when we used the first mode of VMD for time delay estimation, Pr { C V m e a n 0.2 } = 0.8 , Pr { C V t o t a l 0.34 } = 0.8 , and min m Pr { C V ( m ) 0.35 } = 0.8 , indicating that the accuracy of the Δ τ estimation was fairly similar across our normal bone population.
The median of the estimated bone velocity in fractured versus normal clavicles, using v ˜ ( m ) , is plotted for comparison in Figure 5. Velocities significantly differed between groups (p < 0.05) using the rank sum test.
Performance of fracture detection based on a pair-wise comparison between intact and fractured bone for each patient is listed in Table 1. The best performance of fracture detection was 0.80 using the first mode of VMD as the input of the cross correlation technique.
Table 2 shows the estimated QHF based on the time delay and velocity for different estimation methods: constrained velocity and first mode of VMD. We found that the healing process could be detected in 80% of fractured clavicle bones with this method.

4. Discussion

Detection of the fracture and healing process of bone is critical. In this paper, a new noninvasive and quantitative URF method was proposed for the detection of bone fracture and healing. We found that waves generated by ultrasound radiation force excitation vary in speed, depending on the healing of the broken bone to which they are applied. In particular, the velocity and arrival time of the generated waves differed significantly between intact and fractured bone. Using our method, the statistical analysis revealed an 80% detection rate in bones with fracture. In monitoring the healing process, our method, based on an estimation of the time delay, demonstrated a detection rate of 80% in the clavicle bones and around 70% of both clavicle and ulna bones. Conventional ultrasound has also been reported for the detection of clavicle and forearm fractures, with variable sensitivity [7]. The development of quantitative tools for bone healing monitoring may assist physicians in identifying the healing endpoint, thus preventing prolongation of the treatment course or incorrect timing of fixation device removal. Several noninvasive, quantitative techniques to monitor bone fracture healing have been investigated [9,10,11,12,13,14,15,16,18,19]. These studies have shown a correlation between the ultrasound wave velocity and health status of bone, therefore proposing that the ultrasound wave is an accurate method for assessing fracture.
To assess reliability in the estimation of time delay and velocity in normal bone, the coefficient of variation of time delay and velocity was calculated in the normal bone population. All measurements showed that the value was almost 20%, indicating a reliable estimation of time delay and velocity.
There are some limitations to this work. First, the bone was supposed to be an isotropic material. Consequently, the elastic coefficient was assumed to be uniform in all directions. Second, the receiving hydrophone is omnidirectional, as opposed to a directionally focused transducer. Consequently, the received pressure has contributions in all directions; however, the hydrophone is held perpendicular to the bone surface so that the displacement toward the bone (deflection) will be the main contributor to the received pressure. Third, changes in the wave velocity with respect to bone density were not investigated. Velocity is not only dependent on the elastic modulus, but is also linked to density. Fourth, using VMD, we decomposed the signal received by the hydrophone to multiple modes and analyzed the mode with the maximum power for velocity estimation. The decomposed modes, though, could have been different from the waves propagating in the bone tissue material. The results of simulation and experimental ex vivo studies on cortical bone for analyzing the lamb modes using guided waves have also been reported [41,42,43,44].
Further investigation is required to improve the classification of fractured and normal bone and to be able to adequately monitor the healing process. Additional acoustic parameters, such as attenuation, may further enhance the accuracy of classification. A larger dataset can aid in a more accurate performance in identifying fractures and detecting bone healing.

5. Conclusions

Waves generated by ultrasound radiation force excitation of intact bone and bone with fracture demonstrated a significant difference in the velocity and arrival time, and in the clavicle in particular. The healing process of bone with fracture was also detected by identifying the wave velocity in fractured bone over time and comparing it to corresponding intact bone. Results indicated that waves generated by ultrasound radiation force excitation vary in velocity, depending on the health of the bone to which they are applied.

Author Contributions

Conceptualization and methodology implementation: S.G., M.D, and M.B.; Data Analysis S.G.; Investigation: S.G.; Validation: S.G; Supervision: T.A.M., A.N.L., A.A., and M.F.; Patient data collection: S.G., A.G., J.W., and V.K.; Visualization: S.G.; Writing – Original Draft Preparations: S.G.; Writing – Review & Editing: S.G., M.F., and A.A.; Funding acquired by A.A. and M.F.

Funding

This work was supported by the Mayo Clinic Ultrasound Center Grant Award: 1/2017-1/2018.

Acknowledgments

The authors are grateful to Vickie M. Treder for her assistance in patient recruitment.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

The algorithm used for optimization in VMD is summarized below. Here, ^ denotes the Fourier transform of . In this algorithm, for simplicity of notation, superscript ( k , m ) and the subscript of i are omitted, hence x ^ i ( k , m ) ( ω ) and u ^ i , p ( k , m ) ( ω ) are denoted by x ^ ( ω ) and u ^ p ( ω ) , respectively.
Algorithm A1. Steps of optimization for variational mode decomposition.
Initilize   { u ^ p 1 } , { ω p 1 } , λ ^ 1 , n 0 repeat
    n n + 1 for   p = 1 : P   do     Update   u ^ p   for   all   ω 0 :
           u ^ p n + 1 ( ω ) x ^ ( ω ) i < p u ^ i n + 1 ( ω ) i > p u ^ i n ( ω ) + λ ^ n ( ω ) 2 1 + 2 α ( ω ω p n ) 2
        Update   ω p :     ω p n + 1 0 ω | u ^ p n + 1 ( ω ) | 2 d ω 0 | u ^ p n + 1 ( ω ) | 2 d ω
    end   for Dual   ascent   for   all   ω 0 : λ ^ n + 1 ( ω ) λ ^ n ( ω ) + τ ( x ^ ( ω ) p u ^ p n + 1 ( ω ) )
until   convergence : p = 1 P u ^ p n + 1 u ^ p n 2 2 / u ^ p 2 2 < ε

References

  1. Shore, B.J.; Hutchinson, S.; Harris, M.; Bae, D.S.; Kalish, L.A.; Maxwell, W.; Waters, P. Epidemiology and prevention of cast saw injuries: results of a quality improvement program at a single institution. J. Bone Joint Sur. Am. 2014, 96, e31. [Google Scholar] [CrossRef] [PubMed]
  2. You, L.; Temiyasathit, S.; Lee, P.; Kim, C.H.; Tummala, P.; Yao, W.; Kingery, W.; Malone, A.M.; Kwon, R.Y.; Jacobs, C.R. Osteocytes as mechanosensors in the inhibition of bone resorption due to mechanical loading. Bone 2008, 42, 172–179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Blokhuis, T.J.; De Bruine, J.H.D.; Bramer, J.A.M.; Boer, F.C.D.; Bakker, F.C.; Patka, P.; Haarman, H.J.T.M.; Manoliu, R.A. The reliability of plain radiography in experimental fracture healing. Skeletal. Radiol. 2001, 30, 151–156. [Google Scholar] [CrossRef] [PubMed]
  4. E Berg, E.; Chebuhar, C.; Bell, R.M. Pelvic trauma imaging: A blinded comparison of computed tomography and roentgenograms. J. Trauma: Inj. Infect. Crit. Care 1996, 41, 994–998. [Google Scholar] [CrossRef]
  5. Nicholls, P.J.; Berg, E.; Bliven, F.E.; Kling, J.M. X-Ray Diagnosis of Healing Fractures in Rabbits. Clin Orthop Relat Res 1979, 234–236. [Google Scholar] [CrossRef]
  6. Tiedeman, J.J.; Lippiello, L.; Connolly, J.F.; Strates, B.S. Quantitative Roentgenographic Densitometry for Assessing Fracture Healing. Clin Orthop Relat Res 1990, 279–286. [Google Scholar] [CrossRef]
  7. Hamer, D.D.-D.; Blanker, M.H.; Edens, M.A.; Buijteweg, L.N.; Boomsma, M.F.; Van Helden, S.H.; Mauritz, G.-J. Ultrasound for Distal Forearm Fracture: A Systematic Review and Diagnostic Meta-Analysis. PloS ONE 2016, 11, e0155659. [Google Scholar]
  8. Moed, B.R.; Kim, E.C.; Van Holsbeeck, M.; Schaffler, M.B.; Subramanian, S.; Bouffard, J.A.; Craig, J.G. Ultrasound for the Early Diagnosis of Tibial Fracture Healing After Static Interlocked Nailing Without Reaming: Histologic Correlation Using a Canine Model. J. Orthop. Trauma 1998, 12, 200–205. [Google Scholar] [CrossRef] [PubMed]
  9. Protopappas, V.; Vavva, M.; Fotiadis, D.; Malizos, K. Ultrasonic monitoring of bone fracture healing. IEEE Trans. Ultrasonics Ferroelectr. Freq. Control 2008, 55, 1243–1255. [Google Scholar] [CrossRef] [PubMed]
  10. Cunningham, J.L.; Kenwright, J.; Kershaw, C.J. Biomechanical measurement of fracture healing. J. Med. Eng. Technol. 1990, 14, 92–101. [Google Scholar] [CrossRef] [PubMed]
  11. Nakatsuchi, Y.; Tsuchikane, A.; Nomura, A. Assessment of fracture healing in the tibia using the impulse response method. J. Orthop. Trauma 1996, 10, 50–62. [Google Scholar] [CrossRef] [PubMed]
  12. Nikiforidis, G.; Bezerianos, A.; Dimarogonas, A.; Sutherland, C. Monitoring of fracture healing by lateral and axial vibration analysis. J. Biomech. 1990, 23, 323–330. [Google Scholar] [CrossRef]
  13. Hirasawa, Y.; Takai, S.; Kim, W.C.; Takenaka, N.; Yoshino, N.; Watanabe, Y. Biomechanical monitoring of healing bone based on acoustic emission technology. Clin. Orthop. Relat. Res. 2002, 402, 236–244. [Google Scholar] [CrossRef]
  14. Watanabe, Y.; Takai, S.; Arai, Y.; Yoshino, N.; Hirasawa, Y. Prediction of mechanical properties of healing fractures using acoustic emission. J. Orthop. Res. 2001, 19, 548–553. [Google Scholar] [CrossRef] [Green Version]
  15. Lowet, G.; Van der Perre, G. Ultrasound velocity measurement in long bones: measurement method and simulation of ultrasound wave propagation. J. Biomech. 1996, 29, 1255–1262. [Google Scholar] [CrossRef]
  16. Raum, K.; Leguerney, I.; Chandelier, F.; Bossy, E.; Talmant, M.; Saïed, A.; Peyrin, F.; Laugier, P. Bone microstructure and elastic tissue properties are reflected in QUS axial transmission measurements. Ultrasound Med. Biol. 2005, 31, 1225–1235. [Google Scholar] [CrossRef] [PubMed]
  17. Protopappas, V.C.; Baga, D.A.; Fotiadis, D.I.; Likas, A.C.; Papachristos, A.A.; Malizos, K.N. An ultrasound wearable system for the monitoring and acceleration of fracture healing in long bones. IEEE Trans. Bio-Med. Eng. 1597, 52, 1597–1608. [Google Scholar] [CrossRef] [PubMed]
  18. Machado, C.B.; Pereira, W.C.D.A.; Granke, M.; Talmant, M.; Padilla, F.; Laugier, P. Experimental and simulation results on the effect of cortical bone mineralization in ultrasound axial transmission measurements: A model for fracture healing ultrasound monitoring. Bone 2011, 48, 1202–1209. [Google Scholar] [CrossRef] [PubMed]
  19. Machado, C.B.; Pereira, W.C.D.A.; Talmant, M.; Padilla, F.; Laugier, P. Computational Evaluation of the Compositional Factors in Fracture Healing Affecting Ultrasound Axial Transmission Measurements. Ultrasound Med. Biol. 2010, 36, 1314–1326. [Google Scholar] [CrossRef] [PubMed]
  20. Bochud, N.; Vallet, Q.; Minonzio, J.G.; Laugier, P. Predicting bone strength with ultrasonic guided waves. Sci. Rep. 2017, 7, 43628. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Ziegert, J.C.; Lewis, J.L. The Effect of Soft Tissue on Measurements of Vibrational Bone Motion by Skin-Mounted Accelerometers. J. Biomech. Eng. 1979, 101, 218–220. [Google Scholar] [CrossRef]
  22. Fatemi, M.; Greenleaf, J.F. Ultrasound-stimulated vibro-acoustic spectrography. Science 1998, 280, 82–85. [Google Scholar] [CrossRef] [PubMed]
  23. Fatemi, M.; Greenleaf, J.F. Vibro-acoustography: An imaging modality based on ultrasound-stimulated acoustic emission. Proc. Natl. Acad. Sci. USA 1999, 96, 6603–6608. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Alizad, A.; Fatemi, M.; Nishimura, R.A.; Kinnick, R.R.; Rambod, E.; Greenleaf, J.F. Detection of calcium deposits on heart valve leaflets by vibro-acoustography: An in vitro study. J. Am. Soc. Echocardiogr. 2002, 15, 1391–1395. [Google Scholar] [CrossRef] [PubMed]
  25. Urban, M.W.; Silva, G.T.; Fatemi, M.; Greenleaf, J.F. Multifrequency vibro-acoustography. IEEE Trans. Med. Imag. 2006, 25, 1284–1295. [Google Scholar] [CrossRef]
  26. Urban, M.W.; Chalek, C.; Haider, B.; Thomenius, K.E.; Fatemi, M.; Alizad, A. A beamforming study for implementation of vibro-acoustography with a 1.75-D array transducer. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2013, 60, 535–551. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Urban, M.W.; Chalek, C.; Kinnick, R.R.; Kinter, T.M.; Haider, B.; Greenleaf, J.F.; Thomenius, K.E.; Fatemi, M. Implementation of vibro-acoustography on a clinical ultrasound system. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2011, 58, 1169–1181. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Mitri, F.G.; Greenleaf, J.F.; Fatemi, M. Chirp imaging vibro-acoustography for removing the ultrasound standing wave artifact. IEEE Trans. Med. Imag. 2005, 24, 1249–1255. [Google Scholar] [CrossRef] [PubMed]
  29. Fatemi, M.; Greenleaf, J.F. Probing the dynamics of tissue at low frequencies with the radiation force of ultrasound. Phys. Med. Biol. 2000, 45, 1449–1464. [Google Scholar] [CrossRef] [PubMed]
  30. Silva, G.T.; Greenleaf, J.F.; Fatemi, M. Linear arrays for vibro-acoustography: a numerical simulation study. Ultrasonic Imag. 2004, 26, 1–17. [Google Scholar] [CrossRef] [PubMed]
  31. Silva, G.T.; Chen, S.; Frery, A.C.; Greenleaf, J.F.; Fatemi, M. Stress field forming of sector array transducers for vibro-acoustography. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2005, 52, 1943–1951. [Google Scholar] [CrossRef] [PubMed]
  32. Alizad, A.; Whaley, D.H.; Urban, M.W.; Carter, R.E.; Kinnick, R.R.; Greenleaf, J.F.; Fatemi, M. Breast vibro-acoustography: initial results show promise. Breast Cancer Res. 2012, 14, R128. [Google Scholar] [CrossRef] [PubMed]
  33. Alizad, A.; Wold, L.E.; Greenleaf, J.F.; Fatemi, M. Imaging mass lesions by vibro-acoustography: Modeling and experiments. IEEE Trans. Med. Imag. 2004, 23, 1087–1093. [Google Scholar] [CrossRef] [PubMed]
  34. Fatemi, M.; Wold, L.E.; Alizad, A.; Greenleaf, J.F. Vibro-acoustic tissue mammography. IEEE Trans. Med. Imag. 2002, 21, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Kamimura, H.A.; Wang, L.; Carneiro, A.A.; Kinnick, R.R.; An, K.N.; Fatemi, M. Vibroacoustography for the assessment of total hip arthroplasty. Clinics 2013, 68, 463–468. [Google Scholar] [CrossRef]
  36. Mitri, F.G.; Davis, B.J.; Alizad, A.; Greenleaf, J.F.; Wilson, T.M.; Mynderse, L.A.; Fatemi, M. Prostate cryotherapy monitoring using vibroacoustography: Preliminary results of an ex vivo study and technical feasibility. IEEE Trans. Bio-Med. Eng. 2008, 55, 2584–2592. [Google Scholar] [CrossRef] [PubMed]
  37. Denis, M.; Wan, L.; Fatemi, M.; Alizad, A. Ultrasound Characterization of Bone Demineralization Using a Support Vector Machine. Ultrasound Med. Biol. 2018, 44, 714–725. [Google Scholar] [CrossRef] [PubMed]
  38. Alizad, A.; Walch, M.; Greenleaf, J.F.; Fatemi, M. Vibrational characteristics of bone fracture and fracture repair: application to excised rat femur. J. Biomech. Eng. 2006, 128, 300–308. [Google Scholar] [CrossRef] [PubMed]
  39. Dragomiretskiy, K.; Zosso, D. Variational Mode Decomposition. IEEE Trans. Signal Proces. 2014, 62, 531–544. [Google Scholar] [CrossRef]
  40. Bertsekas, D.P. Constrained optimization and Lagrange multiplier methods; Academic Press: New York, NY, USA, 1982. [Google Scholar]
  41. Minonzio, J.G.; Bochud, N.; Vallet, Q.; Bala, Y.; Ramiandrisoa, D.; Follet, H.; Mitton, D.; Laugier, P. Bone cortical thickness and porosity assessment using ultrasound guided waves: An ex vivo validation study. Bone 2018, 116, 111–119. [Google Scholar] [CrossRef] [PubMed]
  42. Xu, K.; Minonzio, J.; Ta, D.; Hu, B.; Wang, W.; Laugier, P. Sparse SVD Method for High-Resolution Extraction of the Dispersion Curves of Ultrasonic Guided Waves. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2016, 63, 1514–1524. [Google Scholar] [CrossRef] [PubMed]
  43. Nauleau, P.; Minonzio, J.G.; Chekroun, M.; Cassereau, D.; Laugier, P.; Prada, C.; Grimal, Q. A method for the measurement of dispersion curves of circumferential guided waves radiating from curved shells: Experimental validation and application to a femoral neck mimicking phantom. Phys. Med. Biol. 2016, 61, 4746–4762. [Google Scholar] [CrossRef] [PubMed]
  44. Potsika, V.T.; Spiridon, I.F.; Protopappas, V.C.; Vavva, M.G.; Lymperopoulos, P.D.; Massalas, C.V.; Polyzos, D.K.; Fotiadis, D.I. Computational study of the influence of callus porosity on ultrasound propagation in healing bones. In Proceedings of the 2014 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Chicago, IL, USA, 26–30 August 2014; pp. 684–687. [Google Scholar]
Figure 1. Experimental setup of URF bone measurement including: (a) Verasonics V-1 machine and probe-hydrophone assembly; (b) ultrasound transducer and hydrophone assembly placement on skin. Red arrow shows the hydrophone and blue arrow shows the ultrasound transducer.
Figure 1. Experimental setup of URF bone measurement including: (a) Verasonics V-1 machine and probe-hydrophone assembly; (b) ultrasound transducer and hydrophone assembly placement on skin. Red arrow shows the hydrophone and blue arrow shows the ultrasound transducer.
Sensors 19 00955 g001
Figure 2. (a). Probe/hydrophone on the skin, where the fracture site is between the probe and hydrophone; (b) Conventional B-mode ultrasound-guided ultrasound with four radiation force excitation beam focuses (red circles) Δ d apart. Time delay measurements are conducted between two excitations. The hydrophone is located at the left side of the image. The arrow shows the bone surface.
Figure 2. (a). Probe/hydrophone on the skin, where the fracture site is between the probe and hydrophone; (b) Conventional B-mode ultrasound-guided ultrasound with four radiation force excitation beam focuses (red circles) Δ d apart. Time delay measurements are conducted between two excitations. The hydrophone is located at the left side of the image. The arrow shows the bone surface.
Sensors 19 00955 g002
Figure 3. A typical time domain signal received by the hydrophone from intact clavicle bone and clavicle bone with fracture; the signal starting point is 14 µs.
Figure 3. A typical time domain signal received by the hydrophone from intact clavicle bone and clavicle bone with fracture; the signal starting point is 14 µs.
Sensors 19 00955 g003
Figure 4. Pr { C V t o t a l T H R } , Pr { C V m e a n T H R } , and min m Pr { C V m T H R } for the first mode, mode with maximum power, and received signal for the clavicle bone with completed visit, respectively.
Figure 4. Pr { C V t o t a l T H R } , Pr { C V m e a n T H R } , and min m Pr { C V m T H R } for the first mode, mode with maximum power, and received signal for the clavicle bone with completed visit, respectively.
Sensors 19 00955 g004
Figure 5. Box-Whisker plots of the median of estimated bone velocity in fractured clavicles versus normal clavicles, p < 0.05.
Figure 5. Box-Whisker plots of the median of estimated bone velocity in fractured clavicles versus normal clavicles, p < 0.05.
Sensors 19 00955 g005
Table 1. Performance of fracture detection based on a pair-wise comparison between intact and fractured bone.
Table 1. Performance of fracture detection based on a pair-wise comparison between intact and fractured bone.
Input of Delay Estimation AlgorithmDecision VariableThreshold on p Value
0.050.1
Constrained velocityTime delay40%60%
Velocity73%73%
First modeTime delay80%80%
Velocity67%80%
Table 2. Estimated quality of healing factor (QHF) based on time delay and velocity for different estimation methods: received signal and first mode of variational mode decomposition (VMD).
Table 2. Estimated quality of healing factor (QHF) based on time delay and velocity for different estimation methods: received signal and first mode of variational mode decomposition (VMD).
Estimation MethodBone TypeDecision VariableQHF
Constrained velocityClavicleTime delay60%
Velocity60%
Clavicle + UlnaTime delay54%
Velocity62%
First mode of VMDClavicleTime delay80%
Velocity70%
Clavicle + UlnaTime delay69%
Velocity62%

Share and Cite

MDPI and ACS Style

Ghavami, S.; Gregory, A.; Webb, J.; Bayat, M.; Denis, M.; Kumar, V.; Milbrand, T.A.; Larson, A.N.; Fatemi, M.; Alizad, A. Ultrasound Radiation Force for the Assessment of Bone Fracture Healing in Children: An In Vivo Pilot Study. Sensors 2019, 19, 955. https://doi.org/10.3390/s19040955

AMA Style

Ghavami S, Gregory A, Webb J, Bayat M, Denis M, Kumar V, Milbrand TA, Larson AN, Fatemi M, Alizad A. Ultrasound Radiation Force for the Assessment of Bone Fracture Healing in Children: An In Vivo Pilot Study. Sensors. 2019; 19(4):955. https://doi.org/10.3390/s19040955

Chicago/Turabian Style

Ghavami, Siavash, Adriana Gregory, Jeremy Webb, Mahdi Bayat, Max Denis, Viksit Kumar, Todd A. Milbrand, A. Noelle Larson, Mostafa Fatemi, and Azra Alizad. 2019. "Ultrasound Radiation Force for the Assessment of Bone Fracture Healing in Children: An In Vivo Pilot Study" Sensors 19, no. 4: 955. https://doi.org/10.3390/s19040955

APA Style

Ghavami, S., Gregory, A., Webb, J., Bayat, M., Denis, M., Kumar, V., Milbrand, T. A., Larson, A. N., Fatemi, M., & Alizad, A. (2019). Ultrasound Radiation Force for the Assessment of Bone Fracture Healing in Children: An In Vivo Pilot Study. Sensors, 19(4), 955. https://doi.org/10.3390/s19040955

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop