Dataset Information

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

ABSTRACT: 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.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.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.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.

SUBMITTER: Lazos D

PROVIDER: S-EPMC3443192 | biostudies-other | 2012 Sep

REPOSITORIES: biostudies-other

ACCESS DATA

Similar Datasets

Project description:To determine the mean and range of location-averaged breast skin thickness using high-resolution dedicated breast CT for use in Monte Carlo-based estimation of normalized glandular dose coefficients.This study retrospectively analyzed image data from a clinical study investigating dedicated breast CT. An algorithm similar to that described by Huang et al. ["The effect of skin thickness determined using breast CT on mammographic dosimetry," Med. Phys. 35(4), 1199-1206 (2008)] was used to determine the skin thickness in 137 dedicated breast CT volumes from 136 women. The location-averaged mean breast skin thickness for each breast was estimated and the study population mean and range were determined. Pathology results were available for 132 women, and were used to investigate if the distribution of location-averaged mean breast skin thickness varied with pathology. The effect of surface fitting to account for breast curvature was also studied.The study mean (± interbreast SD) for breast skin thickness was 1.44 ± 0.25 mm (range: 0.87-2.34 mm), which was in excellent agreement with Huang et al. Based on pathology, pair-wise statistical analysis (Mann-Whitney test) indicated that at the 0.05 significance level, there were no significant difference in the location-averaged mean breast skin thickness distributions between the groups: benign vs malignant (p = 0.223), benign vs hyperplasia (p = 0.651), hyperplasia vs malignant (p = 0.229), and malignant vs nonmalignant (p = 0.172).Considering this study used a different clinical prototype system, and the study participants were from a different geographical location, the observed agreement between the two studies suggests that the choice of 1.45 mm thick skin layer comprising the epidermis and the dermis for breast dosimetry is appropriate. While some benign and malignant conditions could cause skin thickening, in this study cohort the location-averaged mean breast skin thickness distributions did not differ significantly with pathology. The study also underscored the importance of considering breast curvature in estimating breast skin thickness.

Project description:To investigate the effect of the number of projection views on image noise in cone-beam CT (CBCT) with a flat-panel detector.This fairly fundamental consideration in CBCT system design and operation was addressed experimentally (using a phantom presenting a uniform medium as well as statistically motivated "clutter") and theoretically (using a cascaded systems model describing CBCT noise) to elucidate the contributing factors of quantum noise (σ(Q)), electronic noise (σ(E)), and view aliasing (σ(view)). Analysis included investigation of the noise, noise-power spectrum, and modulation transfer function as a function of the number of projections (N(proj)), dose (D(tot)), and voxel size (b(vox)).The results reveal a nonmonotonic relationship between image noise and N(proj) at fixed total dose: for the CBCT system considered, noise decreased with increasing N(proj) due to reduction of view sampling effects in the regime N(proj) <~200, above which noise increased with N(proj) due to increased electronic noise. View sampling effects were shown to depend on the heterogeneity of the object in a direct analytical relationship to power-law anatomical clutter of the form κ/f(β)--and a general model of individual noise components (σ(Q), σ(E), and σ(view)) demonstrated agreement with measurements over a broad range in N(proj), D(tot), and b(vox).The work elucidates fairly basic elements of CBCT noise in a manner that demonstrates the role of distinct noise components (viz., quantum, electronic, and view sampling noise). For configurations fairly typical of CBCT with a flat-panel detector (FPD), the analysis reveals a "sweet spot" (i.e., minimum noise) in the range N(proj) ~ 250-350, nearly an order of magnitude lower in N(proj) than typical of multidetector CT, owing to the relatively high electronic noise in FPDs. The analysis explicitly relates view aliasing and quantum noise in a manner that includes aspects of the object ("clutter") and imaging chain (including nonidealities of detector blur and electronic noise) to provide a more rigorous basis for commonly held intuition and heurism in CBCT system design and operation.

Project description:Cone-beam CT (CBCT) with a flat-panel detector (FPD) is finding application in areas such as breast and musculoskeletal imaging, where dual-energy (DE) capabilities offer potential benefit. The authors investigate the accuracy of material classification in DE CBCT using filtered backprojection (FBP) and penalized likelihood (PL) reconstruction and optimize contrast-enhanced DE CBCT of the joints as a function of dose, material concentration, and detail size.Phantoms consisting of a 15 cm diameter water cylinder with solid calcium inserts (50-200 mg/ml, 3-28.4 mm diameter) and solid iodine inserts (2-10 mg/ml, 3-28.4 mm diameter), as well as a cadaveric knee with intra-articular injection of iodine were imaged on a CBCT bench with a Varian 4343 FPD. The low energy (LE) beam was 70 kVp (+0.2 mm Cu), and the high energy (HE) beam was 120 kVp (+0.2 mm Cu, +0.5 mm Ag). Total dose (LE+HE) was varied from 3.1 to 15.6 mGy with equal dose allocation. Image-based DE classification involved a nearest distance classifier in the space of LE versus HE attenuation values. Recognizing the differences in noise between LE and HE beams, the LE and HE data were differentially filtered (in FBP) or regularized (in PL). Both a quadratic (PLQ) and a total-variation penalty (PLTV) were investigated for PL. The performance of DE CBCT material discrimination was quantified in terms of voxelwise specificity, sensitivity, and accuracy.Noise in the HE image was primarily responsible for classification errors within the contrast inserts, whereas noise in the LE image mainly influenced classification in the surrounding water. For inserts of diameter 28.4 mm, DE CBCT reconstructions were optimized to maximize the total combined accuracy across the range of calcium and iodine concentrations, yielding values of ∼ 88% for FBP and PLQ, and ∼ 95% for PLTV at 3.1 mGy total dose, increasing to ∼ 95% for FBP and PLQ, and ∼ 98% for PLTV at 15.6 mGy total dose. For a fixed iodine concentration of 5 mg/ml and reconstructions maximizing overall accuracy across the range of insert diameters, the minimum diameter classified with accuracy >80% was ∼ 15 mm for FBP and PLQ and ∼ 10 mm for PLTV, improving to ∼ 7 mm for FBP and PLQ and ∼ 3 mm for PLTV at 15.6 mGy. The results indicate similar performance for FBP and PLQ and showed improved classification accuracy with edge-preserving PLTV. A slight preference for increased smoothing of the HE data was found. DE CBCT discrimination of iodine and bone in the knee was demonstrated with FBP and PLTV at 6.2 mGy total dose.For iodine concentrations >5 mg/ml and detail size ∼ 20 mm, material classification accuracy of >90% was achieved in DE CBCT with both FBP and PL at total doses <10 mGy. Optimal performance was attained by selection of reconstruction parameters based on the differences in noise between HE and LE data, typically favoring stronger smoothing of the HE data, and by using penalties matched to the imaging task (e.g., edge-preserving PLTV in areas of uniform enhancement).

Project description:PURPOSE:Scatter is a major factor degrading the image quality of cone beam computed tomography (CBCT). Conventional scatter correction strategies require handcrafted analytical models with ad hoc assumptions, which often leads to less accurate scatter removal. This study aims to develop an effective scatter correction method using a residual convolutional neural network (CNN). METHODS:A U-net based 25-layer CNN was constructed for CBCT scatter correction. The establishment of the model consists of three steps: model training, validation, and testing. For model training, a total of 1800 pairs of x-ray projection and the corresponding scatter-only distribution in nonanthropomorphic phantoms taken in full-fan scan were generated using Monte Carlo simulation of a CBCT scanner installed with a proton therapy system. An end-to-end CNN training was implemented with two major loss functions for 100 epochs with a mini-batch size of 10. Image rotations and flips were randomly applied to augment the training datasets during training. For validation, 200 projections of a digital head phantom were collected. The proposed CNN-based method was compared to a conventional projection-domain scatter correction method named fast adaptive scatter kernel superposition (fASKS) method using 360 projections of an anthropomorphic head phantom. Two different loss functions were applied for the same CNN to evaluate the impact of loss functions on the final results. Furthermore, the CNN model trained with full-fan projections was fine-tuned for scatter correction in half-fan scan by using transfer learning with additional 360 half-fan projection pairs of nonanthropomorphic phantoms. The tuned-CNN model for half-fan scan was compared with the fASKS method as well as the CNN-based method without the fine-tuning using additional lung phantom projections. RESULTS:The CNN-based method provides projections with significantly reduced scatter and CBCT images with more accurate Hounsfield Units (HUs) than that of the fASKS-based method. Root mean squared error of the CNN-corrected projections was improved to 0.0862 compared to 0.278 for uncorrected projections or 0.117 for the fASKS-corrected projections. The CNN-corrected reconstruction provided better HU quantification, especially in regions near the air or bone interfaces. All four image quality measures, which include mean absolute error (MAE), mean squared error (MSE), peak signal-to-noise ratio (PSNR), and structural similarity (SSIM), indicated that the CNN-corrected images were significantly better than that of the fASKS-corrected images. Moreover, the proposed transfer learning technique made it possible for the CNN model trained with full-fan projections to be applicable to remove scatters in half-fan projections after fine-tuning with only a small number of additional half-fan training datasets. SSIM value of the tuned-CNN-corrected images was 0.9993 compared to 0.9984 for the non-tuned-CNN-corrected images or 0.9990 for the fASKS-corrected images. Finally, the CNN-based method is computationally efficient - the correction time for the 360 projections only took less than 5 s in the reported experiments on a PC (4.20 GHz Intel Core-i7 CPU) with a single NVIDIA GTX 1070 GPU. CONCLUSIONS:The proposed deep learning-based method provides an effective tool for CBCT scatter correction and holds significant value for quantitative imaging and image-guided radiation therapy.

Project description:To measure and investigate the improvement of microcalcification (MC) visibility in cone beam breast CT with a high pitch (75 μm), thick (500 μm) scintillator CMOS/CsI flat panel detector (Dexela 2923, Perkin Elmer).Aluminum wires and calcium carbonate grains of various sizes were embedded in a paraffin cylinder to simulate imaging of calcifications in a breast. Phantoms were imaged with a benchtop experimental cone beam CT system at various exposure levels. In addition to the Dexela detector, a high pitch (50 μm), thin (150 μm) scintillator CMOS/CsI flat panel detector (C7921CA-09, Hamamatsu Corporation, Hamamatsu City, Japan) and a widely used low pitch (194 μm), thick (600 μm) scintillator aSi/CsI flat panel detector (PaxScan 4030CB, Varian Medical Systems) were also used in scanning for comparison. The images were independently reviewed by six readers (imaging physicists). The MC visibility was quantified as the fraction of visible MCs and measured as a function of the estimated mean glandular dose (MGD) level for various MC sizes and detectors. The modulation transfer functions (MTFs) and detective quantum efficiencies (DQEs) were also measured and compared for the three detectors used.The authors have demonstrated that the use of a high pitch (75 μm) CMOS detector coupled with a thick (500 μm) CsI scintillator helped make the smaller 150-160, 160-180, and 180-200 μm MC groups more visible at MGDs up to 10.8, 9, and 10.8 mGy, respectively. It also made the larger 200-212 and 212-224 μm MC groups more visible at MGDs up to 7.2 mGy. No performance improvement was observed for 224-250 μm or larger size groups. With the higher spatial resolution of the Dexela detector based system, the apparent dimensions and shapes of MCs were more accurately rendered. The results show that with the aforementioned detector, a 73% visibility could be achieved in imaging 160-180 μm MCs as compared to 28% visibility achieved by the low pitch (194 μm) aSi/CsI flat panel detector. The measurements confirm that the Hamamatsu detector has the highest MTF, followed by the Dexel detector, and then the Varian detector. However, the Dexela detector, with its thick (500 μm) CsI scintillator and low noise level, has the highest DQE at all frequencies, followed by the Varian detector, and then the Hamamatsu detector. The findings on the MC visibility correlated well with the differences in MTFs, noise power spectra, and DQEs measured for these three detectors.The authors have demonstrated that the use of the CMOS type Dexela detector with its high pitch (75 μm) and thick (500 μm) CsI scintillator could help improve the MC visibility. However, the improvement depended on the exposure level and the MC size. For imaging larger MCs or scanning at high exposure levels, there was little advantage in using the Dexela detector as compared to the aSi type Varian detector. These findings correlate well with the higher measured DQEs of the Dexela detector, especially at higher frequencies.

Project description:X-ray scatter is a significant barrier to image quality improvements in cone-beam computed tomography (CBCT). A moving blocker-based strategy was previously proposed to simultaneously estimate scatter and reconstruct the complete volume within the field of view (FOV) from a single CBCT scan. A blocker consisting of lead stripes is inserted between the X-ray source and the imaging object, and moves back and forth along the rotation axis during gantry rotation. While promising results were obtained in our previous studies, the geometric design and moving speed of the blocker were set empirically. The goal of this work is to optimize the geometry and speed of the moving block system.Performance of the blocker was examined through Monte Carlo (MC) simulation and experimental studies with various geometry designs and moving speeds. All hypothetical designs employed an anthropomorphic pelvic phantom. The scatter estimation accuracy was quantified by using lead stripes ranging from 5 to 100 pixels on the detector plane. An iterative reconstruction based on total variation minimization was used to reconstruct CBCT images from unblocked projection data after scatter correction. The reconstructed image was evaluated under various combinations of lead strip width and interspace (ranging from 10 to 60 pixels) and different moving speed (ranging from 1 to 30 pixels per projection).MC simulation showed that the scatter estimation error varied from 0.8% to 5.8%. Phantom experiment showed that CT number error in the reconstructed CBCT images varied from 13 to 35. Highest reconstruction accuracy was achieved when the strip width was 20 pixels and interspace was 60 pixels and the moving speed was 15 pixels per projection.Scatter estimation can be achieved in a large range of lead strip width and interspace combinations. The moving speed does not have a very strong effect on reconstruction result if it is above 5 pixels per projection. Geometry design of the blocker affected image reconstruction accuracy more. The optimal geometry of the blocker has a strip width of 20 pixels and an interspace three times the strip width, which means 25% detector is covered by the blocker, while the optimal moving speed is 15 pixels per projection.

Project description:Mechanical instabilities that occur during gantry rotation of on-board cone-beam computed tomography (CBCT) imaging systems limit the efficacy of image-guided radiotherapy. Various methods for calibrating the CBCT geometry and correcting errors have been proposed, including some that utilize dedicated fiducial phantoms. The purpose of this work was to investigate the role of phantom fabrication imprecision on the accuracy of a particular CT cone-beam geometry estimate and to test a new method to mitigate errors in beam geometry arising from imperfectly fabricated phantoms.The authors implemented a fiducial phantom-based beam geometry estimation following the one described by Cho et al. [Med Phys 32(4), 968-983 (2005)]. The algorithm utilizes as input projection images of the phantom at various gantry angles and provides a full nine parameter beam geometry characterization of the source and detector position and detector orientation versus gantry angle. A method was developed for recalculating the beam geometry in a coordinate system with origin at the source trajectory center and aligned with the axis of gantry rotation, thus making the beam geometry estimation independent of the placement of the phantom. A second CBCT scan with the phantom rotated 180 degrees about its long axis was averaged with the first scan to mitigate errors from phantom imprecision. Computer simulations were performed to assess the effect of 2D fiducial marker positional error on the projections due to image discretization, as well as 3D fiducial marker position error due to phantom fabrication imprecision. Experimental CBCT images of a fiducial phantom were obtained and the algorithm used to measure beam geometry for a Varian Trilogy with an on-board CBCT.Both simulations and experimental results reveal large sinusoidal oscillations in the calculated beam geometry parameters with gantry angle due to displacement of the phantom from CBCT isocenter and misalignment with the gantry axis, which are eliminated by recalculating the beam geometry in the source coordinate system. Simulations and experiments also reveal an additional source of oscillations arising from fiducial marker position error due to phantom fabrication imprecision that are mitigated by averaging the results with those of a second CBCT scan with phantom rotated. With a typical fiducial marker position error of 0.020 mm (0.001 in.), source and detector position are found in simulations to be within 250 microm of the true values, and detector and gantry angles less than 0.2 degrees. Detector offsets are within 100 microm of the known value. Experimental results verify the efficacy of the second scan in mitigating beam geometry errors, as well as large apparent source/detector isocenter offsets arising from phantom fabrication imprecision.The authors have developed and validated a novel fiducial phantom-based CBCT beam geometry estimation algorithm that does not require precise positioning of the phantom at machine isocenter and is insensitive to positional imprecision of fiducial markers within the phantom due to fabrication errors. The method can accurately locate source and detector isocenters even when using an imprecise phantom, which is very important for measurement of isocenter coincidence of the therapy and on-board imaging systems.

Dataset Information

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

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets