Terahertz frequency-wavelet domain deconvolution for stratigraphic and subsurface investigation of art painting.

Terahertz frequency-wavelet deconvolution is utilized specifically for the stratigraphic and subsurface investigation of art paintings with terahertz reflective imaging. In order to resolve the optically thin paint layers, a deconvolution technique is enhanced by the combination of frequency-domain filtering and stationary wavelet shrinkage, and applied to investigate a mid-20th century Italian oil painting on paperboard, After Fishing, by Ausonio Tanda. Based on the deconvolved terahertz data, the stratigraphy of the painting including the paint layers is reconstructed and subsurface features are clearly revealed, demonstrating that terahertz frequency-wavelet deconvolution can be an effective tool to characterize stratified systems with optically thin layers.

Sub-surface terahertz imaging through uneven surfaces: visualizing Neolithic wall paintings in Çatalhöyük.

Pulsed terahertz imaging is being developed as a technique to image obscured mural paintings. Due to significant advances in terahertz technology, portable systems are now capable of operating in unregulated environments and this has prompted their use on archaeological excavations. August 2011 saw the first use of pulsed terahertz imaging at the archaeological site of Çatalhöyük, Turkey, where mural paintings dating from the Neolithic period are continuously being uncovered by archaeologists. In these particular paintings the paint is applied onto an uneven surface, and then covered by an equally uneven surface. Traditional terahertz data analysis has proven unsuccessful at sub-surface imaging of these paintings due to the effect of these uneven surfaces. For the first time, an image processing technique is presented, based around Gaussian beam-mode coupling, which enables the visualization of the obscured painting.

Terahertz deconvolution.

The ability to retrieve information from different layers within a stratified sample using terahertz pulsed reflection imaging and spectroscopy has traditionally been resolution limited by the pulse width available. In this paper, a deconvolution algorithm is presented which circumvents this resolution limit, enabling deep sub-wavelength and sub-pulse width depth resolution. The algorithm is explained through theoretical investigation, and demonstrated by reconstructing signals reflected from boundaries in stratified materials that cannot be resolved directly from the unprocessed time-domain reflection signal. Furthermore, the deconvolution technique has been used to recreate sub-surface images from a stratified sample: imaging the reverse side of a piece of paper.

Optimization of apodization functions in terahertz transient spectrometry.

A quadratic programming optimization procedure for designing asymmetric apodization windows tailored to the shape of time-domain sample waveforms recorded using a terahertz transient spectrometer is proposed. By artificially degrading the waveforms, the performance of the designed window in both the time and the frequency domains is compared with that of conventional rectangular, triangular (Mertz), and Hamming windows. Examples of window optimization assuming Gaussian functions as the building elements of the apodization window are provided. The formulation is sufficiently general to accommodate other basis functions.

Materials for phantoms for terahertz pulsed imaging.

Phantoms are commonly used in medical imaging for quality assurance, calibration, research and teaching. They may include test patterns or simulations of organs, but in either case a tissue substitute medium is an important component of the phantom. The aim of this work was to identify materials suitable for use as tissue substitutes for the relatively new medical imaging modality terahertz pulsed imaging. Samples of different concentrations of the candidate materials TX151 and napthol green dye were prepared, and measurements made of the frequency-dependent absorption coefficient (0.5 to 1.5 THz) and refractive index (0.5 to 1.0 THz). These results were compared qualitatively with measurements made in a similar way on samples of excised human tissue (skin, adipose tissue and striated muscle). Both materials would be suitable for phantoms where the dominant mechanism to be simulated is absorption (approximately 100 cm(-1) at 1 THz) and where simulation of the strength of reflections from boundaries is not important; for example, test patterns for spatial resolution measurements. Only TX151 had a frequency-dependent refractive index close to that of tissue, and could therefore be used to simulate the layered structure of skin, the complexity of microvasculature or to investigate frequency-dependent interference effects that have been noted in terahertz images.

Modelling the propagation of terahertz radiation through a tissue simulating phantom.

Terahertz (THz) frequency radiation, 0.1 THz to 20 THz, is being investigated for biomedical imaging applications following the introduction of pulsed THz sources that produce picosecond pulses and function at room temperature. Owing to the broadband nature of the radiation, spectral and temporal information is available from radiation that has interacted with a sample; this information is exploited in the development of biomedical imaging tools and sensors. In this work, models to aid interpretation of broadband THz spectra were developed and evaluated. THz radiation lies on the boundary between regions best considered using a deterministic electromagnetic approach and those better analysed using a stochastic approach incorporating quantum mechanical effects, so two computational models to simulate the propagation of THz radiation in an absorbing medium were compared. The first was a thin film analysis and the second a stochastic Monte Carlo model. The Cole-Cole model was used to predict the variation with frequency of the physical properties of the sample and scattering was neglected. The two models were compared with measurements from a highly absorbing water-based phantom. The Monte Carlo model gave a prediction closer to experiment over 0.1 to 3 THz. Knowledge of the frequency-dependent physical properties, including the scattering characteristics, of the absorbing media is necessary. The thin film model is computationally simple to implement but is restricted by the geometry of the sample it can describe. The Monte Carlo framework, despite being initially more complex, provides greater flexibility to investigate more complicated sample geometries.