Measurement of ultrasound velocity in yolk and blastula of fish embryo in vivo

An acoustic microscopy method for measuring the velocity of ultrasound in the yolk and blastula of bony fish embryos at early stages of development was proposed. The yolk and blastula were approximated as a sphere and a spherical dome, respectively, consisting of a homogeneous liquid. A theoretical model of ultrasonic wave propagation through a spherical liquid drop located on a solid substrate was developed in the ray approximation. The dependence of the wave propagation time on the speed of sound in the drop, its diameter, and the position of the focus of the ultrasonic transducer has been determined. It was shown that the velocity in the drop can be found by solving the inverse problem by minimizing the discrepancy between the experimental and model spatial distributions of the propagation time, assuming that the velocity in the immersion liquid and the radius of the drop are known. The velocities in the yolk and blastula of the loach (Misgurnus fossilis) embryo at the stage of development of the middle blastula were measured in vivo using a pulsed scanning acoustic microscope operating at a central frequency of 50 MHz. The yolk and blastula radii were determined from ultrasound images of the embryo. Acoustic microscopy measurements conducted with four embryos provide velocities of the acoustic longitudinal wave in the yolk and blastula. They were measured to be 1581 ± 5 m/s and 1525 ± 4 m/s when the temperature of the liquid in the water tank was kept at 22 ± 2 °C.

Towards an integrated radiofrequency safety concept for implant carriers in MRI based on sensor-equipped implants and parallel transmission

To protect implant carriers in MRI from excessive radiofrequency (RF) heating it has previously been suggested to assess that hazard via sensors on the implant. Other work recommended parallel transmission (pTx) to actively mitigate implant-related heating. Here, both ideas are integrated into one comprehensive safety concept where native pTx safety (without implant) is ensured by state-of-the-art field simulations and the implant-specific hazard is quantified in situ using physical sensors. The concept is demonstrated by electromagnetic simulations performed on a human voxel model with a simplified spinal-cord implant in an eight-channel pTx body coil at 3 T . To integrate implant and native safety, the sensor signal must be calibrated in terms of an established safety metric (e.g., specific absorption rate [SAR]). Virtual experiments show that E -field and implant-current sensors are well suited for this purpose, while temperature sensors require some caution, and B 1 probes are inadequate. Based on an implant sensor matrix Q s , constructed in situ from sensor readings, and precomputed native SAR limits, a vector space of safe RF excitations is determined where both global (native) and local (implant-related) safety requirements are satisfied. Within this safe-excitation subspace, the solution with the best image quality in terms of B 1 + magnitude and homogeneity is then found by a straightforward optimization algorithm. In the investigated example, the optimized pTx shim provides a 3-fold higher mean B 1 + magnitude compared with circularly polarized excitation for a maximum implant-related temperature increase ∆ T imp ≤ 1 K . To date, sensor-equipped implants interfaced to a pTx scanner exist as demonstrator items in research labs, but commercial devices are not yet within sight. This paper aims to demonstrate the significant benefits of such an approach and how this could impact implant-related RF safety in MRI. Today, the responsibility for safe implant scanning lies with the implant manufacturer and the MRI operator; within the sensor concept, the MRI manufacturer would assume much of the operator's current responsibility.

Frequency-domain full-waveform inversion-based musculoskeletal ultrasound computed tomography

Recently, full-waveform inversion (FWI) has become a promising tool for ultrasound computed tomography (USCT). However, as a computationally intensive technique, FWI suffers from computational burden, especially in conventional time-domain full-waveform inversion (TDFWI). On the contrary, frequency-domain full-waveform inversion (FDFWI) provides a relatively high computational efficiency as the propagation of discrete frequencies is much cheaper than full time-domain modeling. FDFWI has already been applied in soft tissue imaging, such as breast, but for the musculoskeletal model with high impedance contrast between hard and soft tissues, there is still a lack of an effective source estimation method. In this paper, a water-referenced data calibration method is proposed to address the source estimation challenge in the presence of bones, which achieves consistency between the measured and simulated data before the FDFWI procedure. To avoid the cycle-skipping local minimum effect and facilitate the algorithm convergence, a starting frequency criterion for musculoskeletal FDFWI is further proposed. The feasibility of the proposed method is demonstrated by numerical studies on retrieving the anatomies of the leg models and different musculoskeletal lesions. The study extends the advanced FDFWI method to the musculoskeletal system and provides an alternative solution for musculoskeletal USCT imaging with high computational efficiency.

Development of ex vivo brain hemorrhage phantom for photoacoustic imaging

Brain hemorrhage, specifically intraventricular hemorrhage (IVH), is considered one of the primary and leading causes of cerebral anomalies in neonates. Several imaging modalities including the most popular, cranial ultrasound, are not capable of detecting early stage IVHs. Photoacoustic imaging (PAI) exhibited great potential for detecting cerebral hemorrhage in studies limited to small animal models, but these models are not comparable to neonatal brain morphology. However, hemorrhage detection in large animal models using PAI is rare due to the complexity and cost of inducing hemorrhage in vivo. Moreover, in vitro studies are unable to represent the physiology and environment of the hemorrhagic lesion. Here, we proposed a pseudo hemorrhage implementation method in the sheep brain that allows us to mimic different hemorrhagic lesions ex vivo without compromising the complexity of cerebral imaging. This approach enables a true evaluation of PAI performance for detecting hemorrhages and can be utilized as a reference to optimize the PAI system for in vivo imaging.

A Graph-Space Optimal Transport Approach Based on Kaniadakis κ-Gaussian Distribution for Inverse Problems Related to Wave Propagation

Data-centric inverse problems are a process of inferring physical attributes from indirect measurements. Full-waveform inversion (FWI) is a non-linear inverse problem that attempts to obtain a quantitative physical model by comparing the wave equation solution with observed data, optimizing an objective function. However, the FWI is strenuously dependent on a robust objective function, especially for dealing with cycle-skipping issues and non-Gaussian noises in the dataset. In this work, we present an objective function based on the Kaniadakis κ-Gaussian distribution and the optimal transport (OT) theory to mitigate non-Gaussian noise effects and phase ambiguity concerns that cause cycle skipping. We construct the κ-objective function using the probabilistic maximum likelihood procedure and include it within a well-posed version of the original OT formulation, known as the Kantorovich-Rubinstein metric. We represent the data in the graph space to satisfy the probability axioms required by the Kantorovich-Rubinstein framework. We call our proposal the κ-Graph-Space Optimal Transport FWI (κ-GSOT-FWI). The results suggest that the κ-GSOT-FWI is an effective procedure to circumvent the effects of non-Gaussian noise and cycle-skipping problems. They also show that the Kaniadakis κ-statistics significantly improve the FWI objective function convergence, resulting in higher-resolution models than classical techniques, especially when κ=0.6.

Acute stress induces long-term metabolic, functional, and structural remodeling of the heart

Takotsubo cardiomyopathy is a stress-induced cardiovascular disease with symptoms comparable to those of an acute coronary syndrome but without coronary obstruction. Takotsubo was initially considered spontaneously reversible, but epidemiological studies revealed significant long-term morbidity and mortality, the reason for which is unknown. Here, we show in a female rodent model that a single pharmacological challenge creates a stress-induced cardiomyopathy similar to Takotsubo. The acute response involves changes in blood and tissue biomarkers and in cardiac in vivo imaging acquired with ultrasound, magnetic resonance and positron emission tomography. Longitudinal follow up using in vivo imaging, histochemistry, protein and proteomics analyses evidences a continued metabolic reprogramming of the heart towards metabolic malfunction, eventually leading to irreversible damage in cardiac function and structure. The results combat the supposed reversibility of Takotsubo, point to dysregulation of glucose metabolic pathways as a main cause of long-term cardiac disease and support early therapeutic management of Takotsubo.

Volumetric registration framework for multimodal functional magnetic resonance and optoacoustic tomography of the rodent brain

Optoacoustic tomography (OAT) provides a non-invasive means to characterize cerebral hemodynamics across an entire murine brain while attaining multi-parametric readouts not available with other modalities. This unique capability can massively impact our understanding of brain function. However, OAT largely lacks the soft tissue contrast required for unambiguous identification of brain regions. Hence, its accurate registration to a reference brain atlas is paramount for attaining meaningful functional readings. Herein, we capitalized on the simultaneously acquired bi-modal data from the recently-developed hybrid magnetic resonance optoacoustic tomography (MROT) scanner in order to devise an image coregistration paradigm that facilitates brain parcellation and anatomical referencing. We evaluated the performance of the proposed methodology by coregistering OAT data acquired with a standalone system using different registration methods. The enhanced performance is further demonstrated for functional OAT data analysis and characterization of stimulus-evoked brain responses. The suggested approach enables better consolidation of the research findings thus facilitating wider acceptance of OAT as a powerful neuroimaging tool to study brain functions and diseases.

Full-view LED-based optoacoustic tomography

Optoacoustic tomography is commonly performed with bulky and expensive short-pulsed solid-state lasers providing high per-pulse energies in the millijoule range. Light emitting diodes (LEDs) represent a cost-effective and portable alternative for optoacoustic signal excitation that can additionally provide excellent pulse-to-pulse stability. Herein, we introduce a full-view LED-based optoacoustic tomography (FLOAT) system for deep tissue in vivo imaging. It is based on a custom-made electronic unit driving a stacked array of LEDs, which attains 100 ns pulse width and highly stable (0.62 % standard deviation) total per-pulse energy of 0.48 mJ. The illumination source is integrated into a circular array of cylindrically-focused ultrasound detection elements to result in a full-view tomographic configuration, which plays a critical role in circumventing limited-view effects, enhancing the effective field-of-view and image quality for cross-sectional (2D) imaging. We characterized the FLOAT performance in terms of pulse width, power stability, excitation light distribution, signal-to-noise and penetration depth. FLOAT of the human finger revealed a comparable imaging performance to that achieved with the standard pulsed Nd:YAG laser. It is anticipated that this compact, affordable and versatile illumination technology will facilitate optoacoustic imaging developments in resource-limited settings for biological and clinical applications.

Feasibility of spinal cord imaging at 7 T using rosette trajectory with magnetization transfer preparation and compressed sensing

MRI is a valuable diagnostic tool to investigate spinal cord (SC) pathology. SC MRI can benefit from the increased signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) at ultra-high fields such as 7 T. However, SC MRI acquisitions with routine Cartesian readouts are prone to image artifacts caused by physiological motion. MRI acquisition techniques with non-Cartesian readouts such as rosette can help reduce motion artifacts. The purpose of this study was to demonstrate the feasibility of high-resolution SC imaging using rosette trajectory with magnetization transfer preparation (MT-prep) and compressed sensing (CS) at 7 T. Five healthy volunteers participated in the study. Images acquired with rosette readouts demonstrated reduced motion artifacts compared to the standard Cartesian readouts. The combination of multi-echo rosette-readout images improved the CNR by approximately 50% between the gray matter (GM) and white matter (WM) compared to single-echo images. MT-prep images showed excellent contrast between the GM and WM with magnetization transfer ratio (MTR) and cerebrospinal fluid normalized MT signal (MTCSF) = 0.12 ± 0.017 and 0.74 ± 0.013, respectively, for the GM; and 0.18 ± 0.011 and 0.58 ± 0.009, respectively, for the WM. Under-sampled acquisition using rosette readout with CS reconstruction demonstrated up to 6 times faster scans with comparable image quality as the fully-sampled acquisition.

Synthesized Image Reconstruction for Post-Reconstruction Resolution Recovery

Resolution recovery (RR) techniques in positron emission tomography (PET) imaging aim to mitigate spatial resolution losses and related inaccuracies in quantification by using a model of the system's point spread function (PSF) during reconstruction or post-processing. However, including PSF modeling in fully 3-D image reconstruction is far from trivial as access to the scanner-specific forward and back-projectors is required, along with access to the 3-D sinogram data. Hence, post-reconstruction RR methods, such as the Richardson-Lucy (RL) algorithm, can be more practical. However, the RL method leads to relatively rapid noise amplification in early image iterations, giving inferior image quality compared to iterates obtained by placing the PSF model in the reconstruction algorithm. We propose a post-reconstruction RR method by synthesizing PET data by a forward projection of an initial real data reconstruction (such reconstructions are usually available via a scanner's standard reconstruction software). The synthetic PET data are then used to reconstruct an image, but crucially now including a modeled PSF within the system model used during reconstruction. Results from simulations and real data demonstrate the proposed method improves image quality compared to the RL algorithm, whilst avoiding the need for scanner-specific projectors and raw sinogram data (as required by standard PSF modeling within reconstruction).

Full-Waveform Inversion for Breast Ultrasound Tomography Using Line-Shape Modeled Elements

OBJECTIVE

The objective of the work described here was to incorporate the spatial shapes of the transducer elements into the framework of the full-waveform inversion.

METHODS

An element is treated as its cross-section in the 2-D imaging plane, that is, a line segment. The elements are not simply modeled as a set of point sources on their surface to avoid staircasing artifacts. By use of the Fourier collocation method, an element is spatially represented as the discrete convolution between its spatial distribution and a band-limited delta function. The excitation pulses on the emitters and recorded signals on the receivers are then weighted based on the discrete convolution results. Digital and physical experiments are implemented to validate the method.

DISCUSSION

It is meaningful to model the shapes of the elements if their spatial sizes are similar to or larger than the acoustic wavelengths. It should, however, be noted that because this article focuses on 2-D imaging, the inter-plane effects are not considered.

CONCLUSION

The approach helps reduce the root mean square errors and increase the structural similarity of the reconstructed images. It also helps to improve the stability of convergence and to accelerate the convergence speed.

Advanced Radio Frequency Applicators for Thermal Magnetic Resonance Theranostics of Brain Tumors

Thermal Magnetic Resonance (ThermalMR) is a theranostic concept that combines diagnostic magnetic resonance imaging (MRI) with targeted thermal therapy in the hyperthermia (HT) range using a radiofrequency (RF) applicator in an integrated system. ThermalMR adds a therapeutic dimension to a diagnostic MRI device. Focused, targeted RF heating of deep-seated brain tumors, accurate non-invasive temperature monitoring and high-resolution MRI are specific requirements of ThermalMR that can be addressed with novel concepts in RF applicator design. This work examines hybrid RF applicator arrays combining loop and self-grounded bow-tie (SGBT) dipole antennas for ThermalMR of brain tumors, at magnetic field strengths of 7.0 T, 9.4 T and 10.5 T. These high-density RF arrays improve the feasible transmission channel count, and provide additional degrees of freedom for RF shimming not afforded by using dipole antennas only, for superior thermal therapy and MRI diagnostics. These improvements are especially relevant for ThermalMR theranostics of deep-seated brain tumors because of the small surface area of the head. ThermalMR RF applicators with the hybrid loop+SGBT dipole design outperformed applicators using dipole-only and loop-only designs, with superior MRI performance and targeted RF heating. Array variants with a horse-shoe configuration covering an arc (270°) around the head avoiding the eyes performed better than designs with 360° coverage, with a 1.3 °C higher temperature rise inside the tumor while sparing healthy tissue. Our EMF and temperature simulations performed on a virtual patient with a clinically realistic intracranial tumor provide a technical foundation for implementation of advanced RF applicators tailored for ThermalMR theranostics of brain tumors.

Ultrasound computed tomography based on full waveform inversion with source directivity calibration

Ultrasound computed tomography based on full waveform inversion has the potential to provide high-resolution images of human tissues in a quantitative manner. A successful ultrasound computed tomography system requires the decent knowledge of acquisition array, including the spatial position and the directivity of each transducer, to meet the high-level demand of clinical applications. The conventional full waveform inversion algorithm assumes a point source with the omni-directional emission. Such assumption does not hold when the directivity of emitting transducer is not negligible. For a practical implementation, an efficient and accurate self-checking evaluation of directivity is crucial prior to the reconstruction of images. We propose to measure the directivity of each emitting transducer using the full-matrix captured data obtained with a water-immersed and target-free experiment. We introduce the weighted virtual point-source array to act as the proxy of emitting transducer during the numerical simulation. The weights of different points in the virtual array can be calculated from the observed data using the gradient-based local optimization method. Although the full waveform imaging method relies on the finite-difference solver of wave equation, such directivity estimation benefits from the introduction of analytical solver. The trick significantly reduces the numerical cost, enabling an automatic directivity self-check at boot. We verify the feasibility, efficiency, and accuracy of the virtual array method through simulated and experimental tests. For the experimental test, we also illustrate that full waveform inversion with directivity calibration can reduce the artifacts introduced by the conventional point source assumption, improving the quality of reconstructed images..

Dynamic Speed of Sound Adaptive Transmission-Reflection Ultrasound Computed Tomography

Ultrasound computed tomography (USCT) can visualize a target with multiple imaging contrasts, which were demonstrated individually previously. Here, to improve the imaging quality, the dynamic speed of sound (SoS) map derived from the transmission USCT will be adapted for the correction of the acoustic speed variation in the reflection USCT. The variable SoS map was firstly restored via the optimized simultaneous algebraic reconstruction technique with the time of flights selected from the transmitted ultrasonic signals. Then, the multi-stencils fast marching method was used to calculate the delay time from each element to the grids in the imaging field of view. Finally, the delay time in conventional constant-speed-assumed delay and sum (DAS) beamforming would be replaced by the practical computed delay time to achieve higher delay accuracy in the reflection USCT. The results from the numerical, phantom, and in vivo experiments show that our approach enables multi-modality imaging, accurate target localization, and precise boundary detection with the full-view fast imaging performance. The proposed method and its implementation are of great value for accurate, fast, and multi-modality USCT imaging, particularly suitable for highly acoustic heterogeneous medium.

Synthetic PET via Domain Translation of 3-D MRI

Historically, patient datasets have been used to develop and validate various reconstruction algorithms for PET/MRI and PET/CT. To enable such algorithm development, without the need for acquiring hundreds of patient exams, in this article we demonstrate a deep learning technique to generate synthetic but realistic whole-body PET sinograms from abundantly available whole-body MRI. Specifically, we use a dataset of 56 18F-FDG-PET/MRI exams to train a 3-D residual UNet to predict physiologic PET uptake from whole-body T1-weighted MRI. In training, we implemented a balanced loss function to generate realistic uptake across a large dynamic range and computed losses along tomographic lines of response to mimic the PET acquisition. The predicted PET images are forward projected to produce synthetic PET (sPET) time-of-flight (ToF) sinograms that can be used with vendor-provided PET reconstruction algorithms, including using CT-based attenuation correction (CTAC) and MR-based attenuation correction (MRAC). The resulting synthetic data recapitulates physiologic 18F-FDG uptake, e.g., high uptake localized to the brain and bladder, as well as uptake in liver, kidneys, heart, and muscle. To simulate abnormalities with high uptake, we also insert synthetic lesions. We demonstrate that this sPET data can be used interchangeably with real PET data for the PET quantification task of comparing CTAC and MRAC methods, achieving ≤ 7.6% error in mean-SUV compared to using real data. These results together show that the proposed sPET data pipeline can be reasonably used for development, evaluation, and validation of PET/MRI reconstruction methods.

Head-to-tail imaging of mice with spiral volumetric optoacoustic tomography

Optoacoustic tomography has been established as a powerful modality for preclinical imaging. However, efficient whole-body imaging coverage has not been achieved owing to the arduous requirement for continuous acoustic coupling around the animal. In this work, we introduce panoramic (360⁰) head-to-tail 3D imaging of mice with spiral volumetric optoacoustic tomography (SVOT). The system combines multi-beam illumination and a dedicated head holder enabling uninterrupted acoustic coupling for whole-body scans. Image fidelity is optimized with self-gated respiratory motion rejection and dual speed-of-sound reconstruction algorithms to attain spatial resolution down to 90 µm. The developed system is thus highly suitable for visualizing rapid biodynamics across scales, such as hemodynamic changes in individual organs, responses to treatments and stimuli, perfusion, total body accumulation, or clearance of molecular agents and drugs with unmatched contrast, spatial and temporal resolution.

Application of Machine learning to predict RF heating of cardiac leads during magnetic resonance imaging at 1.5 T and 3 T: A simulation study

Predicting magnetic resonance imaging (MRI)-induced heating of elongated conductive implants, such as leads in cardiovascular implantable electronic devices, is essential to assessing patient safety. Phantom experiments have traditionally been used to estimate radio-frequency (RF) heating of implants, but they are time-consuming. Recently, machine learning has shown promise for fast prediction of RF heating of orthopaedic implants when the implant position within the MRI RF coil was predetermined. We explored whether deep learning could be applied to predict RF heating of conductive leads with variable positions and orientations during MRI at 1.5 T and 3 T. Models of 600 cardiac leads with clinically relevant trajectories were generated, and electromagnetic simulations were performed to calculate the maximum of the 1 g-averaged specific absorption rate (SAR) of RF energy at the tips of lead models during MRI at 1.5 T and 3 T. Neural networks were trained to predict the maximum SAR at the lead tip from the knowledge of the coordinates of points along the lead trajectory. Despite the large range of SAR values (∼230 W/kg to ∼ 3200 W/kg and ∼ 10 W/kg to ∼ 3300 W/kg), the root- mean-square error of the predicted vs ground truth SAR remained at 223 W/kg and 206 W/kg, with the R² scores of 0.89 and 0.85 on the testing set for 1.5 T and 3 T models, respectively. The results suggest that machine learning is a promising approach for fast assessment of RF heating of lead-like implants when only the knowledge of the lead geometry and MRI RF coil features are in hand.

Location-Dependent Spatiotemporal Antialiasing in Photoacoustic Computed Tomography

Photoacoustic computed tomography (PACT) images optical absorption contrast by detecting ultrasonic waves induced by optical energy deposition in materials such as biological tissues. An ultrasonic transducer array or its scanning equivalent is used to detect ultrasonic waves. The spatial distribution of the transducer elements must satisfy the spatial Nyquist criterion; otherwise, spatial aliasing occurs and causes artifacts in reconstructed images. The spatial Nyquist criterion poses different requirements on the transducer elements' distributions for different locations in the image domain, which has not been studied previously. In this research, we elaborate on the location dependency through spatiotemporal analysis and propose a location-dependent spatiotemporal antialiasing method. By applying this method to PACT in full-ring array geometry, we effectively mitigate aliasing artifacts with minimal effects on image resolution in both numerical simulations and in vivo experiments.

MRI of Implantation Sites Using Parallel Transmission of an Optimized Radiofrequency Excitation Vector

Postoperative care of orthopedic implants is aided by imaging to assess the healing process and the implant status. MRI of implantation sites might be compromised by radiofrequency (RF) heating and RF transmission field (B1+) inhomogeneities induced by electrically conducting implants. This study examines the applicability of safe and B1+-distortion-free MRI of implantation sites using optimized parallel RF field transmission (pTx) based on a multi-objective genetic algorithm (GA). Electromagnetic field simulations were performed for eight eight-channel RF array configurations (f = 297.2 MHz), and the most efficient array was manufactured for phantom experiments at 7.0 T. Circular polarization (CP) and orthogonal projection (OP) algorithms were applied for benchmarking the GA-based shimming. B1+ mapping and MR thermometry and imaging were performed using phantoms mimicking muscle containing conductive implants. The local SAR10g of the entire phantom in GA was 12% and 43.8% less than the CP and OP, respectively. Experimental temperature mapping using the CP yielded ΔT = 2.5-3.0 K, whereas the GA induced no extra heating. GA-based shimming eliminated B1+ artefacts at implantation sites and enabled uniform gradient-echo MRI. To conclude, parallel RF transmission with GA-based excitation vectors provides a technical foundation en route to safe and B1+-distortion-free MRI of implantation sites.

Hybrid PET/Compton-camera imaging: an imager for the next generation

Compton cameras can offer advantages over gamma cameras for some applications, since they are well suited for multitracer imaging and for imaging high-energy radiotracers, such as those employed in radionuclide therapy. While in conventional clinical settings state-of-the-art Compton cameras cannot compete with well-established methods such as PET and SPECT, there are specific scenarios in which they can constitute an advantageous alternative. The combination of PET and Compton imaging can benefit from the improved resolution and sensitivity of current PET technology and, at the same time, overcome PET limitations in the use of multiple radiotracers. Such a system can provide simultaneous assessment of different radiotracers under identical conditions and reduce errors associated with physical factors that can change between acquisitions. Advances are being made both in instrumentation developments combining PET and Compton cameras for multimodal or three-gamma imaging systems, and in image reconstruction, addressing the challenges imposed by the combination of the two modalities or the new techniques. This review article summarizes the advances made in Compton cameras for medical imaging and their combination with PET.

Application of photoacoustic computed tomography in biomedical imaging: A literature review

Photoacoustic computed tomography (PACT) is a hybrid imaging modality that combines optical excitation and acoustic detection techniques. It obtains high-resolution deep-tissue images based on the deep penetration of light, the anisotropy of light absorption in objects, and the photoacoustic effect. Hence, PACT shows great potential in biomedical sample imaging. Recently, due to its advantages of high sensitivity to optical absorption and wide scalability of spatial resolution with the desired imaging depth, PACT has received increasing attention in preclinical and clinical practice. To date, there has been a proliferation of PACT systems designed for specific biomedical imaging applications, from small animals to human organs, from ex vivo to in vivo real-time imaging, and from simple structural imaging to functional and molecular imaging with external contrast agents. Therefore, it is of great importance to summarize the previous applications of PACT systems in biomedical imaging and clinical practice. In this review, we searched for studies related to PACT imaging of biomedical tissues and samples over the past two decades; divided the studies into two categories, PACT imaging of preclinical animals and PACT imaging of human organs and body parts; and discussed the significance of the studies. Finally, we pointed out the future directions of PACT in biomedical applications. With the development of exogenous contrast agents and advances of imaging technique, in the future, PACT will enable biomedical imaging from organs to whole bodies, from superficial vasculature to internal organs, from anatomy to functions, and will play an increasingly important role in biomedical research and clinical practice.

Simulation of a novel, non-invasive radiation detector to measure the arterial input function for dynamic positron emission tomography

BACKGROUND

Dynamic positron emission tomography (dPET) is a nuclear medicine imaging technique providing functional images for organs of interest with applications in oncology, cardiology, and drug discovery. This technique requires the acquisition of the time-course arterial plasma activity concentration, called the arterial input function (AIF), which is conventionally acquired via arterial blood sampling.

PURPOSE

The aim of this study was to (A) optimize the geometry for a novel and cost efficient non-invasive detector called NID designed to measure the AIF for dPET scans through Monte Carlo simulations and (B) develop a clinical data analysis chain to successfully separate the arterial component of a simulated AIF signal from the venous component.

METHODS

The NID was optimized by using an in-house Geant4-based software package. The sensitive volume of the NID consists of a band of 10 cm long and 1 mm in diameter scintillating fibers placed over a wrist phantom. The phantom was simulated as a cylinder, 10 cm long and 6.413 cm in diameter comprised of polyethylene with two holes placed through it to simulate the patient's radial artery and vein. This phantom design was chosen to match the wrist phantom used in our previous proof of concept work. Two geometries were simulated with different arrangements of scintillating fibers. The first design used a single layer of 64 fibers. The second used two layers, an inner layer with 29 fibers and an outer layer with 30 fibers. Four positron emitting radioisotopes were simulated: ¹⁸ F, ¹¹ C, ¹⁵ O, and ⁶⁸ Ga with 100 million simulated decay events per run. The total and intrinsic efficiencies of both designs were calculated as well as the full width half maximum (FWHM) of the signal. In addition, contribution by the annihilation photons versus positrons to the signal was investigated. The results obtained from the two simulated detector models were compared. A clinical data analysis chain using an expectation maximization maximum likelihood algorithm was tested. This analysis chain will be used to separate arterial counts from the total signal.

RESULTS

The second NID design with two layers of scintillating fibers had a higher efficiency for all simulations with a maximum increase of 17% total efficiency for ¹¹ C simulation. All simulations had a significant annihilation photon contribution. The signal for ¹⁸ F and ¹¹ C was almost entirely due to photons. The clinical data analysis chain was within 1% of the true value for 434 out of 440 trials. Further experimental studies to validate these simulations will be required.

CONCLUSIONS

The design of the NID was optimized and its efficiency increased through Monte Carlo simulations. A clinical data analysis chain was successfully developed to separate the arterial component of an AIF signal from the venous component. The simulations show that the NID can be used to accurately measure the AIF non-invasively for dPET scans.

8-channel Tx dipole and 20-channel Rx loop coil array for MRI of the cervical spinal cord at 7 Tesla

The quality of cervical spinal cord images can be improved by the use of tailored radiofrequency coil solutions for ultra-high field imaging; however, very few commercial and research 7 Tesla radiofrequency coils currently exist for the spinal cord, and in particular those with parallel transmit capabilities. This work presents the design, testing and validation of a pTx/Rx coil for the human neck and cervical/upper-thoracic spinal cord. The pTx portion is composed of 8 dipoles to ensure high homogeneity over this large region of the spinal cord. The Rx portion is made of 20 semi-adaptable overlapping loops to produce high Signal-to-noise ratio (SNR) across the patient population. The coil housing is designed to facilitate patient positioning and comfort, while being tight fitting to ensure high sensitivity. We demonstrate RF shimming capabilities to optimize B ₁ ⁺ uniformity, power efficiency and/or specific absorption rate (SAR) efficiency. B ₁ ⁺ homogeneity, SNR and g-factor was evaluated in adult volunteers and demonstrated excellent performance from the occipital lobe down to the T4-T5 level. We compared the proposed coil with two state-of-the-art head and head/neck coils, confirming its superiority in the cervical and upper-thoracic regions of the spinal cord. This coil solution therefore provides a convincing platform for producing the high image quality necessary for clinical and research scanning of the upper spinal cord.

Simulation-to-real generalization for deep-learning-based refraction-corrected ultrasound tomography image reconstruction

Objective. The image reconstruction of ultrasound computed tomography is computationally expensive with conventional iterative methods. The fully learned direct deep learning reconstruction is promising to speed up image reconstruction significantly. However, for direct reconstruction from measurement data, due to the lack of real labeled data, the neural network is usually trained on a simulation dataset and shows poor performance on real data because of the simulation-to-real gap.Approach. To improve the simulation-to-real generalization of neural networks, a series of strategies are developed including a Fourier-transform-integrated neural network, measurement-domain data augmentation methods, and a self-supervised-learning-based patch-wise preprocessing neural network. Our strategies are evaluated on both the simulation dataset and real measurement datasets from two different prototype machines.Main results. The experimental results show that our deep learning methods help to improve the neural networks' robustness against noise and the generalizability to real measurement data.Significance. Our methods prove that it is possible for neural networks to achieve superior performance to traditional iterative reconstruction algorithms in imaging quality and allow for real-time 2D-image reconstruction. This study helps pave the path for the application of deep learning methods to practical ultrasound tomography image reconstruction based on simulation datasets.

Ultra-high field MRI: parallel-transmit arrays and RF pulse design

This paper reviews the field of multiple or parallel radiofrequency (RF) transmission for magnetic resonance imaging (MRI). Currently the use of ultra-high field (UHF) MRI at 7 tesla and above is gaining popularity, yet faces challenges with non-uniformity of the RF field and higher RF power deposition. Since its introduction in the early 2000s, parallel transmission (pTx) has been recognized as a powerful tool for accelerating spatially selective RF pulses and combating the challenges associated with RF inhomogeneity at UHF. We provide a survey of the types of dedicated RF coils used commonly for pTx and the important modeling of the coil behavior by electromagnetic (EM) field simulations. We also discuss the additional safety considerations involved with pTx such as the specific absorption rate (SAR) and how to manage them. We then describe the application of pTx with RF pulse design, including a practical guide to popular methods. Finally, we conclude with a description of the current and future prospects for pTx, particularly its potential for routine clinical use.

Recent Advances in Contrast-Enhanced Photoacoustic Imaging: Overcoming the Physical and Practical Challenges

For decades now, photoacoustic imaging (PAI) has been investigated to realize its potential as a niche biomedical imaging modality. Despite its highly desirable optical contrast and ultrasonic spatiotemporal resolution, PAI is challenged by such physical limitations as a low signal-to-noise ratio (SNR), diminished image contrast due to strong optical attenuation, and a lower-bound on spatial resolution in deep tissue. In addition, contrast-enhanced PAI has faced practical limitations such as insufficient cell-specific targeting due to low delivery efficiency and difficulties in developing clinically translatable agents. Identifying these limitations is essential to the continuing expansion of the field, and substantial advances in developing contrast-enhancing agents, complemented by high-performance image acquisition systems, have synergistically dealt with the challenges of conventional PAI. This review covers the past four years of research on pushing the physical and practical challenges of PAI in terms of SNR/contrast, spatial resolution, targeted delivery, and clinical application. Promising strategies for dealing with each challenge are reviewed in detail, and future research directions for next generation contrast-enhanced PAI are discussed.

Fast 3D kernel computation method for positron range correction in PET

Objective. The positron range is a fundamental, detector-independent physical limitation to spatial resolution in positron emission tomography (PET) as it causes a significant blurring of underlying activity distribution in the reconstructed images. A major challenge for positron range correction methods is to provide accurate range kernels that inherently incorporate the generally inhomogeneous stopping power, especially at tissue boundaries. In this work, we propose a novel approach to generate accurate three-dimensional (3D) blurring kernels both in homogenous and heterogeneous media to improve PET spatial resolution.Approach. In the proposed approach, positron energy deposition was approximately tracked along straight paths, depending on the positron stopping power of the underlying material. The positron stopping power was derived from the attenuation coefficient of 511 keV gamma photons according to the available PET attenuation maps. Thus, the history of energy deposition is taken into account within the range of kernels. Special emphasis was placed on facilitating the very fast computation of the positron annihilation probability in each voxel.Results. Positron path distributions of¹⁸F in low-density polyurethane were in high agreement with Geant4 simulation at an annihilation probability larger than 10^-2∼ 10^-3of the maximum annihilation probability. The Geant4 simulation was further validated with measured¹⁸F depth profiles in these polyurethane phantoms. The tissue boundary of water with cortical bone and lung was correctly modeled. Residual artifacts from the numerical computations were in the range of 1%. The calculated annihilation probability in voxels shows an overall difference of less than 20% compared to the Geant4 simulation.Significance. The proposed method is expected to significantly improve spatial resolution for non-standard isotopes by providing sufficiently accurate range kernels, even in the case of significant tissue inhomogeneities.

Reconstruction of thermoacoustic emission sources induced by proton irradiation using numerical time reversal

Objective.Mapping of dose delivery in proton beam therapy can potentially be performed by analyzing thermoacoustic emissions measured by ultrasound arrays. Here, a method is derived and demonstrated for spatial mapping of thermoacoustic sources using numerical time reversal, simulating re-transmission of measured emissions into the medium.Approach.Spatial distributions of thermoacoustic emission sources are shown to be approximated by the analytic-signal form of the time-reversed acoustic field, evaluated at the time of the initial proton pulse. Given calibration of the array sensitivity and knowledge of tissue properties, this approach approximately reconstructs the acoustic source amplitude, equal to the product of the time derivative of the radiation dose rate, mass density, and Grüneisen parameter. This approach was implemented using two models for acoustic fields of the array elements, one modeling elements as line sources and the other as rectangular radiators. Thermoacoustic source reconstructions employed previously reported measurements of emissions from proton energy deposition in tissue-mimicking phantoms. For a phantom incorporating a bone layer, reconstructions accounted for the higher sound speed in bone. Dependence of reconstruction quality on array aperture size and signal-to-noise ratio was consistent with previous acoustic simulation studies.Main results.Thermoacoustic source distributions were successfully reconstructed from acoustic emissions measured by a linear ultrasound array. Spatial resolution of reconstructions was significantly improved in the azimuthal (array) direction by incorporation of array element diffraction. Source localization agreed well with Monte Carlo simulations of energy deposition, and was improved by incorporating effects of inhomogeneous sound speed.Significance.The presented numerical time reversal approach reconstructs thermoacoustic sources from proton beam radiation, based on straightforward processing of acoustic emissions measured by ultrasound arrays. This approach may be useful for ranging and dosimetry of clinical proton beams, if acoustic emissions of sufficient amplitude and bandwidth can be generated by therapeutic proton sources.

Technical opportunities and challenges in developing total-body PET scanners for mice and rats

Positron emission tomography (PET) is the most sensitive in vivo molecular imaging technique available. Small animal PET has been widely used in studying pharmaceutical biodistribution and disease progression over time by imaging a wide range of biological processes. However, it remains true that almost all small animal PET studies using mouse or rat as preclinical models are either limited by the spatial resolution or the sensitivity (especially for dynamic studies), or both, reducing the quantitative accuracy and quantitative precision of the results. Total-body small animal PET scanners, which have axial lengths longer than the nose-to-anus length of the mouse/rat and can provide high sensitivity across the entire body of mouse/rat, can realize new opportunities for small animal PET. This article aims to discuss the technical opportunities and challenges in developing total-body small animal PET scanners for mice and rats.

Optimization of reconstruction time of ultrasound computed tomography with a piecewise homogeneous region-based refract-ray model

In this article, a novel ultrasound computed tomography (USCT) reconstruction algorithm for breast imaging is proposed. This algorithm is based on an ultrasound propagation model, the refract-ray model (RRM). In this model, the field of imaging is assumed as piecewise homogenous and is divided into several regions. The ultrasound propagation paths are considered polylines that only refract at the borders of the regions. The edge information is provided by B-mode imaging. Both simulations and experiments are implemented to validate the proposed algorithm. Compared with the traditional bent-ray model (BRM), the time of reconstructions using RRM decreases by over 90 %. In simulations, the imaging qualities for RRM and BRM are comparable, in terms of the root mean square error, the Tenengrad value, and the deformation of digital phantom. In the experiments, a cylindrical agar phantom is imaged using a customized imaging system. When imaging using RRM, the estimate of the phantom radius is about 0.1 mm in error, while it is about 0.3 mm in error using BRM. Moreover, the Tenengrad value of the result using RRM is much higher than that using BRM (9.76 compared to 0.79). The results show that the proposed algorithm can better delineate the phantom within a water bath. In future work, further experimental work is required to validate the method for improving imaging quality under breast-mimicking imaging conditions.

Non-Invasive 3D Photoacoustic Tomography of Angiographic Anatomy and Hemodynamics of Fatty Livers in Rats

Non-alcoholic fatty liver disease is the most common liver disorder worldwide, which strongly correlates to obesity, diabetes, and metabolic syndromes. Complementary to mainstream liver diagnostic modalities, photoacoustic tomography (PAT) can provide high-speed images with functional optical contrast. However, PAT has not been demonstrated to study fatty liver anatomy with clear volumetric vasculatures. The livers of multiple rats are non-invasively imaged in vivo using the recently developed 3D PAT platform. The system provides isotropically high spatial resolution in 3D space, presenting clear angiographic structures of rat livers without injecting contrast agents. Furthermore, to quantitatively analyze the difference between the livers of lean and obese rats, the authors measured several PAT features and statistical differences between the two groups are observed. In addition to the anatomy, a time-gated strategy is applied to correct respiration-induced motion artifacts and extracted the hemodynamics of major blood vessels during the breathing cycles. This study demonstrates the capabilities of 3D-PAT to reveal both angiographic anatomy and function in rat livers, providing hematogenous information for fatty liver diagnosis. 3D-PAT, as a new tool for preclinical research, warrants further improvements to be transferred to human pediatric liver imaging.

Switchable preamplifier for dual modal photoacoustic and ultrasound imaging

Photoacoustic (PA) imaging is a high-fidelity biomedical imaging technique based on the principle of molecular-specific optical absorption of biological tissue constitute. Because PA imaging shares the same basic principle as that of ultrasound (US) imaging, the use of PA/US dual-modal imaging can be achieved using a single system. However, because PA imaging is limited to a shallower depth than US imaging due to the optical extinction in biological tissue, the PA signal yields a lower signal-to-noise ratio (SNR) than US images. To selectively amplify the PA signal, we propose a switchable preamplifier for acoustic-resolution PA microscopy implemented on an application-specific integrated circuit. Using the preamplifier, we measured the increments in the SNR with both carbon lead and wire phantoms. Furthermore, in vivo whole-body PA/US imaging of a mouse with a preamplifier showed enhancement of SNR in deep tissues, unveiling deeply located organs and vascular networks. By selectively amplifying the PA signal range to a level similar to that of the US signal without contrast agent administration, our switchable amplifier strengthens the mutual complement between PA/US imaging. PA/US imaging is impending toward clinical translation, and we anticipate that this study will help mitigate the imbalance of image depth between the two imaging modalities.

Ionoacoustic application of an optical hydrophone to detect proton beam range in water

BACKGROUND

Proton range uncertainty has been the main factor limiting the ability of proton therapy to concentrate doses to tumors to their full potential. Ionoacoustic (IA) range verification is an approach to reducing this uncertainty by detecting thermoacoustic waves emitted from an irradiated volume immediately following a pulsed proton beam delivery; however, the signal weakness has been an obstacle to its clinical application. To increase the signal-to-noise ratio (SNR) with the conventional piezoelectric hydrophone (PH), the detector-sensitive volume needs to be large, but it could narrow the range of available beam angles and disturb real-time images obtained during beam delivery.

PURPOSE

To prevent this issue, we investigated a millimeter-sized optical hydrophone (OH) that exploits the laser interferometric principle. For two types of IA waves [γ-wave emitted from the Bragg peak (BP) and a spherical IA wave with resonant frequency (SPIRE) emitted from the gold fiducial marker (GM)], comparisons were made with PH in terms of waveforms, SNR, range detection accuracy, and signal intensity robustness against the small detector misalignment, particularly for SPIRE.

METHODS

A 100-MeV proton beam with a 27 ns pulse width and 4 mm beam size was produced using a fixed-field alternating gradient accelerator and was irradiated to the water phantom. The GM was set on the beam's central axis. Acrylic plates of various thicknesses, up to 12 mm, were set in front of the phantoms to shift the proton range. OH was set distal and lateral to the beam, and the range was estimated using the time-of-flight method for γ-wave and by comparing with the calibration data (SPIRE intensity versus the distance between the GM and BP) derived from an IA wave transport simulation for SPIRE. The BP dose per pulse was 0.5-0.6 Gy. To measure the variation in SPIRE amplitude against the hydrophone misalignment, the hydrophone was shifted by ± 2 mm at a maximum in lateral directions.

RESULTS

Despite its small size, OH could detect γ-wave with a higher SNR than the conventional PH (diameter, 29 mm), and a single measurement was sufficient to detect the beam range with a submillimeter accuracy in water. In the SPIRE measurement, OH was far more robust against the detector misalignment than the focused PH (FPH) used in our previous study [5%/mm (OH) versus 80%/mm (FPH)], and the correlation between the measured SPIRE intensity and the distance between the GM and BP agreed well with the simulation results. However, the OH sensitivity was lower than the FPH sensitivity, and about 5.6-Gy dose was required to decrease the intensity variation among measurements to less than 10%.

CONCLUSION

The miniature OH was found to detect weak IA signals produced by proton beams with a BP dose used in hypofractionated regimens. The OH sensitivity improvement at the MHz regime is worth exploring as the next step.

A comparative study of RF heating of deep brain stimulation devices in vertical vs. horizontal MRI systems

The majority of studies that assess magnetic resonance imaging (MRI) induced radiofrequency (RF) heating of the tissue when active electronic implants are present have been performed in horizontal, closed-bore MRI systems. Vertical, open-bore MRI systems have a 90° rotated magnet and a fundamentally different RF coil geometry, thus generating a substantially different RF field distribution inside the body. Little is known about the RF heating of elongated implants such as deep brain stimulation (DBS) devices in this class of scanners. Here, we conducted the first large-scale experimental study investigating whether RF heating was significantly different in a 1.2 T vertical field MRI scanner (Oasis, Fujifilm Healthcare) compared to a 1.5 T horizontal field MRI scanner (Aera, Siemens Healthineers). A commercial DBS device mimicking 30 realistic patient-derived lead trajectories extracted from postoperative computed tomography images of patients who underwent DBS surgery at our institution was implanted in a multi-material, anthropomorphic phantom. RF heating around the DBS lead was measured during four minutes of high-SAR RF exposure. Additionally, we performed electromagnetic simulations with leads of various internal structures to examine this effect on RF heating. When controlling for RMS B1+, the temperature increase around the DBS lead-tip was significantly lower in the vertical scanner compared to the horizontal scanner (0.33 ± 0.24°C vs. 4.19 ± 2.29°C). Electromagnetic simulations demonstrated up to a 17-fold reduction in the maximum of 0.1g-averaged SAR in the tissue surrounding the lead-tip in the vertical scanner compared to the horizontal scanner. Results were consistent across leads with straight and helical internal wires. Radiofrequency heating and power deposition around the DBS lead-tip were substantially lower in the 1.2 T vertical scanner compared to the 1.5 T horizontal scanner. Simulations with different lead structures suggest that the results may extend to leads from other manufacturers.

Amiodarone attenuates cardiac Rubidium-82 in consecutive PET/CT scans in a rodent model

BACKGROUND

Risk stratification and diagnosis using Rubidium-82 (82Rb) positron emission tomography (PET) is a routine clinical approach in coronary artery disease (CAD). Various drugs are used to treat CAD; however, whether any of them change the uptake of 82Rb in the heart has not been investigated. The aim of this study is to determine whether drugs used in treatment of CAD affect the uptake of 82Rb in the heart in healthy rats.

METHODS

Seventy-seven Sprague-Dawley rats were included in the cross-sectional study. All rats underwent baseline 82Rb PET/CT and divided into eleven groups treated with different drugs. One group was control group (no treatment), eight groups were treated with monotherapy (amiodarone, acetylsalicylic acid (ASA), clopidogrel, ticagrelor, atorvastatin, enalapril, amlodipine, metoprolol succinate), and two groups were treated with polypharmacy (ASA, ticagrelor, atorvastatin, amlodipine or ASA, clopidogrel, atorvastatin, amlodipine). Once a day, they were administered pharmacological therapy through oral gavage, and on day seven, follow-up scanned with 82Rb PET/CT.

RESULTS

In the control group without pharmacological treatment, no difference in the standard uptake value (SUV) ratio between heart and muscle from baseline to follow-up (5.8 vs 7.0, P = .3) was found. The group treated with amiodarone had a significantly reduced SUV ratio from baseline to follow-up (5.8 vs 5.1, P = .008). All other drugs investigated had no difference in SUV ratio from baseline to follow-up.

CONCLUSION

In this study, we showed that drugs normally used to treat CAD do not affect the uptake of 82Rb. However, amiodarone result in a significantly lowered 82Rb uptake, compared to control. This information about amiodarone would probably not change the size assessment of a myocardial perfusion defect in a clinical setting. However, it could change the kinetic parameters when assessing absolute myocardial blood flow in patients treated with amiodarone.

Implementation of pencil beam redefinition algorithm (PBRA) for intraoperative electron radiation therapy (IOERT) treatment planning

PURPOSE

Similar to other radiation therapy techniques, intraoperative electron radiation therapy (IOERT) can benefit from an online treatment planning system (TPS). Among all the analytical electron dose calculation algorithms, pencil beam redefinition algorithm (PBRA) has shown an acceptable accuracy in inhomogeneities. The input dataset for PBRA includes electron planar fluence, mean direction and root mean square (RMS) spread about the mean direction which had been introduced based on the conventional linear accelerator geometry in former studies. Herein, three methods for implementing PBRA for IOERT system are presented.

METHODS

The initialization parameters were identified using Monte Carlo (MC) simulation of a dedicated IOERT system equipped with a cylindrical 10 cm applicator, irradiating a water phantom. Phase space distribution of electrons was recorded on a plane below the applicator. The input dataset was extracted for 2 × 2 mm² pixels and energy bin width of 1 MeV.

RESULTS

PBRA was implemented with three initialization methods and compared to MC. The 3D gamma analysis of the algorithm with the Formula method, which was in best agreement with MC in a simple water phantom, showed passing rates of more than 99 % for all nominal energies and it was 97.1 % for 8 MeV in the presence of protecting disk and irregular surface. Implementing PBRA on CUDA C++ resulted in 5 s run time for 8 MeV nominal energy in a water phantom.

CONCLUSIONS

The agreement between PBRA dose calculation and MC is promising for the development of an intraoperative TPS for IOERT.

Hybrid magnetic resonance and optoacoustic tomography (MROT) for preclinical neuroimaging

Multi-modal imaging is essential for advancing our understanding of brain function and unraveling pathophysiological processes underlying neurological and psychiatric disorders. Magnetic resonance (MR) and optoacoustic (OA) imaging have been shown to provide highly complementary contrasts and capabilities for preclinical neuroimaging. True integration between these modalities can thus offer unprecedented capabilities for studying the rodent brain in action. We report on a hybrid magnetic resonance and optoacoustic tomography (MROT) system for concurrent noninvasive structural and functional imaging of the mouse brain. Volumetric OA tomography was designed as an insert into a high-field MR scanner by integrating a customized MR-compatible spherical transducer array, an illumination module, and a dedicated radiofrequency coil. A tailored data processing pipeline has been developed to mitigate signal crosstalk and accurately register image volumes acquired with T1-weighted, angiography, and blood oxygenation level-dependent (BOLD) sequences onto the corresponding vascular and oxygenation data recorded with the OA modality. We demonstrate the concurrent acquisition of dual-mode anatomical and angiographic brain images with the scanner, as well as real-time functional readings of multiple hemodynamic parameters from animals subjected to oxygenation stress. Our approach combines the functional and molecular imaging advantages of OA with the superb soft-tissue contrast of MR, further providing an excellent platform for cross-validation of functional readings by the two modalities.

Proton beam range verification by means of ionoacoustic measurements at clinically relevant doses using a correlation-based evaluation

Purpose

The Bragg peak located at the end of the ion beam range is one of the main advantages of ion beam therapy compared to X-Ray radiotherapy. However, verifying the exact position of the Bragg peak within the patient online is a major challenge. The goal of this work was to achieve submillimeter proton beam range verification for pulsed proton beams of an energy of up to 220 MeV using ionoacoustics for a clinically relevant dose deposition of typically 2 Gy per fraction by i) using optimal proton beam characteristics for ionoacoustic signal generation and ii) improved signal detection by correlating the signal with simulated filter templates.

Methods

A water tank was irradiated with a preclinical 20 MeV proton beam using different pulse durations ranging from 50 ns up to 1 μs in order to maximise the signal-to-noise ratio (SNR) of ionoacoustic signals. The ionoacoustic signals were measured using a piezo-electric ultrasound transducer in the MHz frequency range. The signals were filtered using a cross correlation-based signal processing algorithm utilizing simulated templates, which enhances the SNR of the recorded signals. The range of the protons is evaluated by extracting the time of flight (ToF) of the ionoacoustic signals and compared to simulations from a Monte Carlo dose engine (FLUKA).

Results

Optimised SNR of 28.0 ± 10.6 is obtained at a beam current of 4.5 μA and a pulse duration of 130 ns at a total peak dose deposition of 0.5 Gy. Evaluated ranges coincide with Monte Carlo simulations better than 0.1 mm at an absolute range of 4.21 mm. Higher beam energies require longer proton pulse durations for optimised signal generation. Using the correlation-based post-processing filter a SNR of 17.8 ± 5.5 is obtained for 220 MeV protons at a total peak dose deposition of 1.3 Gy. For this clinically relevant dose deposition and proton beam energy, submillimeter range verification was achieved at an absolute range of 303 mm in water.

Conclusion

Optimal proton pulse durations ensure an ideal trade-off between maximising the ionoacoustic amplitude and minimising dose deposition. In combination with a correlation-based post-processing evaluation algorithm, a reasonable SNR can be achieved at low dose levels putting clinical applications for online proton or ion beam range verification into reach.

3Din vivodose verification in prostate proton therapy with deep learning-based proton-acoustic imaging

Dose delivery uncertainty is a major concern in proton therapy, adversely affecting the treatment precision and outcome. Recently, a promising technique, proton-acoustic (PA) imaging, has been developed to provide real-timein vivo3D dose verification. However, its dosimetry accuracy is limited due to the limited-angle view of the ultrasound transducer. In this study, we developed a deep learning-based method to address the limited-view issue in the PA reconstruction. A deep cascaded convolutional neural network (DC-CNN) was proposed to reconstruct 3D high-quality radiation-induced pressures using PA signals detected by a matrix array, and then derive precise 3D dosimetry from pressures for dose verification in proton therapy. To validate its performance, we collected 81 prostate cancer patients' proton therapy treatment plans. Dose was calculated using the commercial software RayStation and was normalized to the maximum dose. The PA simulation was performed using the open-source k-wave package. A matrix ultrasound array with 64 × 64 sensors and 500 kHz central frequency was simulated near the perineum to acquire radiofrequency (RF) signals during dose delivery. For realistic acoustic simulations, tissue heterogeneity and attenuation were considered, and Gaussian white noise was added to the acquired RF signals. The proposed DC-CNN was trained on 204 samples from 69 patients and tested on 26 samples from 12 other patients. Predicted 3D pressures and dose maps were compared against the ground truth qualitatively and quantitatively using root-mean-squared-error (RMSE), gamma-index (GI), and dice coefficient of isodose lines. Results demonstrated that the proposed method considerably improved the limited-view PA image quality, reconstructing pressures with clear and accurate structures and deriving doses with a high agreement with the ground truth. Quantitatively, the pressure accuracy achieved an RMSE of 0.061, and the dose accuracy achieved an RMSE of 0.044, GI (3%/3 mm) of 93.71%, and 90%-isodose line dice of 0.922. The proposed method demonstrates the feasibility of achieving high-quality quantitative 3D dosimetry in PA imaging using a matrix array, which potentially enables the online 3D dose verification for prostate proton therapy.

Ultrasonic Imaging and Sensors

Ultrasound imaging is a wide research field, covering areas from wave propagation physics, sensors and front-end electronics to image reconstruction algorithms and software [...].

ResViT: Residual Vision Transformers for Multimodal Medical Image Synthesis

Generative adversarial models with convolutional neural network (CNN) backbones have recently been established as state-of-the-art in numerous medical image synthesis tasks. However, CNNs are designed to perform local processing with compact filters, and this inductive bias compromises learning of contextual features. Here, we propose a novel generative adversarial approach for medical image synthesis, ResViT, that leverages the contextual sensitivity of vision transformers along with the precision of convolution operators and realism of adversarial learning. ResViT's generator employs a central bottleneck comprising novel aggregated residual transformer (ART) blocks that synergistically combine residual convolutional and transformer modules. Residual connections in ART blocks promote diversity in captured representations, while a channel compression module distills task-relevant information. A weight sharing strategy is introduced among ART blocks to mitigate computational burden. A unified implementation is introduced to avoid the need to rebuild separate synthesis models for varying source-target modality configurations. Comprehensive demonstrations are performed for synthesizing missing sequences in multi-contrast MRI, and CT images from MRI. Our results indicate superiority of ResViT against competing CNN- and transformer-based methods in terms of qualitative observations and quantitative metrics.

Expediting Image Acquisition in Reflection Ultrasound Computed Tomography

Reflection ultrasound computed tomography (RUCT) attains optimal image quality from objects that can be fully accessed from multiple directions, such as the human breast or small animals. Owing to the full-view tomography approach based on the compounding of images taken from multiple angles, RUCT effectively mitigates several deficiencies afflicting conventional pulse-echo ultrasound (US) systems, such as speckle patterns and interuser variability. On the other hand, the small interelement pitch required to fulfill the spatial sampling criterion in the circular transducer configuration used in RUCT typically implies the use of an excessive number of independent array elements. This increases the system's complexity and costs, and limits the achievable imaging speed. Here, we explore acquisition schemes that enable RUCT imaging with the reduced number of transmit/receive elements. We investigated the influence of the element size in transmission and reception in a ring array geometry. The performance of a sparse acquisition approach based on partial acquisition from a subset of the elements has been further assessed. A larger element size is shown to preserve contrast and resolution at the center of the field of view (FOV), while a reduced number of elements is shown to cause uniform loss of contrast and resolution across the entire FOV. The tradeoffs of achievable FOV, contrast-to-noise ratio, and temporal and spatial resolutions are assessed in phantoms and in vivo mouse experiments. The experimental analysis is expected to aid the development of optimized hardware and image acquisition strategies for RUCT and, thus, result in more affordable imaging systems facilitating wider adoption.

Towards an integrated neonatal brain and cardiac examination capability at 7 T: electromagnetic field simulations and early phantom experiments using an 8-channel dipole array

OBJECTIVE

Neonatal brain and cardiac imaging would benefit from the increased signal-to-noise ratio levels at 7 T compared to lower field. Optimal performance might be achieved using purpose designed RF coil arrays. In this study, we introduce an 8-channel dipole array and investigate, using simulations, its RF performances for neonatal applications at 7 T.

METHODS

The 8-channel dipole array was designed and evaluated for neonatal brain/cardiac configurations in terms of SAR efficiency (ratio between transmit-field and maximum specific-absorption-rate level) using adjusted dielectric properties for neonate. A birdcage coil operating in circularly polarized mode was simulated for comparison. Validation of the simulation model was performed on phantom for the coil array.

RESULTS

The 8-channel dipole array demonstrated up to 46% higher SAR efficiency levels compared to the birdcage coil in neonatal configurations, as the specific-absorption-rate levels were alleviated. An averaged normalized root-mean-square-error of 6.7% was found between measured and simulated transmit field maps on phantom.

CONCLUSION

The 8-channel dipole array design integrated for neonatal brain and cardiac MR was successfully demonstrated, in simulation with coverage of the baby and increased SAR efficiency levels compared to the birdcage. We conclude that the 8Tx-dipole array promises safe operating procedures for MR imaging of neonatal brain and heart at 7 T.

Deep-Learning-Based Ultrasound Sound-Speed Tomography Reconstruction with Tikhonov Pseudo-Inverse Priori

Ultrasound sound-speed tomography (USST) is a promising technology for breast imaging and breast cancer detection. Its reconstruction is a complex non-linear mapping from the projection data to the sound-speed image (SSI). The traditional reconstruction methods include mainly the ray-based methods and the waveform-based methods. The ray-based methods with linear approximation have low computational cost but low reconstruction quality; the full wave-based methods with the complex non-linear model have high quality but high cost. To achieve both high quality and low cost, we introduced traditional linear approximation as prior knowledge into a deep neural network and treated the complex non-linear mapping of USST reconstruction as a combination of linear mapping and non-linear mapping. In the proposed method, the linear mapping was seamlessly implemented with a fully connected layer and initialized using the Tikhonov pseudo-inverse matrix. The non-linear mapping was implemented using a U-shape Net (U-Net). Furthermore, we proposed the Tikhonov U-shape net (TU-Net), in which the linear mapping was done before the non-linear mapping, and the U-shape Tikhonov net (UT-Net), in which the non-linear mapping was done before the linear mapping. Moreover, we conducted simulations and experiments for evaluation. In the numerical simulation, the root-mean-squared error was 6.49 and 4.29 m/s for the UT-Net and TU-Net, the peak signal-to-noise ratio was 49.01 and 52.90 dB, the structural similarity was 0.9436 and 0.9761 and the reconstruction time was 10.8 and 11.3 ms, respectively. In this study, the SSIs obtained with the proposed methods exhibited high sound-speed accuracy. Both the UT-Net and the TU-Net achieved high quality and low computational cost.

In Vivo Measurements of the Bulk Ultrasonic Attenuation Coefficient of Breast Tissue Using a Novel Phase-Insensitive Receiver

This study describes the first in vivo acoustic attenuation measurements of breast tissue undertaken using a novel phase-insensitive detection technique employing a differential pyroelectric sensor. The operation of the sensor is thermal in nature, with its output signal being dictated by the acoustic power integrated over its surface. The particularly novel feature of the sensor lies in its differential principle of operation, which significantly enhances its immunity to background acoustic and vibration noise. A large area variant of the sensor was used to detect ultrasonic energy generated by an array of 14 discrete 3.2-MHz plane piston transducers, transmitted through pendent breasts in water. The transduction and reception capability represent key parts of a prototype Quantitative Ultrasound Computed Tomography Test Facility developed at the National Physical Laboratory to study the efficacy of phase-insensitive ultrasound computed tomography of breast phantoms containing a range of appropriate inclusions, in particular, the measurement uncertainties associated with quantitative reconstructions of the acoustic attenuation coefficient. For this study, attenuation coefficient measurements were made using 1-D projections on 12 nominally healthy study volunteers, whose age ranged from 19 to 65 years. Averaged or bulk attenuation coefficient values were generated in the range 1.7-4.6 dBcm ^-1 at 3.2 MHz and have been compared with existing literature, derived from in vivo and ex vivo studies. Results are encouraging and indicate that the relatively simple technique could be applied as a robust method for assessing the properties of breast tissue, particularly the balance of fatty (adipose) and fibroglandular components.

Positron Range Corrections and Denoising Techniques for Gallium-68 PET Imaging: A Literature Review

Gallium-68 (68Ga) is characterized by relatively high positron energy compared to Fluorine-18 (18F), causing substantial image quality degradation. Furthermore, the presence of statistical noise can further degrade image quality. The aim of this literature review is to identify the recently developed positron range correction techniques for 68Ga, as well as noise reduction methods to enhance the image quality of low count 68Ga PET imaging. The search engines PubMed and Scopus were employed, and we limited our research to published results from January 2010 until 1 August 2022. Positron range correction was achieved by using either deblurring or deep learning approaches. The proposed techniques improved the image quality and, in some cases, achieved an image quality comparable to 18F PET. However, none of these techniques was validated in clinical studies. PET denoising for 68Ga-labeled radiotracers was reported using either reconstruction-based techniques or deep learning approaches. It was demonstrated that both approaches can substantially enhance the image quality by reducing the noise levels of low count 68Ga PET imaging. The combination of 68Ga-specific positron range correction techniques and image denoising approaches may enable the application of low-count, high-quality 68Ga PET imaging in a clinical setting.

Generation of Prior Information in a Dual-Mode Microwave-Ultrasound Breast Imaging System

A new breast imaging system capable of obtaining ultrasound and microwave scattered-field measurements with minimal or no movement of the breast between measurements has recently been reported. In this work, we describe the methodology that has been developed to generate prior information about the internal structures of the breast based on ultrasound data measured with the dual-mode system. This prior information, estimating both the geometry and complex-valued permittivity of tissues within the breast, is incorporated into the microwave inversion algorithm as a means of enhancing image quality. Several techniques to map reconstructed ultrasound speed to complex-valued relative permittivity are investigated. Quantitative images of two simplified dual-mode breast phantoms obtained using experimental data and the various forms of prior information are presented. Though preliminary, the results presented herein provide an understanding of the impacts of different forms of prior information on dual-mode reconstructions of the breast and can be used to inform future work on the subject.

Comparison of image quality and spatial resolution between 18F, 68Ga, and 64Cu phantom measurements using a digital Biograph Vision PET/CT

BACKGROUND

PET nuclides can have a considerable influence on the spatial resolution and image quality of PET/CT scans, which can influence diagnostics in oncology, for example. The individual impact of the positron energy of ¹⁸F, ⁶⁸Ga, and ⁶⁴Cu on spatial resolution and image quality was compared for PET/CT scans acquired using a clinical, digital scanner.

METHODS

A Jaszczak phantom and a NEMA PET body phantom were filled with ¹⁸F-FDG, ⁶⁸Ga-HCl, or ⁶⁴Cu-HCl, and PET/CT scans were performed on a Siemens Biograph Vision. Acquired images were analyzed regarding spatial resolution and image quality (recovery coefficients (RC), coefficient of variation within the background, contrast recovery coefficient (CRC), contrast-noise ratio (CNR), and relative count error in the lung insert). Data were compared between scans with different nuclides.

RESULTS

We found that image quality was comparable between ¹⁸F-FDG and ⁶⁴Cu-HCl PET/CT measurements featuring similar maximal endpoint energies of the positrons. In comparison, RC, CRC, and CNR were degraded in ⁶⁸Ga-HCl data despite similar count rates. In particular, the two smallest spheres of 10 mm and 13 mm diameter revealed lower RC, CRC, and CNR values. The spatial resolution was similar between ¹⁸F-FDG and ⁶⁴Cu-HCl but up to 18% and 23% worse compared with PET/CT images of the NEMA PET body phantom filled with ⁶⁸Ga-HCl.

CONCLUSIONS

The positron energy of the PET nuclide influences the spatial resolution and image quality of a digital PET/CT scan. The image quality and spatial resolution of ⁶⁸Ga-HCl PET/CT images were worse than those of ¹⁸F-FDG or ⁶⁴Cu-HCl despite similar count rates.