Search NASA⌕ Search

SEARCH · Search NASA

Results for “Digital image quality phantom”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-J-27: Shifting Multiple EPID Imager Layers to Improve Image Quality and Resolution in MV CBCT

Purpose: Vertical stacking of four conventional EPID layers can improve DQE for MV-CBCT applications. We hypothesize that shifting each layer laterally by half a pixel relative to the layer above, will improve the contrast-to-noise ratio (CNR) and image resolution. Methods: For CNR assessment, a 20 cm diameter digital phantom with 8 inserts is created. The attenuation coefficient of the phantom is similar to lung at the average energy of a 6 MV photon beam. The inserts have attenuations 1, 2…8 times of lung. One of the inserts is close to soft tissue, resembling the case of a tumor in lung. For resolution assessment, a digital phantom featuring a bar pattern is created. The phantom has an attenuation coefficient similar to soft tissue and the bars have an attenuation coefficient of calcium sulfate. A 2 MeV photon beam is attenuated through these phantoms and hits each of the four stacked detector layers. Each successive layer is shifted by half a pixel in the x only, y only, and x and y (combined) directions, respectively. Blurring and statistical noise are added to the projections. Projections from one, two, three and four layers are used for reconstruction. CNR and image resolution are evaluated and compared. Results: When projections from multiple layers are combined for reconstruction, CNR increases with the number of layers involved. CNR in reconstructions from two, three and four layers are 1.4, 1.7 and 1.99 times that from one layer. The resolution from the shifted four layer detector is also improved from a single layer. In a comparison between one layer versus four layers in this preliminary study, the resolution from four shifted layers is at least 20% better. Conclusion: Layer-shifting in a stacked EPID imager design enhances resolution as well as CNR for half scan MV-CBCT. The project described was supported, in part, by a grant from Varian Medical Systems, Inc., and Award No. R01CA188446-01 from the National Cancer Institute. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Cancer Institute or the National Institutes of Health.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-89: Comparative Analysis of MIM and Velocity’s Image Deformation Algorithm Using Simulated KV-CBCT Images for Quality Assurance

Purpose: Deformable image registration (DIR) is used routinely in the clinic without a formalized quality assurance (QA) process. Using simulated deformations to digitally deform images in a known way and comparing to DIR algorithm predictions is a powerful technique for DIR QA. This technique must also simulate realistic image noise and artifacts, especially between modalities. This study developed an algorithm to create simulated daily kV cone-beam computed-tomography (CBCT) images from CT images for DIR QA between these modalities. Methods: A Catphan and physical head-and-neck phantom, with known deformations, were used. CT and kV-CBCT images of the Catphan were utilized to characterize the changes in Hounsfield units, noise, and image cupping that occur between these imaging modalities. The algorithm then imprinted these changes onto a CT image of the deformed head-and-neck phantom, thereby creating a simulated-CBCT image. CT and kV-CBCT images of the undeformed and deformed head-and-neck phantom were also acquired. The Velocity and MIM DIR algorithms were applied between the undeformed CT image and each of the deformed CT, CBCT, and simulated-CBCT images to obtain predicted deformations. The error between the known and predicted deformations was used as a metric to evaluate the quality of the simulated-CBCT image. Ideally, the simulated-CBCT image registration would produce the same accuracy as the deformed CBCT image registration. Results: For Velocity, the mean error was 1.4 mm for the CT-CT registration, 1.7 mm for the CT-CBCT registration, and 1.4 mm for the CT-simulated-CBCT registration. These same numbers were 1.5, 4.5, and 5.9 mm, respectively, for MIM. Conclusion: All cases produced similar accuracy for Velocity. MIM produced similar values of accuracy for CT-CT registration, but was not as accurate for CT-CBCT registrations. The MIM simulated-CBCT registration followed this same trend, but overestimated MIM DIR errors relative to the CT-CBCT registration.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-C-209-07: Phantoms for Digital Breast Tomosynthesis Imaging System Evaluation

Purpose: Digital Breast Tomosynthesis (DBT) is gaining importance in breast imaging. There is a need for phantoms that can be used for image evaluation and comparison. Existing commercially available phantoms for DBT are expensive and may lack clinically relevant test objects. The purpose of this study is to develop phantoms for DBT evaluation. Methods Four phantoms have been designed and constructed to assess the image quality (IQ) of two DBT systems. The first contains a spiral of 0.3 mm SiC beads in gelatin to measure the tomographic slice thickness profile and uniformity of coverage in a series of tomographic planes. The second contains simulated tumors inclined with respect to the phantom base to assess tomographic image quality. The third has a tilted array of discs with varying contrast and diameter. This phantom was imaged alone and in a stack of TE slabs giving 2 to 10 cm thickness. The fourth has a dual wedge of glandular and adipose simulating materials. One wedge contains discs with varying diameter and thickness; the other supports a mass with six simulated spicules of varying size and a cluster of simulated calcifications. The simulated glandular tissue material varies between 35 and 100% of the total thickness (5.5 cm). Results: All phantoms were scanned successfully. The best IQ comparison was achieved with the dual wedge phantom as demonstrated by the spiculated mass and calcifications. Images were evaluated by two radiologists and one physicist. The projection images and corresponding set of tomographic planes were comparable and the synthesized projection images were inferior to the projection images for both systems. Conclusion: Four phantoms were designed, constructed and imaged on two DBT systems. They successfully demonstrated performance differences between two systems, and between true and synthesized projection images. Future work will incorporate these designs into a single phantom.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-BRB-01: A Novel 3D Printed Patient-Specific Phantom for Spine SBRT Quality Assurance: Comparison of 3D Printing Techniques

Purpose: The novel 3 dimensional (3D)-printed spine quality assurance (QA) phantoms generated by two different 3D-printing technologies, digital light processing (DLP) and Polyjet, were developed and evaluated for spine stereotactic body radiation treatment (SBRT). Methods: The developed 3D-printed spine QA phantom consisted of an acrylic body and a 3D-printed spine phantom. DLP and Polyjet 3D printers using the high-density acrylic polymer were employed to produce spine-shaped phantoms based on CT images. To verify dosimetric effects, the novel phantom was made it enable to insert films between each slabs of acrylic body phantom. Also, for measuring internal dose of spine, 3D-printed spine phantom was designed as divided laterally exactly in half. Image fusion was performed to evaluate the reproducibility of our phantom, and the Hounsfield unit (HU) was measured based on each CT image. Intensity-modulated radiotherapy plans to deliver a fraction of a 16 Gy dose to a planning target volume (PTV) based on the two 3D-printing techniques were compared for target coverage and normal organ-sparing. Results: Image fusion demonstrated good reproducibility of the fabricated spine QA phantom. The HU values of the DLP- and Polyjet-printed spine vertebrae differed by 54.3 on average. The PTV Dmax dose for the DLP-generated phantom was about 1.488 Gy higher than for the Polyjet-generated phantom. The organs at risk received a lower dose when the DLP technique was used than when the Polyjet technique was used. Conclusion: This study confirmed that a novel 3D-printed phantom mimicking a high-density organ can be created based on CT images, and that a developed 3D-printed spine phantom could be utilized in patient-specific QA for SBRT. Despite using the same main material, DLP and Polyjet yielded different HU values. Therefore, the printing technique and materials must be carefully chosen in order to accurately produce a patient-specific QA phantom.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Cerebral perfusion computed tomography deconvolution via structure tensor total variation regularization

Purpose: Cerebral perfusion computed tomography (PCT) imaging as an accurate and fast acute ischemic stroke examination has been widely used in clinic. Meanwhile, a major drawback of PCT imaging is the high radiation dose due to its dynamic scan protocol. The purpose of this work is to develop a robust perfusion deconvolution approach via structure tensor total variation (STV) regularization (PD-STV) for estimating an accurate residue function in PCT imaging with the low-milliampere-seconds (low-mAs) data acquisition. Methods: Besides modeling the spatio-temporal structure information of PCT data, the STV regularization of the present PD-STV approach can utilize the higher order derivatives of the residue function to enhance denoising performance. To minimize the objective function, the authors propose an effective iterative algorithm with a shrinkage/thresholding scheme. A simulation study on a digital brain perfusion phantom and a clinical study on an old infarction patient were conducted to validate and evaluate the performance of the present PD-STV approach. Results: In the digital phantom study, visual inspection and quantitative metrics (i.e., the normalized mean square error, the peak signal-to-noise ratio, and the universal quality index) assessments demonstrated that the PD-STV approach outperformed other existing approaches in terms of the performance of noise-induced artifacts reduction and accurate perfusion hemodynamic maps (PHM) estimation. In the patient data study, the present PD-STV approach could yield accurate PHM estimation with several noticeable gains over other existing approaches in terms of visual inspection and correlation analysis. Conclusions: This study demonstrated the feasibility and efficacy of the present PD-STV approach in utilizing STV regularization to improve the accuracy of residue function estimation of cerebral PCT imaging in the case of low-mAs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-AB-202-05: GPU-Based 4D Deformable Image Registration Using Adaptive Tetrahedral Mesh Modeling

Purpose: Deformable image registration (DIR) has been employed today as an automated and effective segmentation method to transfer tumor or organ contours from the planning image to daily images, instead of manual segmentation. However, the computational time and accuracy of current DIR approaches are still insufficient for online adaptive radiation therapy (ART), which requires real-time and high-quality image segmentation, especially in a large datasets of 4D-CT images. The objective of this work is to propose a new DIR algorithm, with fast computational speed and high accuracy, by using adaptive feature-based tetrahedral meshing and GPU-based parallelization. Methods: The first step is to generate the adaptive tetrahedral mesh based on the image features of a reference phase of 4D-CT, so that the deformation can be well captured and accurately diffused from the mesh vertices to voxels of the image volume. Subsequently, the deformation vector fields (DVF) and other phases of 4D-CT can be obtained by matching each phase of the target 4D-CT images with the corresponding deformed reference phase. The proposed 4D DIR method is implemented on GPU, resulting in significantly increasing the computational efficiency due to its parallel computing ability. Results: A 4D NCAT digital phantom was used to test the efficiency and accuracy of our method. Both the image and DVF results show that the fine structures and shapes of lung are well preserved, and the tumor position is well captured, i.e., 3D distance error is 1.14 mm. Compared to the previous voxel-based CPU implementation of DIR, such as demons, the proposed method is about 160x faster for registering a 10-phase 4D-CT with a phase dimension of 256×256×150. Conclusion: The proposed 4D DIR method uses feature-based mesh and GPU-based parallelism, which demonstrates the capability to compute both high-quality image and motion results, with significant improvement on the computational speed.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Evaluation of the spline reconstruction technique for PET

Purpose: The spline reconstruction technique (SRT), based on the analytic formula for the inverse Radon transform, has been presented earlier in the literature. In this study, the authors present an improved formulation and numerical implementation of this algorithm and evaluate it in comparison to filtered backprojection (FBP). Methods: The SRT is based on the numerical evaluation of the Hilbert transform of the sinogram via an approximation in terms of “custom made” cubic splines. By restricting reconstruction only within object pixels and by utilizing certain mathematical symmetries, the authors achieve a reconstruction time comparable to that of FBP. The authors have implemented SRT in STIR and have evaluated this technique using simulated data from a clinical positron emission tomography (PET) system, as well as real data obtained from clinical and preclinical PET scanners. For the simulation studies, the authors have simulated sinograms of a point-source and three digital phantoms. Using these sinograms, the authors have created realizations of Poisson noise at five noise levels. In addition to visual comparisons of the reconstructed images, the authors have determined contrast and bias for different regions of the phantoms as a function of noise level. For the real-data studies, sinograms of an{sup 18}F-FDG injected mouse, a NEMA NU 4-2008 image quality phantom, and a Derenzo phantom have been acquired from a commercial PET system. The authors have determined: (a) coefficient of variations (COV) and contrast from the NEMA phantom, (b) contrast for the various sections of the Derenzo phantom, and (c) line profiles for the Derenzo phantom. Furthermore, the authors have acquired sinograms from a whole-body PET scan of an {sup 18}F-FDG injected cancer patient, using the GE Discovery ST PET/CT system. SRT and FBP reconstructions of the thorax have been visually evaluated. Results: The results indicate an improvement in FWHM and FWTM in both simulated and real point-source studies. In all simulated phantoms, the SRT exhibits higher contrast and lower bias than FBP at all noise levels, by increasing the COV in the reconstructed images. Finally, in real studies, whereas the contrast of the cold chambers are similar for both algorithms, the SRT reconstructed images of the NEMA phantom exhibit slightly higher COV values than those of FBP. In the Derenzo phantom, SRT resolves the 2-mm separated holes slightly better than FBP. The small-animal and human reconstructions via SRT exhibit slightly higher resolution and contrast than the FBP reconstructions. Conclusions: The SRT provides images of higher resolution, higher contrast, and lower bias than FBP, by increasing slightly the noise in the reconstructed images. Furthermore, it eliminates streak artifacts outside the object boundary. Unlike other analytic algorithms, the reconstruction time of SRT is comparable with that of FBP. The source code for SRT will become available in a future release of STIR.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

WE-EF-303-04: An Advanced Image Processing Method to Improve the Spatial Resolution of Proton Radiographies

Purpose: We present an optimization method to improve the spatial resolution and the water equivalent thickness accuracy of proton radiographies. Methods: The method is designed for imaging systems measuring only the residual range of protons without relying on tracker detectors to determine the beam trajectory before and after the target. Specifically, the method was used for an imaging set-up consisting of a stack of 61 parallel-plate ionization chambers (PPIC) working as a range telescope. The method uses a decomposition approach of the residual range signal measured by the PPIC and constructs subimages with small size pixels geometrically rearranged and appropriately averaged to be merged into a final single radiography. The method was tested using Monte Carlo simulated and experimental proton radiographies of a PMMA step phantom and an anthropomorphic head phantom. Results: For the step phantom, the effective spatial resolution was found to be 4 and 3 times higher than the nominal resolution for the simulated and experimental radiographies, respectively. For the head phantom, a gamma index was calculated to quantify the conformity of the simulated proton radiographies with a digitally reconstructed X-ray radiography convolved with a Gaussian kernel equal to the proton beam spot-size. For DTA=2.5 mm and RD=2.5%, the passing ratio was 100%/85% for the optimized/non-optimized case, respectively. An extension of the method allows reducing the dose given to the patient during radiography acquisition. We show that despite a dose reduction of 25 times (leading to a dose of 0.016 mGy for the current imaging set-up), the image quality of the optimized radiographies remains fairly unaffected for both the simulated and experimental results. Conclusion: The optimization method leads to a significant increase of the spatial resolution allowing recovering image details that are unresolved in non-optimized radiographies. These results represent a major step towards clinical application of proton radiography.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Four-dimensional dose reconstruction through in vivo phase matching of cine images of electronic portal imaging device

Purpose: A method is proposed to reconstruct a four-dimensional (4D) dose distribution using phase matching of measured cine images to precalculated images of electronic portal imaging device (EPID). Methods: (1) A phantom, designed to simulate a tumor in lung (a polystyrene block with a 3 cm diameter embedded in cork), was placed on a sinusoidally moving platform with an amplitude of 1 cm and a period of 4 s. Ten-phase 4D computed tomography (CT) images of the phantom were acquired. A planning target volume (PTV) was created by adding a margin of 1 cm around the internal target volume of the tumor. (2) Three beams were designed, which included a static beam, a theoretical dynamic beam, and a planning-optimized dynamic beam (PODB). While the theoretical beam was made by manually programming a simplistic sliding leaf motion, the planning-optimized beam was obtained from treatment planning. From the three beams, three-dimensional (3D) doses on the phantom were calculated; 4D dose was calculated by means of the ten phase images (integrated over phases afterward); serving as “reference” images, phase-specific EPID dose images under the lung phantom were also calculated for each of the ten phases. (3) Cine EPID images were acquired while the beams were irradiated to the moving phantom. (4) Each cine image was phase-matched to a phase-specific CT image at which common irradiation occurred by intercomparing the cine image with the reference images. (5) Each cine image was used to reconstruct dose in the phase-matched CT image, and the reconstructed doses were summed over all phases. (6) The summation was compared with forwardly calculated 4D and 3D dose distributions. Accounting for realistic situations, intratreatment breathing irregularity was simulated by assuming an amplitude of 0.5 cm for the phantom during a portion of breathing trace in which the phase matching could not be performed. Intertreatment breathing irregularity between the time of treatment and the time of planning CT was considered by utilizing the same reduced amplitude when the phantom was irradiated. To examine the phase matching in a humanoid environment, the matching was also performed in a digital phantom (4D XCAT phantom). Results: For the static, the theoretical, and the planning-optimized dynamic beams, the 4D reconstructed doses showed agreement with the forwardly calculated 4D doses within the gamma pass rates of 92.7%, 100%, and 98.1%, respectively, at the isocenter plane given by 3%/3 mm criteria. Excellent agreement in dose volume histogram of PTV and lung-PTV was also found between the two 4D doses, while substantial differences were found between the 3D and the 4D doses. The significant breathing irregularities modeled in this study were found not to be noticeably affecting the reconstructed dose. The phase matching was performed equally well in a digital phantom. Conclusions: The method of retrospective phase determination of a moving object under irradiation provided successful 4D dose reconstruction. This method will provide accurate quality assurance and facilitate adaptive therapy when distinguishable objects such as well-defined tumors, diaphragm, and organs with markers (pancreas and liver) are covered by treatment beam apertures.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Monte Carlo Simulation and Reconstruction: Assessment of Myocardial Perfusion Imaging of Tracer Dynamics With Cardiac Motion Due to Deformation and Respiration Using Gamma Camera With Continuous Acquisition

Purpose: Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera. Methods: A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K1) and wash-out (k2) rates at each gate were also estimated. Results: The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K1 and k2 were also consistent with the ground truth. Conclusion: The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-DE-207B-05: Measuring Spatial Resolution in Digital Breast Tomosynthesis: Update of AAPM Task Group 245

Purpose: Spatial resolution in digital breast tomosynthesis (DBT) is affected by inherent/binned detector resolution, oblique entry of x-rays, and focal spot size/motion; the limited angular range further limits spatial resolution in the depth-direction. While DBT is being widely adopted clinically, imaging performance metrics and quality control protocols have not been standardized. AAPM Task Group 245 on Tomosynthesis Quality Control has been formed to address this deficiency. Methods: Methods of measuring spatial resolution are evaluated using two prototype quality control phantoms for DBT. Spatial resolution in the detector plane is measured in projection and reconstruction domains using edge-spread function (ESF), point-spread function (PSF) and modulation transfer function (MTF). Spatial resolution in the depth-direction and effective slice thickness are measured in the reconstruction domain using slice sensitivity profile (SSP) and artifact spread function (ASF). An oversampled PSF in the depth-direction is measured using a 50 µm angulated tungsten wire, from which the MTF is computed. Object-dependent PSF is derived and compared with ASF. Sensitivity of these measurements to phantom positioning, imaging conditions and reconstruction algorithms is evaluated. Results are compared from systems of varying acquisition geometry (9–25 projections over 15–60°). Dependence of measurements on feature size is investigated. Results: Measurements of spatial resolution using PSF and LSF are shown to depend on feature size; depth-direction spatial resolution measurements are shown to similarly depend on feature size for ASF, though deconvolution with an object function removes feature size-dependence. A slanted wire may be used to measure oversampled PSFs, from which MTFs may be computed for both in-plane and depth-direction resolution. Conclusion: Spatial resolution measured using PSF is object-independent with sufficiently small object; MTF is object-independent. Depth-direction spatial resolution may be measured directly using MTF or indirectly using ASF or SSP as surrogate measurements. While MTF is object-independent, it is invalid for nonlinear reconstructions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-47: Development of a High-Precision, Image-Guided Radiotherapy, Multi- Purpose Radiation Isocenter Quality-Assurance Calibration and Checking System

Purpose: To develop a system that can define the radiation isocenter and correlate this information with couch coordinates, laser alignment, optical distance indicator (ODI) settings, optical tracking system (OTS) calibrations, and mechanical isocenter walkout. Methods: Our team developed a multi-adapter, multi-purpose quality assurance (QA) and calibration device that uses an electronic portal imaging device (EPID) and in-house image-processing software to define the radiation isocenter, thereby allowing linear accelerator (Linac) components to be verified and calibrated. Motivated by the concept that each Linac component related to patient setup for image-guided radiotherapy based on cone-beam CT should be calibrated with respect to the radiation isocenter, we designed multiple concentric adapters of various materials and shapes to meet the needs of MV and KV radiation isocenter definition, laser alignment, and OTS calibration. The phantom's ability to accurately define the radiation isocenter was validated on 4 Elekta Linacs using a commercial ball bearing (BB) phantom as a reference. Radiation isocenter walkout and the accuracy of couch coordinates, ODI, and OTS were then quantified with the device. Results: The device was able to define the radiation isocenter within 0.3 mm. Radiation isocenter walkout was within ±1 mm at 4 cardinal angles. By switching adapters, we identified that the accuracy of the couch position digital readout, ODI, OTS, and mechanical isocenter walkout was within sub-mm. Conclusion: This multi-adapter, multi-purpose isocenter phantom can be used to accurately define the radiation isocenter and represents a potential paradigm shift in Linac QA. Moreover, multiple concentric adapters allowed for sub-mm accuracy for the other relevant components. This intuitive and user-friendly design is currently patent pending.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2013

Comparative performance of modern digital mammography systems in a large breast screening program

Purpose: To compare physical measures pertaining to image quality among digital mammography systems utilized in a large breast screening program. To examine qualitatively differences in these measures and differences in clinical cancer detection rates between CR and DR among sites within that program. Methods: As part of the routine quality assurance program for screening, field measurements are made of several variables considered to correlate with the diagnostic quality of medical images including: modulation transfer function, noise equivalent quanta, d′ (an index of lesion detectability) and air kerma to allow estimation of mean glandular dose. In addition, images of the mammography accreditation phantom are evaluated. Results: It was found that overall there were marked differences between the performance measures of DR and CR mammography systems. In particular, the modulation transfer functions obtained with the DR systems were found to be higher, even for larger detector element sizes. Similarly, the noise equivalent quanta, d′, and the phantom scores were higher, while the failure rates associated with low signal-to-noise ratio and high dose were lower with DR. These results were consistent with previous findings in the authors’ program that the breast cancer detection rates at sites employing CR technology were, on average, 30.6% lower than those that used DR mammography. Conclusions: While the clinical study was not large enough to allow a statistically powered system-by-system assessment of cancer detection accuracy, the physical measures expressing spatial resolution, and signal-to-noise ratio are consistent with the published finding that sites employing CR systems had lower cancer detection rates than those using DR systems for screening mammography.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-F-18C-06: Prospective Patient Evaluation of Iterative Reconstruction in Radiation Oncology

Purpose: This work incorporates iterative reconstruction (IR) into a dose reduction study to characterize image quality metrics, delineation, and dosimetric assessment, with the goal of reducing imaging dose in Radiation Oncology. Methods: Three-dimensional noise power spectrum (NPS) analysis characterized noise magnitude/texture (120 kVp, 50–200 mAs, IR levels 1–6 yielding noise reduction of 0.89–0.55 compared to filtered backprojection (FBP)). Task-specific Modulation Transfer Functions (MTFtask) were characterized across varied subject contrasts. A prospective dose reduction study (500 to 150 mAs) was conducted for 12 patients (43 inter-fraction CTs) for high-dose rate brachytherapy. Three physicians performed qualitative image assessment between full-dose FBP (FD-FBP, 500 mAs), low-dose FBP (LD-FBP, 150–250 mAs), and low-dose IRL5-6 (LD-IR) scans for image noise, cuff/bladder interface detectability, spatial resolution, texture, and segmentation confidence. Comparisons between LD-FBP and LD-IR were conducted for the following metrics: delineation (bladder and rectum evaluated via overlap indices (OI) and Dice similarity coefficients (DSC)), noise, boundary changes, dose calculation, and digitally reconstructed radiographs (DRRs). Results: NPS showed ∼50% reduction in noise magnitude and ∼0.1 1/mm spatial frequency shift with IRL6. The largest MTFtask decrease between FBP and IR was 0.08 A.U. Qualitative patient image evaluation revealed LD-IR was equivalent or slightly worse than FD-FBP, and superior to LD-FBP for all metrics except low contrast interface and texture. The largest CT number discrepancy from FBP occurred at a bone/tissue interface using IRL6 (−1.2 ± 4.9 HU (range: −17.6 – 12.5 HU)). No significant contour differences (OIs and DSCs = 0.85 – 0.95) and dose calculation discrepancy (<0.02%) were observed. DRRs preserved anatomical detail and demonstrated <2% difference in intensity between LD-FBP and LD-IRL6. Conclusion: While phantom analysis showed slight noise texture differences with IR, patient results revealed that image quality, contouring ability, and dosimetric parameters were not adversely affected, thus support integrating IR into treatment planning. Research supported in part by a grant from Philips HealthCare.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

WE-EF-207-03: Design and Optimization of a CBCT Head Scanner for Detection of Acute Intracranial Hemorrhage

Purpose: To design a dedicated x-ray cone-beam CT (CBCT) system suitable to deployment at the point-of-care and offering reliable detection of acute intracranial hemorrhage (ICH), traumatic brain injury (TBI), stroke, and other head and neck injuries. Methods: A comprehensive task-based image quality model was developed to guide system design and optimization of a prototype head scanner suitable to imaging of acute TBI and ICH. Previously reported models were expanded to include the effects of x-ray scatter correction necessary for detection of low contrast ICH and the contribution of bit depth (digitization noise) to imaging performance. Task-based detectablity index provided the objective function for optimization of system geometry, x-ray source, detector type, anti-scatter grid, and technique at 10–25 mGy dose. Optimal characteristics were experimentally validated using a custom head phantom with 50 HU contrast ICH inserts imaged on a CBCT imaging bench allowing variation of system geometry, focal spot size, detector, grid selection, and x-ray technique. Results: The model guided selection of system geometry with a nominal source-detector distance 1100 mm and optimal magnification of 1.50. Focal spot size ∼0.6 mm was sufficient for spatial resolution requirements in ICH detection. Imaging at 90 kVp yielded the best tradeoff between noise and contrast. The model provided quantitation of tradeoffs between flat-panel and CMOS detectors with respect to electronic noise, field of view, and readout speed required for imaging of ICH. An anti-scatter grid was shown to provide modest benefit in conjunction with post-acquisition scatter correction. Images of the head phantom demonstrate visualization of millimeter-scale simulated ICH. Conclusions: Performance consistent with acute TBI and ICH detection is feasible with model-based system design and robust artifact correction in a dedicated head CBCT system. Further improvements can be achieved with incorporation of model-based iterative reconstruction techniques also within the scope of the task-based optimization framework. David Foos and Xiaohui Wang are employees of Carestream Health.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-FG-207B-05: Iterative Reconstruction Via Prior Image Constrained Total Generalized Variation for Spectral CT

Purpose: To investigate iterative reconstruction via prior image constrained total generalized variation (PICTGV) for spectral computed tomography (CT) using fewer projections while achieving greater image quality. Methods: The proposed PICTGV method is formulated as an optimization problem, which balances the data fidelity and prior image constrained total generalized variation of reconstructed images in one framework. The PICTGV method is based on structure correlations among images in the energy domain and high-quality images to guide the reconstruction of energy-specific images. In PICTGV method, the high-quality image is reconstructed from all detector-collected X-ray signals and is referred as the broad-spectrum image. Distinct from the existing reconstruction methods applied on the images with first order derivative, the higher order derivative of the images is incorporated into the PICTGV method. An alternating optimization algorithm is used to minimize the PICTGV objective function. We evaluate the performance of PICTGV on noise and artifacts suppressing using phantom studies and compare the method with the conventional filtered back-projection method as well as TGV based method without prior image. Results: On the digital phantom, the proposed method outperforms the existing TGV method in terms of the noise reduction, artifacts suppression, and edge detail preservation. Compared to that obtained by the TGV based method without prior image, the relative root mean square error in the images reconstructed by the proposed method is reduced by over 20%. Conclusion: The authors propose an iterative reconstruction via prior image constrained total generalize variation for spectral CT. Also, we have developed an alternating optimization algorithm and numerically demonstrated the merits of our approach. Results show that the proposed PICTGV method outperforms the TGV method for spectral CT.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1977

Computerized tomography using video recorded fluoroscopic images

The use of video-recorded fluoroscopic images as input data for digital reconstruction of objects from their projections is examined. The fluoroscopic and the scanning apparatus used for the experiments are of a commercial type already in existence in most hospitals. It is shown that for beams with divergence up to about 15 deg, one can use a convolution algorithm designed for the parallel radiation case with negligible degradation both quantitatively and from a visual quality standpoint. This convolution algorithm is computationally more efficient than either the algebraic techniques or the convolution algorithms for radially diverging data. Results from studies on Lucite phantoms and a freshly sacrificed rat are included.

Kak, A. C.↗
DOE OSTI2014

Characterization of a commercial hybrid iterative and model-based reconstruction algorithm in radiation oncology

Purpose: Iterative reconstruction (IR) reduces noise, thereby allowing dose reduction in computed tomography (CT) while maintaining comparable image quality to filtered back-projection (FBP). This study sought to characterize image quality metrics, delineation, dosimetric assessment, and other aspects necessary to integrate IR into treatment planning. Methods: CT images (Brilliance Big Bore v3.6, Philips Healthcare) were acquired of several phantoms using 120 kVp and 25–800 mAs. IR was applied at levels corresponding to noise reduction of 0.89–0.55 with respect to FBP. Noise power spectrum (NPS) analysis was used to characterize noise magnitude and texture. CT to electron density (CT-ED) curves were generated over all IR levels. Uniformity as well as spatial and low contrast resolution were quantified using a CATPHAN phantom. Task specific modulation transfer functions (MTF{sub task}) were developed to characterize spatial frequency across objects of varied contrast. A prospective dose reduction study was conducted for 14 patients undergoing interfraction CT scans for high-dose rate brachytherapy. Three physicians performed image quality assessment using a six-point grading scale between the normal-dose FBP (reference), low-dose FBP, and low-dose IR scans for the following metrics: image noise, detectability of the vaginal cuff/bladder interface, spatial resolution, texture, segmentation confidence, and overall image quality. Contouring differences between FBP and IR were quantified for the bladder and rectum via overlap indices (OI) and Dice similarity coefficients (DSC). Line profile and region of interest analyses quantified noise and boundary changes. For two subjects, the impact of IR on external beam dose calculation was assessed via gamma analysis and changes in digitally reconstructed radiographs (DRRs) were quantified. Results: NPS showed large reduction in noise magnitude (50%), and a slight spatial frequency shift (∼0.1 mm{sup −1}) with application of IR at L6. No appreciable changes were observed for CT-ED curves between FBP and IR levels [maximum difference ∼13 HU for bone (∼1% difference)]. For uniformity, differences were ∼1 HU between FBP and IR. Spatial resolution was well conserved; the largest MTF{sub task} decrease between FBP and IR levels was 0.08 A.U. No notable changes in low-contrast detectability were observed and CNR increased substantially with IR. For the patient study, qualitative image grading showed low-dose IR was equivalent to or slightly worse than normal dose FBP, and is superior to low-dose FBP (p < 0.001 for noise), although these did not translate to differences in CT number, contouring ability, or dose calculation. The largest CT number discrepancy from FBP occurred at a bone/tissue interface using the most aggressive IR level [−1.2 ± 4.9 HU (range: −17.6–12.5 HU)]. No clinically significant contour differences were found between IR and FBP, with OIs and DSCs ranging from 0.85 to 0.95. Negligible changes in dose calculation were observed. DRRs preserved anatomical detail with <2% difference in intensity from FBP combined with aggressive IRL6. Conclusions: These results support integrating IR into treatment planning. While slight degradation in edges and shift in texture were observed in phantom, patient results show qualitative image grading, contouring ability, and dosimetric parameters were not adversely affected.

60 APPLIED LIFE SCIENCES↗