Adaptively multi-scale microstructure characterization of cancellous bone via Photoacoustic signal decomposition

Osteoporosis is a systemic disease with a high incidence in the elderly and seriously affects the quality of life of patients. Photoacoustic (PA) technology, which combines the advantages of light and ultrasound, can provide information about the physiological structure and chemical information of biological tissues in a non-invasive and non-radiative way. Due to the complex structural characteristics of bone tissue, PA signals generated by bone tissue are non-stationary and nonlinear. However, conventional PA signal processing methods are not effective for non-stationary signal processing. In this study, an empirical mode decomposition (EMD)-based Hilbert-Huang transform (HHT) PA signal analysis method, called HHT PA signal analysis (HPSA), was developed to assess the microstructure information of bone tissue, which is closely related to bone health. The feasibility of the HPSA method in bone health assessment was proven by numerical simulation and experimental studies on animal samples with different bone volume/total volume (BV/TV) and bone mineral densities. First, based on adaptive EMD, the different modes correlated with multi-scale information were mined from the PA signal, the correlations between different intrinsic mode function (IMF) modes and BV/TVs were analyzed, and the optimal mode for more efficient PA time–frequency analysis was selected. Second, multi-wavelength HPSA was used to assess the changes in the chemical components of the bone tissue. The results demonstrate that the HPSA method can distinguish bones with different BV/TVs and microstructure conditions adaptively with high efficiency. They further emphasize the potential of PA techniques in characterizing biological tissues in bones for early and rapid detection of bone diseases.