An Experimental Study and Concept Evaluation on Tree-Interior Imaging Radar Using Sinusoidal Template-Based Focusing Algorithm.
The detection of voids or cavities within tree-bodies at microwave band has become an emerging technology for the determination of unhealthy and weak trees [1,2]. Every year tens of people are being killed all around the world due to falling tree incidents and related incidents [1-4]. Although most of the trees seem to be healthy and strong when looking from the outside, it is not easy to realize the existence of interior cavities that might occur due to fungus and worm attacks by using conventional methods such as acoustic imaging techniques , resistivity measurement methods , and thermographic inspection . It is known by the experts that decaying wood is expected to have a different moisture content and density; therefore, this part of the tree-body should have different electrical characteristics such as dielectric constant and conductivity . One of the most applied methods is the ultrasonic tomography that can sense decay structures only under special circumstances . This technique uses the mechanical behavior of wood that usually depicts a strong anisotropy such that ultrasonic pulses are being used to sense this discontinuity caused by the cavity. On the other hand, ultrasonic pulses in wood are usually strongly attenuated; therefore, this method requires high signal sensitivity and dynamic range while acquiring and processing the raw data [10,11]. Recently, using obstacle-penetrating radar technologies to detect and image cavities within a tree-body has been increasingly gaining attention thanks to recent advances in radar hardware and radar signal processing techniques [12-15]. It is eminent by the electromagnetic (EM) community that microwave imaging techniques have the advantages of providing high resolution features and requiring relatively smaller amount of power to penetrate such obstacles. Algorithms used in ground-penetrating-radar (GPR) [16-18] and through-wall-imaging-radar (TWIR) [19-21] seem to provide promising solutions to this special problem. Also, circular Radon transform could be an opportunity to focus the scattered electric field data collected from a circular measurement set-up for a different application .
In this work, we offer a practical but also a fast focusing/migration method that can be readily used for tree-interior imaging radar (TIIR) applications. We are presenting a conceptual study on TIIR by suggesting an imaging/focusing method that has been specially tailored to this specific application and also demonstrating a proof-of-concept study by providing and assessing the results of related simulations and measurements from a test-bed.
The paper is organized as follows: in the next section, we present the steps and the details of our imaging algorithm for the TIIR problem. In the third section, the validity of the proposed method is first tested with the help of simulated data that are generated with point scatterer model. In the following section, a measurement set-up for TIIR experiments is constructed and both the simulation of the experimental set-up and measurements conducted with real tree-body with cavity are presented. The results after applying the proposed methods are provided. Resultant radar images of both the simulation and the measurement data validate the effectiveness and the success of our suggested algorithm. In the final section, the work is concluded by discussing the applicability and assessment of the measurement system with the proposed algorithm.
2. Imaging/Focusing Algorithm for TIPR
The geometry of our TIIR problem can be described as illustrated in Figure 1 where the tree-body is circularly scanned by a directive antenna for a total of M distinct look-aspects covering the whole azimuth angles. The radius of the scanning circle is [R.sub.0]. The frequency-diverse backscattered electric field data are collected in the monostatic configuration of the radar. For any point scatterer, P([x.sub.n], [y.sub.n]) was assumed to be on the perimeter of the cavity within the radius [R.sub.0]; therefore, one-dimensional (1D) electric field data can be obtained by the antenna of the ith position as
[E.sup.s.sub.i](k) = [A.sub.i] exp(-j2k[r.sub.i]). (1)
Here, [r.sub.i] = sqrt[[([x.sub.n] - [x.sub.i]).sup.2] + [([y.sub.n] - [y.sub.i]).sup.2]] is the distance between the scatterer and the antenna and k is the wavenumber vector. This result is valid for the following assumptions: first, we assume isotropic radiation of the electromagnetic (EM) wave inside the tree-body such that the EM wave travels along its propagation direction. This assumption is usually true for the decayed trees since such trees usually dry in the inside due to decaying phenomenon. It is reported that the dielectric constant of dried wood material is on the order of 1.4 to 2.9 in practice . Therefore, the moisture content for the decayed trees is usually very low and does not present significant anisotropy in terms of the water content of the tree-body. This situation also helps radar signal to penetrate the tree-body better since the relative electric permittivity contrast between the air and the tree structure is not high. Secondly, the scattering structures within the tree-body are assumed to be represented with point scatterer model that is very effective and commonly used in many radar applications . Therefore, near field boundary conditions for the air void and the tree-body are not imposed. Although such modelling will provide a more realistic representation, the surface modelling of air void would be quite complex and yield a very complex calculation of the EM reflection. Therefore, we have encountered a more basic and simplified EM scattering calculation model to ease the numerical load of the problem. Furthermore, such kind of study is out of this research work. Yet, modelling boundary points such as surface of cavities with perfect point scatters has proven to be an effective way of calculating the scattering from such boundaries in many radar problems [25,26]. Thirdly, any multipath phenomenon is not accounted as all radar imaging algorithms are based on only single-reflection assumption.
Based on the above assumptions, then, (1) can be rewritten as a summation over scatterer and provides a twodimensional (2D) data by considering different look angles as the following for a finite number of point scatterers (M) that form the tree-body:
[E.sup.s](k, [phi]) = [M.summation over (i=1)] [A.sub.i] exp(-j2k[R.sub.i]), (2)
where [R.sub.i] = sqrt[[(X - [x.sub.i]).sup.2] + [(Y - [y.sub.i]).sup.2]] is the distance vector between the scatterer and all the positions of radar for whole azimuth observation angles of [phi] ranging from 1 to M and k is the wave number for the stepped frequencies for a total of N distinct frequencies. Then, X and Y vectors correspond to Cartesian coordinates of M antenna positions, accordingly as given below:
X = [R.sub.0] cos [phi],
Y = [R.sub.0] sin [phi]. (3)
Taking the one-dimensional (1D) inverse Fourier transform (IFT) of (2) along the spatial frequency axis, one can easily get the 2D range-angle data, [E.sup.s](r, [phi]). In this resultant range-angle data, any point P([x.sub.n], [y.sub.n]) exhibits a sinusoidal behavior for the 360[degrees] variation of azimuth angle [phi] because of circular scanning arrangement.
Our detection and imaging algorithm can be briefly summarized as follows: (i) after obtaining the 2D range-angle backscattered electric field data, [E.sup.s](r, [phi]), we select an image window in 2D (range, angle) domain that covers the whole geometry. It is important to note that the range-extent of the image window should be larger than or at least equal to the diameter of scanning, that is, 2[R.sub.0]. (ii) Next, we iteratively pick all the points; P([x.sub.p], [y.sub.p]) within the selected window. (iii) Then, we form a template vector for each look angle of [[phi].sub.i] by writing the following template equation:
[R.sup.n.sub.temp]([[phi].sub.n]) = sqrt[[(X - [x.sub.p]).sup.2] + [(Y - [y.sub.p]).sup.2]]. (4)
(iv) Afterwards, we go back to 2D [E.sup.s](r, [phi]) data and pick the data under ([R.sup.n.sub.temp], [[phi].sub.n]) points. At this point, a precise data interpolation scheme needs to be employed to lessen the error associated with data interpolation. (v) In the final step, we add all the data under this template vector and record the result in a new 2D image matrix at ([x.sub.p], [y.sub.p]) location. The new, reconstructed image data at ([x.sub.p], [y.sub.p]) location becomes a pixel in the final image. The algorithm successively continues until all the points in the selected image window are passed through the algorithm. The resultant I (x, y) matrix is the focused image.
3. Simulation Results
The above-explained algorithm was first tested by a simulation scenario whose illustration is given in Figure 2. The simulation scene was constructed in MATLAB programming environment  with the assumptions that were listed in the previous section. To be able to test and evaluate the performance of this technique with ideal conditions, perimeter of the tree-body (shown as black dots) was assumed to be composed of a total of 60 perfect point scatterers with the identical reflectivity amplitudes of "1," that is, independent of frequency and aspect. The radius of the tree-body is taken as 60 cm. In the simulation, there were two cavity regions; one was centered at (10 cm, 0) while the other one was located at (-20 cm, 5 cm). Both cavities were considered to be in the form of circles with diameters of 8 cm and 6 cm for the first and second void regions, respectively. The circumference of the cavities was also assumed to be composed of perfect point scatterers with equal EM reflectivities of "0.5" that are independent of frequency and angle. Dielectric constant for the tree-body was taken as 1.65 and assumed to be homogeneous within the tree-body. The radar antenna was positioned at 35 cm away from the tree's surface and moved along a circular track for a total of 360 distinct look angles to collect the backscattered data. For each of the look angles, frequency-diverse backscattered electric field was gathered for the stepped frequencies ranging from 1 to 8 GHz for the total of 100 discrete frequency points. Therefore, the backscattered electric field data were collected on the 2D frequency-angle plane such that a matrix of 100 x 360 points was obtained.
We apply the steps of our focusing algorithm that was based on summation of the data beneath the sinusoidal template. In Figure 3(a), the range profile matrix for different look angles, [E.sup.s](r, [phi]), was obtained. We have several observations about the range-angle raw data in Figure 3(a) as follows: (i) the scattering from the nearest point from the circular tree-body geometry shows up as a straight line since both the scanning path and the tree-body are in the form of circles. This line occurred at a range distance of 35 cm as expected. (ii) Scatterers on the perimeter of both cavities come out as sinusoids of different oscillating amplitudes due to different round-trip distances between the antenna and the scatterers due to circular movement of the radar while collecting the data. In fact, the peak amplitudes of these oscillating sinusoids depend on their distances from the center of rotation. For our simulation example, these numbers are 4 cm and 3 cm for the first and the second cavity, respectively. The range-angle image in Figure 3(a) exactly exhibits this behavior as we have anticipated. (iii) The widths of the sinusoids represent the diameters of the cavities. From the figure, the widths of the sinusoidal behavior are 8 cm and 6 cm that are precisely the same values of real diameters for the cavity regions. (iv) The starting point of the sinusoidal behavior actually pinpoints the exact location of the cavity: for example, if the sinusoids start from the maximum value, the cavity is located at the nearest point to the antenna's location. If the sinusoid starts from the minimum value, then, the cavity is at 180[degrees] away from the antenna. (v) Center point of the scanning will show up as a straight line in the range-angle image. Therefore, if there is a small cavity at the center, its image will be a horizontal line at the half wave between the scatterings from near and far points of the tree-body. (vi) Second straight line around 95 cm in range-angle image in Figure 3(a) corresponds to scattering from far-most point of tree-body. As the radar moves around a circle, the distance from the far-most point of the tree-body stays constant at 95 cm as it was obtained in Figure 3(a). After analyzing the range-angle image in Figure 3(a) that certainly tells us many facts, we applied the steps of the algorithm such that the final focused image I(x, y) was obtained as given in Figure 3(b). As it is clearly seen from the figure, both the tree-body and the cavity regions were successfully imaged in the resultant focused image at their correct locations in the Cartesian coordinates. Of course, the image strength of cavities is much smaller than that of tree's surface due to decaying of EM wave with 1/(4[pi][R.sup.2]) term.
The fidelity of the resultant image in our simulation is compared to other methods such as Radon transform  and ultrasonic techniques are as follows : focused image can be obtained from Radon transform simulation method that uses projection data associated with cross-sections of an object but this method requires complex calculations for obtaining the final image. Therefore, its implementation is more difficult and takes much more time. Ultrasonic simulation method also can be used for imaging but the disadvantage of this method is the weakness of the ultrasonic signals to be able to penetrate the tree-body. To ensure the required minimum detectable signal level, many receivers are usually required along the perimeter of tree-body. Our method is more practical and has faster migration method to be applied for TIIR applications.
4. Measurement Set-Up and Experimental Results
To assess the effectiveness of our TIIR detection/imaging algorithm with real experimental data, we have constructed a measurement test-bed at our anechoic chamber of Mersin University's Advanced Technologies Research Center (MEATRC) facility. The geometry of the measurement setup is shown in Figure 4. In this set-up, a Vector Network Analyzer (VNA) [Agilent ENA5071B] was utilized to generate/transmit/receive the radar signal by using an ultrawide band (UWB) double-ridged horn antenna [Geozondas GZ0126DRH] that was used in the monostatic configuration. The antenna frequency range is 1-26 GHz. Beam width of the antenna at E plane is 37.7[degrees] at 1 GHz and 22.4[degrees] 8 GHz. A turn table with automation software was employed to get different look angle profiles of the tree-body. Also, a PC has to be used to synchronously control the VNA and the turn table with a MATLAB  script.
For the presented configuration in Figure 4, a piece of pine tree trunk with outer-diameter of approximately 42 cm was used as the target to be imaged. The tree has a cavity whose center is about 2 cm away from the center of the tree with a radius of about 11 cm as pictured in Figure 5(a). Before starting the experiment, we first constructed an approximate 2D simulation of the measurement that we were going to do to visualize the results out of the experiment and to compare it to real experiment that we were going to conduct later on. In Figure 5(b), the simulation for the experiment was constructed again in MATLAB. Our antenna is assumed to be 58 cm away from the center of the tree or 36 cm away from the body of the tree trunk. The geometry for the tree trunk and the cavity was constructed in accordance with the real values as illustrated in Figure 5(b).
For the simulation, the frequency is altered from 1 to 8 GHz for a total of 100 discrete frequencies while the look angle was also varied from 1[degrees] to 360[degrees] for a total of 361 distinct aspects. The tree-body's frequency-independent reflectivity amplitude was chosen to be "1" while that of cavity was selected as "0.1" to emulate the attenuation of EM wave while travelling inside the tree-body. The 2D range-angle data after running the simulation is presented in Figure 6(a) where we can observe different scattering mechanisms easily: the reflection from the tree-body showed up as a straight line around 36 cm from the antenna as expected. The front point and the back point of the cavity experienced two sinusoids as the look angle varies along the whole azimuth. Of course, the reflection magnitude for the back point of the cavity is less than that of the front point of the cavity due to attenuation of the EM wave with distance. After applying our sinusoidal template-based focusing algorithm, the TIIR image of this simulation for the experiment was gathered as plotted in Figure 6(b). By comparing this figure with the simulation geometry in Figure 5(b), one can easily see that the suggested algorithm successfully forms the TIIR image of the scenario with good fidelity.
As the final study, we have conducted a real through-the-tree-imaging-radar experiment as pictured in Figure 5(c). The experiment was performed with the measurement test-bed at our anechoic chamber of MEATRC facility. During this experiment, the tree trunk whose picture is seen in Figure 5(a) was put on the turn table. The antenna in monostatic configuration was positioned as seen in Figure 5(c) such that the phase center of the measurement set-up is about 58 cm away from the scanning axis of the tree trunk just as similar to the simulation. By using the automation software that was adopted by our MATLAB script, turn table has been turned for the full azimuth coverage of 360[degrees] at 1[degrees] increments. For each look angle measurement, the frequency of the VNA was changed from 1 to 8 GHz for a total of 100 discrete frequency points as in the case of the simulation. Therefore, a 2D scattered electric field matrix, [E.sup.s](k, [phi]), of 100 x 360 measurement points was collected. Theoretical achievable resolution is calculated as [DELTA]r = c/(2 * B * [square root of ([[epsilon].sub.r])]) = 1.66 cm for a bandwidth of B = 7 GHz and [[epsilon].sub.r] = 1.65 (relative electric permittivity of tree trunk). Maximum range is [R.sub.max] = N * [DELTA]r = 166 cm for number of the discrete frequency points N = 100.
After applying the IFT operation along the frequency axis, 2D range-aspect data, [E.sup.s](r, [phi]), were obtained as plotted in Figure 6(c). We can notice the scattering from the tree's body as a wavy line around 36 to 40 cm since the perimeter of the tree was not aperfect circle. We also notice the scattering from cavity that experiences the same behavior of the tree's body as the look angle varies during the measurement. Although the attenuation of the EM wave inside the tree texture looks more intense than we have anticipated in the simulation, we can still sense the disturbance of the scattering energy due to existence of the cavity around 50 to 70 cm in range direction. Next, we have employed our sinusoidal template based focusing algorithm to migrate the data from range-aspect domain to the image domain such that we have obtained the final TIIR image of the tree trunk as depicted in Figure 6(d). Thanks to the focusing ability of our algorithm, the reflections from the tree-body and the cavity's perimeter have been amplified. Therefore, we can easily notice the borders of impedance changes from air to tree-body and tree-body to cavity and vice versa. Visually comparing the results in Figure 6(d) [measurement] and Figure 6(b) [simulation] with the real picture of tree trunk in Figure 5(a), one can clearly realize that resultant TIIR images successfully form the geometry of tree-body and the cavity around their correct locations such that the cavity inside the tree trunk is detected.
In this work, we have presented our recent work on tree-interior imaging radar research. On this scope, we have introduced a new algorithm based on summation of raw data over sinusoidal templates to be effectively used for TIIR or similar applications. The algorithm is specially developed for multiangle multifrequency monostatic radar set-up with circular scanning geometry of object under test. The details of the algorithm were given together with the assumptions. Then, the algorithm is first tested with simulation data and then a real experiment was conducted at Mersin University's MEATRC facility. Both the simulated and the measured results of TIIR images demonstrate the validity and the success of the proposed algorithm in detecting and localizing the cavity region inside the tree-body.
Algorithm can be modified to include the effect caused by change in the velocity of the EM waves inside the tree-body as a future study. This algorithm can also be applied to similar problems such as detecting and imaging breast tumors by employing a very similar geometry of data collection.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This work has been supported by the Scientific and Research Council of Turkey (TUBITAK) under 2211-D Ph.D. scholarship programme and funded by Mersin University Scientific Research Project Unit with Grant no. BAP 2015-TP3-1158. The authors are also thankful to Mersin University MEATRC management for their help and support while conducting the experiments.
 J. Fazio, "How to Recognize and Prevent Hazard Trees," Tree City USA Bulletin No. 15, National Arbor Day Foundation, Nebraska City, NE, USA, 1989.
 J. Albers and E. Hayes, How to Detect, Assess and Correct Hazard Trees in Recreational Areas, DNR, St. Paul, MN, USA, 1993.
 A. Brookes, "Preventing death and serious injury from falling trees and branches," Australian Journal of Outdoor Education, vol. 11, no. 2, pp. 50-59, 2007.
 J. Watt and D. J. Ball, "Trees and the risk of harm, report for the national tree safety group," Forestry Commission Rept., Middlesex University, Scotland, 2009.
 V. Bucur, Acoustics of Wood, CRC Pres, Boca Raton, FL, USA, 1995.
 S. Al Hagrey, "Electrical resistivity imaging of tree trunks," Near Surface Geophysics, vol. 4, no. 12, pp. 177-185, 2006.
 G. Catena, L. Palla, and M. Catalano, "Thermal infrared detection of cavities in trees," European Journal of Forest Pathology, vol. 20, no. 4, pp. 201-210, 1990.
 A. Franchois, "Microwave permittivity measurements of two conifers," IEEE Transactions on Geoscience and Remote Sensing, vol. 36, no. 5, pp. 1384-1395, 1998.
 V. Bucur and I. Bohnke, "Factors affecting ultrasonic measurements in solid wood," Ultrasonics, vol. 32, no. 5, pp. 385-390, 1994.
 M. Hasegawa, M. Takata, J. Matsumura, and K. Oda, "Effect of wood properties on within-tree variation in ultrasonic wave velocity in softwood," Ultrasonics, vol. 51, no. 3, pp. 296-302, 2011.
 H. Berndt, A. P. Schniewind, and G. C. Johnson, "Ultrasonic energy propagation trough wood: where, when, how much," in Proceedings of the Proc. Of 12th Int. Symposium on NDT of Wood, pp. 57-66, Sopron, 2000.
 Y. Zhou, "Microwave imaging based on wideband range profiles," Progress in Electromagnetics Research Letters, vol. 19, pp. 57-65, 2010.
 G. Lawday, J. A. Dolwin, D. Lonsdale, and J. R Barnett, "Development and use of stress wave meter, to detect the presence of decay in wood blocks," in Proceedings of the Proc. of 12th Int. Symposium on NDT of Wood, pp. 187-196, Sopron, September.
 V. Bucur, "Techniques for high resolution imaging of wood structure: A review," Measurement Science and Technology, vol. 14, no. 12, pp. R91-R98, 2003.
 B. Yilmaz, C. Ozdemir, and S. Gokkan, "A focusing algorithm for tree-penetrating radar imaging: An experimental study and concept evaluation," in Proceedings of the 2016 17th International Radar Symposium (IRS), pp. 1-4, Krakow, Poland, May 2016.
 J. R. Butnor, J. A. Doolittle, L. Kress, S. Cohen, and K. H. Johnsen, "Use of ground-penetrating radar to study tree roots in the southeastern United States," Tree Physiology, vol. 21, no. 17, pp. 1269-1278, 2001.
 C. (Ozdemir, [section]. Demirci, E. Yigit, and B. Yilmaz, "A review on migration methods in b-scan ground penetrating radar imaging," Mathematical Problems in Engineering, vol. 2014, Article ID 280738, 2014.
 C. Ozdemir, S. Demirci, E. Yigit, and A. Kavak, "A hyperbolic summation method to focus B-scan ground penetrating radar images: An experimental study with a stepped frequency system," Microwave and Optical Technology Letters, vol. 49, no. 3, pp. 671-676, 2007.
 B. Yilmaz, S. Demirci, E. Yigit, and C. Ozdemir, "An Experimental Study of Through-the-Wall Radar for Life sign Detection," in Proceedings of Progress in Electromagn, Research Symposium (PIERS '13), pp. 1602-1604, 2013.
 M. G. Amin, Through-the-Wall Radar Imaging, Taylor and Francis, Bocaraton, FL, USA, 2010.
 A. Muqaibel, A. Safaai-Jazi, A. Bayram, A. M. Attiya, and S. M. Riad, "Ultrawideband through-the-wall propagation," IEE Proceedings--Microwaves, Antennas and Propagation, vol. 152, no. 6, pp. 581-588, 2005.
 N. J. Redding and G. N. Newsam, "Inverting the Circular Radon Transform," DTSO Research Report DTSO-RR-0211, August 2001.
 http://web.hep.uiuc.edu/home/serrede/P435/Lecture_Notes/ Dielectric_Constants.pdf.
 C. Ozdemi, Inverse Synthetic Aperture Radar Imaging with Matlab Algorithms, John Wiley Sons, Hoboken, NJ, USA, 2012.
 R. Bhalla, J. Moore, and H. Ling, "Global scattering center representation of complex targets using the shooting and bouncing ray technique," IEEE Transactions on Antennas and Propagation, vol. 45, no. 12, pp. 1850-1856, 1997.
 C. Ozdemir, R. Bhalla, and H. Ling, "A radiation center representation of antenna radiation patterns on a complex platform," IEEE Transactions on Antennas and Propagation, vol. 48, no. 6, pp. 992-1000, 2000.
 MATLAB, R2015a, MathWorks Inc., Natick, MA, 2015.
Department of Electrical and Electronics Engineering, Mersin University, Ciftlikkoy 33343, Mersin, Turkey
Correspondence should be addressed to Serhat Gokkan; email@example.com
Received 6 June 2017; Accepted 14 September 2017; Published 1 November 2017
Academic Editor: Ping Feng Pai
Caption: Figure 1: The geometry of TIIR problem.
Caption: Figure 2: Simulation geometry for the TIIR algorithm.
Caption: Figure 3: Simulation results for the geometry in Figure 2: (a) range- angle image; [E.sup.s](r, [phi]), (b) focused image; I(x, y) after applying our algorithm.
Caption: Figure 4: The measurement set-up for the TIIR experiment.
Caption: Figure 5: (a) Picture of the piece of pine tree trunk, (b) simulation scenario for the experiment, and (c) picture from the real TIIR experiment.
Caption: Figure 6: Simulation and measured results for the experiment: (a) simulated range-angle data, (b) simulated TIIR image, (c) measured range-angle data, and (d) measured TIIR image.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Yilmaz, Betul; Gokkan, Serhat; Ozdemir, Caner|
|Publication:||Journal of Electrical and Computer Engineering|
|Date:||Jan 1, 2017|
|Previous Article:||Modified Tang and Pun's Current Comparator and Its Application to Full Flash and Two-Step Flash Current Mode ADCs.|
|Next Article:||Electronically Tunable Quadrature Sinusoidal Oscillator with Equal Output Amplitudes during Frequency Tuning Process.|