Elec 2607 Assignment 1 Answer

Abstract

Positron emission tomography (PET) is a tool for metabolic imaging that has been utilized since the earliest days of nuclear medicine. A key component of such imaging systems is the detector modules—an area of research and development with a long, rich history. Development of detectors for PET has often seen the migration of technologies, originally developed for high energy physics experiments, into prototype PET detectors. Of the many areas explored, some detector designs go on to be incorporated into prototype scanner systems and a few of these may go on to be seen in commercial scanners. There has been a steady, often very diverse development of prototype detectors, and the pace has accelerated with the increased use of PET in clinical studies (currently driven by PET/CT scanners) and the rapid proliferation of pre-clinical PET scanners for academic and commercial research applications. Most of these efforts are focused on scintillator-based detectors, although various alternatives continue to be considered. For example, wire chambers have been investigated many times over the years and more recently various solid-state devices have appeared in PET detector designs for very high spatial resolution applications. But even with scintillators, there have been a wide variety of designs and solutions investigated as developers search for solutions that offer very high spatial resolution, fast timing, high sensitivity and are yet cost effective. In this review, we will explore some of the recent developments in the quest for better PET detector technology.

Positron emission tomography—some of the basics

The use of imaging to obtain metabolic information is an area of considerable development. In a review article by Cherry (2004), the many different imaging modalities being utilized to acquire such data and the kinds of studies being performed are well summarized. One of the major imaging modalities discussed in Cherry’s review article is positron emission tomography (PET), which utilizes radioactive tracers to make images of the distribution of labeled molecules in vivo. Since it is a tracer approach, the introduction of the labeled compound used for imaging does not alter the biological state of the subject. As discussed in Cherry (2004), the notion of in vivo molecular and genomic imaging (or metabolic imaging as used in this review) is not new, but is rapidly growing for both human and pre-clinical applications. Here, we will focus on the basic detector technologies being applied to the area of PET and begin with a brief review of some of the basic concepts of PET imaging.

Excluding single photon emission tomography (SPECT) systems with 511 keV rated parallel hole collimators, all PET systems utilize coincidence detection of the annihilation photons from positron decay (figure 1). Since the paired gamma rays from the annihilation of the positron are anti-parallel, the detection of the gamma rays determines a line of response (LOR) along which the annihilation took place. Typically, scanners are designed around rings of detectors (figure 2). The coincidence is determined by imposing an acceptance window on the time difference between detection of the two events (typically, a window a few nanoseconds wide). There are three major types of events that can cause coincidence events (figure 3). A true coincidence occurs when a single positron annihilates and both of the gamma rays are detected without either of them scattering in the object being scanned. A scattered coincidence occurs when one or both of the 511 keV annihilation photons from a single positron decay scatter in the object being scanned. This is a true coincidence since it came from single positron annihilation, but the resulting LOR will be misplaced, resulting in reduction of spatial resolution and image contrast. The third type of coincidence is a random, or accidental, coincidence. This event occurs when two positrons annihilate and one gamma ray from each annihilation is detected. If the two events occur close enough in time, then the tomograph electronics will register the event as a coincidence (or a prompt event). However, such randoms are distributed uniformly in time (figure 1) and only the portion of random events included in the prompt window ‘contaminates’ the primary data set.

Figure 1

Coincidence detection. When a positron annihilates, two 511 keV gamma rays are emitted. If both gamma rays interact in a pair of detectors, the events will be recorded at approximately the same time, in coincidence. If a time spectrum is generated (the...

Figure 2

Typical PET system geometry. The object being scanned is usually surrounded by rings of detectors. Each detector acts as a single event detector (records single gamma-ray events). When single events occur within a short time of each other, they are considered...

Figure 3

The three types of coincidence events. True coincidences occur when the two gamma rays from single positron annihilation are detected and neither gamma ray undergoes Compton scattering. Scatter coincidences are true coincidences where one or both of the...

The prompt data provide a data set consisting of LORs along which it is assumed that an annihilation event occurred. If the system has very high timing resolution (<1 ns), then we have some information about where along the LOR the event occurred. For example, a time-of-flight (TOF) resolution of 600 ps allows us to constrain the annihilation position to a roughly Gaussian-shaped region along the LOR with a full width at half-maximum of ~9 cm. If all other parameters are equal, a TOF system can provide a higher signal-to-noise ratio (SNR) than a non-TOF system. In addition to determining if two events are within a timing acceptance window (if they are prompt or delayed events), the system electronics also apply energy acceptance criteria, reducing some of the scattered events. During image reconstruction, corrections can be made for both randoms and scatter. Randoms can be determined in a scanner by either sampling the randoms for each LOR by applying a timing acceptance window that is delayed from the prompt window (accepting events that are not true coincidence or scattered events—figure 1) or measuring the single event rates in each detector unit and calculating the randoms. Scatter can be estimated by one of several approaches. Details on PET image reconstruction can be easily found in the literature and text books (e.g. Lewellen and Karp (2004)).

The role of the detector is to stop an emitted gamma ray and produce a signal that the downstream electronics can utilize. This signal needs to carry information about how much energy was deposited in the detector, information about the physical location of the event in the detector array and timing information (when did the event occur). To fulfill this role, the ideal detector would

  1. have high stopping power (high probability that a 511 keV gamma ray will be totally absorbed by the detector),

  2. have high spatial resolution (ability to determine the interaction location of the gamma ray in the detector to a small spatial volume),

  3. have very good energy resolution (to reject scattered events),

  4. have very high timing resolution and

  5. be inexpensive to produce.

To further complicate an ideal detector wish list, the attributes that should be emphasized depend on the final target application. For simplicity, I like to group applications into three main categories: (1) pre-clinical (e.g. small animal) imaging, (2) dedicated neuro-imaging and (3) whole-body imaging. For all three categories, there is a need to strive for higher sensitivity. For pre-clinical imaging, the parameter getting the most attention is spatial resolution with newer designs striving to achieve <1 mm3 isotropic volume resolution. Timing resolution is not as important for pre-clinical imaging since the count rates are such that randoms rates are generally low and there is no current expectation of obtaining TOF resolution sufficiently low to be useful for such small volumes. For neuro-imaging, the spatial resolution target can be relaxed (perhaps ~3 mm3). Many brain imaging applications want both high spatial resolution and good dynamic imaging, and thus sensitivity is particularly important. Sensitivity is also important if effects such as non-collinearity of annihilation gamma rays are to be addressed in the image reconstruction algorithm. Further, timing resolution is important for randoms’ rejection. To date, TOF is not a big issue in such dedicated systems given the size of the objects being scanned and the current state-of-the-art in TOF capabilities (current TOF resolution does not improve the SNR significantly in brain-sized objects). Body imaging is a bit different story. Current commercial systems have crystal cross section dimensions as small as 4 × 4 mm2. The realities of dose limitations and the number of gamma rays detected (and hence the variance due to the photon statistics) make going to smaller detector dimensions for body imaging generally not a goal in most detector designs. Current systems cannot generally reach their spatial resolution limits due to the amount of variance in the images. Thus, one of the current major areas of development for body systems is that of pushing the TOF resolution to achieve an improvement in the final signal-to-noise ratio in the images. It should be no surprise that the ideal detector is yet to be developed for any of these applications and is an area of very active research.

The goal of this review is to highlight many of these areas of development. Given the many laboratories and companies engaged in such development, I will not attempt to catalog all of the many investigators and innovative solutions connected with PET detectors. It is a very active and innovative area of research. Instead, I have chosen specific examples of many of the different technologies and solutions being worked on. Even then, the technologies I have included are not all inclusive but represent those that, in my opinion, have received the most interest in the literature and recent scientific meetings. I apologize in advance to those whose efforts are not directly discussed and those technologies not included. In preparing this narrative, I have assumed that the reader is familiar with detector technologies, but might not be working with PET imaging systems. Perhaps the best place to begin is with scintillator-based systems since this is what leads to the development of the first PET scanners and is still the dominant detector type used in both commercial and laboratory-built systems.

Basic scintillator detectors

Early PET scanners were based on the use of NaI(Tl) scintillators coupled to photomultiplier tubes (PMT). Scintillator-based designs remain the most common primary detector for PET scanners and are a very active area of development. Table 1 lists some of the properties of scintillators found in many PET detector designs. Ideally, a scintillator should be fast, dense, have high light output and be cheap to produce. There are no ideal scintillators as of yet, but it is another active area of development. While table 1 is not exhaustive, it does illustrate the advantages and disadvantages of scintillators found in current detector designs.

Table 1

Properties of some scintillators used in PET detectors. Note that some of these specifications are subject to change as developers change dopants and trace elements in the scintillator growth. For example, the light output, peak wavelength, decay time...

The first systems were single crystals mounted on individual PMTs. However, as designs began to reduce the crystal cross section size (to obtain improved spatial resolution), alternatives to the one-on-one coupling (one crystal per PMT) had to be developed. The limitations of continuing the one-on-one schemes included the cost (numbers of PMTs), the lack of very small PMTs, the number of electronics channels and problems of packing the PMTs together, given the physical size of the glass envelopes of the PMTs. Two basic schemes emerged (figure 4). One was developed as a consequence of the use of BGO in PET scanners. As seen in table 1, BGO has the advantages of high stopping power but has low light output compared to NaI(Tl). To offset the light output, ‘block designs’ were developed. These designs were either an array of individual crystals or a block of scintillator with saw cuts (Casey and Nutt 1986, Wong et al 1993). In either case, the goal was to channel the light to produce the desired light response function (LRF) to allow decoding of the array with a small number of PMTs. Typically, four PMTs are placed over the crystals in a rectangular pattern and ratios are formed from the PMT signals to provide a transverse and axial position signal as indicated in figure 5. An example of plotting numbers of events versus the transverse and axial position values is shown in figure 6—a map of the response from a 6×6 array of BGO crystals viewed by four PMTs. This basic approach is still being utilized, but the details have evolved in a variety of design approaches as will be discussed below.

Figure 4

The major crystal/PMT decoding geometry options developed for current PET scanners. The continuous ‘Anger’ approach uses an array of PMTs to decode a large, continuous crystal or array of crystals. The ‘block detector’...

Figure 5

Example of a typical block detector decoding scheme. In this case, a 6 × 6 array is decoded by four PMTs. The PMT signals (A, B, C and D) are summed to determine the energy signal (E) as well as two position-dependent signals (X, Z). X and Z are...

Figure 6

A block map of a 6 × 6 BGO crystal array viewed by four PMTs. The relative light output is indicated by the iso-count curves. The heavy white lines indicate what regions will be assigned to each crystal position.

The other major decoding scheme was to adapt the concept of large crystals viewed by an array of PMTs (the Anger camera approach) (Karp et al 1990, 2003, Muehllehner and Karp 1986). This approach was originally developed for use with NaI(Tl) scintillators since it requires good light output. However, the low stopping power of NaI(Tl) was a disadvantage and is no longer in use for PET detector development for human systems, although at least one group is looking at NaI(Tl) detectors for small animal PET applications. The basic concepts developed have continued to be utilized. For example, one commercial PET scanner developed a panel consisting of an array of GSO crystals coupled to a light pipe and then viewed by an array of PMTs in a manner very similar to the NaI(Tl) panels developed for PET imaging (Surti and Karp 2004) and then extended that design to use LYSO in place of GSO (Surti et al 2007).

Let us consider some of the complications with either design in terms of gamma-ray interactions within the detector. Figure 7 depicts some of the possible issues that can degrade the spatial resolution in a pixelated detector. Ideally, the spatial resolution of the scanner should be related to approximately half the width of an individual crystal. There are several compounding factors that limit the resolution achievable in a real detector module. One is scatter within the detector. In most scintillators used in PET, the probability of a photoelectric interaction (total absorption) on the first interaction is less than 50% for 511 keV gamma rays (e.g. BGO is 43%, LSO is 34% and NaI(Tl) is 18%). If the crystals are very short, a gamma ray that Compton scatters in the detector will often exit without a second interaction. This design option has been used in some pre-clinical PET scanners to reduce resolution loss (e.g. Miyaoka et al (2005)). Otherwise, to improve overall sensitivity, the crystals are usually long enough so that a large number of gamma rays that Compton scatter within the detector array are subsequently absorbed with a photoelectric interaction. The resulting multiple points of light emission within the array lead to errors in decoding the crystals (the assignment of map partitions as illustrated in figure 6 assumes a single point of light emission). Investigators have developed some solutions for this problem including estimates of the detector response function in the image reconstruction process (e.g. Alessio et al (2004), Staelens et al (2004), Lee et al (2004), Derenzo (1986)). Penetration of a gamma ray from an oblique angle into the array can also lead to loss of spatial resolution. Most systems assign a LOR based on some small distance from the front edge of a given crystal that has been assigned as the source of the event. The fact that the event occurred deeper in the crystal means that the LOR is not correct (figure 7). One solution that is a very active area of development is to add the ability to determine how deep in the crystal an event actually occurs (depth of interaction or DOI). Further discussion of DOI is included later in this review. The general problem of statistical uncertainty due to the amount of light emitted in the gamma-ray interactions and collected by the photosensors further degrades the crystal identification process. In linear positioning algorithms such as those presented in figures 5 and ​6, the result is a loss of resolution due to misidentification of the crystal. However, an alternative is to use statistical estimation techniques to determine the point of interaction in the crystal array rather than linear algorithms (Parra and Barrett 1998, Budinger 1998, Joung et al 2000, 2002). These techniques can include models for all of the degradation effects and can be implemented at either the detector level or as part of the image reconstruction process.

Figure 7

Some of the factors that degrade spatial resolution in a pixelated detector array. Photon penetration results in incorrect assignment of the line of response in a scanner if depth-of-interaction (DOI) is not able to be measured. Scatter within the detector...

In a similar manner, there are several basic detection issues that degrade the spatial resolution achievable with a continuous (i.e. non-pixelated) detector. The light response function can vary in shape dramatically for all three of the reasons shown in figure 8. The statistics of light emission and collection as well as the distortions in the shape of the LRF due to reflections from the sides and back surfaces work to degrade decoding of the position of the event. Just as in the pixelated case, scatter within the detector will also distort the LRF. And, again as in the pixelated detector, depth of interaction plays a role. The shape of the LRF can depend strongly on the DOI if the crystal is thick enough and the light reflections are controlled (e.g. putting a black absorber on the edges to enhance the effect or using a properly designed retroreflector). This DOI dependence on the LRF can be either a hindrance (degrade spatial resolution) or an added parameter to extract and utilize to address the parallax problem. We will discuss DOI approaches in more detail later in this review.

Figure 8

Some of the factors that degrade spatial resolution in a continuous crystal detector design. The light output and response function reflect both light sampling and the statistics of the amount of light collected. Scatter in the detector and depth of interaction...

What are the major areas of development for scintillator-based detector designs? One obvious area is the search for new scintillators. The advent of LSO (followed by other lutetium compounds—see table 1 for a few of them) provided designers with a series of fast, bright, dense scintillators. These lutetium scintillators are more expensive than NaI(Tl) or BGO to grow, but they do offer performance options unobtainable with the older scintillators. One option is to use the increased light output to decode more crystals per PMT. We have already mentioned the panels of GSO and LYSO crystals using an array of PMTs. Even with the ‘conventional’ block design with four effective PMTs, the use of these brighter scintillators has allowed the use of larger arrays of smaller crystals (e.g. Jakoby et al (2007)). An illustration of the impact of these brighter scintillators on block decoding capabilities is provided in figure 9, which depicts a profile through a block map like figure 6 for four different scintillators. The ‘peak-to-valley’ ratio improves with more light, and this allows decoding more crystals per block.

Figure 9

Illustration of the impact of brighter scintillators on the ability to decode crystals in a block. The profiles are through block maps similar to figure 6 and illustrate that the best decoding (valleys between the peaks closer to 0) is achieved with the...

With the advent of improved electronics and better performance PMTs, time-of-flight (TOF) PET has also reappeared as a commercial product (Surti et al 2007). Modern commercial TOF systems run with timing resolutions of 580—700 ps, about the same timing as achieved in the 1980s with BaF2 (Lewellen 1998), but now with much higher stopping power so that the TOF can be utilized without a major compromise of overall detection efficiency (Moses 2007). Work is also being done on yet brighter scintillators such as LaBr3 and CeBr3 as options for TOF PET with timing resolutions on the order of 300 ps or better (Kuhn et al 2006, Glodo et al 2007). Here, the challenge is the lower stopping power of the scintillator which will require longer crystals to achieve the same overall stopping power of the LYSO and LSO scanners. Longer crystals for a fixed ring diameter usually mean degradation of the radial resolution due to parallax. The fix for this problem brings us back to DOI and how to measure it. But first, we should take a look at changes in photodetectors since the new devices now available to investigators play a large role in the ability to make practical DOI detectors.

Advances in photodetectors

The ‘workhorse’ for scintillation detectors has been the photomultiplier tube (PMT). They are very high gain (typically ~×106), low noise, fast response and relatively low cost. They normally require bias voltages of 800–1200 volts. As noted earlier, the size of the PMTs and their glass envelopes limited the ability to make high spatial resolution detectors with one-on-one (one crystal to one PMT) coupling and lead to the evolution of the block detector concept (figure 4). PMTs have improved and, in particular, have made progress on the problems of the bulk of the devices (the overall size and the size of the vacuum envelope around the photocathode). The detector designer can now utilize PMTs that have multiple dynode chains or channels (effectively multiple PMTs) inside a common vacuum envelope. Photonis, Burle and Hamamatsu, for example, offer PMTs with 4–256 channels in a variety of package sizes. Each channel provides an essentially independent photodetector. For example, the Hamamatsu H9500 provides 256 channels, each one providing a 3 × 3 mm2 photodetector. A few of the many multichannel PMTs are shown in figure 10. These multi-channel PMTs as well as a series of other position-sensitive PMTs (e.g. cross anode) opened up new possibilities in detector designs. Now, designers have less dead space between the active areas of the PMT channels, the ability of using a larger number of channels to view either an array of crystals or a crystal ‘slab’, and more compact PMT housing overall. We will see how some of these PMTs have been used in several detector designs. Work continues to improve PMTs, particularly in terms of the timing response for TOF applications with timing resolutions on the order of 300 ps with LSO, and LaBr3 being reported for PET style detectors (Kuhn et al 2006, Glodo et al 2007, Moses 2007).

Figure 10

Examples of multi-channel PMTs from Burle (left), Hamamatsu (center) and Photonis (right). Such devices come in a variety of sizes with 4–256 channels and have been one of the key components in new PET detector designs.

However, much of the current excitement in new detector designs is centered around alternatives to PMTs. The desire to have options for one-on-one coupling with small crystal cross sections, the use of photodetectors on multiple surfaces of a scintillator array or slab, to further reduce the volume of the photosensor and to have detectors that can operate within magnetic fields led investigators to consider various forms of solid-state photodetectors. One of the first to be used were silicon PIN diodes. These devices are compact, have high quantum efficiency (ranging from 10% to 60% for various different devices) and can be used in magnetic fields. The downside is that they have a gain of unity, requiring low noise preamplifiers for signal readout, and they do not have good enough time resolution for TOF PET applications nor many other PET applications where random rates are a concern. The devices continue to improve with lower dark currents and the availability of fine pitched arrays (figure 11) and have been used in some PET detector designs (Frach et al 2004, Safavi-Naeini et al 2007), but the timing resolution and SNR concerns have limited their use by most developers of PET systems. There have been hybrid designs that use both PIN diodes and photomultiplier tubes as we shall see when we discuss depth of interaction.

Figure 11

Example of a silicon PIN diode 64n element array from RMD with 2.45 mm pixels. Such devices are stable, have become cost effective and can be used in magnetic fields. However, they have unity gain and require low noise preamplifiers and do not have good...

Another device that has been popular is the avalanche photodiode (APD). This device has the advantage over silicon PIN diodes that it has higher gain (~102–103) and has faster timing (on the order of 1 ns), well suited for non-TOF PET applications. While the gain is higher than PIN diodes, APDs still require preamplifiers, but they do not need to be ultralow noise devices. The gain is a function of the bias voltage (as in a PMT) and can run from 200 to 2000 volts depending on the device and the application. One complication is that the gain is temperature sensitive, and gain shifts can be significant with a 1 or 2° (Celsius) change. As with PIN diodes APDs can be found in single packages or as arrays (figure 12). More recently, position-sensitive planar APDs (PS-APD) have been produced, again in a variety of sizes (figure 12). These devices use signals from the four corners of the device to determine the position of an event much like the technique used in the original block detectors. One advantage of PS-APDs is that the number of data channels in the electronics (preamplifiers, ADCs, etc) are reduced compared to discrete arrays in most applications. However, the noise in PS-APDs is generally worse than standard APDs. In general, the SNR of current photodiode devices is a potential limitation that has to be considered for any PET detector design.

Figure 12

Examples of avalanche photodiodes. Devices are available in a variety of packages in both single pixel and arrays and are being used in several different PET detector designs. Such devices have reasonable gains (~102) and can be used in magnetic fields....

One of the newest devices to appear is the Geiger-mode APD or silicon photomultiplier (Herbert et al 2007, Britvitch et al 2007, McElroy et al 2007, Musienko et al 2007, Renker 2007). These devices are under active development by a number of companies (more than nine were known to the author when this review was prepared) and go by a variety of names. For this review, we will term all such devices SiPMs. Because there is such rapid development going on for these devices, performance parameters are also rapidly changing. The basic device is built around a series of APD micro-cells (figure 13). Each cell is an independent Geiger-mode detector that is biased such that when a photon interacts in the cell, it discharges. Each cell is connected to the bias voltage by an individual resistor which results in the cell discharge being quenched. The cells are then connected to a common output. Ideally, a SiPM should produce the same size and shape of a current pulse from each cell, making each cell a digital detection device (on or off). When coupled to a scintillator, the hundreds to thousands of light photons that are emitted by the scintillator interact in a large number of the micro-cells. The result is an output pulse similar to that seen in PMTs (figure 14). Gains for typical devices range from 105 to 107 and produce a signal of several millivolts on a 50 ohm load, making it very similar to PMTs. The bias voltage required again varies between different devices, but typically is between 30 and 150 volts. The devices should be capable of good to excellent timing resolution. Time resolutions of 100 ps have been obtained for single photons, and current devices tested with scintillators have achieved 250–500 ps timing resolution with high light output, fast scintillators. Like most semiconductor devices, SiPMs are susceptible to thermal noise. The dark count rate is due to thermally generated electrons, leading to an avalanche process in one of the micro-cells. While the dark count rate can be quite high (i.e. >300 kHz mm−2), because each micro-cell is independent and the noise pulse is very narrow (i.e. <50 ns), the noise signal can be removed by using a simple low energy threshold. Furthermore, for lutetium-based scintillators with signal integration times on the order of 100 ns, noise due to dark counts during signal integration will rarely exceed the single photoelectron value. One limitation of such devices is that current designs have relatively large capacitance per unit area, which leads to problems in scaling up such devices to large detection areas (greater than about 3 × 3 mm2 active areas in a single device with current devices).

Figure 13

Example of a Geiger-mode avalanche photodiode or silicon photomultiplier (SiPM). A series of micro APD cells are connected via the resistance in the surrounding silicon. When a photon interacts in a cell, it discharges and is then quenched via the resistance...

Figure 14

Example of a 511 keV energy spectrum from a SiPM. In this case, a Zecotek MAPD device was coupled to a 2 × 2 × 12 mm3 LFS-3 crystal (also from Zecotek).

Note in figure 13 that there is a dead space (as far as light photon detection) around each micro-cell. The dead space is needed to reduce optical cross talk between micro-cells (preventing neighboring cells from firing spuriously) and to provide electrical isolation between the micro-cells (except for the needed resistance and capacitance to provide the needed quenching of the current once a cell fires). The amount of dead space required depends on the details of the device design (Sadygov et al 2006). To achieve good linearity, a large number of micro-cells is needed to avoid saturation of the device and multiple interactions in the same micro-cell. However, as the numbers of cells increase by using smaller diameter micro-cells, the amount of dead space also increases, reducing the overall quantum efficiency over the area of the device. So, detector designers have to choose between linearity and quantum efficiency and determine the best compromise for their detector. Current devices can offer a range of microcells mm−2 of 100–10 000. These devices can also be made as either p-silicon on an n-substrate (and are sensitive to blue light) or as n-silicon on a p-substrate (and are sensitive to green light). In either case, the process is typically a CMOS fabrication technique much as is used for integrated circuits and offers the possibility of low cost as well as a wide range of pixel sizes (the collection of micro-cells within a device) and fine pitch arrays. The devices can be used in magnetic fields and samples have been tested in fields up to 15 Tesla without showing any degradation. Manufacturers are currently offering both p–n and n–p variants, and some arrays have also been shown. The interest in these devices is high, and many groups are experimenting with detector designs built around the expectation of low cost, high performance SiPM arrays.

DOI approaches with scintillators

A major area of development for the last several years has been that of measuring DOI in order to reduce parallax errors and maintain good sensitivity. DOI also allows smaller detector rings since the parallax corrections reduce the loss of spatial resolution with the radial position. The initial work has been primarily in scintillator-based detectors although we will return to DOI when we discuss recent progress in non-scintillator-based detector systems. One of the first approaches was to utilize pulse shape discrimination (PSD) and scintillators with different decay times, as illustrated in figure 15 (Costa et al 1986, Dahlbom et al 1997, Seidel et al 1999, Chung et al 2004, Braem et al 2007a, Yong et al 2007). Since the different decay times result in different pulse shapes, once can use several techniques to discriminate which scintillator is the source of the light (e.g. simply setting different windows in the time spectra). Many combinations of scintillators have been used including LSO, GSO, LuAP and BGO. Some designs use different dopants in LSO, LYSO or GSO to change the decay times enough to utilize the pulse shape discrimination approach. The PSD approach has been used in several scanners—both University and commercially developed. Most applications of such detectors have been for pre-clinical and dedicated brain scanners.

Figure 15

Pulse shape discrimination based depth-of-interaction detector designs. The concept is to use two or more layers of crystals that have different light decay times. A pulse shape analyzer can then be used to separate the light from each layer based on...

One of the drawbacks with the PSD approach has been the cost of fabrication of the detector modules due to the mechanical operations to assemble and match the needed crystals. Another complication is limitations on the system timing resolution parameters due to the different decay times, making such approaches not attractive for TOF detector designs. Another problem (in fact, a general problem for most DOI designs) is that of scattering within the crystal array resulting in light being produced from more than one of the scintillators, reducing the DOI accuracy. As we will discuss a bit later, there are methods being developed to improve the estimation of the point of first interaction by including appropriate knowledge about the physics of such multiple interactions and the geometry of the scanning system.

Recently, investigators at UC Davis reported on a new approach to the PSD detector design at the 2007 IEEE Nuclear Science Symposium and Medical Imaging Conference. Rather than combining layers of crystals with different decay times, they coated the lower portion of crystals with phosphors designed to absorb light and then re-emit it, with the expected changed effective decay time. Since the amount of light that reaches the phosphor is dependent on the DOI of the event, the resulting PSD becomes depth dependent as well. Another approach presented by Eriksson at the 2007 IEEE Nuclear Science Symposium and Medical Imaging Conference is to use wave shifting sheets between layers of crystals to shift light from one crystal that would be absorbed by the next crystal in the stack (e.g. LaBr3 followed by LSO). While still in the early stages of development with many technical issues to explore, these new investigations illustrate that even ‘old’ techniques such as PSD can often be modified and improved upon for modern PET detector designs. Most of these new efforts are focused on pre-clinical imaging systems and dedicated head scanners since they are not ideal for TOF applications.

Three other approaches currently being used for DOI scintillator-based detector designs are illustrated in figure 16. One approach is to stack two or more layers of crystals such that each layer is offset from the one below it. In doing so, and if there is enough light collected, each crystal can be identified in the crystal map block. The original work was with two layers, but investigators have succeeded in applying this approach to four layers (Nishikido et al 2007, Hasegawa et al 2007) and even combined the offset approach with a PSD approach to decode up to eight layers of crystals (Inadama et al 2007). This is another case where high light output scintillators are needed to improve the peak-to-valley ratios in the crystal maps. As more layers are added, the crystal centers are moved closer together and misidentification becomes more likely.

Figure 16

Three of the many approaches for DOI detector modules. On the left, layers of crystals are offset so that each crystal is visible in a crystal map similar to that of figure 5. In the center, light is collected from each end of a crystal and the ratio...

Another approach that has been under considerable development is to place photosensors at both ends of a crystal and then use the ratio of the light from the two sensors to determine the DOI. Investigators using this technique are able to achieve DOI resolutions on the order of 2 mm for single crystals, illuminated with a collimated beam of photons from the side. This concept has undergone several changes as improvements in photosensors have occurred. One of the first designs was developed at Lawrence Berkeley Laboratories (Moses et al 1993). This design measures light emitted from both ends of an array of crystals. One end is a PMT for fast timing and the other is a PIN diode array for crystal identification as well as measuring light output (figure 17). This detector module design is used in a pre-clinical scanner built by the same group. Another version of this approach is to use a pair of PS-APDs (figure 18—an example from UC Davis and RMD) (Yang et al 2006). Initial tests of single crystals viewed with SiPMs at both ends have also been reported with similar results (Shao et al 2007).

Figure 17

An early DOI module based on a double-ended readout developed at Lawrence Berkley Laboratory. The crystal array is viewed by a PMT at one end (fast timing and light output measurement) and a PIN diode array at the other end (crystal identification and...

Figure 18

Another DOI approach is to use PS-APDs at both ends of a crystal array. The DOI is determined by the ratio of the light detected at both ends of the array. The DOI data shown are from such a detector developed at UC Davis.

The third approach illustrated in figure 16 is to share the light between paired crystals such that the amount of light shared between crystals reflects the DOI (Lewellen et al 2004, Miyaoka et al 1997). By taking the ratio of the light from the two crystals, DOI information is provided. The DOI accuracy achieved to date (~3–4 mm) is not as good as the double-ended approach, but has the advantage of using a single-ended readout, simplifying the detector module. This approach would best be implemented in full detector modules with one-on-one decoding; until recently, it was not a realistic option. But with expectations for fine pitch SiPM arrays, this approach may well become practical.

Other options for detector designs have become available with the new photonics devices. One approach developed by the Stanford group uses PS-APDs with discrete crystals mounted in layers on the PS-APDs as illustrated in figure 19 (Foudray et al 2007). Conventional mapping techniques are used to determine the crystal producing the light. This then is a layered DOI design with discrete crystals. Since the long side of the crystal is mounted on the PS-APD, more light is collected, which is an advantage of this approach. The technique also requires very thin circuit boards and PS-APD elements to maintain a good packing fraction (the percentage of the detector face that is a useful crystal). Another variation of a layered crystal DOI design is one that will use sheets of wavelength shifting fibers to provide crystal identification and then use any one of a number of photodetectors at the ends of the crystal stacks for energy and timing—figure 20 (Braem et al 2007b). A challenge for this design approach (and indeed for any design with many discrete channels required) is the cost and density of the photodetectors. The current development team for the device of figure 20 is planning on using SiPM arrays to provide the needed number of photosensor channels at an affordable cost. While not in the scope of this review, it is also important to note that advances in electronics also allow realistic systems to be considered with a very large number of channels. With modern field-programmable gate arrays and serial analog-to-digital converters, both increasing in power and speed while becoming more affordable, it is well within reality to consider detectors requiring hundreds of digital channels. But even with such advances in these devices, the designer still have to contend with power requirements, heat generation, and the number and density of interconnections to implement large numbers of channels.

Figure 19

A DOI detector design using PS-APDs and crystals mounted in layers on the PS-APD surface. One of the advantages of such an approach is having a large number of the events that scatter within the crystal array and interact in different discrete crystals...

Figure 20

A DOI design that uses stacks of crystals layered with wavelength shifting fibers that are used to determine crystal identification (Braem et al 2007b). Such designs require a large number of photodetectors, which has now become feasible with the recent...

DOI is not restricted to arrays of discrete crystals. Due to the cost of making arrays of individual crystals, several groups have taken a second look at using crystal slabs. Two basic approaches have emerged. One uses layers of crystal slabs to achieve DOI (figure 21). In this design approach, PS-APDs or arrays of APDs or SiPMs are used to view each layer and decode the event position, determine timing and measure the amount of energy deposited. A challenge in using crystal slabs is the ability to extract good event positioning near the edges of the crystal. But here again, advances in photosensors and support electronics allow one to consider photosensors that can sample the light response functions (LRF) in each slab finely enough to determine changes in the LRF as the light source nears the edges of the crystal. This capability coupled with statistical position estimators (Joung et al 2000, Ling et al 2007, van der Laan et al 2007, Milster et al 1985, Bruyndonckx et al 2007, Tavernier et al 2005) (which can be implemented in modern FPGAs) makes such designs practical to consider. An advantage of this approach is that most events that Compton scatter within the detector (and then undergo photoelectric interactions within the detector) will produce events in more than one layer, providing additional information to determine the point of first interaction. These detectors should also have reduced cost compared to layers of discrete crystals since there is less effort required in cutting and polishing the scintillators. A disadvantage is the many readout channels required.

Figure 21

A DOI design using layers of crystal slabs. APD arrays are used to read out light from each layer. The light distribution in each slab is used to determine the point of interaction in each slab. The layers are also useful in handling events that scatter...

Another approach is to use single- or double-ended readout designs of thick crystal slabs and extract the DOI information from the changes in the measured LRFs (Hunter et al 2007, Ling et al 2007, Moore et al 2007). The double-ended approach is essentially the same as used for the double-ended DOI scheme for discrete crystals. One or both of the two photosensors (in this case, arrays or position-sensitive devices) are used to determine the event position centroid as projected on the surface of the detector, and the ratio of the signals between the two arrays is used to determine the DOI. However, one can also derive DOI information from a single-ended readout of a thick crystal slab (Moore et al 2007, Ling et al 2007, Lerche et al 2005) with appropriate modeling of the shape of the LRF as a function of depth. Figure 22 shows one such device being developed for a pre-clinical PET scanner. Currently, the module uses a 50 × 50 × 8 mm3 crystal coupled to a 64-anode PMT to provide <1.4 mm spatial resolution and 1 bit of DOI data. As the technology is improved, one can expect to see much thicker crystals (15–25 mm thick) with improved DOI resolution as well. In both of these approaches, the continued development of fast, bright scintillators is important. Since measurement of the LRF shape is needed, a bright scintillator is needed to reduce the noise of the measured shape (improves the accuracy of the estimation). At the same time, fast scintillators reduce the light pile-up that can distort the LRF shapes with a count rate if the detector needs to handle high single rates. While there have been various approaches to correcting the integrated pulse data from pile-up effects, most PET scanners have not implemented sophisticated pile-up correction schemes. In the future, we are likely to see such schemes used to further reduce variance in LRF shapes for detector designs with statistical position estimators such as those being used in these DOI approaches. One disadvantage of the thick slab approach (as compared to the layered slab designs) is that there is less information from the detector about events that scatter within the detector. Thus, scatter for the thick slab designs must be incorporated in the estimation technique. At the same time, the thick slab approach has the advantage of providing a continuous DOI readout that can be improved as scintillators become brighter. The thick slab designs also have the advantage of requiring fewer electronic channels.

Figure 22

An example of an assembled ‘thick’ slab detector. A 50 × 50 × 8 mm3 LYSO crystal is mounted on a 64-anode PMT. A statistical estimation approach is used to provide a spatial resolution of ~1.4 mm and at least 2 bits of...

Alternate designs

While scintillator-based detector designs still dominate scanner designs, there has always been interest in alternative approaches—generally with the goal of being able to offer much higher spatial resolution with higher packing fractions. Of course, the hope is always that such goals can be accomplished at a lower cost than very high volume resolution (<1 mm3) scintillator-based designs. Some of the earliest such alternatives were based on wire chambers, and one such design has been commercialized for pre-clinical applications and at least one research system has been used for human studies (Ott 1993, Schafers et al 2005). These devices have included chambers with converter foils as well as photosensitive chambers that view scintillators. While there is still ongoing work with wire chambers, most efforts recently on non-scintillator solutions have been focused elsewhere. In a similar fashion, there have been several efforts to utilize plastic scintillator sheets or layers of scintillation fibers as high spatial resolution detectors. Due to the lack of good photoelectric cross section in these types of detectors, they have not made an impact on human whole-body imaging (lack of sensitivity and scatter rejection) compared to the many other solutions. However, much as wire chambers, these devices are being considered for pre-clinical imaging and a small animal system based on layers of fiber optic scintillation detectors has been built (Tsyganov et al 2006). Another area that is being looked at again is the use of liquid xenon as a scintillator, in particular for TOF PET (Doke et al 2006, Gallin-Martel et al 2006). One of the principal advantages of layered fiber or wire chamber designs is that the multiple interactions within the detector system can be kept essentially independent, making the estimation of the point of first interaction a simpler task. Another variation on this concept is that of the resistive plate chamber (a gaseous particle detector: RPC). At least, one group investigates RPC devices for small animal PET (Blanco et al 2006), and there have been discussions about investigating such technology for human whole-body scanners. Prototype devices have obtained spatial resolutions of 0.5 mm with timing resolutions of ~300 ps. A drawback of these devices is that they do not provide energy resolution, raising concerns about high-scattered event rates and the impact on single rates (randoms). However, all of these systems are generally more complex to fabricate into a scanner than scintillator or solid-state-based solutions. This leads us to consider some of the solid-state solutions being investigated.

One area of development is that of silicon detectors. One challenge in considering silicon detectors is that such devices have a relatively low proton number compared to scintillators commonly used in PET designs. For example, if one compares a 4.8 cm thick Si detector to a 1 cm thick BGO detector, the probability of interactions for Si is ~48% and for BGO ~58%, but the photofraction for Si is <1% compared to ~43% for BGO (Cesca et al 2007). However, silicon detectors generate a large number of charge carriers per unit of energy absorbed and these carriers have good mobility. As a result, these devices generally have very good energy resolution. With a prototype silicon pad detector, Park et al (2007a, 2007c) have obtained ~1% energy resolution for 140.5 keV and a 82 ns timing resolution for 511 keV events. The detector used for this experiment is intended for use as a Compton scatter detector since more than 99% of the first interactions will be Compton at 511 keV. In such a design, the silicon detector is placed in front of a more conventional PET detector array (e.g. BGO or LSO blocks). Those events that interact in the pad detector (and are then absorbed in the conventional PET detector ring) will provide the highest spatial resolution with good sensitivity. For the design presented by Park et al, the simulations for a small animal system indicate spatial resolutions of ~1 mm with a 9% sensitivity (for those events that interact in the pad detector) and about 1.7 mm resolution with a 21% sensitivity for a ring of 3 × 3 × 20 mm3 BGO detectors as the conventional PET detector portion of the system. The prototypes in this work were originally developed for Compton imaging and consisted of detectors with 1.4 mm2 pads arranged in a 32 × 16 array (512 pads). Thus, each detector has 512 channels and this means that a very large number of electronics channels will need to be accommodated for such a system. Since silicon detectors are also of interest in high energy physics, there are ASICs available for reading out such systems. However, if timing resolutions of 82 ns such as were obtained by Park et al are normal for such devices, their application will be restricted to relatively low count rate pre-clinical systems.

Another approach with silicon detectors is to implement a row/column scheme similar in concept to what has been proposed for scintillator detector systems by using silicon strip detectors rather than pad detectors (figure 23). Again, such a detector can be used in a Compton imaging arrangement as is being developed by Park et al. However, another group which focused on pre-clinical imaging has proposed stacking many layers of silicon strip detectors to form the entire detector system (Cesca et al 2007). Due to the low photoelectric fraction, the proposal is to only keep those events that occur once in the stack (discard all multiple events due to scatter within the detector stacks). Such an approach will provide very high spatial resolution, but has not been adopted by a large number of PET detector developers due to the low photofraction and the resulting sensitivity issues when using silicon devices. An alternative that has received increasing attention is the use of CdZnTe (CZT) and CdTe room temperature semiconductors that have been under development for a number of years.

Figure 23

An example of a silicon strip detector. The gamma ray produces hole–electron pairs that are collected on the strips and the charge then directed to appropriate amplifiers. On the right is an example of one such detector used at the University...

CZT and CdTe can be made into large detectors with either pixelated anodes or cross strip anodes for reading out resolution elements. One of the main differences between the two materials is the charge transport properties. The actual differences can favor either CZT or CdTe depending on the growth techniques used and the temperature of the devices. Most investigators have reported that the charge mobility of CdTe is somewhat better (typically on the order of 20%) at room temperature. The devices generally have much better energy resolution than scintillator-based detectors and can be made to provide very good spatial resolution. The challenges for PET have been timing resolution (they generally have low charge mobility) and the low photopeak detection efficiency for 511 keV gamma rays (Arnaud et al 2007, Levin et al 2006, Vaska et al 2006). CZT and CdTe detectors have not been developed for any human PET scanners due to the challenges in timing resolution. However, for pre-clinical imaging, the small size of the subjects and the restrictions on dose administered make the timing challenges of CZT/CdTe less critical. One approach for a CZT detector being developed at Stanford is illustrated in figure 24, and a similar approach has been investigated for CdTe (Hadong et al 2007). The investigators at Stanford use 40 × 40 × 5 mm3 slabs of CZT to provide sufficient stopping power and maintain good spatial resolution. In this design, orthogonal anode and cathode strips are sandwiched between the CZT slabs to provide spatial decoding. The slabs are 5 mm wide along one axis, and the orthogonal strips are used to localize the event within each slab. Such a design will have multiple interactions within the detector, so once again a mechanism to measure DOI is needed. With the orthogonal anode/cathode strips and the limited charge mobility in CZT, an event in the slab is only ‘seen’ by the local strips. Much in the same way as a pad detector, this feature of CZT allows multiple events within the slab to be recorded independently. For this design, the investigators plan to use the kinematics of Compton scattering within the detector array along with the fine pitch readout (ability to separate multiple interactions into separate events) to sort out the most likely point of the first interaction (Chinn et al 2007).

Figure 24

A design for a PET detector based on CZT (left) and an image of a 32 × 32 × 2 mm3 CZT prototype detector element being developed at Stanford University.

One might wonder about high purity germanium (HPGe) detectors—long a mainstay of high resolution x-ray and gamma-ray spectroscopy. As with CdTe and CZT, advances have been made in growing HPGe detectors and the fabrication of HPGe strip detectors is available from commercial suppliers. With the use of orthogonal strips, these detectors can perform DOI determination much as is done with the CdTe and CZT approaches we have already mentioned and there is a project exploring such detectors for use in PET (Cooper et al 2007a, 2007b). In addition, HPGe offers very high energy resolution (better than CZT), useful in determining the point of first interaction both within the detector and in rejected events that scatter in the object being scanned. HPGe is a mature technology and is available in much larger volumes than other materials used in solid-state detectors. HPGe detectors are probably limited to pre-clinical applications due to the timing resolution of HPGe detectors and the complications of keeping the detectors cold enough.

The ongoing challenge of detector design

While this review has presented more material on scintillator designs than on alternative approaches (e.g. solid state detectors, wire chambers, etc), the basic challenges for all design directions are similar. With the increasing use of PET/CT scanners in the clinic and PET and PET/CT scanners for pre-clinical imaging and the preliminary developments in MRI/PET systems, there is ample incentive for instrument designers to continue to push the capabilities of detector systems to attack the goals of improved spatial resolution, improved sensitivity and improved timing. In the beginning, we listed the main attributes for a PET detector. As we have seen, there is no one solution. But it is even more complex than we have discussed thus far. As we push the limits of the detectors, we cannot do so without considering other aspects of the potential scanner system, in particular the image reconstruction algorithm. With the everincreasing computing power available and the continued advancement in statistical image reconstruction algorithms that can include the physical effects of the photon transport in both the object and the detector system, it is not always clear where to put the most development effort in the basic detector. For example, DOI can be done at the detector level based on single events, or the basic distribution of data in the detector array can be passed directly to the reconstruction engine and it can use coincidence information to provide more information on the limits of the direction of the incident gamma rays to aid in the point of first interaction determination within the detector array. Thus, modern detector designers really have to start looking at potential system architectures and image reconstruction capabilities early on their efforts.

Из самолета? - повторила.  - Что происходит. С какой стати университетский профессор… Это не университетские дела. Я позвоню и все объясню. Мне в самом деле пора идти, они связи, обещаю.

Comments

Leave a Reply

Your email address will not be published. Required fields are marked *