Imaging of Insect Hole in Living Tree Trunk Based on Joint Driven Algorithm of Electromagnetic Inverse Scattering

Trunk pests have always been one of the most important species of tree pests. Trees eroded by trunk pests will be blocked in the transport of nutrients and water and will wither and die or be broken by strong winds. Most pests are social and distributed in the form of communities inside trees. However, it is difficult to know from the outside if a tree is infected inside. A new method for the non-invasive detecting of tree interiors is proposed to identify trees eroded by trunk pests. The method is based on electromagnetic inverse scattering. The scattered field data are obtained by an electromagnetic wave receiver. A Joint-Driven algorithm is proposed to realize the electromagnetic scattered data imaging to determine the extent and location of pest erosion of the trunk. This imaging method can effectively solve the problem of unclear imaging in the xylem of living trees due to the small area of the pest community. The Joint-Driven algorithm proposed by our group can achieve accurate imaging with a ratio of pest community radius to live tree radius equal to 1:60 under the condition of noise doping. The Joint-Driven algorithm proposed in this paper reduces the time cost and computational complexity of tree internal defect detection and improves the clarity and accuracy of tree internal defect inversion images.


Introduction
According to the Global Forest Resources Assessment report released by FAO in 2020, there are 4.06 billion hectares of forest in the world and, since 2015, an average of 10 million hectares of forestland has disappeared every year. The annual outbreak of forest pests will damage about 35 million hectares of forests. Some insects eat leaves and branches, but some trunk-boring pests such as red fat sacs live in the trunk and cannot be found in a timely manner, causing the deaths of trees and economic losses. Therefore, it is necessary to detect pest communities in the xylems of living trees in order to timely prevent and control pests and diseases and protect forest resources.
Mainstream methods for the detection of internal defects in living tree trunks include stress wave, ultrasonic, and computed tomography (CT) scanning [1][2][3]. However, most of these methods have their drawbacks [4,5]. For example, the stress wave method involves nailing into the trunk of the tree due to its detection characteristics [6]. This detection method will cause damage to the tree and cannot be non-invasive. The ultrasonic inspection process is susceptible to interference from the external environment, and coupling agents may cause environmental pollution [7]. CT equipment is costly and can easily cause radiation hazards to researchers [8,9]. Although nuclear magnetic resonance imaging has proven to be one of the most accurate and powerful characterization techniques in applications such as the diagnosis of wounded branch structures in beech trees [10], its high cost makes it unfeasible in most applications.
At present, electromagnetic inverse scattering is a good method for non-invasive testing. Electromagnetic inverse scattering (EMIS) [11] refers to the use of measured field data to invert the electrical property parameters of a scatterer in the detection region.
A large body of research has been accumulated on the EMIS problem, which can be broadly classified into two categories. One type is the linearization of nonlinear problems, represented by the Born and Rytov approximations [12,13], which have very good imaging results for Low dielectric constants. Another type of inversion method is the nonlinear method [14], which uses an iterative idea to improve the inversion accuracy through continuous iteration; the main methods are the Contrast Source Inversion (CSI) and the Subspace Optimization Method (SOM) [15,16]. Although these nonlinear iterative methods improve accuracy compared to linear approximation methods, they have disadvantages such as sensitivity to initial values and slow convergence [17]. See Table 1 for details. Although electromagnetic inverse scattering has a wide range of applications in engineering measurement and medical fields [18,19], when it is used to detect the interior of trees we find that, with larger trees, the relatively smaller proportion of the pest community area will make inversion imaging more difficult. This is because the electromagnetic inverse scattering is nonlinear, which easily makes the program fall into an error solution that differs from the accurate algorithm. Moreover, the scattering equation is highly illconditioned; that is, a small error in the input data may cause a great change in the output results. The traditional methods, whether iterative or non-iterative, do not have a very strong data processing capacity, so retrieving the model with a small proportion of pest communities is difficult.
The Deep Convolutional network is improved by studying the Contrast Source Inversion algorithm, and the Super-Resolution network is subsequently introduced. A Jointdriven Super-Resolution algorithm with a low computational load and high resolution is proposed to solve the smaller proportion of pest communities in living trees.

Electromagnetic Inverse Scattering Formulation
The algorithm is tested from a two-dimensional (2-D) perspective [20], and Figure 1 shows the model legend of our EMIS, constructed as a living tree-pest model. Our group will start with Maxwell's equations [21] by constructing scattering theoretical equations to obtain the parameter optimization ideas. To find the solution of the integral equations and investigate two major types of model-driven and data-driven inversion algorithms, our group will add Super-Resolution networks based on issues such as noise clutter during imaging to form a model-driven deep learning super-resolution-based inversion algorithm [22][23][24][25].
In this paper, we use a two-dimensional plane wave as the incident wave to detect the unknown domain, denoted as E s , and the scattered field, denoted as E i . The forward propagation of electromagnetic waves is composed of two equations known as the Lippmann-Schwinger equation [26]. The first is the equation of state, which describes the interaction of the wave scatters: where ω = 2π f is the angular frequency, µ 0 is the air magnetic permeability, and G(r, r ) denotes the Green's function. For a two-dimensional plane wave, the G(r, r ) is defined as where H (2) 0 denotes the second class of Hankel functions [27], k 0 denotes the number of waves in free space, r = (x , y ) is the source point of the object domain D,r = (x, y) denotes the position vector at the receiver, E t denotes the total field, expressed as E t = E s + E i , and ε r (r ) is the relative permittivity of the target scatter.  Figure 1. Electromagnetic inverse scattering model.
In this paper, we use a two-dimensional plane wave as the incident wave to detect the unknown domain, denoted as , and the scattered field, denoted as . The forward propagation of electromagnetic waves is composed of two equations known as the Lippmann-Schwinger equation [26]. The first is the equation of state, which describes the interaction of the wave scatters: where = 2 is the angular frequency, is the air magnetic permeability, and ( , ) denotes the Green's function. For a two-dimensional plane wave, the ( , ) is defined as where ( ) denotes the second class of Hankel functions [27], denotes the number of waves in free space, = ( , ) is the source point of the object domain D, = ( , ) denotes the position vector at the receiver, denotes the total field, expressed as = + , and ( ) is the relative permittivity of the target scatter.
The control equations can then be written in two non-permissible forms. The first is called the field-type equation, where the electric field involves two equa- The second equation, known as the data equation, describes the equivalent current radiation process for the scattered field: In (3), the physical meaning of −iωε 0 (ε r (r ) − 1) is the inductive contrast current density, but we define the standardized contrast current density as J(r) = (ε r (r) − 1)E t (r), making variable k 0 = ω √ µ 0 ε 0 , express ε r (r) − 1 as contrast χ(r), and introducing operators G s (·) and G D (·).
The control equations can then be written in two non-permissible forms. The first is called the field-type equation, where the electric field involves two equations: The second is called the source-type equation, in which the current involves two equations: The electromagnetic inverse scattering imaging algorithm is based on the above two types of equations to solve for the target cross-section parameters.

Joint-Driven Algorithm
To overcome the nonlinearity and high morbidity of inverse scattering and to solve the problem that the pest community accounts for a small proportion of the cross-sectional area of living trees, a Joint-Driven algorithm is proposed by combining the data-driven algorithm with the model-driven algorithm. The advantage of the model-driven algorithm in solving problems lies in accurate mathematical modeling based on the mechanism information of the problem itself, while the disadvantage is its excessive dependence on the initial value. The data-driven algorithm relies on the extraction of feature information in the data to give full play to the advantages of big data. It can only solve problems by constantly learning the relationship between data. Therefore, this paper introduces the model-driven mechanism into the data-driven, integrating the advantages of both.
Because the non-destructive testing of trees using electromagnetic waves is almost carried out in the field, it causes the data we obtain to contain noise. To improve the resolution, increase the clarity, and optimize the inversion results, we add a super-resolution network after the inversion and train it in a targeted way so that it can match the optimized inversion results. We call this behavior the Joint-Driven Super-Resolution algorithm.
Neural networks are typically data-driven algorithms, with Deep Convolutional neural networks being the most capable of extracting and classifying data features [28]. The convolutional neural network structure in this paper is shown in Figure 2 and is a typical ConvNets network, consisting of four types of layers: input layer, convolutional layer, pooling layer, and fully connected layer.
The electromagnetic inverse scattering imaging algorithm is based on the above two types of equations to solve for the target cross-section parameters.

Joint-Driven Algorithm
To overcome the nonlinearity and high morbidity of inverse scattering and to solve the problem that the pest community accounts for a small proportion of the cross-sectional area of living trees, a Joint-Driven algorithm is proposed by combining the data-driven algorithm with the model-driven algorithm. The advantage of the model-driven algorithm in solving problems lies in accurate mathematical modeling based on the mechanism information of the problem itself, while the disadvantage is its excessive dependence on the initial value. The data-driven algorithm relies on the extraction of feature information in the data to give full play to the advantages of big data. It can only solve problems by constantly learning the relationship between data. Therefore, this paper introduces the model-driven mechanism into the data-driven, integrating the advantages of both.
Because the non-destructive testing of trees using electromagnetic waves is almost carried out in the field, it causes the data we obtain to contain noise. To improve the resolution, increase the clarity, and optimize the inversion results, we add a super-resolution network after the inversion and train it in a targeted way so that it can match the optimized inversion results. We call this behavior the Joint-Driven Super-Resolution algorithm.
Neural networks are typically data-driven algorithms, with Deep Convolutional neural networks being the most capable of extracting and classifying data features [28]. The convolutional neural network structure in this paper is shown in Figure 2 and is a typical ConvNets network, consisting of four types of layers: input layer, convolutional layer, pooling layer, and fully connected layer.

Convolution
Convolution Convolution Fully-connected E s Figure 2. Structure of the convolutional neural network in this paper.
is fed into the network as a 32 × 32 matrix form as an input into the convolutional network. The perceptual field size in the network is 5 × 5, and the feature data is evaluated in the pooling layer with a 2 × 2 pooling operation. Finally, a 100 × 100 regression value is involved as the output for imaging through the fully connected layer.

Unite the Contrast Source Inversion with a Deep Convolutional Network
To prevent the occurrence of the overfitting phenomenon in the process of neural network training-that is, when the neural network has achieved good detection accuracy on the training data, but the detection accuracy on the test data set is greatly decreased-it E s is fed into the network as a 32 × 32 matrix form as an input into the convolutional network. The perceptual field size in the network is 5 × 5, and the feature data is evaluated in the pooling layer with a 2 × 2 pooling operation. Finally, a 100 × 100 regression value is involved as the output for imaging through the fully connected layer.

Unite the Contrast Source Inversion with a Deep Convolutional Network
To prevent the occurrence of the overfitting phenomenon in the process of neural network training-that is, when the neural network has achieved good detection accuracy on the training data, but the detection accuracy on the test data set is greatly decreased-it is necessary to normalize the cost function. In this paper, L2 normalization is adopted to improve the cost function.
If λ 2n ∑ ω ω 2 is the normalization term and the cost function before improvement is defined as C 0 , then the cost function after specification can be written as The partial derivative of C concerning the weight and bias in the network is ∂C ∂ω ∂C ∂b In (12), it shows that the partial derivative of the cost function L2 concerning bias does not change after specification and the biased gradient descent rule does not change, while the weight learning rule becomes In the learning process, the weight is adjusted by the factor 1 − η m . In the cost function, the learning rate η and parameter λ are hyperparameters. Although they are not directly related to the establishment of the Deep Convolutional neural network, the selection of hyperparameters λ and η has a great relationship with the training speed and imaging effect.
To speed up the training process of a Deep Convolutional neural network and optimize the selection of hyperparameters, the Contrast Source Inversion is combined with the cost function of a Deep Convolutional neural network. The core of the Contrast Source Inversion is to minimize the objective function [29,30], which is in the form of The aim of the Deep Convolutional neural network is to acquire the appropriate weight and bias by minimizing the cost function, so the normalization of the Contrast Source Inversion is introduced into the Deep Convolutional neural network and the value of the hyperparameter λ is determined by the mechanism of information of the Contrast Source Inversion. Formally, the value of λ is changed from a constant value to a dynamic real value determined by the Contrast Source Inversion to form a Jointly-Driven deep learning network. Therefore, the cost-function form of a model-driven deep learning network is: In (15), s represents the j group of data and the sum of the binary norm of scattering field data obtained by K receivers in each group of data. Although E s k is known, s is a dynamic real value. λ 0 is a constant weighting coefficient. In addition to preventing the Deep Convolutional neural network from over-fitting, λ 0 can also prevent abnormal training data input from causing Σ j Σ k E s k 2 s mutation, thus reducing the training effect of Deep Convolutional neural network.

Analysis and Optimization of Weights
To speed up training, weights and bias are optimized. Traditional convolutional networks are selected by Gaussian random variables, and the mean of normalized Gaussian distribution is 0 and the standard deviation is 1 [31]. Since z = ∑ ωx + b, z itself obeys the Gaussian distribution after initializing ω and b using gaussian random variables, as shown in Figure 3.

Analysis and Optimization of Weights
To speed up training, weights and bias are optimized. Traditional convolutional networks are selected by Gaussian random variables, and the mean of normalized Gaussian distribution is 0 and the standard deviation is 1 [31]. Since = ∑ + , z itself obeys the Gaussian distribution after initializing ω and b using gaussian random variables, as shown in Figure 3. As can be seen in Figure 3, the slope of the curve is small and the overly large input values make ( ), the output data, almost saturated. This characteristic makes the weight updating process change very slowly. To improve the learning speed of the network, this paper sets the mean of the Gaussian random distribution to 0 and the standard deviation to 1 , where is the number of input neurons. The improved distribution is shown in Figure 4.

Super-Resolution Network-Assisted Imaging
When the proportion of pest communities in the cross-section of trees is less than 10%, the inversion results will be noisy. If the proportion of the pest community in the cross-section of the tree becomes smaller and smaller, the noise of the inversion image becomes denser and denser. To solve this problem, an image optimization and image reconstruction technology is introduced: Real-ESRGAN [32]. The Super-Resolution algorithm can not only improve the resolution of the inversion image but also optimize the image quality and remove the dryness of the image. Real-ESRGAN uses a U-Net discriminator with spectral normalization [33][34][35], which can identify where the noise is and where the pest community is located. Its network structure is shown in Figure 5. As can be seen in Figure 3, the slope of the curve is small and the overly large input values make σ(z), the output data, almost saturated. This characteristic makes the weight updating process change very slowly. To improve the learning speed of the network, this paper sets the mean of the Gaussian random distribution to 0 and the standard deviation to 1 / √ n in , where n in is the number of input neurons. The improved z distribution is shown in Figure 4.

Analysis and Optimization of Weights
To speed up training, weights and bias are optimized. Traditional convolutional networks are selected by Gaussian random variables, and the mean of normalized Gaussian distribution is 0 and the standard deviation is 1 [31]. Since = ∑ + , z itself obeys the Gaussian distribution after initializing ω and b using gaussian random variables, as shown in Figure 3. As can be seen in Figure 3, the slope of the curve is small and the overly large input values make ( ), the output data, almost saturated. This characteristic makes the weight updating process change very slowly. To improve the learning speed of the network, this paper sets the mean of the Gaussian random distribution to 0 and the standard deviation to 1 , where is the number of input neurons. The improved distribution is shown in Figure 4.

Super-Resolution Network-Assisted Imaging
When the proportion of pest communities in the cross-section of trees is less than 10%, the inversion results will be noisy. If the proportion of the pest community in the cross-section of the tree becomes smaller and smaller, the noise of the inversion image becomes denser and denser. To solve this problem, an image optimization and image reconstruction technology is introduced: Real-ESRGAN [32]. The Super-Resolution algorithm can not only improve the resolution of the inversion image but also optimize the image quality and remove the dryness of the image. Real-ESRGAN uses a U-Net discriminator with spectral normalization [33][34][35], which can identify where the noise is and where the pest community is located. Its network structure is shown in Figure 5.

Super-Resolution Network-Assisted Imaging
When the proportion of pest communities in the cross-section of trees is less than 10%, the inversion results will be noisy. If the proportion of the pest community in the crosssection of the tree becomes smaller and smaller, the noise of the inversion image becomes denser and denser. To solve this problem, an image optimization and image reconstruction technology is introduced: Real-ESRGAN [32]. The Super-Resolution algorithm can not only improve the resolution of the inversion image but also optimize the image quality and remove the dryness of the image. Real-ESRGAN uses a U-Net discriminator with spectral normalization [33][34][35], which can identify where the noise is and where the pest community is located. Its network structure is shown in Figure 5. The U-Net discriminator focuses not only on the image as a whole but also on the characteristics of each part, the most important for inverse imaging being the shape and position of the scatter.

Evaluation Indicators
The results are evaluated using Mean Square Error (MSE) and Image Intersection Over Union (IOU) as a criterion for the accuracy of our single inversion images [36], where IOU is defined as is the area of internal pests in the scatterer in the inversion diagram and is the area of the internal pest community set in the living tree-borers model. It ·is specified that a single pest detection is judged to be accurate when > 0.87; ideally, = 1. To test the performance of the algorithm, this paper sets up multiple groups of experimental data to test the algorithm and calculates the accuracy rate of the algorithm in different detection environments. The accuracy rate formula is as follows.
represents the detection accuracy of the algorithm for all test data, is the number of all test results that are judged to be accurate, and is the total number of test data.

Experiment Condition
Considering that the detection of trees is mostly carried out outdoors, in order to better simulate the noise interference in the field detection in real life, a method of artificially adding interference factors was proposed to simulate the detection of trees in the field. The noise construction equation is shown below.
where is our constructed noise scattered field data, is the scattered field data received by the receiver, is the number of transmitters, is the number of receivers, and is the noise matrix, which is structured as follows: The U-Net discriminator focuses not only on the image as a whole but also on the characteristics of each part, the most important for inverse imaging being the shape and position of the scatter.

Evaluation Indicators
The results are evaluated using Mean Square Error (MSE) and Image Intersection Over Union (IOU) as a criterion for the accuracy of our single inversion images [36], where IOU is defined as S i is the area of internal pests in the scatterer in the inversion diagram and S f is the area of the internal pest community set in the living tree-borers model. It ·is specified that a single pest detection is judged to be accurate when IOU > 0.87; ideally, IOU = 1. To test the performance of the algorithm, this paper sets up multiple groups of experimental data to test the algorithm and calculates the accuracy rate of the algorithm in different detection environments. The accuracy rate formula is as follows.
A cc represents the detection accuracy of the algorithm for all test data, N tp is the number of all test results that are judged to be accurate, and N t is the total number of test data.

Experiment Condition
Considering that the detection of trees is mostly carried out outdoors, in order to better simulate the noise interference in the field detection in real life, a method of artificially adding interference factors was proposed to simulate the detection of trees in the field. The noise construction equation is shown below. where E n is our constructed noise scattered field data, E s 0 is the scattered field data received by the receiver, N i is the number of transmitters,N s is the number of receivers, and n 0 is the noise matrix, which is structured as follows: where n l is the noise factor, which can be changed by changing the value of n l , and n 0 is the N i × N s matrix. Our final input scattered field data for training E s is To better describe the ratio of pest community size to the cross-sectional area of wood, we call it the Inversion Ratio, and its formula is as follows: where r p represents the radius of the pest community in m and r t is the radius of the living tree in meters. This section uses simulated electromagnetic inverse scattering data to validate our proposed method with the simulation model parameters set as shown in the Table 2 below. The relative permittivity of air 1 The internal relative dielectric constant of tree 7 Noise factor n l 0.2 Our group used finite element simulations to obtain 18,000 sets of scattering data in each of the two states. To ensure generalization, we selected a random sample of 15,000 sets of data in the dataset to train our proposed combination of model-driven and data-driven methods. Another 500 sets of data in the remaining 3000 sets were randomly selected for testing.
Within 2 m × 2 m square areas, white indicates air; brown is a two-dimensional cross-section of the xylem, a circle of radius 0.6 m with a relative permittivity of 7; and black indicates an internal trunk pest community with a relative permittivity of 60.

Relationship between Accuracy and InversionRatio
Using the Deep Convolutional algorithm, experiments were carried out on four models with different sizes but the same Inversion Ratio. The results are shown in Figure 6.
The accuracy of the inversion results was evaluated according to Formula (16). The results are shown in Table 3, which demonstrates that under the condition of maintaining the Inversion Ratio, changing the radius will not affect the accuracy value. The accuracy of the inversion results was evaluated according to Formula (16). The results are shown in Table 3, which demonstrates that under the condition of maintaining the , changing the radius will not affect the accuracy value.

Detection of Pest Communities with Different Values
The Contrast source Inversion, Deep Convolution algorithm, Joint-Driven algorithm, and Joint-Driven Super-Resolution algorithm were used to detect pest communities with different . The results are shown in Figures 7 and 8. The is calculated according to Formula (16) to evaluate the accuracy of inversion results, and the results are shown in Table 4.     Table 4.  Table 4.

When
= 1: 60, the Contrast Source Inversion cannot reverse a regular graph, resulting in = error. When the inversion ratio was 1:30 to 1:50, the pest community size retrieved by the Contrast Source Inversion algorithm was the same. When = 1: 10 and = 1: 20, the Contrast Source Inversion algorithm imaging produces a certain "halo." With a threshold value of 80% of the highest relative dielectric constant value in the imaging image, only the red part of the image can

Discussion of Experimental Results
When Inversion Ratio = 1 : 60, the Contrast Source Inversion cannot reverse a regular graph, resulting in IOU = error. When the inversion ratio was 1:30 to 1:50, the pest community size retrieved by the Contrast Source Inversion algorithm was the same. When Inversion Ratio = 1 : 10 and Inversion Ratio = 1 : 20, the Contrast Source Inversion algorithm imaging produces a certain "halo." With a threshold value of 80% of the highest relative dielectric constant value in the imaging image, only the red part of the image can be observed.
Regarding the Deep Convolutional algorithm: when Inversion Ratio = 1 : 20, 1:30, and 1:40, its inversion insect pest community size is unchanged, and its algorithm is insensitive to size change, which is prone to error.
When the Joint-Driven algorithm inverts the model of Inversion Ratio = 1 : 10, 1 : 20, the inversion results are clear and accurate. Moreover, the Joint-Driven algorithm can accurately reverse the size of the pest community. Although there are interference points in the inversion results, the location of the pest community can also be distinguished.
The Joint-Driven Super-Resolution algorithm can not only de-noise inverted images, but also further highlights the position of the pest community, further improves the inversion effect, and is more conducive to judging the internal conditions of trees by further optimizing the image. Figure 10 demonstrates the Inversion Ratio = 1 : 30 single-partial pest detection process because when the Inversion Ratio is between 1:30 and 1:60, Contrast Source Inversion cannot correctly show the shape of the scatter. The Contrast Source Inversion, Deep Convolution, and Joint-Driven methods are compared to analyze the stability of model-driven deep learning networks in the iterative process.  Figure 10a shows that the Contrast Source Inversion tends to stabilize after around 500 iterations, with the number of iterations being higher and more time-consuming. Figure 10b shows that the Deep Convolutional algorithm needs to iterate about 350 times to reach the stable range, and the training process is not smooth enough. Figure 10c shows that the Joint-Driven algorithm only needs about 60 iterations to reach the stable state, which greatly reduces the number of iterations, improves the imaging time, and smooths the training process.
According to Equation (17), the detection accuracy statistics of each inversion algorithm were carried out separately for each of the three radii, and Figure 11 shows the detection accuracy of each algorithm in 300 sets of single-pest inversion tests.  Figure 10a shows that the Contrast Source Inversion tends to stabilize after around 500 iterations, with the number of iterations being higher and more time-consuming. Figure 10b shows that the Deep Convolutional algorithm needs to iterate about 350 times to reach the stable range, and the training process is not smooth enough. Figure 10c shows that the Joint-Driven algorithm only needs about 60 iterations to reach the stable state, which greatly reduces the number of iterations, improves the imaging time, and smooths the training process.
According to Equation (17), the detection accuracy statistics of each inversion algorithm were carried out separately for each of the three radii, and Figure 11 shows the detection accuracy of each algorithm in 300 sets of single-pest inversion tests.  Figure 10a shows that the Contrast Source Inversion tends to stabilize after around 500 iterations, with the number of iterations being higher and more time-consuming. Figure 10b shows that the Deep Convolutional algorithm needs to iterate about 350 times to reach the stable range, and the training process is not smooth enough. Figure 10c shows that the Joint-Driven algorithm only needs about 60 iterations to reach the stable state, which greatly reduces the number of iterations, improves the imaging time, and smooths the training process.
According to Equation (17), the detection accuracy statistics of each inversion algorithm were carried out separately for each of the three radii, and Figure 11 shows the detection accuracy of each algorithm in 300 sets of single-pest inversion tests. Figure 11. Accuracy of four inversion algorithms for detection. Figure 11. Accuracy of four inversion algorithms for detection.
As can be seen from Figure 11, the Contrast Source Inversion can only invert the living tree burrow scatters with an Inversion Ratio greater than 1:30 and its application range is limited. Although the Deep Convolutional network algorithm can invert living tree burrow scatters with an Inversion Ratio = 1 : 60, its accuracy is only about 70%, which cannot meet the actual requirements. The Joint-Driven algorithm proposed in this paper not only solves the problem of imaging the wood cross-sectional area with a small insect colony proportion, but it also has an inversion accuracy of up to 90% after the optimization of a specially-trained super-resolution network.
To sum up, the four methods used are summarized and the results are shown in the following Table 5. When Inversion Ratio = 1 : 10 and Inversion Ratio = 1 : 20, the Contrast Source Inversion algorithm can reverse meaningful results. If the Inversion Ratio continues to change, the results generated by the Contrast Source Inversion algorithm will lose their reference value. Therefore, the calculation error is represented by 'Non'.

Conclusions
In this paper, a model-driven deep learning Super-Resolution inversion algorithm is proposed to solve the problem of high noise and poor imaging in electromagnetic wave detection of tree pest communities. By studying the propagation process of electromagnetic waves in scatters and combining the advantages of the Contrast Source Inversion, Deep Convolutional network, and Super-Resolution network, a Joint-Driven Super-Resolution algorithm is proposed. This algorithm overcomes the problem that the selection of neural network structure and the process of parameter optimization depend too much on the experience of the experimenter by introducing the Contrast Source Inversion, so that the neural network can better fit the nonlinear problem and overcome the ill-posed problem by learning a large amount of data. The experiments are carried out by continuously reducing the radius of the model pest community and using the Contrast Source Inversion, the Convolutional Network algorithm, the Joint-Driven algorithm, and the Joint-Driven Super-Resolution algorithm. Our solution to the problem of imaging tiny high-contrast scatters is provided.