Sparse Acquisition Gamma Cameras

20170285191 · 2017-10-05

    Inventors

    Cpc classification

    International classification

    Abstract

    An imaging method and device are described for improving the performance of a gamma camera by optimizing a figure of merit that depends upon cost, efficiency, and spatial resolution. In a modular gamma camera comprising a tiled array of gamma detector modules, the performance figure of merit can be optimized by sparsely placing gamma detector modules within the gamma camera, optimizing collimation, and providing means for detector and/or collimator motion. Sparse gamma cameras can be constructed as flat or curved panels, and elliptical or circular rings.

    Claims

    1. An imaging method for improving a performance figure of merit of a gamma camera, the performance figure of merit comprising at least one variable related to at least one of cost, efficiency, and spatial resolution and the gamma camera comprising a collimator and at least one tiled array of gamma detectors, the collimator having a performance, the method comprising at least one of the steps of: sparsely placing a gamma detector within the gamma camera, optimizing the collimator performance, and moving at least one of the gamma detector and the collimator.

    2. The imaging method of claim 1, further comprising a step of predetermining at least one performance requirement for at least one application, the at least one performance requirement defining an improvement of the performance figure of merit.

    3. The imaging method of claim 1, wherein the at least one tiled array of gamma detectors comprises at least one of a flat panel, a curved panel, a circular ring, and an elliptical ring.

    4. The imaging method of claim 1, wherein the optimizing the collimator performance step comprises using at least one of a parallel hole collimator, a focused hole collimator, a slit-slat collimator, a rotating slat collimator, a multiple pinhole collimator, a coded aperture collimator, and a Compton scatter collimator.

    5. The imaging method of claim 1, wherein the at least one tiled array of gamma detectors comprise tileable gamma photon detectors having an area less than 400 cm.sup.2.

    6. The imaging method of claim 1, wherein the at least one tiled array of gamma detectors comprises at least one of: a scintillator with an optically-coupled photodetector and a semiconductor direct conversion detector.

    7. The imaging method of claim 1, further comprising a step of iteratively searching for an optimal performance figure of merit by: generating a simulated numerical distribution of a gamma emission radioactive material; modeling an imaging system and a data acquisition sequence, the imaging system comprising a means to move the gamma camera; generating a simulated acquisition of imaging data using the simulated numeral distribution, the modeled imaging system, and the modeled data acquisition sequence; iteratively reconstructing the imaging data to obtain an imaged distribution; calculating the performance figure of merit for each iteration; and incrementally changing at least one of the imaging system and the data acquisition sequence.

    8. The imaging method of claim 1, wherein the sparsely placing a gamma detector within the gamma camera step comprises using at least one scintillator gamma camera and at least one semiconductor gamma camera.

    9. The imaging method of claim 1, wherein the sparsely placing a gamma detector within the gamma camera step comprises removing at least one gamma detector from the at least one tiled array.

    10. The imaging method of claim 1, wherein the sparsely placing a gamma detector within the gamma camera step comprises removing at least one column or row of gamma detectors from the at least one tiled array.

    11. The imaging method of claim 1, wherein the sparsely placing a gamma detector within the gamma camera step comprises removing gamma detectors in a checkerboard pattern from the at least one tiled array.

    12. The imaging method of claim 1, wherein the collimator has design parameters and wherein the optimizing collimator performance step comprises adjusting the collimator design parameters to affect at least one of efficiency and spatial resolution.

    13. The imaging method of claim 1, further comprising a step of acquiring imaging data during a period of time, wherein the moving at least one of the gamma detector and the collimator step comprises at least one of moving the gamma detector from a first detector position to a second detector position and moving the collimator from a first collimator position to a second collimator position during the acquiring imaging data step.

    14. The imaging method of claim 1, wherein the performance figure of merit comprises a ratio of efficiency divided by spatial resolution squared divided by gamma camera cost.

    15. A gamma camera with an improved performance figure of merit, the performance figure of merit comprising at least one variable related to at least one of cost, efficiency, and spatial resolution and the gamma camera comprising an electrical power input, a control input, a data output, a collimator, the collimator having a performance, and at least one tiled array of gamma detectors and further comprising at least one of: a sparsely placed gamma detector within the gamma camera, an optimized performance collimator; and a means for at least one of a gamma detector motion and a collimator motion.

    16. The gamma camera of claim 15, further comprising at least one application with at least one predetermined performance requirement defining the improved performance figure of merit.

    17. The gamma camera of claim 15, wherein the at least one tiled array of gamma detectors comprises at least one of a flat panel, a curved panel, a circular ring, and an elliptical ring.

    18. The gamma camera of claim 15, wherein the optimized performance collimator comprises at least one of a parallel hole collimator, a focused hole collimator, a slit-slat collimator, a rotating slat collimator, a multiple pinhole collimator, a coded aperture collimator, and a Compton scatter collimator.

    19. The gamma camera of claim 15, wherein the at least one tiled array of gamma detectors comprise tileable gamma photon detectors having an area less than 400 cm2.

    20. The gamma camera of claim 15, wherein the at least one tiled array of gamma detectors comprise at least one of: a scintillator with an optically-coupled photodetector, and a semiconductor direct conversion detector.

    21. The gamma camera of claim 15, wherein the improved performance figure of merit is found by iteratively searching for an optimal performance figure of merit by: generating a simulated numerical distribution of a gamma emission radioactive material; modeling an imaging system and a data acquisition sequence, the imaging system comprising a means to move the gamma camera; generating a simulated acquisition of imaging data using the simulated numeral distribution, the modeled imaging system, and the modeled data acquisition sequence; iteratively reconstructing the imaging data to obtain an imaged distribution; calculating the performance figure of merit for each iteration; and incrementally changing at least one of the imaging system and the data acquisition sequence.

    22. The gamma camera of claim 15, wherein the sparsely placed gamma detectors within the gamma camera comprises using at least one scintillator gamma camera and at least one semiconductor gamma camera.

    23. The gamma camera of claim 15, wherein the sparsely placed gamma detectors within the gamma camera are removable or can be left out from the at least one tiled array.

    24. The gamma camera of claim 15, wherein the sparsely placed gamma detectors within the gamma camera comprise at least one column or row of gamma detectors removed or left out from the at least one tiled array.

    25. The gamma camera of claim 15, wherein the sparsely placed gamma detectors within the gamma camera are removed or left out in a checkerboard pattern from the at least one tiled array.

    26. The gamma camera of claim 15, wherein the collimator has design parameters and wherein the optimized performance collimator has adjusted collimator design parameters that affect at least one of efficiency and spatial resolution.

    27. The gamma camera of claim 15, the camera capable of acquiring imaging data during a period of time, wherein the means for at least one of a gamma detector motion and a collimator motion comprises means for at least one of moving the gamma detector from a first detector position to a second detector position and moving the collimator from a first collimator position to a second collimator position during a period of acquiring imaging data.

    28. The gamma camera of claim 15, wherein the performance figure of merit comprises a ratio of efficiency divided by spatial resolution squared divided by gamma camera cost.

    Description

    BRIEF DESCRIPTION OF THE DRAWINGS

    [0019] The organization and manner of the structure and operation of the invention, together with further objects and advantages thereof, may best be understood by reference to the following description, taken in connection with the accompanying non-scale drawings, wherein like reference numerals identify like elements in which:

    [0020] FIG. 1 is a pixelated CZT detector module of the preferred embodiment of the invention.

    [0021] FIG. 2 is an illustration of an array of the CZT detector module of FIG. 1 in a size appropriate for a general purpose flat panel SPECT gamma camera.

    [0022] FIG. 3 illustrates a transaxial view of a human body surrounded by two opposed SPECT scanning CZT gamma cameras in either a flat panel or concave curved panel configuration.

    [0023] FIG. 4 is a graph comparing the planar system resolution for scintillator and CZT gamma cameras using registered square-hole collimator concepts.

    [0024] FIG. 5 is a process flow diagram depicting optimization of a performance figure of merit for a gamma camera based on sparely placed gamma camera modules.

    [0025] FIG. 6 depicts multiple possible configurations of flat panel CZT gamma cameras in which sparse placement of detector modules is achieved in various ways.

    [0026] FIG. 7 depicts a pair of sparsely populated flat panel CZT gamma cameras in which detector motion is employed to compensate for the missing detectors.

    [0027] FIG. 8 depicts a circular or elliptical ring of CZT detectors surrounding a human body, similar to a PET scanning configuration.

    DETAILED DESCRIPTION OF THE ILLUSTRATED EMBODIMENTS

    [0028] While the invention may be subject to embodiment in different forms, there are shown in the drawings, and herein will be described in detail, specific embodiments with the understanding that the present invention is to be considered an exemplification of the principles of the invention, and is not intended to limit the invention to that as illustrated and described herein.

    [0029] FIG. 1 shows a modular solid-state direct conversion CZT gamma detector called D-Matrix™. This aggregator module (AM) 20 is composed of four CZT crystal detectors 21, each designated a gamma module (GM), each measuring about 2.2 cm×2.2 cm×0.5 cm thick. Each GM 21 may comprise a monolithic crystal or may comprise multiple tiles of CZT mounted in an abutting composite with small gaps between the tiles. The top GM 21 surface in FIG. 1 is the metallic cathode, which is typically monolithic but which can also be pixelated. Gamma photons generally enter the GM from the cathode side of the CZT crystal. The lower surface of the GMs 21 has metallic anode pixels with a pitch of about two mm in an 11×11 pixel pattern for each CZT detector. The CZT is bonded to a substrate that provides mechanical stability and electrical connections to the read-out electronics boards, which include ASICs, one or more ADCs, and an FPGA. Negative DC bias voltage (about −500 V for 0.5 cm thick CZT) is provided for the cathodes and an optional capacitive decoupling circuit allows reading the induced cathode signal as well as the anode signals. The overall dimension of the D-Matrix™ aggregator module (AM) 20 with 4 GMs 21 (CZT detectors, 0.5 cm thick) and readout electronics is about 4.4 cm×4.4 cm×5 cm thick. This AM 20 can be tiled to comprise various sized arrays to constitute flat panels, curved panels, rings, ellipses, or other structures as will be disclosed in this invention. Some of the features of the AM 20 are listed in Table 1.

    [0030] One of the advantages of CZT gamma detectors is that they occupy a smaller volume than scintillators and PMTs, so that CZT gamma cameras can be shielded by less volume and weight of heavy metal, such as lead (Pb) or tungsten (W). Thus a SPECT system gantry designed specifically for a CZT gamma camera can be smaller and lighter than a conventional SPECT gantry. It may have a smaller footprint and fit into a smaller examination room.

    [0031] Furthermore, the energy resolution of a CZT gamma camera, as listed in Table 2, is much better (about 4% or less compared to about 9.6% for NaI). Practically, this means that a narrower energy window can be used for the photopeak projection images that are reconstructed into 3D SPECT images. This narrower energy window will discriminate against scattered gamma photons which blur the SPECT image. Thus, the final image contrast should significantly improve in a CZT camera compared to a scintillation gamma camera.

    [0032] FIG. 2 depicts one such array 22 of AMs 20 (only two of 108 are labelled) in a size typical of flat panel gamma cameras used in general purpose medical SPECT. The illustrative array 22 comprises nine rows and 12 columns of AMs 20 forming a flat panel gamma camera 22 with dimensions of about 39.6 cm×52.8 cm. This illustrative arrangement comprises nine×12=108 AMs and (nine×22)×(12×22)=198×264=52,272 pixels, as tabulated in Table 2.

    [0033] The original gamma camera built by Hal Anger in 1957 was a flat panel detector with a circular field of view. All commercial clinical gamma cameras have been derivatives of the original flat panel, although the field of view became rectangular in the 1980s when 3D tomographic image reconstruction from SPECT became practical. There have been several attempts to use curved scintillation cameras for brain or heart imaging applications, but none have been commercially successful. Several small-animal preclinical SPECT systems were built using CZT modules in a ring, but these have not been sold in any successful quantities. Nevertheless, modular detector assemblies, such as the Aggregator Module, AM 20, of FIGS. 1 and 2 are amenable to forming a curved arcuate surface, or ring, or ellipse of detectors.

    [0034] FIGS. 3(a) and 3(b) illustrate this design option for a curved detector surface. FIG. 3(a) illustrates a transaxial cross-section (view from the patient's feet) of a dual-detector SPECT system with a human body 30 as the example object being imaged by gamma emission. Both detector arrays 22 in FIG. 3a, as in FIG. 2, are shown in cross-section with 12 columns 24 illustrated by gray boxes. By reference to FIG. 2, it is apparent that each column 24 comprises nine AMs 20 with length about 39.6 cm. Not shown explicitly are collimators, radiation shielding, and readout electronics, but these elements are well known by those of skill in the art. FIG. 3(b) illustrates curved arcuate surface detectors 23 that could be formed by placing the columns 24 of AMs 20 along a curved surface that roughly follows the contour of a typical human body.

    [0035] A gamma camera could be designed to mechanically transform between the flat panel configuration of FIG. 3(a) and the curved panel configuration of FIG. 3(b). It is anticipated that the curvature could be fixed or variable. The reconstruction would be easier for a fixed curvature, but if the position of each detector is tracked mechanically or by other means, then the iterative reconstruction software could accommodate “body contouring” with variable curve panels. As will be apparent to those skilled in the art, the hinge points for such a transformation would depend upon the collimation scheme employed. For example, if parallel hole collimators were attached to each CZT detector column 24, then the hinge-points would be at the front face of the collimator columns. However, if multiple-pinhole collimation were employed, the hinge points would be at the front face of the CZT detectors 21 of FIG. 1.

    [0036] The smaller square pixels of a CZT gamma camera can also be used to advantage compared to the larger overlapping Gaussian pixels of a scintillation gamma camera. Collimators can be designed (not the subject of the present invention) to optimize the spatial resolution and efficiency of the detector-collimator system, again improving the image contrast. By way of illustration, Table 3 shows a comparison of typical hexagonal parallel-hole collimators, as used in most clinical SPECT systems, and some possible pixel-registered square-hole collimators. These design concepts have not been optimized, but are illustrative of the performance achievable with such pixelated CZT specific collimators.

    TABLE-US-00003 TABLE 3 Comparison of typical SPECT collimators and potential CZT-specific collimators. Type Resolution Resolution Resolution parallel Size Septa Length Penetration Efficiency (mm) (mm) (mm) hole Pb Hole Shape (mm) (mm) (mm) @ 140 keV (cpm/μCi) @ 0 cm @ 10 cm @ 20 cm typical LEGP Hexagonal 1.40 0.18 24.7 2.1% 277 3.9 8.9 14.7 typical LEHR Hexagonal 1.22 0.15 27.0 1.7% 168 3.7 7.4 11.9 LEGP-CZT Square 1.80 0.20 31.0 1.6% 356 (+29%) 2.2 8.9 15.6 LEHR-CZT Square 1.85 0.15 40.9 1.8% 226 (+35%) 2.2 7.4 12.6

    [0037] FIG. 4 shows graphically the comparison of the four collimators shown in Table 3. The dashed curves represent the system resolution for typical LEGP (Low Energy General Purpose, short dash) and LEHR (Low Energy High Resolution, longer dash) hexagonal parallel-hole lead (Pb) collimators. By way of illustration, the solid lines represent the system resolution for two potential designs for pixel-registered square-parallel-hole lead (Pb) collimators. The thinner solid line with open circle markers represents the LEGP-CZT collimator response; the thicker solid line with open square markers represents the LEHR-CZT collimator response. The system resolution at 10 cm distance from the face of the collimator was matched for the two LEGP and two LEHR collimators, as an example of a possible design. Of course the resolution as a function of distance from the collimator can be traded off for system efficiency (sensitivity). For the particular choice illustrated in FIG. 4 and TAB. 3, the registered LEGP-CZT is 29% more efficient than the typical clinical hexagonal hole LEGP collimator. Likewise, the LEHR-CZT is 35% more efficient than the typical clinical hexagonal hole LEHR collimator. The optimization of the design depends on the application task and could be accomplished by Monte Carlo modeling simulation of the system response coupled with iterative reconstruction employing a full physics description of the detector-collimator response.

    [0038] As will be apparent to one skilled in the art, a pixel-registered collimator response has no significant dependence on the detector spatial resolution, hence the system resolution response function is approximately a straight line. In contrast, for a scintillator with Gaussian pixels, the system resolution is a quadrature summation of the detector and collimator resolutions, hence the resolution response function is approximately parabolic. The advantage for closer distances belongs to the CZT gamma camera with pixel-registered collimator. It will be apparent to one skilled in the art that the spatial resolution for a square-hole collimator is not isotropic in the plane. In our experience with such collimators as used in Molecular Breast Imaging (MBI), this off-axis resolution is not a significant issue.

    [0039] FIG. 5 depicts a process flow diagram that can be followed to improve or optimize the design and operation of a CZT gamma camera. We first describe a simple Figure of Merit (FoM) that illustrates one way to optimize the gamma camera by seeking configurations that maximize the FoM. By way of illustration, we choose a simple function that represents performance and cost. A common rule of thumb for gamma camera performance is the ratio of Efficiency (often called sensitivity) and the square of system resolution. We further want to consider cost and make a simplifying assumption that performance and cost are equally weighted. Thus, one possible FoM=Efficiency/(Cost*Resolution.sup.2). The step of calculating this FoM for any particular SPECT system is shown as step 16 in FIG. 5 and the goal (step 17) is to search for an optimal (maximum, in this case) FoM.

    [0040] As shown in FIG. 5, the first step 10 in this optimization process is to choose an application and determine performance requirements. This introduces constraints on the design optimization. The next step 11 is to choose a gamma camera design architecture. In current SPECT systems the gamma cameras are flat panels. As FIG. 3 reminds us, CZT detectors can be configured as curved panels. As we will discuss in FIG. 8, circular or elliptical detector surfaces can also be formed with CZT detectors. The next step 12 is to choose a collimation scheme, which may comprise, for example, parallel or focused holes (hexagonal or registered square), slit-slat or rotating slat, multi-pinhole, coded aperture, or Compton (electronic collimation). The next step 13 is to simulate digitally an application object, typically a portion of a human body or a digital phantom with points, lines, and volumes of radioactivity. Excellent computational human phantoms exist for this purpose. Those skilled in the art will understand that step 13 can be performed before step 11 and/or step 12 without changing the results.

    [0041] The next set of steps 14 through 17 are performed one or more times in an optimization loop. Step 14 considers the choices made in steps 11 and 12 for gamma camera geometry and collimation scheme. This imaging system and an acquisition sequence (for example, helical or stop-and-shoot SPECT) are modeled. Those skilled in the art understand that better results will generally be obtained when more of the physics of gamma emission imaging (such as collimator-detector response) are included in the model. For a particular distribution of radiopharmaceutical tracers in a computational phantom (step 13) the application of the forward projection model (step 14) will produce a simulation of the data that a physical system would acquire. This data set is usually called a sonogram for typical SPECT data acquisition. Step 15 applies iterative reconstruction algorithms, such as OSEM, to the sonogram data to produce a model of the object (step 13). The FoM is calculated in step 16. Typically, the efficiency and resolution can be assessed by using geometric objects as the computational phantom (step 13).

    [0042] The first passage through the optimization loop (steps 14 through 17) establishes a baseline FoM. The next step 17 is to search for an optimal FoM by changing the imaging system and/or the imaging sequence. According to various embodiments of this invention to be discussed below, this optimization may be achieved by sparsely placing CZT gamma detectors 20 in a gamma camera, such as 22 or 23, and optionally adding detector and/or collimator motion to compensate for the “missing” detectors in the sparse arrangement. Steps 14 through 17 may be applied once or iteratively. The optimization process may be driven manually or by computer algorithms that assist the decisions about which gamma detectors 20 to eliminate in a sparse design. Below we will present several embodiments of this invention and will estimate the FoM for each configuration. To be performed rigorously, all the optimization process steps 10 through 17 should be performed. However, to illustrate the invention, steps 13 and 15 have not been performed and steps 14 and 16 have been approximated with simple models.

    Flat Panel CZT Gamma Cameras

    [0043] FIG. 6 depicts pairs of gamma cameras labeled A through J. Each pair is shown side by side in a view from the patient. In clinical practice, the pair of gamma cameras would be mounted in opposition (or at about a right angle) on a SPECT gantry. The pair in FIG. 6A represents the normal flat panel scintillator (NaI & PMT) gamma cameras 24. To compare the FoM for various configurations, we will assert that both scintillator cameras in FIG. 6A are fitted with LEHR collimators (see TAB. 4). The efficiency, resolution, and cost of the following alternative designs B through J will be compared to this baseline reference configuration A. Thus the efficiency, resolution, and cost for design A, two NaI gamma cameras and two LEHR collimators, are each normalized to 1.0. Thus the relative Figure of Merit (FoM)=Efficiency/(Cost*Resolution.sup.2)=1.0 to help compare the various design ideas. Table 4 below lists the comparisons.

    [0044] The pair of flat panel gamma cameras in FIG. 6B represents flat panel CZT gamma cameras 22, as described in Table 2 and FIG. 2. We assert that the relative cost for a fully populated CZT gamma camera is about 2.75 times the cost of the baseline scintillator camera of the same field of view. We also assert that a CZT-LEGP collimator can be used with approximately 2.06 times greater efficiency than hexagonal LEHR, but the same average resolution. The FoM of this approach is 2.06/(2.75*1.00.sup.2)=0.75.

    [0045] FIGS. 6C and 6D depict two potential hybrid approaches to reduce the system cost of deploying CZT gamma cameras. The basic idea is to mix one scintillation camera with one CZT camera. In the pair of flat panel cameras of FIG. 6C “1 NaI+1 CZT” the scintillation camera is coupled to an LEGP collimator and the CZT camera has an LEHR-CZT collimator, resulting in 1.44 times the base efficiency, 1.03 times the resolution, and 1.88 times the cost, for a FoM of 0.73. Another alternative is in FIG. 6D wherein the flat panel scintillator camera has LEGP and the CZT camera has LEGP-CZT collimation. This results in relative efficiency of 1.82, relative resolution of 1.09, relative cost of 1.88, and FoM=0.81. These two hybrid options are for illustration and do not represent an optimized solution, which would require simulation and then confirmation using phantoms.

    [0046] FIG. 6E introduces the concept of sparse acquisition by decimating the number of CZT aggregator modules deployed in the CZT gamma cameras. For purpose of illustration the decimation factor is two: half of the available locations for CZT modules are not filled. Those skilled in the art will recognize that other decimation factors and patterns may be usefully deployed. The pair of flat panel CZT gamma cameras in FIG. 6E “two Moving Slat CZT” gamma cameras shows 6 vertical slats 25 (4.4 cm wide×39.6 cm long) where nine AMs are located, alternating with 6 vertical spaces of white (also 4.4 cm wide) where there are no gamma detectors. Note that in this example embodiment, the upper and lower gamma cameras have a complementary pattern of occupied and vacant detector slats. Thus, in a standard SPECT acquisition where the two cameras are opposed by 180 degrees, all regions of the body are being imaged, so dynamic processes can be reliably sampled.

    [0047] As illustrated in FIG. 7, half of the acquisition time for a planar projection would be spent in the position shown in FIG. 7a, then all the slats would shift quickly right (arrows 26 depicted in the upper camera of FIG. 7b) or left (arrows 26 depicted in the lower camera of FIG. 7b) by 4.4 cm to the blank positions and the second half of the planar acquisition would be acquired. One skilled in the art will recognize various means for achieving this detector motion. By way of illustration, a mechanical frame 27 rigidly supports the sparsely populated gamma camera detector slats (columns) 25. The frame can be moved one column width (4.4 cm in this example) by one or more linear actuators 26. Of course, other components of a functioning gamma camera are not shown, such as electronic readout circuits, power, control circuits, radiation shielding, and collimator. The collimator can be fixed in position and cover the entire gamma camera area, or a slat of collimator can be attached to each slat of detectors and they can shift position together by means of linear actuators in a manner similar to that discussed for the detectors. Thus, 50% of the time is spent in each of two positions and a contiguous planar image is acquired. For this configuration embodiment in FIG. 6E, the relative efficiency and cost are both half of configuration 6B, thus the FoM=0.75 is the same.

    [0048] In FIG. 6F with “1 NaI+1 Moving Slat CZT” gamma cameras, we combine the shifting motion slat concept of FIG. 6E with the hybrid concept shown in FIG. 6D, resulting in relative efficiency of 1.30 and cost of 1.19 with FoM=0.92. It will be apparent to one skilled in the art that the slats could be oriented in the horizontal direction with vertical shifting motion. It is also apparent that different decimation factors could be used to reduce cost. But efficiency and cost both depend linearly on the number of CZT detectors deployed.

    [0049] FIG. 6G shows a variation on the concept of decimated slats (vertical or horizontal). In this variation there is no motion required to obtain projections of the gamma photon emission distribution of radiopharmaceuticals in the human body under examination. Instead, diverging fan (focal point behind detector) pixel-registered square-hole collimators are attached to each slat of gamma detectors. The field of view of adjacent fans will overlap in distal regions of the imaged body, but iterative maximum likelihood estimate reconstruction should have no difficulty in resolving redundant overlapping projections.

    [0050] One perceived disadvantage may be that no planar projection is acquired as with parallel-hole collimators, but synthetic planar projections can be easily derived from the reconstructed 3D SPECT image. Without optimizing the possible diverging fan collimator designs, we estimate that for the embodiment in configuration 6G “two Slat CZT” gamma cameras and “two Fan” collimators, the relative efficiency is 1.5, but the relative resolution suffers from the diverging collimation, and the FoM is only about 0.36. Similarly, in a hybrid configuration embodiment 6H “1 NaI+1 Slat CZT” the FoM is only 0.56. As will be apparent to one skilled in the art, the concept of overlapping fan collimators is an inefficient compensation for decimation of CZT detectors intended to reduce cost.

    [0051] Other embodiments for decimating the number of CZT AMs in a moving detector and collimator concept are contemplated, although they are not depicted here. For example, FIG. 6I illustrates a large checkerboard pattern for the factor of two decimation of CZT gamma detector AMs. Note that complementary detector patterns are shown for the upper and lower cameras, as also preferred for FIGS. 6E and 6G. As shown for illustration, each gray square consists of a 3×3 array of AMs forming a 13.2 cm×13.2 cm subarray 27. Other patterns, including different size checkerboard elements, will be apparent to one skilled in the art. As illustrated in FIGS. 6E and 6G, either motion of the slats 25 or use of diverging fan pixel-registered square-hole collimators could be deployed in this configuration. However, by way of illustration we prefer to deploy multiple-pinhole collimators positioned to illuminate each of the occupied detector subarrays 27. The use of multiple-pinhole collimation with overlapping projections has been proven in commercial products for small animal (“preclinical”) SPECT (e.g., Trifoil Imaging of Chatsworth, Calif. and MiLabs B. V. of The Netherlands) and in human cardiac SPECT (i.e., GE Healthcare, a segment of General Electric Co.). Again, synthetic planar projections can be easily derived from the fully reconstructed SPECT image. Without optimizing the possible multi-pinhole collimator designs, we estimate that for embodiment 6I “two Checker CZT” a relative efficiency of about 2.0 (which could be higher), resolution of 0.87, and cost of 1.38 resulting in an FoM of 1.91. The hybrid variation for this configuration is in FIG. 6J “1 NaI+1 Checker CZT” with FoM of 1.44. These are the first two embodiments illustrated here that have FoM>1, indicating an improvement over the standard scintillator gamma cameras with LEHR collimators, as shown in FIG. 6A.

    [0052] This preference for spatially separated smaller gamma cameras, such as the small 3×3 AM cameras 27 of FIG. 6I, can be generalized. This is one of the advantages of designing SPECT systems with modular CZT detectors. The use of multiple pinhole collimators can also be recommended. We have not discussed technical development risks and time-to-market considerations which may influence the choice of configuration, at least for the early introduction of CZT SPECT gamma cameras.

    [0053] The use of curved detector panels, as indicated in FIG. 3b, would be available as a variation to all the previously discussed examples of preferred flat panel embodiments. We anticipate an improvement factor for relative efficiency of about 1.5 times. Thus, the FoM for a curved panel would be about 1.5 times that for configurations shown in FIG. 6B-J.

    Ring CZT Gamma Cameras

    [0054] FIG. 8 introduces another preferred embodiment of the present invention with a configuration variation that extends the concept of curved panels to the point where the two panels join to form a circular or elliptical ring of AMs 20 encircling the human body 30. In FIG. 8a a ring of detectors 32 about 1.2 m in diameter can be formed from 86 contiguous AMs 20 (86×4.4 cm=378.4 cm=π×120.4 cm=diameter). Recall that the cost of two fully populated flat-panel CZT gamma cameras was estimated to be 2.75 times the reference scintillation camera and that two×nine×12=216 AMs 20 were required. Thus the relative cost of a fully populated single ring of 86 AMs is 2.75×(86/216)=1.09. Of course, this invention contemplates that a tiled ring of CZT or other semiconductor gamma detector modules can be sparsely filled and the effects of missing detectors can be compensated by optimized collimation and potential detector motion.

    [0055] Use of converging fan pixel-registered square-hole collimation would be technically feasible, but the spatial resolution would be poor since the distance from collimator to the interior human body is too large. A much better option is multiple pinhole collimation, for example located on a ring 31 of about 0.6 m diameter with about one pinhole per AM, for example. The density of pinholes could be varied and shutters could be employed to open or close various pinholes to manipulate the SPECT acquisition parameters. This arrangement would provide a magnification of about 1.0 at the central axis and increasing toward the surface of the human body. Based on experience with small animal multi-pinhole SPECT and human cardiac multi-pinhole SPECT, we estimate very conservatively that the relative efficiency of a single circular ring would be about 2.5 and thus an FoM of 3.0. Additional rings can be added, as with PET (Positron Emission Tomography) systems, and they are labeled configuration K, L, and M for one, two and 3 rings. Of course, cost and efficiency both scale linearly with the number of CZT AMs, so the FoM remains about 3.0 for all configurations, assuming fully populated rings. Again synthetic planar projections can be produced from the fully reconstructed 3D SPECT image.

    [0056] Configurations N, O, and P represent one, two, and three elliptical rings as a variation on configurations K, L, and M with circular rings. Ellipses more closely follow the contours of most human bodies 30 (bariatric patients may have a more circular transaxial cross-section) and fewer CZT detectors are required for some cost savings. In the illustrative example shown in FIG. 8b, 79 AMs are deployed in each elliptical detector ring with principle axes of about 1.2 m and 1.0 m. The elliptical ring of multiple pinholes 31 would have proportional principle axes of 0.6 m and 0.5 m. We assume that the relative efficiency of an elliptical ring 32 would be similar to a circular ring 32, but the relative cost would be about 9% less, so the FoM would increase to about 3.27.

    [0057] TAB. 4 lists the relative cost and efficiency (sensitivity) of the flat panel gamma camera configurations B-J depicted in FIG. 6 and the ring configurations K-P depicted in FIG. 8. It is apparent that much higher efficiencies can be achieved with ring CZT gamma cameras. This can be translated into a clinical tradeoff between better image quality and shorter examination times or lower patient doses. Apparently the better technical choice is a fully populated circular or elliptical ring of CZT gamma cameras, scalable to multiple rings (like PET). However, the best business choice may be a sparse flat or curved panel for a first generation CZT SPECT system, considering time to market and investment required to develop the first circular or elliptical CZT SPECT system.

    Ring CZT Compton Gamma Cameras

    [0058] FIG. 8 also illustrates another preferred embodiment of the present invention with another variation on the collimation. The technical risk is higher, so this configuration will likely not be the first introduced commercially, but the potential performance is game changing. The collimation in this preferred embodiment is provided by a Compton scatter (“electronic collimation”) ring 31. Typically, lower atomic number elements have higher Compton scatter cross-sections, so 2D pixelated silicon (Si) gamma detectors are a good choice. To reconstruct a Compton scatter image requires information about the 3D position, energy, and timing of the atomic electron that undergoes a Compton scatter event in the Si detector as well as the 3D position, energy, and timing of the photoelectric interaction in the CZT. If the two events are close in time (“coincidence”) and the energies add to the expected full peak energy of the emitted gamma photon, then the cone of response can be determined for the origin of the gamma photon. Typically, a stack of multiple Si 2D detectors (often pixelated using crossed-strip readout) can be used to enhance the Compton scatter probability. Without detailed modeling, we estimate, based upon general knowledge of one skilled in the art, that the efficiency of a Compton camera ring can be five to 20 times higher than the reference scintillator SPECT system. Also, because a ring of Si detector stacks must be added to the CZT ring, the relative cost is about two. Thus the FoM is about 2.5-10 times higher. We use the more conservative estimates in TAB. 4. Configurations Q, R, and S represent one, two, and 3 circular rings, while T, U, and V represent one, two, and three elliptical rings with FoM of 7.5 or 8.2 for circular or elliptical rings. It is apparent that much higher efficiencies can be achieved with Compton ring CZT gamma cameras. This can be translated into a clinical tradeoff between better image quality and shorter examination times or lower patient doses.

    [0059] As a summary of the various configurations of CZT gamma cameras disclosed as embodiments of the present invention, the following TAB. 4 lists all the configurations and their approximate relative efficiencies, resolutions, costs, and Figures of Merit.

    TABLE-US-00004 TABLE 4 Compilation of relative merits of the various embodiments for flat-panel, ring, and Compton ring SPECT systems using CZT gamma cameras. Note that efficiency and FoM for curved panel detectors of embodiments B through J may be a factor of 1.5 times higher. Effi- Reso- Label Gamma Cameras Collimators ciency lution Cost FoM A 2 NaI 2 LEHR 1.00 1.00 1.00 1.00 B 2 CZT 2 2.06 1.00 2.75 0.75 LEGP-CZT C 1 NaI + LEGP + 1.44 1.03 1.88 0.73 1 CZT LEHR-CZT D 1 NaI + LEGP + 1.82 1.09 1.88 0.81 1 CZT LEGP-CZT E 2 Moving Slat CZT 2 1.03 1.00 1.38 0.75 LEGP-CZT F 1 NaI + LEGP + 1.30 1.09 1.19 0.92 1 Moving Slat CZT LEGP-CZT G 2 Slat CZT 2 Fan 1.50 1.75 1.38 0.36 H 1 NaI + LEHR + 1.25 1.37 1.19 0.56 1 Slat CZT Fan I 2 Checker CZT MPH 2.00 0.87 1.38 1.91 J 1 NaI + LEHR + 1.50 0.94 1.19 1.44 1 Checker CZT MPH K 1 Circular Ring MPH 2.50 0.87 1.09 3.00 L 2 Circular Rings MPH 5.01 0.87 2.19 3.00 M 3 Circular Rings MPH 7.51 0.87 3.28 3.00 N 1 Elliptical Ring MPH 2.50 0.87 1.01 3.27 O 2 Elliptical Rings MPH 5.01 0.87 2.01 3.27 P 3 Elliptical Rings MPH 7.51 0.87 3.02 3.27 Q 1 Circular Ring Compton 12.52 0.87 2.19 7.51 R 2 Circular Rings Compton 25.04 0.87 4.38 7.51 S 3 Circular Rings Compton 37.56 0.87 6.57 7.51 T 1 Elliptical Ring Compton 12.52 0.87 2.01 8.17 U 2 Elliptical Rings Compton 25.04 0.87 4.02 8.17 V 3 Elliptical Rings Compton 37.56 0.87 6.03 8.17

    Higher Energy Isotope SPECT

    [0060] The use of 5 mm thick CZT gives similar stopping power to 9.5 mm (⅜″) NaI of about 85% for Tc-99m (140 keV), the most commonly used medical isotope for SPECT. This thickness is adequate for the most common medical isotopes, such as Tc-99m (140 keV), Tl-201 (70 keV), Xe-133 (81 keV), Ga-67 (90 keV), and I-123 (159 keV). Increasing the CZT thickness to 1.0 cm will increase the stopping power and thus the detection efficiency for higher energy medical isotopes, such as In-111 (171 & 245 keV) and I-131 (364 keV). Increasing the thickness of CZT to 1.0 cm results in higher efficiency for all isotopes and almost doubles the stopping power for 1-131 (364 keV). As will be apparent to one skilled in the art, the stopping power for non-normally incident gamma photons is increased by an interaction path longer by 1/cosine (incident angle). This enhanced efficiency favors somewhat the collimation techniques such as pinhole and fan.

    [0061] With reference once again to FIG. 1 we disclose a method for reducing the cost of a CZT SPECT gamma camera designed to be efficient for both lower and higher energy SPECT isotopes. A D-Matrix™ aggregator module 20 with four gamma modules 21 could be built with two 0.5 cm thick CZT gamma modules 21 and two 1.0 cm thick gamma modules 21, preferably in a checkerboard pattern. In a preferred embodiment, the front (cathode) surface of the AM 20 would be flat, that is, the cathodes would be coplanar but the anodes would be staggered by the 0.5 cm difference in CZT sensor thickness. Iterative reconstruction of the SPECT projections would utilize a full physics description of the collimator-detector response, including the checkerboard CZT thickness pattern. As will be apparent to one skilled in the art, other patterns of mixing detector thicknesses can be devised to accomplish the same goal of providing for more efficient detection of higher energy isotopes, without increasing the price too dramatically.

    [0062] While preferred embodiments of the present invention are shown and described, it is envisioned that those skilled in the art may devise various modifications of the present invention without departing from the spirit and scope of the appended claims.