Method and device for determining canopy interception based on lidar

Abstract

The present disclosure relates to a method for determining canopy interception based on laser radar, a device for determining canopy interception based on laser radar, a storage medium, and an electronic device. The method for determining canopy interception based on laser radar is applied to the field of ecological environment monitoring, including obtaining the leaf area index, canopy density, leaf inclination angle, and retention coefficient of a target canopy during rainfall based on laser radar; obtaining the splash retention coefficient of a target canopy during rainfall based on a high-speed camera; and determining the canopy interception amount of the target canopy within the rainfall duration based on rainfall intensity, rainfall duration, canopy interception capacity of the target canopy, and interception ratio of the target canopy. In this way, the attachment and collision state of raindrops on canopy leaves during rainfall can be directly measured by laser radar and high-speed camera, thereby determining the canopy interception amount of the target canopy within the rainfall duration.

Description

Technical Field

[0001] The present disclosure relates to the technical field of ecological environment monitoring, and in particular to a method for determining canopy interception based on laser radar, a device for determining canopy interception based on laser radar, a storage medium, and an electronic device. Background Art

[0002] Canopy interception significantly impacts watershed hydrology by affecting rainfall redistribution, reducing raindrop kinetic energy, and influencing runoff and erosion. Furthermore, as a pathway for water input and a means of replenishing nutrients and minerals, it plays a crucial role in ecosystem water balance, nutrient cycling, and vegetation growth. Therefore, the role and impact of canopy interception should be fully considered in watershed management and ecosystem protection. Summary of the Invention

[0003] In view of this, the embodiments of the present disclosure hope to provide a method for determining canopy interception amount based on laser radar, a device for determining canopy interception amount based on laser radar, a storage medium and an electronic device.

[0004] The technical solution of the present disclosure is achieved as follows:

[0005] In a first aspect, the present disclosure provides a method for determining canopy interception based on lidar.

[0006] The method for determining canopy interception based on laser radar provided in the embodiments of the present disclosure includes:

[0007] Obtaining point cloud data of the target canopy based on LiDAR, and performing three-dimensional voxel processing on the point cloud data to obtain the leaf area index LAI, canopy density γ, and leaf inclination angle of the target canopy during rainfall;

[0008] The attachment and collision state of raindrops on canopy leaves during rainfall is captured using a high-speed camera, and the captured images are processed to obtain the splash retention coefficient p of the target canopy during rainfall.

[0009] determining a retention coefficient β of the target canopy based on the leaf inclination angle of the target canopy and a distribution function of the leaf inclination angle of the target canopy;

[0010] Determining the interception ratio K of the target canopy based on the leaf area index LAI of the target canopy, the canopy density of the target canopy, the splash retention coefficient p of the target canopy, and the retention coefficient β of the target canopy;

[0011] Based on the rainfall intensity, the rainfall duration, the canopy interception capacity Y of the target canopy, and the interception ratio K of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined.

[0012] In some embodiments, determining the retention coefficient β of the target canopy based on the leaf inclination angle of the target canopy and the distribution function of the leaf inclination angle of the target canopy includes:

[0013] Based on the leaf inclination angle of the target canopy, the critical leaf inclination angle α at which water droplets attach is determined according to the force balance and the physical properties of raindrops. crit ;

[0014] Based on the critical blade inclination angle α crit and the distribution function f(α) of the leaf inclination angle of the target canopy, and determine the retention coefficient β of the target canopy; wherein,

[0015] Among them, α crit is the critical leaf inclination angle, f(α) is the distribution function of the leaf inclination angle of the target canopy, β is the retention coefficient of the target canopy, and α is the leaf inclination angle of the target canopy.

[0016] In some embodiments, determining the interception ratio K of the target canopy based on the leaf area index LAI of the target canopy, the canopy density of the target canopy, the splash retention coefficient p of the target canopy, and the retention coefficient β of the target canopy includes:

[0017] Among them, K is the interception ratio of the target canopy, γ is the canopy density of the target canopy, β is the retention coefficient of the target canopy, p is the splash retention coefficient of the target canopy, and LAI is the leaf area index of the target canopy.

[0018] In some embodiments, determining the canopy interception amount of the target canopy within the rainfall duration based on rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, and the interception ratio K of the target canopy includes:

[0019] Among them, w is the canopy interception amount of the target canopy during the rainfall duration, Y is the canopy interception capacity of the target canopy, K is the interception ratio of the target canopy, I is the rainfall intensity, and t is the rainfall duration.

[0020] In some embodiments, determining the canopy interception amount of the target canopy within the rainfall duration based on rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, and the interception ratio K of the target canopy includes:

[0021] Obtaining the canopy evaporation capacity of the target canopy and the canopy saturation of the target canopy;

[0022] Based on the rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy and the interception ratio K of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined; wherein,

[0023]

[0024]

[0025] Wherein, w is the canopy interception amount of the target canopy during the rainfall duration, Y is the canopy interception capacity of the target canopy, and D r is the canopy saturation of the target canopy, K is the interception ratio of the target canopy, I is the rainfall intensity, t is the rainfall duration, e p is the canopy evaporation capacity of the target canopy, β is the retention coefficient of the target canopy, p is the splash retention coefficient of the target canopy, and LAI is the leaf area index of the target canopy.

[0026] In some embodiments, determining the canopy interception amount of the target canopy within the rainfall duration based on rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, and the interception ratio K of the target canopy includes:

[0027] Obtaining the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the stem flow ratio of the target canopy, and the stem interception capacity of the target canopy;

[0028] Based on the rainfall intensity, the rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the interception ratio K of the target canopy, the stem flow ratio of the target canopy and the stem interception capacity of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined; wherein,

[0029]

[0030]

[0031] Wherein, w is the canopy interception amount of the target canopy during the rainfall duration, Y is the canopy interception capacity of the target canopy, and D r is the canopy saturation of the target canopy, K is the interception ratio of the target canopy, I is the rainfall intensity, t is the rainfall duration, e p is the canopy evaporation capacity of the target canopy, S is the stem interception capacity of the target canopy, pt is the stem flow ratio, β is the retention coefficient of the target canopy, p is the splash retention coefficient of the target canopy, and LAI is the leaf area index of the target canopy.

[0032] In some embodiments, determining the canopy interception amount of the target canopy within the rainfall duration based on rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, and the interception ratio K of the target canopy includes:

[0033] Obtaining the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the stem flow ratio of the target canopy, the stem interception capacity of the target canopy, and the splash drop evaporation of the target canopy;

[0034] Based on the rainfall intensity, the rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the interception ratio K of the target canopy, the stem flow ratio of the target canopy, the stem interception capacity of the target canopy and the splash drop evaporation of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined; wherein,

[0035]

[0036]

[0037]

[0038] E f =I·K ef t;

[0039] in,

[0040] w is the canopy interception amount of the target canopy during the rainfall duration, Y is the canopy interception capacity of the target canopy, D r is the canopy saturation of the target canopy, K is the interception ratio of the target canopy, I is the rainfall intensity, t is the rainfall duration, e p is the canopy evaporation capacity of the target canopy, S is the stem interception capacity of the target canopy, and p t is the stem flow ratio, β is the retention coefficient of the target canopy, p is the splash retention coefficient of the target canopy, LAI is the leaf area index of the target canopy, K ef is the splash drop evaporation amount of the target canopy, and m is the splash drop evaporation ratio.

[0041] In a second aspect, the present disclosure provides a canopy interception determination device based on laser radar, comprising:

[0042] The first information acquisition module is used to obtain point cloud data of the target canopy based on the laser radar, and perform three-dimensional voxel processing on the point cloud data to obtain the leaf area index LAI of the target canopy during rainfall, the canopy density γ of the target canopy, and the leaf inclination angle of the target canopy;

[0043] The second information acquisition module is used to capture the attachment and collision state of raindrops on canopy leaves during rainfall using a high-speed camera, and to process the captured images to obtain a splash retention coefficient p of the target canopy during rainfall;

[0044] A first information determination module is configured to determine a retention coefficient β of a target canopy based on the leaf inclination angle of the target canopy and a distribution function of the leaf inclination angle of the target canopy;

[0045] A second information determination module is configured to determine a retention ratio K of the target canopy based on a leaf area index LAI of the target canopy, a canopy density of the target canopy, a splash retention coefficient p of the target canopy, and a retention coefficient β of the target canopy;

[0046] The canopy interception amount determination module is used to determine the canopy interception amount of the target canopy within the rainfall duration based on the rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy and the interception ratio K of the target canopy.

[0047] In a third aspect, the present disclosure provides a computer-readable storage medium on which a lidar-based canopy interception determination program is stored. When the lidar-based canopy interception determination program is executed by a processor, the lidar-based canopy interception determination method described in the first aspect above is implemented.

[0048] In a fourth aspect, the present disclosure provides an electronic device comprising a memory, a processor, and a lidar-based canopy interception determination program stored in the memory and executable on the processor. When the processor executes the lidar-based canopy interception determination program, the lidar-based canopy interception determination method described in the first aspect above is implemented.

[0049] According to the lidar-based canopy interception amount determination method of the embodiment of the present disclosure, the method includes acquiring point cloud data of the target canopy based on the lidar, and performing three-dimensional voxel processing on the point cloud data to obtain the leaf area index, canopy density and leaf inclination angle of the target canopy during rainfall; photographing the attachment and collision state of raindrops on canopy leaves during rainfall based on a high-speed camera, and performing image processing on the photographed image to obtain the splash retention coefficient of the target canopy during rainfall; determining the retention coefficient of the target canopy based on the leaf inclination angle of the target canopy and the distribution function of the leaf inclination angle of the target canopy; determining the interception ratio of the target canopy based on the leaf area index of the target canopy, the canopy density of the target canopy, the splash retention coefficient of the target canopy and the retention coefficient of the target canopy; and determining the canopy interception amount of the target canopy within the rainfall duration based on the rainfall intensity, rainfall duration, the canopy interception capacity of the target canopy and the interception ratio of the target canopy. In the present application, point cloud data of the target canopy is obtained by using a laser radar, and the point cloud data is subjected to three-dimensional voxel processing to obtain the leaf area index, canopy density and leaf inclination angle of the target canopy during rainfall. The attachment and collision state of raindrops on the canopy leaves during rainfall are photographed by a high-speed camera, and the photographed images are processed to obtain the splash retention coefficient of the target canopy during rainfall. On this basis, the interception ratio of the target canopy is determined based on the leaf area index of the target canopy, the canopy density of the target canopy, the splash retention coefficient of the target canopy and the retention coefficient of the target canopy. The canopy interception amount of the target canopy within the rainfall duration is determined based on the rainfall intensity, rainfall duration, the canopy interception capacity of the target canopy and the interception ratio of the target canopy. In this way, the attachment and collision state of raindrops on the canopy leaves during rainfall can be directly measured by a laser radar and a high-speed camera, so as to determine the canopy interception amount of the target canopy within the rainfall duration.

[0050] Additional aspects and advantages of the present disclosure will be given in part in the following description and in part will be obvious from the following description, or will be learned through practice of the present disclosure. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 is a flow chart of a method for determining canopy interception based on lidar according to an exemplary embodiment;

[0052] Figure 2 It is a schematic structural diagram of a device for determining canopy interception based on lidar according to an exemplary embodiment. DETAILED DESCRIPTION

[0053] The following describes in detail embodiments of the present disclosure, examples of which are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present disclosure, and should not be construed as limiting the present disclosure.

[0054] Canopy interception significantly impacts watershed hydrology by affecting rainfall redistribution, reducing raindrop kinetic energy, and influencing runoff and erosion. Furthermore, as a pathway for water input and a means of replenishing nutrients and minerals, it plays a crucial role in ecosystem water balance, nutrient cycling, and vegetation growth. Therefore, the role and impact of canopy interception should be fully considered in watershed management and ecosystem protection.

[0055] In response to the above situation, the present disclosure provides a method for determining canopy interception based on lidar. Figure 1 This is a process flow of a method for determining canopy interception based on lidar according to an exemplary embodiment. Figure 1 .like Figure 1 As shown, the lidar-based canopy interception determination method includes:

[0056] Step 10: acquiring point cloud data of the target canopy based on the laser radar, and performing three-dimensional voxel processing on the point cloud data to obtain the leaf area index LAI of the target canopy during rainfall, the canopy density γ of the target canopy, and the leaf inclination angle of the target canopy;

[0057] Step 11: Using a high-speed camera, capture the attachment and collision states of raindrops on canopy leaves during rainfall, and perform image processing on the captured images to obtain a splash retention coefficient p of the target canopy during rainfall;

[0058] Step 12: determining a retention coefficient β of the target canopy based on the leaf inclination angle of the target canopy and a distribution function of the leaf inclination angle of the target canopy;

[0059] Step 13: determining the interception ratio K of the target canopy based on the leaf area index LAI of the target canopy, the canopy density of the target canopy, the splash retention coefficient p of the target canopy, and the retention coefficient β of the target canopy;

[0060] Step 14: Determine the canopy interception amount of the target canopy within the rainfall duration based on the rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, and the interception ratio K of the target canopy.

[0061] In this exemplary embodiment, the laser radar may be a ground-based laser radar, a handheld mobile laser radar, an unmanned aerial vehicle (UAV)-mounted laser radar, or the like. In specific applications, any one of these radar devices may be used, or multiple radar devices may be used in combination, such as a ground-based laser radar combined with an unmanned aerial vehicle (UAV) laser radar to acquire point cloud data of the target canopy.

[0062] In this exemplary embodiment, the leaf area index LAI of the target canopy is:

[0063] LAI = -lnP(θ)cosθ / G(θ)·, where P(θ) is the clearance ratio at the zenith angle θ, G(θ) is the average projected area ratio of the leaves at the zenith angle θ, and Ω is the clustering coefficient. The calculation of θ assumes that the rain falls vertically and is set to 0.

[0064] In this exemplary embodiment, a high-speed camera is used to observe the adhesion and collision of water droplets on canopy leaf surfaces during rainfall. The canopy's retention capacity is determined based on the amount and shape of the attached droplets. Based on the collision and splashing of the droplets, the splash retention coefficient is determined—the proportion of water retained on the canopy structure during splashing. Combined with the leaf inclination distribution measured in the first step, the critical leaf inclination angle at which a droplet can adhere is calculated based on force balance and the physical properties of the raindrops, thereby determining the retention coefficient.

[0065] In this exemplary embodiment, vegetation coverage (canopy density, FVC) represents the percentage of canopy projected area within a given area. Assuming uniform rainfall across the area, it can represent the proportion of raindrops that collide with the canopy (non-penetrating rainfall) or the proportion of penetrating rainfall (i.e., penetration rate). Penetrating rainfall is defined as rainfall that is not intercepted.

[0066] In the present application, point cloud data of the target canopy is obtained by using a laser radar, and the point cloud data is subjected to three-dimensional voxel processing to obtain the leaf area index, canopy density and leaf inclination angle of the target canopy during rainfall. The attachment and collision state of raindrops on the canopy leaves during rainfall are photographed by a high-speed camera, and the photographed images are processed to obtain the splash retention coefficient of the target canopy during rainfall. On this basis, the interception ratio of the target canopy is determined based on the leaf area index of the target canopy, the canopy density of the target canopy, the splash retention coefficient of the target canopy and the retention coefficient of the target canopy. The canopy interception amount of the target canopy within the rainfall duration is determined based on the rainfall intensity, rainfall duration, the canopy interception capacity of the target canopy and the interception ratio of the target canopy. In this way, the attachment and collision state of raindrops on the canopy leaves during rainfall can be directly measured by a laser radar and a high-speed camera, so as to determine the canopy interception amount of the target canopy within the rainfall duration.

[0067] In some embodiments, determining the retention coefficient β of the target canopy based on the leaf inclination angle of the target canopy and the distribution function of the leaf inclination angle of the target canopy includes:

[0068] Based on the leaf inclination angle of the target canopy, the critical leaf inclination angle α at which water droplets attach is determined according to the force balance and the physical properties of raindrops. crir ;

[0069] Based on the critical blade inclination angle α crit and the distribution function f(α) of the leaf inclination angle of the target canopy, and determine the retention coefficient β of the target canopy; wherein,

[0070] Among them, α crit is the critical leaf inclination angle, f(α) is the distribution function of the leaf inclination angle of the target canopy, β is the retention coefficient of the target canopy, and α is the leaf inclination angle of the target canopy.

[0071] In this exemplary embodiment,

[0072]

[0073] Among them, θ a is the advancing contact angle, θ r is the receding contact angle, V is the droplet volume, and σ is the surface tension coefficient of water σ=72×10 -3 N / m, g is the acceleration due to gravity, and ρ is the density of water.

[0074] In this exemplary embodiment, the physical meaning of the attachment and retention coefficient β is that, for droplets that temporarily adhere to a leaf surface after splashing, a β-percentage of droplets remain permanently on the leaf surface, while a 1-β-percentage of droplets will immediately slide off after temporary attachment due to an inclination angle greater than the critical angle. α is the canopy leaf inclination angle, and f(α) is the distribution function of the canopy leaf inclination angle, which can be obtained by actual measurement or fitted using the beta function.

[0075] In this exemplary embodiment, splashing drops 1-p formed by raindrop collisions and droplets 1-β sliding off tilted leaves continue to fall, repeating the same collision process. Assuming the canopy is composed of multiple layers of leaves that are completely covered and non-overlapping within a range of canopy density, and ignoring lateral displacement after raindrop collisions—that is, raindrops do not move beyond the canopy density after collision—the average number of raindrop collisions is the approximate number of canopy layers, expressed as LAI / γ. The majority of splashing and sliding raindrops repeat the same collision, splashing, and other motion processes.

[0076] In some embodiments, determining the interception ratio K of the target canopy based on the leaf area index LAI of the target canopy, the canopy density of the target canopy, the splash retention coefficient p of the target canopy, and the retention coefficient β of the target canopy includes:

[0077] Among them, K is the interception ratio of the target canopy, γ is the canopy density of the target canopy, β is the retention coefficient of the target canopy, p is the splash retention coefficient of the target canopy, and LAI is the leaf area index of the target canopy.

[0078] In this exemplary embodiment, determining the canopy interception amount of the target canopy within the rainfall duration based on the rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, and the interception ratio K of the target canopy includes:

[0079] Among them, w is the canopy interception amount of the target canopy during the rainfall duration, Y is the canopy interception capacity of the target canopy, K is the interception ratio of the target canopy, I is the rainfall intensity, and t is the rainfall duration.

[0080] In this exemplary embodiment, when determining the canopy interception amount of the target canopy, when considering evaporation during the rainy period, the canopy interception amount of the target canopy within the rainfall duration is determined based on the rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, and the interception ratio K of the target canopy, including:

[0081] Obtaining the canopy evaporation capacity of the target canopy and the canopy saturation of the target canopy;

[0082] Based on the rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy and the interception ratio K of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined; wherein,

[0083]

[0084]

[0085] Wherein, w is the canopy interception amount of the target canopy during the rainfall duration, Y is the canopy interception capacity of the target canopy, and D r is the canopy saturation of the target canopy, K is the interception ratio of the target canopy, I is the rainfall intensity, t is the rainfall duration, e p is the canopy evaporation capacity of the target canopy, β is the retention coefficient of the target canopy, p is the splash retention coefficient of the target canopy, and LAI is the leaf area index of the target canopy.

[0086] In this exemplary embodiment, when determining the canopy interception capacity of a target canopy, the canopy evaporation capacity and canopy saturation of the target canopy need to be considered when considering evaporation during a rainy period. The canopy interception capacity of the target canopy during the rainy period is determined based on the rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, and the interception ratio K of the target canopy.

[0087] In some embodiments, when determining the canopy interception amount of the target canopy, when considering evaporation and stem flow during the rainy period, the method of determining the canopy interception amount of the target canopy within the rainfall duration based on rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, and the interception ratio K of the target canopy includes:

[0088] Obtaining the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the stem flow ratio of the target canopy, and the stem interception capacity of the target canopy;

[0089] Based on the rainfall intensity, the rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the interception ratio K of the target canopy, the stem flow ratio of the target canopy and the stem interception capacity of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined; wherein,

[0090]

[0091]

[0092] Wherein, w is the canopy interception amount of the target canopy during the rainfall duration, Y is the canopy interception capacity of the target canopy, and D r is the canopy saturation of the target canopy, K is the interception ratio of the target canopy, I is the rainfall intensity, t is the rainfall duration, e p is the canopy evaporation capacity of the target canopy, S is the stem interception capacity of the target canopy, and p t is the stemflow ratio, β is the retention coefficient of the target canopy, p is the splash retention coefficient of the target canopy, and LAI is the leaf area index of the target canopy.

[0093] In this exemplary embodiment, p t is the stem flow ratio,

[0094] Among them, SAI is the stem area index, which can be calculated from the point cloud data obtained by lidar.

[0095] In this exemplary embodiment, when determining the canopy interception amount of the target canopy, when considering evaporation and stem flow during the rainy season, it is necessary to obtain the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the stem flow ratio of the target canopy, and the stem interception capacity of the target canopy;

[0096] Based on the rainfall intensity, the rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the interception ratio K of the target canopy, the stem flow ratio of the target canopy and the stem interception capacity of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined.

[0097] In some embodiments, when determining the canopy interception amount of the target canopy, when considering rainy period evaporation, stem flow, and splash drop evaporation, the canopy interception amount of the target canopy within the rainfall duration is determined based on the rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, and the interception ratio K of the target canopy, including:

[0098] Obtaining the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the stem flow ratio of the target canopy, the stem interception capacity of the target canopy, and the splash drop evaporation of the target canopy;

[0099] Based on the rainfall intensity, the rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the interception ratio K of the target canopy, the stem flow ratio of the target canopy, the stem interception capacity of the target canopy and the splash drop evaporation of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined; wherein,

[0100]

[0101]

[0102]

[0103] E f =I·K ef t;

[0104] Wherein, w is the canopy interception amount of the target canopy during the rainfall duration, Y is the canopy interception capacity of the target canopy, and D r is the canopy saturation of the target canopy, K is the interception ratio of the target canopy, I is the rainfall intensity, t is the rainfall duration, e p is the canopy evaporation capacity of the target canopy, S is the stem interception capacity of the target canopy, and p tis the stem flow ratio, β is the retention coefficient of the target canopy, p is the splash retention coefficient of the target canopy, LAI is the leaf area index of the target canopy, K ef is the splash drop evaporation amount of the target canopy, and m is the splash drop evaporation ratio.

[0105] In this exemplary embodiment, when determining the canopy interception capacity of the target canopy, when considering rainy season evaporation, stem flow, and splash drop evaporation, it is necessary to obtain the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the stem flow ratio of the target canopy, the stem interception capacity of the target canopy, and the splash drop evaporation of the target canopy;

[0106] Based on the rainfall intensity, the rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the interception ratio K of the target canopy, the stem flow ratio of the target canopy, the stem interception capacity of the target canopy and the splash drop evaporation of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined.

[0107] In this application, the layered canopy is formed by superimposing the overall form in a layered form. The impact of the previous rainfall is based on the existing interception of the canopy. The final interception minus the existing interception is the interception during this rainfall period.

[0108] The present disclosure provides a canopy interception determination device based on laser radar. Figure 2 FIG is a schematic diagram of a structure of a device for determining canopy interception based on a laser radar according to an exemplary embodiment. Figure 2 As shown, a canopy interception determination device based on laser radar includes:

[0109] The first information acquisition module 20 is used to acquire point cloud data of the target canopy based on the laser radar, and perform three-dimensional voxel processing on the point cloud data to obtain the leaf area index LAI, the canopy density γ and the leaf inclination angle of the target canopy during rainfall;

[0110] The second information acquisition module 21 is used to capture the attachment and collision state of raindrops on canopy leaves during rainfall using a high-speed camera, and to process the captured images to obtain a splash retention coefficient p of the target canopy during rainfall;

[0111] A first information determination module 22 is configured to determine a retention coefficient β of a target canopy based on the leaf inclination angle of the target canopy and a distribution function of the leaf inclination angle of the target canopy;

[0112] A second information determination module 23 is configured to determine a retention ratio K of the target canopy based on a leaf area index LAI of the target canopy, a canopy density of the target canopy, a splash retention coefficient p of the target canopy, and a retention coefficient β of the target canopy;

[0113] The canopy interception amount determination module 24 is used to determine the canopy interception amount of the target canopy within the rainfall duration based on the rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy and the interception ratio K of the target canopy.

[0114] In this exemplary embodiment, the laser radar may be a ground-based laser radar, a handheld mobile laser radar, an unmanned aerial vehicle (UAV)-mounted laser radar, or the like. In specific applications, any one of these radar devices may be used, or multiple radar devices may be used in combination, such as a ground-based laser radar combined with an unmanned aerial vehicle (UAV) laser radar to acquire point cloud data of the target canopy.

[0115] In this exemplary embodiment, the leaf area index LAI of the target canopy is:

[0116] LAI = -lnP(θ)cosθ / G(θ)·, where P(θ) is the clearance ratio at the zenith angle θ, G(θ) is the average projected area ratio of the leaves at the zenith angle θ, and Ω is the clustering coefficient. The calculation of θ assumes that the rain falls vertically and is set to 0.

[0117] In this exemplary embodiment, a high-speed camera is used to observe the adhesion and collision of water droplets on canopy leaf surfaces during rainfall. The canopy's retention capacity is determined based on the amount and shape of the attached droplets. Based on the collision and splashing of the droplets, the splash retention coefficient is determined—the proportion of water retained on the canopy structure during splashing. Combined with the leaf inclination distribution measured in the first step, the critical leaf inclination angle at which a droplet can adhere is calculated based on force balance and the physical properties of the raindrops, thereby determining the retention coefficient.

[0118] In this exemplary embodiment, vegetation coverage (canopy density, FVC) represents the percentage of canopy projected area within a given area. Assuming uniform rainfall across the area, it can represent the proportion of raindrops that collide with the canopy (non-penetrating rainfall) or the proportion of penetrating rainfall (i.e., penetration rate). Penetrating rainfall is defined as rainfall that is not intercepted.

[0119] In the present application, point cloud data of the target canopy is obtained by using a laser radar, and the point cloud data is subjected to three-dimensional voxel processing to obtain the leaf area index, canopy density and leaf inclination angle of the target canopy during rainfall. The attachment and collision state of raindrops on the canopy leaves during rainfall are photographed by a high-speed camera, and the photographed images are processed to obtain the splash retention coefficient of the target canopy during rainfall. On this basis, the interception ratio of the target canopy is determined based on the leaf area index of the target canopy, the canopy density of the target canopy, the splash retention coefficient of the target canopy and the retention coefficient of the target canopy. The canopy interception amount of the target canopy within the rainfall duration is determined based on the rainfall intensity, rainfall duration, the canopy interception capacity of the target canopy and the interception ratio of the target canopy. In this way, the attachment and collision state of raindrops on the canopy leaves during rainfall can be directly measured by a laser radar and a high-speed camera, so as to determine the canopy interception amount of the target canopy within the rainfall duration.

[0120] In some embodiments, the first information determination module 22 is used to

[0121] Based on the leaf inclination angle of the target canopy, the critical leaf inclination angle α at which water droplets attach is determined according to the force balance and the physical properties of raindrops. crit ;

[0122] Based on the critical blade inclination angle α crit and the distribution function f(α) of the leaf inclination angle of the target canopy, and determine the retention coefficient β of the target canopy; wherein,

[0123] Among them, α crit is the critical leaf inclination angle, f(α) is the distribution function of the leaf inclination angle of the target canopy, β is the retention coefficient of the target canopy, and α is the leaf inclination angle of the target canopy.

[0124] In this exemplary embodiment,

[0125]

[0126] Among them, θ a is the advancing contact angle, θ r is the receding contact angle, V is the droplet volume, and σ is the surface tension coefficient of water σ=72×10 -3 N / m, g is the acceleration due to gravity, and ρ is the density of water.

[0127] In this exemplary embodiment, the physical meaning of the attachment and retention coefficient β is that, for droplets that temporarily adhere to a leaf surface after splashing, a β-percentage of droplets remain permanently on the leaf surface, while a 1-β-percentage of droplets will immediately slide off after temporary attachment due to an inclination angle greater than the critical angle. α is the canopy leaf inclination angle, and f(α) is the distribution function of the canopy leaf inclination angle, which can be obtained by actual measurement or fitted using the beta function.

[0128] In this exemplary embodiment, splashing drops 1-p formed by raindrop collisions and droplets 1-β sliding off tilted leaves continue to fall, repeating the same collision process. Assuming the canopy is composed of multiple layers of leaves that are completely covered and non-overlapping within a range of canopy density, and ignoring lateral displacement after raindrop collisions—that is, raindrops do not move beyond the canopy density after collision—the average number of raindrop collisions is the approximate number of canopy layers, expressed as LAI / γ. The majority of splashing and sliding raindrops repeat the same collision, splashing, and other motion processes.

[0129] In some embodiments, determining the interception ratio K of the target canopy based on the leaf area index LAI of the target canopy, the canopy density of the target canopy, the splash retention coefficient p of the target canopy, and the retention coefficient β of the target canopy includes:

[0130] Among them, K is the interception ratio of the target canopy, γ is the canopy density of the target canopy, β is the retention coefficient of the target canopy, p is the splash retention coefficient of the target canopy, and LAI is the leaf area index of the target canopy.

[0131] In this exemplary embodiment, determining the canopy interception amount of the target canopy within the rainfall duration based on the rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, and the interception ratio K of the target canopy includes:

[0132] Among them, w is the canopy interception amount of the target canopy during the rainfall duration, Y is the canopy interception capacity of the target canopy, K is the interception ratio of the target canopy, I is the rainfall intensity, and t is the rainfall duration.

[0133] In this exemplary embodiment, the canopy interception determination module 24 is used to

[0134] Obtaining the canopy evaporation capacity of the target canopy and the canopy saturation of the target canopy;

[0135] Based on the rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy and the interception ratio K of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined; wherein,

[0136]

[0137]

[0138] Wherein, w is the canopy interception amount of the target canopy during the rainfall duration, Y is the canopy interception capacity of the target canopy, and D r is the canopy saturation of the target canopy, K is the interception ratio of the target canopy, I is the rainfall intensity, t is the rainfall duration, e p is the canopy evaporation capacity of the target canopy, β is the retention coefficient of the target canopy, p is the splash retention coefficient of the target canopy, and LAI is the leaf area index of the target canopy.

[0139] In this exemplary embodiment, when determining the canopy interception capacity of a target canopy, the canopy evaporation capacity and canopy saturation of the target canopy need to be considered when considering evaporation during a rainy period. The canopy interception capacity of the target canopy during the rainy period is determined based on the rainfall intensity, rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, and the interception ratio K of the target canopy.

[0140] In some embodiments, when determining the canopy interception of the target canopy, when considering rainy season evaporation and stem flow, the canopy interception determination module 24 is used to

[0141] Obtaining the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the stem flow ratio of the target canopy, and the stem interception capacity of the target canopy;

[0142] Based on the rainfall intensity, the rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the interception ratio K of the target canopy, the stem flow ratio of the target canopy and the stem interception capacity of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined; wherein,

[0143]

[0144]

[0145] Wherein, w is the canopy interception amount of the target canopy during the rainfall duration, Y is the canopy interception capacity of the target canopy, and D r is the canopy saturation of the target canopy, K is the interception ratio of the target canopy, I is the rainfall intensity, t is the rainfall duration, e pis the canopy evaporation capacity of the target canopy, S is the stem interception capacity of the target canopy, and p t is the stemflow ratio, β is the retention coefficient of the target canopy, p is the splash retention coefficient of the target canopy, and LAI is the leaf area index of the target canopy.

[0146] In this exemplary embodiment, p t is the stem flow ratio,

[0147] Among them, SAI is the stem area index, which can be calculated from the point cloud data obtained by lidar.

[0148] In this exemplary embodiment, when determining the canopy interception amount of the target canopy, when considering evaporation and stem flow during the rainy season, it is necessary to obtain the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the stem flow ratio of the target canopy, and the stem interception capacity of the target canopy;

[0149] Based on the rainfall intensity, the rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the interception ratio K of the target canopy, the stem flow ratio of the target canopy and the stem interception capacity of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined.

[0150] In some embodiments, when determining the canopy interception of the target canopy, when considering rainy season evaporation, stem flow and splash drop evaporation, the canopy interception determination module 24 is used to

[0151] Obtaining the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the stem flow ratio of the target canopy, the stem interception capacity of the target canopy, and the splash drop evaporation of the target canopy;

[0152] Based on the rainfall intensity, the rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the interception ratio K of the target canopy, the stem flow ratio of the target canopy, the stem interception capacity of the target canopy and the splash drop evaporation of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined; wherein,

[0153]

[0154]

[0155]

[0156] E f =I·K ef t;

[0157] in,

[0158] w is the canopy interception amount of the target canopy during the rainfall duration, Y is the canopy interception capacity of the target canopy, D r is the canopy saturation of the target canopy, K is the interception ratio of the target canopy, I is the rainfall intensity, t is the rainfall duration, e p is the canopy evaporation capacity of the target canopy, S is the stem interception capacity of the target canopy, and p t is the stem flow ratio, β is the retention coefficient of the target canopy, p is the splash retention coefficient of the target canopy, LAI is the leaf area index of the target canopy, K ef is the splash drop evaporation amount of the target canopy, and m is the splash drop evaporation ratio.

[0159] In this exemplary embodiment, when determining the canopy interception capacity of the target canopy, when considering rainy season evaporation, stem flow, and splash drop evaporation, it is necessary to obtain the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the stem flow ratio of the target canopy, the stem interception capacity of the target canopy, and the splash drop evaporation of the target canopy;

[0160] Based on the rainfall intensity, the rainfall duration, the canopy interception capacity Y of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the interception ratio K of the target canopy, the stem flow ratio of the target canopy, the stem interception capacity of the target canopy and the splash drop evaporation of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined.

[0161] In this application, the layered canopy is formed by superimposing the overall form in a layered form. The impact of the previous rainfall is based on the existing interception of the canopy. The final interception minus the existing interception is the interception during this rainfall period.

[0162] The present disclosure provides a computer-readable storage medium on which a lidar-based canopy interception determination program is stored. When the lidar-based canopy interception determination program is executed by a processor, the lidar-based canopy interception determination method described in the above embodiments is implemented.

[0163] The present disclosure provides an electronic device, including a memory, a processor, and a lidar-based canopy interception determination program stored in the memory and runnable on the processor. When the processor executes the lidar-based canopy interception determination program, the lidar-based canopy interception determination method described in the above-mentioned embodiments is implemented.

[0164] It should be noted that the logic and / or steps represented in the flowcharts or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing the logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (e.g., a computer-based system, a system including a processor, or other system that can retrieve and execute instructions from an instruction execution system, apparatus, or device). For purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transmit a program for use by, or in conjunction with, an instruction execution system, apparatus, or device. More specific examples (non-exhaustive list) of computer-readable media include the following: an electrical connection portion having one or more wires (electronic device), a portable computer disk cartridge (magnetic device), random access memory (RAM), read-only memory (ROM), erasable and programmable read-only memory (EPROM or flash memory), fiber optic devices, and portable compact disc read-only memory (CDROM). Furthermore, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium and then editing, interpreting or otherwise processing it in a suitable manner if necessary, and then storing it in a computer memory.

[0165] It should be understood that various parts of the present disclosure can be implemented using hardware, software, firmware, or a combination thereof. In the above embodiments, multiple steps or methods can be implemented using software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented using hardware, as in another embodiment, any one of the following technologies known in the art or a combination thereof can be used to implement: a discrete logic circuit having a logic gate circuit for implementing a logic function on a data signal, an application-specific integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0166] Throughout this specification, reference to terms such as "one embodiment," "some embodiments," "examples," "specific examples," or "some examples" means that a specific feature, structure, material, or characteristic described in conjunction with that embodiment or example is included in at least one embodiment or example of the present disclosure. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

[0167] In the description of the present disclosure, it should be understood that the terms "center", "longitudinal", "lateral", "length", "width", "thickness", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside", "clockwise", "counterclockwise", "axial", "radial", "circumferential" and the like indicate orientations or positional relationships based on the orientations or positional relationships shown in the accompanying drawings, and are only for the convenience of describing the present disclosure and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be understood as a limitation to the present disclosure.

[0168] In addition, the terms "first" and "second" used in the embodiments of the present disclosure are only used for descriptive purposes and should not be understood as indicating or implying relative importance, or implicitly indicating the number of technical features indicated in the embodiments. Therefore, the features defined in the embodiments of the present disclosure with terms such as "first" and "second" can explicitly or implicitly indicate that the embodiment includes at least one such feature. In the description of the present disclosure, the word "plurality" means at least two or two or more, such as two, three, four, etc., unless otherwise clearly and specifically defined in the embodiments.

[0169] In this disclosure, unless otherwise clearly specified or limited in the embodiments, the terms "installed," "connected," "connect," and "fixed" appearing in the embodiments should be understood in a broad sense. For example, the connection can be a fixed connection, a detachable connection, or an integral connection. It can also be a mechanical connection, an electrical connection, etc.; of course, it can also be a direct connection, or an indirect connection through an intermediate medium, or it can be the internal communication between two elements, or the interaction between two elements. For those skilled in the art, the specific meanings of the above terms in this disclosure can be understood based on the specific implementation.

[0170] In the present disclosure, unless otherwise expressly specified or limited, when a first feature is "above" or "below" a second feature, it may mean that the first and second features are in direct contact, or the first and second features are in indirect contact through an intermediate medium. Furthermore, when a first feature is "above," "above," or "above" a second feature, it may mean that the first feature is directly above or diagonally above the second feature, or simply means that the first feature is at a higher level than the second feature. When a first feature is "below," "below," or "below" a second feature, it may mean that the first feature is directly below or diagonally below the second feature, or simply means that the first feature is at a lower level than the second feature.

[0171] Although the embodiments of the present disclosure have been shown and described above, it is understood that the above embodiments are illustrative and are not to be construed as limitations on the present disclosure. A person skilled in the art may change, modify, replace and vary the above embodiments within the scope of the present disclosure.

Claims

1. A method for determining canopy interception based on laser radar, characterized in that: include: Obtaining point cloud data of the target canopy based on a laser radar, and performing three-dimensional voxel processing on the point cloud data to obtain the leaf area index, canopy density, and leaf inclination angle of the target canopy during rainfall; The attachment and collision state of raindrops on canopy leaves during rainfall is captured using a high-speed camera, and the captured images are processed to obtain the splash retention coefficient of the target canopy during rainfall. determining a retention coefficient of the target canopy based on the leaf inclination angle of the target canopy and a distribution function of the leaf inclination angle of the target canopy; Determining the interception ratio of the target canopy based on the leaf area index of the target canopy, the canopy density of the target canopy, the splash retention coefficient of the target canopy, and the retention coefficient of the target canopy; The canopy interception amount of the target canopy within the rainfall duration is determined based on the rainfall intensity, the rainfall duration, the canopy interception capacity of the target canopy and the interception ratio of the target canopy.

2. The method for determining canopy interception based on laser radar according to claim 1, characterized in that: The determining of the retention coefficient of the target canopy based on the leaf inclination angle of the target canopy and the distribution function of the leaf inclination angle of the target canopy includes: Based on the leaf inclination angle of the target canopy, determining the critical leaf inclination angle at which water droplets attach according to force balance and physical properties of raindrops; Based on the critical blade inclination angle and the distribution function of the blade inclination angle of the target canopy, the retention coefficient of the target canopy is determined; wherein, ;in, is the critical inclination angle of the blade, is the distribution function of leaf inclination angle of target canopy, is the retention coefficient of the target canopy, is the leaf inclination angle of the target canopy.

3. The method for determining canopy interception based on laser radar according to claim 1, characterized in that: The determining of the interception ratio of the target canopy based on the leaf area index of the target canopy, the canopy density of the target canopy, the splash retention coefficient of the target canopy, and the retention coefficient of the target canopy includes: ;in, is the interception ratio of the target canopy, is the target canopy density, is the retention coefficient of the target canopy, is the splash retention coefficient of the target canopy, is the leaf area index of the target canopy.

4. The method for determining canopy interception based on laser radar according to claim 1, characterized in that: The determining of the canopy interception amount of the target canopy within the rainfall duration based on the rainfall intensity, the rainfall duration, the canopy interception capacity of the target canopy, and the interception ratio of the target canopy includes: ;in, is the canopy interception amount of the target canopy during the rainfall duration, is the canopy interception capacity of the target canopy, is the rainfall intensity, is the interception ratio of the target canopy, The duration of rainfall.

5. The method for determining canopy interception based on laser radar according to claim 2, characterized in that: The determining of the canopy interception amount of the target canopy within the rainfall duration based on the rainfall intensity, the rainfall duration, the canopy interception capacity of the target canopy, and the interception ratio of the target canopy includes: Obtaining the canopy evaporation capacity of the target canopy and the canopy saturation of the target canopy; Based on the rainfall intensity, rainfall duration, the canopy interception capacity of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy and the interception ratio of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined; wherein, ； ； ;in, is the canopy interception amount of the target canopy during the rainfall duration, is the canopy interception capacity of the target canopy, is the target canopy density, is the canopy saturation of the target canopy, is the interception ratio of the target canopy, is the rainfall intensity, For the duration of rainfall, is the canopy evaporation capacity of the target canopy, is the retention coefficient of the target canopy, is the splash retention coefficient of the target canopy, is the leaf area index of the target canopy.

6. The method for determining canopy interception based on laser radar according to claim 2, characterized in that: The determining of the canopy interception amount of the target canopy within the rainfall duration based on the rainfall intensity, the rainfall duration, the canopy interception capacity of the target canopy, and the interception ratio of the target canopy includes: Obtaining the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the stem flow ratio of the target canopy, and the stem interception capacity of the target canopy; Based on the rainfall intensity, the rainfall duration, the canopy interception capacity of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the interception ratio of the target canopy, the stemflow ratio of the target canopy and the stem interception capacity of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined; wherein, ； ； ;in, is the canopy interception amount of the target canopy during the rainfall duration, is the canopy interception capacity of the target canopy, is the target canopy density, is the canopy saturation of the target canopy, is the interception ratio of the target canopy, is the rainfall intensity, For the duration of rainfall, is the canopy evaporation capacity of the target canopy, is the stem interception capacity of the target canopy, is the stem flow ratio, is the retention coefficient of the target canopy, is the splash retention coefficient of the target canopy, is the leaf area index of the target canopy.

7. The method for determining canopy interception based on laser radar according to claim 2, characterized in that: The determining of the canopy interception amount of the target canopy within the rainfall duration based on the rainfall intensity, the rainfall duration, the canopy interception capacity of the target canopy, and the interception ratio of the target canopy includes: Obtaining the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the stem flow ratio of the target canopy, the stem interception capacity of the target canopy, and the splash drop evaporation of the target canopy; Based on the rainfall intensity, the rainfall duration, the canopy interception capacity of the target canopy, the canopy evaporation capacity of the target canopy, the canopy saturation of the target canopy, the interception ratio of the target canopy, the stem flow ratio of the target canopy, the stem interception capacity of the target canopy and the splash drop evaporation of the target canopy, the canopy interception amount of the target canopy within the rainfall duration is determined; wherein, ； ； ； ； ;in, is the canopy interception amount of the target canopy during the rainfall duration, is the canopy interception capacity of the target canopy, is the target canopy density, is the canopy saturation of the target canopy, is the interception ratio of the target canopy, is the rainfall intensity, For the duration of rainfall, is the canopy evaporation capacity of the target canopy, is the stem interception capacity of the target canopy, is the stem flow ratio, is the retention coefficient of the target canopy, is the splash retention coefficient of the target canopy, is the leaf area index of the target canopy, is the splash drop evaporation coefficient of the target canopy, is the splash drop evaporation ratio.

8. A device for determining canopy interception based on laser radar, characterized in that: include: The first information acquisition module is used to obtain point cloud data of the target canopy based on the laser radar, and perform three-dimensional voxel processing on the point cloud data to obtain the leaf area index, canopy density and leaf inclination angle of the target canopy during rainfall; The second information acquisition module is used to capture the attachment and collision state of raindrops on canopy leaves during rainfall using a high-speed camera, and to process the captured images to obtain the splash retention coefficient of the target canopy during rainfall; a first information determination module, configured to determine a retention coefficient of a target canopy based on the leaf inclination angle of the target canopy and a distribution function of the leaf inclination angle of the target canopy; A second information determination module is configured to determine the interception ratio of the target canopy based on the leaf area index of the target canopy, the canopy density of the target canopy, the splash retention coefficient of the target canopy, and the retention coefficient of the target canopy; The canopy interception amount determination module is used to determine the canopy interception amount of the target canopy within the rainfall duration based on the rainfall intensity, rainfall duration, the canopy interception capacity of the target canopy and the interception ratio of the target canopy.

9. A computer-readable storage medium, characterized in that A laser radar-based canopy interception determination program is stored thereon. When the laser radar-based canopy interception determination program is executed by the processor, the laser radar-based canopy interception determination method described in any one of claims 1 to 7 is implemented.

10. An electronic device, characterized in that: It includes a memory, a processor, and a lidar-based canopy interception determination program stored in the memory and executable on the processor. When the processor executes the lidar-based canopy interception determination program, the lidar-based canopy interception determination method described in any one of claims 1 to 7 is implemented.

Images