Multi-frequency characterization of the speed of sound and attenuation coefficient for longitudinal transmission of freshly excised human skulls

Abstract

For medical applications of ultrasound inside the brain, it is necessary to understand the relationship between the apparent density of skull bone and its corresponding speed of sound and attenuation coefficient. Although there have been previous studies exploring this phenomenon, there is still a need to extend the measurements to cover more of the clinically relevant frequency range. The results of measurements of the longitudinal speed of sound and attenuation coefficient are presented for specimens of human calvaria. The study was performed for the frequencies of 0.27, 0.836, 1.402, 1.965 and 2.525 MHz. Specimens were obtained from fresh cadavers through a protocol with the Division of Anatomy of the University of Toronto. The protocol was approved by the Research Ethics Board of Sunnybrook Health Sciences Centre. The specimens were mounted in polycarbonate supports that were marked for stereoscopic positioning. Computer tomography (CT) scans of the skulls mounted on their supports were performed, and a three-dimensional skull surface was reconstructed. This surface was used to guide a positioning system to ensure the normal sound incidence of an acoustic signal. This signal was produced by a focused device with a diameter of 5 cm and a focal length of 10 cm. Measurements of delay in time of flight were carried out using a needle hydrophone. Measurements of effective transmitted energy were carried out using a radiation force method with a 10 μg resolution scale. Preliminary functions of speed of sound and attenuation coefficient, both of which are related to apparent density, were established using a multi-layer propagation model that takes into account speed of sound, density and thickness of the layer. An optimization process was executed from a large set of random functions and the best functions were chosen for those ones that closest reproduced the experimental observations. The final functions were obtained after a second pass of the optimization process was executed, but this time using a finite-difference time-difference solution of the Westervelt equation, which is more precise than the multi-layer model but much more time consuming for computation. For six of seven specimens, measurements were carried out on five locations on the calvaria, and for the other specimen three measurements were made. In total, measurements were carried out on 33 locations. Results indicated the presence of dispersion effects and that these effects are different according to the type of bone in the skull (cortical and trabecular). Additionally, both the speed of sound and attenuation showed dependence on the skull density that varied with the frequency. Using the optimal functions and the information of density from the CT scans, the average values (±s.d.) of the speed of sound for cortical bone were estimated to be 2384(±130), 2471(±90), 2504(±120), 2327(±90) and 2053(±40) m s−1 for the frequencies of 270, 836, 1402, 1965 and 2526 kHz, respectively. For trabecular bone, and in the same order of frequency values, the sp...

Experimental demonstration of passive acoustic imaging in the human skull cavity using CT‐based aberration corrections

“…For a given source frequency of interest, maps of the longitudinal speed of sound and attenuation coefficient were generated by interpolating the empirical relations determined in Ref. 80. Voxels within the computational domain not corresponding to skull bone were assigned material parameters to that of water (ρ = 1000 kg m −3 , α = 0) based on the values presented in Ref.…”

Different magnetic resonance imaging patterns after transcranial magnetic resonance–guided focused ultrasound of the ventral intermediate nucleus of the thalamus and anterior limb of the internal capsule in patients with essential tremor or obsessive-compulsive disorder

“…12 However, heterogeneous aspects of the human skull such as volume, thickness, density, and shape are possible causes of defocusing, which could lead to insufficient temperature rise. 5,21 The transmission of ultrasound through the skull is reduced by the intrinsic absorption of acoustic energy by the bone and the reflection of energy at the inner and outer tables or the trabecular level of the skull. Furthermore, the target location and related acoustic energy beams' characteristics (e.g., ray angle to the skull, distance from transducer to the skull) may also affect energy delivery.…”

Sensitivity of simulated transcranial ultrasound fields to acoustic medium property maps

“…This works by simulating the propagation of ultrasound waves through a map of the individual patient's Figure 1: Relationship between the density of the skull bone, derived from CT, and the compressional sound speed. Values taken from Connor et al [19], Mast et al [20], Aubry et al [5], and Pichardo et al [21].…”