METHOD FOR DETECTING THICKNESS OF TREE CANOPY BASED ON ULTRASONIC ECHO SIGNAL

Abstract

The present invention relates to a method for detecting a thickness of a tree canopy based on an ultrasonic echo signal, and belongs to the technical field of agricultural machinery information sensing and detection. The present invention analyzes and discriminates a first effective peak and a last effective peak of the ultrasonic echo signal, obtains a first effective peak time and a last effective peak time, and calculates the thickness of the tree canopy according to a provided calculation formula, thus realizing the direct detection of the thickness of the tree canopy. This method has characteristics such as a high accuracy and a wide application range, and is applicable to the field of detection of the thickness of a tree canopy in industries such as agricultural machinery and forestry machinery.

Claims

1. A method for detecting a thickness of a tree canopy based on an ultrasonic echo signal, comprising: sending, by an ultrasonic sensor, an ultrasonic emission signal toward the tree canopy; forming the ultrasonic echo signal by the tree canopy reflecting the ultrasonic emission signal; and acquiring the ultrasonic echo signal, the ultrasonic echo signal being composed of voltage data and corresponding time data thereof, wherein the voltage data generates an array V.sub.i=[V.sub.1, V.sub.2, V.sub.3, . . . , V.sub.n], the corresponding time data generates an array T=[T.sub.1, T.sub.2, T.sub.3, . . . , T.sub.n], i=1, 2, 3, . . . , n, and a calculation formula of the thickness L of the tree canopy is as follows: $L=k\frac{v_o\left(T_y-T_x\right)}{2}$ wherein L is the thickness of the tree canopy, in meters, v.sub.o is a local ultrasonic propagation rate, in meters/second, T.sub.x is a first effective peak time of the ultrasonic echo signal, in seconds, T.sub.y is a last effective peak time of the ultrasonic echo signal, in seconds, k is a correction coefficient, when a leaf area index (LAI) of the tree canopy meets 1.5≤LAI≤4.0, k=1, when the leaf area index of the tree canopy is $\mathrm{LAI}<1.5,k=\frac{.Math.\mathrm{LAI}-1.5.Math.}{3}+1,$  and when the leaf area index of the tree canopy is $\mathrm{LAI}>4,k=\frac{.Math.\mathrm{LAI}-4.Math.}{2}+1.$

2. The method for detecting the thickness of the tree canopy based on the ultrasonic echo signal according to claim 1, wherein the first effective peak time of the ultrasonic echo signal is acquired by: discriminating V.sub.i in the array V sequentially starting from i=2, as i gradually increases, when a first V.sub.i meets both V.sub.i>V.sub.i−1 and V.sub.i>MAX(V.sub.i+1, V.sub.i+2 . . . , V.sub.i+m); setting the first V.sub.i as a first effective peak V.sub.x of the ultrasonic echo signal; and setting time data corresponding to the V.sub.x as the first effective peak time T.sub.x, wherein i is 2, 3, 4, . . . , n, MAX(V.sub.i+1, V.sub.i+2, . . . , V.sub.i+m) is a maximum value in an array [V.sub.i+1, V.sub.i+2, . . . , V.sub.i+m], and m is an effective number of the first effective peak.

3. The method for detecting the thickness of the tree canopy based on the ultrasonic echo signal according to claim 2, wherein the effective number of the first effective peak is m=k.sub.1f, k.sub.1 is a sampling duration, and k.sub.1=0.001-0.0005, in seconds, and f is a sampling frequency of the ultrasonic echo signal, in hertz.

4. The method for detecting the thickness of the tree canopy based on the ultrasonic echo signal according to claim 1, wherein the last effective peak time of the ultrasonic echo signal is acquired by: discriminating V.sub.j in the array V sequentially starting from j=n−1, as j gradually decreases, when a first V.sub.j meets both V.sub.j>MAX (V.sub.j−1, V.sub.j−2, . . . , V.sub.j−s) and V.sub.j>MAX(V.sub.j+1, . . . , V.sub.n); setting the first V.sub.j as a last effective peak V.sub.y of the ultrasonic echo signal; and setting time data corresponding to the V.sub.y as the last effective peak time T.sub.y, wherein j is 1, 2, 3, . . . , n−1, MAX(V.sub.j−i, V.sub.j−2, . . . , V.sub.j−s) is a maximum value in an array [V.sub.j−i, V.sub.j−2, . . . , V.sub.j−s], and MAX(V.sub.j+1, . . . , V.sub.n) is a maximum value in an array [V.sub.j+1, . . . , V.sub.n], and s is an effective number of the last effective peak.

5. The method for detecting the thickness of the tree canopy based on the ultrasonic echo signal according to claim 4, wherein the effective number of the last effective peak is s=k.sub.2 f, k.sub.2 is a sampling duration, and k.sub.2=0.0005-0.0001, in seconds, and f is a sampling frequency of the ultrasonic echo signal, in hertz.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0015] The present invention will be further described below with reference to accompanying drawings and embodiments, wherein

[0016] FIG. 1 is a schematic diagram of a solution for detecting the thickness of a tree canopy using an ultrasonic sensor according to an embodiment of the present invention;

[0017] FIG. 2 is a schematic diagram of an ultrasonic emission signal and an ultrasonic echo signal according to an embodiment of the present invention; and

[0018] FIG. 3 is a schematic diagram of an ultrasonic echo signal according to an embodiment of the present invention.

[0019] In the drawings: 1. ultrasonic sensor, 2. thickness L of a tree canopy, 3. tree canopy, 4. ultrasonic emission signal, 5. ultrasonic echo signal.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0020] The technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention.

[0021] As shown in FIG. 1, it is a schematic diagram of a solution for detecting the thickness of a tree canopy using an ultrasonic sensor according to an embodiment of the present invention. The ultrasonic sensor 1 sends an ultrasonic emission signal 4 toward a tree canopy 3, the tree canopy 3 reflects the ultrasonic emission signal to form an ultrasonic echo signal 5, and the thickness 2 of the tree canopy in the detection direction of the ultrasound emitted by the ultrasonic sensor 1 is shown in FIG. 1.

[0022] A typical ultrasonic emission signal 4 and ultrasonic echo signal 5 are shown in FIG. 2. The ultrasonic emission signal 4 and ultrasonic echo signal 5 each are composed of voltage data within a certain period of time, and corresponding time data thereof. The ultrasonic emission signal 4 is composed of voltage data between time T.sub.a and time T.sub.b, and corresponding time data thereof, and the ultrasonic echo signal 5 is composed of voltage data between time T.sub.c and time T.sub.d, and corresponding time data thereof.

[0023] The typical ultrasonic echo signal 5 is shown in FIG. 3, and the ultrasonic echo signal 5 is composed of the voltage data between the time T.sub.c and the time T.sub.d. The voltage data generates an array V=[V.sub.1; V.sub.2; V.sub.3; . . . ; V.sub.n], the time data generates an array T=[T.sub.1; T.sub.2; T.sub.3; . . . ; T.sub.n], and i=1, 2, 3, . . . , n.

[0024] On the basis of the above, the detection and calculation of the thickness 2 of the tree canopy are performed according to the array V and the array T, and specific steps are as follows.

[0025] In step 1, the first effective peak and the first effective peak time are acquired from the ultrasonic echo signal: Starting from i=2, V.sub.i in the array V is sequentially discriminated. As i gradually increases, when a first V.sub.i meets both V.sub.i>V.sub.i−1 and V.sub.i>MAX(V.sub.i+1, V.sub.i+2, . . . , V.sub.i+m), the first V.sub.i is set as the first effective peak V.sub.x of the ultrasonic echo signal, and the time data corresponding to the V.sub.x is set as the first effective peak time T.sub.x, as shown in FIG. 3.

[0026] Where i is 2, 3, 4, . . . , n, MAX(V.sub.i+1, V.sub.i+2, . . . , V.sub.i+m) is the maximum value in the array [V.sub.i+1, V.sub.i+2, . . . , V.sub.i+m], m is the effective number of the first effective peak, m=k.sub.1f, k.sub.1 is a sampling duration, and k.sub.1=0.001-0.0005, in seconds; and f is a sampling frequency of the ultrasonic echo signal, in hertz.

[0027] In step 2, the last effective peak and the last effective peak time are acquired from the ultrasonic echo signal: Starting from j=n−1, V.sub.j in the array V is sequentially discriminated. As j gradually decreases, when a first V.sub.j meets both V.sub.j>MAX(V.sub.j−1, V.sub.j−2, . . . , V.sub.j−s) and V.sub.j>MAX(V.sub.j+1, . . . , V.sub.n), the first V.sub.j is set as the last effective peak V.sub.y of the ultrasonic echo signal, and the time data corresponding to the V.sub.y is set as the last effective peak time T.sub.y, as shown in FIG. 3.

[0028] Where j is 1, 2, 3, . . . , n−1; MAX(V.sub.j−1, V.sub.j−2, . . . , V.sub.j−s) is the maximum value in the array [V.sub.j−i, V.sub.j−2, . . . , V.sub.j−s], and MAX(V.sub.j+1, . . . , V.sub.n) is the maximum value in the array [V.sub.j+1, . . . , V.sub.n], s is the effective number of the last effective peak, s=k.sub.2f, k.sub.2 is a sampling duration, and k.sub.2=0.0005-0.0001, in seconds; and f is the sampling frequency of the ultrasonic echo signal, in hertz.

[0029] In step 3, according to the first effective peak time T.sub.x acquired in step 1 and the last effective peak time T.sub.y acquired in step 2, a calculation formula of the thickness L of the tree canopy is obtained as follows:

[00004]$L=k\frac{v_o\left(T_y-T_x\right)}{2}$

[0030] where L is the thickness of the tree canopy, in meters; v.sub.o is a local ultrasonic propagation rate, in meters/second; T.sub.x is the first effective peak time of the ultrasonic echo signal, in seconds; T.sub.y is the last effective peak time of the ultrasonic echo signal, in seconds; k is a correction coefficient. When a leaf area index LAI of the tree canopy meets 1.5≤LAI≤4.0, k=1. When the leaf area index of the tree canopy is

[00005]$\mathrm{LAI}<1.5,k=\frac{.Math.\mathrm{LAI}-1.5.Math.}{3}+1.$

When the leaf area index of the tree canopy is

[00006]$\mathrm{LAI}>4,k=\frac{.Math.\mathrm{LAI}-4.Math.}{2}+1.$

[0031] When the leaf area index LAI of the tree canopy is 3, the correction coefficient k=1, the local ultrasonic propagation rate v.sub.o=340 meters/second, the sampling frequency of the ultrasonic echo signal is f=400000 Hz, the sampling duration k.sub.1=0.001 seconds, and the sampling duration k.sub.2=0.0005 seconds. According to the method provided in the above embodiment, the thickness L of the tree canopy is detected, and the comparison between the detected thickness L of the tree canopy and the actual thickness LO of the tree canopy is shown in Table 1.

TABLE-US-00001 TABLE 1 Comparison between the embodiment of the present invention and the actual thickness of the tree canopy Thickness L of the tree Actual thickness LO of canopy (m) the tree canopy (m) Relative error δ (%) 0.29 0.30 −3.3 0.42 0.40 5.0 0.52 0.50 4.0 0.57 0.60 −5.0

[0032] As shown in Table 1, the value of the thickness L of the tree canopy obtained in the embodiment of the present invention is very close to the value of the actual thickness LO of the tree canopy, with the relative error δ between 3% and 5%. The result shows that the present invention has a characteristic of high accuracy, and is especially suitable for the detection of various fruit trees in mountainous and hilly areas.

[0033] The relative error δ is defined as

[00007]$\delta =100\times \frac{.Math.L-\mathrm{LO}.Math.}{\mathrm{LO}},$

where δ is the relative error, in %. L is the thickness of the tree canopy, in meters, and LO is the actual thickness of the tree canopy, in meters.