PURPOSE
Detector lag, or residual signal, in a-Si flat-panel (FP) detectors can cause significant shading artifacts in cone-beam computed tomography reconstructions. To date, most correction models have assumed a linear, time-invariant (LTI) model and correct lag by deconvolution with an impulse response function (IRF). However, the lag correction is sensitive to both the exposure intensity and the technique used for determining the IRF. Even when the LTI correction that produces the minimum error is found, residual artifact remains. A new non-LTI method was developed to take into account the IRF measurement technique and exposure dependencies.
METHODS
First, a multiexponential (N = 4) LTI model was implemented for lag correction. Next, a non-LTI lag correction, known as the nonlinear consistent stored charge (NLCSC) method, was developed based on the LTI multiexponential method. It differs from other nonlinear lag correction algorithms in that it maintains a consistent estimate of the amount of charge stored in the FP and it does not require intimate knowledge of the semiconductor parameters specific to the FP. For the NLCSC method, all coefficients of the IRF are functions of exposure intensity. Another nonlinear lag correction method that only used an intensity weighting of the IRF was also compared. The correction algorithms were applied to step-response projection data and CT acquisitions of a large pelvic phantom and an acrylic head phantom. The authors collected rising and falling edge step-response data on a Varian 4030CB a-Si FP detector operating in dynamic gain mode at 15 fps at nine incident exposures (2.0%-92% of the detector saturation exposure). For projection data, 1st and 50th frame lag were measured before and after correction. For the CT reconstructions, five pairs of ROIs were defined and the maximum and mean signal differences within a pair were calculated for the different exposures and step-response edge techniques.
RESULTS
The LTI corrections left residual 1st and 50th frame lag up to 1.4% and 0.48%, while the NLCSC lag correction reduced 1st and 50th frame residual lags to less than 0.29% and 0.0052%. For CT reconstructions, the NLCSC lag correction gave an average error of 11 HU for the pelvic phantom and 3 HU for the head phantom, compared to 14-19 HU and 2-11 HU for the LTI corrections and 15 HU and 9 HU for the intensity weighted non-LTI algorithm. The maximum ROI error was always smallest for the NLCSC correction. The NLCSC correction was also superior to the intensity weighting algorithm.
CONCLUSIONS
The NLCSC lag algorithm corrected for the exposure dependence of lag, provided superior image improvement for the pelvic phantom reconstruction, and gave similar results to the best case LTI results for the head phantom. The blurred ring artifact that is left over in the LTI corrections was better removed by the NLCSC correction in all cases.
[1]
Rebecca Fahrig,et al.
A forward bias method for lag correction of an a-Si flat panel detector.
,
2011,
Medical physics.
[2]
H. Wieczorek.
Effects of Trapping in a-Si:H Diodes
,
1995
.
[3]
Evangelos Matsinos.
Current status of the CBCT project at Varian Medical Systems
,
2005,
SPIE Medical Imaging.
[4]
Torsten Solf,et al.
Temporal artifacts in flat dynamic x-ray detectors
,
2001,
SPIE Medical Imaging.
[5]
Isaias D. Job,et al.
Multidetector-row CT with a 64-row amorphous silicon flat panel detector
,
2007,
SPIE Medical Imaging.
[6]
S. N. Friedman,et al.
A moving slanted-edge method to measure the temporal modulation transfer function of fluoroscopic systems.
,
2008,
Medical physics.
[7]
Josh Star-Lack,et al.
SU‐FF‐I‐04: A Fast Variable‐Intensity Ring Suppression Algorithm
,
2006
.
[8]
J H Siewerdsen,et al.
Cone-beam computed tomography with a flat-panel imager: effects of image lag.
,
1999,
Medical physics.
[9]
G. Virshup,et al.
TH‐E‐330A‐06: Scatter Characterization in Cone‐Beam CT Systems with Offset Flat Panel Imagers
,
2006
.
[10]
D. Jaffray,et al.
A ghost story: spatio-temporal response characteristics of an indirect-detection flat-panel imager.
,
1999,
Medical physics.
[11]
R. Fahrig,et al.
TH‐E‐330A‐04: Investigation Into the Cause of a New Artifact in Cone Beam CT Reconstructions On a Flat Panel Imager
,
2006
.
[12]
Jiang Hsieh,et al.
Recursive correction algorithm for detector decay characteristics in CT
,
2000,
Medical Imaging.
[13]
Carlo Tognina,et al.
Photodiode forward bias to reduce temporal effects in a-Si based flat panel detectors
,
2008,
SPIE Medical Imaging.
[14]
J H Siewerdsen,et al.
An empirical method for lag correction in cone-beam CT.
,
2008,
Medical physics.
[15]
Rebecca Fahrig,et al.
Investigation into the optimal linear time-invariant lag correction for radar artifact removal.
,
2011,
Medical physics.
[16]
Kai Yang,et al.
Computer modeling of the spatial resolution properties of a dedicated breast CT system.
,
2007,
Medical physics.
[17]
S. N. Friedman,et al.
A small-signal approach to temporal modulation transfer functions with exposure-rate dependence and its application to fluoroscopic detective quantum efficiency.
,
2009,
Medical physics.
[18]
Mark A. Hartney,et al.
High-performance amorphous silicon image sensor for x-ray diagnostic medical imaging applications
,
1999,
Medical Imaging.
[19]
L. Feldkamp,et al.
Practical cone-beam algorithm
,
1984
.
[20]
Carlo Tognina,et al.
Parameter investigation and first results from a digital flat panel detector with forward bias capability
,
2008,
SPIE Medical Imaging.