• ### Diffusion time dependence of microstructural parameters in fixed spinal cord(1705.03630)

March 19, 2019 physics.med-ph, physics.bio-ph
Biophysical modelling of diffusion MRI is necessary to provide specific microstructural tissue properties. However, estimating model parameters from data with limited diffusion gradient strength, such as clinical scanners, has proven unreliable due to a shallow optimization landscape. On the other hand, estimation of diffusion kurtosis (DKI) parameters is more robust as the clinical acquisitions typically probe a regime in which the associated 4th order cumulant expansion is adequate; however, its parameters are not microstructurally specific a priori. Given an appropriate biophysical model, its parameters may be connected to DKI parameters, but it was previously shown that at the DKI level, it still does not provide sufficient information to uniquely determine all model parameters. Earlier work has shown that by neglecting axonal dispersion, this parameter degeneracy reduces to the question of whether intra-axonal diffusivity is larger than or smaller than extra-axonal diffusivity. Here we develop a model of diffusion in spinal cord white matter including axonal dispersion and demonstrate stable estimation of all model parameters from DKI. By employing the recently developed fast axisymmetric DKI, we use stimulated echo acquisition mode to collect data over an unprecedented diffusion time range with very narrow diffusion gradient pulses, enabling finely resolved measurements of diffusion time dependence of both net diffusion and kurtosis metrics, as well as model intra- and extra-axonal diffusivities, and axonal dispersion. Our results demonstrate substantial time dependence of all parameters except volume fractions, and the additional time dimension provides support for intra-axonal diffusivity to be larger than extra-axonal diffusivity in spinal cord white matter, although not unambiguously. We compare our findings to predictions from effective medium theory.
• ### A volumetric Modulated Arc Therapy (VMAT) combined with deep inspiration breath hold (DIBH) technique for adjuvant irradiation for left sided breast cancer(1706.06867)

March 13, 2019 physics.med-ph
Background: The deep inspiration breath hold technique (DIBH) is widely used for left-sided breast radiotherapy (RT), in order to reduce the dose to the heart and the risk of cardiotoxicity. The volumetric arc therapy (VMAT) technique increases the dose conformity to the planning target volume (PTV). We hereby describe the procedure of combining VMAT and DIBH without a dedicated equipment, and report a dosimetric study related on our implemented technique. Methods: The DIBH technique is based on voluntary breath hold of the patient which is controlled by a laser on tattoo marks. Patients were selected depending their ability to breath hold for 25s, and VMAT plans were optimized with an arrangement of 4 arcs with a maximum delivery time of 20s each. A retrospective dosimetric study was undertaken on 30 patients treated with this technique: 10 received local (whole breast) RT (L-RT) with a simultaneous integrated boost (SIB), 10 a locoregional (LR)RT of the chestwall, 10 a LR RT of the whole breast with SIB. For these patients, their DIBH-VMAT treatment plans were compared to free breathing (FB) VMAT plans, and dose characteristics of PTV and organs at risk (OAR) were evaluated. Results: The PTV coverage was identical among plans under DIBH or FB. The mean heart dose was reduced from 2.8+/-0.9Gy in FB to 1.8+/-0.6Gy in DIBH for local RT, and reduced from 3.4+/-1.3Gy in FB to 2.2+/-0.5Gy in DIBH for LR-RT. The ipsilateral mean lung dose was reduced by a factor of 12% for LR-RT, but remained unchanged for L-RT with DIBH-VMAT compared to FB treatment. No significant differences were observed in the mean dose of contra-lateral organs with either DIBH or FB VMAT techniques. Conclusion: The VMAT-DIBH technique performed without a dedicated breath hold equipment, provided important dose reductions at the heart, while preserving the ipsilateral lung and contralateral normal tissues.
• ### Wideband Excitation of Microbubbles with Chirps to Maximize the Sonoporation Efficiency and Contrast in Ultrasound Imaging(1606.02626)

Feb. 4, 2019 physics.med-ph
The importance of the excitation bandwidth is well known in diagnostic ultrasound imaging. However, the effect of excitation bandwidth in therapeutic applications of microbubbles has been mostly overlooked. A majority of contrast agent production techniques generate polydisperse microbubble populations, so a wide range of resonance frequencies exist. Therefore, wideband excitation is necessary to fully utilize microbubble resonance behavior and maximize the reradiated energy from a microbubble population, both for imaging and therapy. Oscillations of sixty SonoVue microbubbles in proximity of a rigid boundary were captured on a high speed camera at 3 Mfps, excited with a peak negative pressure of 50 kPa at 1 MHz. Measurements were analyzed according to their peak radiated pressure, radial oscillations, root mean squared pressure, and shear stress generated by microbubbles. Results showed that long duration and wideband excitation at low intensity levels was preferable for sonoporation, where microbubbles can be driven in a stable oscillation state without experiencing inertial cavitation or destruction.
• ### Maximum entropy based non-negative optoacoustic tomographic image reconstruction(1707.08391)

Objective:Optoacoustic (photoacoustic) tomography is aimed at reconstructing maps of the initial pressure rise induced by the absorption of light pulses in tissue. In practice, due to inaccurate assumptions in the forward model, noise and other experimental factors, the images are often afflicted by artifacts, occasionally manifested as negative values. The aim of the work is to develop an inversion method which reduces the occurrence of negative values and improves the quantitative performance of optoacoustic imaging. Methods: We present a novel method for optoacoustic tomography based on an entropy maximization algorithm, which uses logarithmic regularization for attaining non-negative reconstructions. The reconstruction image quality is further improved using structural prior based fluence correction. Results: We report the performance achieved by the entropy maximization scheme on numerical simulation, experimental phantoms and in-vivo samples. Conclusion: The proposed algorithm demonstrates superior reconstruction performance by delivering non-negative pixel values with no visible distortion of anatomical structures. Significance: Our method can enable quantitative optoacoustic imaging, and has the potential to improve pre-clinical and translational imaging applications.
• ### Correction of Blurring due to a Difference in Scanning Direction of Field-Free Line in Projection-Based Magnetic Particle Imaging(1606.05740)

Jan. 7, 2019 physics.med-ph
In projection-based magnetic particle imaging (MPI) with a field-free-line (FFL) encoding scheme, projection data are usually acquired by moving the FFL in a zigzag and a difference in the projection data occurs depending on the scanning direction of FFL, resulting in blurring in the reconstructed images. In this study, we developed a method for correcting the blur by deconvolution using a signal-delay constant ({\xi}). The {\xi} value for correction ({\xi}c) was determined by acquiring projection data in positive and negative directions and searching for the {\xi} value which minimized the 2-norm between the deconvolved projection data in the two directions. We validated our method using a line and A-shaped phantoms for various velocities of FFL (vFFL). The {\xi}c value correlated linearly with vFFL. The full width at half maximum of the line phantom decreased significantly after correction of the blur. The effectiveness of our method was also confirmed by the MPI images of the A-shaped phantom. These results suggest that our method will be useful for enhancing the reliability of projection-based MPI.
• ### Development of a method for measuring blood coagulation using superparamagnetic iron oxide nanoparticles and an alternating magnetic field(1605.00255)

Jan. 5, 2019 physics.med-ph
We developed a method for measuring blood coagulation using superparamagnetic iron oxide nanoparticles (SPIONs) and an alternating magnetic field (AMF). The 3rd and 5th harmonic signals from SPIONs mixed with blood induced by AMF were detected using a gradiometer coil. Blood coagulation was induced artificially by adding CaCl2 solution to whole blood of sheep at various temperatures and hematocrits. We calculated the coagulation rate (k) and normalized signal intensity at infinite time (Sinf) by fitting the time course of the normalized 3rd harmonic signal to S(t)=(1-Sinf)exp(-kt)+Sinf. The k values increased significantly with increasing temperature and decreased significantly with increasing hematocrit. The Sinf values decreased significantly with increasing temperature and tended to increase with increasing hematocrit. Blood anticoagulation was induced by adding heparin to the whole blood sampled from mice. There were significant differences in both the 3rd and 5th harmonic signals between groups with and without heparin at 25 min or more after adding heparin. We also calculated the 3rd and 5th harmonic signals for viscosities ranging from 0.001 to 1 kg/m/s, with an assumption that the magnetization and particle size distribution of SPIONs obey the Langevin theory of paramagnetism and log-normal distribution, respectively. The 3rd and 5th harmonic signals increased slowly with increasing viscosity and had peaks at approximately 0.015 and 0.025 kg/m/s, respectively. After these peaks, they decreased monotonically with increasing viscosity. These results confirm the rationale of our method. In conclusion, our method will be useful for measuring blood coagulation and anticoagulation and for studying their processes.
• ### Optimized Quantification of Spin Relaxation Times in the Hybrid State(1703.00481)

Dec. 3, 2018 physics.med-ph
Purpose: The analysis of optimized spin ensemble trajectories for relaxometry in the hybrid state. Methods: First, we constructed visual representations to elucidate the differential equation that governs spin dynamics in hybrid state. Subsequently, numerical optimizations were performed to find spin ensemble trajectories that minimize the Cram\'er-Rao bound for $T_1$-encoding, $T_2$-encoding, and their weighted sum, respectively, followed by a comparison of the Cram\'er-Rao bounds obtained with our optimized spin-trajectories, as well as Look-Locker and multi-spin-echo methods. Finally, we experimentally tested our optimized spin trajectories with in vivo scans of the human brain. Results: After a nonrecurring inversion segment on the southern hemisphere of the Bloch sphere, all optimized spin trajectories pursue repetitive loops on the northern half of the sphere in which the beginning of the first and the end of the last loop deviate from the others. The numerical results obtained in this work align well with intuitive insights gleaned directly from the governing equation. Our results suggest that hybrid-state sequences outperform traditional methods. Moreover, hybrid-state sequences that balance $T_1$- and $T_2$-encoding still result in near optimal signal-to-noise efficiency. Thus, the second parameter can be encoded at virtually no extra cost. Conclusion: We provide insights regarding the optimal encoding processes of spin relaxation times in order to guide the design of robust and efficient pulse sequences. We find that joint acquisitions of $T_1$ and $T_2$ in the hybrid state are substantially more efficient than sequential encoding techniques.
• ### A feasibility study for predicting optimal radiation therapy dose distributions of prostate cancer patients from patient anatomy using deep learning(1709.09233)

Nov. 29, 2018 cs.AI, cs.CV, physics.med-ph, cs.LG
With the advancement of treatment modalities in radiation therapy for cancer patients, outcomes have improved, but at the cost of increased treatment plan complexity and planning time. The accurate prediction of dose distributions would alleviate this issue by guiding clinical plan optimization to save time and maintain high quality plans. We have modified a convolutional deep network model, U-net (originally designed for segmentation purposes), for predicting dose from patient image contours of the planning target volume (PTV) and organs at risk (OAR). We show that, as an example, we are able to accurately predict the dose of intensity-modulated radiation therapy (IMRT) for prostate cancer patients, where the average Dice similarity coefficient is 0.91 when comparing the predicted vs. true isodose volumes between 0% and 100% of the prescription dose. The average value of the absolute differences in [max, mean] dose is found to be under 5% of the prescription dose, specifically for each structure is [1.80%, 1.03%](PTV), [1.94%, 4.22%](Bladder), [1.80%, 0.48%](Body), [3.87%, 1.79%](L Femoral Head), [5.07%, 2.55%](R Femoral Head), and [1.26%, 1.62%](Rectum) of the prescription dose. We thus managed to map a desired radiation dose distribution from a patient's PTV and OAR contours. As an additional advantage, relatively little data was used in the techniques and models described in this paper.
• ### An analytical Model which Determines the Apparent T1 for Modified Look-Locker Inversion Recovery (MOLLI) -- Analysis of the Longitudinal Relaxation under the Influence of Discontinuous Balanced and Spoiled Gradient Echo Readouts(1704.06898)

Quantitative nuclear magnetic resonance imaging (MRI) shifts more and more into the focus of clinical research. Especially determination of relaxation times without/and with contrast agents becomes the foundation of tissue characterization, e.g. in cardiac MRI for myocardial fibrosis. Techniques which assess longitudinal relaxation times rely on repetitive application of readout modules, which are interrupted by free relaxation periods, e.g. the Modified Look-Locker Inversion Recovery = MOLLI sequence. These discontinuous sequences reveal an apparent relaxation time, and, by techniques extrapolated from continuous readout sequences, the real T1 is determined. What is missing is a rigorous analysis of the dependence of the apparent relaxation time on its real partner, readout sequence parameters and biological parameters as heart rate. This is provided in this paper for the discontinuous balanced steady state free precession (bSSFP) and spoiled gradient echo readouts. It turns out that the apparente longitudinal relaxation rate is the time average of the relaxation rates during the readout module, and free relaxation period. Knowing the heart rate our results vice versa allow to determine the real T1 from its measured apparent partner.
• ### Full modelling of high-intensity focused ultrasound and thermal heating in the kidney using realistic patient models(1708.02860)

Oct. 8, 2018 physics.med-ph
Objective: High-intensity focused ultrasound (HIFU) therapy can be used for non-invasive treatment of kidney (renal) cancer, but the clinical outcomes have been variable. In this study, the efficacy of renal HIFU therapy was studied using nonlinear acoustic and thermal simulations in three patients. Methods: The acoustic simulations were conducted with and without refraction in order to investigate its effect on the shape, size and pressure distribution at the focus. The values for the attenuation, sound speed, perfusion and thermal conductivity of the kidney were varied over the reported ranges to determine the effect of variability on heating. Furthermore, the phase aberration was studied in order to quantify the underlying phase shifts using a second order polynomial function. Results: The ultrasound field intensity was found to drop on average 11.1 dB with refraction and 6.4 dB without refraction. Reflection at tissue interfaces was found to result in a loss less than 0.1 dB. Focal point splitting due to refraction significantly reduced the heating efficacy. Perfusion did not have a large effect on heating during short sonication durations. Small changes in temperature were seen with varying attenuation and thermal conductivity, but no visible changes were present with sound speed variations. The aberration study revealed an underlying trend in the spatial distribution of the phase shifts. Conclusion: The results show that the efficacy of HIFU therapy in the kidney could be improved with aberration correction. Significance: A method is proposed by which patient specific pre-treatment calculations could be used to overcome the aberration and therefore make ultrasound treatment possible.
• ### ENLIVE: An Efficient Nonlinear Method for Calibrationless and Robust Parallel Imaging(1706.09780)

Sept. 17, 2018 physics.med-ph
Robustness against data inconsistencies, imaging artifacts and acquisition speed are crucial factors limiting the possible range of applications for magnetic resonance imaging (MRI). Therefore, we report a novel calibrationless parallel imaging technique which simultaneously estimates coil profiles and image content in a relaxed forward model. Our method is robust against a wide class of data inconsistencies, minimizes imaging artifacts and is comparably fast combining important advantages of many conceptually different state-of-the-art parallel imaging approaches. Depending on the experimental setting, data can be undersampled well below the Nyquist limit. Here, even high acceleration factors yield excellent imaging results while being robust to noise and the occurrence of phase singularities in the image domain, as we show on different data. Moreover, our method successfully reconstructs acquisitions with insufficient field-of-view. We further compare our approach to ESPIRiT and SAKE using spin-echo and gradient echo MRI data from the human head and knee. In addition, we show its applicability to non-Cartesian imaging on radial FLASH cardiac MRI data. Using theoretical considerations, we show that ENLIVE can be related to a low-rank formulation of blind multi-channel deconvolution, explaining why it inherently promotes low-rank solutions.
• ### A Blind Deconvolution Technique Based on Projection Onto Convex Sets for Magnetic Particle Imaging(1705.07506)

Jan. 29, 2020 physics.med-ph
Magnetic Particle Imaging (MPI) is an emerging imaging modality that maps the spatial distribution of magnetic nanoparticles. The x-space reconstruction in MPI results in highly blurry images, where the resolution depends on both system parameters and nanoparticle type. Previous techniques to counteract this blurring rely on the knowledge of the imaging point spread function (PSF), which may not be available or may require additional measurements. This work proposes a blind deconvolution algorithm for MPI to recover the precise spatial distribution of nanoparticles. The proposed algorithm exploits the observation that the imaging PSF in MPI has zero phase in Fourier domain. Thus, even though the reconstructed images are highly blurred, phase remains unaltered. We leverage this powerful property to iteratively enforce consistency of phase and bounded l1 energy information, using an orthogonal Projections Onto Convex Sets (POCS) algorithm. To demonstrate the method, comprehensive simulations were performed without and with nanoparticle relaxation effects, and at various noise levels. In addition, imaging experiments were performed on an in-house MPI scanner using a three-vial phantom that contained different nanoparticle types. Image quality was compared with conventional deconvolution methods, Wiener deconvolution and Lucy-Richardson method, which explicitly rely on the knowledge of PSF. Both the simulation results and experimental imaging results show that the proposed blind deconvolution algorithm outperforms the conventional deconvolution methods. Without utilizing the imaging PSF, the proposed algorithm improves image quality and resolution even in the case of different nanoparticle types, while displaying reliable performance against loss of the fundamental harmonic, nanoparticle relaxation effects, and noise.
• ### Invertibility of spectral x-ray data with pileup--two dimension-two spectrum case(1703.07557)

In the Alvarez-Macovski method, the line integrals of the x-ray basis set coefficients are computed from measurements with multiple spectra. An important question is whether the transformation from measurements to line integrals is invertible. This paper presents a proof that for a system with two spectra and a photon counting detector, pileup does not affect the invertibility of the system. If the system is invertible with no pileup, it will remain invertible with pileup although the reduced Jacobian may lead to increased noise.
• ### A mathematical model for DNA(1707.04283)

Recently, some authors have shown that a DNA molecule produces electromagnetic signals and communicates with other DNA molecules or other molecules. In fact, a DNA acts like a receiver or transmitter of radio waves. In this paper, we suggest a mathematical model for the DNA molecule and use of its communication to cure some diseases like cancer. In this model, first, by using concepts from string theory and M-theory, we calculate the energy of a DNA in terms of interactions between free electrons and bound electrons. We show that when a DNA is damaged, its energy changes and an extra current is produced. This extra current causes the electromagnetic signals of a damaged DNA molecule to be different when compared to the electromagnetic signals of a normal DNA molecule. The electromagnetic signals of a damaged DNA molecule induces an extra current in a normal DNA molecule and leads to its destruction. By sending crafted electromagnetic signals to normal DNA molecules and inducing an opposite current with respect to this extra current, we can prevent the destruction of normal DNA. Finally, we argue that the type of packing of DNA in chromosomes of men and women are different. This causes radiated waves from DNAs of men and women to have opposite signs and cancel the effect of each other in a pair. Using this property, we suggest another mechanism to cancel the effect of extra waves, which are produced by DNAs in cancer cells of a male or a female, by extra waves which are produced by DNAs in similar cells of a female or a male and prevent the progression of the disease.
• ### Quantifying birefringence in the bovine model of early osteoarthritis using polarisation-sensitive optical coherence tomography and mechanical indentation(1709.05129)

Aug. 1, 2018 physics.med-ph
Recent studies have shown potential for using polarisation sensitive optical coherence tomography (PS-OCT) to study cartilage morphology, and to be potentially used as an in-vivo, non-invasive tool for detecting osteoarthritic changes. However, there has been relatively limited ability of this method to quantify the subtle changes that occur in the early stages of cartilage degeneration. An established mechanical indenting technique that has previously been used to examine the microstructural response of articular cartilage was employed to fix the bovine samples in an indented state. The samples were subject to creep loading with a constant compressive stress of 4.5 MPa and, when imaged using PS-OCT, enabled birefringent banding patterns to be observed. The magnitude of the birefringence was quantified using the birefringence coefficient (BRC) and statistical analysis revealed that PS-OCT is able to detect and quantify significant changes between healthy and early osteoarthritic cartilage (p<0.001). This presents a novel utilization of PS-OCT for future development as an in-vivo assessment tool.
• ### High-Performance and Distributed Computing in a Probabilistic Finite Element Comparison Study of the Human Lower Leg Model with Total Knee Replacement(1611.07297)

July 18, 2018 physics.med-ph, cs.CE
Reliability theory is used to assess the sensitivity of a passive flexion and active flexion of the human lower leg Finite Element (FE) models with Total Knee Replacement (TKR) to the variability in the input parameters of the respective FE models. The sensitivity of the active flexion simulating the stair ascent of the human lower leg FE model with TKR was presented before in [1,2] whereas now in this paper a comparison is made with the passive flexion of the human lower leg FE model with TKR. First, with the Monte Carlo Simulation Technique (MCST), a number of randomly generated input data of the FE model(s) are obtained based on the normal standard deviations of the respective input parameters. Then a series of FE simulations are done and the output kinematics and peak contact pressures are obtained for the respective FE models (passive flexion and/or active flexion models). Seven output performance measures are reported for the passive flexion model and one more parameter was reported for the active flexion FE model (patello-femoral peak contact pressure) in [1]. A sensitivity study will be performed based on the Response Surface Method (RSM) to identify the key parameters that influence the kinematics and peak contact pressures of the passive flexion FE model. Another two MCST and RSM-based probabilistic FE analyses will be performed based on a reduced list of 19 key input parameters. In total 4 probabilistic FE analyses will be performed: 2 probabilistic FE analyses (MCST and RSM) based on an extended set of 78 input variables and another 2 probabilistic FE analyses (MCST and RSM) based on a reduced set of 19 input variables. Due to the likely computation cost in order to make hundreds of FE simulations with MCST, a high-performance and distributed computing system will be used for the passive flexion FE model the same as it was used for the active flexion FE model in [1].
• ### Efficiency of magnetic hyperthermia in the presence of rotating and static fields(1706.07426)

July 12, 2018 cond-mat.soft, physics.med-ph
Single-domain ferromagnetic nanoparticle systems can be used to transfer energy from a time-dependent magnetic field into their environment. This local heat generation, i.e., magnetic hyperthermia, receives applications in cancer therapy which requires the enhancement of the energy loss. A possible way to improve the efficiency is to chose a proper type of applied field, e.g., a rotating instead of an oscillating one. The latter case is very well studied and there is an increasing interest in the literature to investigate the former although it is still unclear under which circumstances the rotating applied field can be more favourable than the oscillating one. The goal of this work is to incorporate the presence of a static field and to perform a systematic study of the non-linear dynamics of the magnetisation in the framework of the deterministic Landau-Lifshitz-Gilbert equation in order to calculate energy losses. Two cases are considered: the static field is either assumed to be perpendicular to the plane of rotation or situated in the plane of rotation. In the latter case a significant increase in the energy loss/cycle is observed if the magnitudes of the static and the rotating fields have a certain ratio (e.g. it should be one for isotropic nanoparticles). It can be used to "super-localise" the heat transfer: in case of an inhomogeneous applied static field, tissues are heated up only where the magnitudes of the static and rotating fields reach the required ratio.
• ### Graphene quantum dots prevent alpha-synucleinopathy in Parkinson's disease(1710.07213)

While the emerging evidence indicates that the pathogenesis of Parkinson's disease (PD) is strongly correlated to the accumulation of alpha-synuclein ({\alpha}-syn) aggregates, there has been no clinical success in anti-aggregation agents for the disease to date. Here we show that graphene quantum dots (GQDs) exhibit anti-amyloid activity via direct interaction with {\alpha}-syn. Employing biophysical, biochemical, and cell-based assays as well as molecular dynamics (MD) simulation, we find that GQDs have notable potency in not only inhibiting fibrillization of {\alpha}-syn but also disaggregating mature fibrils in a time-dependent manner. Remarkably, GQDs rescue neuronal death and synaptic loss, reduce Lewy body (LB)/Lewy neurite (LN) formation, ameliorate mitochondrial dysfunctions, and prevent neuron-to-neuron transmission of {\alpha}-syn pathology induced by {\alpha}-syn preformed fibrils (PFFs) in neurons. In addition, in vivo administration of GQDs protects against {\alpha}-syn PFFs-induced loss of dopamine neurons, LB/LN pathology, and behavioural deficits through the penetration of the blood-brain barrier (BBB). The finding that GQDs function as an anti-aggregation agent provides a promising novel therapeutic target for the treatment of PD and related {\alpha}-synucleinopathies.
• ### A Bayesian approach for energy-based estimation of acoustic aberrations in high intensity focused ultrasound treatment(1602.08080)

June 7, 2018 physics.med-ph, stat.AP
High intensity focused ultrasound is a non-invasive method for treatment of diseased tissue that uses a beam of ultrasound to generate heat within a small volume. A common challenge in application of this technique is that heterogeneity of the biological medium can defocus the ultrasound beam. Here we reduce the problem of refocusing the beam to the inverse problem of estimating the acoustic aberration due to the biological tissue from acoustic radiative force imaging data. We solve this inverse problem using a Bayesian framework with a hierarchical prior and solve the inverse problem using a Metropolis-within-Gibbs algorithm. The framework is tested using both synthetic and experimental datasets. We demonstrate that our approach has the ability to estimate the aberrations using small datasets, as little as 32 sonication tests, which can lead to significant speedup in the treatment process. Furthermore, our approach is compatible with a wide range of sonication tests and can be applied to other energy-based measurement techniques.
• ### Simple Auto-Calibrated Gradient Delay Estimation From Few Spokes Using Radial Intersections (RING)(1805.04334)

May 11, 2018 physics.med-ph
• ### Quantitative cone-beam CT reconstruction with polyenergetic scatter model fusion(1805.04411)

May 11, 2018 physics.med-ph
Scatter can account for large errors in cone-beam CT (CBCT) due to its wide field of view, and its complicated nature makes its compensation difficult. Iterative polyenergetic reconstruction algorithms offer the potential to provide quantitative imaging in CT, but they are usually incompatible with scatter contaminated measurements. In this work, we introduce a polyenergetic convolutional scatter model that is directly fused into the reconstruction process, and exploits information readily available at each iteration for a fraction of additional computational cost. We evaluate this method with numerical and real CBCT measurements, and show significantly enhanced electron density estimation and artefact mitigation over pre-calculated fast adaptive scatter kernel superposition (fASKS). We demonstrate our approach has two levels of benefit: reducing the bias introduced by estimating scatter prior to reconstruction; and adapting to the spectral and spatial properties of the specimen.
• ### Liquid Metal Enabled Electrobiology: A Generalized Easy Going Way to Tackle Disease Challenges(1805.04002)

May 10, 2018 physics.med-ph
In this article, a new conceptual biomedical engineering strategy to tackle modern disease challenges, termed as liquid metal enabled electrobiology, is proposed. This generalized and easy going way is based on the physiological fact that specially administrated electricity would induce a series of subsequent desired biological effects, either shortly, transitional or permanently. Owing high compliance within any part of the biological tissues, the liquid metal would aid to mold a pervasive way to treat physiological or psychological diseases. As highly conductive and non-toxic multifunctional flexible materials, such liquid metals (LMs) consist of the core to generate any requested electric treating fields (ETFields) which can adapt to various sites inside the human body. The basic mechanisms of electrobiology in delivering electricity to the target tissues and then inducing expected outputs for disease treatment are interpreted. The methods toward realizing soft and conformable electronics based on liquid metal are illustrated. Further, a group of typical disease challenges were taken to illustrate the basic strategies to perform liquid metal electrobiology therapy, which include but are not limited to: tissue electronics, brain disorder, immunotherapy, neural functional recovery, muscle stimulation, skin rejuvenation, cosmetology and dieting, artificial organs, cardiac pacing and cancer therapy etc. Some practical issues involved in the electrobiology for future disease therapy were discussed. Perspective along this direction to incubate an easy-going biomedical tool for health care was pointed out.
• ### An EEG pre-processing technique for the fast recognition of motor imagery movements(1805.03169)

May 8, 2018 physics.med-ph, eess.SP
In this paper we propose a new pre-processing technique of Electroencephalography (EEG) signals produced by motor imagery movements. This technique results to an accelerated determination of the imagery movement and the command to carry it out, within the time limits imposed by the requirements of brain-based real-time control of rehabilitation devices, making thus feasible to drive these devices according to patient's will. Based on event related desynchronization and synchronization (ERD/ERS) of motor imagery, the received patient signal is first subjected to the removal of environmental, system and interference noise which correspond to normal human activities such as eye-blinking and cardiac motion. Next, power and energy features of the processed signal are compared with the same features of classified signals from an available database and the class to which the processed signal belongs, is identified. The database classification is done off-line by using the SVM algorithm.
• ### Highly Scalable Image Reconstruction using Deep Neural Networks with Bandpass Filtering(1805.03300)

May 8, 2018 cs.CV, physics.med-ph
To increase the flexibility and scalability of deep neural networks for image reconstruction, a framework is proposed based on bandpass filtering. For many applications, sensing measurements are performed indirectly. For example, in magnetic resonance imaging, data are sampled in the frequency domain. The introduction of bandpass filtering enables leveraging known imaging physics while ensuring that the final reconstruction is consistent with actual measurements to maintain reconstruction accuracy. We demonstrate this flexible architecture for reconstructing subsampled datasets of MRI scans. The resulting high subsampling rates increase the speed of MRI acquisitions and enable the visualization rapid hemodynamics.
• ### Three-dimensional Simulation of Quantitative Ultrasound in Cancellous Bone using the Echographic Response of a Metallic Pin(1805.03083)

Degenerative discopathy is a common pathology which may require spine surgery. A metallic cylindrical pin is inserted into the vertebral body to maintain soft tissues and may be used as a reflector of ultrasonic wave to estimate bone density. The first aim of this paper is to validate a 3D model to simulate the ultrasonic propagation in a trabecular bone sample in which a metallic pin has been inserted. We also aim at determining the effect of changes of bone volume fraction (BV/TV) and of positioning errors on the quantitative ultrasound (QUS) parameters in this specific configuration. The approach consists in coupling finite difference time domain simulation with X-ray microcomputed tomography. The correlation coefficient between experimental and simulated speed of sound (SOS) (respectively broadband ultrasonic attenuation (BUA)) was equal to 0.90 (respectively 0.55). The results show a significant correlation of SOS with BV/TV (R = 0.82), while BUA values exhibit a non-linear behavior versus BV/TV. The orientation of the pin should be controlled with an accuracy of around 1{\textdegree} in order to obtain accurate results. The results indicate that using the ultrasonic wave reflected by a pin has a potential to estimate the bone density. SOS is more reliable than BUA due its lower sensitivity to the tilt angle.