Three dimensional image reconstruction for a positron emission tomograph

• US 5,331,553 A
• Filed: 04/15/1992
• Issued: 07/19/1994
• Est. Priority Date: 04/15/1992
• Status: Expired due to Fees

• Current Assignee: Ugm Medical Systems, Inc. (Koninklijke Philips N.V.)
• Original Assignee: Ugm Medical Systems, Inc. (Koninklijke Philips N.V.)
• Inventors: Lewitt, Robert M., Muehllehner, Gerd
• Primary Examiner(s): Thomas, Joseph
• Application Number:
  US07/869,523
• Time in Patent Office: 825 Days
• Field of Search: 364/413.13, 364/413.24, 364/413.14, 364/413.15, 364/413.16, 364/413.19, 364/413.22, 378/24, 378/25, 378/26, 378/27, 395/119, 395/124, 250/363.03, 250/363.04, 250/363.07, 250/363.02, 250/363.1, 250/369
• US Class Current: 250/363.02
• CPC Class Codes: G01T 1/2985 In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis); using external radiation sources A61B6/02
• CPC Class Codes: G06T 17/00 Three dimensional [3D] modelling, e.g. data description of 3D objects

Abstract

Images are reconstructed using multi-slice rebinning of raw data to create a stack of two-dimensional data sets. A three-dimensional image is reconstructed slice-by-slice from the data sets and the rebinned data is axially filtered to reduce blurring resulting from the rebinning, the filtering being performed either before or after reconstruction of the three-dimensional image.

Claims

1. A computer implemented method for reconstructing images, comprising the steps of:

providing an array of detectors, the detectors being disposed around and having therebetween a transversely and axially extending field of view, the detectors having an axial acceptance angle such as to enable sensing of both transverse and oblique rays;

detecting in said field of view both transverse and oblique rays;

converting the detected rays to data in electrical form;

collecting and assembling said data representative of the detected rays;

rebinning said data to create a stack of two-dimensional data sets representative of individual transverse slice images, whereby said slice images signify in the aggregate a three-dimensional image;

reconstructing slice-by-slice from said data sets a three-dimensional image; and

axially filtering said rebinned data to reduce blurring resulting from said rebinning.

2. A computer implemented method in accordance with claim 1, wherein said data sets are representative of sinograms.

3. A computer implemented method in accordance with claim 1, wherein said step of axial filtering is performed on said data sets before said step of reconstructing the three-dimensional image.

4. A computer implemented method in accordance with claim 1, wherein said step of axial filtering is performed on the three-dimensional image after said step of reconstructing the three-dimensional image.

5. A computer implemented method for reconstructing images using a device of the type having an array of detectors, the detectors having therebetween a transversely and axially extending field of view and having an axial acceptance angle such as to enable sensing of both transverse and oblique rays, comprising the steps of:

detecting in said field of view both transverse and oblique rays;

converting the detected rays to data in electrical form;

collecting and assembling said data representative of the detected rays;

rebinning said data to create a stack of two-dimensional data sets representative of individual transverse slice images, whereby said slice images signify in the aggregate a three-dimensional image;

reconstructing slice-by-slice from said data sets a three-dimensional image; and

axially filtering said rebinned data to reduce blurring resulting from said rebinning;

said step of rebinning said data comprising the further steps of;

for each oblique ray, identifying among said transverse data sets the images intersected by said ray; and

applying to each of the data sets associated with the images intersected by said ray an increment representative of said ray.

6. A computer implemented method in accordance with claim 5, wherein the data sets are sinograms, said step of applying an increment to each of the data sets for a given line comprising applying the increment to each of the sinograms at the same two-dimensional location thereon.

7. A computer implemented method for reconstructing images using a device of the type having an array of detectors, the detectors having therebetween a transversely and axially extending field of view and having an axial acceptance angle such as to enable sensing of both transverse and oblique rays, comprising the steps of:

detecting in said field of view both transverse and oblique rays;

converting the detected rays to data in electrical form;

collecting and assembling said data representative of the detected rays;

rebinning said data to create a stack of two-dimensional data sets representative of individual transverse slice images, whereby said slice images signify in the aggregate a three-dimensional image;

reconstructing slice-by-slice from said data sets a three-dimensional image; and

axially filtering said rebinned data to reduce blurring resulting from said rebinning, said step of axial filtering being performed on the three-dimensional image after said step of reconstructing the three-dimensional image and comprising the further steps of;

selecting a plurality of axially spaced locals z'"'"' along the longitudinal axis of the field of view, the spacing of the locales corresponding to the thickness of said slice images;

establishing at each locale z'"'"' a grid defining a plurality of positions (x, y, z'"'"') in the transverse plane through the locale z'"'"';

simulating for each position a point source and rays through said point source;

rebinning and storing data representative of said simulated rays, whereby said rebinned and stored data provides a spread function relating to said simulated point source; and

applying said spread functions to the data representative of said three-dimensional image in order to reduce the axial blurring present in said three-dimensional image.

8. A computer implemented method in accordance with claim 7, wherein the axial field of view has a central transverse plane and the selected locales z'"'"' are disposed from the central transverse plane of the axial field of view to the end of the field of view, the longitudinal axis of the field of view being disposed at a central location within the grid and the grid having a transverse (x,y) dimension approximately equal to the transverse dimension of the field of view.

9. A computer implemented method for reconstructing images using a device of the type having an array of detectors, the detectors having therebetween a transversely and axially extending field of view and having an axial acceptance angle such as to enable sensing of both transverse and oblique rays, comprising the steps of:

detecting in said field of view both transverse and oblique rays;

converting the detected rays to data in electrical form;

collecting and assembling said data representative of the detected rays;

rebinning said data to create a stack of two-dimensional data sets representative of individual transverse slice images, whereby said slice images signify in the aggregate a three-dimensional image;

reconstructing slice-by-slice from said data sets a three-dimensional image; and

axially filtering said rebinned data to reduce blurring resulting from said rebinning, said step of axial filtering being performed on said data sets before said step of reconstructing the three-dimensional image and comprising the further steps of;

selecting a plurality of axially spaced locales z'"'"' along the longitudinal axis of the field of view, the spacing of the locales corresponding to the thickness of said slice images;

simulating at each locale z'"'"' a point source and rays through said point source, said simulated rays being uniformly distributed in the transverse plane through the locale z'"'"' and in oblique planes through the locale z within the range defined by the axial acceptance angle of the detectors;

selecting from among said simulated rays those which intersect a pair of detectors;

rebinning, accumulating and storing data representative of said selected rays, whereby said data provides a spread function relating to the locale z'"'"'; and

applying said spread functions to the data representative of the respective sinograms in order to reduce the axial blurring present in said stack of two-dimensional data sets.

10. A computer implemented method in accordance with claim 9, wherein the axial field of view has a central transverse plane and the selected locales z'"'"' are disposed from the central transverse plane (z'"'"'=o) of the axial field of view to the end of the field of view.

11. A computer implemented method for reconstructing images in positron emission tomography using a device of the type having an array of detectors, the detectors having therebetween a transversely and axially extending field of view and having an axial acceptance angle such as to enable sensing of both transverse and oblique rays, comprising said steps of:

detecting in the field of view both transverse and oblique rays;

converting the detected rays to data in electrical form;

collecting and assembling said data representative of the detected rays;

rebinning said data to create a stack of two-dimensional sinograms representative of individual transverse slice images, whereby said slice images signify in the aggregate a three-dimensional image;

reconstructing slice-by-slice from said sinograms a three-dimensional image; and

axially filtering said rebinned data to reduce blurring resulting from said rebinning.

12. A computer implemented method in accordance with claim 11, wherein said step of axial filtering is performed on said sinograms before said step of reconstructing the three-dimensional image.

13. A computer implemented method in accordance with claim 11, wherein said step of axial filtering is performed on the three-dimensional image after said step of reconstructing the three-dimensional image.

14. A computer implemented method for reconstructing images in positron emission tomography using a device of the type having an array of detectors, the detectors having therebetween a transversely and axially extending field of view and having an axial acceptance angle such as to enable sensing of both transverse and oblique rays, comprising the steps of;

detecting in said field of view both transverse and oblique rays;

converting the detected rays to data in electrical form;

collecting and assembling said data representative of the detected rays;

rebinning said data to create a stack of two-dimensional sinograms representative of individual transverse slice images, whereby slice images signify in the aggregate a three-dimensional image;

reconstructing slice-by-slice from said sinograms a three-dimensional image; and

axially filtering said rebinned data to reduce blurring resulting from said rebinning, said step of

axially filtering said rebinned data comprising the further steps of;

selecting a plurality of axially spaced locales z'"'"' along the longitudinal axis of the field of view, the spacing of the locales corresponding to the thickness of said slice images;

establishing at each locale z'"'"' a grid defining a plurality of positions (x, y, z'"'"') in the transverse plane through the locale z'"'"';

simulating for each position a point source and rays through said point source;

rebinning and storing data representative of said simulated rays, whereby said rebinned and stored data provides a spread function relating to said simulated point source; and

applying said spread functions to the data representative of said three-dimensional image in order to reduce the axial blurring present in said three-dimensional image.

15. A computer implemented method in accordance with claim 14, wherein the axial field of view has a central transverse plane and the selected locales z'"'"' are disposed from the central transverse plane of the axial field of view to the end of the field of view, the longitudinal axis of the field of view being disposed at a central location within the grid and the grid having a transverse (x,y) dimension approximately equal to the transverse dimension of the field of view.

16. A computer implemented method for reconstructing images in positron emission tomography using a device of the type having an array of detectors, the detectors having therebetween a transversely and axially extending field of view and having an axial acceptance angle such as to enable sensing of both transverse and oblique rays, comprising the steps of:

detecting in said field of view both transverse and oblique rays;

converting the detected rays to data in electrical form;

collecting and assembling said data representative of the detected rays;

rebinning said data to create a stack of two-dimensional sinograms representative of individual transverse slice images, whereby said slice images signify in the aggregate a dimensional image;

reconstructing slice-by-slice form said sinograms a three-dimensional image; and

axially filtering said rebinned data to reduce blurring resulting from said rebinning, said step of

axially filtering said rebinned data comprising the further steps of;

selecting a plurality of axially spaced locales z'"'"' along the longitudinal axis of the field of view, the spacing of the locales corresponding to the thickness of said slice images;

simulating at each locale z'"'"' a point source and rays through said point source, said simulated rays being uniformly distributed in the transverse plane through the locale z'"'"' and in oblique planes through the locale z'"'"' within the range defined by the axial acceptance angle of the detectors;

selecting from among said simulated rays those which intersect a pair of detectors;

rebinning, accumulating said storing data representative of said selected rays, whereby said data provides a spread function relating to the locale z'"'"'; and

applying said spread functions to the data representative of the respective sinograms in order to reduce the axial blurring present in said stack of sinograms.

17. A computer implemented method in accordance with claim 16, wherein the axial field of view has a central transverse plane and the selected locales z'"'"' are disposed form the central transverse plane (z'"'"'=o) of the axial field of view to the end of the field of view.