Showing papers in "arXiv: Medical Physics in 2019"

SANDI: a compartment-based model for non-invasive apparent soma and neurite imaging by diffusion MRI

Abstract: This work introduces a compartment-based model for apparent soma and neurite density imaging (SANDI) using non-invasive diffusion-weighted MRI (DW-MRI). The existing conjecture in brain microstructure imaging trough DW-MRI presents water diffusion in white (WM) and grey (GM) matter as restricted diffusion in neurites, modelled by infinite cylinders of null radius embedded in the hindered extra-neurite water. The extra-neurite pool in WM corresponds to water in the extra-axonal space, but in GM it combines water in the extra-cellular space with water in soma. While several studies showed that this microstructure model successfully describe DW-MRI data in WM and GM at b >3 ms/{\mum^2}). Here we hypothesize that the unmodelled soma compartment may be responsible for this failure and propose SANDI as a new model of brain microstructure where soma (i.e. cell body of any brain cell type: from neuroglia to neurons) is explicitly included. We assess the effects of size and density of soma on the direction-averaged DW-MRI signal at high b values and the regime of validity of the model using numerical simulations and comparison with experimental data from mouse (bmax = 40 ms/{/mum^2}) and human (bmax = 10 ms/{\mum^2}) brain. We show that SANDI defines new contrasts representing new complementary information on the brain cyto- and myelo-architecture. Indeed, we show for the first-time maps from 25 healthy human subjects of MR soma and neurite signal fractions, that remarkably mirror contrasts of histological images of brain cyto- and myelo-architecture. Although still under validation, SANDI might provide new insight into tissue architecture by introducing a new set of biomarkers of potential great value for biomedical applications and pure neuroscience.

Convolutional Sparse Coding for Compressed Sensing CT Reconstruction

Abstract: Over the past few years, dictionary learning (DL)-based methods have been successfully used in various image reconstruction problems. However, traditional DL-based computed tomography (CT) reconstruction methods are patch-based and ignore the consistency of pixels in overlapped patches. In addition, the features learned by these methods always contain shifted versions of the same features. In recent years, convolutional sparse coding (CSC) has been developed to address these problems. In this paper, inspired by several successful applications of CSC in the field of signal processing, we explore the potential of CSC in sparse-view CT reconstruction. By directly working on the whole image, without the necessity of dividing the image into overlapped patches in DL-based methods, the proposed methods can maintain more details and avoid artifacts caused by patch aggregation. With predetermined filters, an alternating scheme is developed to optimize the objective function. Extensive experiments with simulated and real CT data were performed to validate the effectiveness of the proposed methods. Qualitative and quantitative results demonstrate that the proposed methods achieve better performance than several existing state-of-the-art methods.

Next Generation Radiogenomics Sequencing for Prediction of EGFR and KRAS Mutation Status in NSCLC Patients Using Multimodal Imaging and Machine Learning Approaches

Abstract: Aim: In the present work, we aimed to evaluate a comprehensive radiomics framework that enabled prediction of EGFR and KRAS mutation status in NSCLC cancer patients based on PET and CT multi-modalities radiomic features and machine learning (ML) algorithms. Methods: Our study involved 211 NSCLC cancer patient with PET and CTD images. More than twenty thousand radiomic features from different image-feature sets were extracted Feature value was normalized to obtain Z-scores, followed by student t-test students for comparison, high correlated features were eliminated and the False discovery rate (FDR) correction were performed Six feature selection methods and twelve classifiers were used to predict gene status in patient and model evaluation was reported on independent validation sets (68 patients). Results: The best predictive power of conventional PET parameters was achieved by SUVpeak (AUC: 0.69, P-value = 0.0002) and MTV (AUC: 0.55, P-value = 0.0011) for EGFR and KRAS, respectively. Univariate analysis of radiomics features improved prediction power up to AUC: 75 (q-value: 0.003, Short Run Emphasis feature of GLRLM from LOG preprocessed image of PET with sigma value 1.5) and AUC: 0.71 (q-value 0.00005, The Large Dependence Low Gray Level Emphasis from GLDM in LOG preprocessed image of CTD sigma value 5) for EGFR and KRAS, respectively. Furthermore, the machine learning algorithm improved the perdition power up to AUC: 0.82 for EGFR (LOG preprocessed of PET image set with sigma 3 with VT feature selector and SGD classifier) and AUC: 0.83 for KRAS (CT image set with sigma 3.5 with SM feature selector and SGD classifier). Conclusion: We demonstrated that radiomic features extracted from different image-feature sets could be used for EGFR and KRAS mutation status prediction in NSCLC patients, and showed that they have more predictive power than conventional imaging parameters.

Review of Cost Reduction Methods in Photoacoustic Computed Tomography

Abstract: Photoacoustic Computed Tomography (PACT) is a major configuration of photoacoustic imaging, a hybrid noninvasive modality for both functional and molecular imaging. PACT has rapidly gained importance in the field of biomedical imaging due to superior performance as compared to conventional optical imaging counterparts. However, the overall cost of developing a PACT system is one of the challenges towards clinical translation of this novel technique. The cost of a typical commercial PACT system originates from optical source, ultrasound detector, and data acquisition unit. With growing applications of photoacoustic imaging, there is a tremendous demand towards reducing its cost. In this review article, we have discussed various approaches to reduce the overall cost of a PACT system, and provided a cost estimation to build a low-cost PACT system.

Accelerated Nuclear Magnetic Resonance Spectroscopy with Deep Learning

Abstract: Nuclear magnetic resonance (NMR) spectroscopy serves as an indispensable tool in chemistry and biology but often suffers from long experimental time. We present a proof-of-concept of application of deep learning and neural network for high-quality, reliable, and very fast NMR spectra reconstruction from limited experimental data. We show that the neural network training can be achieved using solely synthetic NMR signal, which lifts the prohibiting demand for a large volume of realistic training data usually required in the deep learning approach.

Objective Study of Sensor Relevance for Automatic Cough Detection

Abstract: The development of a system for the automatic, objective and reliable detection of cough events is a need underlined by the medical literature for years. The benefit of such a tool is clear as it would allow the assessment of pathology severity in chronic cough diseases. Even though some approaches have recently reported solutions achieving this task with a relative success, there is still no standardization about the method to adopt or the sensors to use. The goal of this paper is to study objectively the performance of several sensors for cough detection: ECG, thermistor, chest belt, accelerometer, contact and audio microphones. Experiments are carried out on a database of 32 healthy subjects producing, in a confined room and in three situations, voluntary cough at various volumes as well as other event categories which can possibly lead to some detection errors: background noise, forced expiration, throat clearing, speech and laugh. The relevance of each sensor is evaluated at three stages: mutual information conveyed by the features, ability to discriminate at the frame level cough from these latter other sources of ambiguity, and ability to detect cough events. In this latter experiment, with both an averaged sensitivity and specificity of about 94.5%, the proposed approach is shown to clearly outperform the commercial Karmelsonix system which achieved a specificity of 95.3% and a sensitivity of 64.9%.

Overview of image-to-image translation by use of deep neural networks: denoising, super-resolution, modality conversion, and reconstruction in medical imaging.

Abstract: Since the advent of deep convolutional neural networks (DNNs), computer vision has seen an extremely rapid progress that has led to huge advances in medical imaging. This article does not aim to cover all aspects of the field but focuses on a particular topic, image-to-image translation. Although the topic may not sound familiar, it turns out that many seemingly irrelevant applications can be understood as instances of image-to-image translation. Such applications include (1) noise reduction, (2) super-resolution, (3) image synthesis, and (4) reconstruction. The same underlying principles and algorithms work for various tasks. Our aim is to introduce some of the key ideas on this topic from a uniform point of view. We introduce core ideas and jargon that are specific to image processing by use of DNNs. Having an intuitive grasp of the core ideas of and a knowledge of technical terms would be of great help to the reader for understanding the existing and future applications. Most of the recent applications which build on image-to-image translation are based on one of two fundamental architectures, called pix2pix and CycleGAN, depending on whether the available training data are paired or unpaired. We provide computer codes which implement these two architectures with various enhancements. Our codes are available online with use of the very permissive MIT license. We provide a hands-on tutorial for training a model for denoising based on our codes. We hope that this article, together with the codes, will provide both an overview and the details of the key algorithms, and that it will serve as a basis for the development of new applications.

A Method of Rapid Quantification of Patient-Specific Organ Dose for CT Using Coupled Deep-Learning based Multi-Organ Segmentation and GPU-accelerated Monte Carlo Dose Computing

Abstract: Purpose: This paper describes a new method to apply deep-learning algorithms for automatic segmentation of radiosensitive organs from 3D tomographic CT images before computing organ doses using a GPU-based Monte Carlo code. Methods: A deep convolutional neural network (CNN) for organ segmentation is trained to automatically delineate radiosensitive organs from CT. With a GPU-based Monte Carlo dose engine (ARCHER) to derive CT dose of a phantom made from a subject's CT scan, we are then able to compute the patient-specific CT dose for each of the segmented organs. The developed tool is validated by using Relative Dose Error (RDE) against the organ doses calculated by ARCHER with manual segmentation performed by radiologists. The dose computation results are also compared against organ doses from population-average phantoms to demonstrate the improvement achieved by using the developed tool. In this study, two datasets were used: The Lung CT Segmentation Challenge 2017 (LCTSC) dataset, which contains 60 thoracic CT scan patients each with 5 segmented organs, and the Pancreas-CT (PCT) dataset, which contains 43 abdominal CT scan patients each with 8 segmented organs. Five-fold cross-validation of the new method is performed on both datasets. Results: Comparing with the traditional organ dose evaluation method that based on population-average phantom, our proposed method achieved the smaller RDE range on all organs with -4.3%~1.5% vs -31.5%~33.9% (lung), -7.0%~2.3% vs -15.2%~125.1% (heart), -18.8%~40.2% vs -10.3%~124.1% (esophagus) in the LCTSC dataset and -5.6%~1.6% vs -20.3%~57.4% (spleen), -4.5%~4.6% vs -19.5%~61.0% (pancreas), -2.3%~4.4% vs -37.8%~75.8% (left kidney), -14.9%~5.4% vs -39.9% ~14.6% (gall bladder), -0.9%~1.6% vs -30.1%~72.5% (liver), and -23.0%~11.1% vs -52.5%~-1.3% (stomach) in the PCT dataset.

A computational model of radiolytic oxygen depletion during FLASH irradiation and its effect on the oxygen enhancement ratio

Abstract: Recent results from animal irradiation studies have rekindled interest in the potential of ultra-high dose rate irradiation (also known as FLASH) for reducing normal tissue toxicity. However, despite mounting evidence of a "FLASH effect", a mechanism has yet to be elucidated. This article hypothesizes that the radioprotecting effect of FLASH irradiation could be due to the specific sparing of hypoxic stem cell niches, which have been identified in several organs including the bone marrow and the brain. To explore this hypothesis, a new computational model is presented that frames transient radiolytic oxygen depletion (ROD) during FLASH irradiation in terms of its effect on the oxygen enhancement ratio (OER). The model takes into consideration oxygen diffusion through the tissue, its consumption by metabolic cells, and its radiolytic depletion to estimate the relative decrease in radiosensitivity of cells receiving FLASH irradiation. Based on this model, several predictions are made that could be tested in future experiments: (1) the FLASH effect should gradually disappear as the radiation pulse duration is increased from <1s to 10 s; (2) dose should be deposited using the smallest number of radiation pulses to achieve the greatest FLASH effect; (3) a FLASH effect should only be observed in cells that are already hypoxic at the time of irradiation; and (4) changes in capillary oxygen tension (increase or decrease) should diminish the FLASH effect.

Non-Invasive Fuhrman Grading of Clear Cell Renal Cell Carcinoma Using Computed Tomography Radiomics Features and Machine Learning

Abstract: Purpose: To identify optimal classification methods for computed tomography (CT) radiomics-based preoperative prediction of clear cells renal cell carcinoma (ccRCC) grade. Methods and material: Seventy one ccRCC patients were included in the study. Three image preprocessing techniques (Laplacian of Gaussian, wavelet filter, and discretization of the intensity values) were applied on tumor volumes. In total, 2530 radiomics features (tumor shape and size, intensity statistics, and texture) were extracted from each segmented tumor volume. Univariate analysis was performed to assess the association of each feature with the histological condition. In the case of multivariate analysis, the following was implemented: three feature selection including the least absolute shrinkage and selection operator (LASSO), students t-test and minimum Redundancy Maximum Relevance (mRMR) algorithms. These selected features were then used to construct three classification models (SVM, random forest, and logistic regression) to discriminate the high from low-grade ccRCC at nephrectomy. Lastly, multivariate model performance was evaluated on the bootstrapped validation cohort using the area under receiver operating characteristic curve (AUC). Results: Univariate analysis demonstrated that among different image sets, 128 bin discretized images have statistically significant different (q-value < 0.05) texture parameters with a mean of AUC 0.74 (q-value < 0.05). The three ML-based classifier shows proficient discrimination of the high from low-grade ccRCC. The AUC was 0.78 in logistic regression, 0.62 in random forest, and 0.83 in SVM model, respectively. Conclusion: Radiomics features can be a useful and promising non-invasive method for preoperative evaluation of ccRCC Fuhrman grades. Key words: RCC, Radiomics, Machine Learning, Computed Tomography

CBCT-to-CT synthesis with a single neural network for head-and-neck, lung and breast cancer adaptive radiotherapy

Abstract: Purpose: CBCT-based adaptive radiotherapy requires daily images for accurate dose calculations. This study investigates the feasibility of applying a single convolutional network to facilitate CBCT-to-CT synthesis for head-and-neck, lung, and breast cancer patients. Methods: Ninety-nine patients diagnosed with head-and-neck, lung or breast cancer undergoing radiotherapy with CBCT-based position verification were included in this study. CBCTs were registered to planning CTs according to clinical procedures. Three cycle-consistent generative adversarial networks (cycle-GANs) were trained in an unpaired manner on 15 patients per anatomical site generating synthetic-CTs (sCTs). Another network was trained with all the anatomical sites together. Performances of all four networks were compared and evaluated for image similarity against rescan CT (rCT). Clinical plans were recalculated on CT and sCT and analysed through voxel-based dose differences and {\gamma}-analysis. Results: A sCT was generated in 10 seconds. Image similarity was comparable between models trained on different anatomical sites and a single model for all sites. Mean dose differences 95% were achieved for all sites. Conclusions: Cycle-GAN reduced CBCT artefacts and increased HU similarity to CT, enabling sCT-based dose calculations. The speed of the network can facilitate on-line adaptive radiotherapy using a single network for head-and-neck, lung and breast cancer patients.

A Feasibility Study on Deep Learning-Based Radiotherapy Dose Calculation

Abstract: Purpose: Various dose calculation algorithms are available for radiation therapy for cancer patients. However, these algorithms are faced with the tradeoff between efficiency and accuracy. The fast algorithms are generally less accurate, while the accurate dose engines are often time consuming. In this work, we try to resolve this dilemma by exploring deep learning (DL) for dose calculation. Methods: We developed a new radiotherapy dose calculation engine based on a modified Hierarchically Densely Connected U-net (HD U-net) model and tested its feasibility with prostate intensity-modulated radiation therapy (IMRT) cases. Mapping from an IMRT fluence map domain to a 3D dose domain requires a deep neural network of complicated architecture and a huge training dataset. To solve this problem, we first project the fluence maps to the dose domain using a modified ray-tracing algorithm, and then we use the HD U-net to map the ray-tracing dose distribution into an accurate dose distribution calculated using a collapsed cone convolution/superposition (CS) algorithm. Results: It takes about one second to compute a 3D dose distribution for a typical 7-field prostate IMRT plan, which can be further reduced to achieve real-time dose calculation by optimizing the network. For all eight testing patients, evaluation with Gamma Index and various clinical goals for IMRT optimization shows that the DL dose distributions are clinically identical to the CS dose distributions. Conclusions: We have shown the feasibility of using DL for calculating radiotherapy dose distribution with high accuracy and efficiency.

Deep Model with Siamese Network for Viability and Necrosis Tumor Assessment in Osteosarcoma

Abstract: Osteosarcoma is the most common primary malignant bone tumor, which has high mortality due to easy lung metastasis. Osteosarcoma is a highly anaplastic, pleomorphic tumor with a variety of tumor cell morphology, including fusiform, oval, epithelial, lymphocyte like, small round, transparent cells, etc. Due to the multiple patterns of osteosarcoma cell morphology, pathologists have differences in the classification (viable tumor, necrotic tumor, non-tumor) of osteosarcoma. Therefore, automatic and accurate recognition algorithms can help pathologists greatly reduce time and improve diagnostic accuracy. In recent years, deep learning technology has made great progress in the field of natural images and medical images, and has achieved excellent results beyond human performance in classification. In this paper, we propose a Deep Model with Siamese Network (DS-Net) for automatic classification in Hematoxylin and Eosin (H&E) stained osteosarcoma histology images.

Estimation of blood oxygenation with learned spectral decoloring for quantitative photoacoustic imaging (LSD-qPAI).

Abstract: One of the main applications of photoacoustic (PA) imaging is the recovery of functional tissue properties, such as blood oxygenation (sO2). This is typically achieved by linear spectral unmixing of relevant chromophores from multispectral photoacoustic images. Despite the progress that has been made towards quantitative PA imaging (qPAI), most sO2 estimation methods yield poor results in realistic settings. In this work, we tackle the challenge by employing learned spectral decoloring for quantitative photoacoustic imaging (LSD-qPAI) to obtain quantitative estimates for blood oxygenation. LSD-qPAI computes sO2 directly from pixel-wise initial pressure spectra Sp0, which are vectors comprised of the initial pressure at the same spatial location over all recorded wavelengths. Initial results suggest that LSD-qPAI is able to obtain accurate sO2 estimates directly from multispectral photoacoustic measurements in silico and plausible estimates in vivo.

Maxwell-compensated design of asymmetric gradient waveforms for tensor-valued diffusion encoding

Abstract: Purpose: Asymmetric gradient waveforms are attractive for diffusion encoding due to their superior efficiency, however, the asymmetry may cause a residual gradient moment at the end of the encoding. Depending on the experiment setup, this residual moment may cause significant signal bias and image artifacts. The purpose of this study was to develop an asymmetric gradient waveform design for tensor-valued diffusion encoding that is not affected by concomitant gradient. Methods: The Maxwell index was proposed as a scalar invariant that captures the effect of concomitant gradients and was constrained in the numerical optimization to 100 (mT/m)$^2$ms to yield Maxwell-compensated waveforms. The efficacy of this design was tested in an oil phantom, and in a healthy human brain. For reference, waveforms from literature were included in the analysis. Simulations were performed to investigate if the design was valid for a wide range of experiments and if it could predict the signal bias. Results: Maxwell-compensated waveforms showed no signal bias in oil or in the brain. By contrast, several waveforms from literature showed gross signal bias. In the brain, the bias was large enough to markedly affect both signal and parameter maps, and the bias could be accurately predicted by theory. Conclusion: Constraining the Maxwell index in the optimization of asymmetric gradient waveforms yields efficient tensor-valued encoding with concomitant gradients that have a negligible effect on the signal. This waveform design is especially relevant in combination with strong gradients, long encoding times, thick slices, simultaneous multi-slice acquisition and large/oblique FOVs.

Patient-specific predictions of aneurysm growth and remodeling in the ascending thoracic aorta using the homogenized constrained mixture model

Abstract: In its permanent quest of mechanobiological homeostasis, our vascula-ture significantly adapts across multiple length and time scales in various physiological and pathological conditions. Computational modeling of vascular growth and remodeling (G\&R) has significantly improved our insights of the mechanobio-logical processes of diseases such as hypertension or aneurysms. However, patient-specific computational modeling of ascending thoracic aortic aneurysm (ATAA) evolution, based on finite-element models (FEM), remains a challenging scientific problem with rare contributions, despite the major significance of this topic of research. Challenges are related to complex boundary conditions and geometries combined with layer-specific G\&R responses. To address these challenges, in the current paper, we employed the constrained mixture model (CMM) to model the arterial wall as a mixture of different constituents such as elastin, collagen fiber families and smooth muscle cells (SMCs). Implemented in Abaqus as a UMAT, this first patient-specific CMM-based FEM of G\&R in human ATAA was first validated for canonical problems such as single-layer thick-wall cylindrical and bi-layer thick-wall toric arterial geometries. Then it was used to predict ATAA evolution for a patient-specific aortic geometry, showing that the typical shape of an ATAA can be simply produced by elastin proteolysis localized in regions of deranged hemodymanics. The results indicate a transfer of stress to the adventitia by elastin loss and continuous adaptation of the stress distribution due to change of ATAA shape. Moreover, stress redistribution leads to collagen deposition where the maximum elastin mass is lost, which in turn leads to stiffening of the arterial wall. As future work, the predictions of this G\&R framework will be validated on datasets of patient-specific ATAA geometries followed up over a significant number of years.

A time-dependent diffusion MRI signature of axon caliber variations and beading

Abstract: MRI provides a unique non-invasive window into the brain, yet is limited to millimeter resolution, orders of magnitude coarser than cell dimensions. Here we show that diffusion MRI is sensitive to the micrometer-scale variations in axon caliber or pathological beading, by identifying a signature power-law diffusion time-dependence of the along-fiber diffusion coefficient. We observe this signature in human brain white matter, and uncover its origins by Monte Carlo simulations in realistic substrates from 3d electron microscopy of mouse corpus callosum. Simulations reveal that the time-dependence originates from axon caliber variation, rather than from mitochondria or axonal undulations. We report a decreased amplitude of time-dependence in multiple sclerosis lesions, illustrating the sensitivity of our method to axonal beading in a plethora of neurodegenerative disorders. This specificity to microstructure offers an exciting possibility of bridging across scales to image cellular-level pathology with a clinically feasible MRI technique.

Flexible and Comprehensive Patient-Specific Mitral Valve Silicone Models with Chordae Tendinae Made From 3D-Printable Molds.

Abstract: Given the multitude of challenges surgeons face during mitral valve repair surgery, they should have a high confidence in handling of instruments and in the application of surgical techniques before they enter the operating room. Unfortunately, opportunities for surgical training of minimally-invasive repair are very limited, leading to a situation where most surgeons undergo a steep learning curve while operating the first patients. In order to provide a realistic tool for surgical training, a commercial simulator was augmented by flexible patient-specific mitral valve replica. In an elaborated production pipeline, finalized after many optimization cycles, models were segmented from 3D ultrasound and then 3D-printable molds were computed automatically and printed in rigid material, the lower part being water-soluble. After silicone injection, the silicone model was dissolved from the mold and anchored in the simulator. To our knowledge, our models are the first to comprise the full mitral valve apparatus, i.e. the annulus, leaflets, chordae tendineae and papillary muscles. Nine different valve molds were automatically created according to the proposed workflow (seven prolapsed valves and two valves with functional mitral insufficiency). From these mold geometries, 16 replica were manufactured. A material test revealed that Ecoflex\textsuperscript{TM} 00-30 is the most suitable material for leaflet-mimicking tissue out of seven mixtures. Production time was around 36h per valve. Twelve surgeons performed various surgical techniques, e.g. annuloplasty, neo-chordae implantation, triangular leaflet resection and assessed the realism of the valves very positively. The standardized production process guarantees a high anatomical recapitulation of the silicone valves to the segmented models and the ultrasound data...

Dry electrodes for bioimpedance measurements -- design, characterization and comparison

Abstract: Objective: Bioimpedance measurements are mostly performed utilizing gel electrodes to decrease the occurring electrode-skin impedance. Since in many measurement environments this kind of electrode is not appropriate, the usability of dry electrodes is analysed. Approach: The development of five different kinds of dry electrodes, including gold, stainless steel, carbon rubber and metallized textile as contact materials are proposed. All test electrodes are based on a circular printed circuit board as carrier and have the same contact surface dimensions. To compare the electrodes' characteristics, the occurring electrode-skin impedances are measured under variation of signal frequency, contact duration, contact pressure, placement position and subjects. Additionally, all measurements are performed with silver/silver chloride (Ag/AgCl) hydrogel electrodes for comparison purposes. Main results: The analysed parameters play a significant role regarding the electrode-skin impedance. Choosing a wise setup of these parameters can decrease the electrode-skin impedance of dry electrodes down to ranges of hydrogel electrodes and even below. Significance: The usage of dry electrodes is one of the most difficult challenges when transferring scientific measurement techniques to clinical environments or commercial products but it is indispensable for many applications like body composition measurements or prosthesis control.

Evaluation of principal component analysis image denoising on multi-exponential MRI relaxometry

Abstract: PURPOSE: Multi-exponential relaxometry is a powerful tool for characterizing tissue, but generally requires high image signal-to-noise ratio (SNR). This work evaluates the use of principal-component-analysis (PCA) denoising to mitigate these SNR demands and improve the precision of relaxometry measures. METHODS: PCA denoising was evaluated using both simulated and experimental MRI data. Bi-exponential transverse relaxation signals were simulated for a wide range of acquisition and sample parameters, and experimental data were acquired from three excised and fixed mouse brains. In both cases, standard relaxometry analysis was performed on both original and denoised image data, and resulting estimated signal parameters were compared. RESULTS: Denoising reduced the root-mean-square-error of parameters estimated from multi-exponential relaxometry by factors of 2 to 4x, for typical acquisition and sample parameters. Denoised images and subsequent parameter maps showed little or no signs of spatial artifact or loss of resolution. CONCLUSION: Experimental studies and simulations demonstrate that PCA denoising of MRI relaxometry data is an effective method of improving parameter precision without sacrificing image resolution. This simple yet important processing step thus paves the way for broader applicability of multi-exponential MRI relaxometry.

Speed-of-Sound Imaging using Diverging Waves

Abstract: Recent ultrasound imaging modalities based on ultrasound computed tomography indicate a huge potential to detect pathologies is tissue due to altered biomechanical properties. Especially the imaging of speed-of-sound (SoS) distribution in tissue has shown clinical promise and thus gained increasing attention in the field -- with several methods proposed based on transmission mode tomography. SoS imaging using conventional ultrasound (US) systems would be convenient and easy for clinical translation, but this requires using conventional US probes with single-sides tissue access and thus pulse-echo imaging sequences. Recent pulse-echo SoS imaging methods rely on plane wave (PW) insonifications, which is prone to strong aberration effects for non-homogeneous tissue composition. In this paper we propose to use diverging waves (DW) for SoS imaging and thus substantially improve the reconstruction of SoS distributions. We study this proposition by first plane wavefront aberrations compared to DW. We then present the sensitivity of both approaches to major parameterization choices on a set of simulated phantoms. Using the optimum parameter combination for each method for a given transducer model and imaging sequence, we analyze the SoS imaging performance comparatively between the two approaches. Results indicate that using DW instead of PW, the reconstruction accuracy improves substantially, by over 22% in reconstruction error (RMSE) and by 55% in contrast (CNR). We also demonstrate improvements in SoS reconstructions from an actual US acquisition of a breast phantom with tumor- and cyst-representative inclusions, with high and low SoS contrast, respectively.

Improvement on corrosion resistance and biocompability of ZK60 magnesium alloy by carboxyl ion implantation

Abstract: Magnesium alloys have been considered to be potential biocompatible metallic materials. Further improvement on the anti-corrosion is expected to make this type of materials more suitable for biomedical applications in the fields of orthopedics, cardiovascular surgery and others. In this paper, we introduce a method of carboxyl ion (COOH+) implantation to reduce the degradation of ZK60 Mg alloy and improve its functionality in physiological environment. X-ray photoelectron spectroscopy (XPS) and atomic force microscopy (AFM) experiments show the formation of a smooth layer containing carbaxylic group, carbonate, metal oxides and hydroxides on the ion implanted alloy surface. Corrosion experiments and in vitro cytotoxicity tests demonstrate that the ion implantation treatment can both reduce the corrosion rate and improve the biocompatibility of the alloy. The promising results indicate that organic functional group ion implantation may be a practical method of improving the biological and corrosion properties of magnesium alloys.

A survey of computational frameworks for solving the acoustic inverse problem in three-dimensional photoacoustic computed tomography

Abstract: Photoacoustic computed tomography (PACT), also known as optoacoustic tomography, is an emerging imaging technique that holds great promise for biomedical imaging. PACT is a hybrid imaging method that can exploit the strong endogenous contrast of optical methods along with the high spatial resolution of ultrasound methods. In its canonical form that is addressed in this article, PACT seeks to estimate the photoacoustically-induced initial pressure distribution within the object. Image reconstruction methods are employed to solve the acoustic inverse problem associated with the image formation process. When an idealized imaging scenario is considered, analytic solutions to the PACT inverse problem are available; however, in practice, numerous challenges exist that are more readily addressed within an optimization-based, or iterative, image reconstruction framework. In this article, the PACT image reconstruction problem is reviewed within the context of modern optimization-based image reconstruction methodologies. Imaging models that relate the measured photoacoustic wavefields to the sought-after object function are described in their continuous and discrete forms. The basic principles of optimization-based image reconstruction from discrete PACT measurement data are presented, which includes a review of methods for modeling the PACT measurement system response and other important physical factors. Non-conventional formulations of the PACT image reconstruction problem, in which acoustic parameters of the medium are concurrently estimated along with the PACT image, are also introduced and reviewed.

MD-Recon-Net: A Parallel Dual-Domain Convolutional Neural Network for Compressed Sensing MRI

Abstract: Compressed sensing magnetic resonance imaging (CS-MRI) is a theoretical framework that can accurately reconstruct images from undersampled k-space data with a much lower sampling rate than the one set by the classical Nyquist-Shannon sampling theorem. Therefore, CS-MRI can efficiently accelerate acquisition time and relieve the psychological burden on patients while maintaining high imaging quality. The problems with traditional CS-MRI reconstruction are solved by iterative numerical solvers, which usually suffer from expensive computational cost and the lack of accurate handcrafted priori. In this paper, inspired by deep learning's (DL's) fast inference and excellent end-to-end performance, we propose a novel cascaded convolutional neural network called MD-Recon-Net to facilitate fast and accurate MRI reconstruction. Especially, different from existing DL-based methods, which operate on single domain data or both domains in a certain order, our proposed MD-Recon-Net contains two parallel and interactive branches that simultaneously perform on k-space and spatial-domain data, exploring the latent relationship between k-space and the spatial domain. The simulated experimental results show that the proposed method not only achieves competitive visual effects to several state-of-the-art methods, but also outperforms other DL-based methods in terms of model scale and computational cost.

Towards Real-Time Respiratory Motion Prediction based on Long Short-Term Memory Neural Networks

Abstract: Radiation therapy of thoracic and abdominal tumors requires incorporating the respiratory motion into treatments. To precisely account for the patient respiratory motions and predict the respiratory signals, a generalized model for predictions of different types of respiratory motions is desired. The aim of this study is to explore the feasibility of developing a Long Short-Term Memory (LSTM)-based model for the respiratory signal prediction. To achieve that, 1703 sets of Real-Time Position Management data were collected from retrospective studies across three clinical institutions. These datasets were separated as the training, internal validity and external validity groups. 1187 datasets were used for model development and the remaining 516 datasets were used to test the generality power of the model. Furthermore, an exhaustive grid search was implemented to find the optimal hyper-parameters of the LSTM model. The hyper-parameters are the number of LSTM layers, the number of hidden units, the optimizer, the learning rate, the number of epochs, and the length of time lags. Our model achieved superior accuracy over conventional artificial neural network models: with a prediction window of 500ms, the LSTM model achieved an average relative Mean Absolute Error (MAE) of 0.037, an average Root Mean Square Error (RMSE) of 0.048, and a Maximum Error (ME) of 1.687 in the internal validity data, and an average relative MAE of 0.112, an average RMSE of 0.139 and an ME of 1.811 in the external validity data. Compared to the LSTM model trained with default hyper-parameters, the MAE of the optimized model results decreased by 20%, indicating the importance of tuning the hyper-parameters of LSTM models to obtain superior accuracy. This study demonstrates the potential of deep LSTM models for the respiratory signal prediction and illustrates the impacts of major hyper-parameters in LSTM models.

Artificial Intelligence for Pediatric Ophthalmology.

Abstract: PURPOSE OF REVIEW: Despite the impressive results of recent artificial intelligence (AI) applications to general ophthalmology, comparatively less progress has been made toward solving problems in pediatric ophthalmology using similar techniques. This article discusses the unique needs of pediatric ophthalmology patients and how AI techniques can address these challenges, surveys recent applications of AI to pediatric ophthalmology, and discusses future directions in the field. RECENT FINDINGS: The most significant advances involve the automated detection of retinopathy of prematurity (ROP), yielding results that rival experts. Machine learning (ML) has also been successfully applied to the classification of pediatric cataracts, prediction of post-operative complications following cataract surgery, detection of strabismus and refractive error, prediction of future high myopia, and diagnosis of reading disability via eye tracking. In addition, ML techniques have been used for the study of visual development, vessel segmentation in pediatric fundus images, and ophthalmic image synthesis. SUMMARY: AI applications could significantly benefit clinical care for pediatric ophthalmology patients by optimizing disease detection and grading, broadening access to care, furthering scientific discovery, and improving clinical efficiency. These methods need to match or surpass physician performance in clinical trials before deployment with patients. Due to widespread use of closed-access data sets and software implementations, it is difficult to directly compare the performance of these approaches, and reproducibility is poor. Open-access data sets and software implementations could alleviate these issues, and encourage further AI applications to pediatric ophthalmology. KEYWORDS: pediatric ophthalmology, machine learning, artificial intelligence, deep learning

Flexible water excitation for fat-free MRI at 3 Tesla using lipid insensitive binomial off-resonant RF excitation (LIBRE) pulses

Abstract: Purpose: To develop a robust and flexible low power water excitation pulse that enables effective fat suppression at high magnetic field strength. Methods: A water excitation method that uses spatially non-selective pulses was optimized in numerical simulations, and implemented and tested in phantoms and healthy volunteers at 3T. The lipid insensitive binomial off-resonant excitation (LIBRE) pulse comprises two low power rectangular sub-pulses that have a variable frequency offset, phase offset and duration. The capability and extent of LIBRE fat suppression was quantitatively compared with conventional fat saturation (FS) and water excitation (WE) techniques. Results: LIBRE enables simultaneous water excitation and near complete fat suppression in large volumes at 3T as demonstrated by numerical simulations, and experiments. In phantoms and in human subjects, the frequency responses matched well with those from the numerical simulation. Comparing FS and WE, LIBRE demonstrated an improved robustness to magnetic field inhomogeneities, and a much more effectively suppressed fat signal. This applied for a range of pulse durations and pulses as short as 1.4 ms. Conclusion: A flexible water excitation method was developed that shows robust, near complete fat suppression at 3T.

Monitoring dynamic collagen reorganization during skin stretching with fast polarization-resolved second harmonic generation imaging

Abstract: The mechanical properties of biological tissues are strongly correlated to the specific distribution of their collagen fibers. Monitoring the dynamic reorganization of the collagen network during mechanical stretching is however a technical challenge because it requires mapping orientation of collagen fibers in a thick and deforming sample. In this work, a fast polarization-resolved SHG microscope is implemented to map collagen orientation during mechanical assays. This system is based on line-to-line switching of polarization using an electro-optical modulator and works in epidetection geometry. After proper calibration, it successfully highlights the collagen dynamic alignment along the traction direction in ex vivo murine skin dermis. This microstructure reorganization is quantified by the entropy of the collagen orientation distribution as a function of the stretch ratio. It exhibits a linear behavior, whose slope is measured with a good accuracy. This approach can be generalized to probe a variety of dynamic processes in thick tissues.

Single-shot x-ray speckle-based imaging of a single-material object

Abstract: We develop a means for speckle-based phase imaging of the projected thickness of a single-material object, under the assumption of illumination by spatially random time-independent x-ray speckles. These speckles are generated by passing x rays through a suitable spatially random mask. The method makes use of a single image obtained in the presence of the object, which serves to deform the illuminating speckle field relative to a reference speckle field (which only needs to be measured once) obtained in the presence of the mask and the absence of the object. The method implicitly rather than explicitly tracks speckles, and utilizes the transport-of-intensity equation to give a closed-form solution to the inverse problem of determining the complex transmission function of the object. Implementation using x-ray synchrotron data shows the method to be robust and efficient with respect to noise. Applications include x-ray phase--amplitude radiography and tomography, as well as time-dependent imaging of dynamic and radiation-sensitive samples using low-flux sources.

Generation of abdominal synthetic CTs from 0.35T MR images using generative adversarial networks for MR-only liver radiotherapy

Abstract: Electron density maps must be accurately estimated to achieve valid dose calculation in MR-only radiotherapy. The goal of this study is to assess whether two deep learning models, the conditional generative adversarial network (cGAN) and the cycle-consistent generative adversarial network (cycleGAN), can generate accurate abdominal synthetic CT (sCT) images from 0.35T MR images for MR-only liver radiotherapy. A retrospective study was performed using CT images and 0.35T MR images of 12 patients with liver (n=8) and non-liver abdominal (n=4) cancer. CT images were deformably registered to the corresponding MR images to generate deformed CT (dCT) images for treatment planning. Both cGAN and cycleGAN were trained using MR and dCT transverse slices. Four-fold cross-validation testing was conducted to generate sCT images for all patients. The HU prediction accuracy was evaluated by voxel-wise similarity metric between each dCT and sCT image for all 12 patients. dCT-based and sCT-based dose distributions were compared using gamma and dose-volume histogram (DVH) metric analysis for 8 liver patients. sCTcycleGAN achieved the average mean absolute error (MAE) of 94.1 HU, while sCTcGAN achieved 89.8 HU. In both models, the average gamma passing rates within all volumes of interest were higher than 95% using a 2%, 2 mm criterion, and 99% using a 3%, 3 mm criterion. The average differences in the mean dose and DVH metrics were within +/-0.6% for the planning target volume and within +/-0.15% for evaluated organs in both models. Results demonstrated that abdominal sCT images generated by both cGAN and cycleGAN achieved accurate dose calculation for 8 liver radiotherapy plans. sCTcGAN images had smaller average MAE and achieved better dose calculation accuracy than sCTcyleGAN images. More abdominal patients will be enrolled in the future to further evaluate two models.