Cancellous bone analysis with modified least squares Prony’s method and chirp filter: Phantom experiments and simulation

Abstract

The presence of two longitudinal waves in porous media is predicted by Biot's theory and has been confirmed experimentally in cancellous bone. When cancellous bone samples are interrogated in through-transmission, these two waves can overlap in time. Previously, the Modified Least-Squares Prony's (MLSP) method was validated for estimation of amplitudes, attenuation coefficients, and phase velocities of fast and slow waves, but tended to overestimate phase velocities by up to about 5%. In the present paper, a pre-processing chirp filter to mitigate the phase velocity bias is derived. The MLSP/chirp filter (MLSPCF) method was tested for decomposition of a 500 kHz-center-frequency signal containing two overlapping components: one passing through a low-density-polyethylene plate (fast wave) and another passing through a cancellous-bone-mimicking phantom material (slow wave). The chirp filter reduced phase velocity bias from 100 m/s (5.1%) to 69 m/s (3.5%) (fast wave) and from 29 m/s (1.9%) to 10 m/s (0.7%) (slow wave). Similar improvements were found for 1) measurements in polycarbonate (fast wave) and a cancellous-bone-mimicking phantom (slow wave), and 2) a simulation based on parameters mimicking bovine cancellous bone. The MLSPCF method did not offer consistent improvement in estimates of attenuation coefficient or amplitude.