Imaging system and method using scattered and diffused radiation

• US 5,070,455 A
• Filed: 11/22/1989
• Issued: 12/03/1991
• Est. Priority Date: 11/22/1989
• Status: Expired due to Fees

• Current Assignee: Singer Imaging Incorporated
• Original Assignee: Singer Imaging Incorporated
• Inventors: Singer, Jerome R., Zubelli, Jorge P., Grunbaum, Francisco A., Naparst, Harold L., Latham, Geoffrey, Kohn, Philip D., Couch, John L.
• Primary Examiner(s): Jablon, Clark A.
• Assistant Examiner(s): Brutman, Laura
• Application Number:
  US07/440,088
• Time in Patent Office: 741 Days
• Field of Search: 364/413.19, 364/413.14
• US Class Current: 378/6
• CPC Class Codes: G01N 2021/1785 Three dimensional
• CPC Class Codes: G01N 21/4795 spatially resolved investigating of object in scattering medium in vivo A61B
• CPC Class Codes: G06T 11/006 Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
• CPC Class Codes: G06T 2211/424 Iterative

Abstract

An imaging system generates images of the interior of an object using radiation which is attenuated and scattered by the interior of the object. The object is radiated at a number of points near the exterior of the object, and the radiation emerging from the object is measured at an array of exit points near the exterior of the object. The interior of the object is modelled as an array of volume elements called voxels. Each voxel is modelled as having scattering and attenuation properties, represented by numerical coefficients. The system computes the intensity of radiation that would emerge from the object at the exit points if the interior of the object were characterized by a currently assigned set of values for the scattering and attenuation coefficients. Then, the differences between the measured intensities and the computer intensities are used to compute an error function relates to the magnitude of the errors in the reconstruction. This error function is then minimized using a gradient descent methodology, i.e., the value of the coefficients are modified so as to reduce the value of the error function. The coefficients are repeatedly updated in this manner until the error function falls below a specified threshold value. The final values of the scattering and attenuation coefficients from this process are mapped so as to generate a series of images of the interior of the object depicting the attenuation and scattering characteristcs of the object'"'"'s interior.

Claims

1. A system for generating images corresponding to the internal characteristics of an object, comprising:

radiating means for generating a directed beam of radiation, said beam being directed at distinct times at each of a multiplicity of input positions around an object;

sensor means arrayed around at least a portion of said object for measuring the intensity of radiation emerging from said object at a multiplicity of sensor positions;

data storage means, coupled to said sensor means, for storing a multiplicity of sets of measured intensity values, each set of intensity values corresponding to the intensities of radiation measured by said sensor means while said radiating means is generating a beam directed at distinct ones of said input positions;

reconstruction means, coupled to said data storage means, for reconstructing a composite image of a portion of the interior of said object, said reconstructing means computing, using said stored intensity values, an array of image data values comprising radiation attenuation and radiation scattering coefficients which characterize both radiation and scattering characteristics of said portion of the interior of said object.

2. The imaging system of claim 1, wherein said reconstruction means includes means for repeatedly (A) computing intensity values comprising the intensity of radiation that would emerge from said object at said sensor positions while said object is radiated at at least selected ones of said input positions, were said object characterized by said computed attenuation and scattering coefficients, (B) evaluating differences between said computed and measured intensity values, and (C) adjusting said computed attenuation and scattering coefficients, until the difference between said computed and measured intensity values meet predefined criteria.

3. A method of imaging the interior of an object, the steps of the method comprising:

generating a directed beam of radiation, said beam being directed at distinct times at each of a multiplicity of input positions around an object;

measuring the intensity of radiation emerging from said object at a multiplicity of sensor positions;

reconstructing a composite image of a portion of the interior of said object, including computing, using said measured intensity values, an array of image data values comprising radiation attenuation and radiation scattering coefficients which characterize both radiation and scattering characteristics of said portion of the interior of said object.

4. The imaging method of claim 3, wherein said reconstructing step repeatedly (A) computing intensity values comprising the intensity of radiation that would emerge from said object at said sensor positions while said object is radiated at at least selected ones of said input positions, were said object characterized by said computed attenuation and scattering coefficients, (B) evaluating differences between said computed and measured intensity values, and (C) adjusting said computed attenuation and scattering coefficients, until the differences between said computed and measured intensity values meet predefined criteria.

5. A system for generating images corresponding to the internal characteristics of an object, comprising:

radiating means for generating a directed beam of radiation, said beam being directed at distinct times at each of a multiplicity of input positions around an object;

sensor means arrayed around at least a portion of said object for measuring the intensity of radiation emerging from said object at a multiplicity of sensor positions;

data storage means, coupled to said sensor means, for storing a multiplicity of sets of measured intensity values, each set of intensity values corresponding to the intensities of radiation measured by said sensor means while said radiating means is generating a beam directed at distinct ones of said input positions;

reconstruction means, coupled to said data storage means, for reconstructing a composite image of at least a portion of the interior of said object, said reconstructing means including;

(A) modeling means for modeling the interior of said object as an array of voxels, each voxel being characterized by an attenuation coefficient and a plurality of scattering coefficients, said coefficients corresponding to portions of radiation entering said voxel which are absorbed and scattered by said voxel;

said modeling means including;

coefficient denoting means for denoting current values for said attenuation and scattering coefficients for all of said voxels in said model of the interior of said object, including means for assigning initial values to said coefficients;

forward computing means for computing a set of intensity values comprising the intensity of radiation that would emerge from said object at said sensor positions while said object is radiated at at least selected ones of said input positions, were said object characterized by said current values for said attenuation and scattering coefficients as denoted by said coefficient denoting means; and

updating means, coupled to said forward computing means, said data storage means and said coefficient denoting means, for evaluating differences between said computed and measured intensity values, and for generating new values for said attenuation and scattering coefficients which reduce the differences between said computed and measured intensity values; and

(B) control means, coupled to said forward computing and updating means, including means for performing a successive approximation process in which said forward computing and updating means repeatedly compute intensity values and then update said attenuation and scattering coefficients until the differences between said computed and measured intensity values meet predefined criteria.

6. An imaging system as set forth in claim 5, said reconstruction means including image generation means for generating an image corresponding to at least a subset of said attenuation and scattering coefficients.

7. An imaging system as set forth in claim 5, said updating means including:

gradient means for computing a predefined function of said differences between said computed and measured intensity values and for computing a multi-dimensional slope corresponding to partial differentials of said predefined function with respect to each of said coefficients;

means for generating new values for said coefficients by adjusting the current values of said coefficients by amounts corresponding to said multi-dimensional slope.

8. An imaging system as set forth in claim 5, said updating means including:

gradient means for computing a predefined function corresponding to the sum of the squares of said differences between said computed and measured intensity values, and for computing a multi-dimensional slope corresponding to partial differentials of said predefined function with respect to each of said coefficients;

means for generating new values for said coefficients by adjusting the current values of said coefficients by amounts corresponding to said multi-dimensional slope.

9. A system for generating images corresponding to the internal characteristics of an object, comprising:

light generating means for generating a directed beam of light, said beam being directed at distinct times at each of a multiplicity of input positions around an object;

sensor means arrayed around at least a portion of said object for measuring the intensity of light emerging from said object at a multiplicity of sensor positions;

data storage means, coupled to said sensor means, for storing a multiplicity of sets of measured intensity values, each set of intensity values corresponding to the intensities of light measured by said sensor means while said light generating means is generating a beam directed at distinct one of said input positions;

reconstruction means, coupled to said data storage means, for reconstructing a composite image of at least a portion of the interior of said object, said reconstructing means including;

(A) modeling means for modeling the interior of said object as an array of voxels, each voxel being characterized by an attenuation coefficient and a plurality of scattering coefficients, said coefficients corresponding to portions of light entering said voxel which are absorbed and scattered by said voxel;

said modeling means including;

coefficient denoting means for denoting current values for said attenuation and scattering coefficients for all of said voxels in said model of the interior of said object, including means for assigning initial values to said coefficients;

forward computing means for computing a set of intensity values comprising the intensity of radiation that would emerge from said object at said sensor positions while said object is radiated with light at at least selected ones of said input positions, were said object characterized by said current values for said attenuation and scattering coefficients as denoted by said coefficient denoting means; and

updating means, coupled to said forward computing means, said data storage means and said coefficient denoting means, for evaluating differences between said computed and measured intensity values, and for generating new values for said attenuation and scattering coefficients which reduce the differences between said computed and measured intensity values; and

(B) control means, coupled to said forward computing and updating means, including means for performing a successive approximation process in which said forward computing and updating means repeatedly compute intensity values and then update said attenuation and scattering coefficients until the differences between said computed and measured intensity values meet predefined criteria.

10. An imaging system as set forth in claim 9, said reconstruction means including image generation means for generating an image corresponding to at least a subset of said transmission, reflection and scattering coefficients.

11. An imaging system as set forth in claim 9, said updating means including:

gradient means for computing a predefined function of said differences between said computed and measured intensity values and for computing a multi-dimensional slope corresponding to partial differentials of said predefined function with respect to each of said coefficients;

means for generating new values for said coefficients by adjusting the current values of said coefficients by amounts corresponding to said multi-dimensional slope.

12. An imaging system as set forth in claim 9, said updating means including:

gradient means for computing a predefined function corresponding to the sum of the squares of said differences between said computed and measured intensity values, and for computing a multi-dimensional slope corresponding to partial differentials of said predefined function with respect to each of said coefficients;

means for generating new values for said coefficients by adjusting the current values of said coefficients by amounts corresponding to said multi-dimensional slope.

15. A method of imaging the interior of an object as set forth in claim 10, said updating step including the steps of:

computing a predefined function of the differences between said computed and measured intensity values,

computing a multi-dimensional slope corresponding to partial differentials of said predefined function with respect to each of said coefficients; and

computing new values for said coefficients by adjusting the current values of said coefficients by amounts corresponding to said multi-dimensional slope.

13. A method of imaging the interior of an object, the steps of the method comprising:

generating a directed beam of radiation directed at distinct times at each of a multiplicity of input positions around an object;

measuring the intensity of radiation emerging from said object at a multiplicity of sensor positions and storing sets of measured intensity values, each set of intensity values corresponding to said measured intensities while said object is being radiated by a beam directed at distinct ones of said input positions;

reconstructing a composite image of at least a portion of the interior of said object, said reconstructing step including;

(A) modeling the interior of said object as an array of voxels, each voxel being characterized by an attenuation coefficient and a plurality of scattering coefficients, said coefficients corresponding to portions of radiation entering said voxel which are absorbed and scattered by said voxel;

said modeling step including the steps of;

denoting current values for said attenuation and scattering coefficients for all of said voxels;

computing a set of intensity values comprising the intensity of radiation that would emerge from said object at said sensor positions while said object is radiated at at least selected ones of said input positions, were said object characterized by said current values for said attenuation and scattering coefficients;

evaluating differences between said computed and measured intensity values, and

updating said current values of said coefficients by computing new values for said attenuation and scattering coefficients which reduce the differences between said computed and measured intensity values; and

(B) performing a successive approximation process by repeatedly performing said computing, evaluating and updating steps until the differences between said computed and measured intensity values meet predefined criteria.

14. A method of imaging the interior of an object as set forth in claim 13, said reconstructing step including generating an image corresponding to at least a subset of said attenuation and scattering coefficients.

16. A method of imaging the interior of an object as set forth in claim 19, said updating step including the steps of:

computing a predefined function corresponding to the sum of the squares of the differences between said computed and measured intensity values;

computing a multi-dimensional slope corresponding to partial differentials of said predefined function with respect to each of said coefficients; and

computing new values for said coefficients by adjusting the current values of said coefficients by amounts corresponding to said multi-dimensional slope.