Search NASA⌕ Search

SEARCH · Search NASA

Results for “modulation transfer function”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

A forward bias method for lag correction of an a-Si flat panel detector

Purpose: Digital a-Si flat panel (FP) x-ray detectors can exhibit detector lag, or residual signal, of several percent that can cause ghosting in projection images or severe shading artifacts, known as the radar artifact, in cone-beam computed tomography (CBCT) reconstructions. A major contributor to detector lag is believed to be defect states, or traps, in the a-Si layer of the FP. Software methods to characterize and correct for the detector lag exist, but they may make assumptions such as system linearity and time invariance, which may not be true. The purpose of this work is to investigate a new hardware based method to reduce lag in an a-Si FP and to evaluate its effectiveness at removing shading artifacts in CBCT reconstructions. The feasibility of a novel, partially hardware based solution is also examined. Methods: The proposed hardware solution for lag reduction requires only a minor change to the FP. For pulsed irradiation, the proposed method inserts a new operation step between the readout and data collection stages. During this new stage the photodiode is operated in a forward bias mode, which fills the defect states with charge. A Varian 4030CB panel was modified to allow for operation in the forward bias mode. The contrast of residual lag ghosts was measured for lag frames 2 and 100 after irradiation ceased for standard and forward bias modes. Detector step response, lag, SNR, modulation transfer function (MTF), and detective quantum efficiency (DQE) measurements were made with standard and forward bias firmware. CBCT data of pelvic and head phantoms were also collected. Results: Overall, the 2nd and 100th detector lag frame residual signals were reduced 70%-88% using the new method. SNR, MTF, and DQE measurements show a small decrease in collected signal and a small increase in noise. The forward bias hardware successfully reduced the radar artifact in the CBCT reconstruction of the pelvic and head phantoms by 48%-81%. Conclusions: Overall, the forward bias method has been found to greatly reduce detector lag ghosts in projection data and the radar artifact in CBCT reconstructions. The method is limited to improvements of the a-Si photodiode response only. A future hybrid mode may overcome any limitations of this method.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

Physical evaluation of a needle photostimulable phosphor based CR mammography system

Purpose: Needle phosphor based computed radiography (CR) systems promise improved image quality compared to powder phosphor based CR units for x-ray screening mammography. This paper compares the imaging performance of needle CR cassettes, powder based CR cassettes and a well established amorphous selenium (a-Se) based flat panel based mammography system, using consistent beam qualities. Methods: Detector performance was assessed using modulation transfer function (MTF), normalized noise power spectrum (NNPS), and detective quantum efficiency (DQE). Mammography system performance was assessed against levels from the European Guidelines, including threshold gold thickness (c-d), relative signal difference to noise (SdNR) and mean glandular dose, for automatic exposure control settings suggested by the manufacturers. The needle based Agfa HM5.0 CR detector was compared against the single sided readout Agfa MM3.0R and dual sided readout Fuji Profect CS powder CR plates using a 28 kV Mo/Rh spectrum, while a 28 kV W/Rh spectrum was used to compare the Agfa HM5.0 against the Siemens MAMMOMAT Inspiration a-Se based system. Results: MTF at 5 mm{sup -1} was 0.16 and 0.24 for the needle CR detector in the fast and slow scan directions, respectively, indicating a slight improvement ({approx}20%) over the two powder CR systems but remained 50% lower than the result at 5 mm{sup -1} for the a-Se detector ({approx}0.55). Structured screen noise was lower for the needle phosphor compared to the powder plates. CR system gain, estimated from the measured absorption fraction and NNPS results, was 6.3 for the (single sided) needle phosphor and 5.1 and 7.2 for the single sided and dual sided powder phosphor systems. Peak DQE at {approx}100 {mu}Gy was 0.47 for the needle system compared to peak DQE figures of 0.33 and 0.46 for the single sided readout powder plates and dual sided readout plates. The high frequency DQE (at 5 mm{sup -1}) was 0.19 for the needle CR plates, a factor of approximately 3 greater than for the powder CR plates. At 28 kV W/Rh, 2 mm Al, peak DQE for the needle CR system was 0.45 against a value of 0.50 for the a-Se detector. The needle CR detector reached the Acceptable limit for 0.1 mm details in the European Guidelines at a mean glandular dose (MGD) of approximately 1.31 mGy imaged at 28 kV Mo/Rh, compared to figures of 2.19 and 1.43 mGy for the single sided and dual sided readout powder CR systems. The a-Se detector could reach the limit at 0.65 mGy using a 28 kV W/Rh spectrum, while the needle CR system required 1.09 mGy for the same spectrum. Conclusions: Imaging performance for the needle CR phosphor technology, characterized using MTF and DQE and threshold gold thickness demonstrated a clear improvement compared to both single and dual sided reading powder phosphor based CR systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

Monte Carlo simulation of a novel water-equivalent electronic portal imaging device using plastic scintillating fibers

Purpose: Most electronic portal imaging devices (EPIDs) developed so far use a thin Cu plate/phosphor screen to convert x-ray energies into light photons, while maintaining a high spatial resolution. This results in a low x-ray absorption and thus a low quantum efficiency (QE) of approximately 2-4% for megavoltage (MV) x-rays. A significant increase of QE is desirable for applications such as MV cone-beam computed tomography (MV-CBCT). Furthermore, the Cu plate/phosphor screen contains high atomic number (high-Z) materials, resulting in an undesirable over-response to low energy x-rays (due to photoelectric effect) as well as high energy x-rays (due to pair production) when used for dosimetric verification. Our goal is to develop a new MV x-ray detector that has a high QE and uses low-Z materials to overcome the obstacles faced by current MV x-ray imaging technologies. Methods: A new high QE and low-Z EPID is proposed. It consists of a matrix of plastic scintillating fibers embedded in a water-equivalent medium and coupled to an optically sensitive 2D active matrix flat panel imager (AMFPI) for image readout. It differs from the previous approach that uses segmented crystalline scintillators made of higher density and higher atomic number materials to detect MV x-rays. The plastic scintillating fibers are focused toward the x-ray source to avoid image blurring due to oblique incidence of off-axis x-rays. When MV x-rays interact with the scintillating fibers in the detector, scintillation light will be produced. The light photons produced in a fiber core and emitted within the acceptance angle of the fiber will be guided toward the AMFPI by total internal reflection. A Monte Carlo simulation has been used to investigate imaging and dosimetric characteristics of the proposed detector under irradiation of MV x-rays. Results: Properties, such as detection efficiency, modulation transfer function, detective quantum efficiency (DQE), energy dependence of detector response, and water-equivalence of dose response have been investigated. It has been found that the zero frequency DQE of the proposed detector can be up to 37% at 6 MV. The detector, also, is water-equivalent with a relatively uniform response to different energy x-rays as compared to current EPIDs. Conclusions: The results of our simulations show that, using plastic scintillating fibers, it is possible to construct a water-equivalent EPID that has a better energy response and a higher detection efficiency than current flat panel based EPIDs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

High resolution stationary digital breast tomosynthesis using distributed carbon nanotube x-ray source array

Purpose: The purpose of this study is to investigate the feasibility of increasing the system spatial resolution and scanning speed of Hologic Selenia Dimensions digital breast tomosynthesis (DBT) scanner by replacing the rotating mammography x-ray tube with a specially designed carbon nanotube (CNT) x-ray source array, which generates all the projection images needed for tomosynthesis reconstruction by electronically activating individual x-ray sources without any mechanical motion. The stationary digital breast tomosynthesis (s-DBT) design aims to (i) increase the system spatial resolution by eliminating image blurring due to x-ray tube motion and (ii) reduce the scanning time. Low spatial resolution and long scanning time are the two main technical limitations of current DBT technology. Methods: A CNT x-ray source array was designed and evaluated against a set of targeted system performance parameters. Simulations were performed to determine the maximum anode heat load at the desired focal spot size and to design the electron focusing optics. Field emission current from CNT cathode was measured for an extended period of time to determine the stable life time of CNT cathode for an expected clinical operation scenario. The source array was manufactured, tested, and integrated with a Selenia scanner. An electronic control unit was developed to interface the source array with the detection system and to scan and regulate x-ray beams. The performance of the s-DBT system was evaluated using physical phantoms. Results: The spatially distributed CNT x-ray source array comprised 31 individually addressable x-ray sources covering a 30 angular span with 1 pitch and an isotropic focal spot size of 0.6 mm at full width at half-maximum. Stable operation at 28 kV(peak) anode voltage and 38 mA tube current was demonstrated with extended lifetime and good source-to-source consistency. For the standard imaging protocol of 15 views over 14, 100 mAs dose, and 2 x 2 detector binning, the projection resolution along the scanning direction increased from 4.0 cycles/mm [at 10% modulation-transfer-function (MTF)] in DBT to 5.1 cycles/mm in s-DBT at magnification factor of 1.08. The improvement is more pronounced for faster scanning speeds, wider angular coverage, and smaller detector pixel sizes. The scanning speed depends on the detector, the number of views, and the imaging dose. With 240 ms detector readout time, the s-DBT system scanning time is 6.3 s for a 15-view, 100 mAs scan regardless of the angular coverage. The scanning speed can be reduced to less than 4 s when detectors become faster. Initial phantom studies showed good quality reconstructed images. Conclusions: A prototype s-DBT scanner has been developed and evaluated by retrofitting the Selenia rotating gantry DBT scanner with a spatially distributed CNT x-ray source array. Preliminary results show that it improves system spatial resolution substantially by eliminating image blur due to x-ray focal spot motion. The scanner speed of s-DBT system is independent of angular coverage and can be increased with faster detector without image degration. The accelerated lifetime measurement demonstrated the long term stability of CNT x-ray source array with typical clinical operation lifetime over 3 years.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

The relationship between anatomic noise and volumetric breast density for digital mammography

Purpose: The appearance of parenchymal/stromal patterns in mammography have been characterized as having a Wiener power spectrum with an inverse power-law shape described by the exponential parameter, {beta}. The amount of fibroglandular tissue, which can be quantified in terms of volumetric breast density (VBD), influences the texture and appearance of the patterns formed in a mammogram. Here, a large study is performed to investigate the variations in {beta} in a clinical population and to indicate the relationship between {beta} and breast density. Methods: From a set of 2686 cranio-caudal normal screening mammograms, the parameter {beta} was extracted from log-log fits to the Wiener spectrum over the range 0.15-1 mm{sup -1}. The Wiener spectrum was calculated from regions of interest in the compression paddle contact region of the breast. An in-house computer program, Cumulus V, was used to extract the volumetric breast density and identify the compression paddle contact regions of the breast. The Wiener spectra were calculated with and without modulation transfer function (MTF) correction to determine the impact of VBD on the intrinsic anatomic noise. Results: The mean volumetric breast density was 25.5% ({+-}12.6%) over all images. The mean {beta} following a MTF correction which decreased the {beta} slightly ( Almost-Equal-To -0.08) was found to be 2.87. Varying the maximum of the spatial frequency range of the fits from 0.7 to 1.0, 1.25 or 1.5 mm{sup -1} showing small decreases in the result, although the effect of the quantum noise power component on reducing {beta} was clearly observed at 1.5 mm{sup -1}. Conclusions: The texture parameter, {beta}, was found to increase with VBD at low volumetric breast densities with an apparent leveling off at higher densities. The relationship between {beta} and VBD measured here can be used to create probabilistic models for computer simulations of detectability. As breast density is a known risk predictor for breast cancer, the correlation between {beta} and VBD suggests that {beta} may provide predictive information and this will be investigated in the future.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Ideal-observer detectability in photon-counting differential phase-contrast imaging using a linear-systems approach

Purpose: To provide a cascaded-systems framework based on the noise-power spectrum (NPS), modulation transfer function (MTF), and noise-equivalent number of quanta (NEQ) for quantitative evaluation of differential phase-contrast imaging (Talbot interferometry) in relation to conventional absorption contrast under equal-dose, equal-geometry, and, to some extent, equal-photon-economy constraints. The focus is a geometry for photon-counting mammography. Methods: Phase-contrast imaging is a promising technology that may emerge as an alternative or adjunct to conventional absorption contrast. In particular, phase contrast may increase the signal-difference-to-noise ratio compared to absorption contrast because the difference in phase shift between soft-tissue structures is often substantially larger than the absorption difference. We have developed a comprehensive cascaded-systems framework to investigate Talbot interferometry, which is a technique for differential phase-contrast imaging. Analytical expressions for the MTF and NPS were derived to calculate the NEQ and a task-specific ideal-observer detectability index under assumptions of linearity and shift invariance. Talbot interferometry was compared to absorption contrast at equal dose, and using either a plane wave or a spherical wave in a conceivable mammography geometry. The impact of source size and spectrum bandwidth was included in the framework, and the trade-off with photon economy was investigated in some detail. Wave-propagation simulations were used to verify the analytical expressions and to generate example images. Results: Talbot interferometry inherently detects the differential of the phase, which led to a maximum in NEQ at high spatial frequencies, whereas the absorption-contrast NEQ decreased monotonically with frequency. Further, phase contrast detects differences in density rather than atomic number, and the optimal imaging energy was found to be a factor of 1.7 higher than for absorption contrast. Talbot interferometry with a plane wave increased detectability for 0.1-mm tumor and glandular structures by a factor of 3-4 at equal dose, whereas absorption contrast was the preferred method for structures larger than {approx}0.5 mm. Microcalcifications are small, but differ from soft tissue in atomic number more than density, which is favored by absorption contrast, and Talbot interferometry was barely beneficial at all within the resolution limit of the system. Further, Talbot interferometry favored detection of 'sharp' as opposed to 'smooth' structures, and discrimination tasks by about 50% compared to detection tasks. The technique was relatively insensitive to spectrum bandwidth, whereas the projected source size was more important. If equal photon economy was added as a restriction, phase-contrast efficiency was reduced so that the benefit for detection tasks almost vanished compared to absorption contrast, but discrimination tasks were still improved close to a factor of 2 at the resolution limit. Conclusions: Cascaded-systems analysis enables comprehensive and intuitive evaluation of phase-contrast efficiency in relation to absorption contrast under requirements of equal dose, equal geometry, and equal photon economy. The benefit of Talbot interferometry was highly dependent on task, in particular detection versus discrimination tasks, and target size, shape, and material. Requiring equal photon economy weakened the benefit of Talbot interferometry in mammography.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

A robust geometry estimation method for spiral, sequential and circular cone-beam micro-CT

Purpose: The authors propose a novel method for misalignment estimation of micro-CT scanners using an adaptive genetic algorithm. Methods: The proposed algorithm is able to estimate the rotational geometry, the direction vector of table movement and the displacement between different imaging threads of a dual source or even multisource scanner. The calibration procedure does not rely on dedicated calibration phantoms and a sequence scan of a single metal bead is sufficient to geometrically calibrate the whole imaging system for spiral, sequential, and circular scan protocols. Dual source spiral and sequential scan protocols in micro-computed tomography result in projection data that-besides the source and detector positions and orientations-also require a precise knowledge of the table direction vector to be reconstructed properly. If those geometric parameters are not known accurately severe artifacts and a loss in spatial resolution appear in the reconstructed images as long as no geometry calibration is performed. The table direction vector is further required to ensure that consecutive volumes of a sequence scan can be stitched together and to allow the reconstruction of spiral data at all. Results: The algorithm's performance is evaluated using simulations of a micro-CT system with known geometry and misalignment. To assess the quality of the algorithm in a real world scenario the calibration of a micro-CT scanner is performed and several reconstructions with and without geometry estimation are presented. Conclusions: The results indicate that the algorithm successfully estimates all geometry parameters, misalignment artifacts in the reconstructed volumes vanish, and the spatial resolution is increased as can be shown by the evaluation of modulation transfer function measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Impact of flat panel-imager veiling glare on scatter-estimation accuracy and image quality of a commercial on-board cone-beam CT imaging system

Purpose: The purposes of this study is to measure the low frequency drop (LFD) of the modulation transfer function (MTF), associated with the long tails of the detector point spread function (PSF) of an on-board flat panel imager and study its impact on cone-beam CT (CBCT) image quality and scatter measurement accuracy. Methods: Two different experimental methods were used to characterize LFD and its associated PSF of a Varian OBI flat-panel detector system: the edge response function (ERF) method and the disk transfer function (DTF) method. PSF was estimated by fitting parametric models to these measurements for four values of the applied voltage (kVp). The resultant PSF was used to demonstrate the effect of LFD on image contrast and CT number accuracy in CBCT images reconstructed from synthetic datasets, as well as, accuracy of scatter measurements with the beam-stop method. Results: The MTFs derived from the measured ERF data revealed LFDs varying from 8% (at 60 kVp) to 10.5% (at 120 kVp), while the intensity of the long PSF tails was found to increase with increasing kVp. The veiling glare line spread functions derived from the ERF and DTF methods were in excellent agreement. Uncorrected veiling glare reduced contrast and the image intensity in CBCT reconstruction, near the phantom periphery (by 67 Hounsfield units in a 20 cm-in-diameter water phantom) and (to a smaller degree) near inhomogeneities. Use of the bow-tie filter mitigated these effects. Veiling glare also resulted in about 10%-15% overestimation of the scatter-to-primary ratio when measured with the beam-stop or beam-stop array method. Conclusions: The long tails of the detector PSF were found to have a modest dependence of beam spectrum, which is reflected on the MTF curve LFD. Our findings show that uncorrected veiling glare can affect quantitative accuracy and contrast in CBCT imaging, based on flat panel imager. In addition, it results in overestimation of the scatter-to-primary ratio, measured with the beam-stop methods.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Comparison of signal to noise ratios from spatial and frequency domain formulations of nonprewhitening model observers in digital mammography

Purpose: Image quality indices based upon model observers are promising alternatives to laborious human readings of contrast-detail images. This is especially appealing in digital mammography as limiting values for contrast thresholds determine, according to some international protocols, the acceptability of these systems in the radiological practice. The objective of the present study was to compare the signal to noise ratios (SNR) obtained with two nonprewhitening matched filter model observer approaches, one in the spatial domain and the other in the frequency domain, and with both of them worked out for disks as present in the CDMAM phantom. Methods: The analysis was performed using images acquired with the Siemens Novation and Inspiration digital mammography systems. The spatial domain formulation uses a series of high dose CDMAM images as the signal and a routine exposure of two flood images to calculate the covariance matrix. The frequency domain approach uses the mathematical description of a disk and modulation transfer function (MTF) and noise power spectrum (NPS) calculated from images. Results: For both systems most of the SNR values calculated in the frequency domain were in very good agreement with the SNR values calculated in the spatial domain. Both the formulations in the frequency domain and in the spatial domain show a linear relationship between SNR and the diameter of the CDMAM discs. Conclusions: The results suggest that both formulations of the model observer lead to very similar figures of merit. This is a step forward in the adoption of figures of merit based on NPS and MTF for the acceptance testing of mammography systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

The effect of averaging adjacent planes for artifact reduction in matrix inversion tomosynthesis

Purpose: Matrix inversion tomosynthesis (MITS) uses linear systems theory and knowledge of the imaging geometry to remove tomographic blur that is present in conventional backprojection tomosynthesis reconstructions, leaving in-plane detail rendered clearly. The use of partial-pixel interpolation during the backprojection process introduces imprecision in the MITS modeling of tomographic blur, and creates low-contrast artifacts in some MITS planes. This paper examines the use of MITS slabs, created by averaging several adjacent MITS planes, as a method for suppressing partial-pixel artifacts. Methods: Human chest tomosynthesis projection data, acquired as part of an IRB-approved pilot study, were used to generate MITS planes, three-plane MITS slabs (MITSa3), five-plane MITS slabs (MITSa5), and seven-plane MITS slabs (MITSa7). These were qualitatively examined for partial-pixel artifacts and the visibility of normal and abnormal anatomy. Additionally, small (5 mm) subtle pulmonary nodules were simulated and digitally superimposed upon human chest tomosynthesis projection images, and their visibility was qualitatively assessed in the different reconstruction techniques. Simulated images of a thin wire were used to generate modulation transfer function (MTF) and slice-sensitivity profile curves for the different MITS and MITS slab techniques, and these were examined for indications of partial-pixel artifacts and frequency response uniformity. Finally, mean-subtracted, exposure-normalized noise power spectra (ENNPS) estimates were computed and compared for MITS and MITS slab reconstructions, generated from 10 sets of tomosynthesis projection data of an acrylic slab. The simulated in-plane MTF response of each technique was also combined with the square root of the ENNPS estimate to yield stochastic signal-to-noise ratio (SNR) information about the different reconstruction techniques. Results: For scan angles of 20 Degree-Sign and 5 mm plane separation, seven MITS planes must be averaged to sufficiently remove partial-pixel artifacts. MITSa7 does appear to subtly reduce the contrast of high-frequency 'edge' information, but the removal of partial-pixel artifacts makes the appearance of low-contrast, fine-detail anatomy even more conspicuous in MITSa7 slices. MITSa7 also appears to render simulated subtle 5 mm pulmonary nodules with greater visibility than MITS alone, in both the open lung and regions overlying the mediastinum. Finally, the MITSa7 technique reduces stochastic image variance, though the in-plane stochastic SNR (for very thin objects which do not span multiple MITS planes) is only improved at spatial frequencies between 0.05 and 0.20 cycles/mm. Conclusions: The MITSa7 method is an improvement over traditional single-plane MITS for thoracic imaging and the pulmonary nodule detection task, and thus the authors plan to use the MITSa7 approach for all future MITS research at the authors' institution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Validation of columnar CsI x-ray detector responses obtained with hybridMANTIS, a CPU-GPU Monte Carlo code for coupled x-ray, electron, and optical transport

Purpose: hybridMANTIS is a Monte Carlo package for modeling indirect x-ray imagers using columnar geometry based on a hybrid concept that maximizes the utilization of available CPU and graphics processing unit processors in a workstation. Methods: The authors compare hybridMANTIS x-ray response simulations to previously published MANTIS and experimental data for four cesium iodide scintillator screens. These screens have a variety of reflective and absorptive surfaces with different thicknesses. The authors analyze hybridMANTIS results in terms of modulation transfer function and calculate the root mean square difference and Swank factors from simulated and experimental results. Results: The comparison suggests that hybridMANTIS better matches the experimental data as compared to MANTIS, especially at high spatial frequencies and for the thicker screens. hybridMANTIS simulations are much faster than MANTIS with speed-ups up to 5260. Conclusions: hybridMANTIS is a useful tool for improved description and optimization of image acquisition stages in medical imaging systems and for modeling the forward problem in iterative reconstruction algorithms.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

Effective detective quantum efficiency for two mammography systems: Measurement and comparison against established metrics

Purpose: The aim of this paper was to illustrate the value of the new metric effective detective quantum efficiency (eDQE) in relation to more established measures in the optimization process of two digital mammography systems. The following metrics were included for comparison against eDQE: detective quantum efficiency (DQE) of the detector, signal difference to noise ratio (SdNR), and detectability index (d′) calculated using a standard nonprewhitened observer with eye filter.Methods: The two systems investigated were the Siemens MAMMOMAT Inspiration and the Hologic Selenia Dimensions. The presampling modulation transfer function (MTF) required for the eDQE was measured using two geometries: a geometry containing scattered radiation and a low scatter geometry. The eDQE, SdNR, and d′ were measured for poly(methyl methacrylate) (PMMA) thicknesses of 20, 40, 60, and 70 mm, with and without the antiscatter grid and for a selection of clinically relevant target/filter (T/F) combinations. Figures of merit (FOMs) were then formed from SdNR and d′ using the mean glandular dose as the factor to express detriment. Detector DQE was measured at energies covering the range of typical clinically used spectra.Results: The MTF measured in the presence of scattered radiation showed a large drop at low spatial frequency compared to the low scatter method and led to a corresponding reduction in eDQE. The eDQE for the Siemens system at 1 mm{sup −1} ranged between 0.15 and 0.27, depending on T/F and grid setting. For the Hologic system, eDQE at 1 mm{sup −1} varied from 0.15 to 0.32, again depending on T/F and grid setting. The eDQE results for both systems showed that the grid increased the system efficiency for PMMA thicknesses of 40 mm and above but showed only small sensitivity to T/F setting. While results of the SdNR and d′ based FOMs confirmed the eDQE grid position results, they were also more specific in terms of T/F selection. For the Siemens system at 20 mm PMMA, the FOMs indicated Mo/Mo (grid out) as optimal while W/Rh (grid in) was the optimal configuration at 40, 60, and 70 mm PMMA. For the Hologic, the FOMs pointed to W/Rh (grid in) at 20 and 40 mm of PMMA while W/Ag (grid in) gave the highest FOM at 60 and 70 mm PMMA. Finally, DQE at 1 mm{sup −1} averaged for the four beam qualities studied was 0.44 ± 0.02 and 0.55 ± 0.03 for the Siemens and Hologic detectors, respectively, indicating only a small influence of energy on detector DQE.Conclusions: Both the DQE and eDQE data showed only a small sensitivity to T/F setting for these two systems. The eDQE showed clear preferences in terms of scatter reduction, being highest for the grid-in geometry for PMMA thicknesses of 40 mm and above. The SdNR and d′ based figures of merit, which contain additional weighting for contrast and dose, pointed to specific T/F settings for both systems.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Characterization of a clinical unit for digital radiography based on irradiation side sampling technology

Purpose: A characterization of a clinical unit for digital radiography (FUJIFILM FDR D-EVO) is presented. This system is based on the irradiation side sampling (ISS) technology and can be equipped with two different scintillators: one traditional gadolinium-oxysulphide phosphor (GOS) and a needle structured cesium iodide (CsI) phosphor panel.Methods: The characterization was achieved in terms of response curve, modulation transfer function (MTF), noise power spectra (NPS), detective quantum efficiency (DQE), and psychophysical parameters (contrast-detail analysis with an automatic reading of CDRAD images). For both scintillation screens the authors accomplished the measurements with four standard beam conditions: RAQ3, RQA5, RQA7, and RQA9.Results: At the Nyquist frequency (3.33 lp/mm) the MTF is about 35% and 25% for CsI and GOS detectors, respectively. The CsI scintillator has better noise properties than the GOS screen in almost all the conditions. This is particularly true for low-energy beams, where the noise for the GOS system can go up to a factor 2 greater than that found for CsI. The DQE of the CsI detector reaches a peak of 60%, 60%, 58%, and 50% for the RQA3, RQA5, RQA7, and RQA9 beams, respectively, whereas for the GOS screen the maximum DQE is 40%, 44%, 44%, and 35%. The contrast-detail analysis confirms that in the majority of cases the CsI scintillator is able to provide improved outcomes to those obtained with the GOS screen.Conclusions: The limited diffusion of light produced by the ISS reading makes possible the achievement of very good spatial resolution. In fact, the MTF of the unit with the CsI panel is only slightly lower to that achieved with direct conversion detectors. The combination of very good spatial resolution, together with the good noise properties reached with the CsI screen, allows achieving DQE on average about 1.5 times greater than that obtained with GOS. In fact, the DQE of unit equipped with CsI is comparable to the best alternative methods available which are based on the same technology, and similar to others based on an a-Se direct conversion detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

A quality assurance framework for the fully automated and objective evaluation of image quality in cone-beam computed tomography

Purpose: Thousands of cone-beam computed tomography (CBCT) scanners for vascular, maxillofacial, neurological, and body imaging are in clinical use today, but there is no consensus on uniform acceptance and constancy testing for image quality (IQ) and dose yet. The authors developed a quality assurance (QA) framework for fully automated and time-efficient performance evaluation of these systems. In addition, the dependence of objective Fourier-based IQ metrics on direction and position in 3D volumes was investigated for CBCT. Methods: The authors designed a dedicated QA phantom 10 cm in length consisting of five compartments, each with a diameter of 10 cm, and an optional extension ring 16 cm in diameter. A homogeneous section of water-equivalent material allows measuring CT value accuracy, image noise and uniformity, and multidimensional global and local noise power spectra (NPS). For the quantitative determination of 3D high-contrast spatial resolution, the modulation transfer function (MTF) of centrally and peripherally positioned aluminum spheres was computed from edge profiles. Additional in-plane and axial resolution patterns were used to assess resolution qualitatively. The characterization of low-contrast detectability as well as CT value linearity and artifact behavior was tested by utilizing sections with soft-tissue-equivalent and metallic inserts. For an automated QA procedure, a phantom detection algorithm was implemented. All tests used in the dedicated QA program were initially verified in simulation studies and experimentally confirmed on a clinical dental CBCT system. Results: The automated IQ evaluation of volume data sets of the dental CBCT system was achieved with the proposed phantom requiring only one scan for the determination of all desired parameters. Typically, less than 5 min were needed for phantom set-up, scanning, and data analysis. Quantitative evaluation of system performance over time by comparison to previous examinations was also verified. The maximum percentage interscan variation of repeated measurements was less than 4% and 1.7% on average for all investigated quality criteria. The NPS-based image noise differed by less than 5% from the conventional standard deviation approach and spatially selective 10% MTF values were well comparable to subjective results obtained with 3D resolution pattern. Determining only transverse spatial resolution and global noise behavior in the central field of measurement turned out to be insufficient. Conclusions: The proposed framework transfers QA routines employed in conventional CT in an advanced version to CBCT for fully automated and time-efficient evaluation of technical equipment. With the modular phantom design, a routine as well as an expert version for assessing IQ is provided. The QA program can be used for arbitrary CT units to evaluate 3D imaging characteristics automatically.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Improved proton computed tomography by dual modality image reconstruction

Purpose: Proton computed tomography (CT) is a promising image modality for improving the stopping power estimates and dose calculations for particle therapy. However, the finite range of about 33 cm of water of most commercial proton therapy systems limits the sites that can be scanned from a full 360° rotation. In this paper the authors propose a method to overcome the problem using a dual modality reconstruction (DMR) combining the proton data with a cone-beam x-ray prior. Methods: A Catphan 600 phantom was scanned using a cone beam x-ray CT scanner. A digital replica of the phantom was created in the Monte Carlo code Geant4 and a 360° proton CT scan was simulated, storing the entrance and exit position and momentum vector of every proton. Proton CT images were reconstructed using a varying number of angles from the scan. The proton CT images were reconstructed using a constrained nonlinear conjugate gradient algorithm, minimizing total variation and the x-ray CT prior while remaining consistent with the proton projection data. The proton histories were reconstructed along curved cubic-spline paths. Results: The spatial resolution of the cone beam CT prior was retained for the fully sampled case and the 90° interval case, with the MTF = 0.5 (modulation transfer function) ranging from 5.22 to 5.65 linepairs/cm. In the 45° interval case, the MTF = 0.5 dropped to 3.91 linepairs/cm For the fully sampled DMR, the maximal root mean square (RMS) error was 0.006 in units of relative stopping power. For the limited angle cases the maximal RMS error was 0.18, an almost five-fold improvement over the cone beam CT estimate. Conclusions: Dual modality reconstruction yields the high spatial resolution of cone beam x-ray CT while maintaining the improved stopping power estimation of proton CT. In the case of limited angles, the use of prior image proton CT greatly improves the resolution and stopping power estimate, but does not fully achieve the quality of a 360° proton CT scan.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Free software for performing physical analysis of systems for digital radiography and mammography

Purpose: In this paper, the authors present a free software for assisting users in achieving the physical characterization of x-ray digital systems and image quality checks. Methods: The program was developed as a plugin of a well-known public-domain suite ImageJ. The software can assist users in calculating various physical parameters such as the response curve (also termed signal transfer property), modulation transfer function (MTF), noise power spectra (NPS), and detective quantum efficiency (DQE). It also includes the computation of some image quality checks: defective pixel analysis, uniformity, dark analysis, and lag. Results: The software was made available in 2009 and has been used during the last couple of years by many users who gave us valuable feedback for improving its usability. It was tested for achieving the physical characterization of several clinical systems for digital radiography and mammography. Various published papers made use of the outcomes of the plugin. Conclusions: This software is potentially beneficial to a variety of users: physicists working in hospitals, staff working in radiological departments, such as medical physicists, physicians, engineers. The plugin, together with a brief user manual, are freely available and can be found online ( http://www.medphys.it/downloads.htm ). With our plugin users can estimate all three most important parameters used for physical characterization (MTF, NPS, and also DQE). The plugin can run on any operating system equipped with ImageJ suite. The authors validated the software by comparing MTF and NPS curves on a common set of images with those obtained with other dedicated programs, achieving a very good agreement.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Task-based detectability in CT image reconstruction by filtered backprojection and penalized likelihood estimation

Purpose: Nonstationarity is an important aspect of imaging performance in CT and cone-beam CT (CBCT), especially for systems employing iterative reconstruction. This work presents a theoretical framework for both filtered-backprojection (FBP) and penalized-likelihood (PL) reconstruction that includes explicit descriptions of nonstationary noise, spatial resolution, and task-based detectability index. Potential utility of the model was demonstrated in the optimal selection of regularization parameters in PL reconstruction. Methods: Analytical models for local modulation transfer function (MTF) and noise-power spectrum (NPS) were investigated for both FBP and PL reconstruction, including explicit dependence on the object and spatial location. For FBP, a cascaded systems analysis framework was adapted to account for nonstationarity by separately calculating fluence and system gains for each ray passing through any given voxel. For PL, the point-spread function and covariance were derived using the implicit function theorem and first-order Taylor expansion according toFessler [“Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): Applications to tomography,” IEEE Trans. Image Process. 5(3), 493–506 (1996)]. Detectability index was calculated for a variety of simple tasks. The model for PL was used in selecting the regularization strength parameter to optimize task-based performance, with both a constant and a spatially varying regularization map. Results: Theoretical models of FBP and PL were validated in 2D simulated fan-beam data and found to yield accurate predictions of local MTF and NPS as a function of the object and the spatial location. The NPS for both FBP and PL exhibit similar anisotropic nature depending on the pathlength (and therefore, the object and spatial location within the object) traversed by each ray, with the PL NPS experiencing greater smoothing along directions with higher noise. The MTF of FBP is isotropic and independent of location to a first order approximation, whereas the MTF of PL is anisotropic in a manner complementary to the NPS. Task-based detectability demonstrates dependence on the task, object, spatial location, and smoothing parameters. A spatially varying regularization “map” designed from locally optimal regularization can improve overall detectability beyond that achievable with the commonly used constant regularization parameter. Conclusions: Analytical models for task-based FBP and PL reconstruction are predictive of nonstationary noise and resolution characteristics, providing a valuable framework for understanding and optimizing system performance in CT and CBCT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
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↗