This paper is devoted to a method of extraction of guided waves phase velocities from experimental signals. Measurements are performed using an axial transmission device consisting of a linear arrangement of emitters and receivers placed on the surface of the inspected specimen. The technique takes benefit of using both multiple emitters and receivers and is validated on a reference wave guide. The guided mode phase velocities are obtained using a projection in the singular vectors basis. The singular vectors are determined by the singular values decomposition (SVD) of the response matrix between the two arrays in the frequency domain. This technique enables to recover accurately guided wave phase velocity dispersion curves. The SVD based approach was designed to overcome limitations of spatio-temporal Fourier transform for receiver array of limited spatial extent as in the case of clinical assessment of cortical bone in axial transmission.

This paper reports for the first time on inverse estimation of several bone properties from guided-wave measurements in human bone samples. Previously, related approaches have focused on ultrasonic estimation of a single bone property at a time. The method is based on two steps: the multi-Lamb mode response is analyzed using the singular value decomposition signal processing method recently introduced in the field, then an identification procedure is run to find thickness and anisotropic elastic properties of the considered specimen. Prior to the measurements on bone, the method is validated on cortical bone-mimicking phantoms. The repeatability and the trueness of the estimated parameters on bone-mimicking phantoms were found around a few percent. Estimation of cortical thickness on bone samples was in good agreement with cortical thickness derived from high-resolution peripheral quantitative computed tomography data analysis of the samples.

Recent bone quantitative ultrasound approaches exploit the multimode waveguide response of long bones for assessing properties such as cortical thickness and stiffness. Clinical applications remain, however, challenging, as the impact of soft tissue on guided waves characteristics is not fully understood yet. In particular, it must be clarified whether soft tissue must be incorporated in waveguide models needed to infer reliable cortical bone properties. We hypothesize that an inverse procedure using a free plate model can be applied to retrieve the thickness and stiffness of cortical bone from experimental data. This approach is first validated on a series of laboratory-controlled measurements performed on assemblies of bone- and soft tissue mimicking phantoms and then onin vivomeasurements. The accuracy of the estimates is evaluated by comparison with reference values. To further support our hypothesis, these estimates are subsequently inserted into a bilayer model to test its accuracy. Our results show that the free plate model allows retrieving reliable waveguide properties, despite the presence of soft tissue. They also suggest that the more sophisticated bilayer model, although it is more precise to predict experimental data in the forward problem, could turn out to be hardly manageable for solving the inverse problem.

The 2-D Fourier transform analysis of multichannel signals is a straightforward method to extract the dispersion curves of guided modes. Basically, the time signals recorded at several positions along the waveguide are converted to the wavenumber-frequency space, so that the dispersion curves (i.e., the frequency-dependent wavenumbers) of the guided modes can be extracted by detecting peaks of energy trajectories. In order to improve the dispersion curve extraction of low-amplitude modes propagating in a cortical bone, a multiemitter and multireceiver transducer array has been developed together with an effective singular vector decomposition (SVD)-based signal processing method. However, in practice, the limited number of positions where these signals are recorded results in a much lower resolution in the wavenumber axis than in the frequency axis. This prevents a clear identification of overlapping dispersion curves. In this paper, a sparse SVD (S-SVD) method, which combines the signal-to-noise ratio improvement of the SVD-based approach with the high wavenumber resolution advantage of the sparse optimization, is presented to overcome the above-mentioned limitation. Different penalty constraints, i.e., $l_{1}$ -norm, Frobenius norm, and revised Cauchy norm, are compared with the sparse characteristics. The regularization parameters are investigated with respect to the convergence property and wavenumber resolution. The proposed S-SVD method is investigated using synthetic wideband signals and experimental data obtained from a bone-mimicking phantom and from an ex-vivo human radius. The analysis of the results suggests that the S-SVD method has the potential to significantly enhance the wavenumber resolution and to improve the extraction of the dispersion curves.

Robust signal processing methods adapted to clinical measurements of guided modes are required to assess bone properties such as cortical thickness and porosity. Recently, an approach based on the singular value decomposition (SVD) of multidimensional signals recorded with an axial transmission array of emitters and receivers has been proposed for materials with negligible absorption, see Minonzio et al. [J. Acoust. Soc. Am. 127, 2913–2919 (2010)]. In presence of absorption, the ability to extract guided mode degrades. The objective of the present study is to extend the method to the case of absorbing media, considering attenuated plane waves (complex wavenumber). The guided mode wavenumber extraction is enhanced and the order of magnitude of the attenuation of the guided mode is estimated. Experiments have been carried out on 2 mm thick plates in the 0.2–2 MHz bandwidth. Two materials are inspected: polymethylacrylate (PMMA) (isotropic with absorption) and artificial composite bones (Sawbones, Pacific Research Laboratory Inc, Vashon, WA) which is a transverse isotropic absorbing medium. Bulk wave velocities and bulk attenuation have been evaluated from transmission measurements. These values were used to compute theoretical Lamb mode wavenumbers which are consistent with the experimental ones obtained with the SVD-based approach.