Method for eliminating scattered gamma-rays and reconstructing image, and gamma camera apparatus

Abstract: A method of eliminating scattered gamma rays and collecting an image comprises the steps of detecting a radioisotope (RI) given to a living body, to form an RI distribution image, and eliminating scattered gamma rays generated in the living body or a gamma camera (1) from the RI distribution image. In this case, the energy spectrums of gamma rays incident to individual positions in the gamma camera (1) in association with pixels of the RI distribution image are collected at the same time as or during the collection of the RI distribution image. A ratio of the scattered gamma rays in respect to the RI distribution image is acquired from the energy spectrum for each position, and a scattered gamma ray component is eliminated from the RI distribution image for each pixel, by multiplying the ratio of the scattered gamma rays to the RI distribution image.

Description: The present invention relates to a scattered γ-ray eliminating and image a gamma camera apparatus which detects a radioisotope given to a living body and accumulated or deposited in an organ, to form an RI distribution image, and eliminates scattered γ-rays in the living body and scattered γ-rays in a gamma camera from this RI distribution image, by setting a window. Nuclear medical apparatuses have been used to detect a radioisotope given to a living body and accumulated or deposited in an organ using a gamma camera and form an image of the two-dimensional distribution of the radioisotope for diagnostic use. In this system, scattered γ-rays (scattered gamma-rays) are generated in the living body or in the gamma camera (e.g., collimator or NaI scintillator). Since the scattered γ-rays are not necessary for diagnostic information, they should be eliminated. There are two methods known of eliminating scattered γ-rays from an image acquired by a gamma camera (scintillation camera) or the like as disclosed in J. Nucl. Med. 14; 67-72, 1972, J. Nucl. Med. 25; 490-494, 1984, J. Nucl. Med. 29; 195-202, 1988, and IEEE. Tran. Nucl. Science. NS32. 786-793, 1985. According to the first method, a window a0 is set for a photoelectric peak Pl in the relation of the amplitude versus an energy spectrum E as shown in Fig. 10. Images within the window a0 are collected and a window b0 is set for a Compton scattered component C0 at the same time or as the next sequence. Based on a photoelectric image A(x, y) and a scattered-ray image S(x, y) acquired from the windows a0 and b0, a process of A(x, y)-R·S(x, y) is executed as a method of eliminating the scattered rays, where radioisotope is a constant representing a predicted ratio of scattered rays included in a photoelectric absorbing peak P1. According to the second scattered-ray eliminating method, as scattered γ-rays have a distribution dependent on a position (x, y), they are taken as an image more accurately than in the first method. In the relation of the amplitude versus the energy spectrum E as shown in Fig. 11, therefore, a window having a sufficiently narrow window width ΔE is sequentially shifted from El to Ep having a peak. At this time, images E(x, y) in individual steps El-Ep are collected and their respective energy images dependent on the γ detected positions, at which γ-rays are detected by the gamma camera, are formed, whereby scattered components for the individual positions are acquired. Those conventional methods of eliminating scattered rays, however, have the following shortcomings. The first method eliminates scattered γ-rays using only one value of the predicted constant radioisotope. If the distribution of scattered rays varies depending on γ detected positions, however, this distribution differs from the actual physical phenomena. Therefore, the proper scattered-ray elimination can not be performed for those positions, thus providing less accurate images. According to the second method, improving the accuracy in eliminating scattered rays requires making ΔE smaller and collection of a large number of images. This takes a significant amount of time in collecting the images. Further, in order to eliminate nuclide scattered rays having two or more photoelectric peaks position signal, a greater number of images should be collected, thus requiring a longer time for collecting images. Accordingly, it is an object of this invention to provide a scattered γ-ray eliminating method and an image collecting method and a gamma camera apparatus which can accurately and simply eliminate scattered γ-rays in a living body or a gamma camera even if these scattered γ-rays depend on the detection positions in the gamma camera. According to the present invention, there is provided a scattered γ-ray eliminating method and an image collecting method of detecting γ-rays irradiated from a radioisotope given to a living body to form an RI distribution image and eliminating scattered γ-rays generated in the living body or a gamma camera from this RI distribution image by setting a window, in which RI distribution images are collected and energy spectrums of γ-rays incident to individual positions in the gamma camera at the same time or during the collection are collected in association with the RI distribution image, a ratio of scattered γ-rays is obtained from the energy spectrum for each position, and a scattered γ-ray component is eliminated from the RI distribution image for each pixel based on the ratio. A gamma camera apparatus comprises a gamma camera, having a scintillator and a plurality of photomultipliers, for outputting γ-rays irradiated from a radioisotope given to a living body as position signals x and y and an energy signal Z proportional to the energy of γ-rays, a first memory for storing an RI image at an address designated by the position signals x and y when the energy signal Z from the gamma camera is within a predetermined value range, a second memory for storing the energy signal Z from the gamma camera assembly for each position as an energy spectrum, and a controller for controlling writing of the energy signal Z into the second memory and writing of the RI image into the first memory. The gamma camera collects the energy spectrum E depending on the position (x, y) at the same time as or following collection of RI images whereby the whole position-dependent spectrum image is accurately processed, the ratio of scattered rays is acquired from this spectrum, and the collected images are multiplied by this ratio. As a result, scattered γ-rays can be accurately and simply eliminated for each detected position in the gamma camera. Therefore, the speed for processing the RI image is increased and the resolution of the collected images is improved. In a case where nuclides that generate γ-rays of two or more energies or two types of more of nuclides of different energies are simultaneously collected and scattered rays are eliminated from the collected image, in particular, this scattered-ray elimination can be executed accurately and simply. This invention can be more fully understood from the following detailed description when taken in conjunction with the accompanying drawings, in which: Fig. 1 is a schematic block diagram of a gamma camera apparatus to which a method of eliminating scattered γ-rays and collecting an image according to one embodiment of the present invention is applied; Fig. 2 is a block diagram of a gamma camera according to another embodiment of this invention; Fig. 3 is a diagram illustrating an RI image obtained by a gamma camera; Figs. 4A through 4C are schematic diagrams illustrating energy spectrums at individual positions of the RI image shown in Fig. 3; Fig. 5 is a flowchart for scattered γ-ray elimination and image collection executed by the gamma camera apparatus shown in Fig. 1; Fig. 6 is a schematic diagram illustrating the relation between the energy and amplitude per channel according to this invention; Fig. 7 is a diagram illustrating an energy spectrum E and a count value of photoelectrons; Fig. 8 is a diagram illustrating a flowchart for acquiring the area NPA of a photoelectric peak; Fig. 9 is a flowchart for scattered γ-ray elimination and image collection executed by the gamma camera apparatus shown in Fig. 2; and Figs. 10 and 11 are diagrams for explaining a conventional method of eliminating scattered γ-rays and collecting an image. Referring to Fig. 1, a scintillation camera 1 (called a gamma camera) has a scintillator and a plurality of photomultipliers, and is connected to an A/D converter 2. This A/D converter 2 converts position signals x and y and an energy signal Z from the scintillation camera 1 into digital signals. The A/D converter 2 has its output terminal connected to input terminal of an image memory controller 3, a window circuit 5, an x-y address selector 6 and a peak discriminator 7. The window circuit 5 issues a write enable command S1 to the image memory controller 3 when the energy signal Z from the A/D converter 2 falls within the width between the upper limit WU and lower limit WL of a given window set by a main CPU 9. The x-y address selector 6 selects x and y addresses in accordance with the position signals x and y from the A/D converter 2. Upon reception of the write enable command S1 from the window circuit 5, the image memory controller 3 serving as a control circuit, adds "1" to the content of an x-y corresponding memory address on an image data memory 4 in response to the address signal from the x-y address selector 6 and writes image data in the memory 4. The position (x, y) for the energy signal Z input to the window circuit 5 is discriminated by the x-y address selector 6. Further, the peak of the signal Z is discriminated by the peak discriminator 7, "1" is added to the content of the memory address corresponding to an energy associated with the magnitude of the spectrum signal Z corresponding to the position (x, y) of a spectrum data memory 8, and an energy spectrum (X, Y, e) corresponding to the position is collected at the same time as collection of images. In other words, the energy spectrum Z at the position (X, Y) in the scintillation camera is stored in the form of a combination of energy spectrums (X, Y, e) in the data memory 8. The window circuit 5 sets the window width corresponding to about 20-30% of the energy of a photoelectric peak according to a control signal from the main CPU 9. The image memory 4 stores an RI distribution image i.e., γ-ray distribution image P (x, y) obtained at the coordinates (x, y) of the gamma camera to which there is input the input energy within the window width set by the window circuit 5. The spectrum data memory 8 stores an energy spectrum for each position in the field of view of the gamma camera, and stores spectrums E (x, y, e) collected within a sufficiently wider window width in order to presume a ratio of the scattered components in respect to the γ-rays within the photoelectric peak to stored in the image memory 4. Fig. 3 illustrates a collected γ-ray distribution image (hereinafter referred to as collected image or P (x, y)). Figs. 4A through 4C illustrate the relation between the energy spectrum collected in association with each pixel position (i-1, j;i, j;i+1, j) of a collected image in Fig. 3 (hereinafter referred to as energy spectrum or E (x, y, e)) and the collected γ-ray distribution image (collected image or P (x, y)). Referring to Figs. 1 through 4, there will now be described the scattered γ-ray eliminating method, the image collecting method and the gamma camera. An RI image signal including scattered γ-rays from the gamma camera 1, i.e., the position signals x and y and energy signal Z, are converted into digital signals by the A/D converter 2. The position signals x and y are input to the x-y address selector 6 while the energy signal Z is input to the window circuit 5 and peak discriminator 7. The window circuit 5, which receives a command from the main CPU 9, collects an image P (x, y) in a certain window width, e.g., a width of 20-30% of the photoelectric peak energy. At the same time, the window circuit 5 collects the energy spectrums E (x, y, e) corresponding to the position on an image with the full window width (step A). Then, the main CPU 9 determines whether or not the count of the collected energy spectrum is large (step B). If the count is large, the position of the photoelectric peak is detected for E (x, y, e) (step C). Further, the area of the photoelectric peak, NPA (Net Peak Area), within WL ≦ e ≦ WU with respect to the window position (upper limit position WU, lower limit position WL) set at the time of collecting P (x, y) is obtained from the following equation (1). NPA (x, y) = ∫ WU WL E (x, y, e)de - B (x, y) where B (x, y) = {E (x, y, WU) + E (x, y, wl)} × (WU-WL)/2. That is, the photoelectric peak NPA (x, y) and scattered ray component B (x, y) are obtained (step D). Then, the RI image P (x, y) is multiplied by the ratio of the scattered rays as expressed by the following equation (2). P (x, y) × NPA (x, y)/{NPA (x, y) + B (x, y)} Through execution of the equation (2), an image P′ having scattered rays eliminated can be acquired (step E). In step B, if the count of the energy spectrum is smaller than a given value, the image memory controller 3 executes a filtering process using the following equation (3) with the filtering coefficient al to reduce the statistic noise of the amplitude. Then, after a more accurate energy spectrum is obtained by shading off the position information (step F), the process following the step C is executed. Therefore, the collected image P (x, y) is replaced with the image P′ (x, y) free of a scattered component which deteriorates the quantization of the image. According to this embodiment, as described above, the energy spectrum E dependent on the position (x, y) is accurately processed as the whole position-dependent image at the same time as the image collection, and the ratio of the scattered rays is obtained from this spectrum. The collected image P (x, y) is multiplied by this ratio. As a result, the scattered rays can be accurately and simply eliminated, thus increasing the image processing speed and improving the resolution of the diagnostic image. Further, nuclide that irradiates γ-rays of two or more energies or two or more types of nuclides of different energies can be simultaneously collected as images from which scattered rays can be eliminated in particularly accurate and simple manner. In the steps C-E of eliminating scattered rays from the collected image as shown in Fig. 5, since the scattered components within the energy spectrum are considered only as the shaded portion B of the energy spectrum E (i-l, j, k) in Fig. 4A, or the base background of the energy spectrum E, they are eliminated. However, in clinical examinations, since γ-rays irradiated from RI (radioisotope) distributed in a living body are scattered in the body, the shape of the photoelectric peak is deformed as compared with the one attained when no scattering material exists in the living body. Utilizing this phenomenon, the amount of deformation is measured, given that the measured value is S (x, y), a more accurate NPA can be acquired from the following equation (4). NPA' (x, y) = NPA (x, y) - S (x, y) The above processing will be described more specifically referring to Fig. 6. Given that E (x, y, e) is the energy spectrum in clinic and Einair (x, y, e) is the energy spectrum in the air where no scattering material exists, i.e., the response function of the scintillation camera, then S (x, y) can be obtained from the equation (5) below. S (x, y) = ∫ WU WL [E (x, y, e) - Einair (x, y, e)]de where Einair is integer-multiplied so that, for example, the peak and the corner of a curve, coincide with each other in the range Pc < e < WU. In other words, the above equation can be rewritten as the following equation (6). S (x, y) = ∫ WU WL E (x, y, e) - A·Einair (x, y, e) de Then, A (real number) to minimize the result of the equation (6). A more detailed algorithm for squiring the area NPA of the photoelectric peak in the above manner by this embodiment will be described below, referring to Fig. 7. In Fig. 7, the horizontal scale indicates the energy spectrum E, and K is a main peak energy. The peak discriminator 7 shown in Fig. 7 is designed to collect energy spectrums within the same energy range as the energy window width set by the window circuit 5. In the window circuit 5 at this time is set a window of K ± (Ws + We)/2 for the main peak energy K, where Ws is a preset energy window width and We is an arbitrary energy width set such that it partially overlaps the edge portion of Ws. For instance, there two values are so set as to make (Ws/K) x 100 = approximately 20% and We ≒ approximately 2 to 4 KeV. Further, the coefficients expressed by the following equations are computed. Then, a diagnostic image P (x, y) is multiplied by the ratio of scattered rays using the individual coefficients to predict NPA (x, y), as expressed by the following equation. P (x, y) × Ct - (Ch + Cℓ) × Ws 2We Ct Based on the result of the above computation, P' free of scattered rays is acquired. Fig. 8 illustrates a flowchart for acquiring P' in the above manner. First, a window width of K ± (Ws + We)/2 is set as described above (step a). Then, the diagnostic image P (x, y) and energy spectrum E (x, y, e) are computed (step b) followed by the computation of the individual coefficients Ct (x, y), Ch (x, y) and cℓ (x, y) (step c). Finally, after NAP (x, y) is computed using the equation (10), the scattered ray component is eliminated on the basis of NPA (x, y) to provide the image P' free of the scattered rays (step d). With NPA acquired through the above method, if the amplitude is low, the scattered-ray component B (x, y) can accurately be predicted by specifying the individual coefficients, Ct, Ch and Cℓ. Accordingly, the scattered component can accurately be eliminated and the collecting time can be shortened. In other words, although it is difficult to accurately obtain the scattered-ray component B (x, y), the component B (x, y) can be stably determined since the term (Ch + Cℓ)Ws/2We in the equation (10) can be expressed as B (x, y) ≒ (Ch + Cℓ)Ws/2We. It should however be noted that the above is applicable to a system having a energy resolution of 10 to 15% in terms of a half-width. The second embodiment of the present invention will be described below. There will now be described a case where a gamma camera is applied to a single-photon emission (hereinafter referred to as SPECT). To begin with, P (x, y) represents an collected image by the aforementioned scintillation camera. The gamma camera is rotated 360° or 180° around an object to be examined, and collects a projection image for each n degrees to acquire a projection image P (x, y, ϑ). For instance, if n = 6°, the energy spectrum is collected for each projection image to acquire E (x, y, ϑ) where ϑ represents the angle at which the image collection is done. Accordingly, the scattered ray elimination for each angle ϑ is performed in the above-described procedures shown in Fig. 5 for a projection image P (x, y, ϑ) for each angle ϑ. The same effect as described above can be obtained in an apparatus using a tomographic image. If the amplitude of the energy spectrum is insufficient to provide a sufficient accuracy, it may be assumed that a change distribution of scattered rays does not depend so much on the angle compared with a projection image, and the energy spectrum for each 2·n or 3·n degrees can be used as a mean value. The third embodiment will be described below referring to Fig. 2. The same reference numerals as used in Fig. 1 are used to denote the identical or corresponding elements. The split circuit 15 output, to the image memory controller 3, only position signals x and y which falls in a specific range of x0 < X < x1, y0 < Y < y1 with respect to the A/D converted position signals. The circuit also outputs the energy signal Z to the window circuit 5. Under certain conditions, the position signals x and y in the specific range are stored by the image memory controller 3, at the corresponding X-Y address in the image data memory 4. Under other conditions, the energy spectrum E (x, y, e) of γ-rays set by the split circuit 15 and entered in the range (x, y) (hereinafter referred to as "split"), is collected by a main CPU memory 17. The collection of the scintigram images and the energy spectrum in this device will be described referring to the flow chart illustrated in Fig. 7. First, x and y are set to be 0 < (x, y) < (xmax, ymax). The split circuit 15 is driven by the main CPU/memory 17 for full window with 0 < x < xmax, 0 < y < ymax (step a). The CPU/memory 17 sets the required window width for example, 20% with respect to the main peak in the window circuit 5 (step b), then collects image P (x, y) (step c). The CPU/memory 17 initializes m and n (step d), and sets the window circuit 15 for full window (step e). Further, the CPU/memory 17 sets the split circuit 15 within a predetermined range in accordance with m·Δx ≦ x < (m + 1)·Δx, m·Δy ≦ y < (m + 1)·Δy (step f). Then, the energy spectrum E (x, y, e) is collected and stored in the main CPU/memory 17 (step g). Next m and n are set by the main CPU 17 (step h). When these m and n are set, those m and n within the range which exceeds a certain count on the already collected image P (x, y) are set (step k). If the next m and n are not set, the sequence is terminated (step j). Even if the energy spectrum E (x, y, e) and collected image P (x, y) are not simultaneously collected, therefore, the scattered rays can be eliminated from the final information acquired in the same manner as done in the case of simultaneous collection, except for those of earlier movement tests. In the actual clinic or, SPECT collection, the time can be shortens by setting ϑ greater than the collecting angle n by the energy spectrum E (x, y, e) as per the second embodiment. However, it takes a significant time for collection if the energy spectrum E over the entire image is collected. Alternatively, in a region of interest (hereinafter referred to as ROI) having a count more than a predetermined value may be set in the collected image P (x, y, ϑ) before the energy spectrum is collected, and the write enable command s1 may be issued so as to collect the energy spectrum E only within the ROI. According to the present invention, the energy spectrum E depending on the position (x, y) is collected at the same time as or following the image collection to accurately process the whole image of the spectrum, the ratio of the scattered rays can be acquired from this spectrum, and the collected image is multiplied by this ration, thus ensuring accurate and simple elimination of the scattered rays. This increases the processing speed and improves the resolution of the diagnostic image. Further, if images of a nuclide irradiating γ-rays of two or more energies or two or more types of nuclides with different energies, for example, are collected simultaneously for elimination of scattered rays, the scattered ray eliminating method and image collecting method can be executed particularly in accurate and simple manner.

Claims: A method of eliminating scattered rays and collecting an image, comprising the steps of: detecting γ-rays irradiated from a radioisotope (RI) given to a living body and collecting an RI image in a predetermined energy window; collecting energy spectrums of γ-rays incident to individual positions in a gamma camera (1) following collection of said RI image in association with pixels of said RI image; characterized by collecting said energy spectrums of γ-rays alternative- ly at the time of collection of the RI image; acquiring a ratio of said scattered rays from said energy spectrum for each of said individual positions given by: P (x, y) x NPA (x, y)/{NPA(x, y) + B (x, y)} whereby P (x, y) represents said RI image, NPA (x, y) represents the net peak area of the photoelectric peak, and B (x, y) represents the scattered ray components; and eliminating a scattered ray component from said RI image for each of said pixels using said ratio. A method according to claim 1, characterized in that said RI image collecting step is for collecting said RI image in accordance with a window width corresponding to 20 to 30 % of a peak energy of detected γ-rays. A method according to claim 1, characterized in that said energy spectrum collecting step is for collecting said energy spectrum in a full window state. A gamma camera apparatus comprising: a gamma camera (1), having a scintillator and a plurality of photomultipliers, for detecting γ-rays irradiating from a radioisotope given to a living body and outputting a position signal and an energy signal proportional to energy of said γ-rays, first memory means (4) for storing image data at an address corresponding to said position signal and collecting an RI image when said energy signal from said gamma camera is within a predetermined energy window; characterized by second memory means (8) for storing said energy signal from said gamma camera as an energy spectrum for each position; and control means (9) for acquiring a ratio of said scattered rays from said energy spectrum for said each position and eliminating a scattered ray component from said RI image for each pixel using said ratio given by: P (x, y) x NPA (x, y)/{NPA(x, y) + B (x, y)} whereby P (x, y) represents said RI image, NPA (x, y) represents the net peak area of the photoelectric peak, and B (x, y) represents the scattered ray components. A gamma camera apparatus according to claim 4, characterized in that said first memory means includes window means (5) for subjecting said energy signal to window processing with a predetermined window width. A gamma camera apparatus according to claim 5, characterized in that said window means (5) includes means for subjecting said energy signal to window processing in accordance with a window width corresponding to 20 to 30 % of peak energy of γ-rays. A gamma camera apparatus according to claim 4, 5 or 6, characterized in that said first memory means includes address selecting means (6) for selecting an address according to said position signal and outputting an address signal. A gamma camera apparatus according to anyone of claims 3 to 7, characterized in that said first memory means includes a first memory (4) for storing image data in response to said address signal. A gamma camera apparatus according to anyone of claims 4 to 8, characterized in that said first memory means (4) stores a γ-ray distribution image corresponding to an input energy signal in a predetermined window width, and said second memory means (8) stores an energy spectrum for each position in a field of view of said gamma camera (1). A gamma camera apparatus according to anyone of claims 4 to 9, further comprising means (3) for filtering said position signal and energy signal in accordance with a level of an energy spectrum. A gamma camera apparatus according to anyone of claims 4 to 10, characterized in that said gamma camera (1) is rotated around an object to be examined by a predetermined angle m and outputs a position signal and an energy signal corresponding to a projection image for each n (n < m).

Inventors: ICHIHARA TAKASHI C O INTELLECT, ICHIHARA, TAKASHI, ICHIHARA, TAKASHI, C/O INTELLECTUAL PROPERTY DIV.
Applicants: TOSHIBA KK, KABUSHIKI KAISHA TOSHIBA
Application year: 1990
Main classification: G01T1
IPC codes: G01T1