Abstract
This article deals with pulsed excitation terahertz computed tomography (THz CT). Opposite to x-ray CT, where just a single value (pixel) is obtained, in case of pulsed THz CT the time signal is acquired for each position. Recorded waveform can be parametrized - many features carrying various information about examined structure can be calculated. Based on this, multiparametric reconstruction algorithm was proposed: inverse Radon transform based reconstruction is applied for each parameter and then fusion of results is utilized. Performance of the proposed imaging scheme was experimentally verified using dielectric phantoms.
1 Introduction
Electromagnetic waves in terahertz frequency range (T-rays) enable non-invasive, non-contact and non-ionizing inspection of dielectric materials and semiconductors. T-rays are sensitive to complex permittivity changes, therefore any defect or structure detail causing its noticeable disturbance can be detected. The most common method of defects localization is transmission or reflection imaging based on pulsed terahertz Time Domain Spectroscopy (THz TDS) [1,2,3,4,5].
Tomography is defined as the cross-sectional imaging of an object under test (OUT) by measuring scattered wave (rays). There are several terahertz imaging techniques [6,7,8,9,10]:
time of flight reflection tomography (THz ToFRT),
digital holography,
diffraction tomography (THz DT),
tomography with binary lens,
computed tomography (THz CT).
Terahertz CT enables analogous to conventional x-ray CT reconstruction of the objects physical properties distribution based on sectional slices (Figure 1). However, in case of terahertz tomography both amplitude and phase data are collected (opposite to intensity of radiation for x-ray CT). The above fact causes that THz CT provides more information about the examined object [1]. For each slice object under test is rotated by angle Θ and shifted by distance r. The terahertz radiation passing through the object is recorded in the measuring plane, where the projection R (Θ, r) is obtained. In Ref. [7] it has been shown that the change in phase of the electromagnetic wave in the measuring plane can be used to determine the distribution of the refractive index n (and consequently permittivity ∊) in the cross-section of the object. Forward problem (projection R (Θ, r) calculation) mathematically is described by Radon transform [7, 9, 11]:
where: g(x,y) – spatial distribution of the selected electromagnetic parameter of the imaged object, e.g. refractive index n(x, y), δ (•) – Dirac delta function, Θ – projection angle, r – ray distance from the axis of rotation.
The inverse problem – the reconstruction of the spatial distribution of the OUT’s electromagnetic parameters – is mathematically described by the inverse Radon transform [9, 11]:
In this paper terahertz computed tomography system with pulsed excitation is shown. Moreover, the reconstruction algorithm based on standard inverse Radon transform (2) and artificial neural network is proposed and verified using tomographic inspections of phantoms made of foamed polystyrene.
2 Measuring system
The proposed tomographic system is based on TRay4000 pulsed terahertz spectroscope of Picometrix (Figure 2). Picosecond order pulses are generated and detected by transmitter and receiver heads using photoconductive antennas (PCA). Terahertz beam focused by HDPE lenses is transmitted through examined object. The OUT was mounted on a (Θ, r, z) rotation/translation stage and scanned in order to obtain a slice image. For each position r and projection angle Θ, a transmission measurement is performed and a time-domain signal s(t) is recorded. Exemplary waveform obtained during tomographic inspection as well as some basic features of THz pulse are shown in Figure 3. Whole system was controlled by a software implemented in Matlab environment.
3 Reconstruction algorithm
Other (than phase shift) parameters of the signal measured in the measurement plane may be used for reconstruction. In the case of pulsed excitation (TDS), the pulse energy, amplitude or time delay introduced by the presence of the object is determined. Time signal acquisition makes it possible to design a very large number of time domain (TD), frequency (FD) and joint time-frequency domain (TFD) parameters. Depending on the parameter used during reconstruction, the resulting image (slice) will have a different character, i.e. other features will be exposed. For example, the use of frequency parameters (the amplitudes of different harmonics) will increase the sensitivity of detection of materials with different frequency dependent refractive indices or absorption coefficients. The use of multiple parameters reconstruction and subsequent fusion of the resulting distributions enables to achieve better results, which will be shown later in this paper. Generally, the multiparametric approach is becoming more popular in non-destructive testing [12].
Based on acquired s(t) time-domain waveform, a set of features ζn was calculated both in time and frequency domain. Each obtained parameter ζn can be utilized in order to form a sinogram (Figure 1). Sinograms acquired by measurements as well as reconstruction results using inverse Radon transform [11] in case of selected parameters ζn of waveform s(t) are shown in Figure 4. Based on observable differences in obtained images, the following conclusion can be drawn – selected parameters carry various information content about spatial distribution of dielectric properties in object under test. Therefore, multiparametric reconstruction algorithm using Radon’s inverse transformation and artificial neural networks in the process of combining information derived from different parameters is proposed. The block diagram of the proposed procedure is shown in Figure 5. The first step is to perform an inspection using the developed tomographic system. Then a preprocessing in form of filtration and deconvolution is carried out. Filtration may be implemented in the form of a one-dimensional median filter applied directly to the measurement signal s(t) or a two-dimensional median filter applied on the sinogram RM(Θ, r). In both cases, instead of median filtration, one- or two-dimensional wavelet denoising may be used. Another task implemented at this stage is the deconvolution of the sinogram RM(Θ, r) – blur limiting operation due to utilization of terahertz Gaussian beam instead of the straight ray beam. The data processed in this way are used in the next step – multi-parametric reconstruction. It is schematically presented in Figure 6.
Twelve two-dimensional distributions of selected parameters ζn (x,y) - obtained by reconstruction using the inverse Radon transform - were scanned by the sliding window as shown in Figure 6. A set of 12 parameters from 27 initial ones was selected based on a preliminary analysis of the signal to noise ratio (SNR). Their definitions are presented in Table 1. Both time domain features (ζ1 – ζ9) and frequency domain (ζ10 – ζ12) were used. Reconstructions obtained using the selected parameters enable identification of both the edges of the examined phantoms and the objects inside the structure, such as artificial or real defects. They contain a variety of information about the OUT. Data from the window that moves over 12 concurrent reconstructions (based on 12 different parameters) is given to the input of a previously trained artificial neural network (ANN). The neural model estimates the value of the final reconstruction OR for the position determined by the center of the sliding window (x, y). Three-layer network was used: 2 hidden layers with sigmoidal activation function and 1 output layer with linear activation function. The neural network was trained using the Levenberg-Marquardt method.
No. | Parameters | |
---|---|---|
Definition | Description | |
1. | ζ1 = max [s (t)] = maxVal | maximum value of s(t) signal (maxVal in Figure 3) |
2. | ζ2 = td max | time delay at which the maximum value of the signal occurs (maxPos in Figure 3) |
3. | ζ3 = min [s (t)] = min Val | minimum value of s(t) signal (minVal in Figure 3) |
4. | ζ4 = td min | time delay at which the minimum value of the signal occurs (minPos in Figure 3) |
5. | ζ5 = max [s (t)] – min [s (t)] | signal peak-to-peak value |
6. | area under absolute value of signal in time domain | |
7. | ζ7 = mean [s (t)] | mean value of the signal |
8. | ζ8 = med [s (t)] | median of sampled signal |
9. | ζ9 = sid [s (t)] | standard deviation of the sampled signal |
10. | ζ10 = 10log10 |S (f1)| | amplitude of selected harmonic f1 = 70GHz |
11. | ζ11 = 10log10 |S (f3)| | amplitude of selected harmonic f2= 250GHz |
12. | ζ12 = 10log10 |S (f5)| | amplitude of selected harmonic f3 = 1.3THz |
4 Results
In order to determine the parameters of proposed reconstruction procedure and neural model, the following criterion of reconstructed image evaluation was proposed:
where: KR – criterion for evaluation of the reconstruction result by the proposed procedure, wms – size of the sliding window,NA – number of neurons in the neural model, kζ – number of ζn parameters used for reconstruction, KFN – sum of false negative indications, KFP – sum of false positive indications.
Such simple criterion can be utilized because of quasi-homogenous character of utilized phantoms. For all considered criteria (KFN, KFP and KR), the smaller the value, the better is the result of the reconstruction. Thus, the minimum value of KR will be sought to determine the parameters of the utilized neural model. In further analysis, the above criterion will be used in a standardized form:
An experiment was conducted, in which the influence of the number of parameters used in the reconstruction on its quality was investigated. The size of sliding window ranged from 3 to 11 and total number of neurons in ANN ranged from 10 to 60. For each combination of (wms, NA, kζ), 100 networks were trained based on T-shape phantom presented in Figure 7. The validation was performed using several phantoms of different shape. Exemplary results of reconstruction using standard inverse Radon transform and proposed procedure are presented in Figure 8. In case of proposed method, one can observe blur reduction and better defined phantom edges (no oscillation at the air-phantom boundary), compared to standard procedure. This allows a more accurate determination of the size of the examined object and its details. The disadvantage of the proposed algorithm is the need to train a neural model for combining data from various parameters. An influence of number of parameters on reconstruction quality is shown in Figure 9. The order of the features used in this experiment was chosen randomly. One can observe, the reconstruction error (KRN) decreases nonlinearly with increase of kζ. Utilization of 4 parameters causes noticeably smaller KRN than in case of one parameter. Further increase of kζ doesn’t affect the reconstruction error noticeably.
5 Conclusions
The proposed THz CT system and multiparametric algorithm enable proper reconstruction of the considered objects cross-section (as shown in Figure 8) and cause blur reduction of the resulting image. Moreover, the algorithm enables reduction of reconstruction error in comparison to the case of single parameter based inversion (as shown in Figure 9). That proves correctness of proposed multiparametric approach.
The drawback of proposed algorithm is a necessity of training the neural model, thus it is not as widely applicable as transforms based inversions.
References
[1] Mittelman D., Sensing with terahertz radiation, Springer, Berlin, 2010Search in Google Scholar
[2] Palka N., Miedzińska D., Detailed non-destructive evaluation of UHMWPE composites in the terahertz range, Optical and Quantum Electronics, 2014, 46, 515-525.10.1007/s11082-013-9836-4Search in Google Scholar
[3] Yakovlev E., Zaytsev K., Dolganova I., Yurchenko S., Nondestructive evaluation of polymer composite materials at the manufacturing stage using terahertz pulsed spectroscopy, IEEE Trans. Terahertz Sci. Technol., 2015, 5, 5, 810-816.10.1109/TTHZ.2015.2460671Search in Google Scholar
[4] Ahi K., Anwar M., Advanced terahertz techniques for quality control and counterfeit detection, Proc. SPIE 9856, Terahertz Physics, Devices, and Systems X: Advanced Applications in Industry and Defense, 98560G, 2016, 10.1117/12.2228684Search in Google Scholar
[5] Lopato P., Automatic defect recognition for pulsed terahertz inspection of basalt fiber reinforced composites, Compel, 2016, 35, 4, 1346-1359.10.1108/COMPEL-09-2015-0351Search in Google Scholar
[6] Ferguson B., Wang S., Gray D., Abbott D., Zhang X.C., Towards functional 3D T-ray imaging, Physics in Medicine and Biology, 2002, 47, 3735-3742.10.1088/0031-9155/47/21/309Search in Google Scholar PubMed
[7] Wang S., Zhang X.C., Pulsed terahertz tomography, Journal of Physics D: Applied Physics, 2004, 37, R1-R36.10.1088/0022-3727/37/4/R01Search in Google Scholar
[8] Zimdars D., Fichter G., Chemovski A., Time domain terahertz computed axial tomography for nondestructive evaluation, Review of Quantitative Nondestructive Evaluation, 2009, 28, 426-433.10.1063/1.3114269Search in Google Scholar
[9] Choi H., Son J.-H., Terahertz imaging and tomography techniques, Terahertz Biomedical Science and Technology, Joo-Hiuk Son (ed.), CRC Press, 2014, 47-66.Search in Google Scholar
[10] Guillet J.P., Recur B., Frederique L., Bousquet B., Canioni L., Manek-Honninger I., Desbarats P., Mounaix P., Review of terahertz tomography techniques, Journal of Infrared, Millimeter and Terahertz Waves, 2014, 35, 382-411.10.1007/s10762-014-0057-0Search in Google Scholar
[11] Kak A.C., Slaney M., Principles of computerized tomographic imaging, IEEE Press, New York, 1988Search in Google Scholar
[12] Psuj G., Fusion of multiple parameters of signals obtained by vector magnetic flux observation forevaluation of stress loaded steel samples, International Journal of Applied Electromagnetics and Mechanics, 2015, 49, 1, 1-7.10.3233/JAE-152003Search in Google Scholar
© 2018 P. Lopato, published by De Gruyter
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.