Sparse Acquisition Gamma Cameras
20170285191 · 2017-10-05
Inventors
Cpc classification
G01T1/243
PHYSICS
A61B6/4275
HUMAN NECESSITIES
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]
[0021]
[0022]
[0023]
[0024]
[0025]
[0026]
[0027]
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]
[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]
[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
[0034]
[0035] A gamma camera could be designed to mechanically transform between the flat panel configuration of
[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]
[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]
[0040] As shown in
[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]
[0044] The pair of flat panel gamma cameras in
[0045]
[0046]
[0047] As illustrated in
[0048] In
[0049]
[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,
[0052] This preference for spatially separated smaller gamma cameras, such as the small 3×3 AM cameras 27 of
[0053] The use of curved detector panels, as indicated in
Ring CZT Gamma Cameras
[0054]
[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
[0057] TAB. 4 lists the relative cost and efficiency (sensitivity) of the flat panel gamma camera configurations B-J depicted in
Ring CZT Compton Gamma Cameras
[0058]
[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
[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.