SPECT gamma camera

• US 6,943,355 B2
• Filed: 08/08/2001
• Issued: 09/13/2005
• Est. Priority Date: 08/10/2000
• Status: Expired due to Fees

• Current Assignee: Ultraspect Ltd.
• Original Assignee: Ume Technologies Co., Ltd.
• Inventors: Ohana, Israel, Shwartz, Shoulamit C.
• Primary Examiner(s): Hannaher, Constantine
• Application Number:
  US10/333,947
• Publication Number:
  US 20030208117A1
• Time in Patent Office: 1,497 Days
• Field of Search: 250/363.02, 250/363.04, 378/4, 378/21
• US Class Current: 250/363.04
• CPC Class Codes: A61B 6/037 Emission tomography
• CPC Class Codes: G01T 1/1644 using an array of optically separate scintillation elements permitting direct location of scintillations G01T1/1645 takes precedence

Abstract

A method and apparatus of obtaining and reconstructing an image of a portion of a body, administered by a radiopharmaceutical substance, by using Single-photon emission computerized tomography (SPECT) for determination of functional information thereon. The method comprises (a) acquiring gamma ray photons emitted from said portion by means of a detector capable of converting the photons into electric signals, the detector having at least one crystal and allowing said gamma rays having incident angles essentially exceeding 5 degrees and, preferably, exceeding 10 degrees to be detected; (b) processing said electric signals by a position logic circuitry and thereby transforming them into data indicative of positions on said photon detector crystal, where the photons have impinged the detector; and (c) reconstructing an image of a spatial distribution of the pharmaceutical substance within the portion of the body by processing said data and taking into consideration weight values which are functions of angles and, possibly, distances between different elements of the portion of the body and corresponding elements of this position'"'"'s projection on the detector.

Claims

1. A method of obtaining and reconstructing an image of a portion of a body, administered by radiopharmaceutical substance radiating gamma rays, by using SPECT (single photon emission computerized tomography), for determination of functional information thereon, comprising the steps of:

(a) acquiring photons emitted from said portion of the body, by means of a detector capable of converting the photons into electric signals, the detector having at least one crystal and adapted to detect emitted photons having incident angles in the range og 0 to 5 degrees;

(b) processing said electric signals by a position logic circuitry and thereby deriving therefrom data indicative of positions on said photon detector crystal, where the photons have impinged the detector; and

(c) reconstructing an image of a spatial distribution of the pharmaceutical substance within the portion of the body by processing said data in conjunction with weight values, derived from functions of either solid angles or solid angles and distances between different discrete elements of the portion of the body and corresponding discrete elements of the projection of the portion of the body on the detector;

wherein said reconstructing an image by processing said data comprises the steps of;

(a) dividing an area of the detector facing the body into M bins;

(b) dividing the portion of the body into N voxels;

(c) providing a set of values D_i (wherein i=1, . . . , M) reflective of the number of photons acquired by each bin;

(d) constructing a matrix P having matrix elements P_ij of weight values of the voxels of the portion of the body (wherein i=1, . . . , M and j=1, . . . , N), the matrix P setting a relation between each bin of the detector each voxel of the portion of the body;

(e) modeling a relation between said set of values D_i and a set of voxel values V_j of said image and deriving said set of voxel V_j of said image, whereby said spatial distribution of the pharmaceutical substance indicating the functional information on said portion of the body is obtained.

2. The method of claim 1 wherein the step of modeling a relation between said set of values D_i and a set of voxel values V_j of said image and deriving said set of voxel values V_j of said image comprises the step of solving a set of equations $D_i=\sum _{}^{}$

j=1N⁢

 

⁢

Pij⁢

Vj

with respect to each value V_j, (wherein i=1, . . . , M).

3. The method of claim 1 wherein the step of modeling a relation between said set of values D_i and a set of voxel values V_j of said image and deriving said set of voxel values V_j of said image comprises the step of solving a set of equations $D_i=\sum _{}^{}$

j=1N⁢

 

⁢

Pij⁢

Vj+Ei

with respect to each value V_j, wherein E_i is a set of measurement errors.

4. The method of claim 1 wherein the step of modeling a relation between said set of values D_i and a set of voxel values V_j of said image and deriving said set of voxel values V_j of said image comprises the step of solving an optimization problem for estimating a mean value λ

(V_j) of Poisson random variables V_j, using the D_i values measured by the detector.

5. The method of claim 4, wherein estimating the value λ

(V_j) of the Poisson random variables V_j, is carried out by calculating a maximum of the likelihood function $L<br><br>\left[\mathrm{\lambda <br><br>}<br><br>\left(D_i\right)\right]=\underset{i=1}{\overset{M}{\prod <br><br>}}<br><br>\text{\hspace{1em}<br><br>}<br><br>\frac{{\mathrm{e<br><br>}}^{-\mathrm{\lambda <br><br>}<br><br>\left(D_i\right)}<br><br>{\mathrm{\lambda <br><br>}<br><br>\left(D_i\right)}^{D_i}}{D_i!},$

with respect to the unknowns V_j.

6. The method of claim 1 wherein the step of modeling a relation between said set of values D_i and a set of voxel values V_j of said image and deriving said set of voxel values V_j of said image comprises the step of solving a Bayesian optimization problem utilizing a likelihood and penalty functions.

7. The method of claim 6, wherein the Bayesian optimization problem has the general form V=argmax{L[λ

(D_i)]+α

F(V_j,V_k)}, where α

is the weight that is given to the penalty function F.

8. The method of claim 1 wherein the matrix P is matrix in which each of the matrix elements P_ij is a function of an average distance and solid angle at which detector bin having an index i is view from the voxel having an index j.

9. The method of claim 1 wherein the matrix P is a matrix in which each of the matrix elements P_ij is a function of an average solid angle at which a detector bin having an index i is viewed from the voxel having an index j.

10. The method of claim 1 wherein the matrix P is a matrix in which each of the matrix elements P_ij is a function of an average distance and solid angle at which a voxel having an index j is viewed from a detector bin having an index i.

11. The method of claim 1 wherein the matrix P is a matrix in which each of the matrix elements P_ij is a function of an average solid angle at which a voxel having an index j is viewed from a detector bin having an index i.

12. The method of claim 1 wherein the matrix P is a matrix in which each of the matrix elements P_ij is a function of a solid angle and distance at which a detector bin having an index i is viewed from a center of a voxel having an index j.

13. The method of claim 1 wherein the matrix P is a matrix in which each of the matrix elements P_ij is a function of a solid angle at which a detector bin having an index i is viewed from a center of a voxel having an index j.

14. The method of claim 1 wherein the matrix P is a matrix in which each of the matrix elements is a function of a solid angle and distance at which a center of a voxel having an index j is viewed from a detector bin having an index i.

15. The method of claim 1 wherein the matrix P is a matrix in which each of the matrix elements P_ij is a function of a solid angle at which a center of a voxel having an index j is viewed from a detector bin having an index i.

16. The method of claim 1 wherein the modeling of the relation between said set of values D_i and a set of voxel values V_j of said image, is biased by an attenuation effect of the patient body.

17. The method of claim 16 wherein the matrix P is a matrix in which each of the matrix elements P_ij is influenced by the attenuation density that exists on the path of the photons that emanate at voxel j and arrive at bin i, for each detector position k.

18. The method of claim 16 wherein the matrix P is a matrix in which each of the matrix elements P_ij is influenced by the attenuation density that exists on the path of the photons that emanate at voxel j and arrive at bin i, for each detector position k and for each energy peak of the radio-pharmaceutical substance that is used.

19. The method of claim 1, further comprising the step of providing a a collimator between said portion of body and said detector, and wherein the matrix P is a matrix in which each of the matrix elements P_ij is multiplied by the relative effective area of bin i associated with the solid angle between voxel j and bin I.