A method for delineation of bone surfaces in photoacoustic computed tomography of the finger

Photoacoustic imaging of interphalangeal peripheral joints is of interest in the context of using the synovial membrane as a surrogate marker of rheumatoid arthritis. Previous work has shown that ultrasound produced by absorption of light at the epidermis reflects on the bone surfaces within the finger. When the reflected signals are backprojected in the region of interest, artifacts are produced, confounding interpretation of the images. In this work, we present an approach where the photoacoustic signals known to originate from the epidermis, are treated as virtual ultrasound transmitters, and a separate reconstruction is performed as in ultrasound reflection imaging. This allows us to identify the bone surfaces. Further, the identification of the joint space is important as this provides a landmark to localize a region-of-interest in seeking the inflamed synovial membrane. The ability to delineate bone surfaces allows us not only to identify the artifacts, but also to identify the interphalangeal joint space without recourse to new US hardware or a new measurement. We test the approach on phantoms and on a healthy human finger

Attenuation Imaging with Pulse-Echo Ultrasound based on an Acoustic Reflector

Ultrasound attenuation is caused by absorption and scattering in tissue and is thus a function of tissue composition, hence its imaging offers great potential for screening and differential diagnosis. In this paper we propose a novel method that allows to reconstruct spatial attenuation distribution in tissue based on computed tomography, using reflections from a passive acoustic reflector. This requires a standard ultrasound transducer operating in pulse-echo mode, thus it can be implemented on conventional ultrasound systems with minor modifications. We use calibration with water measurements in order to normalize measurements for quantitative imaging of attenuation. In contrast to earlier techniques, we herein show that attenuation reconstructions are possible without any geometric prior on the inclusion location or shape. We present a quantitative evaluation of reconstructions based on simulations, gelatin phantoms, and ex-vivo bovine skeletal muscle tissue, achieving contrast-to-noise ratio of up to 2.3 for an inclusion in ex-vivo tissue.Comment: Accepted at MICCAI 2019 (International Conference on Medical Image Computing and Computer Assisted Intervention

Toward Speed-of-Sound Anisotropy Quantification in Muscle With Pulse-Echo Ultrasound

The velocity of ultrasound longitudinal waves (speed of sound) is emerging as a valuable biomarker for a wide range of diseases, including musculoskeletal disorders. Muscles are fiber-rich tissues that exhibit anisotropic behavior, meaning that velocities vary with the wave-propagation direction. Therefore, quantifying anisotropy is essential to improve velocity estimates while providing a new metric related to muscle composition and architecture. For the first time, this work presents a method to estimate speed-of-sound anisotropy in transversely isotropic tissues using pulse-echo ultrasound. We assume elliptical anisotropy and consider an experimental setup with a flat reflector parallel to the linear probe, with the muscle in between. This setup allows us to measure first-arrival reflection traveltimes using multistatic operation. Unknown muscle parameters are the orientation angle of the anisotropy symmetry axis and the velocities along and across this axis. We derive analytical expressions for the nonlinear relationship between traveltimes and anisotropy parameters, including reflector inclinations. These equations are exact for homogeneous media and are useful to estimate the effective average anisotropy in muscles. To analyze the structure of this forward problem, we formulate the inversion statistically using the Bayesian framework. We demonstrate that anisotropy parameters can be uniquely constrained by combining traveltimes from different reflector inclinations. Numerical results from wide-ranging acquisition and anisotropy properties show that uncertainties in velocity estimates are substantially lower than expected velocity differences in the muscle. Thus, our approach could provide meaningful muscle anisotropy estimates in future clinical applications

Frequency-Dependent Attenuation Reconstruction with an Acoustic Reflector

Attenuation of ultrasound waves varies with tissue composition, hence its estimation offers great potential for tissue characterization and diagnosis and staging of pathology. We recently proposed a method that allows to spatially reconstruct the distribution of the overall ultrasound attenuation in tissue based on computed tomography, using reflections from a passive acoustic reflector. This requires a standard ultrasound transducer operating in pulse-echo mode and a calibration protocol using water measurements, thus it can be implemented on conventional ultrasound systems with minor adaptations. Herein, we extend this method by additionally estimating and imaging the frequency-dependent nature of local ultrasound attenuation for the first time. Spatial distributions of attenuation coefficient and exponent are reconstructed, enabling an elaborate and expressive tissue-specific characterization. With simulations, we demonstrate that our proposed method yields a low reconstruction error of 0.04dB/cm at 1MHz for attenuation coefficient and 0.08 for the frequency exponent. With tissue-mimicking phantoms and ex-vivo bovine muscle samples, a high reconstruction contrast as well as reproducibility are demonstrated. Attenuation exponents of a gelatin-cellulose mixture and an ex-vivo bovine muscle sample were found to be, respectively, 1.4 and 0.5 on average, from images of their heterogeneous compositions. Such frequency-dependent parametrization could enable novel imaging and diagnostic techniques, as well as help attenuation compensation other ultrasound-based imaging techniques

Learning the Imaging Model of Speed-of-Sound Reconstruction via a Convolutional Formulation

Speed-of-sound (SoS) is an emerging ultrasound contrast modality, where pulse-echo techniques using conventional transducers offer multiple benefits. For estimating tissue SoS distributions, spatial domain reconstruction from relative speckle shifts between different beamforming sequences is a promising approach. This operates based on a forward model that relates the sought local values of SoS to observed speckle shifts, for which the associated image reconstruction inverse problem is solved. The reconstruction accuracy thus highly depends on the hand-crafted forward imaging model. In this work, we propose to learn the SoS imaging model based on data. We introduce a convolutional formulation of the pulse-echo SoS imaging problem such that the entire field-of-view requires a single unified kernel, the learning of which is then tractable and robust. We present least-squares estimation of such convolutional kernel, which can further be constrained and regularized for numerical stability. In experiments, we show that a forward model learned from k-Wave simulations improves the median contrast of SoS reconstructions by 63%, compared to a conventional hand-crafted line-based wave-path model. This simulation-learned model generalizes successfully to acquired phantom data, nearly doubling the SoS contrast compared to the conventional hand-crafted alternative. We demonstrate equipment-specific and small-data regime feasibility by learning a forward model from a single phantom image, where our learned model quadruples the SoS contrast compared to the conventional hand-crafted model. On in-vivo data, the simulation- and phantom-learned models respectively exhibit impressive 7 and 10 folds contrast improvements over the conventional model

Breast Cancer Assessment With Pulse-Echo Speed of Sound Ultrasound From Intrinsic Tissue Reflections: Proof-of-Concept

PURPOSE The aim of this study was to differentiate malignant and benign solid breast lesions with a novel ultrasound (US) technique, which measures speed of sound (SoS) using standard US transducers and intrinsic tissue reflections and scattering (speckles) as internal reference. MATERIALS AND METHODS This prospective, institutional review board-approved, Health Insurance Portability and Accountability Act-compliant prospective comparison study was performed with prior written informed consent from 20 women. Ten women with histological proven breast cancer and 10 with fibroadenoma were measured. A conventional US system with a linear probe was used for SoS-US (SonixTouch; Ultrasonix, Richmond, British Columbia, Canada). Tissue speckle reflections served as a timing reference for the US signals transmitted through the breasts. Relative phase inconsistencies were detected using plane wave measurements from different angular directions, and SoS images with 0.5-mm resolution were generated using a spatial domain reconstruction algorithm. The SoS of tumors were compared with the breast density of a larger cohort of 106 healthy women. RESULTS Breast lesions show focal increments ΔSoS (meters per second) with respect to the tissue background. Peak ΔSoS values were evaluated. Breast carcinoma showed significantly higher ΔSoS than fibroadenomas ([INCREMENT]SoS > 41.64 m/s: sensitivity, 90%; specificity, 80%; area under curve, 0.910) and healthy breast tissue of different densities (area under curve, 0.938; sensitivity, 90%; specificity, 96.5%). The lesion localization in SoS-US images was consistent with B-mode imaging and repeated SoS-US measurements were reproducible. CONCLUSIONS Using SoS-US, based on conventional US and tissue speckles as timing reference, breast carcinoma showed significantly higher SoS values than fibroadenoma and healthy breast tissue of different densities. The SoS presents a promising technique for differentiating solid breast lesions

A new approach to depict bone surfaces in finger imaging using photoacoustic tomography

Imaging the vasculature close around the finger joints is of interest in the field of rheumatology. Locally increased vasculature in the synovial membrane of these joints can be a marker for rheumatoid arthritis. In previous work we showed that part of the photoacoustically induced ultrasound from the epidermis reflects on the bone surface within the finger. These reflected signals could be wrongly interpreted as new photoacoustic sources. In this work we show that a conventional ultrasound reconstruction algorithm, that considers the skin as a collection of ultrasound transmitters and the PA tomography probe as the detector array, can be used to delineate bone surfaces of a finger. This can in the future assist in the localization of the joint gaps. This can provide us with a landmark to localize the region of the inflamed synovial membrane. We test the approach on finger mimicking phantoms