Search

Background: A common rule of thumb for object detection is the Rose criterion, which states that a signal must be five standard deviations above background to be detectable to a human observer. The validity of the Rose criterion in CT imaging is limited due to the presence of correlated noise. Recent reconstruction and denoising methodologies are also able to restore apparent image quality in very noisy conditions, and the ultimate limits of these methodologies are not yet known. Purpose: To establish a lower bound on the minimum achievable signal‐to‐noise ratio (SNR) for object detection, below which detection performance is poor regardless of reconstruction or denoising methodology. Methods: We consider a numerical observer that operates on projection data and has perfect knowledge of the background and the objects to be detected, and determine the minimum projection SNR that is necessary to achieve predetermined lesion‐level sensitivity and case‐level specificity targets. We define a set of discrete signal objects O$\mathcal{O}$ that encompasses any lesion of interest and could include lesions of different sizes, shapes, and locations. The task is to determine which object of O$\mathcal{O}$ is present, or to state the null hypothesis that no object is present. We constrain each object in O$\mathcal{O}$ to have equivalent projection SNR and use Monte Carlo methods to calculate the required projection SNR necessary. Because our calculations are performed in projection space, they impose an upper limit on the performance possible from reconstructed images. We chose O$\mathcal{O}$ to be a collection of elliptical or circular low contrast metastases and simulated detection of these objects in a parallel beam system with Gaussian statistics. Unless otherwise stated, we assume a target of 80% lesion‐level sensitivity and 80% case‐level specificity and a search field of view that is 6 cm by 6 cm by 10 slices. Results: When O$\mathcal{O}$ contains only a single object, our problem is equivalent to two‐alternative forced choice (2AFC) and the required projection SNR is 1.7. When O$\mathcal{O}$ consists of circular 6‐mm lesions at different locations in space, the required projection SNR is 5.1. When O$\mathcal{O}$ is extended to include ellipses and circles of different sizes, the required projection SNR increases to 5.3. The required SNR increases if the sensitivity target, specificity target, or search field of view increases. Conclusions: Even with perfect knowledge of the background and target objects, the ideal observer still requires an SNR of approximately 5. This is a lower bound on the SNR that would be required in real conditions, where the background and target objects are not known perfectly. Algorithms that denoise lesions with less than 5 projection SNR, regardless of the denoising methodology, are expected to show vanishing effects or false positive lesions. [ABSTRACT FROM AUTHOR]

Purpose: Eye‐tracking approaches have been used to understand the visual search process in radiology. However, previous eye‐tracking work in computer tomography (CT) has been limited largely to single cross‐sectional images or video playback of the reconstructed volume, which do not accurately reflect radiologists' visual search activities and their interactivity with three‐dimensional image data at a computer workstation (e.g., scroll, pan, and zoom) for visual evaluation of diagnostic imaging targets. We have developed a platform that integrates eye‐tracking hardware with in‐house‐developed reader workstation software to allow monitoring of the visual search process and reader‐image interactions in clinically relevant reader tasks. The purpose of this work is to validate the spatial accuracy of eye‐tracking data using this platform for different eye‐tracking data acquisition modes. Methods: An eye‐tracker was integrated with a previously developed workstation designed for reader performance studies. The integrated system captured real‐time eye movement and workstation events at 1000 Hz sampling frequency. The eye‐tracker was operated either in head‐stabilized mode or in free‐movement mode. In head‐stabilized mode, the reader positioned their head on a manufacturer‐provided chinrest. In free‐movement mode, a biofeedback tool emitted an audio cue when the head position was outside the data collection range (general biofeedback) or outside a narrower range of positions near the calibration position (strict biofeedback). Four radiologists and one resident were invited to participate in three studies to determine eye‐tracking spatial accuracy under three constraint conditions: head‐stabilized mode (i.e., with use of a chin rest), free movement with general biofeedback, and free movement with strict biofeedback. Study 1 evaluated the impact of head stabilization versus general or strict biofeedback using a cross‐hair target prior to the integration of the eye‐tracker with the image viewing workstation. In Study 2, after integration of the eye‐tracker and reader workstation, readers were asked to fixate on targets that were randomly distributed within a volumetric digital phantom. In Study 3, readers used the integrated system to scroll through volumetric patient CT angiographic images while fixating on the centerline of designated blood vessels (from the left coronary artery to dorsalis pedis artery). Spatial accuracy was quantified as the offset between the center of the intended target and the detected fixation using units of image pixels and the degree of visual angle. Results: The three head position constraint conditions yielded comparable accuracy in the studies using digital phantoms. For Study 1 involving the digital crosshairs, the median ± the standard deviation of offset values among readers were 15.2 ± 7.0 image pixels with the chinrest, 14.2 ± 3.6 image pixels with strict biofeedback, and 19.1 ± 6.5 image pixels with general biofeedback. For Study 2 using the random dot phantom, the median ± standard deviation offset values were 16.7 ± 28.8 pixels with use of a chinrest, 16.5 ± 24.6 pixels using strict biofeedback, and 18.0 ± 22.4 pixels using general biofeedback, which translated to a visual angle of about 0.8° for all three conditions. We found no obvious association between eye‐tracking accuracy and target size or view time. In Study 3 viewing patient images, use of the chinrest and strict biofeedback demonstrated comparable accuracy, while the use of general biofeedback demonstrated a slightly worse accuracy. The median ± standard deviation of offset values were 14.8 ± 11.4 pixels with use of a chinrest, 21.0 ± 16.2 pixels using strict biofeedback, and 29.7 ± 20.9 image pixels using general biofeedback. These corresponded to visual angles ranging from 0.7° to 1.3°. Conclusions: An integrated eye‐tracker system to assess reader eye movement and interactive viewing in relation to imaging targets demonstrated reasonable spatial accuracy for assessment of visual fixation. The head‐free movement condition with audio biofeedback performed similarly to head‐stabilized mode. [ABSTRACT FROM AUTHOR]

Objective: In X‐raycomputed tomography (CT), many important clinical applications may benefit from a fast acquisition speed. The helical scan is the most widely used acquisition mode in clinical CT, where a fast helical pitch can improve the acquisition speed. However, on a typical single‐source helical CT (SSCT) system, the helical pitch p typically cannot exceed 1.5; otherwise, reconstruction artifacts will result from data insufficiency. The purpose of this work is to develop a deep convolutional neural network (CNN) to correct for artifacts caused by an ultra‐fast pitch, which can enable faster acquisition speed than what is currently achievable. Methods: A customized CNN (denoted as ultra‐fast‐pitch network (UFP‐net)) was developed to restore the underlying anatomical structure from the artifact‐corrupted post‐reconstruction data acquired from SSCT with ultra‐fast pitch (i.e., p ≥ 2). UFP‐net employed residual learning to capture the features of image artifacts. UFP‐net further deployed in‐house‐customized functional blocks with spatial‐domain local operators and frequency‐domain non‐local operators, to explore multi‐scale feature representation. Images of contrast‐enhanced patient exams (n = 83) with routine pitch setting (i.e., p < 1) were retrospectively collected, which were used as training and testing datasets. This patient cohort involved CT exams over different scan ranges of anatomy (chest, abdomen, and pelvis) and CT systems (Siemens Definition, Definition Flash, Definition AS+, Siemens Healthcare, Inc.), and the corresponding base CT scanning protocols used consistent settings of major scan parameters (e.g., collimation and pitch). Forward projection of the original images was calculated to synthesize helical CT scans with one regular pitch setting (p = 1) and two ultra‐fast‐pitch setting (p = 2 and 3). All patient images were reconstructed using the standard filtered‐back‐projection (FBP) algorithm. A customized multi‐stage training scheme was developed to incrementally optimize the parameters of UFP‐net, using ultra‐fast‐pitch images as network inputs and regular pitch images as labels. Visual inspection was conducted to evaluate image quality. Structural similarity index (SSIM) and relative root‐mean‐square error (rRMSE) were used as quantitative quality metrics. Results: The UFP‐net dramatically improved image quality over standard FBP at both ultra‐fast‐pitch settings. At p = 2, UFP‐net yielded higher mean SSIM (> 0.98) with lower mean rRMSE (< 2.9%), compared to FBP (mean SSIM < 0.93; mean rRMSE > 9.1%). Quantitative metrics at p = 3: UFP‐net—mean SSIM [0.86, 0.94] and mean rRMSE [5.0%, 8.2%]; FBP—mean SSIM [0.36, 0.61] and mean rRMSE [36.0%, 58.6%]. Conclusion: The proposed UFP‐net has the potential to enable ultra‐fast data acquisition in clinical CT without sacrificing image quality. This method has demonstrated reasonable generalizability over different body parts when the corresponding CT exams involved consistent base scan parameters. [ABSTRACT FROM AUTHOR]

Purpose: Low contrast detectability (LCD) is a metric of fundamental importance in computed tomography (CT) imaging. In spite of this, its measurement is challenging in the context of nonlinear data processing. We introduce a new framework for objectively characterizing LCD with a single scan of a special‐purpose phantom and automated analysis software. The output of the analysis software is a "machine LCD" metric which is more representative of LCD than contrast‐noise ratio (CNR). It is not intended to replace human observer or model observer studies. Methods: Following preliminary simulations, we fabricated a phantom containing hundreds of low‐contrast beads. These beads are acrylic spheres (1.6 mm, net contrast ~10 HU) suspended and randomly dispersed in a background matrix of nylon pellets and isoattenuating saline. The task was to search for and localize the beads. A modified matched filter was used to automatically scan the reconstruction and select candidate bead localizations of varying confidence. These were compared to bead locations as determined from a high‐dose reference scan to produce free‐response ROC curves. We compared iterative reconstruction (IR) and filtered backpropagation (FBP) at multiple dose levels between 40 and 240 mAs. The scans at 60, 120, and 180 mAs were performed three times each to estimate uncertainty. Results: Experimental scans demonstrated the feasibility of our technique. Our metric for machine LCD was the area under the exponential transform of the FROC curve (AUC). AUC increased monotonically from 0.21 at 40 mAs to 0.84 at 240 mAs. The sample standard deviation of AUC was approximately 0.02. This measurement uncertainty in AUC corresponded to a change in tube current of 4% to 8%. Surprisingly, we found that AUCs for IR were slightly worse than AUCs for FBP. While the phantom was sufficient for these experiments, it contained small air bubbles and alternative fabrication methods will be necessary for widespread utilization. Conclusions: It is feasible to measure machine LCD using a search task on a phantom with hundreds of beads and to obtain tight error bars using only a single scan. Our method could facilitate routine quality assurance or possibly enable comparisons between different protocols and scanners. [ABSTRACT FROM AUTHOR]

Photon counting detectors (PCDs) are classically described as being either paralyzable or nonparalyzable. When the PCD is paralyzed, it is no longer sensitive to the detection of additional flux. A recent strategy in PCD design has been to compensate for detector paralysis by embedding specialized paralysis compensation electronics into the application-specific integrated circuit (ASIC). One such compensation mechanism is the pileup trigger, which places an additional energy bin at very high energy that is triggered only during pileup. Another compensation mechanism is the retrigger architecture, which converts a paralyzable PCD into a nonparalyzable PCD. We propose a third mechanism that modifies the retrigger architecture using dedicated secondary counters. We studied the incremental benefit of these three paralysis compensation mechanisms in simulation. We modeled the spectral response using Monte Carlo simulations and then estimated the variance in basis material decomposition of a single pixel using the Cramér-Rao lower bound (CRLB). In the absence of paralysis compensation, noise in basis material images shows sharp increases at moderate flux (near the characteristic count rate) due to contrast inversion and again at high flux. The pileup trigger reduces noise at high flux but does not eliminate contrast inversion. The retrigger architecture eliminates contrast inversion but does not reduce noise at high flux. Our proposed retrigger architecture with dedicated secondary counters reduce noise at both moderate and high flux. [ABSTRACT FROM AUTHOR]

The performance of X-ray photon counting detectors (PCDs), especially on spectral tasks, is compromised by charge sharing. Existing mechanisms to compensate for charge sharing, such as charge summing circuitry or larger pixel sizes, increase and aggravate pileup effects. We propose a new mechanism, the coincidence counting bin (CCB), which does not increase pileup and which has implementation similarities to existing energy bins. The CCB is triggered by coincident events in adjacent pixels and provides an estimate of the double counts arising from charge sharing. Unlike charge summing, the CCB does not directly restore corrupted events. Nonetheless, knowledge of the number of coincident counts can be used by the estimator to reduce noise. We simulated a PCD with and without the CCB using Monte Carlo simulations, modeling PCD pixels as instantaneous charge collectors and X-ray energy deposition as producing a Gaussian charge cloud with 75 micron FWHM, independent of energy. With typical operating conditions and at low flux (120 kVp, incident count rate 1% of characteristic count rate, 30 cm object thickness, five energy bins, pixel pitch of 300 microns), the CCB improved dose efficiency of iodine and water basis material decomposition by 70% and 50%, respectively. An improvement of 20% was also seen in an iodine CNR task. These improvements are attenuated as incident flux increases and show moderate dependence on filtration and pixel size. At high flux, the CCB does not provide useful information and is discarded by the estimator. The CCB may be an effective and practical mechanism for charge sharing compensation in PCDs. [ABSTRACT FROM AUTHOR]

Charge sharing, scatter, and fluorescence events in a photon counting detector can result in counting of a single incident photon in multiple neighboring pixels, each at a fraction of the true energy. This causes energy distortion and correlation of data across energy bins in neighboring pixels (spatio-energy correlation), with the severity depending on the detector pixel size and detector material. If a “macro-pixel” is formed by combining the counts from multiple adjacent small pixels, it will exhibit correlations across its energy bins. Understanding these effects can be crucial for detector design and for model-based imaging applications. This paper investigates the impact of these effects in basis material and effective monoenergetic estimates using the Cramér-Rao Lower Bound. To do so, we derive a correlation model for the multi-counting events. CdTe detectors with grids of pixels with side length of $250~\mu \text{m}$ , $500~\mu \text{m}$ , and 1 mm were compared, with binning of $4\times4$ , $2\times2$ , and $1\times1$ pixels, respectively, to keep the same net 1 mm2 aperture constant. The same flux was applied to each. The mean and covariance matrix of measured photon counts were derived analytically using spatio-energy response functions precomputed from Monte Carlo simulations. Our results show that a 1 mm2 macro-pixel with $250\times 250\,\,\mu \text{m}^{\textsf {2}}$ sub-pixels shows 35% higher standard deviation than a single 1 mm2 pixel for material-specific imaging, while the penalty for effective monoenergetic imaging is <10% compared with a single 1 mm $^{\textsf {2}}$ pixel. Potential benefits of sub-pixels (higher spatial resolution and lower pulse pile-up effects) are important but were not investigated here. [ABSTRACT FROM AUTHOR]

Purpose: To facilitate investigations into the impacts of acquisition and reconstruction parameters on quantitative imaging, radiomics and CAD using CT imaging, we previously released an open‐source implementation of a conventional weighted filtered backprojection reconstruction called FreeCT_wFBP. Our purpose was to extend that work by providing an open‐source implementation of a model‐based iterative reconstruction method using coordinate descent optimization, called FreeCT_ICD. Methods: Model‐based iterative reconstruction offers the potential for substantial radiation dose reduction, but can impose substantial computational processing and storage requirements. FreeCT_ICD is an open‐source implementation of a model‐based iterative reconstruction method that provides a reasonable tradeoff between these requirements. This was accomplished by adapting a previously proposed method that allows the system matrix to be stored with a reasonable memory requirement. The method amounts to describing the attenuation coefficient using rotating slices that follow the helical geometry. In the initially proposed version, the rotating slices are themselves described using blobs. We have replaced this description by a unique model that relies on trilinear interpolation together with the principles of Joseph's method. This model offers an improvement in memory requirement while still allowing highly accurate reconstruction for conventional CT geometries. The system matrix is stored column‐wise and combined with an iterative coordinate descent (ICD) optimization. The result is FreeCT_ICD, which is a reconstruction program developed on the Linux platform using C++ libraries and released under the open‐source GNU GPL v2.0 license. The software is capable of reconstructing raw projection data of helical CT scans. In this work, the software has been described and evaluated by reconstructing datasets exported from a clinical scanner which consisted of an ACR accreditation phantom dataset and a clinical pediatric thoracic scan. Results: For the ACR phantom, image quality was comparable to clinical reconstructions as well as reconstructions using open‐source FreeCT_wFBP software. The pediatric thoracic scan also yielded acceptable results. In addition, we did not observe any deleterious impact in image quality associated with the utilization of rotating slices. These evaluations also demonstrated reasonable tradeoffs in storage requirements and computational demands. Conclusion: FreeCT_ICD is an open‐source implementation of a model‐based iterative reconstruction method that extends the capabilities of previously released open‐source reconstruction software and provides the ability to perform vendor‐independent reconstructions of clinically acquired raw projection data. This implementation represents a reasonable tradeoff between storage and computational requirements and has demonstrated acceptable image quality in both simulated and clinical image datasets. [ABSTRACT FROM AUTHOR]

Purpose: Photon‐counting detectors using CdTe or CZT substrates are promising candidates for future CT systems but suffer from a number of nonidealities, including charge sharing and pulse pileup. By increasing the pixel size of the detector, the system can improve charge sharing characteristics at the expense of increasing pileup. The purpose of this work is to describe these considerations in the optimization of the detector pixel pitch. Methods: The transport of x rays through the CdTe substrate was simulated in a Monte Carlo fashion using GEANT4. Deposited energy was converted into charges distributed as a Gaussian function with size dependent on interaction depth to capture spreading from diffusion and Coulomb repulsion. The charges were then collected in a pixelated fashion. Pulse pileup was incorporated separately with Monte Carlo simulation. The Cramér–Rao lower bound (CRLB) of the measurement variance was numerically estimated for the basis material projections. Noise in these estimates was propagated into CT images. We simulated pixel pitches of 250, 350, and 450 microns and compared the results to a photon counting detector with pileup but otherwise ideal energy response and an ideal dual‐energy system (80/140 kVp with tin filtration). The modeled CdTe thickness was 2 mm, the incident spectrum was 140 kVp and 500 mA, and the effective dead time was 67 ns. Charge summing circuitry was not modeled. We restricted our simulations to objects of uniform thickness and did not consider the potential advantage of smaller pixels at high spatial frequencies. Results: At very high x‐ray flux, pulse pileup dominates and small pixel sizes perform best. At low flux or for thick objects, charge sharing dominates and large pixel sizes perform best. At low flux and depending on the beam hardness, the CRLB of variance in basis material projections tasks can be 32%–55% higher with a 250 micron pixel pitch compared to a 450 micron pixel pitch. However, both are about four times worse in variance than the ideal photon counting detector. The optimal pixel size depends on a number of factors such as x‐ray technique and object size. At high technique (140 kVp/500 mA), the ratio of variance for a 450 micron pixel compared to a 250 micron pixel size is 2126%, 200%, 97%, and 78% when imaging 10, 15, 20, and 25 cm of water, respectively. If 300 mg/cm2 of iodine is also added to the object, the variance ratio is 117%, 91%, 74%, and 72%, respectively. Nonspectral tasks, such as equivalent monoenergetic imaging, are less sensitive to spectral distortion. Conclusions: The detector pixel size is an important design consideration in CdTe detectors. Smaller pixels allow for improved capabilities at high flux but increase charge sharing, which in turn compromises spectral performance. The optimal pixel size will depend on the specific task and on the charge shaping time. [ABSTRACT FROM AUTHOR]

Purpose: To find an upper bound on the maximum dose reduction possible for any reconstruction algorithm, analytic or iterative, that result from the inclusion of the data statistics. The authors do not analyze noise reduction possible from prior knowledge or assumptions about the object. Methods: The authors examined the task of estimating the density of a circular lesion in a cross section. Raw data were simulated by forward projection of existing images and numerical phantoms. To assess an upper bound on the achievable dose reduction by any algorithm, the authors assume that both the background and the shape of the lesion are completely known. Under these conditions, the best possible estimate of the density can be determined by solving a weighted least squares problem directly in the raw data domain. Any possible reconstruction algorithm that does not use prior knowledge or make assumptions about the object, including filtered backprojection (FBP) or iterative reconstruction methods with this constraint, must be no better than this least squares solution. The authors simulated 10 000 sets of noisy data and compared the variance in density from the least squares solution with those from FBP. Density was estimated from FBP images using either averaging within a ROI, or streak-adaptive averaging with better noise performance. Results: The bound on the possible dose reduction depends on the degree to which the observer can read through the possibly streaky noise. For the described low contrast detection task with the signal shape and background known exactly, the average dose reduction possible compared to FBP with streak-adaptive averaging was 42% and it was 64% if only the ROI average is used with FBP. The exact amount of dose reduction also depends on the background anatomy, with statistically inhomogeneous backgrounds showing greater benefits. Conclusions: The dose reductions from new, statistical reconstruction methods can be bounded. Larger dose reductions in the density estimation task studied here are only possible with the introduction of prior knowledge, which can introduce bias. [ABSTRACT FROM AUTHOR]

Energy-discriminating, photon counting (EDPC) detectors have high potential in spectral imaging applications but exhibit degraded performance when the incident count rate approaches or exceeds the characteristic count rate of the detector. In order to reduce the requirements on the detector, we explore the strategy of modulating the X-ray flux field using a recently proposed dynamic, piecewise-linear attenuator. A previous paper studied this modulation for photon counting detectors but did not explore the impact on spectral applications. In this work, we modeled detection with a bipolar triangular pulse shape (Taguchi , 2011) and estimated the Cramer-Rao lower bound (CRLB) of the variance of material selective and equivalent monoenergetic images, assuming deterministic errors at high flux could be corrected. We compared different materials for the dynamic attenuator and found that rare earth elements, such as erbium, outperformed previously proposed materials such as iron in spectral imaging. The redistribution of flux reduces the variance or dose, consistent with previous studies on benefits with conventional detectors. Numerical simulations based on DICOM datasets were used to assess the impact of the dynamic attenuator for detectors with several different characteristic count rates. The dynamic attenuator reduced the peak incident count rate by a factor of 4 in the thorax and 44 in the pelvis, and a 10 Mcps/mm^2 EDPC detector with dynamic attenuator provided generally superior image quality to a 100 Mcps/mm^2 detector with reference bowtie filter for the same dose. The improvement is more pronounced in the material images. [ABSTRACT FROM PUBLISHER]

Purpose: Truncation artifacts in CT occur if the object to be imaged extends past the scanner field of view (SFOV). These artifacts impede diagnosis and could possibly introduce errors in dose plans for radiation therapy. Several approaches exist for correcting truncation artifacts, but existing correction algorithms do not accurately recover the skin line (or support) of the patient, which is important in some dose planning methods. The purpose of this paper was to develop an iterative algorithm that recovers the support of the object. Methods: The authors assume that the truncated portion of the image is made up of soft tissue of uniform CT number and attempt to find a shape consistent with the measured data. Each known measurement in the sinogram is interpreted as an estimate of missing mass along a line. An initial estimate of the object support is generated by thresholding a reconstruction made using a previous truncation artifact correction algorithm (e.g., water cylinder extrapolation). This object support is iteratively deformed to reduce the inconsistency with the measured data. The missing data are estimated using this object support to complete the dataset. The method was tested on simulated and experimentally truncated CT data. Results: The proposed algorithm produces a better defined skin line than water cylinder extrapolation. On the experimental data, the RMS error of the skin line is reduced by about 60%. For moderately truncated images, some soft tissue contrast is retained near the SFOV. As the extent of truncation increases, the soft tissue contrast outside the SFOV becomes unusable although the skin line remains clearly defined, and in reformatted images it varies smoothly from slice to slice as expected. Conclusions: The support recovery algorithm provides a more accurate estimate of the patient outline than thresholded, basic water cylinder extrapolation, and may be preferred in some radiation therapy applications. [ABSTRACT FROM AUTHOR]

Purpose: The authors describe algorithms to control dynamic attenuators in CT and compare their performance using simulated scans. Dynamic attenuators are prepatient beam shaping filters that modulate the distribution of x-ray fluence incident on the patient on a view-by-view basis. These attenuators can reduce dose while improving key image quality metrics such as peak or mean variance. In each view, the attenuator presents several degrees of freedom which may be individually adjusted. The total number of degrees of freedom across all views is very large, making many optimization techniques impractical. The authors develop a theory for optimally controlling these attenuators. Special attention is paid to a theoretically perfect attenuator which controls the fluence for each ray individually, but the authors also investigate and compare three other, practical attenuator designs which have been previously proposed: the piecewise-linear attenuator, the translating attenuator, and the double wedge attenuator. Methods: The authors pose and solve the optimization problems of minimizing the mean and peak variance subject to a fixed dose limit. For a perfect attenuator and mean variance minimization, this problem can be solved in simple, closed form. For other attenuator designs, the problem can be decomposed into separate problems for each view to greatly reduce the computational complexity. Peak variance minimization can be approximately solved using iterated, weighted mean variance (WMV) minimization. Also, the authors develop heuristics for the perfect and piecewise-linear attenuators which do not require a priori knowledge of the patient anatomy. The authors compare these control algorithms on different types of dynamic attenuators using simulated raw data from forward projected DICOM files of a thorax and an abdomen. Results: The translating and double wedge attenuators reduce dose by an average of 30% relative to current techniques (bowtie filter with tube current modulation) without increasing peak variance. The 15-element piecewise-linear dynamic attenuator reduces dose by an average of 42%, and the perfect attenuator reduces dose by an average of 50%. Improvements in peak variance are several times larger than improvements in mean variance. Heuristic control eliminates the need for a prescan. For the piecewise-linear attenuator, the cost of heuristic control is an increase in dose of 9%. The proposed iterated WMV minimization produces results that are within a few percent of the true solution. Conclusions: Dynamic attenuators show potential for significant dose reduction. A wide class of dynamic attenuators can be accurately controlled using the described methods. [ABSTRACT FROM AUTHOR]

Purpose: The authors recently proposed a dynamic, prepatient x-ray attenuator capable of producing a piecewise-linear attenuation profile customized to each patient and viewing angle. This attenuator was intended to reduce scatter-to-primary ratio (SPR), dynamic range, and dose by redistributing flux. In this work the authors tested the ability of the attenuator to reduce dose and SPR in simulations. Methods: The authors selected four clinical applications, including routine full field-of-view scans of the thorax and abdomen, and targeted reconstruction tasks for an abdominal aortic aneurysm and the pancreas. Raw data were estimated by forward projection of the image volume datasets. The dynamic attenuator was controlled to reduce dose while maintaining peak variance by solving a convex optimization problem, assuming a priori knowledge of the patient anatomy. In targeted reconstruction tasks, the noise in specific regions was given increased weighting. A system with a standard attenuator (or 'bowtie filter') was used as a reference, and used either convex optimized tube current modulation (TCM) or a standard TCM heuristic. The noise of the scan was determined analytically while the dose was estimated using Monte Carlo simulations. Scatter was also estimated using Monte Carlo simulations. The sensitivity of the dynamic attenuator to patient centering was also examined by shifting the abdomen in 2 cm intervals. Results: Compared to a reference system with optimized TCM, use of the dynamic attenuator reduced dose by about 30% in routine scans and 50% in targeted scans. Compared to the TCM heuristics which are typically used without a priori knowledge, the dose reduction is about 50% for routine scans. The dynamic attenuator gives the ability to redistribute noise and variance and produces more uniform noise profiles than systems with a conventional bowtie filter. The SPR was also modestly reduced by 10% in the thorax and 24% in the abdomen. Imaging with the dynamic attenuator was relatively insensitive to patient centering, showing a 17% increase in peak variance for a 6 cm shift of the abdomen, instead of an 82% increase in peak variance for a fixed bowtie filter. Conclusions: A dynamic prepatient x-ray attenuator consisting of multiple wedges is capable of achieving substantial dose reductions and modest SPR reductions. [ABSTRACT FROM AUTHOR]