Diffraction tomography systems and methods with arbitrary source and detector array surfaces

• US 4,598,366 A
• Filed: 01/25/1984
• Issued: 07/01/1986
• Est. Priority Date: 11/12/1982
• Status: Expired due to Term

• Current Assignee: Schlumberger Technology Corporation (Schlumberger NV)
• Original Assignee: Schlumberger Technology Corporation (Schlumberger NV)
• Inventors: Devaney, Anthony J.
• Primary Examiner(s): Smith, Jerry
• Assistant Examiner(s): Jablon, Clark A.
• Application Number:
  US06/573,540
• Time in Patent Office: 888 Days
• Field of Search: 364/414, 364/420, 364/421, 364/422, 364/400, 73/602, 250/253, 250/256, 250/269, 128/660, 367/14, 378/4, 378/901
• US Class Current: 700/90
• CPC Class Codes: G01S 15/8977 using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
• CPC Class Codes: G06T 11/005 Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
• 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/421 Filtered back projection [FBP]
• CPC Class Codes: Y10S 128/916 Ultrasound 3-D imaging
• CPC Class Codes: Y10S 378/901 Computer tomography program or processor

Abstract

Diffraction tomography systems and methods are disclosed for the reconstruction of physical properties of two and three-dimensional objects from information collected by detectors located along an arbitrary surface in space. The sources of energy investigating the object may likewise be located along an arbitrary surface. Optimal filtered backpropagation techniques and preprocessing required for such reconstructions are disclosed for systems utilizing plane wave sources, and for systems using fan beam cylindrical or spherical wave insonification. The preprocessing and filtered backpropagation techniques are shown to simplify for the cases of circular boundaries, and other separable boundaries within the Helmholtz equation.

Claims

1. A method for the tomographic reconstruction of an object which has diffracted waves of energy propagating through it, said waves of energy being selected from a group consisting of sonic and electromagnetic waves, comprising:

obtaining signals which are a function of at least one of the phase and amplitude of the diffracted propagating waves after said waves have interacted with said object, wherein the sources of said energy waves are located along an arbitrary curvilinear surface, other than a single straight line and the detectors of said diffracted waves are located along an arbitrary curvilinear surface other than a single straight line; and

converting said signals by means of a filtered backpropagation technique into an array representing a partial reconstruction of the object.

2. A method according to claim 1 further comprising:

preprocessing said obtained signals before said converting step wherein said signals converted by said converting step are preprocessed signals, and wherein said preprocessing includes transforming said obtained signals with a general transform.

3. A method according to claim 2 wherein:

said sources of energy are plane wave sources; and

said general transform transforms a determined scattered field into a plane wave scattering amplitude for input into said filtered backpropagation technique.

4. A method according to claim 3 wherein:

said preprocessed signals are filtered according to the filtered backpropagation technique ##EQU30## wherein s_o is a unit propagation vector of the incident plane waves, s is a unit vector which may assume all directions, r is the vector from a coordinate system origin to a particular point of reconstruction in said coordinate system, γ

(s;

s_o) is the plane wave scattering amplitude, χ and

χ

_o are respectively the angles formed by s and s_o with a fixed reference direction, and O(r;

χ

_o) is the partial reconstruction of the object.

5. A method according to claim 4 wherein said object is two-dimensional further comprising:

summing said partial reconstructions over a plurality of view angles to obtain a full two-dimensional reconstruction of said object.

6. A method according to claim 5 wherein:

said detectors are located along a circular boundary.

7. A method according to claim 6 wherein:

said energy sources are located along a circular boundary.

8. A method according to claim 5 wherein:

said detectors are located along a boundary of weak curvature.

9. A method according to claim 5 wherein:

said detectors are located along a boundary which approximates two infinite parallel lines.

10. A method according to claim 4 wherein said preprocessing further comprises:

transforming said obtained signals from analog to digital form;

Fourier transforming said digital signals in time; and

determining the scattered field from said Fourier tansformed signals.

11. A method according to claim 2 wherein:

said sources of energy are cylindrical or spherical wave sources;

said general transform transforms a scattered field into a cylindrical wave response; and

said cylindrcial wave response is slant stacked to provide the plane wave scattering amplitude for input into said filtered backpropagation technique.

12. A method according to claim 11 wherein:

said preprocessed signals are filtered according to the filtered backpropagation technique ##EQU31## wherein s_o is a unit propagation vector of the incident waves, s is a unit vector which may assume all directions, r is the vector from the origin of a coordinate system to a particular point of reconstruction in said coordinate system, γ

(s;

s_o) is the plane wave scattering amplitude, χ and

χ

_o are respectively the angles formed by s and s_o with a fixed reference direction, and O(r;

χ

_o) is the partial reconstruction of the object.

13. A method according to claim 12 wherein said object is two-dimensional further comprising:

summing said partial reconstructions over a plurality of view angles to obtain a full two-dimensional reconstruction of said object.

14. A method according to claim 13 wherein said preprocessing further comprises:

transforming said obtained signals from analog to digital form;

Fourier transforming said digital signals in time; and

determining the scattered field from said Fourier tansformed signals.

15. A method according to claim 14 wherein:

said detectors are located along a circular boundary.

16. A method according to claim 15 wherein:

said energy sources are located along a circular boundary.

17. A method according to claim 14 wherein:

said detectors are located along a boundary of weak curvature.

18. A method according to claim 14 wherein:

said detectors are located along a boundary which approximates two infinite parallel lines.

19. A method according to claim 1 further comprising:

preprocessing said obtained signals before said converting step wherein said signals converted by said converting step are preprocessed signals.

20. A method according to claim 19 wherein said sources of energy are cylindrical or spherical wave sources and said proprocessing comprises:

transforming said obtained signals from analog to digital form;

Fourier transforming said digital signals in time; and

determining the scattered field from said Fourier transformed signals.

21. A method according to claim 20 further comprising:

summing said partial reconstructions to obtain a full reconstruction by integrating over the point locations of said sources.

22. A method according to claim 20 wherein:

said object is two-dimensional;

said sources and detectors are located along circular boundaries; and

said scattered field is filtered according to a filtered backpropagation technique ##EQU32## where ##EQU33## and wherein O(r,β

) is a partial reconstruction of said object, r is the vector from the origin of a coordinate system to a particular point of reconstruction in the coordinate system, β

is the angle formed by R_o with an arbitrary reference where R_o is the cylindrical or spherical wave source location, σ

is the angle made by r'"'"' with an arbitrary reference direction where r'"'"' is a detector point location, Ψ

(σ

,β

) is the scattered field, χ and

χ

_o are respectively the angles formed by s and s_o with fixed reference directions where s is a unit vector which may assume all directions and s_o is a unit propagation vector of the incident waves, and F_R is a function given by ##EQU34## where H_n is the n'"'"'th order Hankel function of the first kind, and k is the wave numbeer, and R is the radius of the detector boundary and R_o is the radius of the source boundary.

23. A method according to claim 20 wherein:

said object is two-dimensional;

said detectors are located along a boundary which approximates two infinite parallel lines; and

said scattered field is filtered according to a filtered backpropagation technique ##EQU35## where ##EQU36## and wherein O(r,l_o) is a partial reconstruction of said object, r is the vector from the origin of a coordinate system to a particular point of reconstruction in the coordinate system, l_o is the location of R_o along a straight boundary line where R_o is the cylindrical or spherical wave source location, l'"'"' is the location of r'"'"' along the detector surface boundary where r'"'"' is the detector point location, Ψ

(l'"'"';

l_o) is the scattered field, χ and

χ

_o are respectively the angles formed by s and s_o with fixed reference directions where s is a unit vector which may assume all directions and s_o is a unit propagation vector of the incident waves, n'"'"' is the unit outward normal vector to the detecting surface, and n_o is the unit outward normal vector to the source surface.

24. A method according to claim 20 wherein:

said object is two-dimensional;

said detectors are located along a boundary of weak curvature; and

said scattered field is filtered according to a filtered backpropagation technique ##EQU37## where ##EQU38## and wherein O(r,l_o) is a partial reconstruction of said object, r is the vector from the origin of a coordinate system to a particular point in the coordinate system, l_o is the location of R_o along a straight boundary line where R_o is the cylindrical or spherical wave source location, l'"'"' is the location of r'"'"' along the detector surface boundary where r'"'"' is the detector point location, Ψ

(r'"'"';

R_o) is the scattered field, χ and

χ

_o are respectively the angles formed by s and s_o with fixed reference directions where s is a unit vector which may assume all directions and s_o is a unit propagation vector of the incident waves, n'"'"' is the unit outward normal vector to the detecting surface, and n_o is the unit outward normal vector to the source surface.

25. A method for the tomographic reconstruction of an object comprising:

directing, from energy sources located along an arbitrary curvilinear surface, waves of energy toward said object;

detecting as a function of time, with detectors located along an arbitrary curvilinear surface other than a single straight line, the fields produced by said waves;

producing signals which are functions of said detected fields;

preprocessing said produced signals; and

filtering said preprocessed signals according to a filtered backpropagation technique to generate an array representing a partial reconstruction of said object.

26. A method according to claim 25 wherein:

said preprocessed signals are filtered according to the filtered backpropagation technique ##EQU39## wherein s_o is a unit propagation vector of the incident waves, s is a unit vector which may assume all directions, r is the vector from the origin of a coordinate system to a particular point of reconstruction in said coordinate system, γ

(s;

s_o) is the plane wave scattering amplitude, χ and

χ

_o are respectively the angles formed by s and s_o with a fixed reference direction, and O(r;

χ

_o) is the partial reconstruction of the object.

27. A method according to claim 26 wherein said object is a two-dimensional object, further comprising:

summing said partial reconstructions over a plurality of view angles to obtain a full two-dimensional reconstruction of said object.

28. A method according to claim 26 wherein said sources of energy are plane wave sources, and said preprocessing step comprises:

transforming said produced signals from analog to digital form;

Fourier transforming said digital signals in time;

determining the scattered field from said Fourier transformed signals; and

determining the plane wave scattering amplitude by transforming said scattered field with a general transform.

29. A method according to claim 28 wherein:

said detectors and sources are located along circular boundaries.

30. A method according to claim 28 wherein:

said detectors and sources are located along boundaries of weak curvature.

31. A method according to claim 28 wherein:

said detectors and sources are located along boundaries which approximate two infinite parallel lines.

32. A method according to claim 26 wherein said sources of energy are cylindrical or spherical wave sources and said preprocessing step comprises:

transforming said produced signals from analog to digital form;

Fourier transforming said digital signals in time;

determining the scattered field from said Fourier transformed signals;

transforming the scattered field into a cylindrical wave response according to a general transform; and

determining the plane wave scattering amplitude by slant stacking said cylindrical wave response.

33. A method according to claim 32 wherein:

said detectors and sources are located along circular boundaries.

34. A method according to claim 32 wherein:

said detectors and sources are located along boundaries of weak curvature.

35. A method according to claim 32 wherein

said detectors and sources are located along boundaries which approximate two infinite parallel lines.

36. A method according to claim 25 wherein said sources of energy are cylindrical or spherical wave sources and said preprocessing comprises:

transforming said produced signals from analog to digital form;

Fourier transforming said digital signals in time; and

determining the scattered field from said Fourier transformed signals.

37. A method according to claim 36 further comprising:

summing said partial reconstructions by integrating over the point locations of said sources of energy to obtain a full reconstruction of said object.

38. A method according to claim 36 wherein:

said sources and detectors are located along circular boundaries; and

said scattered field is filtered according to the filtered backpropagation technique ##EQU40## where ##EQU41## and wherein O(r,β

) is a partial reconstruction of said object, r is the vector from the origin of a coordinate system to a particular point located in said coordinate system, β

is the angle formed by R_o with an arbitrary reference where R_o is the cylindrical or spherical wave source location, σ

is the angle made by r'"'"' with arbitrary reference direction where r'"'"' is a detector point location, Ψ

(σ

,β

) is the scattered field, χ and

χ

_o are respectively the angles formed by s and s_o with fixed reference directions where s is a unit vector which may assume all directions and s_o is a unit propagation vector of the incident waves, and F_R is a function given by ##EQU42## where H_n is the n'"'"'th order Hankel function of the first kind, and k is the wavenumber. and R is the radius of the detector boundary and R is the radius of the source boundary.

39. A method according to claim 36 wherein:

said detectors are located along a boundary which approximates two infinite parallel lines; and

said scattered field is filtered according to the filtered backpropagation technique ##EQU43## where ##EQU44## and wherein O(r,l_o) is a partial reconstruction of said object, r is the vector from the origin of a coordinate system to a particular point in the coordinate system, l_o is the location of R_o along a straight boundary line where R_o is the cylindrical or spherical wave source location, l'"'"' is the location of r'"'"' along the detector surface boundary where r'"'"' is the detector point location, Ψ

(l_o ;

l'"'"') is the scattered field, χ and

χ

_o are respectively the angles formed by s and s_o with fixed reference directions where s is a unit vector which may assume all directions and s_o is a unit propagation vector of the incident waves, n'"'"' is the unit outward normal vector to the detecting surface, and n_o is the unit outward normal vector to the source surface.

40. A method according to claim 36 wherein:

said detectors are located along a boundary of weak curvature; and

said scattered field is filtered according to a filtered backpropagation technique ##EQU45## where ##EQU46## and wherein O(r,l_o) is a partial reconstruction of said object, r is the vector from the origin of a coordinate system to a particular point in said coordinate system, l_o is the location of R_o along a straight boundary line where R_o is the cylindrical or spherical wave source location, l'"'"' is the location of r'"'"' along the detector surface boundary where r'"'"' is the detector point location, Ψ

(r'"'"';

R_o) is the scattered field, χ and

χ

_o are respectively the angles formed by s and s_o with fixed reference directions where s is a unit vector which may assume all directions and s_o is a unit propagation vector of the incident waves, n'"'"' is the unit outward normal vector to the detecting surface, and n_o is the unit outward normal vector to the source surface.

41. A system for the tomographic reconstruction of an object comprising:

energy sources located along an arbitrary curvilinear surface for insonifying said object with waves of energy;

detectors of energy located along an arbitrary curvilinear surface other than a single straight line, for detecting the fields produced by the interaction of said source energy with said object and for producing signals which are functions of said detected fields;

preprocessing means for preprocessing said produced signals; and

processing means for filtering and backpropagating said preprocessed signals according to a filtered backpropagation technique to generate an array representing a partial reconstruction of said object.

42. A system according to claim 41 wherein:

said processing means includes means for filtering and backpropagating said preprocessed signals according to the filtered backpropagation technique ##EQU47## wherein s_o is a unit propagation vector of the incident waves, s is a unit vector which may assume all directions, r is the vector from the origin of a coordinate system to a particular point of reconstruction in said coordinate system, γ

(s;

s_o) is the plane wave scattering amplitude, χ and

χ

_o are respectively the angles formed by s and s_o with a fixed reference direction, and O(r;

χ

_o) is the partial reconstruction of the object.