Results for “SLICING” — Page 23

DOE OSTI2015

TU-EF-BRA-03: Free Induction Decay (without the Decay) and Spin-Echo Imaging

NMR, and Proton Density MRI of the 1D Patient - Anthony Wolbarst Net Voxel Magnetization, m(x,t). T1-MRI; The MRI Device - Lisa Lemen ‘Classical’ NMR; FID Imaging in 1D via k-Space - Nathan Yanasak Spin-Echo; S-E/Spin Warp in a 2D Slice - Ronald Price Magnetic resonance imaging not only reveals the structural, anatomic details of the body, as does CT, but also it can provide information on the physiological status and pathologies of its tissues, like nuclear medicine. It can display high-quality slice and 3D images of organs and vessels viewed from any perspective, with resolution better than 1 mm. MRI is perhaps most extraordinary and notable for the plethora of ways in which it can create unique forms of image contrast, reflective of fundamentally different biophysical phenomena. As with ultrasound, there is no risk from ionizing radiation to the patient or staff, since no X-rays or radioactive nuclei are involved. Instead, MRI harnesses magnetic fields and radio waves to probe the stable nuclei of the ordinary hydrogen atoms (isolated protons) occurring in water and lipid molecules within and around cells. MRI consists, in essence, of creating spatial maps of the electromagnetic environments around these hydrogen nuclei. Spatial variations in the proton milieus can be related to clinical differences in the biochemical and physiological properties and conditions of the associated tissues. Imaging of proton density (PD), and of the tissue proton spin relaxation times known as T1 and T2, all can reveal important clinical information, but they do so with approaches so dissimilar from one another that each is chosen for only certain clinical situations. T1 and T2 in a voxel are determined by different aspects of the rotations and other motions of the water and lipid molecules involved, as constrained by the local biophysical surroundings within and between its cells – and they, in turn, depend on the type of tissue and its state of health. Three other common applications of MRI exploit its capability to detect and image distinct movements of fluids: MR angiography (MRA), which rivals CT angiography but often requires no contrast medium, monitors the bulk flow of blood; functional MRI ( f MRI), distinguishes the perfusion of oxygenated blood from that of de-oxygenated, and lights up parts of the brain that are activated by a stimulus, rather like PET; and diffusion tensor imaging (DTI) indicates the diffusion of free water along tracts of axons, thereby bringing nerve trunks into view. There are variants on all of these themes, and on others as well. Magnetic Resonance Spectroscopy (MRS), for example, can perform non-invasive ‘virtual biopsies’ that allow identification of certain cancers and other lesions. And an MRI-guided needle biopsy can sample brain tissue from a region only millimeters in dimensions. MRI, however, involves deeper and more complex aspects of physics, technology, and biology than do most other imaging modalities, and it is widely considered to be correspondingly more difficult to learn. We could probably cover all of this rather comprehensively if we had 50 hours available rather than 2 ̶ but, to paraphrase a former Secretary of Defense, you tell your story in the time you have allotted. The four presenters and another physicist, Kevin King from GE, have combined their efforts to co-author a single slide show that describes essentials of MRI as simply as possible. It is obviously far from thorough, but hopefully it will succeed in explaining some of the basics in a simplified but still valid fashion; in providing a taste of the numerous capabilities and complexities of the modality; and in whetting your appetite to learn more. Part I. NMR, and Proton Density MRI of the 1D Patient (Wolbarst), begins with an introductory case study that illustrates a half dozen ways in which MRI provides valuable clinical information. It then explores the nuclear magnetic resonance (NMR) phenomenon, which underlies MRI. NMR can be introduced with either of two approaches. In the first, one thinks (loosely) of the nuclei of hydrogen atoms as (rotating and charged and therefore) magnetic objects, whose spin-axes tend to align in a strong external magnetic field, much like a compass needle. As with the Bohr atom, this spin-up/spin-down picture is a highly abridged version of the full quantum mechanical treatment, but still it leads to some useful, legitimate pictures of the NMR process occurring within a voxel: When RF photons of the correct (Larmor) frequency elevate protons in a fixed magnetic field out of their lower-energy spin state into the upper, the NMR phenomenon is indicated by the detectable absorption of RF power. With the addition of a linear gradient field along a multi-voxel, one-dimensional patient/phantom, as well, we can determine the water content of each compartment – an example of a real MRI study, albeit in 1D. Part I concludes with a discussion of the net magnetization at position x, m0(x), under conditions of dynamic thermal equilibrium, which leads into: Part II. Net Voxel Magnetization, m(x,t); T1-MRI; The MRI Device (Lemen), investigates the biophysics of the form of proton spin relaxation process characterized by the time T1. It then moves on to the creation of an MR image that displays the spatial variation in the values of this clinically relevant parameter, again in 1D. Finally, the design and workings of a clinical MRI machine are sketched, in preparation for: Part III. ‘Classical’ NMR; FID Imaging in 1D via k-Space (Yanasak) presents the second standard approach to NMR and MRI, the classical model. It focuses on the time dependence of the net nuclear magnetization, m(x,t), the overall magnetic field generated by the cohort of protons in the voxel at position x. Quite remarkably, this nuclear net magnetization itself acts in a strong magnetic field like a gyroscope in a gravitational field. This tack is better for explaining Free Induction Decay (FID), which involves a brief introduction to the Fourier transform and k-space. This leads to conventional Spin-Echo (S-E) reconstruction techniques for creating clinical images from raw data, and sets the stage for: Part IV. Spin-Echo; S-E / Spin Warp in a 2D Slice (Price) discusses application of the S- E sequence of radiofrequency pulses and gradient magnetic fields to the 1D patient. T2 is introduced but not explained. This Part also considers how to manipulate the image acquisition parameters so as to generate clinical pictures with contrast dominated by spatial variations in PD, T1, or T2. We conclude by demonstrating the spin-warp approach to imaging in 2D with a simple 2×2, 4-voxel example. Much of this material is presented in more detail in the chapter “MRI of the One-dimensional Patient, Part I”, in Advances in Medical Physics, Vol 5 (2014). Copies are on display at the Medical Physics Publishing booth. Learning Objectives: The participant will learn about the processes of NMR and T1 spin relaxation in a tissue voxel in a uniform magnetic field. The participant will learn about combining spin-up/spin-down NMR and T1 processes with a linear gradient to effect frequency-encoding of voxel spatial position. This approach can be used to create proton density and T1 MRI maps, respectively, of the contents of multi-voxel 1D phantoms. The participant will learn about how the ‘classical’ model of NMR it can generate Free Induction Decay (FID) images of 1D phantoms, which involves the use of the Fourier transform in k-space. This can lead simply into standard Spin-Echo images. The participant will learn about extending Spin-Echo imaging into 2 and more dimensions.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

TU-EF-BRA-04: Into 2, 3, and 4 Dimensions

NMR, and Proton Density MRI of the 1D Patient - Anthony Wolbarst Net Voxel Magnetization, m(x,t). T1-MRI; The MRI Device - Lisa Lemen ‘Classical’ NMR; FID Imaging in 1D via k-Space - Nathan Yanasak Spin-Echo; S-E/Spin Warp in a 2D Slice - Ronald Price Magnetic resonance imaging not only reveals the structural, anatomic details of the body, as does CT, but also it can provide information on the physiological status and pathologies of its tissues, like nuclear medicine. It can display high-quality slice and 3D images of organs and vessels viewed from any perspective, with resolution better than 1 mm. MRI is perhaps most extraordinary and notable for the plethora of ways in which it can create unique forms of image contrast, reflective of fundamentally different biophysical phenomena. As with ultrasound, there is no risk from ionizing radiation to the patient or staff, since no X-rays or radioactive nuclei are involved. Instead, MRI harnesses magnetic fields and radio waves to probe the stable nuclei of the ordinary hydrogen atoms (isolated protons) occurring in water and lipid molecules within and around cells. MRI consists, in essence, of creating spatial maps of the electromagnetic environments around these hydrogen nuclei. Spatial variations in the proton milieus can be related to clinical differences in the biochemical and physiological properties and conditions of the associated tissues. Imaging of proton density (PD), and of the tissue proton spin relaxation times known as T1 and T2, all can reveal important clinical information, but they do so with approaches so dissimilar from one another that each is chosen for only certain clinical situations. T1 and T2 in a voxel are determined by different aspects of the rotations and other motions of the water and lipid molecules involved, as constrained by the local biophysical surroundings within and between its cells – and they, in turn, depend on the type of tissue and its state of health. Three other common applications of MRI exploit its capability to detect and image distinct movements of fluids: MR angiography (MRA), which rivals CT angiography but often requires no contrast medium, monitors the bulk flow of blood; functional MRI ( f MRI), distinguishes the perfusion of oxygenated blood from that of de-oxygenated, and lights up parts of the brain that are activated by a stimulus, rather like PET; and diffusion tensor imaging (DTI) indicates the diffusion of free water along tracts of axons, thereby bringing nerve trunks into view. There are variants on all of these themes, and on others as well. Magnetic Resonance Spectroscopy (MRS), for example, can perform non-invasive ‘virtual biopsies’ that allow identification of certain cancers and other lesions. And an MRI-guided needle biopsy can sample brain tissue from a region only millimeters in dimensions. MRI, however, involves deeper and more complex aspects of physics, technology, and biology than do most other imaging modalities, and it is widely considered to be correspondingly more difficult to learn. We could probably cover all of this rather comprehensively if we had 50 hours available rather than 2 ̶ but, to paraphrase a former Secretary of Defense, you tell your story in the time you have allotted. The four presenters and another physicist, Kevin King from GE, have combined their efforts to co-author a single slide show that describes essentials of MRI as simply as possible. It is obviously far from thorough, but hopefully it will succeed in explaining some of the basics in a simplified but still valid fashion; in providing a taste of the numerous capabilities and complexities of the modality; and in whetting your appetite to learn more. Part I. NMR, and Proton Density MRI of the 1D Patient (Wolbarst), begins with an introductory case study that illustrates a half dozen ways in which MRI provides valuable clinical information. It then explores the nuclear magnetic resonance (NMR) phenomenon, which underlies MRI. NMR can be introduced with either of two approaches. In the first, one thinks (loosely) of the nuclei of hydrogen atoms as (rotating and charged and therefore) magnetic objects, whose spin-axes tend to align in a strong external magnetic field, much like a compass needle. As with the Bohr atom, this spin-up/spin-down picture is a highly abridged version of the full quantum mechanical treatment, but still it leads to some useful, legitimate pictures of the NMR process occurring within a voxel: When RF photons of the correct (Larmor) frequency elevate protons in a fixed magnetic field out of their lower-energy spin state into the upper, the NMR phenomenon is indicated by the detectable absorption of RF power. With the addition of a linear gradient field along a multi-voxel, one-dimensional patient/phantom, as well, we can determine the water content of each compartment – an example of a real MRI study, albeit in 1D. Part I concludes with a discussion of the net magnetization at position x, m0(x), under conditions of dynamic thermal equilibrium, which leads into: Part II. Net Voxel Magnetization, m(x,t); T1-MRI; The MRI Device (Lemen), investigates the biophysics of the form of proton spin relaxation process characterized by the time T1. It then moves on to the creation of an MR image that displays the spatial variation in the values of this clinically relevant parameter, again in 1D. Finally, the design and workings of a clinical MRI machine are sketched, in preparation for: Part III. ‘Classical’ NMR; FID Imaging in 1D via k-Space (Yanasak) presents the second standard approach to NMR and MRI, the classical model. It focuses on the time dependence of the net nuclear magnetization, m(x,t), the overall magnetic field generated by the cohort of protons in the voxel at position x. Quite remarkably, this nuclear net magnetization itself acts in a strong magnetic field like a gyroscope in a gravitational field. This tack is better for explaining Free Induction Decay (FID), which involves a brief introduction to the Fourier transform and k-space. This leads to conventional Spin-Echo (S-E) reconstruction techniques for creating clinical images from raw data, and sets the stage for: Part IV. Spin-Echo; S-E / Spin Warp in a 2D Slice (Price) discusses application of the S- E sequence of radiofrequency pulses and gradient magnetic fields to the 1D patient. T2 is introduced but not explained. This Part also considers how to manipulate the image acquisition parameters so as to generate clinical pictures with contrast dominated by spatial variations in PD, T1, or T2. We conclude by demonstrating the spin-warp approach to imaging in 2D with a simple 2×2, 4-voxel example. Much of this material is presented in more detail in the chapter “MRI of the One-dimensional Patient, Part I”, in Advances in Medical Physics, Vol 5 (2014). Copies are on display at the Medical Physics Publishing booth. Learning Objectives: The participant will learn about the processes of NMR and T1 spin relaxation in a tissue voxel in a uniform magnetic field. The participant will learn about combining spin-up/spin-down NMR and T1 processes with a linear gradient to effect frequency-encoding of voxel spatial position. This approach can be used to create proton density and T1 MRI maps, respectively, of the contents of multi-voxel 1D phantoms. The participant will learn about how the ‘classical’ model of NMR it can generate Free Induction Decay (FID) images of 1D phantoms, which involves the use of the Fourier transform in k-space. This can lead simply into standard Spin-Echo images. The participant will learn about extending Spin-Echo imaging into 2 and more dimensions.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

WE-G-BRD-01: A Data-Driven 4D-MRI Motion Model to Estimate Full Field-Of-View Abdominal Motion From 2D Image Navigators During MR-Linac Treatment

Purpose: To estimate full field-of-view abdominal respiratory motion from fast 2D image navigators using a 4D-MRI based motion model. This will allow for radiation dose accumulation mapping during MR-Linac treatment. Methods: Experiments were conducted on a Philips Ingenia 1.5T MRI. First, a retrospectively ordered 4D-MRI was constructed using 3D transient-bSSFP with radial in-plane sampling. Motion fields were calculated through 3D non-rigid registration. From these motion fields a PCA-based abdominal motion model was constructed and used to warp a 3D reference volume to fast 2D cine-MR image navigators that can be used for real-time tracking. To test this procedure, a time-series consisting of two interleaved orthogonal slices (sagittal and coronal), positioned on the pancreas or kidneys, were acquired for 1m38s (dynamic scan-time=0.196ms), during normal, shallow, or deep breathing. The coronal slices were used to update the optimal weights for the first two PCA components, in order to warp the 3D reference image and construct a dynamic 4D-MRI time-series. The interleaved sagittal slices served as an independent measure to test the model’s accuracy and fit. Spatial maps of the root-mean-squared error (RMSE) and histograms of the motion differences within the pancreas and kidneys were used to evaluate the method. Results: Cranio-caudal motion was accurately calculated within the pancreas using the model for normal and shallow breathing with an RMSE of 1.6mm and 1.5mm and a histogram median and standard deviation below 0.2 and 1.7mm, respectively. For deep-breathing an underestimation of the inhale amplitude was observed (RMSE=4.1mm). Respiratory-induced antero-posterior and lateral motion were correctly mapped (RMSE=0.6/0.5mm). Kidney motion demonstrated good motion estimation with RMSE-values of 0.95 and 2.4mm for the right and left kidney, respectively. Conclusion: We have demonstrated a method that can calculate dynamic 3D abdominal motion in a large volume, while acquiring real-time cine-MR images for MR-guided radiotherapy.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

Image-based reconstruction of three-dimensional myocardial infarct geometry for patient-specific modeling of cardiac electrophysiology

Purpose: Accurate three-dimensional (3D) reconstruction of myocardial infarct geometry is crucial to patient-specific modeling of the heart aimed at providing therapeutic guidance in ischemic cardiomyopathy. However, myocardial infarct imaging is clinically performed using two-dimensional (2D) late-gadolinium enhanced cardiac magnetic resonance (LGE-CMR) techniques, and a method to build accurate 3D infarct reconstructions from the 2D LGE-CMR images has been lacking. The purpose of this study was to address this need. Methods: The authors developed a novel methodology to reconstruct 3D infarct geometry from segmented low-resolution (Lo-res) clinical LGE-CMR images. Their methodology employed the so-called logarithm of odds (LogOdds) function to implicitly represent the shape of the infarct in segmented image slices as LogOdds maps. These 2D maps were then interpolated into a 3D image, and the result transformed via the inverse of LogOdds to a binary image representing the 3D infarct geometry. To assess the efficacy of this method, the authors utilized 39 high-resolution (Hi-res) LGE-CMR images, including 36 in vivo acquisitions of human subjects with prior myocardial infarction and 3 ex vivo scans of canine hearts following coronary ligation to induce infarction. The infarct was manually segmented by trained experts in each slice of the Hi-res images, and the segmented data were downsampled to typical clinical resolution. The proposed method was then used to reconstruct 3D infarct geometry from the downsampled images, and the resulting reconstructions were compared with the manually segmented data. The method was extensively evaluated using metrics based on geometry as well as results of electrophysiological simulations of cardiac sinus rhythm and ventricular tachycardia in individual hearts. Several alternative reconstruction techniques were also implemented and compared with the proposed method. Results: The accuracy of the LogOdds method in reconstructing 3D infarct geometry, as measured by the Dice similarity coefficient, was 82.10% ± 6.58%, a significantly higher value than those of the alternative reconstruction methods. Among outcomes of electrophysiological simulations with infarct reconstructions generated by various methods, the simulation results corresponding to the LogOdds method showed the smallest deviation from those corresponding to the manual reconstructions, as measured by metrics based on both activation maps and pseudo-ECGs. Conclusions: The authors have developed a novel method for reconstructing 3D infarct geometry from segmented slices of Lo-res clinical 2D LGE-CMR images. This method outperformed alternative approaches in reproducing expert manual 3D reconstructions and in electrophysiological simulations.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

A dual-view digital tomosynthesis imaging technique for improved chest imaging

Purpose: Digital tomosynthesis (DTS) has been shown to be useful for reducing the overlapping of abnormalities with anatomical structures at various depth levels along the posterior–anterior (PA) direction in chest radiography. However, DTS provides crude three-dimensional (3D) images that have poor resolution in the lateral view and can only be displayed with reasonable quality in the PA view. Furthermore, the spillover of high-contrast objects from off-fulcrum planes generates artifacts that may impede the diagnostic use of the DTS images. In this paper, the authors describe and demonstrate the use of a dual-view DTS technique to improve the accuracy of the reconstructed volume image data for more accurate rendition of the anatomy and slice images with improved resolution and reduced artifacts, thus allowing the 3D image data to be viewed in views other than the PA view. Methods: With the dual-view DTS technique, limited angle scans are performed and projection images are acquired in two orthogonal views: PA and lateral. The dual-view projection data are used together to reconstruct 3D images using the maximum likelihood expectation maximization iterative algorithm. In this study, projection images were simulated or experimentally acquired over 360° using the scanning geometry for cone beam computed tomography (CBCT). While all projections were used to reconstruct CBCT images, selected projections were extracted and used to reconstruct single- and dual-view DTS images for comparison with the CBCT images. For realistic demonstration and comparison, a digital chest phantom derived from clinical CT images was used for the simulation study. An anthropomorphic chest phantom was imaged for the experimental study. The resultant dual-view DTS images were visually compared with the single-view DTS images and CBCT images for the presence of image artifacts and accuracy of CT numbers and anatomy and quantitatively compared with root-mean-square-deviation (RMSD) values computed using the digital chest phantom or the CBCT images as the reference in the simulation and experimental study, respectively. High-contrast wires with vertical, oblique, and horizontal orientations in a PA view plane were also imaged to investigate the spatial resolutions and how the wire signals spread in the PA view and lateral view slice images. Results: Both the digital phantom images (simulated) and the anthropomorphic phantom images (experimentally generated) demonstrated that the dual-view DTS technique resulted in improved spatial resolution in the depth (PA) direction, more accurate representation of the anatomy, and significantly reduced artifacts. The RMSD values corroborate well with visual observations with substantially lower RMSD values measured for the dual-view DTS images as compared to those measured for the single-view DTS images. The imaging experiment with the high-contrast wires shows that while the vertical and oblique wires could be resolved in the lateral view in both single- and dual-view DTS images, the horizontal wire could only be resolved in the dual-view DTS images. This indicates that with single-view DTS, the wire signals spread liberally to off-fulcrum planes and generated wire shadow there. Conclusions: The authors have demonstrated both visually and quantitatively that the dual-view DTS technique can be used to achieve more accurate rendition of the anatomy and to obtain slice images with improved resolution and reduced artifacts as compared to the single-view DTS technique, thus allowing the 3D image data to be viewed in views other than the PA view. These advantages could make the dual-view DTS technique useful in situations where better separation of the objects-of-interest from the off-fulcrum structures or more accurate 3D rendition of the anatomy are required while a regular CT examination is undesirable due to radiation dose considerations.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

Automated selection of the optimal cardiac phase for single-beat coronary CT angiography reconstruction

Purpose: Reconstructing a low-motion cardiac phase is expected to improve coronary artery visualization in coronary computed tomography angiography (CCTA) exams. This study developed an automated algorithm for selecting the optimal cardiac phase for CCTA reconstruction. The algorithm uses prospectively gated, single-beat, multiphase data made possible by wide cone-beam imaging. The proposed algorithm differs from previous approaches because the optimal phase is identified based on vessel image quality (IQ) directly, compared to previous approaches that included motion estimation and interphase processing. Because there is no processing of interphase information, the algorithm can be applied to any sampling of image phases, making it suited for prospectively gated studies where only a subset of phases are available. Methods: An automated algorithm was developed to select the optimal phase based on quantitative IQ metrics. For each reconstructed slice at each reconstructed phase, an image quality metric was calculated based on measures of circularity and edge strength of through-plane vessels. The image quality metric was aggregated across slices, while a metric of vessel-location consistency was used to ignore slices that did not contain through-plane vessels. The algorithm performance was evaluated using two observer studies. Fourteen single-beat cardiac CT exams (Revolution CT, GE Healthcare, Chalfont St. Giles, UK) reconstructed at 2% intervals were evaluated for best systolic (1), diastolic (6), or systolic and diastolic phases (7) by three readers and the algorithm. Pairwise inter-reader and reader-algorithm agreement was evaluated using the mean absolute difference (MAD) and concordance correlation coefficient (CCC) between the reader and algorithm-selected phases. A reader-consensus best phase was determined and compared to the algorithm selected phase. In cases where the algorithm and consensus best phases differed by more than 2%, IQ was scored by three readers using a five point Likert scale. Results: There was no statistically significant difference between inter-reader and reader-algorithm agreement for either MAD or CCC metrics (p > 0.1). The algorithm phase was within 2% of the consensus phase in 15/21 of cases. The average absolute difference between consensus and algorithm best phases was 2.29% ± 2.47%, with a maximum difference of 8%. Average image quality scores for the algorithm chosen best phase were 4.01 ± 0.65 overall, 3.33 ± 1.27 for right coronary artery (RCA), 4.50 ± 0.35 for left anterior descending (LAD) artery, and 4.50 ± 0.35 for left circumflex artery (LCX). Average image quality scores for the consensus best phase were 4.11 ± 0.54 overall, 3.44 ± 1.03 for RCA, 4.39 ± 0.39 for LAD, and 4.50 ± 0.18 for LCX. There was no statistically significant difference (p > 0.1) between the image quality scores of the algorithm phase and the consensus phase. Conclusions: The proposed algorithm was statistically equivalent to a reader in selecting an optimal cardiac phase for CCTA exams. When reader and algorithm phases differed by >2%, image quality as rated by blinded readers was statistically equivalent. By detecting the optimal phase for CCTA reconstruction, the proposed algorithm is expected to improve coronary artery visualization in CCTA exams.

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2016

SU-C-207B-04: Automated Segmentation of Pectoral Muscle in MR Images of Dense Breasts

Purpose: To develop and evaluate a fully automated method for segmentation of the pectoral muscle boundary in Magnetic Resonance Imaging (MRI) of dense breasts. Methods: Segmentation of the pectoral muscle is an important part of automatic breast image analysis methods. Current methods for segmenting the pectoral muscle in breast MRI have difficulties delineating the muscle border correctly in breasts with a large proportion of fibroglandular tissue (i.e., dense breasts). Hence, an automated method based on dynamic programming was developed, incorporating heuristics aimed at shape, location and gradient features.To assess the method, the pectoral muscle was segmented in 91 randomly selected participants (mean age 56.6 years, range 49.5–75.2 years) from a large MRI screening trial in women with dense breasts (ACR BI-RADS category 4). Each MR dataset consisted of 178 or 179 T1-weighted images with voxel size 0.64 × 0.64 × 1.00 mm3. All images (n=16,287) were reviewed and scored by a radiologist. In contrast to volume overlap coefficients, such as DICE, the radiologist detected deviations in the segmented muscle border and determined whether the result would impact the ability to accurately determine the volume of fibroglandular tissue and detection of breast lesions. Results: According to the radiologist’s scores, 95.5% of the slices did not mask breast tissue in such way that it could affect detection of breast lesions or volume measurements. In 13.1% of the slices a deviation in the segmented muscle border was present which would not impact breast lesion detection. In 70 datasets (78%) at least 95% of the slices were segmented in such a way it would not affect detection of breast lesions, and in 60 (66%) datasets this was 100%. Conclusion: Dynamic programming with dedicated heuristics shows promising potential to segment the pectoral muscle in women with dense breasts.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-D-BRA-03: Analysis of Systematic Errors with 2D/3D Image Registration for Target Localization and Treatment Delivery in Stereotactic Radiosurgery

Purpose: Determine systematic deviations between 2D/3D and 3D/3D image registrations with six degrees of freedom (6DOF) for various imaging modalities and registration algorithms on the Varian Edge Linac. Methods: The 6DOF systematic errors were assessed by comparing automated 2D/3D (kV/MV vs. CT) with 3D/3D (CBCT vs. CT) image registrations from different imaging pairs, CT slice thicknesses, couch angles, similarity measures, etc., using a Rando head and a pelvic phantom. The 2D/3D image registration accuracy was evaluated at different treatment sites (intra-cranial and extra-cranial) by statistically analyzing 2D/3D pre-treatment verification against 3D/3D localization of 192 Stereotactic Radiosurgery/Stereotactic Body Radiation Therapy treatment fractions for 88 patients. Results: The systematic errors of 2D/3D image registration using kV-kV, MV-kV and MV-MV image pairs using 0.8 mm slice thickness CT images were within 0.3 mm and 0.3° for translations and rotations with a 95% confidence interval (CI). No significant difference between 2D/3D and 3D/3D image registrations (P>0.05) was observed for target localization at various CT slice thicknesses ranging from 0.8 to 3 mm. Couch angles (30, 45, 60 degree) did not impact the accuracy of 2D/3D image registration. Using pattern intensity with content image filtering was recommended for 2D/3D image registration to achieve the best accuracy. For the patient study, translational error was within 2 mm and rotational error was within 0.6 degrees in terms of 95% CI for 2D/3D image registration. For intra-cranial sites, means and std. deviations of translational errors were −0.2±0.7, 0.04±0.5, 0.1±0.4 mm for LNG, LAT, VRT directions, respectively. For extra-cranial sites, means and std. deviations of translational errors were - 0.04±1, 0.2±1, 0.1±1 mm for LNG, LAT, VRT directions, respectively. 2D/3D image registration uncertainties for intra-cranial and extra-cranial sites were comparable. Conclusion: The Varian Edge radiosurgery 6DOF-based system, can perform 2D/3D image registration with high accuracy for target localization in image-guided stereotactic radiosurgery. The work was supported by a Research Scholar Grant, RSG-15-137-01-CCE from the American Cancer Society.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-I-31: Reproducibility of An Automatic Exposure Control Technique in the Low-Dose CT Scan of Cardiac PET/CT Exams

Purpose: Two separate low-dose CT scans are usually performed for attenuation correction of rest and stress N-13 ammonia PET/CT myocardial perfusion imaging (PET/CT). We utilize an automatic exposure control (AEC) technique to reduce CT radiation dose while maintaining perfusion image quality. Our goal is to assess the reproducibility of displayed CT dose index (CTDI) on same-day repeat CT scans (CT1 and CT2). Methods: Retrospectively, we reviewed CT images of PET/CT studies performed on the same day. Low-dose CT utilized AEC technique based on tube current modulation called Smart-mA. The scan parameters were 64 × 0.625mm collimation, 5mm slice thickness, 0.984 pitch, 1-sec rotation time, 120 kVp, and noise index 50 with a range of 10–200 mA. The scan length matched with PET field of view (FOV) with the heart near the middle of axial FOV. We identified the reference slice number (RS) for an anatomical landmark (carina) and used it to estimate axial shift between two CTs. For patient size, we measured an effective diameter on the reference slice. The effect of patient positioning to CTDI was evaluated using the table height. We calculated the absolute percent difference of the CTDI (%diff) for estimation of the reproducibility. Results: The study included 168 adults with an average body-mass index of 31.72 ± 9.10 (kg/m{sup 2}) and effective diameter was 32.72 ± 4.60 cm. The average CTDI was 1.95 ± 1.40 mGy for CT1 and 1.97 ± 1.42mGy for CT2. The mean %diff was 7.8 ± 6.8%. Linear regression analysis showed a significant correlation between the table height and %diff CTDI. (r=0.82, p<0.001) Conclusion: We have shown for the first time in human subjects, using two same-day CT images, that the AEC technique in low-dose CT is reproducible within 10% and significantly depends on the patient centering.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-J-156: The Feasibility of MR-Only IMRT Planning for Prostate Anatomy

Purpose: For prostate anatomy, previous investigations have shown that simulated CT (sCT) generated from MR images can be used for accurate dose computation. In this study, we demonstrate the feasibility of MR-only IMRT planning for prostate case. Methods: Regular CT (rCT) and MR images of the same patient were acquired for prostate anatomy. Regions-of-interest (ROIs) i.e. target and risk structures are delineated on the rCT. A simulated CT (sCT) is generated from the MR image using the method described by Schadewaldt N et al. Their work establishes the clinical acceptability of dose calculation results on the sCT when compared to rCT. rCT and sCT are rigidly registered to ensure proper alignment between the two images. rCT and sCT are overlaid on each other and slice-wise visual inspection confirms excellent agreement between the two images. ROIs on the rCT are copied over to sCT. Philips AutoPlanning solution is used for generating treatment plans. The same treatment technique protocol (plan parameters and clinical goals) is used to generate AutoPlan-rCT and AutoPlan-sCT respectively for rCT and and sCT. DVH comparison on ROIs and slice-wise evaluation of dose is performed between AutoPlan-rCT and AutoPlan-sCT. Delivery parameters i.e. beam and corresponding segments from the AutoPlan-sCT are copied over to rCT and dose is computed to get AutoPlan-sCT-on-rCT. Results: Plan evaluation is done based on Dose Volume Histogram (DVH) of ROIs and manual slice-wise inspection of dose distribution. Both AutoPlan-rCT and AutoPlan-sCT provide a clinically acceptable plan. Also, AutoPlan-sCT-on-rCT shows excellent agreement with AutoPlan-sCT. Conclusion: The study demonstrates that it is feasible to do IMRT planning on the simulated CT image obtained from MR image for prostate anatomy. The research is supported by Philips India Ltd.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-470: Implementation of a Daily Quality Assurance Phantom for Six Degrees-Of-Freedom Couch Testing

Purpose: To develop a daily quality assurance (QA) device that can test the 6DoF (degrees of freedom) couch repositioning accuracy, prior to SBRT treatment deliveries, with an accuracy of ±0.3 degrees and ±0.3 mm. Methods: A daily QA phantom is designed with a focus on the derived center of projections of its markers, rather than tracking its individual markers one at a time. This approach can be the most favorable to address the intended machining accuracy of the QA phantom and the CBCT spatial resolution limitations, primarily 1 mm min slice thickness, simultaneously. With the current design, ±0.1 mm congruence of the resultant center of gravity of the markers with reference CT (0.6 mm minimum slice thickness) vs CBCT (1.0 mm minimum slice thickness) can be achieved. If successful, the QA device should be qualified to test 6DoF couch performance with a gauged accuracy of ±0.3 degrees/±0.3 mm. Testing is performed for the Varian True Beam 2.0 6DoF system. Results: Once the QA phantom is constructed and tested, agreement of the center of gravity of the reference CT scan and the CBCT scan of ±0.1 mm is achieved. This has translated into a consistent 3D-3D match on the treatment machine, CT vs CBCT, with a repetitive ±0.1 mm variation, thus exceeding our expectations. We have deployed the phantom for daily QA on one of our accelerators, and found that the QA time has increased by only 10 minutes. Conclusion: A 6DoF phantom has been designed (patent pending) and built with a realistic work flow in mind where the daily couch accuracy QA checks taking less than 10 minutes. Current developments include integration with the Varian’s Machine Performance Check consistency module.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-G-IeP4-11: Monitoring Tumor Growth in Subcutaneous Murine Tumor Model in Vivo: A Comparison Between MRI and Small Animal CT

Purpose: The purpose of the study is to compare the volume measurement of subcutaneous tumors in mice with different imaging platforms, namely a GE MRI and a Sofie-Biosciences small animal CT scanner. Methods: A549 human lung carcinoma cells and FaDu human head and neck squamous cell carcinoma cells were implanted subcutaneously into flanks of nude mice. Three FaDu tumors and three A549 tumors were included in this study. The MRI scans were done with a GE Signa 1.5 Tesla MR scanner using a fast T2-weighted sequence (70mm FOV and 1.2mm slice thickness), while the CT scans were done with the CT scanner on a Sofie-Biosciences G8 PET/CT platform dedicated for small animal studies (48mm FOV and 0.2mm slice thickness). Imaging contrast agent was not used in this study. Based on the DICOM images from MRI and CT scans, the tumors were contoured with Philips DICOM Viewer and the tumor volumes were obtained by summing up the contoured area and multiplied by the slice thickness. Results: The volume measurements based on the CT scans agree reasonably with that obtained with MR images for the subcutaneous tumors. The mean difference in the absolute tumor volumes between MRI- and CT-based measurements was found to be −6.2% ± 1.0%, with the difference defined as (VMR – VCT)*100%/VMR. Furthermore, we evaluated the normalized tumor volumes, which were defined for each tumor as V/V{sub 0} where V{sub 0} stands for the volume from the first MR or CT scan. The mean difference in the normalized tumor volumes was found to be 0.10% ± 0.96%. Conclusion: Despite the fact that the difference between normal and abnormal tissues is often less clear on small animal CT images than on MR images, one can still obtain reasonable tumor volume information with the small animal CT scans for subcutaneous murine xenograft models.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TU-AB-207A-03: Image Quality, Dose, and Clinical Applications

Practicing medical physicists are often time charged with the tasks of evaluating and troubleshooting complex image quality issues related to CT scanners. This course will equip them with a solid and practical understanding of common CT imaging chain and its major components with emphasis on acquisition physics and hardware, reconstruction, artifacts, image quality, dose, and advanced clinical applications. The core objective is to explain the effects of these major system components on the image quality. This course will not focus on the rapid-changing advanced technologies given the two-hour time limit, but the fundamental principles discussed in this course may facilitate better understanding of those more complicated technologies. The course will begin with an overview of CT acquisition physics and geometry. X-ray tube and CT detector are important acquisition hardware critical to the overall image quality. Each of these two subsystems consists of several major components. An in-depth description of the function and failure modes of these components will be provided. Examples of artifacts related to these failure modes will be presented: off-focal radiation, tube arcing, heel effect, oil bubble, offset drift effect, cross-talk effect, and bad pixels. The fundamentals of CT image reconstruction will first be discussed on an intuitive level. Approaches that do not require rigorous derivation of mathematical formulations will be presented. This is followed by a detailed derivation of the Fourier slice theorem: the foundation of the FBP algorithm. FBP for parallel-beam, fan-beam, and cone-beam geometries will be discussed. To address the issue of radiation dose related to x-ray CT, recent advances in iterative reconstruction, their advantages, and clinical applications will also be described. Because of the nature of fundamental physics and mathematics, limitations in data acquisition, and non-ideal conditions of major system components, image artifact often arise in the reconstructed images. Because of the limited scope of this course, only major imaging artifacts, their appearance, and possible mitigation and corrections will be discussed. Assessment of the performance of a CT scanner is a complicated subject. Procedures to measure common image quality metrics such as high contrast spatial resolution, low contrast detectability, and slice profile will be described. The reason why these metrics used for FBP may not be sufficient for statistical iterative reconstruction will be explained. Optimizing radiation dose requires comprehension of CT dose metrics. This course will briefly describe various dose metrics, and interaction with acquisition parameters and patient habitus. CT is among the most frequently used imaging tools due to its superior image quality, easy to operate, and a broad range of applications. This course will present several interesting CT applications such as a mobile CT unit on an ambulance for stroke patients, low dose lung cancer screening, and single heartbeat cardiac CT. Learning Objectives: Understand the function and impact of major components of X-ray tube on the image quality. Understand the function and impact of major components of CT detector on the image quality. Be familiar with the basic procedure of CT image reconstruction. Understand the effect of image reconstruction on CT image quality and artifacts. Understand the root causes of common CT image artifacts. Be familiar with image quality metrics especially high and low contrast resolution, noise power spectrum, slice sensitivity profile, etc. Understand why basic image quality metrics used for FBP may not be sufficient to characterize the performance of advanced iterative reconstruction. Be familiar with various CT dose metrics and their interaction with acquisition parameters. New development in advanced CT clinical applications. JH: Employee of GE Healthcare. FD: No disclosure.; J. Hsieh, Jiang Hsieh is an employee of GE Healthcare.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

TU-H-206-02: Novel Linearly-Filled Derenzo PET Phantom Design

Purpose: To design a linearly-filled Derenzo positron emission tomography (PET) phantom, eliminating the extraneous radioisotope volumes in a conventional reservoir-type design. This activity reduction combined with the elimination of bubbles in smaller phantom channels would significantly reduce personnel dose, radioisotope cost, and would improve image quality by reducing out-of-slice activity scatter. Methods: A computer-aided design (CAD) was created of a modular Derenzo phantom consisting of three phantom layers with gaskets between the layers. The central piece contains the active pattern volume and channels connecting adjacent rods in a serpentine pattern. The two end-pieces contained an inlet and an outlet for filling purposes. Phantom prototypes were 3D printed on a Viper Si2 stereolithography machine. The two gaskets were fabricated from silicon sheets using a PLS 6.75 laser cutter. Phantoms were held together by pass-through glass-filled nylon bolts and nuts. Phantoms were filled with {sup 52}Mn, {sup 64}Cu, {sup 74}Br, and {sup 124}I for testing, and were imaged on a Siemens Inveon MicroPET scanner. Results: Four phantom prototypes were constructed using male Leur Lock fittings for inlet/outlet ports. 3D printed layers were sanded to ensure proper coupling to the silicon gaskets. The filling volume for each prototype was approximately 2.4 mL. The filling process was found to be rapid, leak-tight, and with minimal back-pressure. PET images were reconstructed by OSEM3D, and axial slices along the phantom pattern length were averaged to provide final images. Image distortion was isotope dependent with {sup 52}Mn and {sup 64}Cu having the least distortion and {sup 124}I having the most distortion. Conclusion: These results indicate that the linearlyfilled Derenzo design improves on conventional reservoir-type designs by eliminating potential bubbles in small channels and by reducing activity level, radioisotope volume, radioisotope cost, personnel dose, filling time, and out-of-slice activity scatter The method described in this abstract has been filed as a patent application to the US Patent and Trade Office by the Wisconsin Alumni Research Foundation (WARF).

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

Characterizing Corrosion Mechanisms Of Structural Alloys In Actinide-based Molten Chloride Salt

This is a dataset consisting of 33 GB of .png image files, each featuring a two-dimensional distribution of different element species captured an individual slice of the specimen, and .ang files that capture microstructural information captured from each respective slice using a technique called electron backscatter diffraction. The .ang files can be processed to visualize alloy grains captured in a respective slice.

Copeland, Trishelle M.↗

DOE OSTI2025

Machine Learning Atom Probe Tomography Tool For Automatic And Fast Clustering

The software uses a YOLO11 segmentation model trained on synthetic data to analyze APT datasets. The workflow operates as follows: 1. Data Slicing: The APT dataset is divided into multiple 2D cross-sections of a specified thickness. 2. Segmentation: The model identifies point-dense regions within each 2D slice. 3. 3D Reconstruction: Detected regions (masks) from all slices are combined and reconstructed back into the original 3D space, forming clusters. The integration with HPC resources enables the software to process large-scale APT datasets efficiently. This combination of automation and scalability reduces manual intervention, improves reproducibility, and accelerates the clustering workflow.

Tang, Yalei [Idaho National Laboratory (INL), Idah↗

DOE OSTI2017

Small-angle X-ray scattering characteristics of mouse brain: Planar imaging measurements and tomographic imaging simulations

Small-angle x-ray scattering (SAXS) imaging can differentiate tissue types based on their nanoscale molecular structure. However, characterization of the coherent scattering cross-section profile of relevant tissues is needed to optimally design SAXS imaging techniques for a variety of biomedical applications. Reported measured nervous tissue x-ray scattering cross sections under a synchrotron source have had limited agreement. We report a set of x-ray cross-section measurements obtained from planar SAXS imaging of 1 mm thick mouse brain (APP/PS1 wild-type) coronal slices using an 8 keV laboratory x-ray source. Two characteristic peaks were found at 0.96 and 1.60 nm -1 attributed to myelin. The peak intensities varied by location in the slice. We found that regions of gray matter, white matter, and corpus callosum could be segmented by their increasing intensities of myelin peaks respectively. Measured small-angle x-ray scattering cross sections were then used to define brain tissue scattering properties in a GPU-accelerated Monte Carlo simulation of SAXS computed tomography (CT) using a higher monochromatic x-ray energy (20 keV) to study design trade-offs for noninvasive in vivo SAXS imaging on a small-animal head including radiation dose, signal-to-noise ratio (SNR), and the effect of skull presence on the previous two metrics. Simulation results show the estimated total dose to the mouse head for a single SAXS-CT slice was 149.4 mGy. The pixel SNR was approximately 30.8 for white matter material whether or not a skull was present. In this early-stage proof-of-principle work, we have demonstrated our brain cross-section data and simulation tools can be used to assess optimal instrument parameters for dedicated small-animal SAXS-CT prototypes.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2012

Local Longitudinal Microwave Instability Limits During Bunch Rotation

Bunch width compression can be accomplished by rf rotating an elongated bunch with minimal energy spread. The formation of tails at the two bunch ends produces disconnected regions in many time-advance slices. The stability limits of longitudinal microwave growth for each time-advance slice of the beam will be affected. This effect is studied and the Keil-Schnell stability limits for such time-advance slices are derived. Application is made to the bunch-width compression in the Fermilab Compressor Ring, destined for pion and subsequently muon production

43 PARTICLE ACCELERATORS↗