-
RadRotator: 3D Rotation of Radiographs with Diffusion Models
Authors:
Pouria Rouzrokh,
Bardia Khosravi,
Shahriar Faghani,
Kellen L. Mulford,
Michael J. Taunton,
Bradley J. Erickson,
Cody C. Wyles
Abstract:
Transforming two-dimensional (2D) images into three-dimensional (3D) volumes is a well-known yet challenging problem for the computer vision community. In the medical domain, a few previous studies attempted to convert two or more input radiographs into computed tomography (CT) volumes. Following their effort, we introduce a diffusion model-based technology that can rotate the anatomical content o…
▽ More
Transforming two-dimensional (2D) images into three-dimensional (3D) volumes is a well-known yet challenging problem for the computer vision community. In the medical domain, a few previous studies attempted to convert two or more input radiographs into computed tomography (CT) volumes. Following their effort, we introduce a diffusion model-based technology that can rotate the anatomical content of any input radiograph in 3D space, potentially enabling the visualization of the entire anatomical content of the radiograph from any viewpoint in 3D. Similar to previous studies, we used CT volumes to create Digitally Reconstructed Radiographs (DRRs) as the training data for our model. However, we addressed two significant limitations encountered in previous studies: 1. We utilized conditional diffusion models with classifier-free guidance instead of Generative Adversarial Networks (GANs) to achieve higher mode coverage and improved output image quality, with the only trade-off being slower inference time, which is often less critical in medical applications; and 2. We demonstrated that the unreliable output of style transfer deep learning (DL) models, such as Cycle-GAN, to transfer the style of actual radiographs to DRRs could be replaced with a simple yet effective training transformation that randomly changes the pixel intensity histograms of the input and ground-truth imaging data during training. This transformation makes the diffusion model agnostic to any distribution variations of the input data pixel intensity, enabling the reliable training of a DL model on input DRRs and applying the exact same model to conventional radiographs (or DRRs) during inference.
△ Less
Submitted 19 April, 2024;
originally announced April 2024.
-
CONFLARE: CONFormal LArge language model REtrieval
Authors:
Pouria Rouzrokh,
Shahriar Faghani,
Cooper U. Gamble,
Moein Shariatnia,
Bradley J. Erickson
Abstract:
Retrieval-augmented generation (RAG) frameworks enable large language models (LLMs) to retrieve relevant information from a knowledge base and incorporate it into the context for generating responses. This mitigates hallucinations and allows for the updating of knowledge without retraining the LLM. However, RAG does not guarantee valid responses if retrieval fails to identify the necessary informa…
▽ More
Retrieval-augmented generation (RAG) frameworks enable large language models (LLMs) to retrieve relevant information from a knowledge base and incorporate it into the context for generating responses. This mitigates hallucinations and allows for the updating of knowledge without retraining the LLM. However, RAG does not guarantee valid responses if retrieval fails to identify the necessary information as the context for response generation. Also, if there is contradictory content, the RAG response will likely reflect only one of the two possible responses. Therefore, quantifying uncertainty in the retrieval process is crucial for ensuring RAG trustworthiness. In this report, we introduce a four-step framework for applying conformal prediction to quantify retrieval uncertainty in RAG frameworks. First, a calibration set of questions answerable from the knowledge base is constructed. Each question's embedding is compared against document embeddings to identify the most relevant document chunks containing the answer and record their similarity scores. Given a user-specified error rate (α), these similarity scores are then analyzed to determine a similarity score cutoff threshold. During inference, all chunks with similarity exceeding this threshold are retrieved to provide context to the LLM, ensuring the true answer is captured in the context with a (1-α) confidence level. We provide a Python package that enables users to implement the entire workflow proposed in our work, only using LLMs and without human intervention.
△ Less
Submitted 3 April, 2024;
originally announced April 2024.
-
Toward Clinically Trustworthy Deep Learning: Applying Conformal Prediction to Intracranial Hemorrhage Detection
Authors:
Cooper Gamble,
Shahriar Faghani,
Bradley J. Erickson
Abstract:
As deep learning (DL) continues to demonstrate its ability in radiological tasks, it is critical that we optimize clinical DL solutions to include safety. One of the principal concerns in the clinical adoption of DL tools is trust. This study aims to apply conformal prediction as a step toward trustworthiness for DL in radiology. This is a retrospective study of 491 non-contrast head CTs from the…
▽ More
As deep learning (DL) continues to demonstrate its ability in radiological tasks, it is critical that we optimize clinical DL solutions to include safety. One of the principal concerns in the clinical adoption of DL tools is trust. This study aims to apply conformal prediction as a step toward trustworthiness for DL in radiology. This is a retrospective study of 491 non-contrast head CTs from the CQ500 dataset, in which three senior radiologists annotated slices containing intracranial hemorrhage (ICH). The dataset was split into definite and challenging subsets, where challenging images were defined to those in which there was disagreement among readers. A DL model was trained on 146 patients (10,815 slices) from the definite data (training dataset) to perform ICH localization and classification for five classes of ICH. To develop an uncertainty-aware DL model, 1,546 cases of the definite data (calibration dataset) was used for Mondrian conformal prediction (MCP). The uncertainty-aware DL model was tested on 8,401 definite and challenging cases to assess its ability to identify challenging cases. After the MCP procedure, the model achieved an F1 score of 0.920 for ICH classification on the test dataset. Additionally, it correctly identified 6,837 of the 6,856 total challenging cases as challenging (99.7% accuracy). It did not incorrectly label any definite cases as challenging. The uncertainty-aware ICH detector performs on par with state-of-the-art models. MCP's performance in detecting challenging cases demonstrates that it is useful in automated ICH detection and promising for trustworthiness in radiological DL.
△ Less
Submitted 15 January, 2024;
originally announced January 2024.
-
Physics-Informed Deep Learning of Rate-and-State Fault Friction
Authors:
Cody Rucker,
Brittany A. Erickson
Abstract:
Direct observations of earthquake nucleation and propagation are few and yet the next decade will likely see an unprecedented increase in indirect, surface observations that must be integrated into modeling efforts. Machine learning (ML) excels in the presence of large data and is an actively growing field in seismology. However, not all ML methods incorporate rigorous physics, and purely data-dri…
▽ More
Direct observations of earthquake nucleation and propagation are few and yet the next decade will likely see an unprecedented increase in indirect, surface observations that must be integrated into modeling efforts. Machine learning (ML) excels in the presence of large data and is an actively growing field in seismology. However, not all ML methods incorporate rigorous physics, and purely data-driven models can predict physically unrealistic outcomes due to observational bias or extrapolation. Our work focuses on the recently emergent Physics-Informed Neural Network (PINN), which seamlessly integrates data while ensuring that model outcomes satisfy rigorous physical constraints. In this work we develop a multi-network PINN for both the forward problem as well as for direct inversion of nonlinear fault friction parameters, constrained by the physics of motion in the solid Earth, which have direct implications for assessing seismic hazard. We present the computational PINN framework for strike-slip faults in 1D and 2D subject to rate-and-state friction. Initial and boundary conditions define the data on which the PINN is trained. While the PINN is capable of approximating the solution to the governing equations to low-errors, our primary interest lies in the network's capacity to infer friction parameters during the training loop. We find that the network for the parameter inversion at the fault performs much better than the network for material displacements to which it is coupled. Additional training iterations and model tuning resolves this discrepancy, enabling a robust surrogate model for solving both forward and inverse problems relevant to seismic faulting.
△ Less
Submitted 14 December, 2023;
originally announced December 2023.
-
MedYOLO: A Medical Image Object Detection Framework
Authors:
Joseph Sobek,
Jose R. Medina Inojosa,
Betsy J. Medina Inojosa,
S. M. Rassoulinejad-Mousavi,
Gian Marco Conte,
Francisco Lopez-Jimenez,
Bradley J. Erickson
Abstract:
Artificial intelligence-enhanced identification of organs, lesions, and other structures in medical imaging is typically done using convolutional neural networks (CNNs) designed to make voxel-accurate segmentations of the region of interest. However, the labels required to train these CNNs are time-consuming to generate and require attention from subject matter experts to ensure quality. For tasks…
▽ More
Artificial intelligence-enhanced identification of organs, lesions, and other structures in medical imaging is typically done using convolutional neural networks (CNNs) designed to make voxel-accurate segmentations of the region of interest. However, the labels required to train these CNNs are time-consuming to generate and require attention from subject matter experts to ensure quality. For tasks where voxel-level precision is not required, object detection models offer a viable alternative that can reduce annotation effort. Despite this potential application, there are few options for general purpose object detection frameworks available for 3-D medical imaging. We report on MedYOLO, a 3-D object detection framework using the one-shot detection method of the YOLO family of models and designed for use with medical imaging. We tested this model on four different datasets: BRaTS, LIDC, an abdominal organ Computed Tomography (CT) dataset, and an ECG-gated heart CT dataset. We found our models achieve high performance on commonly present medium and large-sized structures such as the heart, liver, and pancreas even without hyperparameter tuning. However, the models struggle with very small or rarely present structures.
△ Less
Submitted 7 June, 2024; v1 submitted 12 December, 2023;
originally announced December 2023.
-
Synthetically Enhanced: Unveiling Synthetic Data's Potential in Medical Imaging Research
Authors:
Bardia Khosravi,
Frank Li,
Theo Dapamede,
Pouria Rouzrokh,
Cooper U. Gamble,
Hari M. Trivedi,
Cody C. Wyles,
Andrew B. Sellergren,
Saptarshi Purkayastha,
Bradley J. Erickson,
Judy W. Gichoya
Abstract:
Chest X-rays (CXR) are essential for diagnosing a variety of conditions, but when used on new populations, model generalizability issues limit their efficacy. Generative AI, particularly denoising diffusion probabilistic models (DDPMs), offers a promising approach to generating synthetic images, enhancing dataset diversity. This study investigates the impact of synthetic data supplementation on th…
▽ More
Chest X-rays (CXR) are essential for diagnosing a variety of conditions, but when used on new populations, model generalizability issues limit their efficacy. Generative AI, particularly denoising diffusion probabilistic models (DDPMs), offers a promising approach to generating synthetic images, enhancing dataset diversity. This study investigates the impact of synthetic data supplementation on the performance and generalizability of medical imaging research. The study employed DDPMs to create synthetic CXRs conditioned on demographic and pathological characteristics from the CheXpert dataset. These synthetic images were used to supplement training datasets for pathology classifiers, with the aim of improving their performance. The evaluation involved three datasets (CheXpert, MIMIC-CXR, and Emory Chest X-ray) and various experiments, including supplementing real data with synthetic data, training with purely synthetic data, and mixing synthetic data with external datasets. Performance was assessed using the area under the receiver operating curve (AUROC). Adding synthetic data to real datasets resulted in a notable increase in AUROC values (up to 0.02 in internal and external test sets with 1000% supplementation, p-value less than 0.01 in all instances). When classifiers were trained exclusively on synthetic data, they achieved performance levels comparable to those trained on real data with 200%-300% data supplementation. The combination of real and synthetic data from different sources demonstrated enhanced model generalizability, increasing model AUROC from 0.76 to 0.80 on the internal test set (p-value less than 0.01). In conclusion, synthetic data supplementation significantly improves the performance and generalizability of pathology classifiers in medical imaging.
△ Less
Submitted 7 July, 2024; v1 submitted 15 November, 2023;
originally announced November 2023.
-
Coronal Heating as Determined by the Solar Flare Frequency Distribution Obtained by Aggregating Case Studies
Authors:
James Paul Mason,
Alexandra Werth,
Colin G. West,
Allison A. Youngblood,
Donald L. Woodraska,
Courtney Peck,
Kevin Lacjak,
Florian G. Frick,
Moutamen Gabir,
Reema A. Alsinan,
Thomas Jacobsen,
Mohammad Alrubaie,
Kayla M. Chizmar,
Benjamin P. Lau,
Lizbeth Montoya Dominguez,
David Price,
Dylan R. Butler,
Connor J. Biron,
Nikita Feoktistov,
Kai Dewey,
N. E. Loomis,
Michal Bodzianowski,
Connor Kuybus,
Henry Dietrick,
Aubrey M. Wolfe
, et al. (977 additional authors not shown)
Abstract:
Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms th…
▽ More
Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfvén waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, $α=2$ as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed $>$600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: pre-flare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that $α= 1.63 \pm 0.03$. This is below the critical threshold, suggesting that Alfvén waves are an important driver of coronal heating.
△ Less
Submitted 9 May, 2023;
originally announced May 2023.
-
Report of the Medical Image De-Identification (MIDI) Task Group -- Best Practices and Recommendations
Authors:
David A. Clunie,
Adam Flanders,
Adam Taylor,
Brad Erickson,
Brian Bialecki,
David Brundage,
David Gutman,
Fred Prior,
J Anthony Seibert,
John Perry,
Judy Wawira Gichoya,
Justin Kirby,
Katherine Andriole,
Luke Geneslaw,
Steve Moore,
TJ Fitzgerald,
Wyatt Tellis,
Ying Xiao,
Keyvan Farahani
Abstract:
This report addresses the technical aspects of de-identification of medical images of human subjects and biospecimens, such that re-identification risk of ethical, moral, and legal concern is sufficiently reduced to allow unrestricted public sharing for any purpose, regardless of the jurisdiction of the source and distribution sites. All medical images, regardless of the mode of acquisition, are c…
▽ More
This report addresses the technical aspects of de-identification of medical images of human subjects and biospecimens, such that re-identification risk of ethical, moral, and legal concern is sufficiently reduced to allow unrestricted public sharing for any purpose, regardless of the jurisdiction of the source and distribution sites. All medical images, regardless of the mode of acquisition, are considered, though the primary emphasis is on those with accompanying data elements, especially those encoded in formats in which the data elements are embedded, particularly Digital Imaging and Communications in Medicine (DICOM). These images include image-like objects such as Segmentations, Parametric Maps, and Radiotherapy (RT) Dose objects. The scope also includes related non-image objects, such as RT Structure Sets, Plans and Dose Volume Histograms, Structured Reports, and Presentation States. Only de-identification of publicly released data is considered, and alternative approaches to privacy preservation, such as federated learning for artificial intelligence (AI) model development, are out of scope, as are issues of privacy leakage from AI model sharing. Only technical issues of public sharing are addressed.
△ Less
Submitted 1 April, 2023; v1 submitted 18 March, 2023;
originally announced March 2023.
-
Multitask Brain Tumor Inpainting with Diffusion Models: A Methodological Report
Authors:
Pouria Rouzrokh,
Bardia Khosravi,
Shahriar Faghani,
Mana Moassefi,
Sanaz Vahdati,
Bradley J. Erickson
Abstract:
Despite the ever-increasing interest in applying deep learning (DL) models to medical imaging, the typical scarcity and imbalance of medical datasets can severely impact the performance of DL models. The generation of synthetic data that might be freely shared without compromising patient privacy is a well-known technique for addressing these difficulties. Inpainting algorithms are a subset of DL…
▽ More
Despite the ever-increasing interest in applying deep learning (DL) models to medical imaging, the typical scarcity and imbalance of medical datasets can severely impact the performance of DL models. The generation of synthetic data that might be freely shared without compromising patient privacy is a well-known technique for addressing these difficulties. Inpainting algorithms are a subset of DL generative models that can alter one or more regions of an input image while matching its surrounding context and, in certain cases, non-imaging input conditions. Although the majority of inpainting techniques for medical imaging data use generative adversarial networks (GANs), the performance of these algorithms is frequently suboptimal due to their limited output variety, a problem that is already well-known for GANs. Denoising diffusion probabilistic models (DDPMs) are a recently introduced family of generative networks that can generate results of comparable quality to GANs, but with diverse outputs. In this paper, we describe a DDPM to execute multiple inpainting tasks on 2D axial slices of brain MRI with various sequences, and present proof-of-concept examples of its performance in a variety of evaluation scenarios. Our model and a public online interface to try our tool are available at: https://github.com/Mayo-Radiology-Informatics-Lab/MBTI
△ Less
Submitted 30 March, 2023; v1 submitted 21 October, 2022;
originally announced October 2022.
-
A Rokhlin Lemma for Noninvertible Totally-Ordered Measure-Preserving Dynamical Systems
Authors:
Adam R. B. Erickson
Abstract:
Let $(X,\mathcal{F},μ,T)$ be a not necessarily invertible non-atomic measure-preserving dynamical system where the $σ$-algebra $\mathcal{F}$ is generated by the intervals according to some total order. The main result is that the classical Rokhlin lemma may be adapted to such a situation assuming a slight extension of aperiodicity. This result is compared to previous noninvertible versions of the…
▽ More
Let $(X,\mathcal{F},μ,T)$ be a not necessarily invertible non-atomic measure-preserving dynamical system where the $σ$-algebra $\mathcal{F}$ is generated by the intervals according to some total order. The main result is that the classical Rokhlin lemma may be adapted to such a situation assuming a slight extension of aperiodicity. This result is compared to previous noninvertible versions of the Rokhlin lemma.
△ Less
Submitted 25 March, 2022; v1 submitted 22 March, 2022;
originally announced March 2022.
-
SOUP-GAN: Super-Resolution MRI Using Generative Adversarial Networks
Authors:
Kuan Zhang,
Haoji Hu,
Kenneth Philbrick,
Gian Marco Conte,
Joseph D. Sobek,
Pouria Rouzrokh,
Bradley J. Erickson
Abstract:
There is a growing demand for high-resolution (HR) medical images in both the clinical and research applications. Image quality is inevitably traded off with the acquisition time for better patient comfort, lower examination costs, dose, and fewer motion-induced artifacts. For many image-based tasks, increasing the apparent resolution in the perpendicular plane to produce multi-planar reformats or…
▽ More
There is a growing demand for high-resolution (HR) medical images in both the clinical and research applications. Image quality is inevitably traded off with the acquisition time for better patient comfort, lower examination costs, dose, and fewer motion-induced artifacts. For many image-based tasks, increasing the apparent resolution in the perpendicular plane to produce multi-planar reformats or 3D images is commonly used. Single image super-resolution (SR) is a promising technique to provide HR images based on unsupervised learning to increase resolution of a 2D image, but there are few reports on 3D SR. Further, perceptual loss is proposed in the literature to better capture the textual details and edges than using pixel-wise loss functions, by comparing the semantic distances in the high-dimensional feature space of a pre-trained 2D network (e.g., VGG). However, it is not clear how one should generalize it to 3D medical images, and the attendant implications are still unclear. In this paper, we propose a framework called SOUP-GAN: Super-resolution Optimized Using Perceptual-tuned Generative Adversarial Network (GAN), in order to produce thinner slice (e.g., high resolution in the 'Z' plane) medical images with anti-aliasing and deblurring. The proposed method outperforms other conventional resolution-enhancement methods and previous SR work on medical images upon both qualitative and quantitative comparisons. Specifically, we examine the model in terms of its generalization for various SR ratios and imaging modalities. By addressing those limitations, our model shows promise as a novel 3D SR interpolation technique, providing potential applications in both clinical and research settings.
△ Less
Submitted 4 June, 2021;
originally announced June 2021.
-
A Non-stiff Summation-By-Parts Finite Difference Method for the Scalar Wave Equation in Second Order Form: Characteristic Boundary Conditions and Nonlinear Interfaces
Authors:
Brittany A Erickson,
Jeremy E Kozdon,
Tobias W Harvey
Abstract:
Curvilinear, multiblock summation-by-parts finite difference operators with the simultaneous approximation term method provide a stable and accurate framework for solving the wave equation in second order form. That said, the standard method can become arbitrarily stiff when characteristic boundary conditions and nonlinear interface conditions are used. Here we propose a new technique that avoids…
▽ More
Curvilinear, multiblock summation-by-parts finite difference operators with the simultaneous approximation term method provide a stable and accurate framework for solving the wave equation in second order form. That said, the standard method can become arbitrarily stiff when characteristic boundary conditions and nonlinear interface conditions are used. Here we propose a new technique that avoids this stiffness by using characteristic variables to "upwind" the boundary and interface treatment. This is done through the introduction of an additional block boundary displacement variable. Using a unified energy, which expresses both the standard as well as characteristic boundary and interface treatment, we show that the resulting scheme has semidiscrete energy stability for the scalar anisotropic wave equation. The theoretical stability results are confirmed with numerical experiments that also demonstrate the accuracy and robustness of the proposed scheme. The numerical results also show that the characteristic scheme has a time step restriction based on standard wave propagation considerations and not the boundary closure.
△ Less
Submitted 28 June, 2022; v1 submitted 1 June, 2021;
originally announced June 2021.
-
Hybridized Summation-By-Parts Finite Difference Methods
Authors:
Jeremy E. Kozdon,
Brittany A. Erickson,
Lucas C. Wilcox
Abstract:
We present a hybridization technique for summation-by-parts finite difference methods with weak enforcement of interface and boundary conditions for second order, linear elliptic partial differential equations. The method is based on techniques from the hybridized discontinuous Galerkin literature where local and global problems are defined for the volume and trace grid points, respectively. By us…
▽ More
We present a hybridization technique for summation-by-parts finite difference methods with weak enforcement of interface and boundary conditions for second order, linear elliptic partial differential equations. The method is based on techniques from the hybridized discontinuous Galerkin literature where local and global problems are defined for the volume and trace grid points, respectively. By using a Schur complement technique the volume points can be eliminated, which drastically reduces the system size. We derive both the local and global problems, and show that the linear systems that must be solved are symmetric positive definite. The theoretical stability results are confirmed with numerical experiments as is the accuracy of the method.
△ Less
Submitted 1 June, 2021; v1 submitted 31 January, 2020;
originally announced February 2020.
-
A New Cryogenic Apparatus to Search for the Neutron Electric Dipole Moment
Authors:
M. W. Ahmed,
R. Alarcon,
A. Aleksandrova,
S. Baessler,
L. Barron-Palos,
L. M. Bartoszek,
D. H. Beck,
M. Behzadipour,
I. Berkutov,
J. Bessuille,
M. Blatnik,
M. Broering,
L. J. Broussard,
M. Busch,
R. Carr,
V. Cianciolo,
S. M. Clayton,
M. D. Cooper,
C. Crawford,
S. A. Currie,
C. Daurer,
R. Dipert,
K. Dow,
D. Dutta,
Y. Efremenko
, et al. (69 additional authors not shown)
Abstract:
A cryogenic apparatus is described that enables a new experiment, nEDM@SNS, with a major improvement in sensitivity compared to the existing limit in the search for a neutron Electric Dipole Moment (EDM). It uses superfluid $^4$He to produce a high density of Ultra-Cold Neutrons (UCN) which are contained in a suitably coated pair of measurement cells. The experiment, to be operated at the Spallati…
▽ More
A cryogenic apparatus is described that enables a new experiment, nEDM@SNS, with a major improvement in sensitivity compared to the existing limit in the search for a neutron Electric Dipole Moment (EDM). It uses superfluid $^4$He to produce a high density of Ultra-Cold Neutrons (UCN) which are contained in a suitably coated pair of measurement cells. The experiment, to be operated at the Spallation Neutron Source at Oak Ridge National Laboratory, uses polarized $^3$He from an Atomic Beam Source injected into the superfluid $^4$He and transported to the measurement cells as a co-magnetometer. The superfluid $^4$He is also used as an insulating medium allowing significantly higher electric fields, compared to previous experiments, to be maintained across the measurement cells. These features provide an ultimate statistical uncertainty for the EDM of $2-3\times 10^{-28}$ e-cm, with anticipated systematic uncertainties below this level.
△ Less
Submitted 20 November, 2019; v1 submitted 26 August, 2019;
originally announced August 2019.
-
Interactive segmentation of medical images through fully convolutional neural networks
Authors:
Tomas Sakinis,
Fausto Milletari,
Holger Roth,
Panagiotis Korfiatis,
Petro Kostandy,
Kenneth Philbrick,
Zeynettin Akkus,
Ziyue Xu,
Daguang Xu,
Bradley J. Erickson
Abstract:
Image segmentation plays an essential role in medicine for both diagnostic and interventional tasks. Segmentation approaches are either manual, semi-automated or fully-automated. Manual segmentation offers full control over the quality of the results, but is tedious, time consuming and prone to operator bias. Fully automated methods require no human effort, but often deliver sub-optimal results wi…
▽ More
Image segmentation plays an essential role in medicine for both diagnostic and interventional tasks. Segmentation approaches are either manual, semi-automated or fully-automated. Manual segmentation offers full control over the quality of the results, but is tedious, time consuming and prone to operator bias. Fully automated methods require no human effort, but often deliver sub-optimal results without providing users with the means to make corrections. Semi-automated approaches keep users in control of the results by providing means for interaction, but the main challenge is to offer a good trade-off between precision and required interaction. In this paper we present a deep learning (DL) based semi-automated segmentation approach that aims to be a "smart" interactive tool for region of interest delineation in medical images. We demonstrate its use for segmenting multiple organs on computed tomography (CT) of the abdomen. Our approach solves some of the most pressing clinical challenges: (i) it requires only one to a few user clicks to deliver excellent 2D segmentations in a fast and reliable fashion; (ii) it can generalize to previously unseen structures and "corner cases"; (iii) it delivers results that can be corrected quickly in a smart and intuitive way up to an arbitrary degree of precision chosen by the user and (iv) ensures high accuracy. We present our approach and compare it to other techniques and previous work to show the advantages brought by our method.
△ Less
Submitted 19 March, 2019;
originally announced March 2019.
-
The neutron electric dipole moment experiment at the Spallation Neutron Source
Authors:
K. K. H. Leung,
M. Ahmed,
R. Alarcon,
A. Aleksandrova,
S. Baeßler,
L. Barrón-Palos,
L. Bartoszek,
D. H. Beck,
M. Behzadipour,
J. Bessuille,
M. A. Blatnik,
M. Broering,
L. J. Broussard,
M. Busch,
R. Carr,
P. -H. Chu,
V. Cianciolo,
S. M. Clayton,
M. D. Cooper,
C. Crawford,
S. A. Currie,
C. Daurer,
R. Dipert,
K. Dow,
D. Dutta
, et al. (68 additional authors not shown)
Abstract:
Novel experimental techniques are required to make the next big leap in neutron electric dipole moment experimental sensitivity, both in terms of statistics and systematic error control. The nEDM experiment at the Spallation Neutron Source (nEDM@SNS) will implement the scheme of Golub & Lamoreaux [Phys. Rep., 237, 1 (1994)]. The unique properties of combining polarized ultracold neutrons, polarize…
▽ More
Novel experimental techniques are required to make the next big leap in neutron electric dipole moment experimental sensitivity, both in terms of statistics and systematic error control. The nEDM experiment at the Spallation Neutron Source (nEDM@SNS) will implement the scheme of Golub & Lamoreaux [Phys. Rep., 237, 1 (1994)]. The unique properties of combining polarized ultracold neutrons, polarized $^3$He, and superfluid $^4$He will be exploited to provide a sensitivity to $\sim 10^{-28}\,e{\rm \,\cdot\, cm}$. Our cryogenic apparatus will deploy two small ($3\,{\rm L}$) measurement cells with a high density of ultracold neutrons produced and spin analyzed in situ. The electric field strength, precession time, magnetic shielding, and detected UCN number will all be enhanced compared to previous room temperature Ramsey measurements. Our $^3$He co-magnetometer offers unique control of systematic effects, in particular the Bloch-Siegert induced false EDM. Furthermore, there will be two distinct measurement modes: free precession and dressed spin. This will provide an important self-check of our results. Following five years of "critical component demonstration," our collaboration transitioned to a "large scale integration" phase in 2018. An overview of our measurement techniques, experimental design, and brief updates are described in these proceedings.
△ Less
Submitted 4 October, 2019; v1 submitted 6 March, 2019;
originally announced March 2019.
-
The Buzzard Flock: Dark Energy Survey Synthetic Sky Catalogs
Authors:
Joseph DeRose,
Risa H. Wechsler,
Matthew R. Becker,
Michael T. Busha,
Eli S. Rykoff,
Niall MacCrann,
Brandon Erickson,
August E. Evrard,
Andrey Kravtsov,
Daniel Gruen,
Sahar Allam,
Santiago Avila,
Sarah Bridle,
David Brooks,
Elizabeth Buckley-Geer,
Aurelio Carnero Rosell,
Matias Carrasco Kind,
Jorge Carretero,
Francisco J. Castander,
Ross Cawthon,
Martin Crocce,
Luiz N. da Costa,
Christopher Davis,
Juan De Vicente,
Jörg P. Dietrich
, et al. (30 additional authors not shown)
Abstract:
We present a suite of 18 synthetic sky catalogs designed to support science analysis of galaxies in the Dark Energy Survey Year 1 (DES Y1) data. For each catalog, we use a computationally efficient empirical approach, ADDGALS, to embed galaxies within light-cone outputs of three dark matter simulations that resolve halos with masses above ~5x10^12 h^-1 m_sun at z <= 0.32 and 10^13 h^-1 m_sun at z~…
▽ More
We present a suite of 18 synthetic sky catalogs designed to support science analysis of galaxies in the Dark Energy Survey Year 1 (DES Y1) data. For each catalog, we use a computationally efficient empirical approach, ADDGALS, to embed galaxies within light-cone outputs of three dark matter simulations that resolve halos with masses above ~5x10^12 h^-1 m_sun at z <= 0.32 and 10^13 h^-1 m_sun at z~2. The embedding method is tuned to match the observed evolution of galaxy counts at different luminosities as well as the spatial clustering of the galaxy population. Galaxies are lensed by matter along the line of sight --- including magnification, shear, and multiple images --- using CALCLENS, an algorithm that calculates shear with 0.42 arcmin resolution at galaxy positions in the full catalog. The catalogs presented here, each with the same LCDM cosmology (denoted Buzzard), contain on average 820 million galaxies over an area of 1120 square degrees with positions, magnitudes, shapes, photometric errors, and photometric redshift estimates. We show that the weak-lensing shear catalog, redMaGiC galaxy catalogs and redMaPPer cluster catalogs provide plausible realizations of the same catalogs in the DES Y1 data by comparing their magnitude, color and redshift distributions, angular clustering, and mass-observable relations, making them useful for testing analyses that use these samples. We make public the galaxy samples appropriate for the DES Y1 data, as well as the data vectors used for cosmology analyses on these simulations.
△ Less
Submitted 8 January, 2019;
originally announced January 2019.
-
Comparison of Gaussian process modeling software
Authors:
Collin B. Erickson,
Bruce E. Ankenman,
Susan M. Sanchez
Abstract:
Gaussian process fitting, or kriging, is often used to create a model from a set of data. Many available software packages do this, but we show that very different results can be obtained from different packages even when using the same data and model. We describe the parameterization, features, and optimization used by eight different fitting packages that run on four different platforms. We then…
▽ More
Gaussian process fitting, or kriging, is often used to create a model from a set of data. Many available software packages do this, but we show that very different results can be obtained from different packages even when using the same data and model. We describe the parameterization, features, and optimization used by eight different fitting packages that run on four different platforms. We then compare these eight packages using various data functions and data sets, revealing that there are stark differences between the packages. In addition to comparing the prediction accuracy, the predictive variance--which is important for evaluating precision of predictions and is often used in stopping criteria--is also evaluated.
△ Less
Submitted 9 October, 2017;
originally announced October 2017.
-
Predicting 1p19q Chromosomal Deletion of Low-Grade Gliomas from MR Images using Deep Learning
Authors:
Zeynettin Akkus,
Issa Ali,
Jiri Sedlar,
Timothy L. Kline,
Jay P. Agrawal,
Ian F. Parney,
Caterina Giannini,
Bradley J. Erickson
Abstract:
Objective: Several studies have associated codeletion of chromosome arms 1p/19q in low-grade gliomas (LGG) with positive response to treatment and longer progression free survival. Therefore, predicting 1p/19q status is crucial for effective treatment planning of LGG. In this study, we predict the 1p/19q status from MR images using convolutional neural networks (CNN), which could be a noninvasive…
▽ More
Objective: Several studies have associated codeletion of chromosome arms 1p/19q in low-grade gliomas (LGG) with positive response to treatment and longer progression free survival. Therefore, predicting 1p/19q status is crucial for effective treatment planning of LGG. In this study, we predict the 1p/19q status from MR images using convolutional neural networks (CNN), which could be a noninvasive alternative to surgical biopsy and histopathological analysis. Method: Our method consists of three main steps: image registration, tumor segmentation, and classification of 1p/19q status using CNN. We included a total of 159 LGG with 3 image slices each who had biopsy-proven 1p/19q status (57 nondeleted and 102 codeleted) and preoperative postcontrast-T1 (T1C) and T2 images. We divided our data into training, validation, and test sets. The training data was balanced for equal class probability and then augmented with iterations of random translational shift, rotation, and horizontal and vertical flips to increase the size of the training set. We shuffled and augmented the training data to counter overfitting in each epoch. Finally, we evaluated several configurations of a multi-scale CNN architecture until training and validation accuracies became consistent. Results: The results of the best performing configuration on the unseen test set were 93.3% (sensitivity), 82.22% (specificity), and 87.7% (accuracy). Conclusion: Multi-scale CNN with their self-learning capability provides promising results for predicting 1p/19q status noninvasively based on T1C and T2 images. Significance: Predicting 1p/19q status noninvasively from MR images would allow selecting effective treatment strategies for LGG patients without the need for surgical biopsy.
△ Less
Submitted 21 November, 2016;
originally announced November 2016.
-
Cosmic Shear Measurements with DES Science Verification Data
Authors:
M. R. Becker,
M. A. Troxel,
N. MacCrann,
E. Krause,
T. F. Eifler,
O. Friedrich,
A. Nicola,
A. Refregier,
A. Amara,
D. Bacon,
G. M. Bernstein,
C. Bonnett,
S. L. Bridle,
M. T. Busha,
C. Chang,
S. Dodelson,
B. Erickson,
A. E. Evrard,
J. Frieman,
E. Gaztanaga,
D. Gruen,
W. Hartley,
B. Jain,
M. Jarvis,
T. Kacprzak
, et al. (80 additional authors not shown)
Abstract:
We present measurements of weak gravitational lensing cosmic shear two-point statistics using Dark Energy Survey Science Verification data. We demonstrate that our results are robust to the choice of shear measurement pipeline, either ngmix or im3shape, and robust to the choice of two-point statistic, including both real and Fourier-space statistics. Our results pass a suite of null tests includin…
▽ More
We present measurements of weak gravitational lensing cosmic shear two-point statistics using Dark Energy Survey Science Verification data. We demonstrate that our results are robust to the choice of shear measurement pipeline, either ngmix or im3shape, and robust to the choice of two-point statistic, including both real and Fourier-space statistics. Our results pass a suite of null tests including tests for B-mode contamination and direct tests for any dependence of the two-point functions on a set of 16 observing conditions and galaxy properties, such as seeing, airmass, galaxy color, galaxy magnitude, etc. We furthermore use a large suite of simulations to compute the covariance matrix of the cosmic shear measurements and assign statistical significance to our null tests. We find that our covariance matrix is consistent with the halo model prediction, indicating that it has the appropriate level of halo sample variance. We compare the same jackknife procedure applied to the data and the simulations in order to search for additional sources of noise not captured by the simulations. We find no statistically significant extra sources of noise in the data. The overall detection significance with tomography for our highest source density catalog is 9.7sigma. Cosmological constraints from the measurements in this work are presented in a companion paper (DES et al. 2015).
△ Less
Submitted 27 July, 2016; v1 submitted 20 July, 2015;
originally announced July 2015.
-
Cosmology from Cosmic Shear with DES Science Verification Data
Authors:
The Dark Energy Survey Collaboration,
T. Abbott,
F. B. Abdalla,
S. Allam,
A. Amara,
J. Annis,
R. Armstrong,
D. Bacon,
M. Banerji,
A. H. Bauer,
E. Baxter,
M. R. Becker,
A. Benoit-Lévy,
R. A. Bernstein,
G. M. Bernstein,
E. Bertin,
J. Blazek,
C. Bonnett,
S. L. Bridle,
D. Brooks,
C. Bruderer,
E. Buckley-Geer,
D. L. Burke,
M. T. Busha,
D. Capozzi
, et al. (104 additional authors not shown)
Abstract:
We present the first constraints on cosmology from the Dark Energy Survey (DES), using weak lensing measurements from the preliminary Science Verification (SV) data. We use 139 square degrees of SV data, which is less than 3\% of the full DES survey area. Using cosmic shear 2-point measurements over three redshift bins we find $σ_8 (Ω_{\rm m}/0.3)^{0.5} = 0.81 \pm 0.06$ (68\% confidence), after ma…
▽ More
We present the first constraints on cosmology from the Dark Energy Survey (DES), using weak lensing measurements from the preliminary Science Verification (SV) data. We use 139 square degrees of SV data, which is less than 3\% of the full DES survey area. Using cosmic shear 2-point measurements over three redshift bins we find $σ_8 (Ω_{\rm m}/0.3)^{0.5} = 0.81 \pm 0.06$ (68\% confidence), after marginalising over 7 systematics parameters and 3 other cosmological parameters. We examine the robustness of our results to the choice of data vector and systematics assumed, and find them to be stable. About $20$\% of our error bar comes from marginalising over shear and photometric redshift calibration uncertainties. The current state-of-the-art cosmic shear measurements from CFHTLenS are mildly discrepant with the cosmological constraints from Planck CMB data; our results are consistent with both datasets. Our uncertainties are $\sim$30\% larger than those from CFHTLenS when we carry out a comparable analysis of the two datasets, which we attribute largely to the lower number density of our shear catalogue. We investigate constraints on dark energy and find that, with this small fraction of the full survey, the DES SV constraints make negligible impact on the Planck constraints. The moderate disagreement between the CFHTLenS and Planck values of $σ_8 (Ω_{\rm m}/0.3)^{0.5}$ is present regardless of the value of $w$.
△ Less
Submitted 3 May, 2017; v1 submitted 20 July, 2015;
originally announced July 2015.
-
A High Throughput Workflow Environment for Cosmological Simulations
Authors:
Brandon M. S. Erickson,
Raminderjeet Singh,
August E. Evrard,
Matthew R. Becker,
Michael T. Busha,
Andrey V. Kravtsov,
Suresh Marru,
Marlon Pierce,
Risa H. Wechsler
Abstract:
The next generation of wide-area sky surveys offer the power to place extremely precise constraints on cosmological parameters and to test the source of cosmic acceleration. These observational programs will employ multiple techniques based on a variety of statistical signatures of galaxies and large-scale structure. These techniques have sources of systematic error that need to be understood at t…
▽ More
The next generation of wide-area sky surveys offer the power to place extremely precise constraints on cosmological parameters and to test the source of cosmic acceleration. These observational programs will employ multiple techniques based on a variety of statistical signatures of galaxies and large-scale structure. These techniques have sources of systematic error that need to be understood at the percent-level in order to fully leverage the power of next-generation catalogs. Simulations of large-scale structure provide the means to characterize these uncertainties. We are using XSEDE resources to produce multiple synthetic sky surveys of galaxies and large-scale structure in support of science analysis for the Dark Energy Survey. In order to scale up our production to the level of fifty 10^10-particle simulations, we are working to embed production control within the Apache Airavata workflow environment. We explain our methods and report how the workflow has reduced production time by 40% compared to manual management.
△ Less
Submitted 13 November, 2012; v1 submitted 11 October, 2012;
originally announced October 2012.
-
Influence of Projection in Cluster Cosmology Studies
Authors:
Brandon M. S. Erickson,
Carlos E. Cunha,
August E. Evrard
Abstract:
Projection tends to skew the mass-observable relation of galaxy clusters by creating a small fraction of severely blended systems, those for which the measured observable property of a cluster is strongly boosted relative to the value of its primary host halo. We examine the bias in cosmological parameter estimates caused by incorrectly assuming a Gaussian (projection-free) mass-observable relatio…
▽ More
Projection tends to skew the mass-observable relation of galaxy clusters by creating a small fraction of severely blended systems, those for which the measured observable property of a cluster is strongly boosted relative to the value of its primary host halo. We examine the bias in cosmological parameter estimates caused by incorrectly assuming a Gaussian (projection-free) mass-observable relation when the true relation is non-Gaussian due to projection. We introduce a mixture model for projection and explore Fisher forecasts for a survey of 5000 sq. deg. to z=1.1 and an equivalent mass threshold of 10^13.7 h^-1 solar masses. Using a blended fraction motivated by optical cluster finding applied to the Millennium Simulation and applying Planck and otherwise weak priors, we find that the biases in Omega_DE and w are significant, being factors of 2.8 and 2.4, respectively, times previous forecast uncertainties. Incorporating eight new degrees of freedom to describe cluster selection with projection increases the forecast uncertainty in Omega_DE and w by similar factors. Knowledge of these additional parameters at the 5% level limits degradation in dark energy constraints to <10% relative to projection-free forecasts. We discuss strategies for using simulations and complementary observations to characterize the fraction of blended clusters and their mass selection properties.
△ Less
Submitted 24 June, 2011; v1 submitted 15 June, 2011;
originally announced June 2011.