Patent Application
20090072155
The present invention relates generally to nuclear medicine and systems for obtaining nuclear medicine images. In particular, the present invention relates to reconstructing nuclear medicine images from time-of-flight (TOF) positron emission tomography (PET) data that has been time-truncated.
Nuclear medicine is a unique specialty wherein radiation is used to acquire images which show the function and anatomy of organs, bones or tissues of the body. The technique of acquiring nuclear medicine images entails first introducing radiopharmaceuticals into the body--either by injection or ingestion. These radiopharmaceuticals are attracted to specific organs, bones or tissues of interest (These exemplary organs, bones, or tissues are also more generally referred to herein using the term “objects”.). Upon arriving at their specified area of interest, the radiopharmaceuticals produce gamma photon emissions which emanate from the body and are then captured by a scintillation crystal. The interaction of the gamma photons with the scintillation crystal produces flashes of light which are referred to as “events.” Events are detected by an array of photo detectors (such as photomultiplier tubes) and their spatial locations or positions are then 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.
One particular nuclear medicine imaging technique is known as positron emission tomography, or PET. PET is used to produce images for diagnosing the biochemistry or physiology of a specific organ, tumor or other metabolically active site. The measurement of tissue concentration using a positron emitting radionuclide is based on coincidence detection of the two gamma photons arising from a positron annihilation. When a positron is annihilated by an electron, two 511 keV gamma photons are simultaneously produced and travel in approximately opposite directions. Gamma photons produced by an annihilation event can be detected by a pair of oppositely disposed radiation detectors capable of producing a signal in response to the interaction of the gamma photons with a scintillation crystal. Annihilation events are typically identified by a time coincidence between the detection of the two 511 keV gamma photons in the two oppositely disposed detectors; i.e., the gamma photon emissions are detected virtually simultaneously by each detector. When two oppositely disposed gamma photons each strike an oppositely disposed detector to produce a time coincidence event, they also identify a line(s)-of-response (LOR) along which the annihilation event has occurred. An example of a PET method and apparatus is described in U.S. Pat. No. 6,858,847, which patent is incorporated herein by reference in its entirety.
After being sorted into parallel projections, the LOR defined by the coincidence events are used to reconstruct a three-dimensional distribution of the positron-emitting radionuclide within the patient. In two-dimensional PET, each 2D transverse section or “slice” of the radionuclide distribution is reconstructed independently of adjacent sections. In fully three-dimensional PET, the data are sorted into sets of LOR, where each set is parallel to a particular detector angle, and therefore represents a two dimensional parallel projection p(s, φ) of the three dimensional radionuclide distribution within the patient—where “s” corresponds to the distance of the LOR from the center of the detector and “φ” corresponds to the angle of the detector plane with respect to the x axis in (x, y) coordinate space (in other words, φ corresponds to a particular LOR direction).
Coincidence events are integrated or collected for each LOR and stored in a sinogram. In this format, a single fixed point in f(x, y) traces a sinusoid in the sinogram. In each sinogram, there is one row containing the LOR for a particular azimuthal angle φ; each such row corresponds to a one-dimensional parallel projection of the tracer distribution at a different coordinate along the scanner axis. This is shown conceptually in
An event is registered if both crystals detect an annihilation photon within a coincidence time window τ (e.g, on the order of 4-5 nsec), depending on the timing properties of the scintillator and the field of view (FOV). The FOV is defined as the volume between the detectors; and a pair of detectors is sensitive only to coincidence events occurring in the FOV. Therefore, the need for physical collimation is eliminated and sensitivity is significantly increased. Accurate corrections (for example, attenuation correction) can be made for the self-absorption of photons within the patient so that accurate measurements of tracer concentration can be made.
The number of time coincidences detected per second within a FOV of a detector is the count rate of the detector. The count rate at each of two oppositely disposed detectors, A and B, can be referred to as singles counts or SA and SB, respectively. The time required for a gamma photon to travel from its point of origin to a point of detection is referred to as the time-of-flight (TOF) of the gamma photon. TOF is dependent upon the speed of light c and the distance traveled. A time coincidence or coincidence event is identified if the time difference between the arrivals of signals in a pair of oppositely disposed detectors is within the coincidence time window τ. In conventional PET, the coincidence detection time window τ is wide enough so that an annihilation event occurring anywhere within the object will produce annihilation gamma photons reaching their respective detectors within the coincidence window. Coincidence time windows of 4.5-12 nsec are common for conventional PET, and are largely determined by the time resolution capabilities of the detectors and electronics.
As illustrated in
Time-of-flight (TOF) positron emission tomography (PET) (“TOF-PET”) is based on the measurement of the difference Δt between the detection times of the two gamma photons arising from the positron annihilation event. This measurement allows the annihilation event to be localized along the LOR with a resolution of about 75-120 mm FWHM, assuming a time resolution of 500-800 ps (picoseconds). Though less accurate than the spatial resolution of the scanner, this approximate localization is effective in reducing the random coincidence rate and in improving both the stability of the reconstruction and the signal-to-noise ratio (SNR), especially when imaging large objects. Thus, in TOF-PET, the “TOF” coordinate, Δt, is stored together with s and φ.
Local tomography refers to the ability to reconstruct an image from truncated or incomplete projection data. Such cases typically arise when the object to be imaged is relatively large, the entire object is not sufficiently measured, and the projection data are truncated due to a small detector size and resulting small field of view (FOV). Here, the truncated projection data is in the radial direction as the small FOV cannot directly detect gamma photons from outside the FOV. In this case, because the distribution of activity extends beyond a region of interest (ROI) that is targeted by the FOV, some of this activity from outside the FOV is received by the detector and represents artifact data. However, in such cases, a small ROI can be sufficiently measured, and the small ROI is able to be exactly reconstructed.
Because TOF-PET data includes Δt information in addition to LOR information, the present inventors have investigated the possibility of truncating data acquisition in the time direction as opposed to the radial direction which would be useful in significantly high count rate acquisitions, such as 82Rb PET imaging scans.
The present invention provides a method of TOF-PET image reconstruction using time-truncated TOF-PET projection data. Such time-truncated TOF-PET data is obtained by narrowing the scanner time window to a smaller “time window field of view.” This results in a lower list mode stream transfer rate, which can be useful in high count rate data acquisitions.
The present invention will now be described and disclosed in greater 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.
The present invention provides a method of image reconstruction from TOF-PET projection data that has been truncated in the time direction. Such time truncation is unique to TOF-PET data because only TOF-PET includes a time parameter. The truncation of acquired TOF-PET projection data is obtained by time window narrowing. This results in a lower transfer rate of a list mode data stream, which can be useful in high count rate data acquisitions, such as 82Rb cardiac studies. An advantage of time data truncation as related to radial data truncation is that radial tails information is preserved, thus providing reliable source information for scaling in scatter correction estimation calculations.
It is not apparent from consideration of image reconstruction of radially truncated projection data that local tomography reconstruction of TOF-PET data truncated in the time direction would be possible. However, one variation of a 2D TOF FBP (Filtered Back-Projection) reconstruction algorithm proves this local tomography property. TOF projections are defined as:
p(φ, r, t)=∫{right arrow over (dx)}f({right arrow over (x)})δ(x{right arrow over (x)}{right arrow over (u)}⊥−r)gσ(t−{right arrow over (x)}{right arrow over (u)})
where {right arrow over (u)}⊥=(cos φ, sin φ): {right arrow over (u)}=(−sin φ, cos φ), φ is the azimuthal angle, r is the radial coordinate, and t is the time coordinate. The TOF profile is Gaussian:
The FBP formula (wherein projections are convolved in both the time and radial directions) states:
From this, it can be seen that local tomography with respect to t is possible, in that image reconstruction of a specific local region requires only the TOF bins that contribute to that specific local region. For example, uniform truncation of TOF data will result in exact reconstruction of a central circular region. The local tomography property is also extended to 3D reconstruction, as there is exact TOF rebinning applied for each t separately.
Thus, in an iterative reconstruction algorithm such as OS-EM, one can expect that a local tomography region will be reconstructed exactly. Further, while the extent of reconstruction outside of the local tomography region at the moment is not totally clear, it can be expected that an image outside of the time window FOV can be reconstructed to some degree.
The following references are relevant to this description and are fully incorporated herein by reference:
This application claims benefit under 35 U.S.C. §119(e) to copending provisional application Ser. No. 60/972,894 filed Sep. 17, 2007, which is incorporated herein by reference in its entirety.
| Number | Date | Country | |
|---|---|---|---|
| 60972894 | Sep 2007 | US |
