Terahertz tomographic imaging of XVIIIth Dynasty Egyptian sealed pottery

600-GHz Fourier imaging based on heterodyne detection at the 2nd sub-harmonic

Fourier imaging is an indirect imaging method which records the diffraction pattern of the object scene coherently in the focal plane of the imaging system and reconstructs the image using computational resources. The spatial resolution, which can be reached, depends on one hand on the wavelength of the radiation, but also on the capability to measure - in the focal plane - Fourier components with high spatial wave-vectors. This leads to a conflicting situation at THz frequencies, because choosing a shorter wavelength for better resolution usually comes at the cost of less radiation power, concomitant with a loss of dynamic range, which limits the detection of higher Fourier components. Here, aiming at maintaining a high dynamic range and limiting the system costs, we adopt heterodyne detection at the 2nd sub-harmonic, working with continuous-wave (CW) radiation for object illumination at 600 GHz and local-oscillator (LO) radiation at 300 GHz. The detector is a single-pixel broad-band Si CMOS TeraFET equipped with substrate lenses on both the front- and backside for separate in-coupling of the waves. The entire scene is illuminated by the object wave, and the Fourier spectrum is recorded by raster scanning of the single-detector unit through the focal plane. With only 56 µW of power of the 600-GHz radiation, a dynamic range of 60 dB is reached, sufficient to detect the entire accessible Fourier space spectrum in the test measurements. We present a detailed comparison between plane-to-plane imaging and Fourier imaging, and show that, with both, a lateral spatial resolution of better than 0.5 mm, at the diffraction limit, is reached.

Efficient iterative reconstruction with beam shape compensation for THz computed tomography

Terahertz (THz) computed tomography is an emerging nondestructive and non-ionizing imaging method. Most THz reconstruction methods rely on the Radon transform, originating from x-ray imaging, in which x rays propagate in straight lines. However, a THz beam has a finite width, and ignoring its shape results in blurred reconstructed images. Moreover, accounting for the THz beam model in a straightforward way in an iterative reconstruction method results in extreme demands in memory and in slow convergence. In this paper, we propose an efficient iterative reconstruction that incorporates the THz beam shape, while avoiding the above disadvantages. Both simulation and real experiments show that our approach results in improved resolution recovery in the reconstructed image. Furthermore, we propose a suitable preconditioner to improve the convergence speed of our reconstruction.

3D image reconstruction of terahertz computed tomography at sparse angles by total variation minimization

Continuous-wave terahertz computed tomography (THz-CT) is an important three-dimensional imaging method for probing the profile and inner properties of a sample's structure. We applied the total variation (TV) minimization iterative algorithm to squeeze 75% data acquisition time of THz-CT without the loss of reconstruction fidelity. The imaging system is built based on a 278.6 GHz avalanche diode source. A zero-order Bessel beam is generated by an axicon, for which the intensity profile remains essentially propagation invariant within the non-diffracting zone. The effectiveness of the proposed method is verified by using three optically opaque objects. The reconstruction results show that the TV-minimization algorithm can effectively suppress noise, artefacts, and shape distortion created in sparse angle reconstruction.

Layer-Resolving Terahertz Light-Field Imaging Based on Angular Intensity Filtering Method

Terahertz focal plane array imaging methods, direct camera imaging and conventional light field imaging methods are incapable of resolving and separating layers of multilayer objects. In this paper, for the purpose of fast, high-resolution and layer-resolving imaging of multilayer structures with different reflection characteristics, a novel angular intensity filtering (AIF) method based on terahertz light-field imaging is purposed. The method utilizes the extra dimensional information from the 4D light field and the reflection characteristics of the imaging object, and the method is capable to resolve and reconstruct layers individually. The feasibility of the method is validated by experiment on both "idealized" and "practical" multilayer samples, and the advantages in performance of the method are proven by quantitative comparison with conventional methods.

Twenty years of terahertz imaging [Invited]

The birth of terahertz imaging approximately coincides with the birth of the journal Optics Express. The 20^th anniversary of the journal is therefore an opportune moment to consider the state of progress in the field of terahertz imaging. This article discusses some of the compelling reasons that one may wish to form images in the THz range, in order to provide a perspective of how far the field has come since the early demonstrations of the mid-1990's. It then focuses on a few of the more prominent frontiers of current research, highlighting their impacts on both fundamental science and applications.

Ordered subsets convex algorithm for 3D terahertz transmission tomography

We investigate in this paper a new reconstruction method in order to perform 3D Terahertz (THz) tomography using a continuous wave acquisition setup in transmission mode. This method is based on the Maximum Likelihood for TRansmission tomography (ML-TR) first developed for X-ray imaging. We optimize the Ordered Subsets Convex (OSC) implementation of the ML-TR by including the Gaussian propagation model of THz waves and take into account the intensity distributions of both blank calibration scan and dark-field measured on THz detectors. THz ML-TR reconstruction quality and accuracy are discussed and compared to other tomographic reconstructions.

Methodology for materials analysis using swept-frequency feedback interferometry with terahertz frequency quantum cascade lasers

Recently, we demonstrated an interferometric materials analysis scheme at terahertz frequencies based on the self-mixing effect in terahertz quantum cascade lasers. Here, we examine the impact of variations in laser operating parameters, target characteristics, laser-target system properties, and the quality calibration standards on our scheme. We show that our coherent scheme is intrinsically most sensitive to fluctuations in interferometric phase, arising primarily from variations in external cavity length. Moreover we demonstrate that the smallest experimental uncertainties in the determination of extinction coefficients are expected for lossy materials.

Fast three-dimensional terahertz computed tomography using real-time line projection of intense terahertz pulse

We demonstrated fast three-dimensional transmission terahertz computed tomography by using real-time line projection of intense terahertz beam generated by optical rectification in lithium niobate crystal. After emphasizing the advantage of intense terahertz pulse generation for two-dimensional spatio-temporal terahertz imaging, peak-to-peak amplitudes of pulsed terahertz electric field have been used to obtain a series of projection images at different rotation angles. Then a standard reconstruction algorithm has been employed to perform final three-dimensional reconstruction. Test samples including a medicine capsule have been investigated with a total acquisition time to only 6 minutes.

Three-dimensional terahertz computed tomography of human bones

Three-dimensional terahertz computed tomography has been used to investigate dried human bones such as a lumbar vertebra, a coxal bone, and a skull, with a direct comparison with standard radiography. In spite of lower spatial resolution compared with x-ray, terahertz imaging clearly discerns a compact bone from a spongy one, with strong terahertz absorption as shown by additional terahertz time-domain transmission spectroscopy.