AMENDED CLAIMS
[received by the International Bureau on 11 June 2003 (11.06.03); original claims 1-39 replaced by new claims 1-37 (6 pages)]
We claim: 1. A method for obtaining a compositional image of an object comprising obtaining a series of images of said object by transmitting pulsed electromagnetic radiation in a frequency range between about 100 GHz to 10 THz through the entire object from a plurality of angles, detecting changes in the electromagnetic radiation following transmission through said object for each of said plurality of angles using an optical probe beam comprising a chirped laser pulse, and constructing a plurality of imaged slices of the three dimensional object using the detected changes in the transmitted THz radiation.
2. The method according to claim 1 wherein the step of detecting comprises optically detecting said changes in said electromagnetic radiation.
3. The method according to claim 1 wherein said detecting changes in the electromagnetic radiation comprises receiving both transmitted and scattered electromagnetic radiation pulses following transmitting said electromagnetic pulse through said object for each of said plurality of angles.
4. The method according to claim 1 wherein said detecting changes in the electromagnetic radiation following transmission through said object for each of said plurality of angles comprises detecting spectroscopic dielectric information.
5. The method according to claim 4 further comprising using said spectroscopic dielectric information to identify specific materials in said objects.
6. The method according to claim 1 wherein said object is a biological tissue.
7. A method for performing T-ray computed tomography on an object, the method comprising the steps of: (a) emitting a THz pulse and a chirped laser optical probe pulse; (b) transmitting the THz pulse in a first path through the object and transmitting the optical probe pulse in a second path not through the object; (c) modulating the optical probe pulse with the transmitted THz pulse in an electro- optic crystal to create a modulated optical pulse; (d) detecting the modulated optical pulse with a two-dimensional CCD imaging system and storing information corresponding thereto; (e) repeating steps (a)-(d) for a plurality of delay times between the THz pulse and the optical pulse sufficient to characterize a full waveform of the THz pulse; (f) repeating steps (a)-(e) for a plurality of projection angles relative to the object. ; (g) constructing at least one image of the object using data collected during steps (a)-(f).
8. The method according to claim 7 wherein the step of constructing at least one image of the object comprises constructing said at least one image representing a cross sectional slice of said object.
9. The method according to claim 7 wherein in step (g) the image constructed is a 3D representation of the object.
10. The method of claim 9 comprising using a filtered backprojection algorithm in step (g) to construct the three-dimensional image.
11. The method of claim 7 wherein step (g) comprises measuring a phase and amplitude of the transmitted THz pulses and reconstructing the complex refractive index of the object using said measured phase and amplitude as the input to an algorithm which inverts the Radon transform to reconstruct the complex refractive index of the object.
12. The method of claim 11 wherein there are transmitted a plurality of THz pulses and wherein the phase measurement comprises a phase delay and is estimated by interpolating the transmitted THz pulses and interpolating a reference THz pulse, cross-correlating the interpolated reference pulse with the transmitted pulses and using a timing of a peak of the cross correlation as the estimated phase delay.
13. The method of claim 11 where the amplitude and phase are estimated at a plurality of frequencies by Fourier transforming the time domain THz pulses and calculating a phase and amplitude of the complex Fourier coefficients at each given frequency.
14. The method of claim 7 comprising detecting the modulated optical pulse using a spectrometer and a CCD camera.
15. The method of claim 9 wherein said three-dimensional image comprises voxels and wherein a Fourier frequency-dependent refractive index and absorption coefficient for each voxel in the three-dimensional image is constructed using a Fourier transform of the information.
16. The method of claim 7, wherein the THz pulse comprises an expanded THz pulse, the optical probe pulse comprises an expanded optical probe pulse, and detecting the modulated optical probe pulse comprises using a two-dimensional charge coupled device (CCD) camera and the information corresponding thereto includes diffraction information.
17. The method of claim 16, wherein step (e) comprises repeating steps (a)-(d) for a plurality of delay times between the THz pulse and the optical pulse sufficient to characterize a full wavelength of the THz pulse. 22
18. The method of claim 17, wherein step (g) comprises constructing the three-dimensional image using a mathematical algorithm selected from the group consisting of: an algorithm based on a linearization of the wave equation, an algorithm that inverts the non-linear wave equation using iterative finite difference techniques, an iterative algorithm based on the contrast source inversion algorithm, or an algorithm based upon reconstruction of 2D profiles of the object using Fresnel Diffraction.
19. The method of claim 18, wherein the linearization algorithm comprises a Born or Rytov approximation.
20. The method of claim 18, wherein the iterative finite difference technique comprises a propagation-back propagation (PBP) algorithm.
21. The method of claim 7, further comprising using a dynamic subtraction technique to increase the signal to noise ratio.
22. A system for performing T-ray tomography, the system comprising: means for providing a THz pulse and an optical pulse; means for directing the THz pulse in a first path through an object; means for directing the optical pulse through a second path not through the object; an electro-optic crystal at a point of conversion of the first path and second path, for modulating the optical pulse with the THz pulse to create a modulated optical pulse; means for detecting the modulated optical pulse; means for storing information including diffraction information relating to the detected modulated optical pulse; means for obtaining the stored information for the temporal profile of the THz pulse for a plurality of pixels along an x-axis of the object and a y-axis of the object for a plurality of projection angles; means for rotating the object relative to the THz pulse along a z-axis to obtain the plurality of projecting angles; means for constructing a three-dimensional image of the object using the information obtained for the plurality of pixels for the plurality of projection angles.
23. The system of claim 22 wherein the means for obtaining the stored information for a THz pulse comprises means for providing a chirped optical pulse.
24. The system of claim 22 wherein the means for detecting the modulated optical pulse comprises a spectrometer and a CCD camera.
25. The system of claim 23 wherein the means for obtaining the stored information for a full wavelength comprises a delay stage for providing a plurality of delays between the THz pulse and the optical probe pulse and the means for detecting the modulated optical pulse comprises one of a two-dimensional CCD imaging system or single photo-detector.
26. The system of claim 25, wherein the means for providing a THz pulse and an optical pulse comprises means for providing an expanded THz pulse and an expanded optical pulse, the means for obtaining the stored information for the temporal profile comprises a delay stage for providing a plurality of delays between the THz pulse and the optical probe pulse, and the means for detecting the modulated optical pulse comprises a two-dimensional CCD imaging system.
27. The system according to claim 22 further comprising an optical system placed in said first path at a point between said object and said point of conversion, said optical system comprising at least a focusing lens and a limiting aperture means at a focal point of the focusing lens to attenuate scattered radiation components of said expanded THz pulse.
28. A system for performing T-ray tomography, the system comprising: means for simultaneously providing a THz pulse and an optical pulse; means for directing the THz pulse in a first path through an object; means for directing the optical pulse through a second path not through the object; a conversion point for said first and said second paths a photoconductive antenna at said point of conversion of the first path and second path for generating a photocurrent proportional to the instantaneous THz field; means for detecting the generated photocurrent; means for storing information relating to the generated photocurrent; means for obtaining the stored information for the temporal profile of the THz pulse for a plurality of pixels along an x-axis of the object and a y-axis of the object for a plurality of projection angles; means for rotating the object relative to the THz pulse along a z-axis to obtain the plurality of projecting angles; means for constructing a three-dimensional image of the object using the information obtained for the plurality of pixels for the plurality of projection angles.
29. A method for performing T-ray imaging on an object, the method comprising the steps of: (a) emitting a THz pulse and an optical probe pulse; (b) transmitting the THz pulse in a first path through the object and transmitting the optical probe pulse in a second path not through the object; (c) modulating the optical probe pulse with the transmitted THz pulse in an electro-optic crystal to create a modulated optical pulse; (d) detecting the modulated optical pulse and storing information corresponding thereto using a spectrometer and a CCD camera.; (e) repeating steps (a)-(d) from a plurality of locations relative to the object sufficient to characterize the object; (f) constructing a three-dimensional image of the object using the information collected during steps (a)-(e).
30. The method of claim 29 comprising using a filtered backprojection algorithm in step (f) to construct the three-dimensional image.
31. The method of claim 29 wherein said three dimensional image comprises voxels the method further comprising using a Fourier transform of the information to construct a frequency-dependent refractive index and absorption coefficient for each voxel in the three-dimensional image.
32. The method of claim 31 further comprising means for performing a dynamic subtraction technique to improve signal-to-noise ratio in the detected information relating to the modulated expanded optical pulse.
33. A method for providing a two-dimensional profile image of an object, the method comprising the steps of: (a) emitting an expanded THz pulse and an expanded optical probe pulse; (b) transmitting the expanded THz pulse in a first path through the object and transmitting the expanded optical probe pulse in a second path not through the object; (c) modulating the expanded optical probe pulse with the transmitted expanded THz pulse in an electro-optic crystal to create a modulated expanded optical pulse; (d) detecting the modulated expanded optical pulse with a charge coupled device (CCD) camera and storing information corresponding thereto, including diffraction information; (e) constructing a two-dimensional profile image of the object using a mathematical algorithm based upon the time-reversal of the Huygen-Fresnel Diffraction integral.
34. A method for obtaining a compositional image of an object comprising obtaining a series of images of said object by transmitting pulsed electromagnetic radiation in a frequency range between about 100 GHz to 10 THz through the entire object from a plurality of angles, detecting changes in the electromagnetic radiation following transmission through said object for each of said plurality of angles using an optical probe beam, and constructing a three dimensional image of the three dimensional object using the detected changes in the transmitted THz radiation wherein said constructing of said three dimensional image comprises assuming a linearized propagation model and applying a filtered backprojection algorithm.
35. The method of claim 34 wherein the filtered backprojection algorithm is used to perform an inverse Radon equation to construct the three dimensional image.
36. The method of claim 34 wherein the linearized propagation model comprises an approximation relating the transmitted electromagnetic radiation and the detected changes in the electromagnetic radiation following transmission through said object by a line integral.
37. The method of claim 36 wherein the line integral has the form:
r- iω n{r)
Pd (ω,θ, l) = Pi(ω) ex ■dr
L(θ,l) where: Pd is the Fourier transform of the detected THz radiation at frequency ω, a projection angle θ, and a horizontal offset from the axis of rotation /, Pt is the Fourier component of the transmitted THz radiation at the same frequency, L is a straight line between the source and detector, c is the speed of light, and n(τ) = n(r) + ik(r) is the unknown complex refractive index of the sample.