The present invention relates generally to scintillation detectors for imaging devices, and more particularly, to scintillation crystal materials used in scintillation detectors
Nuclear medicine is a unique medical specialty wherein radiation is used to acquire images which show the function and anatomy of organs, bones or tissues of the body. Radiopharmaceuticals are introduced into the body, either by injection or ingestion, and are attracted to specific organs, bones or tissues of interest. Such radiopharmaceuticals produce gamma photon emissions which emanate from the body and are captured by a scintillation crystal, with which the photons interact to produce flashes of light or “events.” Events are detected by an array of photodetectors, such as photomultiplier tubes, and their spatial locations or positions are calculated and stored. In this way, an image of the organ or tissue under study is created from detection of the distribution of the radioisotopes in the body.
A number of scintillation crystals are known for use in nuclear medicine and each has particular advantages and disadvantages depending upon the type of nuclear imaging to be performed, the required image quality, cost of manufacture, light sharing scheme, type/size of a photomultiplier tube (PMT), detector/gantry geometry, performance in terms of sensitivity and spatial resolution, etc. Similarly, there are a number of known crystal surface treatments, which can affect image quality. Prior art approaches generally employ a series of channels cut or otherwise formed in the light output surface of the scintillation crystal. The purpose of the channels is to guide scintillation photons through the crystal and to minimize the lateral movement (“spreading”) of photons within the crystal. By minimizing spreading, the sensors are more likely to receive photons emanating from scintillation events located directly below them. The channels therefore function somewhat as collimators to direct or channel the photons in a direction substantially perpendicular to the plane of the scintillation crystal and/or perpendicular to the emission face. In one such prior art approach, as described, e.g., in U.S. Pat. No. 6,881,960, a scintillation crystal is formed with two sets of substantially orthogonal channels on its light output side which form a plurality of rectangular solid portions. The channels of this scintillation crystal are identical and have a uniform depth. Another prior art approach is disclosed in U.S. Pat. No. 6,841,783, which is incorporated by reference herein in its entirety. In the '783 patent, a scintillation crystal (Na(Tl)) is formed with two sets of substantially orthogonal channels on its light output side, however, the channels have varying depths.
While the above-identified, and other, surface treatments for output surfaces of scintillation crystals are known, there are few surface treatment methods for the input surfaces of scintillation crystals to prevent spreading, particularly with respect to relatively thick NaI(Tl) scintillation crystals that can be used for PET, or high energy SPECT applications. Indeed, a significant challenge in the development NaI based PET, and high energy SPECT, systems is that relatively thick scintillation crystals must be utilized in order to compensate for the low stopping power of NaI(Tl) crystals receiving 300 keV or 511 keV gamma photons. In sum, employing relatively thick crystals increases the scintillation light spread within the detector and degrades the spatial resolution
What is needed then are surface treatments for the input sides of scintillation crystals, particularly relatively thick NaI(Tl) scintillation crystals, that address the above-identified deficiencies.
The present invention addresses the above identified deficiencies by providing a method and surface treatment for scintillation crystals used in radiographic imaging devices, particularly imaging devices utilizing high energy particles, e.g., 300, keV, 511 keV, etc. gamma photons, such as PET, and high energy SPECT, imaging devices. In one embodiment, the invention comprises a method for enhancing the light output efficiency of a scintillation crystal having an input face and an output face. The input face is configured for receiving a radiation particle and the output face is configured for emitting light photons produced in response to absorption of a radiation particle in the crystal. The method comprises providing the input face with a first set of channels having a first channel depth formed in a portion of the input face, wherein each channel extends in a first direction along the portion of the input face in a substantially parallel, spaced apart relationship with other channels in said first set. The method also comprises providing a second set of channels having a second channel depth formed in a portion of the input face wherein each channel of the second set of channels extends in a second direction along a portion of the input face in a substantially parallel, spaced apart relationship with other channels in the second set. The first and second directions are non-parallel to each other.
In one embodiment, the invention comprises a method for enhancing the light output efficiency of a scintillation crystal having an input face and an output face. In some embodiments of the invention, the first set of channels and the second set of channels intersect with one another. In some embodiments, the first set of channels and the second set of channels are orthogonal. In some embodiments, the channels comprise a pair of sides and a base side, and the pair of sides are parallel to one another. In some embodiments, the base side is substantially orthogonal with respect to the pair of sides. In some embodiments, the first set of channels has a depth that is equal to the depth of the second set of channels within a smallest obtainable manufacturing tolerance. In some embodiments, the depth of the channels is less than 10 millimeters. In some embodiments, the scintillation crystal has a thickness less than 5 centimeters and preferably, between 2 and 3 centimeters. In some embodiments, the scintillation crystal comprises NaI(Tl).
In some embodiments, the invention comprises a method for enhancing the light output efficiency of a NaI(Tl) scintillation crystal configured for use in PET or SPECT. In such embodiments, the method comprises providing an input face of the crystal with a first set of channels having a first channel depth formed in a portion of the input face, wherein each channel extends in a first direction along said portion of the input face in a substantially parallel, spaced apart relationship with other channels in said first set, and providing a second set of channels having a second channel depth formed in a portion of the input face, wherein each channel of the second set of channels extends in a second direction along a portion of the input face in a substantially parallel, spaced apart relationship with other channels in the second set, and wherein the second direction is non-parallel with the first direction.
In some embodiments wherein the crystal is a NaI(Tl) scintillation crystal, the crystal has thickness that is between 2 and 3 centimeters. In some embodiments wherein the crystal is a NaI(Tl) scintillation crystal, the first set of channels intersect and are substantially orthogonal with respect to the second set of channels. In some embodiments wherein the crystal is a NaI(Tl) scintillation crystal, a channel comprises a pair of sides and a base side and the pair of sides are parallel to one another. In some embodiments wherein the crystal is a NaI(Tl) scintillation crystal, the base side is substantially orthogonally disposed with respect to the pair of sides. In some embodiments wherein the crystal is a NaI(Tl) scintillation crystal, the first set of channels has a depth that is equal to the second set of channels within a smallest obtainable manufacturing tolerance. In some embodiments wherein the crystal is a NaI(Tl) scintillation crystal, the depth is less than 10 millimeters, and preferably between 2-8 centimeters.
In some embodiments according to the invention, a NaI(Tl) scintillation crystal having enhanced light output efficiency comprises an input face and an output face, wherein the input face is configured for receiving a radiation particle and the output face is configured for emitting a photon in response to absorption of the radiation particle. A first set of channels is formed in the input face and extend in a first direction along the input face in a substantially parallel, spaced apart relationship. A second set of channels is formed in the input face and extend in a second direction along the input face in a substantially parallel, spaced apart relationship. The first and second directions are non-parallel with one another.
The invention will now be more fully described by way of example with reference to the accompanying drawings in which:
a-c are illustrations of a scintillation crystal and detector arrangement according to the invention;
a-2c are illustrations of a plurality scintillation crystal input face surface treatments;
a-d are graphical representations of DSI simulator validation results (dashed line=detector, (a) paint surface case: representing Lambertian model, (b) metal coated surface case: representing specular reflector, (c) rough surface (alpha=45), (d) 4 holes f=5 cm, approximately (±10, ±10 cm) from a point of light source;
a-d are graphical representations of validation with experimental cases (dashed line (a) Ecam (⅜ inch) Detect, (b) Ecam (⅜ inch) DSI, (c) Ecam+ (⅝ inch) Detect, (d) Ecam+ (⅝ inch) DSI;
a-d are graphical representations of the expected light response function for a plurality of crystal surface treatments ((a) 1″ standard (b) holes (c) entrance grooves (d) retro-reflector; dashed line represents ⅜ inch LRF);
a-d are graphical representations of examples of CR bound curve as a function of position ((a) entrance side grooves (b) holes (c) retro-reflector (c) duet (e) 1″ standard (f) ⅝ inch E.cam plus; dashed line represents CR bound for ⅜ inch E.cam as a reference CR);
a-b are graphical representations of the expected performance of crystal surface treatments ((a) CR lower bound (b) light collection efficiency; all values are normalized against 1″ standard case performance);
a-d are graphical representations of Depth of Interaction dependency of Light Response Functions (LRFs) ((a) 1″ standard (b) retro-reflector (c) duet (c) entrance grooves; BLK: shallow events, BLU: middle events, R: deep events);
a-b are graphical representations of LRF control effects on spatial resolution and Light Output efficiency (LOE); and
a-f are graphical representations of Mean Point Spread Function (MPSF) for a plurality of crystal surface treatments (MPSFfwhm (a) 1″ standard (b) retro-reflector (c) ⅜ inch E.cam (d) Duet (e) holes (f) entrance grooves).
The present invention will now be described and disclosed in further detail. It is to be understood, however, that the disclosed embodiments are merely exemplary of the invention and that the invention may be embodied in various and alternative forms. Therefore, specific structural and functional details disclosed herein are not to be interpreted as limiting the scope of the claims, but are merely provided as an example to teach one having ordinary skill in the art to make and use the invention.
Adverting now to the figures,
Referring now to
Holed surface 28 can be formed by drilling one or more holes 34 into the surface of the input face by suitable means, such as by means of a drill, and in the embodiment illustrated each hole 34 has a diameter between 1-2 mm and is separated from another by a distance of 1-4 mm. Holes 34 can be arranged in rows, offset rows, or can be randomly disposed about the input face of the crystal 18. In some embodiments, the height, or depth of holes 18 can be at least 25 mm.
Retro-reflector surface 30 can be formed from a plurality of pyramid-shaped structures 36 having a height between 1-2 mm and a width of base between 2-4 mm. The pyramid-shaped structures 36 of the retro-reflector surface 30 can be formed by suitable means such as milling.
Scintillation crystal 23 comprising grooved surface 32 is shown as comprising a first set 38 of channels 42 and a second set 40 of channels 42. In the embodiment illustrated, channels 42 of the first set 38 are formed in a spaced apart parallel relationship with respect to one another and extend in a first direction. Similarly, channels 42 of the second set 40 are formed in a spaced apart parallel relationship with respect to one another and extend in a second direction. In the embodiment illustrated, first set 38 is formed substantially perpendicular to the second set 40. The two sets of channels 42, thus, are illustrated as forming a plurality of substantially rectangular solid portions 43, each which comprises a portion of the input face 24. In the embodiment illustrated, channels 42 have a generally square, or rectangular, cross sectional shape and have a width of approximately 0.5 millimeters (mm), a depth of between 2-8 millimeters (mm), and form a rectangular solid portion 43 that has a width of approximately 5 millimeters (mm).
In the embodiment illustrated, channels 42 of each set are preferably substantially uniform. That is, each is preferably similar in width, depth, and cross-sectional shape. In a preferred embodiment, the channels include a pair of parallel sides 44 and a base side 46, which is orthogonally disposed with respect to the parallel sides 44. It should be appreciated that while first and second sets 38, 40 of channels 42 are illustrated as being orthogonally disposed with respect to one another, other angular configurations may be desirable. Similarly, while sides 44 are described as being parallel with respect to one another and orthogonal with respect to base side 46, one or more of the sides 44 may be disposed at other angles. Similarly, base side 46, while being described as being essentially planar by virtue of its being orthogonally disposed with respect to sides 44, may be angled or arced. Preferably, however, channel 42 has a generally square or rectangular cross sectional shape. Channels 42 may be formed by milling crystal 18, or by other suitable process. Channels 42, while illustrated as comprising air gaps, may be configured to comprise other materials, such as reflective powder, or other solid material, etc. While in a preferred embodiment, scintillation crystal 18 is less than 5 centimeters thick from input face to output face, and preferably between 2-3 centimeters, for other crystal thicknesses, the ratios of the depths of the channels to the thickness of the crystal may vary as a function of crystal thickness as the relationship of channel depth to crystal thickness is non-linear, but may be determined readily through experimentation and testing. It should be understood that the numbers described above are merely exemplary, and the spacing, width, and depth may be changed and still remain within the scope of the invention.
Experiments
Conditions
Referring now to
The experiments were conducted on the detector geometries, optical conditions and PMT allocations described in
The surface modeling, sampling and user interface of the DSI simulator was modeled after the Detect Simulator and was validated with the Detect simulator. The two simulators were then validated with known experimental cases of E.cam and E.Cam+ systems, commercially available from Siemens Medical Solutions USA, Inc. of Hoffman Estates, Ill.
Quality Factors
The statistical fluctuation is generally dependant upon the Light Output Efficiency (LOE) of a crystal, the quantum efficiency of PMT photocathodes and the noise characteristics of a photomultiplier tube. The conversion efficiency and noise characteristics are intrinsic properties of a PMT. Consequently, LOE and, therefore, high signal to noise ratio, is important to the design of an optimal detector.
Anger type positioning algorithms are the predominant positioning estimators for scintillation cameras because of their relatively simple hardware realization and reasonable performance. However, such methods, e.g., linear least square type, are known to be biased and not efficient estimators. Alternative positioning algorithms utilizing statistical properties include maximum likelihood method, transfer function approach and statistics based positioning algorithms. Although such type estimators typically require additional electronics and computational resources, it can be worthwhile to utilize optimum positioning algorithms to improve the system performance, including edge effect artifacts.
Generally, the shape of the Light Response Function (LRF) determines the fundamental performance of a given system and it is most sensitive to thickness, the shape and surface treatment, and light guide configuration of a crystal. The quality factor information was extracted from the LRF, the quality factors are:
1.) Light output efficiency: LOE(%)=(detected photons/total simulated photons)×100
2.) Cramer Rao Lower bound (unbiased case): The CR lower bound, derived in equation (1), was used to determine the best standard deviation, i.e., spatial resolution obtainable from the given LRF of detector geometry. Note that the CR lower bound is a function of the LRF, in other words, the intrinsic shape (slope) of the LRF of a given system determines the spatial resolution of the system. Also, CR provides a scalar performance index that is useful for selecting the optimum configuration of an imaging system.
3) Mean of Point Spread Function FWHM (MPSFfwhm): Although CR lower bound is useful when to examine the performance potential of different imaging system geometries, the equation (1) is derived based on the unbiased assumption so that it is inapplicable where bias is unavoidable situation. On the other hand, the FWHM of PSF of a camera is a hallmark of performance assessment. The individual FWHM of PSF at x and y direction is separately calculated at each test point defined in
Photon transport and interaction models at the boundary of optical components and random sampling models of DSI stimulator were examined against that of the Detect simulator. The validation of surface modeling including Lambertian, specular and rough case are shown in
Further validation of the two simulators was performed with known experimental cases, i.e., E Cam (⅜) and E Cam+(⅝).
Light Output Efficiency
Light Output Efficiency (LOE %) was tested for all proposed manipulation methods and the E Cam series. LOE generally represents the expected energy resolution performance of a given method. Table 1 summarizes LOE results. It should be appreciated that manipulating the entrance side of the crystal surface enhanced its light output efficiency.
Two Dimensional Light Distribution
Light Response Function (LRF)
In practice, Light Response Function (LRF) is obtained by moving a highly collimated point source on a known grid with desired precision. In the simulation, one dimensional light response function was derived from the two dimensional light distribution shown in
Cramer Rao Lower Bound
Cramer Rao (CR) Bound was derived assuming a total of 5 PMTs in a row were utilized in estimating event positions for one direction.
Performance
Expected performance of the crystals is summarized in
Depth of Interaction Decodeability
Depth of Interaction (DOI) dependent or independent LRF is an important property of a scintillation system, especially for positron imaging and single photon imaging with pinhole collimation. For example, DOI independent LRF is desired in SPECT systems for a uniform performance regardless of depth of event interaction while in PET systems DOI dependent LRF is desirable to allow easy DOI decoding features to minimize parallax errors,
Although some manipulation techniques provide DOI dependent LRF, conventional Anger style positioning algorithms are not suitable to decode depth depend events since they are designed to find a centroid of a given LRF. To have a DOI decoding feature based on LRF, sophisticated positioning algorithms based on statistics and/or detector characteristics may be employed.
LRF Control
The ease of LRF control property is an important quality factor.
Mean Point Spread Function
As previously discussed, CR lower bound is a useful tool for assessing system performance. However, point spread function or modular transfer function is a good objective intrinsic spatial resolution assessment tool for a given system. Although MPSF was conducted based on Monte Carlo simulation, it provides expected spatial resolution performance of a given system with a high confidence level. One thousand events were simulated at each test point described in
In view of the foregoing, it is seen that a NaI scintillation crystal having an input face comprising a grooved or channeled surface enhances a number of factors that affects crystal output, particularly light output efficiency, spatial and energy resolution, DOI dependency and uniformity. Consequently, a method of enhancing light output efficiency of relatively thick NaI scintillation crystals generally includes providing a plurality of grooved or channeled surfaces on the input face of the crystal.
It should be appreciated by those having ordinary skill in the art that while the present invention has been illustrated and described in what is deemed to be the preferred embodiments, various changes and modifications may be made to the invention without departing from the spirit and scope of the invention. Therefore, it should be understood that the present invention is not limited to the particular embodiments disclosed herein.