Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57

Imaging outside the box: Resolution enhancement in X-ray coherent diffraction imaging by extrapolation of diffraction patterns

Coherent diffraction imaging is a high-resolution imaging technique whose potential can be greatly enhanced by applying the extrapolation method presented here. We demonstrate the enhancement in resolution of a non-periodical object reconstructed from an experimental X-ray diffraction record which contains about 10% missing information, including the pixels in the center of the diffraction pattern. A diffraction pattern is extrapolated beyond the detector area and as a result, the object is reconstructed at an enhanced resolution and better agreement with experimental amplitudes is achieved. The optimal parameters for the iterative routine and the limits of the extrapolation procedure are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU-F-303-02: Achieving 4D MRI in Regular Breathing Cycle with Extended Acquisition Time of Dynamic MR Images

Purpose: The purpose of this study is to develop a technique to obtain four-dimensional (4D) magnetic resonance (MR) images that are more representative of a patient’s typical breathing cycle by utilizing an extended acquisition time while minimizing the image artifacts. Methods: The 4D MR data were acquired with the balanced steady state free precession in two-dimensional sagittal plane of view. Each slice was acquired repeatedly for about 15 s, thereby obtaining multiple images at each of the 10 phases in the respiratory cycle. This improves the probability that at least one of the images were acquired at the desired phase during a regular breathing cycle. To create optimal 4D MR images, an iterative approach was used to identify the set of images that yielded the highest slice-to-slice similarity. To assess the effectiveness of the approach, the data set was truncated into periods of 7 s (50 time points), 11 s (75 time points) and the full 15 s (100 time points). The 4D MR images were then sorted with data of the three different acquisition periods for comparison. Results: In general, the 4D MR images sorted using data from longer acquisition periods showed less mismatched artifacts. In addition, the normalized cross correlation (NCC) between slices of a 4D volume increases with increased acquisition period. The average NCC was 0.791 from the 7 s period, 0.794 from the 11 s period and 0.796 from the 15 s period. Conclusion: Our preliminary study showed that extending the acquisition time with the proposed sorting technique can improve image quality and reduce artifact presence in the 4D MR images. Data acquisition over two breathing cycles is a good trade-off between artifact reduction and scan time. This research was partially funded by the the Center for Radiation Oncology Research from UT MD Anderson Cancer Center.

60 APPLIED LIFE SCIENCES↗

Atomic Radius and Charge Parameter Uncertainty in Biomolecular Solvation Energy Calculations

Atomic radii and charges are two major parameters used in implicit solvent electrostatics and energy calculations. The optimization problem for charges and radii is under-determined, leading to uncertainty in the values of these parameters and in the results of solvation energy calculations using these parameters. This paper presents a method for quantifying this uncertainty in solvation energies using surrogate models based on generalized polynomial chaos (gPC) expansions. There are relatively few atom types used to specify radii parameters in implicit solvation calculations; therefore, surrogate models for these low-dimensional spaces could be constructed using least-squares fitting. However, there are many more types of atomic charges; therefore, construction of surrogate models for the charge parameter space required compressed sensing combined with an iterative rotation method to enhance problem sparsity. We present results for the uncertainty in small molecule solvation energies based on these approaches. Additionally, we explore the correlation between uncertainties due to radii and charges which motivates the need for future work in uncertainty quantification methods for high-dimensional parameter spaces.

Yang, Xiu↗

ISR-1 Seminar: Nuclear Cloud Lofting

This report goes into detail about motivation, background, lofting methodologies (empirical, parcel methods, Navier-Stokes), DELFIC Model, validation, DIORAMA integration, and has the following conclusions: A parcel methodology was applied to develop the DIORAMA cloud lofting module. It is based on the DELFIC model and treats the cloud as a homogenous unit. It solves the set of 8 ODEs for cloud properties. It outputs the time history of the cloud height, radius, and other parameters. The cloud lofting module was tuned with 54 test cases. An iterative brute force search was carried out to find the best fit tuning parameters. The best fit parameters yielded average cloud height errors of 12.9% and 10.3% for the bottom and top, respectively. This all allows for more accurate modeling of the propagation of delayed gamma rays in DIORAMA.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

SU-E-J-115: Correlation of Displacement Vector Fields Calculated by Deformable Image Registration Algorithms with Motion Parameters of CT Images with Well-Defined Targets and Controlled-Motion

Purpose: To investigate correlation of displacement vector fields (DVF) calculated by deformable image registration algorithms with motion parameters in helical axial and cone-beam CT images with motion artifacts. Methods: A mobile thorax phantom with well-known targets with different sizes that were made from water-equivalent material and inserted in foam to simulate lung lesions. The thorax phantom was imaged with helical, axial and cone-beam CT. The phantom was moved with a cyclic motion with different motion amplitudes and frequencies along the superior-inferior direction. Different deformable image registration algorithms including demons, fast demons, Horn-Shunck and iterative-optical-flow from the DIRART software were used to deform CT images for the phantom with different motion patterns. The CT images of the mobile phantom were deformed to CT images of the stationary phantom. Results: The values of displacement vectors calculated by deformable image registration algorithm correlated strongly with motion amplitude where large displacement vectors were calculated for CT images with large motion amplitudes. For example, the maximal displacement vectors were nearly equal to the motion amplitudes (5mm, 10mm or 20mm) at interfaces between the mobile targets lung tissue, while the minimal displacement vectors were nearly equal to negative the motion amplitudes. The maximal and minimal displacement vectors matched with edges of the blurred targets along the Z-axis (motion-direction), while DVF’s were small in the other directions. This indicates that the blurred edges by phantom motion were shifted largely to match with the actual target edge. These shifts were nearly equal to the motion amplitude. Conclusions: The DVF from deformable-image registration algorithms correlated well with motion amplitude of well-defined mobile targets. This can be used to extract motion parameters such as amplitude. However, as motion amplitudes increased, image artifacts increased significantly and that limited image quality and poor correlation between the motion amplitude and DVF was obtained.

60 APPLIED LIFE SCIENCES↗

Outer Wall Segmentation of Abdominal Aortic Aneurysm by Variable Neighborhood Search Through Intensity and Gradient Spaces

Aortic aneurysm segmentation remains a challenge. Manual segmentation is a time-consuming process which is not practical for routine use. To address this limitation, several automated segmentation techniques for aortic aneurysm have been developed, such as edge detection-based methods, partial differential equation methods, and graph partitioning methods. However, automatic segmentation of aortic aneurysm is difficult due to high pixel similarity to adjacent tissue and a lack of color information in the medical image, preventing previous work from being applicable to difficult cases. This paper uses uses a variable neighborhood search that alternates between intensity-based and gradient-based segmentation techniques. By alternating between intensity and gradient spaces, the search can escape from local optima of each space. The experimental results demonstrate that the proposed method outperforms the other existing segmentation methods in the literature, based on measurements of dice similarity coefficient and jaccard similarity coefficient at the pixel level. In addition, it is shown to perform well for cases that are difficult to segment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-E-T-128: Dosimetric Evaluation of MLC Modeling in Pinnacle V9.2 for Varian TrueBeam STx

Purpose: To evaluate the effects of MLC modeling after commissioning the Varian TrueBeam LINAC in Pinnacle version 9.2. Methods: Stepand-shoot IMRT QAs were investigated when we observed our measured absolute dose results using ion chamber (Capintec PR-05P) were uncharacteristically low; about 4–5% compared to doses calculated by Pinnacle{sup 3} (Phillips, Madison, WI). This problem was predominant for large and highly modulated head and neck (HN) treatments. Intuitively we knew this had to be related to shortcomings in the MLC modeling in Pinnacle. Using film QA we were able to iteratively adjust the MLC parameters. We confirmed results by re-testing five failed IMRT QA patients; and ion chamber measurements were verified in Quasar anthropomorphic phantom. Results: After commissioning the LINAC in Pinnacle version 9.2, the MLC transmission for 6X, 10X and 15X were 2.0%, 1.7% and 2.0%, respectively, and additional Interleaf leakage for all three energies was 0.5%. These parameters were obtained from profiles scanned with an Edge detector (Sun Nuclear, Melbourne, FL) during machine commissioning. A Verification testing with radiographic EDR2 film (Kodak, Rochester, NY) measurement was performed by creating a closed MLC leaf pattern and analyzing using RIT software (RIT, Colorado Springs, CO). This reduced MLC transmission for 6X, 10X and 15X to 0.7%, 0.9% and 0.9%, respectively; while increasing additional Interleaf leakage for all three energies to 1.0%. Conclusion: Radiographic film measurements were used to correct MLC transmission values for step and shoot IMRT fields used in Pinnacle version 9.2. After adjusting the MLC parameters to correlate with the film QA, there was still very good agreement between the Pinnacle model and commissioning data. Using the same QA methodology, we were also able to improve the beam models for the Varian C-series linacs, Novalis-Tx, and TrueBeam M-120 linacs.

07 ISOTOPE AND RADIATION SOURCES↗

SU-E-T-478: IMRT Delivery Parameter Dependence of Dose-Mass Optimization

Purpose: To compare DMH and DVH optimization sensitivity to changes in IMRT delivery parameters. Methods: Two lung and two head and neck (HN) cases were retrospectively optimized using DVH and DMH optimization. For both optimization approaches, changes to two parameters were studied: number of IMRT segments (5 and 10 per beam) and the minimum segment area (2 and 6 cm2). The number of beams, beam angles, and minimum MUs per segment were the same for both optimizations approaches for each patient. During optimization, doses to the organs at risk (OARs) were iteratively lowered until the standard deviation across the PTV was above ∼3.0%. For each patient DVH and DMH plans were normalized such that 95% of the PTV received the same dose. Plan quality was evaluated by dose indices (DIs), which represent the dose delivered to a certain anatomical structure volume. For the lung cases, DIs assessed included: 1% cord, 33% heart, both lungs 20% and 30%, and 50% esophagus. In the HN cases: 1% cord, 1% brainstem, left/right parotids 50%, 50% larynx, and 50% esophagus. Results: When increasing the number of segments, while keeping a small segment area (2cm2), the average percent change of all DIs for DVH/DMH optimizations for each patient were: −4.66/4.71, 3.21/3.46, −9.62/21.69 and −3.28/−7.62. For a large segment area (6cm2): −0.26/−1.46, −5.04/−1.92, −5.23/−2.19 and 4.12/19.63. Results from increasing segment area while keeping a small number of segments (5segments/beam) were: 1.41/7.90, 8.17/11.66, 0.09/33.58 and −4.83/−11.60 for each case. For large number of segments (10 segments/beam): 8.35/1.30, −0.91/5.77, 6.29/7.08 and 2.62/5.16. Conclusion: This preliminary study showed case dependent results. Changes in IMRT parameters did not show consistent DI changes for either optimization approach. A larger population of patients is warranted for such comparison.

60 APPLIED LIFE SCIENCES↗

TH-AB-202-12: The First Clinical Implementation of a Real-Time Six Degree of Freedom Tracking System During Radiation Therapy

Purpose: In current practice, imaging is typically performed prior to treatment; the cancer target motion during treatment is unknown. We present the first clinical implementation of real time Kilovoltage Intrafraction Monitoring (KIM) system which tracks the cancer target translational and rotational motions during treatment. Methods: KIM technology: KIM estimates the 3D position of the target tumour based on segmented 2D positions of the three implanted fiducials in each of the kV images (125 kV, 10 mA at 11 fps) taken continuously during the treatment arc. The 2D-3D target estimation is based on a probability distribution function, obtained during pre-treatment CBCT. Rotations about each axis with the centroid of the markers as the pivot were calculated using the iterative closest point algorithm in real time. Patient: A patient with prostate adenocarcinoma undergoing stereotactic body radiotherapy (SBRT) with 36.25 Gy delivered in 5 fractions (Varian TrueBeam, 6X, VMAT) was enrolled in the study. The trial complies with Australian ethical and regulatory standards. Results: Of the 5 fractions of treatment the patient received, KIM was utilised successfully in 4 fractions with 3 couch shifts due to large persistent prostate movements (>2mm for more than 5 seconds). KIM translational accuracy and precision in comparison with post treatment kV-MV triangulation are 0.28±0.59 mm, −0.19±0.25 mm and 0.23±0.69 mm in the Left-Right, Superior-Inferior and Anterior-Posterior directions, respectively. KIM rotational accuracy as compared with triangulation is: 0.429°±2.22°, −0.44°±4.7° and 0.06°±1.08° in the roll, pitch and yaw direction, respectively. Conclusion: The first six degree of freedom KIM system was successfully implemented clinically. The presented KIM system has sub-millimeters accuracy and precision in all three translational axes, and less than 1° of mean error in all three rotational axes. Acknowledgement: This work is supported by Cancer Australia grants APPXXX, APPYYY.

60 APPLIED LIFE SCIENCES↗

Integral Criticality Estimators in MCATK

The Monte Carlo Application ToolKit (MCATK) is a component-based software toolset for delivering customized particle transport solutions using the Monte Carlo method. Currently under development in the XCP Monte Carlo group at Los Alamos National Laboratory, the toolkit has the ability to estimate the ke f f and a eigenvalues for static geometries. This paper presents a description of the estimators and variance reduction techniques available in the toolkit and includes a preview of those slated for future releases. Along with the description of the underlying algorithms is a description of the available user inputs for controlling the iterations. The paper concludes with a comparison of the MCATK results with those provided by analytic solutions. The results match within expected statistical uncertainties and demonstrate MCATK’s usefulness in estimating these important quantities.

97 MATHEMATICS AND COMPUTING↗

Numerical analysis of flows of rarefied gases in long channels with octagonal cross section shapes

Isothermal, pressure driven rarefied gas flows through long channels with octagonal cross section shapes are analyzed computationally. The capillary is between inlet and outlet reservoirs. The cross section is constant along the axial direction. The boundary condition at the solid-gas interface is assumed to be diffuse reflection. Since the channel is long, the gaseous velocity is small compared to the average molecular speed. Consequently, a linearized description can be used. The flow is described by the linearized Bhatnagar-Gross-Krook kinetic model. The solution of the problem is divided into two stages. First, the local flow field is determined by assuming the local pressure gradient. Secondly, the global flow behavior is deduced by the consideration of the conservation of the mass along the axis of the capillary. The kinetic equation is solved by the discrete velocity method on the cross section. Both spatial and velocity spaces are discretized. A body fitted rectangular grid is used for the spatial space. Near the boundary, first-order, while in the interior part of the flow domain, second-order finite-differences are applied to approximate the spatial derivatives. This combination results into an efficient and straightforward numerical treatment. The velocity space is represented by a Gauss-Legendre quadrature. The kinetic equation is solved in an iterative manner. The local dimensionless flow rate is calculated and tabulated for a wide range of the gaseous rarefaction for octagonal cross sections with various geometrical parameters. It exhibits the Knudsen minimum phenomenon. The flow rates in the octagonal channel are compared to those through capillaries with circular and square cross sections. Typical velocity profiles are also shown. The mass flow rate and the distribution of the pressure are determined and presented for global pressure driven flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A multilevel method for coupling the neutron kinetics and heat transfer equations

We present a computational method for adequate and efficient coupling the neutron transport equation with the precursor and heat transfer equations. It is based on the multilevel nonlinear quasidiffusion (QD) method for solving the multigroup transport equation. The system of equations includes the time-dependent high-order transport equation and time-dependent multigroup and effective one-group low-order QD equations. We also formulate a method applying the {alpha}-approximation for the time-dependent high-order transport equation. This approach enables one to avoid storing the angular flux from the previous time step. Numerical results for a model transient problem are presented. (authors)

97 MATHEMATICS AND COMPUTING↗

A DATA-DRIVEN MODEL FOR SPECTRA: FINDING DOUBLE REDSHIFTS IN THE SLOAN DIGITAL SKY SURVEY

We present a data-driven method-heteroscedastic matrix factorization, a kind of probabilistic factor analysis-for modeling or performing dimensionality reduction on observed spectra or other high-dimensional data with known but non-uniform observational uncertainties. The method uses an iterative inverse-variance-weighted least-squares minimization procedure to generate a best set of basis functions. The method is similar to principal components analysis (PCA), but with the substantial advantage that it uses measurement uncertainties in a responsible way and accounts naturally for poorly measured and missing data; it models the variance in the noise-deconvolved data space. A regularization can be applied, in the form of a smoothness prior (inspired by Gaussian processes) or a non-negative constraint, without making the method prohibitively slow. Because the method optimizes a justified scalar (related to the likelihood), the basis provides a better fit to the data in a probabilistic sense than any PCA basis. We test the method on Sloan Digital Sky Survey (SDSS) spectra, concentrating on spectra known to contain two redshift components: these are spectra of gravitational lens candidates and massive black hole binaries. We apply a hypothesis test to compare one-redshift and two-redshift models for these spectra, utilizing the data-driven model trained on a random subset of all SDSS spectra. This test confirms 129 of the 131 lens candidates in our sample and all of the known binary candidates, and turns up very few false positives.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparison of analytic models of instability of rarefied gas flow in a channel

Numerical and analytical results are compared concerning the limit properties of the trajectories, attractors and bifurcations of rarefied gas flows in channels. The cascade of bifurcations obtained in our previous analytical and numerical investigations is simulated numerically for different scattering functions V generalizing the ray-diffuse reflection of gas particles from the surface. The main purpose of numerical simulation by Monte Carlo method is the investigation of the properties of different analytic nonlinear dynamic systems corresponding to rarefied gas flow in a channel. The results are compared as well for the models suggested originally by R. N. Miroshin, as for the approximations considered for the first time or for studied in our subsequent papers. Analytical solutions we obtained earlier for the ray reflection which means only one determined velocity of scattered from the walls gas atoms, generally different from the specular reflection. The nonlinear iterative equation describing a rarefied gas flow in a long channel becomes unstable in some regions of corresponding parameters of V (it means the sensitivity to boundary conditions). The values of the parameters are found from analytical approximations in these regions. Numerical results show that the chaotic behavior of the nonlinear dynamic system corresponds to strange attractors and distinguishes clearly from Maxwellian distribution and from the equilibrium on the whole. In the regions of instability (as the dimension of the attractor increases) the search for a corresponding state requires a lot more computation time and a lot of data (the amount of data required increases exponentially with embedding dimension). Therefore the main complication in the computation is reducing as well the computing time as the amount of data to find a suitably close solution. To reduce the computing time our analytical results are applied. Flow conditions satisfying the requirements to the experiment are indicated where the instability of considered type can be detected. The advantages and the drawbacks of considered approximations are detected and the features of obtained solutions are indicated. The recommendations are given for applying the results in practical applications and in numerical calculations of rarefied gas flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A digitally reconstructed radiograph algorithm calculated from first principles

Purpose: To develop an algorithm for computing realistic digitally reconstructed radiographs (DRRs) that match real cone-beam CT (CBCT) projections with no artificial adjustments. Methods: The authors used measured attenuation data from cone-beam CT projection radiographs of different materials to obtain a function to convert CT number to linear attenuation coefficient (LAC). The effects of scatter, beam hardening, and veiling glare were first removed from the attenuation data. Using this conversion function the authors calculated the line integral of LAC through a CT along rays connecting the radiation source and detector pixels with a ray-tracing algorithm, producing raw DRRs. The effects of scatter, beam hardening, and veiling glare were then included in the DRRs through postprocessing. Results: The authors compared actual CBCT projections to DRRs produced with all corrections (scatter, beam hardening, and veiling glare) and to uncorrected DRRs. Algorithm accuracy was assessed through visual comparison of projections and DRRs, pixel intensity comparisons, intensity histogram comparisons, and correlation plots of DRR-to-projection pixel intensities. In general, the fully corrected algorithm provided a small but nontrivial improvement in accuracy over the uncorrected algorithm. The authors also investigated both measurement- and computation-based methods for determining the beam hardening correction, and found the computation-based method to be superior, as it accounted for nonuniform bowtie filter thickness. The authors benchmarked the algorithm for speed and found that it produced DRRs in about 0.35 s for full detector and CT resolution at a ray step-size of 0.5 mm. Conclusions: The authors have demonstrated a DRR algorithm calculated from first principles that accounts for scatter, beam hardening, and veiling glare in order to produce accurate DRRs. The algorithm is computationally efficient, making it a good candidate for iterative CT reconstruction techniques that require a data fidelity term based on the matching of DRRs and projections.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU-F-T-195: Systematic Constraining of Contralateral Parotid Gland Led to Improved Dosimetric Outcomes for Multi-Field Optimization with Scanning Beam Proton Therapy: Promising Results From a Pilot Study in Patients with Base of Tongue Carcinoma

Purpose: Treatment planning for Intensity Modulated Proton Therapy (IMPT) for head and neck cancer is time-consuming due to the large number of organs-at-risk (OAR) to be considered. As there are many competing objectives and also wide range of acceptable OAR constraints, the final approved plan may not be most optimal for the given structures. We evaluated the dose reduction to the contralateral parotid by implementing standardized constraints during optimization for scanning beam proton therapy planning. Methods: Twenty-four (24) consecutive patients previously treated for base of tongue carcinoma were retrospectively selected. The doses were 70Gy, 63Gy and 57Gy (SIB in 33 fractions) for high-, intermediate-, and standard-risk clinical target volumes (CTV), respectively; the treatment included bilateral neck. Scanning beams using MFO with standardized bilateral anterior oblique and PA fields were applied. New plans where then developed and optimized by employing additional contralateral parotid constraints at multiple defined dose levels. Using a step-wise iterative process, the volume-based constraints at each level were then further reduced until known target coverages were compromised. The newly developed plans were then compared to the original clinically approved plans using paired student t-testing. Results: All 24 newly optimized treatment plans maintained initial plan quality as compared to the approved plans, and the 98% prescription dose coverage to the CTV’s were not compromised. Representative DVH comparison is shown in FIGURE 1. The contralateral parotid doses were reduced at all levels of interest when systematic constraints were applied to V10, V20, V30 and V40Gy (All P<0.0001; TABLE 1). Overall, the mean contralateral parotid doses were reduced by 2.26 Gy on average, a ∼13% relative improvement. Conclusion: Applying systematic and volume-based contralateral parotid constraints for IMPT planning significantly reduced the dose at all dosimetric levels for patients with base of tongue cancer.

60 APPLIED LIFE SCIENCES↗

Methodology for Sensitivity Analysis, Approximate Analysis, and Design Optimization in CFD for Multidisciplinary Applications

An incremental iterative formulation together with the well-known spatially split approximate-factorization algorithm, is presented for solving the large, sparse systems of linear equations that are associated with aerodynamic sensitivity analysis. This formulation is also known as the 'delta' or 'correction' form. For the smaller two dimensional problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. However, iterative methods are needed for larger two-dimensional and three dimensional applications because direct methods require more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioned coefficient matrix; this problem is overcome when these equations are cast in the incremental form. The methodology is successfully implemented and tested using an upwind cell-centered finite-volume formulation applied in two dimensions to the thin-layer Navier-Stokes equations for external flow over an airfoil. In three dimensions this methodology is demonstrated with a marching-solution algorithm for the Euler equations to calculate supersonic flow over the High-Speed Civil Transport configuration (HSCT 24E). The sensitivity derivatives obtained with the incremental iterative method from a marching Euler code are used in a design-improvement study of the HSCT configuration that involves thickness. camber, and planform design variables.

Taylor, Arthur C., III↗

Assessment of volumetric noise and resolution performance for linear and nonlinear CT reconstruction methods

Purpose: For nonlinear iterative image reconstructions (IR), the computed tomography (CT) noise and resolution properties can depend on the specific imaging conditions, such as lesion contrast and image noise level. Therefore, it is imperative to develop a reliable method to measure the noise and resolution properties under clinically relevant conditions. This study aimed to develop a robust methodology to measure the three-dimensional CT noise and resolution properties under such conditions and to provide guidelines to achieve desirable levels of accuracy and precision. Methods: The methodology was developed based on a previously reported CT image quality phantom. In this methodology, CT noise properties are measured in the uniform region of the phantom in terms of a task-based 3D noise-power spectrum (NPS{sub task}). The in-plane resolution properties are measured in terms of the task transfer function (TTF) by applying a radial edge technique to the rod inserts in the phantom. The z-direction resolution properties are measured from a supplemental phantom, also in terms of the TTF. To account for the possible nonlinearity of IR, the NPS{sub task} is measured with respect to the noise magnitude, and the TTF with respect to noise magnitude and edge contrast. To determine the accuracy and precision of the methodology, images of known noise and resolution properties were simulated. The NPS{sub task} and TTF were measured on the simulated images and compared to the truth, with criteria established to achieve NPS{sub task} and TTF measurements with <10% error. To demonstrate the utility of this methodology, measurements were performed on a commercial CT system using five dose levels, two slice thicknesses, and three reconstruction algorithms (filtered backprojection, FBP; iterative reconstruction in imaging space, IRIS; and sinogram affirmed iterative reconstruction with strengths of 5, SAFIRE5). Results: To achieve NPS{sub task} measurements with <10% error, the number of regions of interest needed to be greater than 65. To achieve TTF measurements with <10% error, the contrast-to-noise ratio of the edge needed to be ≥15, achievable by averaging multiple slices across the same edge. The NPS{sub task} measured on a commercial CT system showed IR's reduced noise (IRIS, 30% and SAFIRE5, 55%) and “waxier” texture (peak frequencies: FBP, 0.25 mm{sup −1}; IRIS, 0.23 mm{sup −1}; and SAFIRE5, 0.16 mm{sup −1}). The TTF measured within the axial plane showed improved in-plane resolution with SAFIRE5 at the TTF 50% frequency, f{sub 50} (FBP, 0.36–0.41 mm{sup −1}; SAFIRE5, 0.37–0.46 mm{sup −1}). The TTF measured along the axial direction showed improved z-direction resolution with thinner slice thickness (f{sub 50}: 0.6 mm, 0.35–0.79 mm{sup −1}; 1.5 mm, 0.22–0.3 mm{sup −1}) and with SAFIRE5 (f{sub 50}: FBP, 0.35–0.52 mm{sup −1}; SAFIRE5, 0.42–0.79 mm{sup −1}). Both in-plane and z-direction resolution of SAFIRE5 showed strong dependency on contrast, reflecting SAFIRE5's nonlinearity. Conclusions: A methodology was developed to measure three-dimensional CT noise and resolution properties for iterative reconstruction, especially at challenging measurement conditions with low contrast and high image noise. The methodology also demonstrated its utility for evaluating commercial CT systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗