Process for obtaining three-dimensional images of an object and its application to the tomography of an organ

• US 4,539,639 A
• Filed: 02/04/1983
• Issued: 09/03/1985
• Est. Priority Date: 02/17/1982
• Status: Expired due to Fees

• Current Assignee: Commissariat a L'Energie Atomique (AES Corporation)
• Original Assignee: Commissariat a L'Energie Atomique (AES Corporation)
• Inventors: Soussaline, Francoise, Le Coq, Gerard
• Primary Examiner(s): Smith, Jerry
• Assistant Examiner(s): Woo, Lonis
• Application Number:
  US06/463,811
• Time in Patent Office: 942 Days
• Field of Search: 364/414, 250/363 S, 378/901
• US Class Current: 378/13
• CPC Class Codes: A61B 6/037 Emission tomography
• 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: Y10S 128/922 including image analysis
• CPC Class Codes: Y10S 378/901 Computer tomography program or processor

Abstract

A process and appparatus for obtaining three-dimensional images of an object with an application of the process and the device to the tomography of an organ. The device includes a detector or mobile rays which rotates about the object and mechanisms which provide references for the angular position of the detector with respect to the reference plane. The detector contains detection areas which receive the radiation and are identified by their number and emit signals, the amplitude of which depends on the intensity of the received rays. Also included in the device is an iterative processing machine which operates by convergent processing of the signal amplitude values with the processing mechanism being connecting to visualization devices in order to display the images of the object in various cross-section planes after each iteration.

Claims

1. A process for obtaining three-dimensional images of an object from radiation emitted by a tracer contained in said object, said radiation being received by certain detection areas of a planar detector, designated by numbers n in the plane of said detector, the latter performing a rotary movement around the object, so that it constantly remains tangential to an imaginary cylinder of revolution (C), whose axis corresponds to an axis of the object, the angular positions of the planar detector being designated by angles θ

_p between a reference plane (PR) passing through the axis of the object and a plane passing through the axis of the object and perpendicular to the plane of the detector, the object being assumed as broken down into elementary cells for the emission of radiation designated by their coordinates in a fixed reference mark (i, j, k), axis k of said mark being parallel to the axis of the object, the radiation received by certain detection areas being those which are emitted by the elementary cells of the object in a direction perpendicular to the plane of the detector, the detection areas supplying, for each angular position θ

_p, electrical signals whose amplitudes have values dependent on the intensities of the radiation received, said amplitude values being stored in order to obtain the images of the object, in the sectional planes, of ordinates k perpendicular to the axis of the object, wherein, for obtaining images of the object in the sectional planes, it comprises;

producing on the basis of the stored amplitude values, for each sectional plane of ordinates k, for all the angular positions θ

_p, and for all the detection areas of number n receiving the radiation, correction values of said amplitude values, whereby these correction values are obtained on initialization, on the basis of simulated amplitude values, corresponding to a simulated image of the object in which it is assumed that the distribution of the tracer is uniform, and from stored correction factors relating to the attenuation (u) of the radiation in the object, as a function of coordinates (i, j) of each elementary cell and relating to the geometry of each detector;

determining and storing for each sectional plane, for all the angular positions (θ

_p) and for all the detection areas of numbers n which receive the radiation, the values of the respective differences between the values of the amplitudes of the signals supplied by these areas and the correction values, the images corresponding to the values of the amplitudes of the signals supplied by the detection areas being available for the sectional planes for which the values of the differences comply with a predetermined convergence criterion, then, when the values of the differences do not comply with the predetermined criterion, the process then comprises;

producing other correction values of the signal amplitude values, for the sectional planes of ordinates k, for all the angular positions θ

_p and for all the detection areas of numbers n which receive the radiation, said other correction values being obtained from values of stored differences, on the basis of said correction factors and on the basis of stored filtering coefficients (λ

) of said other correction values;

performing for each sectional plane of ordinates k, for all the angular positions (θ

p) and for all the detection areas receiving the radiation, additions between respectively the stored amplitude values and the other correction values, the values of the results of these additions being stored, in order to supply the corresponding images of the object in the sectional planes or ordinates k or to obtain new correction values on the basis of the results of said additions until these results converge.

2. A process according to claim 1, wherein the correction values are obtained for each sectional plane of ordinates k, for all the angular positions θ

p and for all the areas receiving radiation, by multiplication of the values of the simulated amplitudes or resulting from additions, with the correction factor values.

3. A process according to claim 2, wherein the other correction values are obtained for each sectional plane of ordinates k, for all the angular positions θ

p and for all the detection areas receiving the radiation, by the multiplication of said differences with said correction factors and with said filtering coefficients.

4. A process according to claim 3, wherein it comprises, for each sectional plane of ordinates k, for all the angular positions θ

p to be determined and stored, prior to obtaining the correction values, of the coordinates (i, j) of a mask of the space external of the images of the desired sections, of determining the correspondences between the numbers of the detection areas n and the coordinates (i, j) of the elementary cells located in the mask and then comparing the numbers of the detection areas defined by the mask with the numbers of the detection areas receiving the radiation.

5. A process according to claim 4, wherein it comprises using a microcomputer connected to at least one memory and to image display means, for controlling the rotation of the detector, the marking of the angular positions θ

p around the object, the markings of the detection areas and the coordinates i, j, k and for carrying out all the operations of storing, correcting and displaying images.

6. An apparatus for obtaining three-dimensional images of an object on the basis of radiation emitted by a tracer contained in said object, comprising a planar detector having detection areas designated by numbers n and able to supply signals respectively dependent on the intensities of the radiation received in a direction perpendicular to the plane (P) of the detector, said detector performing a rotary movement around the object, so that it remains constantly tangential to an imaginary cylinder (C), whose axis corresponds to an axis of the object, means for controlling and marking the angular position (θ

p) of the detector, said positions being marked between a reference plane (P_R) passing through the axis of the object and a plane passing through the axis of the object and perpendicular to the plane of the detector, the object being fictitiously broken down into elementary cells for emitting the radiation and designated by their coordinates in a fixed reference marking (i, j, k), the axis k of said marking being parallel to the axis of the object, the radiation received by certain detection areas designated by their numbers n being those emitted by the elementary cells of the object in a direction perpendicular to the plane of the detector, the detection areas supplying, for each angular position θ

p, electrical signals whose amplitudes have values which are dependent on the intensities of the radiation received, said amplitude values making it possible to obtain images of the object in sectional planes of ordinates k, perpendicular to the axis of the object, wherein it comprises iterative processing means by convergence of the values of the amplitudes of the signals from the detection areas, which have received radiation, said processing means being connected to display means for displaying the images corresponding to each sectional plane of ordinates k after each iteration.

7. An apparatus according to claim 6, wherein the processing means comprise:

storage means which, for all the ordinates k of the sectional planes and all the angular positions θ

p of the detector, make it possible to store;

the values of the amplitudes of the signals supplied by the detection areas of numbers n, which have received radiation,

the coordinates (i, j) of the elementary cells which may have emitted radiation to said detection areas of numbers n,

reference values of the signal amplitudes, said values being simulated or obtained at the end of a preceding processing operation by the device of signal values corresponding to the detection areas,

correction values to be applied to the reference values,

at least one amplitude value filtering coefficient;

first means for multiplying the reference amplitude values by the correction factors, the values of the results of these multiplications being recorded by storage means;

means for comparing the signal amplitude values supplied by the detection areas with reference values, the values of the results of these differences being recorded in the storage means;

means for testing the convergence of the values of the differences, as a function of a predetermined convergence criterion, display means being connected to the testing means for displaying the images of the sectional planes of ordinates k, when the convergence criterion is respected;

second means for the multiplication of the values of the differences with the correction factors and the filtering coefficient, said means being used when the convergence criterion is not respected, the values of the results of these multiplications being recorded in the storage means;

means for the addition of the values of the multiplication results to the reference amplitude values, the testing means being connected to the display means for controlling the display of images of the sectional planes of ordinates k on the basis of the results of the additions, when a convergence criterion of the values of the addition results is respected, said addition means being connected to the first multiplication means and to the storage means, so that when the convergence criterion of the values of the results of the additions is not respected, the reference values supplied to said second multiplication means are values of the results of the additions.

8. An apparatus according to claim 7, wherein it also comprises means for controlling the storage means so that the latter record, for each sectional plane ordinate k, and for all the angular positions θ

p, masking coordinates (i, j) defining a limiting contour for the reconstruction of the object, and means for comparing the masking coordinates with the coordinates (i, j) of the elementary cells which have emitted radiation, so that account is only taken of the detection areas, whose numbers n correspond to elementary cells located in the limiting masking content.

9. An apparatus according to claim 8, wherein the correction factors recorded in the storage means relate to the geometry of the detection areas, at the solid angles under which are seen the detection areas receiving radiation, on the basis of corresponding elementary cells of the object, and to the attenuation of the radiation as a function of the distance separating an elementary cell from the masking contour.