Fast Detector Colbeth2005
Fast Detector Colbeth2005
ABSTRACT
This paper explores the potential of flat panel detectors in sub-second CT scanning applications. Using a PaxScan
4030CB with 600um thick CsI(Tl), a central section of the panel (16 to 32 rows), was scanned at frame rates up to
1000fps. Using this platform, fundamental issues related to high speed scanning were characterized. The offset drift of
the imager over 60 seconds was found to be less than 0.014 ppm/sec relative to full scale. The gain stability over a 10
hour period is better than +/- .45%, which is at the resolution limit of the measurement. Two different types of lag
measurements were performed in order to separate the photodiode array lag from the CsI afterglow. The panel lag was
found to be 0.41% 1st frame and 0.054% 25th frame at 1000fps. The CsI(Tl) afterglow, however, is roughly an order of
magnitude higher, dominating the lag for sub-second scans. At 1000fps the 1st frame lag due to afterglow was 3.3% and
the 25th frame lag was 0.34%. Both the lag and afterglow are independent of signal level and each follows a simple
power law evolution versus time. Reconstructions of anatomical phantoms and the CATPHAN 500 phantom are
presented. With a 2 second, 1200 projection scan of the CATPHAN phantom at 600fps in 32 slice mode, using
120kVp and CTDI100 of 43.2mGy, 0.3% contrast resolution for a 6mm diameter target, can be visualized. In addition,
15lp/cm spatial resolution was achieved with a 2mm slice and a central CTDI100 of 10.8mGy.
Keywords: Flat Panel Detector, FPD, CT, Cone Beam CT, amorphous silicon, high frame rate
1. INTRODUCTION
The purpose of this work is to understand the potential of flat panel technology in sub-second volumetric CT. Since
much of the experimental work is done with a commercially available FPD (the 4030CB) [1], this work also
demonstrates the high frame rate capability of the imager and will hopefully stimulate research on new cone beam CT
(CBCT) applications for FPDs.
Traditional 3rd generation, 64 slice CT scanners are available from multiple manufacturers and a prototype 256 slice
scanner is in clinical trials [2]. The detectors used in these scanners are highly optimized to their task, having dynamic
range on the order of 20 bits and highly absorbing scintillators with low afterglow. However, traditional CT detectors
are also relatively expensive, since each pixel has its own electronics channel, the scintillators are pixelated through a
mechanical process and the underlying visible photon detectors are made from specialized arrays of crystalline Si
photodiodes [3].
In contrast, Flat Panel Detectors (FPDs) utilize multiplexing inside the array, so that the electronics channels are one per
data column. FPDs have typically used grown, columnar CsI, but mechanically pixelated scintillators have been used
in MV applications [4]. A single flat panel detector can be used for multiple applications, including fluoroscopy,
radiography, cine, angiography, DSA and more recently Cone Beam CT. Flat Panel Detectors (FPDs) generally have
many more rows of detectors than even the most advanced multi-slice CT detector and simultaneously, resolution well
beyond that available in CT scanners today. However, FPDs began in relatively low frame rate, low dynamic range
applications. Recently, we have demonstrated expanded dynamic range flat panel detectors in Cone Beam applications,
which approach the contrast resolution of traditional CT scanners [1]. And although the total time required to scan an
organ is roughly the same, the frame rate of flat panel CBCT is generally so low (30fps) that motion artifacts are a
significant issue.
Medical Imaging 2005: Physics of Medical Imaging, edited by Michael J. Flynn, 387
Proceedings of SPIE Vol. 5745 (SPIE, Bellingham, WA, 2005)
1605-7422/05/$15 · doi: 10.1117/12.595631
4030CB
1.24 x 40 cm
(16 x 1024)
to
2.48 x 40cm
(32 x 1024)
Rotation
Stage
Collimator
&
G242
x-ray tube
Figure 1. Cone Beam CT test stand. The high speed central strip location is marked.
The experimental work utilized a commercially available 4030CB with custom configuration files designed to scan a
narrow strip in the center of the detector containing a reduced number of rows. This is a straightforward extension of
the zoom capability of the imager and it is possible to switch between these high speed strip modes and normal full field
of view (FOV) modes in just a few frames. It is desirable to continuously discharge the unused pixels to some level, so
that they don’t saturate from dark current build up or scattered x-rays. Thus there is some overhead time associated
with the read time of the central strip. For our setup, it was possible to scan a central 16 row section of the 4030CB at
frame rates up to 1000fps.
Figure 1 shows the panel, x-ray tube and rotation stage. Dynamic Gain Switching (DGS) was used in order to achieve
the highest frame rates with extended bit depth [1]. DGS modes have low and high sensitivity ranges. The panel
switches to the appropriate signal sensitivity automatically, on a pixel-by-pixel basis, depending on the signal level.
The resulting linear range of the detector is in excess of 80,000 counts, or 16.5 bits. For most of the measurements the
pixels were 2x4 binned, resulting in a 388 x 776um pixels, in rows 1024 pixels across. The number of rows (slices) was
varied from 16 to 32 depending on the frame rate and the particular experiment. The panel was externally triggered by
a pulse generator and the choice of frame rate was continuously variable up to 1000fps. The source-to-imager distance
(SID) was 150cm and the source-to-axial distance (SAD) was typically 100cm. The Varian G242 x-ray tube has 0.4
and 0.8mm spot sizes.
1.2
1.0
0.8
Signal
0.0
-0.2
-5 0 5 10 15
time (msec)
Figure 2. X-ray beam signal compared the average signal from the panel at 1000fps.
The above measurement certainly contains both the array lag and the signal from the CsI afterglow. From the data sheet
for our CsI coating, the afterglow is expected to decay to approximately 0.1% after 10msec. With this in mind a special
version of firmware was constructed in which a single row could be held in accumulation mode, until the afterglow
decayed sufficiently, and then scanned repeatedly at high frame rate. Figure 3 shows the experimental timing for the
2000fps (0.5msec) line scanning. Note that the 1st reading of the line contains the x-ray signal plus the accumulated
signal from the CsI afterglow. This 1st reading, in theory, does not contain the lag of the array, since the amount of
trapped charge should remain in steady state while the signal population on the photodiode remains relatively constant.
The trapped charge depopulation and lag decay start when the scanning begins. The array lag was calculated from
residual signal starting with the 2nd frame. Using this pulsed fluoroscopy method it was possible to use extensive
averaging to reduce the noise in the measurement.
Figure 3. Timing of the experimental procedure used to measure the array lag.
The gain stability is of interest over a period of a day, since CT scanners are typically gain calibrated once a day. The
gain stability was measured with the 4030CB in full FOV mode. A flat, open field image with 90% recursive filtering
was periodically acquired in pulsed fluoroscopy mode, with an x-ray technique of 120kVp, 50mA, 5msec and 1.25mm
Cu plus 1mm Al filtration. A dosimeter was mounted on the front surface of the imager and allowed to accumulate
greater than 20mR for each reading. Using the dosimeter readings and the average counts from the panel, the sensitivity
in counts/uR is calculated and tracked over time. The estimated standard deviation in the measurement, due to
dosimetry error and noise in the imager flat file, is +/- 0.45%.
2.4 DQE
The DQE is a measure of the efficiency with which a system utilizes incoming quanta. It requires knowledge of
the Modulation Transfer Function (MTF), the Noise Power Spectrum (NPS) and the average x-ray flux. A useful
practical definition applicable for linear systems is [5]:
d * MTF 2 (u )
( )
2
DQE q, u =
q * NPS q, u ( ) (1)
where the DQE is expressed as a function of the spatial frequency u and the photon flux q , which is expressed
2
as phtotons / mm / mR ; d is the average pixel signal expressed in the same units as those used in the
calculation of NPS.
We calculated the photon flux using the x-ray spectral modeling methodology of Boone and Seibert [6] for the
actual conditions used in our experiments, and also cross-checked the calculation method against one of the
values defined in ISO 92361 obtaining agreement within ~1.5%.
The NPS was obtained by a method similar to that described by many authors including Dobbins [7]. Flat fields
taken at known x-ray exposures, were processed in regions, typically 128x128 pixels. For each individual region,
the 2D FFT is calculated and the results averaged for many regions. The normalization of the FFT is checked to
ensure that the squared Fourier amplitudes sum to the pixel variance. The 2D NPS is calculated from the 2D FFT
2
by scaling by 1 where f inc is the frequency increment given by f inc = 1 /( N * pixel _ side) and N
f inc
is the size of the region used (128 typically). Providing there is no directional preference indicated in the data the
1D NPS is obtained using all of the points in the 2D plane, excluding those on-axis, which are most susceptible to
artifacts.
Finally, it is necessary to correct for any lag in the imager, as discussed by Granfors and Aufrichtig [8]. The lag
compensation was done by scaling the DQE by (1 − lag1 ) , where lag 1 is the first frame lag.
2
In the data presented below, a pulsed fluoroscopy source was used. The tube voltage was 70kVp with a beam
hardness of 7.1mm Al half value layer (RQA5). For a given DQE vs spatial resolution curve, every data point is
plotted and the scatter in the data indicates the uncertainty in the measurements. Because of the size of our
ionization chamber, the DQE was measured in the full FOV, 2x2 binning, analog gain=1, DGS mode at 30fps.
Table 1. CTDI measurements at a technique of 120 kV @ 200 mAs, slice thickness of 22.5 mm, beam
filtration of 4.8 mm Al equivalents and SAD of 100 cm.
Phantom CTDI100 Center CTDI100 Periphery CTDIW
(mGy/100 mAs) (mGy/100 mAs) (mGy/100 mAs)
Head acrylic 4.5 5.8 5.4
Body acrylic 1.5 4.9 3.8
Air 7.8
3. RESULTS
3.1 DQE in Dynamic Gain Switching Mode
As discussed briefly above, in Dynamic Gain Switching (DGS) mode the panel switches between high and low
sensitivity on a pixel-by-pixel basis. The output data contains a 14bit measurement plus a gain bit for each pixel.
The offset & gain normalization is then applied post image capture. Figure 4 shows the Signal and SNR vs.
mAs/frame for the DGS mode running with 16 rows, 2x4 binning, analog gain equal 1 and 1000fps. Figure 5 is
the measured DQE with beam quality RQA5, for the DGS mode with 2x2 binning and analog voltage gain
= 1 in full FOV mode at 30fps.
1000
80,000
70,000
60,000
100
50,000
Signal
SNR
40,000
30,000
10
20,000
10,000
0
1
0.000 0.050 0.100 0.150 0.200 0.250 0.300 0.350
0.001 0.010 0.100 1.000
mAs/frame mAs/frame
Figure 4. Signal and SNR vs. mAs/frame for 2x4 binning,, 16 x 1024 DGS mode running at 1000fps.
0.7
0.6
5.1 nGy
0.5
13.1 nGy
DQE
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
Spatial Frequency (lp/mm)
Figure 5. DQE of the DGS mode with full FOV 2x2, Gain=1 at 30fps. The beam quality was RQA5.
1713
y = 0.0011x + 1711.7
1712.5
Average Offset
1712
(adu's)
1711.5
1711
1710.5
0 10 20 30 40 50 60
Time (seconds)
Figure 6. Offset drift versus time at 1000fps.
The gain stability measurements seen in Figure 7, show that the peak-to-peak variation in sensitivity over a 10 hour
period was 0.9%. However the standard deviation of the measurement is +/- 0.45%, so the measured gain variations
appear to be limited by the uncertainty in the experiment.
Sensitivity (counts/uR)
98.00
average = 96.9
97.00
96.00
95.00
1 10 100 1000
Time (minutes)
100.00%
10.00%
Residual Signal
(%)
1.00%
0.10%
0.01%
0.001 0.01 0.1
Time (seconds)
Figure 8. CsI afterglow vs. time for input levels from 1% to 75% of the maximum linear input dose.
0.10%
0.01%
1 10 100
Frame #
10.00%
1000fps
1.00%
Residual Signal
y = 0.0651x-0.9025
(%)
0.10%
0.01%
0.10 1.00 10.00 100.00 1000.00
Time (msec)
Figure 10. Residual signal due to photodiode array lag (Method 2) and CsI afterglow (Method 1), at 1000fps.
Figure 11 shows an axial slice from a high frame rate scan of a chest phantom, and for comparison, an image of the
same phantom taken on a commercial CT scanner. In order to generate a comparable level of dose, the flat panel
data is the sum of four individual scans. The flat panel scan has not been corrected for beam hardening, the effects
of afterglow, or scatter.
(a) (b)
Figure 11. (a) Chest phantom scanned over 2 sec, 1200 views, taken at 120kVp, 2x4 binning, 32x1024
matrix. An additional 2x1 binning during reconstruction, created pixels 0.776mm square. The
CTDIW = 57mGy. (b) A scan for the same phantom taken on a commercial scanner is shown on
the right for comparison.
In Figures 12 and 13 are the results of scans of the CATPHAN 500 phantom. In the resolution phantom 15lp/cm can be
resolved. This is less than expected, since scans of the same phantom at 30fps have achieved greater than 20lp/cm.
Since no lag correction has been applied, the degradation in resolution may be attributable to the CsI afterglow, but this
is still under investigation. The scan of the low contrast phantom revealed that the 0.3% contrast objects could be
visualized for target diameters as low as 6mm. Increasing the re-binning in the axial plane did reduce the high
frequency noise observed, but had no significant affect on the observable contrast resolution.
(c)
Figure 12. (a) + (c) 2 sec, 1200 view scan of CATPHAN 500 line pair per cm high resolution gauge taken at 120kVp,
2x4 binning, 32x1024 matrix. The slice thickness was 2mm and the effective patient dose at the center was 10.8mGy
CTDI100. The final pixel size for projection images was 0.388 (tangential) and 0.776 (axial). (b) The calibration target
on the right was also scanned over 2 seconds, 1200 views, but with a 0.776 x 0.776 pixel size. The effective patient
dose was 10.8mGy CTDI100 at patient center. The slice thickness was 5mm.
Figure 13. CATPHAN low contrast phantom, 10mm slice acquired at 600fps, 2x4 pixels, 32 x 1024 matrix, 120kVp
and CTDI100 (central) = 43.2 mGy. The image on the right is reconstructed with additional 4x1 binning, pixel size
1.552 (tangential) x 0.776(axial) mm at the detector. The image on the left is reconstructed with an additional 2x1
binning for a pixel size of 0.776 x 0.776 mm.
4. CONCLUSIONS
Research and applications of flat panel based cone beam CT have increased significantly over the last few years;
however, all this work has, of necessity, been limited to relatively slow scans at low imager frame rates. In this
investigation we have characterized the behavior of a flat panel detector at frame rates compatible with sub-second
scanning. More broadly, this work explores the potential for multiplexed detectors in CT scanning.
The results are very promising. One of the early concerns was the lag behavior of the amorphous silicon photodiode
array, since at low frame rates the short term lag can be near 10% and the long-term lag on the order of 1% for many
seconds. At high frame rates, the integration time for capturing the lag is so reduced that the issue of residual signal
becomes dominated by the scintillator afterglow, which is essentially the situation in state-of-the-art CT scanners.
Stability was also a concern. Although flat panel imagers have always used offset and gain normalization like CT
detectors, the level of stability required is much higher in CT applications. Our results show that the offset stability is
extremely good and that the gain stability appears adequate over the time scales required for CT. Using the extended
dynamic range capability of the 4030CB panel, we have demonstrated 0.3% contrast resolution for a 6mm diameter
target.
The authors would like to acknowledge the valuable discussion and insight of M. Grasruck, Siemens Medical
Solutions, in applying the dynamic gain switching mode to CT.
REFERENCES
1. P.G. Roos et al., “Multiple gain ranging readout method to extend the dynamic range of amorphous silicon flat
panel imagers,” Proceedings of SPIE Medical Imaging, 5368-16, Feb. 2004, San Diego, CA.
2. 6th Annual International Symposium on Multidetector-Row CT, Ritz-Carlton, San Francisco, CA, June 23-26,
2004.
3. S. Utrup et al. “Design and performance of a 32-slice CT detector system using back-illuminated
photodiodes,” Proceedings of SPIE Medical Imaging, 5368-06, Feb. 2004, San Diego, CA.
4. E.J. Seppi et al., “Megavoltage cone-beam CT using a high efficiency image receptor,” Int. J. Radiation
Oncology Biol. Phys., Vol.55, No.3, pp. 793-803, 2003.
5. I.A. Cunningham in Handbook of Medical Imaging (HMI), Volume 1, p124, Edited by J. Beutel, H.L. Kundel
and R.L. Van Metter. SPIE press copyright 2000.
6. J.M. Boone and J.A.Seibert, Medical Physics 24, 1661 (1997)
7. Dobbins in HMI p197 et seq.
8. P.R. Granfors, R. Aufrichtig, “DQE(f) of an amorphous-silicon flat-panel x-ray detector: detector parameter
influences and measurement methodology,” Proc. SPIE Med.Img.2000, vol.3977, pg.2-13, 13-15 Feb 2000.