Patent Grant
7326930
This application claims the benefit of Great Britain application GBP290114 filed on 19 Jan. 2004.
Not Applicable
Not Applicable
This invention relates to apparatus and methods for sensing terahertz radiation, in particular over an area, and to terahertz radiation imaging systems.
The terahertz (THz) region of the electromagnetic spectrum spans the frequency range between the mid-infrared and the millimeter/microwave portion of the spectrum. In this specification, terahertz radiation is considered to be that lying within a frequency range of from 0.01 THz to 100 THz. In this specification we will also refer to light and to optical paths; it is to be understood that this includes non-visible light such as infrared and ultra violet light.
Terahertz systems have shown much promise in recent years (see, for example M. C. Beard, G. M. Turner, and C. A. Schmuttenmaer, J. Phys. Chem. B 106, 7146(2002), X.-C. Zhang, Phys. Med. Biol. 47, 3667 (2002)), however commercialization of this technology has been hampered by a lack of rugged compact imaging systems that allow rapid acquisition of terahertz images.
It has traditionally been very difficult to detect radiation in the terahertz part of the electromagnetic spectrum. Suitable detectors are typically very bulky, they often require cooling with cryogenic fluids (as with bolometric detection), and they are rarely robust and insensitive to vibration. Further more there are significant problems with background noise. Most objects even down to cryogenic temperatures will emit large amounts of thermal radiation in this wavelength band. It is difficult to exclude this noise in detection schemes.
Recently new schemes to detect radiation in the terahertz wavelength region have shown promise—see, for example, Hu BB, Nuss MC Imaging with Terahertz Waves, Optics Letters, 20(16), 1716 (1995). Radiation in this wavelength region induces small changes in the visible wavelength birefringence of electro-optic crystals. If terahertz radiation is provided in to an electro-optic crystal it is possible to detect its field strength by measuring the change in birefringence of the electro-optic crystal. This is known as electro-optic sampling (or EOS). Typically this birefringence change is measured by passing an optical probe pulse through the electro-optic crystal collinearly with a terahertz pulse and measuring the change in polarization of the optical pulse with a balanced diode detection scheme. It is possible by this method to coherently detect the terahertz amplitude and phase of any given terahertz pulse. This method has been employed in both biological and non-biological terahertz imaging.
If one uses ultrashort laser pulses to both generate and detect terahertz pulses this facilitates gating the EOS detection such that the detector is sensitive for only the short space of time when the sample is illuminated by the terahertz pulse. This means that the vast majority of the background radiation is excluded.
Further where the detection is coherent (i.e. the time of arrival of the terahertz pulse can be accurately measured and the entire pulse can be plotted in time) it is possible, in a pixel imaging system, to conduct spectroscopic analysis of each pixel measured in a manner analogous to FTIR.
Generally these schemes measure a single point at a time and scan the terahertz spot over the sample to measure all parts of it. However, there has been a demonstration of a system where the terahertz pulse is spread over a large area electro-optic crystal and is probed by passing an optical probe pulse over the large area of the electro-optic crystal; Lu Z G, Campbell P, Zhang X C, Free-space electro-optic sampling with a high repetition-rate regenerative amplified laser, Appl. Phys. Lett 71(5) 4, Aug. 1997. Detection of the probe pulse occurs by passing that pulse through a polarizer and then on to a CCD chip. Where the optical pulse is detected indicates where terahertz was present on the electro optic crystal.
Terahertz radiation sensors are sometimes known as terahertz probes. One design of probe is described in GB 2,352,512A, together with a number of variants of the basic design. Use of the probe to image a tooth is also described. The contents of the specification are hereby incorporated by reference in their entirety. Broadly speaking the probe of GB'512 uses electro-optic sampling along the lines described above but the device exhibits some drawbacks. One drawback relates to the sensitivity of the probe and more particularly to the size of the electro-optic response to terahertz radiation, which is very small. This, coupled with the relative lack of efficiency of terahertz emitters can result in a need to measure a change in polarization of the order of 1 part in 104 to 1 part in 106. The imaging device described in GB'512 (that is the device of
A system for 3-dimensional imaging of an object using terahertz radiation is described in WO 03/042670. This describes the conventional approach to 3D imaging, that is capturing a plurality of 2D images using a point source of terahertz radiation and a terahertz-imaging sensor, rotating the object through a small angle between each (page 5, lines 1 to 6). All the terahertz probes described in GB'512 also include a terahertz lens at the front end of the probe to, essentially, image the terahertz radiation onto a sensor. Delay control is employed in GB'512 but this is merely to align the timing of the (optically) gated terahertz detector with that of the terahertz emitter.
Other background prior art is described in EP 0841548, EP 0727671, EP 0864857, U.S. Pat. No. 5,894,125, US 2001/038074 (WO 01/77646), US 2003/178584, GB 2350673, and EP 1271115. For example GB 2350673A describes a system for terahertz imaging by measuring the time of flight of a pulse through a sample, optionally Fourier transforming temporal data into the frequency domain in order to obtain a terahertz power spectrum to determine a plurality of images at different frequencies, for example to determine the composition of a structure.
Details of the time of flight technique (which assumes no diffraction) are described in GB'512 with reference to
We will describe a probe system which uses a polarizing beam splitter. This beam splitter is preferably constructed from a birefringent crystal such as calcite or a synthetic material, although other types of polarizing beam splitter may be used.
A large spot size probe beam is inserted in to the beam splitter. The beam then reflects off an internal surface of the beam splitter and passes to the electro-optic crystal. The electro optic crystal is oriented to the polarization of the probe pulse (further details of the alignment of EOS crystals for THz generation and detection are discussed in “Generation of Terahertz Radiation from semiconductors” Thesis—University of Cambridge UK by Alessandra Corchia hereby incorporated by reference). When the terahertz pulse arrives in the EOS crystal it induces a change in the birefringence of the crystal. As a result the polarization of the optical pulse incident on the EOS crystal changes in the locations where a terahertz field is present. The beam is then reflected from a front surface of the EOS crystal; this can occur due to a change in refractive index at the surface of the EOS crystal or due to dielectric coatings in a wavelength region that allow for reflection of the probe pulse (the technology of dielectric coatings is well known to and understood by those skilled in the art). Having reflected off the front surface of the EOS crystal the probe pulse passes back through the beam splitter so that unwanted polarizations are rejected. The remaining light is detected on a photosensitive device. This device may be a CMOS or CCD chip although it is expected that a CMOS chip would provide a higher performance, lower cost system. Any kind of photosensitive array would be suitable. The areas where the visible or near infrared pulse is detected will correspond to the areas where terahertz is received on the EOS crystal.
According to a first aspect of the present invention there is therefore provided a terahertz radiation sensor, the sensor comprising an optical beam input to receive an optical probe beam a detector to modulate said probe beam responsive to terahertz radiation, and a photosensitive detector to provide an output responsive to said probe beam modulation, the sensor being configured to provide a first optical path between said optical beam input and said electro-optic detector and to provide a second optical path between said electro-optic detector and said photosensitive detector; and wherein said sensor further comprises a polarizer, said polarizer being located in both said first and said second optical paths.
If the probe pulse is already highly polarized there is no need for a polarizer in the first optical path.
In a related aspect the invention provides a terahertz radiation sensor, the sensor comprising an optical beam input to receive an optical probe beam, a detector to modulate said probe beam responsive to terahertz radiation, and a photosensitive detector to provide an output responsive to said probe beam modulation, the sensor being configured to provide a first optical path between said optical beam input and said electro-optic detector and to provide a second optical path between said electro-optic detector and said photosensitive detector; and wherein said sensor further comprises a beam splitter located in both said first and said second optical paths.
Preferably the first and second optical paths traverse the detector, which may comprise an electro-optic crystal, and are preferably substantially counter propagating within the detector. Preferably the polarizer/beam splitter comprises a polarizing beam splitter. Using a single polarizing beam splitter to both direct light onto and filter light from a detector such as an electro-optic crystal facilitates the construction of a probe head with a high extinction ratio, for example of the order of 1:106. In embodiments such a configuration also facilitates fabrication of a compact and rugged probe head, as well as potentially reducing the bill of materials for the device. In embodiments the configuration also allows objects to be placed directly in front of the terahertz detector, allowing for near field imaging (Section 10.1.2 on page 437 of Optics (third edition) by Eugene Hecht describes near field (or Fresnel) diffraction and far field (or Fraunhofer) diffraction; the regimes are the same for imaging), for example where back-reflection of the probe beam by the detector is employed. Furthermore situating a polarizing beam splitter, for example comprising two trapezoid prisms, such that it both directs a probe pulse onto an electro-optic sampling crystal and filters the output from the EOS crystal makes the sensor relatively insensitive to vibration, and in embodiments, in effect makes the sensor system self-aligning.
Preferably the sensor also includes a birefringence corrector, preferably positioned so that it lies in both the first and second optical paths. This can be used to compensate for intrinsic birefringence in the detector.
The sensor is preferably configured so that the probe beam is internally reflected off a front face of the detector. This face of the detector may be provided with an optical filter to filter out ambient light which might otherwise reach the photosensitive detector and/or a multi layer stack to enhance transmission of terahertz radiation into the detector.
Preferably the photosensitive detector comprises an imaging sensor such as a CCD (charge coupled device) or CMOS (complimentary metal oxide semiconductor) sensor. Employing a photosensitive array facilitates the detection of terahertz radiation at a plurality of pixels simultaneously (area detection). This has a number of advantages in embodiments, for example, allowing faster acquisition times, allowing for single shot measurement of terahertz pulses making an imaging system less susceptible to interpulse variations and laser instability, and avoiding the need for a mechanical scanning system, thus reducing cost and complexity. A CMOS image sensing device is particularly preferred as this enables the use of higher probe beam chop rates (as described further later) thus reducing noise and improving sensitivity, as well as reducing cost and providing a more compact probe head (relative to use of a CCD sensor). Moreover in embodiments the use of CMOS detection as compared with conventional balanced diode detection can avoid the need for expensive lock-in device.
In some embodiments of the terahertz radiation sensor the sensor (front face) lacks a terahertz focusing system such as a lens or mirror and is a non-terahertz-imaging sensor (that is, the incoming terahertz radiation is not imaged, for example onto the detector). Thus the sensor may be configured that terahertz radiation falling on the detector is substantially unfocussed. This can further reduce cost and increase ruggedness and compactness of the sensor, as well as facilitating near field imaging in potentially hostile environments such as medical/veterinary endoscopic imaging. Near field imaging also potentially provides higher image resolution than available by other techniques.
Thus in another aspect the invention provides a terahertz radiation sensor comprising a terahertz detector configured to receive terahertz radiation and to provide an optical signal, and an optical image sensor coupled to said terahertz detector to provide a data output responsive to said optical signal, and wherein said detector is not in a terahertz radiation image plane.
In this context it will be understood from the description later that it is not necessary for an optical image to be formed on the optical image sensor. Thus the optical image sensor need not lie in an optical image plane compensation for this being applied, for example, by post-processing applied to the sensor output.
The invention also provides a terahertz radiation sensor, the sensor comprising an optical beam input to receive an optical probe beam, a detector to modulate said probe beam responsive to terahertz radiation, and a photosensitive detector to provide an output responsive to said probe beam modulation, the sensor being configured to provide a first optical path between said optical beam input and said electro-optic detector and to provide a second optical path between said electro-optic detector and said photosensitive detector; and wherein said sensor further comprises one or more dielectric layers on a surface of said electro-optic detector configured to enhance transmission of said terahertz radiation into said detector, a said dielectric layer comprising a material transmissive at a terahertz wavelength and having a thickness substantially equal to an odd integral number of quarter wavelengths of said terahertz radiation in said material.
There are also dielectric stack configurations that do not involve quarter wavelength thicknesses, dependent upon the number of layers deposited and their refractive index.
In a further aspect the invention provides a terahertz radiation sensor, the sensor comprising an optical beam input to receive an optical probe beam, a detector to modulate said probe beam responsive to terahertz radiation, and a photosensitive detector to provide an output responsive to said probe beam modulation, the sensor being configured to provide a first optical path between said optical beam input and said electro-optic detector and to provide a second optical path between said electro-optic detector and said photosensitive detector, and wherein said sensor further comprises a polarizing beam splitter located in both said first and said second optical paths, wherein both said first and said second optical paths traverse said detector, and wherein between said polarizing beam splitter and said detector said first and second optical paths substantially coincide.
The invention also provides terahertz imaging systems, that is systems to image an object in 2 or 3 dimensions using terahertz radiation, without necessarily imaging detected terahertz radiation. The systems include the above sensors and, in particular, may further comprise a terahertz emitter, for example controllable by an optical pulse, and means to adjust a relative timing of emission by the terahertz emitter and detection by a sensor as described above in order to determine a 2-or 3-dimensional image of an object in reflection and/or transmission.
The invention further provides a method of sensing terahertz radiation, the method comprising providing an optical probe beam to an electro-optic detector, modulating a polarization of said optical beam at said electro-optic detector using said terahertz radiation, and sensing said terahertz radiation by detecting said polarization modulation, the method further comprising polarizing said optical beam and detecting said polarization modulation using a common polarizing beam-splitter.
In a yet further aspect the invention provides a terahertz imaging system comprising a terahertz source to provide terahertz radiation for imaging an object, a terahertz sensor to detect said terahertz radiation, a timing adjust device to adjust a relative timing of emission of terahertz radiation from said source and detection of said terahertz radiation by said detector, and an image processor to receive a plurality of sensed signals from said terahertz sensor for a plurality of different said relative timings, and to output data, representing a three-dimensional image of said object determined from said plurality of sensed signals.
Such a system images or forms a representation of an object using terahertz radiation but need not (and in preferred embodiments does not) image the terahertz radiation used for imaging.
Preferably the terahertz sensor provides a 2-dimensional image output for the image processor although in other embodiments an x-y scanning system may be employed. The imaging may be performed in transmission and/or reflection; preferably the terahertz source and sensor comprise optically-controlled devices. The source may comprise a photoconductive emitter (see, for example, EP 0828143A, WO 01/38929, and the co-pending UK patent application no. GB 0307096.8 filed on 27, Mar. 2003). The sensor may comprise an electro-optic crystal, which may conveniently be driven by, for example, a common pulsed laser system with a pair of output beams and an optical timing adjust mechanism.
Broadly speaking the inventor has recognized that using successive sampling techniques, sometimes known as equivalent time sample in the context of fast digital oscilloscopes, the terahertz radiation transmitted through and/or reflected from an imaged object may not only be sensed but also characterized to, in effect, determine a wave shape of a transmitted and/or reflected E-field component of the radiation. With this information, and knowing the E-field waveshape of the terahertz radiation illuminating the imaged object, a 3-dimensional configuration of the object may be determined, by determining what 3-dimensional configuration, when applied to the illuminating radiation, gives the observed result. To perform this calculation it is not necessary that the terahertz radiation be imaged onto the radiation sensor—in effect the lack of such imaging may be taken into account by the calculation. This can be understood by thinking of the image processing as replacing the physical imaging of the terahertz radiation although, in practice, the image processor will not generally explicitly perform such a calculation.
Preferably, therefore, the terahertz sensor comprises a non-terahertz radiation imaging sensor—the terahertz detection plane (for example, the front surface of the sensor) need not be at or close to a terahertz imaging plane. The inventor has recognized that it is nonetheless possible to image an object without a physical terahertz imaging system. For example, as alluded to above, terahertz detection may be performed by phase-sensitive electro-optic sampling. As previously described in such an arrangement the probe radiation may be provided to the back of a detector crystal such an EOS crystal in order to create a phase sensitive imaging plane that can easily be sampled at multiple pixels, for example by means of a CMOS optical imaging array.
Preferably the terahertz source comprises an extended source rather than a point source. This facilitates perfusion of an object with terahertz radiation that will reach the detection plane. For some objects, such as a lens like object, perfusion of an object with terahertz radiation may be best achieved with a source that is point like, although in general a diffuse source is normally preferred.
The imaging system can be thought of as generating an interference pattern or hologram of the object (although preferably a substantially complete E-field waveform is available rather than merely amplitude and phase information) and the image processor can be thought of as reconstructing a 3-dimensional image of the object by decoding the interference pattern or hologram. Because the THz radiation tends to diffract strongly from objects with similar length scales to the terahertz radiation wavelength this method is particularly advantageous for imaging systems which contain detail on the same length scale as the terahertz wavelength.
This can be done by Fourier transforming, inverting, and then inverse Fourier transforming captured image data (although in other embodiments other transforms and inverse transforms may be employed). The captured data comprises effectively 3-dimensional data, having two spatial dimensions (x and y) and a third, time dimension (which effectively corresponds to a spatial dimension. Thus the transforms are preferably performed in three dimensions. In order to compensate for effects of the transmitted pulse shape a set of data is captured without an object present; this may then be deconvolved (in three dimensions) from data captured when an object is present. Optionally a reconstructed 3-dimensional image may be filtered or otherwise post-processed as desired.
It will be recognized that in embodiments of the above described technique a voxel of a 3D image is dependent, in general, on a complete set of captured image data (in x, y and t) rather than, say, just on a single time slice through the object. Thus preferably data representing a region of a calculated 3D image is determined from sensed signals for more than one and preferably substantially all of the plurality of different relative timings (unlike prior art techniques such as those described in GB 2350673A where, in effect, a method analogous to sonar is employed).
In a related aspect the invention provides a method of processing data from a terahertz imaging system. The method comprising inputting a plurality of sets of two-dimensional data each representing a pattern of terahertz radiation at one of a succession of sampling times of a terahertz imaging waveshape illuminating an imaged object, and generating image data representing a substantially three-dimensional image of said imaged object from said input data, wherein said generating comprises determining data for a voxel of said three-dimensional data from a plurality of said sets of two-dimensional data.
By shifting the sampling point by a few picoseconds or a few femtoseconds each time a set of 2-dimensional data is acquired a picture of the imaging waveshape as modified by the imaged object can gradually be built up and this data can then be transformed to determine a set of 3-dimensional image voxels. Optionally a priori information relating to the refractive index or refractive index range of the imaged object may be taken into account, and/or refractive index variations may be taken into account by the procedure, including absorption, that is the imaginary part of a complex refractive index.
In another aspect the invention provides a method of determining three-dimensional image data for an object using terahertz radiation. The method comprising illuminating the object with terahertz radiation having a repeated waveform, sampling, over an area, terahertz radiation transmitted and/or reflected by the object at a succession of time offsets to build up an image of said waveform after transmission through and/or reflection from said object, and determining image data for a three-dimensional image of the object which, when applied to said illuminating waveform, results in said waveform after transmission and/or reflection.
Preferably this method further comprises determining the illuminating waveform by a similar illuminating and sampling procedure without the object present.
The invention further provides systems for implementing the above described methods. Thus, in particular, the invention provides a system for determining three-dimensional image data for an object using terahertz radiation, the system comprising means for illuminating the object with terahertz radiation having a repeated waveform. Means for sampling, over an area, terahertz radiation transmitted and/or reflected by the object at a succession of time offsets to build up an image of said waveform after transmission through and/or reflection from said object, and means for determining image data for a three-dimensional image of the object which, when applied to said illuminating waveform, results in said waveform after transmission and/or reflection.
The invention also provides processor control code to implement the above described methods, in particular on a data carrier such as a disk, CD or DVD-ROM, programmed memory such as read only memory (firmware), or on a data carrier such as an optical signal carrier. Embodiments of the above described methods may be implemented using a general purpose computer system or a digital signal processor (DSP), or an application specific integrated circuit (ASIC) or field programmable gate array (FPGA). Thus code to implement embodiments of aspects of the invention may comprise code in a conventional programming language such as C, or microcode, or code for setting up or controlling an ASIC or FPGA, or code for a hardware description language such as Verilog (trademark). As the skilled person will appreciate relevant image processing code and/or data may be distributed between a plurality of coupled components in communication with one another, for example across a network.
Thus in particular the invention also provides computer program code for a terahertz imaging system comprising a terahertz source to provide terahertz radiation for imaging an object, a terahertz sensor to detect said terahertz radiation, and a timing adjust device to adjust a relative timing of emission of terahertz radiation from said source and detection of said terahertz radiation by said detector, said terahertz sensor providing a two-dimensional data output, the code being configured to, when running: capture image data from a plurality of sets of two-dimensional images of said terahertz radiation for a plurality of said relative timings, said captured image data representing an interference pattern generated by said object, and construct a three-dimensional image of said object from said captured image data by decoding said interference pattern.
In still further aspects of the invention, the above described sampling at a succession of times and means to perform this/adjust a relative timing may be replaced by sampling at a succession of frequencies and means to perform this/adjust a terahertz radiation frequency. In such a configuration, however (for example, where relatively narrowband terahertz radiation is generated, say by photomixing) a sampling time or time delay may still be altered to find the maximum (or maxima) of the detector's output.
The invention also provides for a means to image objects placed substantially adjacent to a plane that is sensitive to terahertz. For electro-optic sampling to be effective it is essential that the terahertz radiation and the probe radiation be co-propagating in the electro-optic sensor for this reason conventional teaching was that systems that measure terahertz by electro-optic sampling or antenna detection required the collinear insertion of probe radiation in to the terahertz sensor. The requirement that the probe beam be inserted collinearly with the terahertz pulse prevents objects to be imaged, said objects being substantially transparent to terahertz but opaque to the probe radiation, being placed in the probe beam path, and thus adjacent to the detector, see for example Lu Z G, Campbell P, Zhang X C, Free-space electro-optic sampling with a high repetition-rate regenerative amplified laser, Appl. Phys. Lett 71(5) 4, Aug. 1997. The inventor has recognized that it is in fact possible to insert the probe radiation in to the rear of the detector with the terahertz radiation configured to arrive from the front. In this case the detector is configured such that the probe radiation reflects off the front surface of the detector and from that moment on is co-propagating with the terahertz radiation. Thus the inventor has further recognized that it is possible to place objects adjacent to the detector surface. The inventor has further recognized that the ability to place objects close to the surface of the detector allows a substantial solid angle of the radiation that is reflected, scattered, or transmitted through the object to be measured. In other word the inventor has recognized that such a sensor would be able to sense terahertz in the fresnel region. The inventor has further recognized that capturing a substantial solid angle of the radiation from the object is an important pre-requisite to being able to effectively conduct imaging or inverse problem resolution using the holographic techniques that the inventor has recognized could be applied to imaging with terahertz. The inventor has further recognized that given such a sensor it may be possible to engage in other methods to reconstruct the nature of an arbitrary object placed close to the sensor that does not make strict use of the holography based techniques described herein.
These and other aspects of the invention will now be further described, by way of example only, with reference to the accompanying figures in which:
a to 1c show a simplified embodiment of terahertz radiation sensor according to an aspect of the present invention illustrating stages in the passage of a probe beam through the sensor;
a to 10c show examples of terahertz imaging systems using a short pulse laser system, and difference frequency mixing, respectively and successive sampling of terahertz pulses;
a and 12b show, respectively, a general purpose computer system programmed for implementing a terahertz imaging system control and image processing procedure embodying an aspect of the present invention, and a terahertz imaging system incorporating the computer system of
Referring to
In
A polarizing beam splitter 214 receives incoming probe beam light and directs this towards detector 206 and also analyses the polarization of light reflected back from detector 206, providing this to a light detection system comprising, in this embodiment, a lens 216, an aperture 218, and a photosensitive array 220. Polarizing beam splitters are well known to those skilled in the art and may comprise, for example, a pair of calcite crystals 214a, b. Lens 216 focuses an optical image from detector 206 onto a photosensitive array 220 and aperture 218 is preferably included to reduce the detection of stray light. The photosensitive array 220 preferably comprises a CMOS device although CCD and other photosensitive arrays may alternatively be employed. Preferably photosensitive device 220 provides 2-dimensional image data but in other embodiments a single photosensitive diode or detector may be employed to provide a single reading of terahertz radiation detected at detector 206. In such an arrangement a 2-dimensional image may be captured by scanning the detector.
Radiation from the probe laser beam may conveniently be coupled into the terahertz probe head 200 by means of an optical fiber 222, optionally employing one or more lenses 224 to collimate and/or expand the output of the fiber. Preferably fiber 222 comprises a polarization-preserving single mode fiber although other single mode fibers may be employed as well as photonic band gap fibers (zero dispersion photonic fibers are sold by Crystal Fibre A/S—Denmark) and potentially multi-mode fibers for 2D-imaging. A mirror 226 serves to direct the collimated probe light into the polarizing beam splitter 214. Preferably mirror 226 is aligned slightly off-axis so that light that is rejected back from the probe head is not substantially coupled back into the optical fiber 222 delivering the probe radiation.
Broadly speaking the lenses in the arrangements of
Referring now to
The terahertz generator 1010 may comprise any terahertz generator that can be switched on a time scale which is similar to that of the pump arm pulse. Example terahertz generators include Austin switches, surface field emission devices, electro-optic crystals, difference frequency mixing devices, quantum cascade lasers, high frequency klystrons, and high frequency and multiplied gunn diodes and other terahertz emitters, optionally provided with a terahertz amplifier. One particularly advantageous form of emitter is described in the co-pending UK patent application no. GB 0307096.8 filed on 27, Mar. 2003; this describes an emitter including a beam expander to allow an object to be illuminated with terahertz radiation from what is effectively an extended rather than a point source. This is advantageous in the imaging techniques which are described below. Other techniques which can be employed to illuminate a sample such that it is perfused with terahertz radiation include the use of conventional beam expander arrangements with terahertz lenses or reflectors and optimized for THz wavelengths, and cassegrain, Newtonian and other systems. Techniques for the construction of wide aperture emitters are well established and known to those skilled in the art.
b shows an alternative example of a terahertz imaging system 1050 (in which like elements are indicated by like reference numerals) in which the short pulse laser system 1002 is replaced by a laser system 1052 comprising a pair of lasers 1052a, b emitting laser beams with a difference in frequency that lies in the terahertz range, these beams interfering when they come together. The two laser beams may be provided by two different laser systems, as shown in
In order to synchronize terahertz generation and terahertz detection a variable optical delay 1014 can be placed in either or both of the probe arm 1008 and pump 1006. This provides a basic imaging system.
An ultrafast pulse generally stretches when passed down an optical fiber; a small amount of such stretching does not significantly affect the detection/imaging described herein but to replace stretching a pulse can be prechirped; Coherent Inc of California USA sell systems to do this. Additionally or alternatively sapphire optical fiber may be employed or any fiber that preserves the pulse length of the pulse as it is transmitted down the fiber including photonic bandgap fibres.
We will next describe a more sophisticated technique in which such an optical delay is altered to sample successive terahertz pulses at steadily increasing offsets to build up an image or template of the shape of the electric field amplitude of a terahertz pulse received by the probe head.
Having captured a template of the received terahertz pulse a 3-dimensional image of the object 1012 (including terahertz refractive index variations) can be determined by solving an inverse problem to determine a representation of the object which, when applied to the transmitted terahertz radiation pulse, gives rise to the received signal. There are potentially many ways of solving this problem; the mathematics of one suitable technique are described in “Optical scanning holography as a technique for high-resolution three-dimensional biological microscopy”; Jim Swoger, Manuel Martínez-Corral, Jan Huisken, & Ernst H. K. Stelzer, J. Opt. Soc. Am. A. Vol. 19, No. 9, September 2002, 1910-1918, hereby incorporated by reference. An outline description of this procedure is given later. It will be appreciated that the inverse problem may be solved to determine a (2 or) 3-dimensional configuration of an object in either transmission or reflection or both.
It will be appreciated that with this technique it is important to make an accurate determination of the transmitted terahertz pulse shape and this can be done by making a similar set of measurements, varying the optical delay, without an object present in the arrangements of
In order to enhance the signal-to-noise ratio the pump arm 1006 or terahertz generator 1010 of the system of
To gather a complete set of data a number of output frames is collected at a plurality of different optical delays. Optionally (e.g. in a photomixing system) data may also be collected at a plurality of different terahertz frequencies, which facilitates compositional analysis of an object. This can be achieved, for example, using the imaging system of
The above described probe heads are able to provide a data set which measures terahertz intensity in two dimensions (x and y) and scanning through a set of different time delays provides a 3-dimensional data set including the time dimension. Fourier transforms can be used to determine a 3-dimensional image of the interior of an object from this 3-dimensional data set. Employing the terminology of Swoger et al. (ibid) a hologram generation process is describable by a correlation of a sample distribution with an encoding function:
H(r)=ƒsample(r){circle around (×)}ƒenc(r),
where H(r) is the hologram, ƒsample is the characteristic of the sample that we want to image, ƒenc is an encoding function, r is a lateral position vector (x, y, t), and {circle around (×)} donates correlation of two scalars.
A reconstructed image of ƒsample can be generated by
Irec(r)=H(r){circle around (×)}Rƒ(r),
where Rƒ is a reconstruction function, for an ideal, one-to-one reconstruction of the sample given by
Rƒ(r)−FT−1[1/FT(ƒenc(r)))],
where FT and FT−1 denote Fourier and inverse Fourier transforms, respectively.
In practice the Fourier transforms may be applied numerically to a 3-dimensional data set captured using a probe head by means of a fast Fourier transform (FFT) to reconstruct a 3-dimensional (ideal) image of an object using these “holographic” techniques. The skilled person will appreciate that, in a conventional manner, such discrete FFTs should be implemented in such a way as to avoid division by zero. The skilled person will further appreciate that once an image of an object has been constructed the image may be filtered in many ways, for example to add a bandwidth limitation or to provide edge detection, for example by means of a Gabor filter.
In the above H comprises the effect of the sample convoluted with the encoding function. Thus if we take an image without a sample present this will give us a value of H for blank space, a suitable null value. When trying to reconstruct an image of the sample the acquisition of an image has two phases, one is taking a reference image (with no sample present), one is taking a sample image (with the sample present). In
Irec(r)=H(r){circle around (×)}Rƒ(r),
H(r) is simply the reference image, (the same reference image may be used for multiple samples). Having thus acquired Irec(r) this may be viewed with 3 dimensional visualization software such as OpenDX.
Referring next to
At Step S1100 a user inputs a set of parameters for a scan into the procedure and at Step S1102 a time delay t (and optionally a terahertz frequency f) are set and x-y image data are input from the terahertz probe, optionally averaged over a plurality of laser pulses (S1104). The captured 2-dimensional data is then stored, for example in an array (S1106) and then the time delay t (and optionally frequency f, depending on the scan parameters) are adjusted (Step S1108) and the procedure loops back to Step S1102 until the scan in time (and optionally frequency) is completed. Steps S1100a to S1108a may be employed to calibrate the imaging apparatus to determine the shape of a transmitted terahertz radiation pulse; Steps S1102a to S1108a correspond to Steps S1102 to S1108.
The image processing procedure begins at Step S1110 and is described in outline only; for further details reference may be made, for example, to the Swoger et al. paper. Thus at Step S1110 the procedure reads a complete set of {x, y, t} data which determines the encoding function ƒenc. Then, at Step S1112, the procedure calculates the sample reconstruction function Rƒ and from this a 3-dimensional reconstruction of the imaged sample or object; the procedure then ends at Step S1114.
b shows a terahertz imaging system 1250 incorporating the general purpose computer system 1200 of
No doubt many effective alternatives will occur to the skilled person. For example although the described embodiments of the terahertz probe lack a terahertz radiation imaging system, such a system, for example a terahertz lens, may be included in the terahertz probe prior to detector 206. A coherent detection system (as described above) may also be used, for example, for spectroscopy, analogously to FTIR (Fourier Transform Infrared) spectroscopy.
It will be understood that the invention is not limited to the described embodiments and encompasses modifications apparent to those skilled in the art lying within the spirit and scope of the claims appended hereto.
| Number | Date | Country | Kind |
|---|---|---|---|
| 0401276.1 | Jan 2004 | GB | national |
| Number | Name | Date | Kind |
|---|---|---|---|
| 5894125 | Brener et al. | Apr 1999 | A |
| 6239866 | Bromage et al. | May 2001 | B1 |
| 6690001 | Jiang et al. | Feb 2004 | B2 |
| 6828558 | Arnone et al. | Dec 2004 | B1 |
| 20010038074 | Zhang et al. | Nov 2001 | A1 |
| 20010038471 | Zhang et al. | Nov 2001 | A1 |
| 20020153874 | Jiang et al. | Oct 2002 | A1 |
| 20030155512 | Arnone et al. | Aug 2003 | A1 |
| 20030165003 | Ciesla et al. | Sep 2003 | A1 |
| 20030178584 | Arnone et al. | Sep 2003 | A1 |
| Number | Date | Country |
|---|---|---|
| 0841 548 | Oct 1997 | EP |
| 0864857 | Sep 1998 | EP |
| 1271 115 | Jan 2003 | EP |
| 0727671 | Sep 2003 | EP |
| 2350 673 | Dec 2000 | GB |
| 2350673 | Dec 2000 | GB |
| 2352512 | Jan 2001 | GB |
| 2410081 | Feb 2007 | GB |
| 2425833 | Feb 2007 | GB |
| WO 0177646 | Oct 2001 | WO |
| WO03042070 | May 2003 | WO |
| WO 03095991 | Nov 2003 | WO |
| Number | Date | Country | |
|---|---|---|---|
| 20050156110 A1 | Jul 2005 | US |
