Decay diagnosis of tree trunks using 3D point cloud and reverse time migration of GPR data

Health monitoring and disease mitigation of trees are essential to ensure the sustainability of wood industry, safety of ecosystems, and maintenance of climatic conditions. Several non-destructive testing methods have been applied to monitor and detect the decays inside the trunks. Among them, ground penetrating radar (GPR) has gained recognition due to its high efficiency and good resolution. However, due to the wide beam width of the antenna pattern and the complicated scattering caused by the trunk structure, the recorded GPR profile is far from the actual geometry of the tree trunk. Moreover, the irregular contour of the tree trunk makes traditional data processing algorithms difficult to be performed. Therefore, an efficient migration algorithm with high resolution, as well as a high accuracy survey-line positioning method for curved contour of the trunk should be developed.

In this paper, a combined approach is proposed to image the inner structures inside the irregular-shaped trunks. In the first step, the 3D contour of the targeted tree trunk is built up by a 3D point cloud technique via photographing around the trunk at various angles. Subsequently, the 2D irregular contour of the cross-section of trunk at the position of the GPR survey line is extracted by the Canny edge detection method to locate the accurate position of each GPR A-scans [1]. Thirdly, the raw GPR profile is pre-processed to suppress undesired noise and clutters. Then, an RTM algorithm based on the zero-time imaging condition is applied for image reconstruction using the extracted 2D contour [2]. Lastly, a denoising method based on the total variation (TV) regularization is applied for artifact suppression in the reconstructed images [3].

Numerical, laboratory and field experiments are carried out to validate the applicability of the proposed approach. Both numerical and laboratory experimental results show that the RTM can yield more accurate and higher resolution images of the inner structures of the tree cross section than the BP algorithm. The proposed approach is further applied to a diseased camphor tree, and an elliptical decay defect is found the in the migrated GPR image. The results are validated by a visual inspection after the tree trunk was sawed down.

Fig. 1 Field experiment. (a) Geometric reconstruction result using point cloud data, (b) migrated result by the RTM algorithm and (c) bottom view of the tree trunk after sawing down. The red and yellow ellipses indicate the cavity and the decay region in the trunk, respectively.

References:

[1] Canny, "A Computational Approach to edge detection," IEEE Transactions on Pattern Analysis and Machine Interllgent, vol. PAMI-8, no. 6, pp. 679-698, 1986, doi: 10.1109/TPAMI.1986.4767851.

[2] S. Chattopadhyay and G. A. McMechan, "Imaging conditions for prestack reverse-time migration," Geophysics, vol. 73, no. 3, pp. S81-S89, 2008, doi: 10.1190/1.2903822.

[3] L. I. Rudin, S. Osher, and E. Fatemi, "Nonlinear total variation based noise removal algorithms," Physica D, vol. 60, pp. 259-268, 1992, doi: 10.1016/0167-2789(92)90242-F.