Showing papers in "arXiv: Medical Physics in 2021"

Commissioning of a clinical pencil beam scanning proton therapy unit for ultrahigh dose rates (FLASH)

Abstract: Purpose: The purpose of this work was to provide a flexible platform for FLASH research with protons by adapting a former clinical pencil beam scanning gantry to irradiations with ultrahigh dose rates. Methods: PSI Gantry 1 treated patients until December 2018. We optimized the beamline parameters to transport the 250 MeV beam extracted from the PSI COMET accelerator to the treatment room, maximizing the transmission of beam intensity to the sample. We characterized a dose monitor on the gantry to ensure good control of the dose, delivered in spot-scanning mode. We characterized the beam for different dose rates and field sizes for transmission irradiations. We explored scanning possibilities in order to enable conformal irradiations or transmission irradiations of large targets (with transverse scanning). Results: We achieved a transmission of 86 % from the cyclotron to the treatment room. We reached a peak dose rate of 9000 Gy/s at 3 mm water equivalent depth, along the central axis of a single pencil beam. Field sizes of up to 5x5 mm$^{2}$ were achieved for single spot FLASH irradiations. Fast transverse scanning allowed to cover a field of 16x1.2 cm$^{2}$. With the use of a nozzle-mounted range shifter we are able to span depths in water ranging from 19.6 to 37.9 cm. Various dose levels were delivered with a precision within less than 1 %. Conclusions: We have realized a proton FLASH irradiation setup able to investigate continuously a wide dose rate spectrum, from 1 to 9000 Gy/s in a single spot irradiation as well as in the pencil beam scanning mode. As such, we have developed a versatile test bench for FLASH research.

Neurite Exchange Imaging (NEXI): A minimal model of diffusion in gray matter with inter-compartment water exchange

Abstract: Biophysical models of diffusion in white matter have been center-stage over the past two decades and are essentially based on what is now commonly referred to as the "Standard Model" (SM) of non-exchanging compartments with Gaussian diffusion. In this work, we focus on diffusion MRI in cortical gray matter, which requires rethinking basic microstructure modeling blocks. In particular, at least three contributions beyond the SM need to be considered for gray matter: water exchange across the cell membrane - between neurites and the extracellular space; non-Gaussian diffusion along cellular processes - resulting from structural disorder; and signal contribution from soma. For the first contribution, we propose Neurite Exchange Imaging (NEXI) as an extension of the SM of diffusion, which builds on the anisotropic Karger model of two exchanging compartments. Using multi-shell multi-t datasets acquired in the rat brain in vivo, we investigate the suitability of NEXI to describe the diffusion signal in the cortex, compared to the other two possible contributions. In the light of minimal diffusion time-dependence, kurtosis decaying with diffusion time as approximately 1/t and lower signal for longer diffusion times at high b, we identify exchange as the mechanism that best explains these signal signatures in both low-b and high-b regime, and thereby propose NEXI as the minimal model for cortex microstructure mapping. We finally highlight multi-shell multi-t acquisitions protocols as being best suited to estimate NEXI model parameters reliably. Using this approach, we estimate the inter-compartment water exchange time to be 20 - 60 ms in the rat cortex in vivo, which is of the same order or shorter than the diffusion time in typical diffusion MRI acquisitions.

Low Temperature Plasma for Biology, Hygiene, and Medicine: Perspective and Roadmap

Abstract: Plasma, the fourth and most pervasive state of matter in the visible universe, is a fascinating medium that is connected to the beginning of our universe itself. Man-made plasmas are at the core of many technological advances that include the fabrication of semiconductor devices, which enabled the modern computer and communication revolutions. The introduction of low temperature, atmospheric pressure plasmas to the biomedical field has ushered a new revolution in the healthcare arena that promises to introduce plasma-based therapies to combat some thorny and long-standing medical challenges. This paper presents an overview of where research is at today and discusses innovative concepts and approaches to overcome present challenges and take the field to the next level. It is written by a team of experts who took an in-depth look at the various biomedical applications, made critical analysis, and proposed ideas and concepts that should help the research community focus their efforts on clear and practical steps necessary to keep the field advancing for decades to come.

Diffusion time dependence, power-law scaling, and exchange in gray matter

Abstract: Diffusion MRI (dMRI) provides contrast that reflect diffusing spins' interactions with microstructural features of biological systems, but its specificity remains limited due to the ambiguity of its relation to the underlying microstructure. To improve specificity, biophysical models of white matter (WM) typically express dMRI signals according to the Standard Model (SM) and have more recently in gray matter (GM) attempted to incorporate cell soma (the SANDI model). The validity of the assumptions underlying these models, however, remains largely undetermined, especially in GM. Observing the models' unique, functional properties, such as the $b^{-1/2}$ power-law associated with 1d diffusion, has emerged as a fruitful strategy for experimental validation. The absence of this signature in GM has been explained by neurite water exchange, non-linear morphology, and/or obscuring soma signal contributions. Here, we present simulations in realistic neurons demonstrating that curvature and branching does not destroy the stick power-law in impermeable neurites, but that their signal is drowned by the soma under typical experimental conditions: Nevertheless, we identify an attainable experimental regime in which the neurite signal dominates. Furthermore, we find that exchange-driven time dependence produces a behavior opposite to that expected from restricted diffusion, thereby providing a functional signature disambiguating the two effects. We present data from dMRI experiments in ex vivo rat brain at ultrahigh field and observe a time dependence consistent with substantial exchange and a GM stick power-law. The first finding suggests significant water exchange while the second suggests a small sub-population of impermeable neurites. To quantify our observations, we harness the Karger exchange model and incorporate the corresponding signal time dependence in SM and SANDI.

Deep Learning for fully automatic detection, segmentation, and Gleason Grade estimation of prostate cancer in multiparametric Magnetic Resonance Images.

Abstract: The emergence of multi-parametric magnetic resonance imaging (mpMRI) has had a profound impact on the diagnosis of prostate cancers (PCa), which is the most prevalent malignancy in males in the western world, enabling a better selection of patients for confirmation biopsy. However, analyzing these images is complex even for experts, hence opening an opportunity for computer-aided diagnosis systems to seize. This paper proposes a fully automatic system based on Deep Learning that takes a prostate mpMRI from a PCa-suspect patient and, by leveraging the Retina U-Net detection framework, locates PCa lesions, segments them, and predicts their most likely Gleason grade group (GGG). It uses 490 mpMRIs for training/validation, and 75 patients for testing from two different datasets: ProstateX and IVO (Valencia Oncology Institute Foundation). In the test set, it achieves an excellent lesion-level AUC/sensitivity/specificity for the GGG$\geq$2 significance criterion of 0.96/1.00/0.79 for the ProstateX dataset, and 0.95/1.00/0.80 for the IVO dataset. Evaluated at a patient level, the results are 0.87/1.00/0.375 in ProstateX, and 0.91/1.00/0.762 in IVO. Furthermore, on the online ProstateX grand challenge, the model obtained an AUC of 0.85 (0.87 when trained only on the ProstateX data, tying up with the original winner of the challenge). For expert comparison, IVO radiologist's PI-RADS 4 sensitivity/specificity were 0.88/0.56 at a lesion level, and 0.85/0.58 at a patient level. Additional subsystems for automatic prostate zonal segmentation and mpMRI non-rigid sequence registration were also employed to produce the final fully automated system. The code for the ProstateX-trained system has been made openly available at this https URL. We hope that this will represent a landmark for future research to use, compare and improve upon.

Respiratory aerosols and droplets in the transmission of infectious diseases

Abstract: Knowing the physicochemical properties of exhaled droplets and aerosol particles is a prerequisite for a detailed mechanistic understanding and effective prevention of the airborne transmission of infectious human diseases. This article provides a critical review and synthesis of scientific knowledge on the number concentrations, size distributions, composition, mixing state, and related properties of respiratory particles emitted upon breathing, speaking, singing, coughing, and sneezing. We derive and present a parameterization of respiratory particle size distributions based on five lognormal modes related to different origins in the respiratory tract, which can be used to trace and localize the sources of infectious particles. This approach may support the medical treatment as well as the risk assessment for aerosol and droplet transmission of infectious diseases. It was applied to analyze which respiratory activities may drive the spread of specific pathogens, such as Mycobacterium tuberculosis, influenza viruses, and SARS-CoV-2 viruses. The results confirm the high relevance of vocalization for the transmission of SARS-CoV-2 as well as the usefulness of physical distancing, face masks, room ventilation, and air filtration as preventive measures against COVID-19 and other airborne infectious diseases.

Time Window of Head Impact Kinematics Measurement for Calculation of Brain Strain and Strain Rate in American Football

Abstract: Wearable devices have been shown to effectively measure the head movement during impacts in sports like American football. When a head impact occurs, the device is triggered to collect and save the kinematic measurements during a predefined time window. Then, based on the collected kinematics, finite element (FE) head models can calculate brain strain, which is used to evaluate the risk of mild traumatic brain injury. To find a time window that can provide a sufficient duration of kinematics for FE analysis, we investigated 118 on-field video-confirmed head impacts collected by the Stanford Instrumented Mouthguard. Because the individual differences in brain geometry influence these calculations, we included six representative brain geometries and found that larger brains need a longer time window of kinematics for accurate calculation. Among the different sizes of brains, a pre-trigger time of 20 ms and a post-trigger time of 70 ms were found to yield calculations of brain strain and strain rate that were not significantly different from calculations using the original 200 ms time window recorded by the mouthguard.

Link between SARS-CoV-2 emissions and airborne concentrations: closing the gap in understanding

Abstract: The question of how SARS-CoV-2 is transmitted remains surprisingly controversial today, especially with reference to airborne transmission. In fact, despite a large body of scientific evidence, health and regulatory authorities still require direct proof of this mode of transmission. To close this gap, we measured the saliva viral load of SARS-CoV-2 of an infected subject located in a hospital room, as well as the airborne SARS-CoV-2 concentration in the room resulting from the person breathing and speaking. As the next step, we simulated the same scenarios to estimate the concentration of RNA copies in the air through a novel predictive theoretical approach. Finally, we conducted a comparative analysis (i.e. a metrological compatibility analysis) of the differences between the experimental and theoretical results by estimating the uncertainties of these two approaches. Our results showed that for an infected subject's saliva load ranging between 2.4x106 and 5.5x106 RNA copies mL-1, the corresponding airborne SARS-CoV-2 concentration was not detectable when the person was breathing, but was 16.1 (with an uncertainty of 32.8) RNA copies m-3 when speaking. The application of the novel predictive estimation approach provided average concentrations of 3.2 (uncertainty range of 0.2-8.3) and 18.5 (uncertainty range of 4.5-43.0) RNA copies m-3 for breathing and speaking scenarios, respectively, thus confirming that for the breathing scenario, the airborne RNA concentration would be undetectable, being below the minimum detection threshold of the experimental apparatus (< 2 RNA copies m-3).

Dual-energy CT imaging with limited-angular-range data

Abstract: In dual-energy computed tomography (DECT), low- and high- kVp data are collected often over a full-angular range (FAR) of $360^\circ$. While there exists strong interest in DECT with low- and high-kVp data acquired over limited-angular ranges (LARs), there remains little investigation of image reconstruction in DECT with LAR data. Objective: We investigate image reconstruction with minimized LAR artifacts from low- and high-kVp data over LARs of $\le 180^\circ$ by using a directional-total-variation (DTV) algorithm. Methods: Image reconstruction from LAR data is formulated as a convex optimization problem in which data-$\ell_2$ is minimized with constraints on image's DTVs along orthogonal axes. We then achieve image reconstruction by applying the DTV algorithm to solve the optimization problem. We conduct numerical studies from data generated over arcs of LARs, ranging from $14^\circ$ to $180^\circ$, and perform visual inspection and quantitative analysis of images reconstructed. Results: Monochromatic images of interest obtained with the DTV algorithm from LAR data show substantially reduced artifacts that are observed often in images obtained with existing algorithms. The improved image quality also leads to accurate estimation of physical quantities of interest, such as effective atomic number and iodine-contrast concentration. Conclusion: Our study reveals that from LAR data of low- and high-kVp, monochromatic images can be obtained that are visually, and physical quantities can be estimated that are quantitatively, comparable to those obtained in FAR DECT. Significance: As LAR DECT is of high practical application interest, the results acquired in the work may engender insights into the design of DECT with LAR scanning configurations of practical application significance.

In Vivo Assessment of Hypertensive Nephrosclerosis Using Ultrasound Localization Microscopy.

Abstract: Purpose: As a typical chronic kidney disease (CKD), hypertensive nephrosclerosis (HN) is a common syndrome of hypertension, characterized by chronic kidney microvascular damage. Early diagnosis of microvascular damage using conventional ultrasound imaging encounters challenges in sensitivity and specificity owing to the inherent diffraction limit. Ultrasound localization microscopy (ULM) has been developed to obtain microvasculature and microvascular hemodynamics within the kidney, and would be a promising tool for early diagnosis of CKD. Methods: In this study, the advantage of ULM over conventional clinical inspection (serum and urine tests) and ultrasound imaging (Doppler and contrast-enhanced ultrasound imaging) for early diagnosis of HN was investigated. Examinations were carried out on 6 spontaneously hypertensive rats (SHR) and 5 normal Wistar-Kyoto (WKY) rats at the age of 10 weeks. Results: The experimental results showed that the indicators derived from conventional clinical inspection and ultrasound imaging (PSV, EDV, RI, RT, IMAX, mTT and AUC) did not show significant difference between hypertensive and healthy rats (p > 0.05), while the TTP derived from CEUS (p < 0.05) and the mean blood flow speed in artery of SHR derived from ULM is significantly higher than that of WKY rats (p < 0.01). Conclusion: The quantitative results showed that ULM has higher sensitivity than conventional clinical inspection and ultrasound imaging. ULM may promise a reliable solution for early diagnosis of HN.

Optical and magnetic resonance imaging approaches for investigating the tumour microenvironment: state-of-the-art review and future trends

Abstract: The tumour microenvironment (TME) strongly influences tumorigenesis and metastasis. Two of the most characterized properties of the TME are acidosis and hypoxia, both of which are considered hallmarks of tumours as well as critical factors in response to anticancer treatments. Currently, various imaging approaches exist to measure acidosis and hypoxia in the TME, including magnetic resonance imaging (MRI), positron emission tomography and optical imaging. In this review, we will focus on the latest fluorescent-based methods for optical sensing of cell metabolism and MRI as diagnostic imaging tools applied both in vitro and in vivo. The primary emphasis will be on describing the current and future uses of systems that can measure intra- and extra-cellular pH and oxygen changes at high spatial and temporal resolution. In addition, the suitability of these approaches for mapping tumour heterogeneity, and assessing response or failure to therapeutics will also be covered.

Numerical study of COVID-19 spatial-temporal spreading in London

Abstract: Recent study reported that an aerosolised virus (COVID-19) can survive in the air for a few hours. It is highly possible that people get infected with the disease by breathing and contact with items contaminated by the aerosolised virus. However, the aerosolised virus transmission and trajectories in various meteorological environments remain unclear. This paper has investigated the movement of aerosolised viruses from a high concentration source across a dense urban area. The case study looks at the highly air polluted areas of London: University College Hospital (UCH) and King Cross and St Pancras International Station (KCSPI). We explored the spread and decay of COVID-19 released from the hospital and railway stations with the prescribed meteorological conditions. The study has three key findings: the primary result is that it is possible for the virus to travel from meters up to hundred meters from the source location. The secondary finding shows viruses released into the atmosphere from entry and exit points at KCSPI remain trapped within a small radial distance of < 50m. This strengthens the case for the use of face coverings to reduce the infection rate. The final finding shows that there are different levels of risk at various door locations for UCH, depending on which door is used there can be a higher concentration of COVID-19. Although our results are based on London, since the fundamental knowledge processes are the same, our study can be further extended to other locations (especially the highly air polluted areas) in the world.

DAN-Net: Dual-Domain Adaptive-Scaling Non-local Network for CT Metal Artifact Reduction

Abstract: Metal implants can heavily attenuate X-rays in computed tomography (CT) scans, leading to severe artifacts in reconstructed images, which significantly jeopardize image quality and negatively impact subsequent diagnoses and treatment planning. With the rapid development of deep learning in the field of medical imaging, several network models have been proposed for metal artifact reduction (MAR) in CT. Despite the encouraging results achieved by these methods, there is still much room to further improve performance. In this paper, a novel Dual-domain Adaptive-scaling Non-local network (DAN-Net) for MAR. We correct the corrupted sinogram using adaptive scaling first to preserve more tissue and bone details as a more informative input. Then, an end-to-end dual-domain network is adopted to successively process the sinogram and its corresponding reconstructed image generated by the analytical reconstruction layer. In addition, to better suppress the existing artifacts and restrain the potential secondary artifacts caused by inaccurate results of the sinogram-domain network, a novel residual sinogram learning strategy and nonlocal module are leveraged in the proposed network model. In the experiments, the proposed DAN-Net demonstrates performance competitive with several state-of-the-art MAR methods in both qualitative and quantitative aspects.

Fully automated 3D segmentation of dopamine transporter SPECT images using an estimation-based approach

Abstract: Quantitative measures of uptake in caudate, putamen, and globus pallidus in dopamine transporter (DaT) brain SPECT have potential as biomarkers for the severity of Parkinson disease. Reliable quantification of uptake requires accurate segmentation of these regions. However, segmentation is challenging in DaT SPECT due to partial-volume effects, system noise, physiological variability, and the small size of these regions. To address these challenges, we propose an estimation-based approach to segmentation. This approach estimates the posterior mean of the fractional volume occupied by caudate, putamen, and globus pallidus within each voxel of a 3D SPECT image. The estimate is obtained by minimizing a cost function based on the binary cross-entropy loss between the true and estimated fractional volumes over a population of SPECT images, where the distribution of the true fractional volumes is obtained from magnetic resonance images from clinical populations. The proposed method accounts for both the sources of partial-volume effects in SPECT, namely the limited system resolution and tissue-fraction effects. The method was implemented using an encoder-decoder network and evaluated using realistic clinically guided SPECT simulation studies, where the ground-truth fractional volumes were known. The method significantly outperformed all other considered segmentation methods and yielded accurate segmentation with dice similarity coefficients of ~ 0.80 for all regions. The method was relatively insensitive to changes in voxel size. Further, the method was relatively robust up to +/- 10 degrees of patient head tilt along transaxial, sagittal, and coronal planes. Overall, the results demonstrate the efficacy of the proposed method to yield accurate fully automated segmentation of caudate, putamen, and globus pallidus in 3D DaT-SPECT images.

THz detection of biomolecules in aqueous environments -- status and perspectives for analysis under physiological conditions and clinical use

Abstract: Bioanalytical THz sensing techniques have proven to be an interesting and viable tool for the label-free detection and analysis of biomolecules. However, a major challenge for THz bioanalytics is to perform investigations in the native aqueous environments of the analytes. This review recapitulates the status and future requirements for establishing THz biosensing as a complementary toolbox in the repertoire of standard bioanalytic methods. The potential use in medical research and clinical diagnosis is discussed. Under these considerations, this article presents a comprehensive categorization of biochemically relevant analytes that have been investigated by THz sensing techniques in aqueous media. The detectable concentration levels of ions, carbohydrates, (poly-)nucleotides, active agents, proteins and different biomacromolecules from THz experiments are compared to characteristic physiological concentrations and lower detection limits of state-of-the-art bioanalytical methods. Finally, recent experimental developments and achievements are discussed, which potentially pave the way for THz analysis of biomolecules under clinically relevant conditions.

A physiologically realistic virtual patient database for the study of arterial haemodynamics

Abstract: This study creates a physiologically realistic virtual patient database (VPD), representing the human arterial system, for the primary purpose of studying the affects of arterial disease on haemodynamics. A low dimensional representation of an anatomically detailed arterial network is outlined, and a physiologically realistic posterior distribution for its parameters is constructed through a Bayesian approach. This approach combines both physiological/geometrical constraints and the available measurements reported in the literature. A key contribution of this work is to present a framework for including all such available information for the creation of virtual patients (VPs). The Markov Chain Monte Carlo (MCMC) method is used to sample random VPs from this posterior distribution, and the pressure and flow-rate profiles associated with the VPs are computed through a model of pulse wave propagation. This combination of the arterial network parameters (representing the VPs) and the haemodynamics waveforms of pressure and flow-rates at various locations (representing functional response of the VPs) makes up the VPD. While 75,000 VPs are sampled from the posterior distribution, 10,000 are discarded as the initial burn-in period. A further 12,857 VPs are subsequently removed due to the presence of negative average flow-rate. Due to an undesirable behaviour observed in some VPs -- asymmetric under- and over-damped pressure and flow-rate profiles in the left and right sides of the arterial system -- a filter is proposed for their removal. The final VPD has 28,868 subjects. It is shown that the methodology is appropriate by comparing the VPD statistics to those reported in literature across real populations. A good agreement between the two is found while respecting physiological/geometrical constraints. The pre-filter database is made available at this https URL.

Auditory-visual scenes for hearing research

Abstract: While experimentation with synthetic stimuli in abstracted listening situations has a long standing and successful history in hearing research, an increased interest exists on closing the remaining gap towards real-life listening by replicating situations with high ecological validity in the lab. This is important for understanding the underlying auditory mechanisms and their relevance in real-life situations as well as for developing and evaluating increasingly sophisticated algorithms for hearing assistance. A range of 'classical' stimuli and paradigms have evolved to de-facto standards in psychoacoustics, which are simplistic and can be easily reproduced across laboratories. While they ideally allow for across laboratory comparisons and reproducible research, they, however, lack the acoustic stimulus complexity and the availability of visual information as observed in everyday life communication and listening situations. This contribution aims to provide and establish an extendable set of complex auditory-visual scenes for hearing research that allow for ecologically valid testing in realistic scenes while also supporting reproducibility and comparability of scientific results. Three virtual environments are provided (underground station, pub, living room), consisting of a detailed visual model, an acoustic geometry model with acoustic surface properties as well as a set of acoustic measurements in the respective real-world environments. The current data set enables i) audio-visual research in a reproducible set of environments, ii) comparison of room acoustic simulation methods with "ground truth" acoustic measurements, iii) a condensation point for future extensions and contributions for developments towards standardized test cases for ecologically valid hearing research in complex scenes.

Failure Tolerance of the Human Lumbar Spine in Combined Compression and Flexion Loading

Abstract: Vehicle safety systems have substantially decreased motor vehicle crash-related injuries and fatalities, but injuries to the lumbar spine still have been reported Experimental and computational analyses of upright and, particularly, reclined occupants in frontal crashes have shown that the lumbar spine can be subjected to axial compression followed by combined compression-flexion loading Lumbar spine failure tolerance in combined compression-flexion has not been widely explored in the literature Therefore, the goal of this study was to measure the failure tolerance of the lumbar spine in combined compression and flexion Forty 3-vertebra lumbar spine segments were pre-loaded with axial compression and then subjected to dynamic flexion bending until failure Clinically relevant middle vertebra fractures were observed in twenty-one of the specimens, including compression and burst fractures The remaining nineteen specimens experienced failure at the potting grip interface Since specimen characteristics and pre-test axial load varied widely within the sample, failure forces (mean 34 kN, range 16-51 kN) and moments (mean 73 Nm, range 0-181 Nm) also varied widely Tobit univariate regressions were performed to determine the relationship between censored failure tolerance and specimen sex, segment type (upper/lower), age, and cross-sectional area Age, sex, and cross-sectional area significantly affected failure force and moment individually (p<00024) These data can be used to develop injury prediction tools for lumbar spine fractures and further research in future safety systems

Choroidal vasculature imaging with laser Doppler holography

Abstract: The choroid is a highly vascularized tissues supplying the retinal pigment epithelium and photoreceptors. Its implication in retinal diseases is gaining increasing interest. However, investigating the anatomy and flow of the choroid remains challenging. Here we show that laser Doppler holography provides high contrast imaging of choroidal vessels in humans, with a spatial resolution comparable to state of the art indocyanine green angiography and optical coherence tomography. Additionally, laser Doppler holography contributes to sort out choroidal arteries and veins by using a power Doppler spectral analysis. We thus demonstrate the potential of laser Doppler holography to improve our understanding of the anatomy and flow of the choroidal vascular network.

Deep learning-based attenuation correction in the image domain for myocardial perfusion SPECT imaging

Abstract: Objective: In this work, we set out to investigate the accuracy of direct attenuation correction (AC) in the image domain for the myocardial perfusion SPECT imaging (MPI-SPECT) using two residual (ResNet) and UNet deep convolutional neural networks. Methods: The MPI-SPECT 99mTc-sestamibi images of 99 participants were retrospectively examined. UNet and ResNet networks were trained using SPECT non-attenuation corrected images as input and CT-based attenuation corrected SPECT images (CT-AC) as reference. The Chang AC approach, considering a uniform attenuation coefficient within the body contour, was also implemented. Quantitative and clinical evaluation of the proposed methods were performed considering SPECT CT-AC images of 19 subjects as reference using the mean absolute error (MAE), structural similarity index (SSIM) metrics, as well as relevant clinical indices such as perfusion deficit (TPD). Results: Overall, the deep learning solution exhibited good agreement with the CT-based AC, noticeably outperforming the Chang method. The ResNet and UNet models resulted in the ME (count) of ${-6.99\pm16.72}$ and ${-4.41\pm11.8}$ and SSIM of ${0.99\pm0.04}$ and ${0.98\pm0.05}$, respectively. While the Change approach led to ME and SSIM of ${25.52\pm33.98}$ and ${0.93\pm0.09}$, respectively. Similarly, the clinical evaluation revealed a mean TPD of ${12.78\pm9.22}$ and ${12.57\pm8.93}$ for the ResNet and UNet models, respectively, compared to ${12.84\pm8.63}$ obtained from the reference SPECT CT-AC images. On the other hand, the Chang approach led to a mean TPD of ${16.68\pm11.24}$. Conclusion: We evaluated two deep convolutional neural networks to estimate SPECT-AC images directly from the non-attenuation corrected images. The deep learning solutions exhibited the promising potential to generate reliable attenuation corrected SPECT images without the use of transmission scanning.

Quantitative analysis of image quality in low-dose CT imaging for Covid-19 patients

Abstract: We set out to simulate four reduced dose-levels (60%-dose, 40%-dose, 20%-dose, and 10%-dose) of standard CT imaging using Beer-Lambert's law across 49 patients infected with COVID-19. Then, three denoising filters, namely Gaussian, Bilateral, and Median, were applied to the different low-dose CT images, the quality of which was assessed prior to and after the application of the various filters via calculation of peak signal-to-noise ratio (PSNR), root mean square error (RMSE), structural similarity index measure (SSIM), and relative CT-value bias, separately for the lung tissue and whole-body. The quantitative evaluation indicated that 10%-dose CT images have inferior quality (with RMSE=322.1-+104.0 HU and bias=11.44-+4.49% in the lung) even after the application of the denoising filters. The bilateral filter exhibited superior performance to suppress the noise and recover the underlying signals in low-dose CT images compared to the other denoising techniques. The bilateral filter led to RMSE and bias of 100.21-+16.47 HU, -0.21-+1.20%, respectively in the lung regions for 20%-dose CT images compared to the Gaussian filter with RMSE=103.46-+15.70 HU and bias=1.02-+1.68%, median filter with RMSE=129.60-+18.09 HU and bias=-6.15-+2.24%, and the nonfiltered 20%-dose CT with RMSE=217.37-+64.66 HU and bias=4.30-+1.85%. In conclusion, the 20%-dose CT imaging followed by the bilateral filtering introduced a reasonable compromise between image quality and patient dose reduction.

Bragg-Curve Simulation of Carbon-Ion Beams for Particle-therapy Applications: a study with the GEANT4 toolkit

Abstract: We used the GEANT4 Monte Carlo MC Toolkit to simulate carbon ion beams incident on water, tissue, and bone, taking into account nuclear fragmentation reactions. Upon increasing the energy of the primary beam, the position of the Bragg-Peak transfers to a location deeper inside the phantom. For different materials, the peak is located at a shallower depth along the beam direction and becomes sharper with increasing electron density NZ. Subsequently, the generated depth dose of the Bragg curve is then benchmarked with experimental data from GSI in Germany. The results exhibit a reasonable correlation with GSI experimental data with an accuracy of between 0.02 and 0.08 cm, thus establishing the basis to adopt MC in heavy-ion treatment planning. The Kolmogorov-Smirnov K-S test further ascertained from a statistical point of view that the simulation data matched the experimentally measured data very well. The two-dimensional isodose contours at the entrance were compared to those around the peak position and in the tail region beyond the peak, showing that bone produces more dose, in comparison to both water and tissue, due to secondary doses. In the water, the results show that the maximum energy deposited per fragment is mainly attributed to secondary carbon ions, followed by secondary boron and beryllium. Furthermore, the number of protons produced is the highest, thus making the maximum contribution to the total dose deposition in the tail region. Finally, the associated spectra of neutrons and photons were analyzed. The mean neutron energy value was found to be 16.29 MeV, and 1.03 MeV for the secondary gamma. However, the neutron dose was found to be negligible as compared to the total dose due to their longer range.

Continental generalization of an AI system for clinical seizure recognition

Abstract: Electroencephalogram (EEG) monitoring and objective seizure identification is an essential clinical investigation for some patients with epilepsy. Accurate annotation is done through a time-consuming process by EEG specialists. Computer-assisted systems for seizure detection currently lack extensive clinical utility due to retrospective, patient-specific, and/or irreproducible studies that result in low sensitivity or high false positives in clinical tests. We aim to significantly reduce the time and resources on data annotation by demonstrating a continental generalization of seizure detection that balances sensitivity and specificity. This is a prospective inference test of artificial intelligence on nearly 14,590 hours of adult EEG data from patients with epilepsy between 2011 and 2019 in a hospital in Sydney, Australia. The inference set includes patients with different types and frequencies of seizures across a wide range of ages and EEG recording hours. We validated our inference model in an AI-assisted mode with a human expert arbiter and a result review panel of expert neurologists and EEG specialists on 66 sessions to demonstrate achievement of the same performance with over an order-of-magnitude reduction in time. Our inference on 1,006 EEG recording sessions on the Australian dataset achieved 76.68% with nearly 56 [0, 115] false alarms per 24 hours on average, against legacy ground-truth annotations by human experts, conducted independently over nine years. Our pilot test of 66 sessions with a human arbiter, and reviewed ground truth by a panel of experts, confirmed an identical human performance of 92.19% with an AI-assisted system, while the time requirements reduce significantly from 90 to 7.62 minutes on average.

Real-time non-rigid 3D respiratory motion estimation for MR-guided radiotherapy using MR-MOTUS

Abstract: The MR-Linac is a combination of an MR-scanner and radiotherapy linear accelerator (Linac) which holds the promise to increase the precision of radiotherapy treatments with MR-guided radiotherapy by monitoring motion during radiotherapy with MRI, and adjusting the radiotherapy plan accordingly. Optimal MR-guidance for respiratory motion during radiotherapy requires MR-based 3D motion estimation with a latency of 200-500 ms. Currently this is still challenging since typical methods rely on MR-images, and are therefore limited to the 3D MR-imaging latency. In this work, we present a method to perform non-rigid 3D respiratory motion estimation with 170 ms latency, including both acquisition and reconstruction. The proposed method called real-time low-rank MR-MOTUS reconstructs motion-fields directly from k-space data, and leverages an explicit low-rank decomposition of motion-fields to split the large scale 3D+t motion-field reconstruction problem posed in our previous work into two parts: (I) a medium-scale offline preparation phase and (II) a small-scale online inference phase which exploits the results of the offline phase for real-time computations. The method was validated on free-breathing data of five volunteers, acquired with a 1.5T Elekta Unity MR-Linac. Results show that the reconstructed 3D motion-fields are anatomically plausible, highly correlated with a self-navigation motion surrogate (R = 0.975 +/- 0.0110), and can be reconstructed with a total latency of 170 ms that is sufficient for real-time MR-guided abdominal radiotherapy.

Radioactive nuclei for $\beta^+\gamma$ PET and theranostics: selected candidates

Abstract: PET is an established medical diagnostic imaging method. Continuous improvements are aimed at refining image reconstruction, reducing the amount of radioactive tracer and combining with targeted therapy. TOF-PET provides the localization of the tracer through improved time resolution, nuclear physics may contribute to this goal via selection of radioactive nuclei emitting additional $\gamma$-rays. This additional radiation, when properly detected, localizes the decay of the tracer at the line of response determined by two detected 511 keV quanta. Selected candidates are presented. Some are particularly interesting, as they are strong candidates for theranostic applications.

A Stroke Detection and Discrimination Framework using Broadband Microwave Scattering on Stochastic Models With Deep Learning

Abstract: Stroke poses an immense public health burden and remains among the primary causes of death and disability worldwide. Emergent therapy is often precluded by late or indeterminate times of onset before initial clinical presentation. Rapid, mobile, safe and low-cost stroke detection technology remains a deeply unmet clinical need. Past studies have explored the use of microwave and other small form-factor strategies for rapid stroke detection; however, widespread clinical adoption remains unrealized. Here, we investigated the use of microwave scattering perturbations from ultra-wide-band antenna arrays to learn dielectric signatures of disease. Two deep neural networks (DNNs) were used for: 1) stroke detection ("classification network"), and 2) characterization of the hemorrhage location and size ("discrimination network"). Dielectric signatures were learned on a simulated cohort of 666 hemorrhagic stroke and control subjects using 2D stochastic head models. The classification network yielded a stratified K-fold stroke detection accuracy of 94.60% with standard deviation of 2.86% and an AUC of 0.996, while the discrimination network resulted in a mean squared error of <0.004 cm and <0.02 cm, for the stroke localization and size estimation, respectively. We report a novel approach to intelligent diagnostics using microwave wide-band scattering information thus circumventing conventional image-formation.

Prepolarized MRI of Hard Tissues and Solid-State Matter

Abstract: Prepolarized Magnetic Resonance Imaging (PMRI) is a long-established technique conceived to counteract the loss in signal-to-noise ratio (SNR) inherent to low-field MRI systems. When it comes to hard biological tissues and solid-state matter, PMRI is severely restricted by their ultra-short characteristic relaxation times. Here we demonstrate that efficient hard tissue prepolarization is within reach with a special-purpose 0.26 T scanner designed for dental MRI and equipped with suitable high-power electronics. We have characterized the performance of a 0.5 T prepolarizer module which can be switched on and off in just 200 us. To that end, we have used resin, dental and bone samples, all with T1 times in the order of 20 ms at our field strength. The measured SNR enhancement is in good agreement with a simple theoretical model, and small deviations in extreme regimes can be attributed to mechanical vibrations due to the magnetic interaction between the prepolarization and main magnets. Finally, we argue that these results can be applied to clinical dental imaging, opening the door to replacing hazardous X-ray systems with low-field PMRI scanners.

Probabilistic feature extraction, dose statistic prediction and dose mimicking for automated radiation therapy treatment planning

Abstract: Purpose: We propose a general framework for quantifying predictive uncertainties of dose-related quantities and leveraging this information in a dose mimicking problem in the context of automated radiation therapy treatment planning. Methods: A three-step pipeline, comprising feature extraction, dose statistic prediction and dose mimicking, is employed. In particular, the features are produced by a convolutional variational autoencoder and used as inputs in a previously developed nonparametric Bayesian statistical method, estimating the multivariate predictive distribution of a collection of predefined dose statistics. Specially developed objective functions are then used to construct a probabilistic dose mimicking problem based on the produced distributions, creating deliverable treatment plans. Results: The numerical experiments are performed using a dataset of 94 retrospective treatment plans of prostate cancer patients. We show that the features extracted by the variational autoencoder capture geometric information of substantial relevance to the dose statistic prediction problem and are related to dose statistics in a more regularized fashion than hand-crafted features. The estimated predictive distributions are reasonable and outperforms a non-input-dependent benchmark method, and the deliverable plans produced by the probabilistic dose mimicking agree better with their clinical counterparts than for a non-probabilistic formulation. Conclusions: We demonstrate that prediction of dose-related quantities may be extended to include uncertainty estimation and that such probabilistic information may be leveraged in a dose mimicking problem. The treatment plans produced by the proposed pipeline resemble their original counterparts well, illustrating the merits of a holistic approach to automated planning based on probabilistic modeling.

Skin marker-based subject-specific spinal alignment modeling: A feasibility study

Abstract: Musculoskeletal models have the potential to improve diagnosis and optimize clinical treatment by predicting accurate outcomes on an individual basis. However, the subject-specific modeling of spinal alignment is often strongly simplified or is based on radiographic assessments, exposing subjects to unnecessary radiation. We therefore developed a novel skin marker-based approach for modeling subject-specific spinal alignment and evaluated its feasibility by comparing the predicted with the actual intervertebral joint (IVJ) locations/orientations (ground truth) using lateral-view radiographic images. Moreover, the predictive performance of the subject-specific models was evaluated by comparing the predicted L1/L2 spinal loads during various functional activities with in vivo measured data obtained from the OrthoLoad database. IVJ locations/orientations were predicted closer to ground truth as opposed to standard model scaling, with average location prediction errors of 0.99+/-0.68 cm on the frontal and 1.21+/-0.97 cm on the transverse axis as well as an average orientation prediction error of 4.74°+/-2.80°. Simulated spinal loads showed similar curve patterns but considerably larger values as compared to in vivo measured data. Differences in spinal loads between generic and subject-specific models become only apparent on an individual subject level. These results underline the feasibility of the proposed method and associated workflow for inter- and intra-subject investigations using musculoskeletal simulations. When implemented into standard model scaling workflows, it is expected to improve the accuracy of muscle activity and joint loading simulations, which is crucial for investigations of treatment effects or pathology-dependent deviations.

Performance of an Enhanced Afterloader with Electromagnetic Tracking Capabilities for Channel Reconstruction and Error Detection

Abstract: The reconstruction and error detection performance of a Flexitron research unit equipped with a special check-cable in which is integrated an EM sensor (NDI Aurora V3) that enables tracking and reconstruction capability was assessed. Reconstructions of a 24 cm long catheter were performed using two methods: continuous fixed speed check cable backward stepping and stepping through each dwell position every 1mm. The ability of the system to differentiate between two closely located (parallel) catheters was investigated by connecting catheters to the afterloader and moving it from its axis with an increment of 1mm. A robotic arm was used to move the catheter between each reconstruction. Reconstructions were obtained with a locally weighted scatter-plot smoothing algorithm. To quantify the reconstruction accuracy, distances between two catheters were computed along the reconstruction track with a 5mm step. Reconstructions of curve catheter paths were assessed through parallel and perpendicular phantom configuration to the EM field generator. Indexer length and lateral error were simulated and ROC analysis was made. Using a 50cm/s check cable speed does not allow accurate reconstruction. A slower check cable speed results in better reconstruction performance and smaller standard deviations. The optimum operating distance from the field generator (50-300mm) resulted in a lower absolute mean deviation from the expected value (0.2pm0.1mm) versus being positioned on the edge of the electromagnetic sensitive detection volume (0.6pm0.3mm). Reconstructions of curved and circular paths with a check cable speed under 5cm/s gave a 0.8pm0.5mm error, or better. All indexer and lateral shifts of 1mm were detected with a check cable speed of 2.5cm/s or under. A speed under 5cm/s is recommended for straight and curved catheter reconstructions. The check cable can also be used to detect common shift errors.