Development of a Single Detector Ring Micro Crystal Element Scanner: QuickPET II

MiCE Detector Modules

We have previously reported on the development of different MiCE detector designs [5,11]. The MiCE detectors used in QuickPET II consist of a 22 × 22 array of polished 0.8 × 0.8 × 10 mm mixed lutetium silicate crystals. The crystals are placed within a grid made of a highly reflective polymer film material [13] following the technique described by Miyaoka et al. [5,11]. The grid serves three purposes: (1) it optically isolates the crystals; (2) it functions as a reflective wrap; and (3) it provides structural support for the crystal array. The center-to-center spacing of the crystals within the grid is 0.9 mm. The crystal array is directly coupled to a 6 + 6 cross-anode position-sensitive photomultiplier tube (PMT) (R8520-00-C12, Hamamatsu Photonics K.K., Japan) using Histomount (National Diagnostics, Atlanta, Georgia) glue [14]. A sample crystal array and position-sensitive PMT is pictured in Figure 1.

OPEN IN VIEWER

Figure 1.

Picture of MiCE crystal array and Hamamatsu position-sensitive PMT.

An image, referred to as a crystal map (see Figures 4A and 10), illustrating the decoding performance of a MiCE detector module obtained using the techniques described by Miyaoka et al. [15] was obtained. The crystal map is a two-dimensional (i.e., x-axis and y-axis) histogram produced by mapping the ratios of the electronic signals acquired from a uniform photon flux. Each dot in the crystal map represents a crystal element. The crystal map is used to generate a look up table (LUT) that designates the crystal of interaction for a given set of position signals. The decoding characteristics of a detector module are given by the peak (i.e., the maximum value) to valley (i.e., the minimum value) ratios between neighboring crystal elements in the crystal map.

To maximize the decoding capabilities of the detector module, positioning of the crystal array on the face of the PMT is a multistep process. The initial positioning is done using a specially designed jig that positions the crystal array with respect to the external structure of the PMT. The long curing time of Histomount allows crystal maps to be collected before the glue fully sets. If the initial crystal maps indicate a slight misalignment between the crystal array and the PMT optics, the position of the crystal array is manually adjusted. This process can be repeated if necessary. Because each crystal array is individually positioned with respect to the external structure of the PMT, a crystal array alignment correction that is applied to the acquired data was developed. This procedure is described below in the subsection, “Module Alignment Correction.”

A second performance characteristic of a detector module is the energy resolution of the individual crystal elements. After the crystal boundary map (i.e., the LUT) is determined, an energy histogram (i.e., energy spectrum) for each crystal element is created from the collected data. The energy resolution is determined from the energy histogram. Sample energy spectra are illustrated in Figure 4C. The peak in the energy spectrum corresponds to 511 keV photons that are photoelectrically absorbed by the detector array and is referred to as the photopeak. The energy resolution of a detector is described by the full width at half maximum (FWHM) of the photopeak divided by the bin value corresponding to the highest point of the photopeak. The lower energy events seen in the energy spectrum correspond to Compton scattered photons. The energy resolution of the crystals can affect the choice of the lower energy threshold for the detector system. However, for a scintillator-based PET imaging system optimized to image mice this generally is not the case.

OPEN IN VIEWER

Figure 4.

(A) An image, referred to as a crystal map, illustrating the decoding characteristics for a sample MiCE detector module. Each dot in the image represents a crystal element. (B) A profile of (A) through the center row helps to further elucidate the decoding performance of a detector module. The average peak to valley ratio is ~7. (C) Sample energy spectra from individual MiCE crystals. The average energy resolution for a detector module is ~22%.

QuickPET II System Electronics and Gantry

A block diagram of the QuickPET II system is shown in Figure 2. The system is made up of 18 detector modules. The detector ring is divided into three zones. Coincidences are only allowed between modules in different zones. Each module has its own “front-end” electronics card that produces four position signals (X+, X-, Y+, Y-) and a summed signal. The four position signals are used to determine which crystal within the 22 × 22 array detected the annihilation photon associated with the injected radiotracer. The summed signals are individually fed to constant fraction discriminators (CFDs): one CAEN N843, 16 channel CFD, and two TENNELEC TC453 single-channel CFDs. The CFDs produce trigger signals that are fed to a coincidence processor and an analog multiplexer board, which was developed in-house. The coincidence processor unit sends control signals (i.e., module identifier information) to the analog multiplexers informing them which set of signals to pass on. The control signals are also routed to digital-to-analog converters (DACs) that produce analog detector module identifier signals, which in turn are collected by the acquisition electronics. The acquisition electronics consists mainly of nuclear instrumentation modules and computer automated measurement and control (CAMAC) modules and has been fully described in Ref. [16].

OPEN IN VIEWER

Figure 2.

Block diagram of data and control signal flow for the single ring evaluation system (QuickPET II).

The QuickPET II detectors, gantry and imaging table are pictured in Figure 3. The detector system consists of 18 detector modules. It has an inner ring diameter of 12.8 cm. The in-plane FOV is 5.76 cm and the axial FOV is 1.98 cm. The detector ring is mounted on a bearing that allows partial rotation of the detectors (± ~20°). The purpose of adding a partial rotation was to facilitate sampling the lines of response associated with the gaps between modules. To minimize the number of rotations required per scan, the length of the animal is scanned with the detectors in one position; the detectors are rotated 10° (one half the step between modules); and the animal is scanned in the reverse direction. Due to limitations in the acquisition electronics, the rotation angle of the gantry is not monitored during data acquisition. Therefore, the detectors are not rotated during data acquisition. The imaging table (Summit Medical Equipment, Oregon) was incorporated in such a manner as to allow delivery and exhaust of gas anesthesia from a single end of the table. The tube labeled in Figure 3 is for gas delivery while channels in the table allow the exhaust to be collected from the black tube at the end of the table. Three linear stages and one vertical stage allow for accurate table placement. The table has a linear travel range of ~9 cm, enough for whole-body scanning of a mouse. For whole-body imaging studies, the table is translated, under computer control, in 1.44 cm axial FOV steps. A slight overlap of adjacent axial FOVs is required because of the fully 3-D acquisition geometry.

OPEN IN VIEWER

Figure 3.

Photograph of QuickPET II evaluation system, including animal imaging table.

MiCE Detector Module Intrinsic Spatial Resolution

The intrinsic spatial resolution of a pair of MiCE detector modules was measured by stepping a line source through the axial imaging FOV of the system. The line source was a very small diameter glass capillary tube, outer diameter 0.30 mm and inner diameter unknown, filled with ¹⁸F. The length of the line source was 3.25 cm and it had a slight bow to it. For protection the line source was placed within a 22-gauge needle (outer diameter 0.711 mm; inner diameter 0.394 mm) and then mounted to the edge of the system's imaging table. The height of the table was adjusted so the line source was approximately centered between two opposing detectors. From a point just outside the system's axial imaging FOV, the line source was stepped in 0.1-mm increments across the system's imaging FOV. Profile plots of the coincidence counts for each of the 22 direct slices (opposing rows of crystals) were created. A Gaussian was fit to each of the plots to determine the FWHM for each of the 22 direct slices. Ideally, a point source would be used for this measurement; however, in order to get enough activity to do the measurement, a line source was required. Any misalignment between the line source and the crystal arrays will lead to a broadening of the FWHM of the measurement.

Module Alignment Correction

A procedure was developed to estimate the crystal array offsets before binning the data [17]. The offsets were estimated by comparing simulations and measurements of two axially oriented line sources located 0 and 14 mm radially from the center of the FOV. The 18 detector modules were considered as 9 pairs of opposing detectors. For each pair of detectors, the measured data were compared to the simulated data that represents truth. The detector modules were translated in 0.225-mm steps on a 5 × 5 grid and the position resulting in the minimum root-mean-square-error (RMSE) between the measured and simulated data was determined. The transverse offsets corresponding to the minimum RMSE were chosen as the locations of the detector blocks. The corrections do not account for rotational and z-axis offsets (i.e., detector ring diameter).

Line Phantom Images

Line source measurements were made using 0.3 mm inner diameter and 0.77 mm outer diameter glass capillary tubes. Two tubes were each filled with ~16 μCi (~592 KBq) of ¹⁸F. The tubes were placed parallel to one another and 14 mm apart on a thin, stiff board. Two tubes were imaged simultaneously to halve the time required for data acquisition. A tube separation of 14 mm was used so that the tails of the reconstructed line sources would not interfere with each other.

The capillary tubes were initially positioned at 0 mm (Tube 1) and 14 mm (Tube 2) radial distance from FOV center, respectively. The tubes were stepped outward in 2-mm increments until Tube 1 overlapped the initial position of Tube 2. A total 600,000 events (2000 counts/sec × 5 min) were scanned at each source position. The data were collected in listmode (i.e., event by event) format, which includes information on each event such as the detector coordinates and the energy values of the detected photons.

The detector ring was then rotated 10° counterclockwise. Data were collected as the capillary tubes were stepped back toward the center of the detector ring. By rotating the detector ring, lines of response associated with the gaps between adjacent detector modules were sampled. All scan times were adjusted to a decay corrected length corresponding to the initial acquisition.

Mouse Studies

A C3-1 TAg [18] mouse (weight, 35 g) that spontaneously develops mammary tumors was imaged twice using F-18-2-fluoro-2-deoxyglucose (FDG). The first time it was imaged, the mouse was injected via the tail vein with 340 μCi (12.58 MBq) of FDG and placed back into its cage during the uptake phase of the FDG. Imaging for the mouse began 3 hr 40 min after injection. The long delay between injection and imaging was due to two factors: (1) to let the activity decay to a level more optimal for the QuickPET II scanner and (2) a second mouse was imaged. The mouse was imaged again 10 days later. For the second study, the mouse was injected with 200 μCi (7.4 MBq) FDG via the tail vein. After injection, it was kept in the induction box under a low level of anesthesia (~0.5% isoflurane) for 50 min before imaging began. Four axial FOV scans were acquired for both imaging procedures.

A p53 heterozygote female mouse (weight, 30 g) on an NIH background [19] with chemically induced skin tumors was imaged with FDG. The mouse was injected with 240 μCi (8.88 MBq) of FDG via tail vein injection. Imaging began approximately 1 hr postinjection. To reduce uptake in the background tissues, the mouse was kept under light anesthesia (0.5–2.0% isoflurane) during radiopharmaceutical uptake. A 5 axial FOV scan was acquired over ~90 min; 21.3 million events were acquired during the study. Due to a handling error, the mouse was fed just before imaging and there was intense cardiac uptake as well as the expected tumor uptake.

Image Reconstruction

Fully 3-D sinogram data were rebinned into a stacked set of contiguous 2-D sinograms using either the single slice rebinning [20] or Fourier rebinning (FORE) [21] algorithms. The rebinning process accelerates the image reconstruction process by reducing the fully 3-D raytracing problem to a stack of much faster 2-D raytracings. The module alignment corrections were applied before rebinning of the data. After rebinning, data were corrected for sampling density. No other corrections (e.g., attenuation, scatter, randoms, or other normalization corrections) are currently being applied. Reconstruction algorithm options include 2-D versions of filtered-backprojection (FBP) [22] and ordered-subsets expectation-maximization (OSEM) [23,24] using the ASPIRE library [25]. In addition, we have also developed an iterative reconstruction utilizing a factorized system model including detector blurring (DB). This method is referred to as OSEM(DB) [17]. The FBP images were reconstructed with a Hamming filter at the Nyquist frequency (i.e., 0.9 mm FWHM). A 0.9-mm Gaussian filter was applied to the OSEM images. No postfiltering was applied to the OSEM(DB) images. For both OSEM reconstructions, eight iterations and eight subsets were used.

MiCE Detector Modules

A detector crystal map image is illustrated in Figure 4A. In one dimension, the peak and valley values are best visualized by looking at a profile through a row of crystals in the crystal map, as shown in Figure 4B. The average peak-to-valley ratio for our detector modules was ~7. Sample energy spectra from both a crystal with good light collection efficiency and low light collection efficiency are shown in Figure 4C. The average energy resolution for the individual crystals was ~22% (<20% if edge crystals are not included).

The production MiCE detectors used for QuickPET II have slightly degraded decoding performance along the edges of the crystal arrays compared to the initial prototype detectors [10]. An example of a detector with poor decoding along one of its edges is illustrated in Figure 5. The dark dots on the right side of the crystal map correspond to two crystals. The gray tails along the edge of the crystal map correspond to low energy events and are generally not included in the individual crystal look up tables (i.e., they are excluded once energy thresholds are applied on a crystal by crystal basis). The main reason for the slightly degraded performance is that the center-to-center crystal spacing was increased from 0.88 mm for the prototype crystal arrays to 0.90 mm for the production crystal arrays. This 0.02-mm increase in the crystal spacing translates into a 0.42-mm increase in the linear size of the crystal array that leads to decreased light collection efficiency for the edge crystals. To improve the decoding for crystals along the edge of the PMT, we are redesigning the front-end amplifier board to include threshold amplifiers. The use of threshold amplifiers has been shown to improve the edge decoding performance for detectors using Anger style positioning [26]. Additional reasons for the slightly degraded performance are that many of the current crystal arrays have been glued and unglued multiple times as the fabrication technique was being optimized and the use of analog multiplexers in the data acquisition electronics.

OPEN IN VIEWER

Figure 5.

An example of a crystal map with poor decoding of crystal elements along the right edge of the map.

MiCE Detector Module Intrinsic Spatial Resolution

The FWHM for the 22 direct slices of the detector module ranged from 0.99 to 1.15 mm, with a mean of 1.05 ± 0.07 mm, without correction for the line source (see Figure 6). The FWHM values were calculated from Gaussian fits to the measured data. The differences in the peak heights of the edge rows of crystals are caused by poor decoding of the edge crystals and the lower sensitivity of the edge crystals. The results were measured with the detectors in the gantry and the signals processed using all of the QuickPET II data acquisition electronics. Again, any rotational misalignments between the line source and the detector modules or the detector modules themselves will lead to a broadening of the measured intrinsic spatial resolution. Still these results are slightly better than the 1.21 mm average FWHM reported by Tai et al. [3] for the MicroPET II detectors measured in the system using a 0.5-mm diameter ²²Na point source and equal to the 1.05-mm FWHM reported by Chatziioannou et al [27] for a pair of prototype detectors and a 0.33-mm diameter ¹⁸F needle source.

OPEN IN VIEWER

Figure 6.

Intrinsic spatial resolution of the system measured with an ¹⁸F line source less than 0.4 mm in diameter. Average FWHM is 1.05 mm for the 22 direct slices without correction for the source size. FWHM values were determined from a Gaussian fit to the measured data.

Line Phantom Images

Sinograms of the multiline source phantom are illustrated in Figure 7. The data were binned using single slice rebinning with three plane differences. The sinograms include data with the gantry rotated to sample the gaps between the detectors. Figure 7A is a sinogram, without any alignment adjustment, for the line sources centered and at a 14-mm radial offset. Figure 7B-D contains sinograms after application of our detector alignment adjustment technique. The data were reconstructed using FBP and a ramp filter at Nyquist frequency. The reconstructed images are illustrated in Figure 8. The average FWHM and the average full width at tenth maximum (FWTM) values for the different line source positions are listed in Table 1. The values listed are for all 22 direct imaging planes and are not adjusted for the size of the source. Applying our current alignment correction led to a slight improvement in our image resolution. Methods to estimate residual misalignments for each of the 6 degrees of freedom associated with the location and orientation of each detector module are under development. We believe this is a critical component to fully realize the intrinsic spatial resolution characteristics of our MiCE detectors.

OPEN IN VIEWER

Figure 7.

(A) Sinogram of line sources at center and offset by 14 mm without alignment correction. (B) Sinogram of line sources at center and offset by 14 mm with alignment correction. (C) Sinogram of line sources offset by 4 and 18 mm. (D) Sinogram of line sources offset by 8 and 22 mm. All sinograms were rebinned using single slice rebinning and ±3 slices.

OPEN IN VIEWER

Figure 8.

(A) Reconstructed image of line sources at center and with a transverse offset of 14 mm without alignment correction. (B) Reconstructed image of line sources at center and offset by 14 mm with alignment correction. (C) Reconstructed image of line sources offset by 4 and 18 mm. (D) Reconstructed image of line sources offset by 8 and 22 mm. All images reconstructed with FBP and a ramp filter at Nyquist frequency.

OPEN IN VIEWER

Table 1.

Line Source Image Resolution

	Radial Profilesa		Tangential Profilesa
Position of Line Source	FWHM (mm)	FWTM (mm)	FWHM (mm)	FWTM (mm)
0 mm (no correction)	1.3 ± 0.05	3.0 ± 0.26	1.5 ± 0.38	2.9 ± 0.99
14 mm (no correction)	1.6 ± 0.05	4.1 ± 0.07	1.8 ± 0.36	4.5 ± 0.84
0 mm (with correction)	1.2 ± 0.06	3.0 ± 0.39	1.2 ± 0.05	2.9 ± 0.29
4 mm (with correction)	1.3 ± 0.05	2.9 ± 0.30	1.2 ± 0.04	2.7 ± 0.18
8 mm (with correction)	1.5 ± 0.09	3.8 ± 0.60	1.3 ± 0.08	3.7 ± 0.52
14 mm (with correction)	1.5 ± 0.10	3.8 ± 0.53	1.5 ± 0.16	4.4 ± 0.40
18 mm (with correction)	2.2 ± 0.13	4.7 ± 0.64	1.5 ± 0.12	4.8 ± 0.54
22 mm (with correction)	2.2 ± 0.11	5.0 ± 0.39	1.5 ± 0.16	5.2 ± 0.44

a

No correction for the size of line source was included.

Mouse Studies

The images of the C3-1 TAg mouse in Figure 9 illustrate how sedation of the mouse improves the uptake of FDG in tumors. These results are consistent with other FDG tumor imaging protocols we have conducted on tumor bearing mice and not due to difference in the uptake period for the two studies. In fact for human FDG imaging, lengthening the uptake period usually improves the tumor-to-background contrast [28]. Both datasets were binned using single slice rebinning with nine ring differences and reconstructed using FBP with a Hamming filter at 0.8 the Nyquist frequency (1.13 mm FWHM).

OPEN IN VIEWER

Figure 9.

(A) C3-1 TAg mouse imaged with FDG with no anesthesia during uptake phase. (B) Same mouse imaged 10 days later. For the second scan the mouse was under anesthesia during uptake phase of FDG.

The intense uptake in the heart of the mouse bearing the skin tumors provided an ideal dataset to examine differences between FBP, OSEM, and our OSEM(DB) image reconstructions. Mouse heart images were reconstructed with each of the three different reconstruction techniques (after FORE rebinning). Transaxial and coronal slices through the heart (non-reoriented) are illustrated in Figure 10. The images reconstructed using OSEM are slightly sharper than the FBP images. The OSEM(DB) images were significantly sharper due to the resolution recovery associated with modeling the 2-D detector response function within the system matrix. Hot spots are visible in each of the image sets. Because they appear in the FBP images, we believe they are either a normalization artifact or physiologic uptake differences in the heart rather than an artifact caused by the iterative reconstruction methods. Figure 11 shows a maximum intensity pixel display through the whole-body mouse image volume reconstructed by the FORE + OSEM(DB) approach. FDG uptake of skin tumors are well recognized with improved resolution by the FORE + OSEM(DB).

OPEN IN VIEWER

Figure 10.

Heart images of p53 heterozygote mouse. FORE + FBP, 0.9 mm Hamming filter: (A) Coronal slice; (B) Transaxial slice. FORE + OSEM, 0.9 mm Gaussian postfilter: (C) Coronal slice; (D) Transaxial slice. FORE + OSEM(DB): (E) Coronal slice; (F) Transaxial slice.

OPEN IN VIEWER

Figure 11.

Maximum intensity pixel projection image of p53 heterozygote mouse reconstructed by FORE + OSEM(DB). The same mouse data as Figure 10 was used for the reconstruction.