Terahertz sparse deconvolution based on an iterative shrinkage algorithm is presented in this study to characterize multilayered structures. With an upsampling approach, sparse deconvolution with superresolution is developed to overcome the time resolution limited by the sampling period in the measurement and increase the precision of the estimation of echo arrival times. A simple but effective time-domain model for describing the temporal pulse spreading due to the frequency-dependent loss is also designed and introduced into the algorithm, which greatly improves the performance of sparse deconvolution in processing time-varying pulses during the propagation of terahertz waves in materials. Numerical simulations and experimental measurements verify the algorithms and show that sparse deconvolution can be considered as an effective tool for terahertz nondestructive characterization of multilayered structures.

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.

The process by which art paintings are produced typically involves the successive applications of preparatory and paint layers to a canvas or other support; however, there is an absence of nondestructive modalities to provide a global mapping of the stratigraphy, information that is crucial for evaluation of its authenticity and attribution, for insights into historical or artist-specific techniques, as well as for conservation. We demonstrate sparsity-based terahertz reflectometry can be applied to extract a detailed 3D mapping of the layer structure of the 17th century easel painting Madonna in Preghiera by the workshop of Giovanni Battista Salvi da Sassoferrato, in which the structure of the canvas support, the ground, imprimatura, underpainting, pictorial, and varnish layers are identified quantitatively. In addition, a hitherto unidentified restoration of the varnish has been found. Our approach unlocks the full promise of terahertz reflectometry to provide a global and detailed account of an easel painting’s stratigraphy by exploiting the sparse deconvolution, without which terahertz reflectometry in the past has only provided a meager tool for the characterization of paintings with paint-layer thicknesses smaller than 50 μm. The proposed modality can also be employed across a broad range of applications in nondestructive testing and biomedical imaging.