Plant root system distribution detection method and device and plant root system distribution detection system

Abstract

The invention provides a plant root system distribution detection method and device and a plant root system distribution detection system.The method comprises the steps that a wave field function of soil where a plant root system is located is obtained, dip angle decomposition functions at different migration speeds are determined according to the wave field function, and the amplitude of a target signal is determined according to the dip angle decomposition functions; calculating the focusing degree of the target signal according to the amplitude to obtain a first weighting factor, and calculating the similarity of the target signal according to the amplitude to obtain a second weighting factor; local features of the target signal are calculated according to the amplitude, a third weighting factor is obtained, a plant root system section function is obtained through calculation according to the first weighting factor, the second weighting factor, the third weighting factor and the dip angle decomposition function, and the plant root system section function is used for spatial distribution imaging of soil where the plant root system is located. According to the method, the problem that in the plant root system detection process in the prior art, the imaging precision of plant root system space distribution in a complex underground environment is not high is solved.

Description

Technical Field

[0001] The present invention relates to the field of plant root system detection, and in particular to a plant root system distribution detection method, device, computer-readable storage medium, computer program product and plant root system distribution detection system. Background Art

[0002] The development of the root system directly affects key ecological processes such as soil structure, nutrient cycling and water conservation. Therefore, accurate detection of the spatial distribution of plant roots is of great significance for studying the health and sustainability of ecosystems. However, traditionally, the detection and research of plant roots mainly rely on excavation and direct observation. This method is not only highly destructive and costly, but also difficult to obtain three-dimensional distribution information of roots in the soil. In recent years, ground penetrating radar, as a non-invasive detection technology, can obtain wave field signals of underground non-uniform geological bodies. Plant roots have obvious heterogeneous characteristics, and radar wave field signals can be used to locate the position of plant roots. However, traditional methods may be difficult to provide sufficient resolution and positioning accuracy when facing complex underground environments.

[0003] Traditional ground penetrating radar imaging methods mostly use migration technology based on a single velocity model. In complex soil environments where physical properties such as soil conductivity and moisture content vary greatly, the variability of soil conductivity and moisture content will cause significant fluctuations in the propagation speed of electromagnetic waves. The single velocity model cannot accurately reflect this non-uniformity, thereby introducing systematic errors in the imaging process, reducing the accuracy and reliability of root imaging, and easily leading to distortion and inaccuracy of imaging results. Summary of the invention

[0004] The main purpose of the present application is to provide a plant root distribution detection method, device, computer-readable storage medium, computer program product and plant root distribution detection system, so as to at least solve the problem in the prior art of low accuracy of imaging the spatial distribution of plant roots in the complex underground environment during plant root detection.

[0005] In order to achieve the above-mentioned purpose, according to one aspect of the present application, a method for detecting plant root distribution is provided, comprising: obtaining a wave field function of the soil where the plant roots are located, the wave field function being a function of the physical characteristics of electromagnetic wave propagation in the soil where the plant roots are located with respect to time, and the physical characteristics including the propagation distance, amplitude and phase of the electromagnetic wave in the soil; determining a dip angle decomposition function at different offset speeds according to the wave field function, the dip angle decomposition function being a function of the physical characteristics of electromagnetic wave propagation at different incident angles separated by the wave field function according to different offset speeds with respect to time; determining the amplitude of a target signal according to the dip angle decomposition function, calculating the focusing degree of the target signal according to the amplitude, and obtaining a first weighting factor, wherein the target signal is the electromagnetic During the propagation of the magnetic wave, the electromagnetic wave signal is reflected back after interacting with the plant root system, the amplitude is the maximum amplitude of the target signal, and the focusing degree is the clarity of the target signal; the similarity of the target signal is calculated according to the amplitude to obtain a second weighting factor, and the similarity is the consistency of the waveform change of the target signal at different offset speeds; the local feature of the target signal is calculated according to the amplitude to obtain a third weighting factor, and the local feature is the change of the amplitude of the target signal in the local area; the plant root profile function is calculated according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function, and the plant root profile function is used for spatial distribution imaging of the soil where the plant root system is located.

[0006] Optionally, determining the dip decomposition function at different migration speeds according to the wave field function includes: performing dip decomposition on the wave field function to decompose wave components with different dips; processing the wave components using different migration speeds to obtain the dip decomposition function at different migration speeds, wherein the dip decomposition function is I(x, v, p), wherein x is the spatial position of the plant root system, v is the migration speed, and p is the formation dip.

[0007] Optionally, the amplitude of the target signal is determined according to the inclination decomposition function, and the focusing degree of the target signal is calculated according to the amplitude to obtain a first weighting factor, including: obtaining the amplitude A(x, v) of the plant root system when the spatial position is x and the offset speed is v based on the inclination decomposition function; determining the energy intensity of the target signal when the plant root system is at the spatial position x and the offset speed v according to the amplitude. Among them, v i is the value of the i-th offset speed among the different offset speeds, and the energy intensity is the total energy of the target signal at the different offset speeds; the focusing degree of the target signal is calculated according to the energy intensity of the target signal The focusing degree is determined as a first weighting factor.

[0008] Optionally, calculating the similarity of the target signal according to the amplitude to obtain a second weighting factor includes: calculating a similarity measure according to the amplitude The similarity measure is used to evaluate the similarity of the amplitude of the target signal at different offset speeds; the variance of the offset speed is calculated based on the similarity measure. is the value of the ith similarity measure in the different similarity measures; the similarity of the target signal is calculated based on the similarity measure and the variance The similarity is determined as a second weighting factor.

[0009] Optionally, calculating the local feature of the target signal according to the amplitude to obtain a third weighting factor includes: calculating the local feature of the target signal according to the amplitude and a Laplace operator. Wherein, Δ is the Laplace operator, λ is the balance parameter, and the local feature is determined as the third weighting factor.

[0010] Optionally, calculating the plant root profile function according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function includes: calculating the plant root profile function according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function Among them, k is the nonlinear adjustment factor and m is the number of weighting factors.

[0011] In order to achieve the above-mentioned purpose, according to one aspect of the present application, a plant root distribution detection device is provided, including: a first acquisition unit, used to acquire the wave field function of the soil where the plant roots are located, the wave field function is the physical property of the propagation of electromagnetic waves in the soil where the plant roots are located, and the physical property includes the propagation speed, amplitude change and phase change of the electromagnetic wave in the soil; a first determination unit, used to determine the dip angle decomposition function at different offset speeds according to the wave field function, the dip angle decomposition function is the physical property of the propagation of the electromagnetic wave at different incident angles separated by the wave field function according to different offset speeds; a first calculation unit, used to determine the amplitude of the target signal according to the dip angle decomposition function, calculate the focusing degree of the target signal according to the amplitude, and obtain a first weighting factor, the target signal is the focusing degree of the electromagnetic wave during propagation. During the process, the amplitude of the electromagnetic wave signal reflected back after interacting with the plant root system is the maximum amplitude of the target signal, and the focusing degree is the clarity of the target signal; the second calculation unit is used to calculate the similarity of the target signal according to the amplitude to obtain a second weighting factor, and the similarity is the consistency of the waveform change of the target signal at different offset speeds; the third calculation unit is used to calculate the local characteristics of the target signal according to the amplitude to obtain a third weighting factor, and the local characteristics are the changes in the amplitude of the target signal in the local area; the fourth calculation unit is used to calculate the plant root profile function according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function, and the plant root profile function is the spatial distribution characteristic of the soil where the plant root system is located.

[0012] According to another aspect of the present application, a computer-readable storage medium is provided, wherein the computer-readable storage medium includes a stored program, wherein when the program is executed, the method described in any one of the devices where the computer-readable storage medium is located is controlled.

[0013] According to another aspect of the present application, a computer program product is provided, comprising a computer program, wherein when the computer program is executed by a processor, any one of the methods described above is implemented.

[0014] According to another aspect of the present application, a plant root distribution detection system is provided, characterized in that it includes: one or more processors, a memory, and one or more programs, wherein the one or more programs are stored in the memory and are configured to be executed by the one or more processors, and the one or more programs include methods for executing any one of the described methods.

[0015] By applying the technical solution of the present application, in the above-mentioned method for detecting the distribution of plant roots, a wave field function of the soil where the above-mentioned plant roots are located is obtained, the above-mentioned wave field function is a function of the physical characteristics of electromagnetic wave propagation in the soil where the above-mentioned plant roots are located with respect to time, and the above-mentioned physical characteristics include the propagation distance, amplitude and phase of the above-mentioned electromagnetic wave in the above-mentioned soil; the dip angle decomposition function at different offset speeds is determined according to the above-mentioned wave field function, and the above-mentioned dip angle decomposition function is a function of the above-mentioned physical characteristics of electromagnetic wave propagation at different incident angles separated by the above-mentioned wave field function according to the above-mentioned offset speeds with respect to time; the amplitude of the target signal is determined according to the above-mentioned dip angle decomposition function, and the focusing degree of the above-mentioned target signal is calculated according to the above-mentioned amplitude to obtain a first weighting factor, and the above-mentioned target signal is the focusing degree of the above-mentioned target signal during the propagation process. In the process, the electromagnetic wave signal reflected back after interacting with the above-mentioned plant root system, the above-mentioned amplitude is the maximum amplitude value of the above-mentioned target signal, and the above-mentioned focusing degree is the clarity of the above-mentioned target signal; the similarity of the above-mentioned target signal is calculated according to the above-mentioned amplitude to obtain the second weighting factor, and the above-mentioned similarity is the consistency of the waveform change of the above-mentioned target signal at different above-mentioned offset speeds; the local feature of the above-mentioned target signal is calculated according to the above-mentioned amplitude to obtain the third weighting factor, and the above-mentioned local feature is the change of the above-mentioned amplitude of the above-mentioned target signal in the local area; the plant root profile function is calculated according to the above-mentioned first weighting factor, the above-mentioned second weighting factor, the above-mentioned third weighting factor and the above-mentioned inclination decomposition function, and the above-mentioned plant root profile function is used for spatial distribution imaging of the soil where the above-mentioned plant root system is located. The present application obtains the wave field function of the soil where the plant roots are located, uses the wave field function to determine the inclination decomposition function at different offset speeds, and then comprehensively processes the inclination function at different offset speeds by calculating the three weighting factors of the focus, similarity and local characteristics of the target signal to obtain the plant root profile function. This function reflects the spatial distribution imaging of the soil where the plant roots are located, realizes high-precision imaging of the plant roots in a complex soil environment, and solves the problem of low accuracy of spatial distribution imaging of plant roots in the complex underground environment during plant root detection in the prior art. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] Figure 1 A hardware structure block diagram of a mobile terminal for detecting plant root distribution according to an embodiment of the present application is shown;

[0017] Figure 2 A schematic diagram of a process of detecting plant root distribution provided according to an embodiment of the present application is shown;

[0018] Figure 3 A structural block diagram of a plant root distribution detection device provided according to an embodiment of the present application is shown.

[0019] The above drawings include the following reference numerals:

[0020] 102, processor; 104, memory; 106, transmission device; 108, input and output devices. DETAILED DESCRIPTION

[0021] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0022] In order to enable those skilled in the art to better understand the solution of the present application, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present application.

[0023] It should be noted that the terms "first", "second", etc. in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequential order. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present application described here. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device comprising a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0024] As introduced in the background technology, in the prior art, when faced with a complex underground environment, the results of plant root imaging are not accurate enough. To solve this technical problem, the embodiments of the present application provide a plant root distribution detection method, device, computer-readable storage medium, computer program product and plant root distribution detection system.

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

[0026] The method embodiments provided in the embodiments of the present application can be executed in a mobile terminal, a computer terminal or a similar computing device. Taking running on a mobile terminal as an example, Figure 1 FIG. 1 is a hardware structure block diagram of a mobile terminal of a plant root distribution detection method according to an embodiment of the present invention. Figure 1 As shown, the mobile terminal may include one or more ( Figure 1Only one is shown in the figure) a processor 102 (the processor 102 may include but is not limited to a processing device such as a microprocessor MCU or a programmable logic device FPGA) and a memory 104 for storing data, wherein the mobile terminal may also include a transmission device 106 and an input / output device 108 for communication functions. It can be understood by those skilled in the art that Figure 1 The structure shown is only for illustration and does not limit the structure of the mobile terminal. Figure 1 More or fewer components as shown, or with Figure 1 Different configurations are shown.

[0027] The memory 104 can be used to store computer programs, for example, software programs and modules of application software, such as a computer program corresponding to a plant root distribution detection method in an embodiment of the present invention. The processor 102 executes various functional applications and data processing by running the computer program stored in the memory 104, that is, the above method is implemented. The memory 104 may include a high-speed random access memory, and may also include a non-volatile memory, such as one or more magnetic storage devices, flash memory, or other non-volatile solid-state memory. In some examples, the memory 104 may further include a memory remotely arranged relative to the processor 102, and these remote memories may be connected to the mobile terminal via a network. Examples of the above-mentioned network include but are not limited to the Internet, an intranet, a local area network, a mobile communication network, and a combination thereof. The transmission device 106 is used to receive or send data via a network. The above-mentioned network specific example may include a wireless network provided by a communication provider of the mobile terminal. In one example, the transmission device 106 includes a network adapter (Network Interface Controller, referred to as NIC), which can be connected to other network devices through a base station so as to communicate with the Internet. In one example, the transmission device 106 may be a radio frequency (RF) module, which is used to communicate with the Internet wirelessly.

[0028] In this embodiment, a plant root distribution detection method running on a mobile terminal, a computer terminal or a similar computing device is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and although the logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0029] Figure 2 1 is a flow chart of a method for detecting plant root distribution according to an embodiment of the present application. Figure 2 As shown, the method comprises the following steps:

[0030] Step S201, obtaining a wave field function of the soil where the plant root system is located, wherein the wave field function is a function of the physical characteristics of electromagnetic wave propagation in the soil where the plant root system is located with respect to time, wherein the physical characteristics include the propagation distance, amplitude and phase of the electromagnetic wave in the soil;

[0031] Specifically, we obtain the wave field data of electromagnetic wave propagation in the soil where the plant roots are located. These data describe the propagation characteristics of electromagnetic waves in the soil medium that change over time, including the distance traveled, amplitude, and phase changes of the waves. By analyzing the wave field function, we can understand the propagation laws of electromagnetic waves in the soil and infer the distribution and status of plant roots based on this.

[0032] Step S202, determining a dip decomposition function at different migration speeds according to the wave field function, wherein the dip decomposition function is a function of the physical characteristics of the electromagnetic wave propagation at different incident angles separated by the wave field function according to different migration speeds with respect to time;

[0033] Specifically, based on the wave field function in the soil, the dip angle decomposition function of the electromagnetic wave propagation characteristics at different incident angles is calculated for different offset velocities. The dip angle decomposition function describes the changes in the physical characteristics of the electromagnetic wave propagation distance, amplitude and phase in the soil over time under different offset velocities. This method enables us to analyze the propagation behavior of electromagnetic waves in the soil from multiple angles, thereby more comprehensively understanding the influence of plant roots on the wave field.

[0034] Step S203, determining the amplitude of the target signal according to the inclination decomposition function, calculating the focusing degree of the target signal according to the amplitude, and obtaining a first weighting factor, wherein the target signal is the electromagnetic wave signal reflected back after the electromagnetic wave interacts with the plant root system during propagation, the amplitude is the maximum value of the amplitude of the target signal, and the focusing degree is the clarity of the target signal;

[0035] Specifically, according to the electromagnetic wave propagation characteristics at different incident angles in the wave field, the amplitude of the electromagnetic wave signal reflected after interacting with the plant roots is extracted, and the focusing degree of the signal is evaluated based on these amplitudes, thereby generating a first weighting factor to measure the signal clarity. In this process, the target signal refers to the electromagnetic wave propagating in the soil and reflected back by the plant roots. The peak value of the amplitude represents the maximum intensity of the target signal, and the focusing degree is a measure of the spatial concentration of the signal, that is, the clarity of the signal.

[0036] Step S204, calculating the similarity of the target signal according to the amplitude to obtain a second weighting factor, wherein the similarity is the consistency of the waveform change of the target signal at different offset speeds;

[0037] Specifically, the amplitude information of the target signal is used to evaluate the consistency of the waveform changes of these signals at different offset speeds, so that the similarity of the target signal is calculated and used as the second weighting factor. The similarity index reflects whether the waveform changes of the reflected wave signal remain consistent under different speed conditions, that is, the stability and reliability of the signal. Similarity analysis can help distinguish between real root signals and other possible interference sources, such as soil structure changes or non-target objects.

[0038] Step S205, calculating the local feature of the target signal according to the amplitude to obtain a third weighting factor, wherein the local feature is a change in the amplitude of the target signal in a local area;

[0039] Specifically, by analyzing the local changes in the target signal amplitude, the amplitude characteristics of the signal in a specific area are extracted, and the local eigenvalue of the signal is calculated and used as the third weighting factor. This local feature reflects the changes in the target signal amplitude within the local range. It is an in-depth exploration of the signal details, which helps us understand the propagation characteristics of electromagnetic wave signals in the soil more carefully, especially in the areas near the plant roots. By analyzing the changes in local amplitude, we can capture subtle differences in the signal, which is crucial to revealing the microstructure and dynamic changes of the root system.

[0040] Step S206, calculating a plant root profile function according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function, wherein the plant root profile function is used for spatial distribution imaging of the soil where the plant roots are located.

[0041] Specifically, the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function are combined to calculate the profile function of the plant root system. This profile function is used to generate a spatial distribution image of the plant root system in the soil, which more accurately depicts the distribution of the plant root system in the soil, improves the accuracy of root detection, and can reveal the detailed characteristics of the root system structure.

[0042] In order to determine the inclination decomposition function, in an optional implementation, the above step S202 includes:

[0043] Step S2021, performing inclination decomposition on the wave field function to decompose wave components with different inclination angles;

[0044] Specifically, the inclination decomposition operation is performed on the wave field function. This step decomposes the wave field function into multiple wave components with different inclinations. In this way, the reflected wave signal after interaction with the plant roots can be more accurately determined. This inclination decomposition method can help us understand the wave propagation path in the soil in detail, which is crucial to improving the accuracy and depth of root detection.

[0045] Step S2022, using different above-mentioned migration speeds to process the above-mentioned wave components, to obtain the above-mentioned inclination decomposition function under different above-mentioned migration speeds, the above-mentioned inclination decomposition function is I(x,v,p), wherein x is the spatial position of the above-mentioned plant root system, v is the migration speed, and p is the formation inclination.

[0046] Specifically, by applying different offset velocities to the wave components of the wave field function, the inclination decomposition functions corresponding to these velocities are obtained. These inclination decomposition functions describe the specific paths of electromagnetic waves propagating in the soil to the plant root position under different offset velocities. The spatial position parameter x, offset velocity v and formation inclination p involved are the key variables affecting the inclination decomposition function.

[0047] In order to obtain the first weighting factor, in an optional implementation manner, the above step S203 includes:

[0048] Step S2031, obtaining the amplitude A(x, v) of the plant root when the spatial position is x and the displacement speed is v based on the inclination angle decomposition function;

[0049] Specifically, the amplitude of the plant root system under a specific spatial position x and a specific offset velocity v is determined according to the inclination decomposition function. Through this method, we can capture the details of the interaction between the plant root system and the electromagnetic waves, and then obtain the precise amplitude value of the root system, which helps to deeply understand the spatial structure and dynamic changes of the root system.

[0050] Step S2032, determining the energy intensity of the target signal when the plant root is at the spatial position x and the displacement speed v according to the amplitude Among them, v i is the value of the i-th offset speed among the different offset speeds, and the energy intensity is the total energy of the target signal at the different offset speeds;

[0051] Specifically, the energy intensity of the target signal of the plant root at a specific spatial position x and a specific offset speed v is calculated according to the amplitude value. The energy intensity is the sum of the energy values ​​of the target signal at all different offset speeds, reflecting the overall energy characteristics of the target signal under multiple speed conditions.

[0052] Step S2033: Calculate the focus of the target signal according to the energy intensity of the target signal. The above-mentioned focusing degree is determined as a first weighting factor.

[0053] Specifically, by analyzing the energy intensity of the target signal, we obtain the focus of the target signal and use it as the first weighting factor. The focus measures the spatial concentration of the signal, that is, the clarity of the signal. Through this method, we can quantify the spatial concentration of the target signal and use it as an important indicator for evaluating the quality of the target signal.

[0054] In order to obtain the second weighting factor, in an optional implementation manner, the above step S204 includes:

[0055] Step S2041, calculate the similarity measure based on the above amplitude The similarity measure is used to evaluate the similarity of the amplitude of the target signal at different migration speeds;

[0056] Specifically, the similarity measure is calculated based on the amplitude value of the target signal to measure the similarity of the target signal amplitude under different offset speeds. The similarity of the signal amplitude under various speed conditions is evaluated, which provides us with an indicator to judge the stability and reliability of the signal.

[0057] Step S2042: Calculate the variance of the above-mentioned offset velocity based on the above-mentioned similarity measurement is the value of the i-th similarity measure among the different similarity measures;

[0058] Specifically, based on the similarity measurement of the target signal at different offset speeds, the variance of these similarity values ​​was calculated. The similarity measurement reflects the degree of discreteness of the amplitude similarity at different speeds, thereby evaluating the performance of the signal at different offset speeds. A high variance indicates that the similarity value varies greatly at different speeds, which means that the signal is disturbed or unstable. A low variance indicates that the similarity values ​​are more concentrated, indicating that the signal has good consistency and stability at different speeds.

[0059] Step S2043, obtaining the similarity of the target signal according to the similarity measurement and the variance calculation The above similarity is determined as a second weighting factor.

[0060] Specifically, the similarity of the target signals was comprehensively evaluated based on the similarity measure of the target signals and their variance. The similarity index reflects whether the waveform change of the target signal remains consistent under different speed conditions, that is, the stability of the signal. It is set as the second weighting factor, providing an indicator to measure the consistency of the signal under speed changes.

[0061] In order to obtain the third weighting factor, in an optional implementation manner, the above step S205 includes:

[0062] Step S2051, calculating the local features of the target signal based on the amplitude and Laplace operator Wherein, Δ is the Laplace operator, λ is the balance parameter, and the above local feature is determined as the third weighting factor.

[0063] Specifically, based on the amplitude information and Laplace operator of the target signal and combined with the balance parameter, the local characteristics of the target signal are calculated and used as the third weighting factor. The local characteristics reflect the local changes of the signal amplitude in space. The Laplace operator is used to capture these changes, and the balance parameter is used to adjust the sensitivity of change capture, which provides more detailed spatial information for target signal analysis and helps to reveal the microstructure and dynamic characteristics of the target signal.

[0064] In order to obtain the plant root profile function, in an optional implementation, the above step S206 includes:

[0065] Step S2061: Calculate the plant root profile function according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function. Among them, k is the nonlinear adjustment factor and m is the number of weighting factors.

[0066] Specifically, the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function are integrated, and the nonlinear adjustment factor k is introduced to calculate the profile function of the plant root system. This profile function integrates the influence of multiple factors on root imaging and provides a comprehensive imaging result.

[0067] It should be noted that the steps shown in the flowcharts of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and that, although a logical order is shown in the flowcharts, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0068] The embodiment of the present application also provides a plant root distribution detection device. It should be noted that a plant root distribution detection device of the embodiment of the present application can be used to execute a plant root distribution detection method provided by the embodiment of the present application. The device is used to implement the above-mentioned embodiments and preferred implementation modes, and the descriptions that have been made will not be repeated. As used below, the term "module" can implement a combination of software and / or hardware for a predetermined function. Although the device described in the following embodiments is preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceivable.

[0069] The following is an introduction to a plant root distribution detection device provided in an embodiment of the present application.

[0070] Figure 3 is a structural block diagram of a plant root distribution detection device according to an embodiment of the present application. Figure 3 As shown, the device comprises:

[0071] A first acquisition unit 10 is used to acquire a wave field function of the soil where the plant root system is located, where the wave field function is a physical property of electromagnetic wave propagation in the soil where the plant root system is located, and the physical property includes a propagation speed, amplitude change, and phase change of the electromagnetic wave in the soil;

[0072] Specifically, we obtain the wave field data of electromagnetic wave propagation in the soil where the plant roots are located. These data describe the propagation characteristics of electromagnetic waves in the soil medium that change over time, including the distance traveled, amplitude, and phase changes of the waves. By analyzing the wave field function, we can understand the propagation laws of electromagnetic waves in the soil and infer the distribution and status of plant roots based on this.

[0073] A first determining unit 20 is used to determine a dip decomposition function at different migration speeds according to the wave field function, wherein the dip decomposition function is the physical characteristic of the electromagnetic wave propagation at different incident angles separated by the wave field function according to different migration speeds;

[0074] Specifically, based on the wave field function in the soil, the dip angle decomposition function of the electromagnetic wave propagation characteristics at different incident angles is calculated for different offset velocities. The dip angle decomposition function describes the changes in the physical characteristics of the electromagnetic wave propagation distance, amplitude and phase in the soil over time under different offset velocities. This method enables us to analyze the propagation behavior of electromagnetic waves in the soil from multiple angles, thereby more comprehensively understanding the influence of plant roots on the wave field.

[0075] A first calculation unit 30 is used to determine the amplitude of the target signal according to the inclination decomposition function, calculate the focusing degree of the target signal according to the amplitude, and obtain a first weighting factor, wherein the target signal is the electromagnetic wave signal reflected back after the electromagnetic wave interacts with the plant root system during propagation, the amplitude is the maximum value of the amplitude of the target signal, and the focusing degree is the clarity of the target signal;

[0076] Specifically, according to the electromagnetic wave propagation characteristics at different incident angles in the wave field, the amplitude of the electromagnetic wave signal reflected after interacting with the plant roots is extracted, and the focusing degree of the signal is evaluated based on these amplitudes, thereby generating a first weighting factor to measure the signal clarity. In this process, the target signal refers to the electromagnetic wave propagating in the soil and reflected back by the plant roots. The peak value of the amplitude represents the maximum intensity of the target signal, and the focusing degree is a measure of the spatial concentration of the signal, that is, the clarity of the signal.

[0077] A second calculation unit 40 is used to calculate the similarity of the target signal according to the amplitude to obtain a second weighting factor, wherein the similarity is the consistency of the waveform change of the target signal at different offset speeds;

[0078] Specifically, the amplitude information of the target signal is used to evaluate the consistency of the waveform changes of these signals at different offset speeds, so that the similarity of the target signal is calculated and used as the second weighting factor. The similarity index reflects whether the waveform changes of the reflected wave signal remain consistent under different speed conditions, that is, the stability and reliability of the signal. Similarity analysis can help distinguish between real root signals and other possible interference sources, such as soil structure changes or non-target objects.

[0079] A third calculation unit 50 is used to calculate the local feature of the target signal according to the amplitude to obtain a third weighting factor, wherein the local feature is a change in the amplitude of the target signal in a local area;

[0080] Specifically, by analyzing the local changes in the target signal amplitude, the amplitude characteristics of the signal in a specific area are extracted, and the local eigenvalue of the signal is calculated and used as the third weighting factor. This local feature reflects the changes in the target signal amplitude within the local range. It is an in-depth exploration of the signal details, which helps us understand the propagation characteristics of electromagnetic wave signals in the soil more carefully, especially in the areas near the plant roots. By analyzing the changes in local amplitude, we can capture subtle differences in the signal, which is crucial to revealing the microstructure and dynamic changes of the root system.

[0081] The fourth calculation unit 60 is used to calculate the plant root profile function according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function, wherein the plant root profile function is the spatial distribution characteristic of the soil where the plant roots are located.

[0082] Specifically, the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function are combined to calculate the profile function of the plant root system. This profile function is used to generate a spatial distribution image of the plant root system in the soil, which more accurately depicts the distribution of the plant root system in the soil, improves the accuracy of root detection, and can reveal the detailed characteristics of the root system structure.

[0083] In order to determine the inclination decomposition function, in an optional implementation manner, the first determining unit includes:

[0084] The first determination module is used to perform inclination decomposition on the wave field function to decompose wave components with different inclination angles;

[0085] Specifically, the inclination decomposition operation is performed on the wave field function. This step decomposes the wave field function into multiple wave components with different inclinations. In this way, the reflected wave signal after interaction with the plant roots can be more accurately determined. This inclination decomposition method can help us understand the wave propagation path in the soil in detail, which is crucial to improving the accuracy and depth of root detection.

[0086] The second determination module is used to process the above-mentioned wave components using different above-mentioned migration speeds to obtain the above-mentioned inclination decomposition function under different above-mentioned migration speeds. The above-mentioned inclination decomposition function is I(x,v,p), wherein x is the spatial position of the above-mentioned plant root system, v is the migration speed, and p is the formation inclination.

[0087] Specifically, by applying different offset velocities to the wave components of the wave field function, the inclination decomposition functions corresponding to these velocities are obtained. These inclination decomposition functions describe the specific paths of electromagnetic waves propagating in the soil to the plant root position under different offset velocities. The spatial position parameter x, offset velocity v and formation inclination p involved are the key variables affecting the inclination decomposition function.

[0088] In order to obtain the first weighting factor, in an optional implementation manner, the first calculation unit includes:

[0089] A first calculation module is used to obtain the amplitude A(x, v) of the plant root when the spatial position is x and the displacement speed is v based on the inclination angle decomposition function;

[0090] Specifically, the amplitude of the plant root system under a specific spatial position x and a specific offset velocity v is determined according to the inclination decomposition function. Through this method, we can capture the details of the interaction between the plant root system and the electromagnetic waves, and then obtain the precise amplitude value of the root system, which helps to deeply understand the spatial structure and dynamic changes of the root system.

[0091] The second calculation module is used to determine the energy intensity of the target signal when the plant root is at the spatial position x and the displacement speed v according to the amplitude. Among them, v i is the value of the i-th offset speed among the different offset speeds, and the energy intensity is the total energy of the target signal at the different offset speeds;

[0092] Specifically, the energy intensity of the target signal of the plant root at a specific spatial position x and a specific offset speed v is calculated according to the amplitude value. The energy intensity is the sum of the energy values ​​of the target signal at all different offset speeds, reflecting the overall energy characteristics of the target signal under multiple speed conditions.

[0093] The third calculation module is used to calculate the focus degree of the target signal according to the energy intensity of the target signal. The above-mentioned focusing degree is determined as a first weighting factor.

[0094] Specifically, by analyzing the energy intensity of the target signal, we obtain the focus of the target signal and use it as the first weighting factor. The focus measures the spatial concentration of the signal, that is, the clarity of the signal. Through this method, we can quantify the spatial concentration of the target signal and use it as an important indicator for evaluating the quality of the target signal.

[0095] In order to obtain the second weighting factor, in an optional implementation manner, the second calculation unit includes:

[0096] The fourth calculation module is used to calculate the similarity measure based on the above amplitude The similarity measure is used to evaluate the similarity of the amplitude of the target signal at different migration speeds;

[0097] Specifically, the similarity measure is calculated based on the amplitude value of the target signal to measure the similarity of the target signal amplitude under different offset speeds. The similarity of the signal amplitude under various speed conditions is evaluated, which provides us with an indicator to judge the stability and reliability of the signal.

[0098] The fifth calculation module is used to calculate the variance of the above-mentioned offset speed according to the above-mentioned similarity measurement is the value of the i-th similarity measure among the different similarity measures;

[0099] Specifically, based on the similarity measurement of the target signal at different offset speeds, the variance of these similarity values ​​was calculated. The similarity measurement reflects the degree of discreteness of the amplitude similarity at different speeds, thereby evaluating the performance of the signal at different offset speeds. A high variance indicates that the similarity value varies greatly at different speeds, which means that the signal is disturbed or unstable. A low variance indicates that the similarity values ​​are more concentrated, indicating that the signal has good consistency and stability at different speeds.

[0100] A sixth calculation module, configured to obtain the similarity of the target signal according to the similarity measurement and the variance calculation. The above similarity is determined as a second weighting factor.

[0101] Specifically, the similarity of the target signals was comprehensively evaluated based on the similarity measure of the target signals and their variance. The similarity index reflects whether the waveform change of the target signal remains consistent under different speed conditions, that is, the stability of the signal. It is set as the second weighting factor, providing an indicator to measure the consistency of the signal under speed changes.

[0102] In order to obtain the third weighting factor, in an optional implementation manner, the third calculation unit includes:

[0103] A seventh calculation module, used to calculate the local features of the target signal according to the amplitude and the Laplace operator Wherein, Δ is the Laplace operator, λ is the balance parameter, and the above local feature is determined as the third weighting factor.

[0104] Specifically, based on the amplitude information and Laplace operator of the target signal and combined with the balance parameter, the local characteristics of the target signal are calculated and used as the third weighting factor. The local characteristics reflect the local changes of the signal amplitude in space. The Laplace operator is used to capture these changes, and the balance parameter is used to adjust the sensitivity of change capture, which provides more detailed spatial information for target signal analysis and helps to reveal the microstructure and dynamic characteristics of the target signal.

[0105] In order to obtain the plant root profile function, in an optional implementation, the fourth calculation unit includes:

[0106] An eighth calculation module is used to calculate the plant root profile function according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function. Among them, k is the nonlinear adjustment factor and m is the number of weighting factors.

[0107] Specifically, the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function are integrated, and the nonlinear adjustment factor k is introduced to calculate the profile function of the plant root system. This profile function integrates the influence of multiple factors on root imaging and provides a comprehensive imaging result.

[0108] The plant root distribution detection device includes a processor and a memory, wherein the first acquisition unit, the first determination unit, the first calculation unit, the second calculation unit, the third calculation unit and the fourth calculation unit are stored in the memory as program units, and the processor executes the program units stored in the memory to implement corresponding functions. The modules are all located in the same processor; or, the modules are located in different processors in any combination.

[0109] The processor includes a kernel, which retrieves the corresponding program unit from the memory. One or more kernels can be set, and the imaging accuracy of the spatial distribution of plant roots can be improved by adjusting kernel parameters.

[0110] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM, and the memory includes at least one memory chip.

[0111] An embodiment of the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium includes a stored program, wherein when the program is executed, the device where the computer-readable storage medium is located is controlled to execute the plant root distribution detection method.

[0112] Specifically, a plant root distribution detection method includes:

[0113] Step S201, obtaining a wave field function of the soil where the plant root system is located, wherein the wave field function is a function of the physical characteristics of electromagnetic wave propagation in the soil where the plant root system is located with respect to time, wherein the physical characteristics include the propagation distance, amplitude and phase of the electromagnetic wave in the soil;

[0114] Step S202, determining a dip decomposition function at different migration speeds according to the wave field function, wherein the dip decomposition function is a function of the physical characteristics of the electromagnetic wave propagation at different incident angles separated by the wave field function according to different migration speeds with respect to time;

[0115] Step S203, determining the amplitude of the target signal according to the inclination decomposition function, calculating the focusing degree of the target signal according to the amplitude, and obtaining a first weighting factor, wherein the target signal is the electromagnetic wave signal reflected back after the electromagnetic wave interacts with the plant root system during propagation, the amplitude is the maximum value of the amplitude of the target signal, and the focusing degree is the clarity of the target signal;

[0116] Step S204, calculating the similarity of the target signal according to the amplitude to obtain a second weighting factor, wherein the similarity is the consistency of the waveform change of the target signal at different offset speeds;

[0117] Step S205, calculating the local feature of the target signal according to the amplitude to obtain a third weighting factor, wherein the local feature is a change in the amplitude of the target signal in a local area;

[0118] Step S206, calculating a plant root profile function according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function, wherein the plant root profile function is used for spatial distribution imaging of the soil where the plant roots are located.

[0119] An embodiment of the present invention provides a processor, and the processor is used to run a program, wherein the program executes the plant root distribution detection method when running.

[0120] Specifically, a plant root distribution detection method includes:

[0121] Step S201, obtaining a wave field function of the soil where the plant root system is located, wherein the wave field function is a function of the physical characteristics of electromagnetic wave propagation in the soil where the plant root system is located with respect to time, wherein the physical characteristics include the propagation distance, amplitude and phase of the electromagnetic wave in the soil;

[0122] Step S202, determining a dip decomposition function at different migration speeds according to the wave field function, wherein the dip decomposition function is a function of the physical characteristics of the electromagnetic wave propagation at different incident angles separated by the wave field function according to different migration speeds with respect to time;

[0123] Step S203, determining the amplitude of the target signal according to the inclination decomposition function, calculating the focusing degree of the target signal according to the amplitude, and obtaining a first weighting factor, wherein the target signal is the electromagnetic wave signal reflected back after the electromagnetic wave interacts with the plant root system during propagation, the amplitude is the maximum value of the amplitude of the target signal, and the focusing degree is the clarity of the target signal;

[0124] Step S204, calculating the similarity of the target signal according to the amplitude to obtain a second weighting factor, wherein the similarity is the consistency of the waveform change of the target signal at different offset speeds;

[0125] Step S205, calculating the local feature of the target signal according to the amplitude to obtain a third weighting factor, wherein the local feature is a change in the amplitude of the target signal in a local area;

[0126] Step S206, calculating a plant root profile function according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function, wherein the plant root profile function is used for spatial distribution imaging of the soil where the plant roots are located.

[0127] The present application also provides a computer program product, which, when executed on a data processing device, is suitable for executing a program for initializing at least the following method steps:

[0128] Step S201, obtaining a wave field function of the soil where the plant root system is located, wherein the wave field function is a function of the physical characteristics of electromagnetic wave propagation in the soil where the plant root system is located with respect to time, wherein the physical characteristics include the propagation distance, amplitude and phase of the electromagnetic wave in the soil;

[0129] Step S202, determining a dip decomposition function at different migration speeds according to the wave field function, wherein the dip decomposition function is a function of the physical characteristics of the electromagnetic wave propagation at different incident angles separated by the wave field function according to different migration speeds with respect to time;

[0130] Step S203, determining the amplitude of the target signal according to the inclination decomposition function, calculating the focusing degree of the target signal according to the amplitude, and obtaining a first weighting factor, wherein the target signal is the electromagnetic wave signal reflected back after the electromagnetic wave interacts with the plant root system during propagation, the amplitude is the maximum value of the amplitude of the target signal, and the focusing degree is the clarity of the target signal;

[0131] Step S204, calculating the similarity of the target signal according to the amplitude to obtain a second weighting factor, wherein the similarity is the consistency of the waveform change of the target signal at different offset speeds;

[0132] Step S205, calculating the local feature of the target signal according to the amplitude to obtain a third weighting factor, wherein the local feature is a change in the amplitude of the target signal in a local area;

[0133] Step S206, calculating a plant root profile function according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function, wherein the plant root profile function is used for spatial distribution imaging of the soil where the plant roots are located.

[0134] An embodiment of the present application also provides a plant root distribution detection system, comprising: one or more processors, a memory, and one or more programs, wherein the one or more programs are stored in the memory and are configured to be executed by the one or more processors, including executing any one of the above-mentioned methods for evaluating the power supply capacity of a feeder group.

[0135] Specifically, a plant root distribution detection method includes:

[0136] Step S201, obtaining a wave field function of the soil where the plant root system is located, wherein the wave field function is a function of the physical characteristics of electromagnetic wave propagation in the soil where the plant root system is located with respect to time, wherein the physical characteristics include the propagation distance, amplitude and phase of the electromagnetic wave in the soil;

[0137] Step S202, determining a dip decomposition function at different migration speeds according to the wave field function, wherein the dip decomposition function is a function of the physical characteristics of the electromagnetic wave propagation at different incident angles separated by the wave field function according to different migration speeds with respect to time;

[0138] Step S203, determining the amplitude of the target signal according to the inclination decomposition function, calculating the focusing degree of the target signal according to the amplitude, and obtaining a first weighting factor, wherein the target signal is the electromagnetic wave signal reflected back after the electromagnetic wave interacts with the plant root system during propagation, the amplitude is the maximum value of the amplitude of the target signal, and the focusing degree is the clarity of the target signal;

[0139] Step S204, calculating the similarity of the target signal according to the amplitude to obtain a second weighting factor, wherein the similarity is the consistency of the waveform change of the target signal at different offset speeds;

[0140] Step S205, calculating the local feature of the target signal according to the amplitude to obtain a third weighting factor, wherein the local feature is a change in the amplitude of the target signal in a local area;

[0141] Step S206, calculating a plant root profile function according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function, wherein the plant root profile function is used for spatial distribution imaging of the soil where the plant roots are located.

[0142] Obviously, those skilled in the art should understand that the above modules or steps of the present invention can be implemented by a general computing device, they can be concentrated on a single computing device, or distributed on a network composed of multiple computing devices, they can be implemented by a program code executable by a computing device, so that they can be stored in a storage device and executed by the computing device, and in some cases, the steps shown or described can be executed in a different order than here, or they can be made into individual integrated circuit modules, or multiple modules or steps therein can be made into a single integrated circuit module for implementation. Thus, the present invention is not limited to any specific combination of hardware and software.

[0143] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0144] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0145] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0146] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process in the computer or other programmable device. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0147] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0148] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0149] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0150] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

[0151] From the above description, it can be seen that the above embodiments of the present application achieve the following technical effects:

[0152] 1) In a method for detecting the distribution of plant roots of the present application, a wave field function of the soil where the plant roots are located is obtained, the wave field function is a function of the physical characteristics of electromagnetic wave propagation in the soil where the plant roots are located with respect to time, and the physical characteristics include the propagation distance, amplitude and phase of the electromagnetic wave in the soil; a dip angle decomposition function at different offset speeds is determined according to the wave field function, the dip angle decomposition function is a function of the physical characteristics of electromagnetic wave propagation at different incident angles separated by the wave field function according to different offset speeds with respect to time; the amplitude of the target signal is determined according to the dip angle decomposition function, the focusing degree of the target signal is calculated according to the amplitude, and a first weighting factor is obtained, the target signal is the electromagnetic wave in the propagation process, The electromagnetic wave signal reflected after interacting with the above-mentioned plant root system, the above-mentioned amplitude is the maximum amplitude of the above-mentioned target signal, and the above-mentioned focusing degree is the clarity of the above-mentioned target signal; the similarity of the above-mentioned target signal is calculated according to the above-mentioned amplitude to obtain the second weighting factor, and the above-mentioned similarity is the consistency of the waveform change of the above-mentioned target signal at different above-mentioned offset speeds; the local characteristics of the above-mentioned target signal are calculated according to the above-mentioned amplitude to obtain the third weighting factor, and the above-mentioned local characteristics are the changes in the above-mentioned amplitude of the above-mentioned target signal in the local area; the plant root profile function is calculated according to the above-mentioned first weighting factor, the above-mentioned second weighting factor, the above-mentioned third weighting factor and the above-mentioned inclination decomposition function, and the above-mentioned plant root profile function is used for spatial distribution imaging of the soil where the above-mentioned plant roots are located. The present application obtains the wave field function of the soil where the plant roots are located, uses the wave field function to determine the inclination decomposition function at different offset speeds, and then comprehensively processes the inclination function at different offset speeds by calculating the three weighting factors of the focus, similarity and local characteristics of the target signal to obtain the plant root profile function. This function reflects the spatial distribution imaging of the soil where the plant roots are located, realizes high-precision imaging of the plant roots in a complex soil environment, and solves the problem of low accuracy of spatial distribution imaging of plant roots in the complex underground environment during plant root detection in the prior art.

[0153] 2) In a plant root distribution detection device of the present application, a first acquisition unit is used to obtain a wave field function of the soil where the plant roots are located, the wave field function is a physical characteristic of electromagnetic wave propagation in the soil where the plant roots are located, and the physical characteristic includes the propagation speed, amplitude change and phase change of the electromagnetic wave in the soil; a first determination unit is used to determine a dip angle decomposition function at different offset speeds according to the wave field function, the dip angle decomposition function is the physical characteristic of electromagnetic wave propagation at different incident angles separated by the wave field function according to the different offset speeds; a first calculation unit is used to determine the amplitude of a target signal according to the dip angle decomposition function, calculate the focusing degree of the target signal according to the amplitude, and obtain a first weighting factor, the target signal is the amplitude of the electromagnetic wave during propagation, and the focusing degree of the target signal is calculated according to the amplitude, and a first weighting factor is obtained. The electromagnetic wave signal reflected back after the root interaction, the above amplitude is the maximum amplitude of the above target signal, and the above focusing degree is the clarity of the above target signal; the second calculation unit is used to calculate the similarity of the above target signal according to the above amplitude, and obtain the second weighting factor, and the above similarity is the consistency of the waveform change of the above target signal at different above offset speeds; the third calculation unit is used to calculate the local characteristics of the above target signal according to the above amplitude, and obtain the third weighting factor, and the above local characteristics are the changes in the above amplitude of the above target signal in the local area; the fourth calculation unit is used to calculate the plant root profile function according to the above first weighting factor, the above second weighting factor, the above third weighting factor and the above inclination decomposition function, and the above plant root profile function is the spatial distribution characteristic of the soil where the above plant roots are located. The present application obtains the wave field function of the soil where the plant roots are located, uses the wave field function to determine the inclination decomposition function at different offset speeds, and then comprehensively processes the inclination function at different offset speeds by calculating the three weighting factors of the focus, similarity and local characteristics of the target signal to obtain the plant root profile function. This function reflects the spatial distribution imaging of the soil where the plant roots are located, realizes high-precision imaging of the plant roots in a complex soil environment, and solves the problem of low accuracy of spatial distribution imaging of plant roots in the complex underground environment during plant root detection in the prior art.

[0154] The above description is only the preferred embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method for detecting plant root distribution, characterized in that: include: Acquire a wave field function of the soil where the plant root system is located, where the wave field function is a function of the physical characteristics of electromagnetic wave propagation in the soil where the plant root system is located with respect to time, where the physical characteristics include the propagation distance, amplitude and phase of the electromagnetic wave in the soil; Determining a dip angle decomposition function at different migration speeds according to the wave field function, wherein the dip angle decomposition function is a function of the physical characteristics of the electromagnetic wave propagation at different incident angles separated by the wave field function according to different migration speeds with respect to time; Determine the amplitude of the target signal according to the inclination decomposition function, calculate the focusing degree of the target signal according to the amplitude, and obtain a first weighting factor, wherein the target signal is the electromagnetic wave signal reflected back after the electromagnetic wave interacts with the plant root system during propagation, the amplitude is the maximum value of the amplitude of the target signal, and the focusing degree is the clarity of the target signal; Calculating the similarity of the target signal according to the amplitude to obtain a second weighting factor, wherein the similarity is the consistency of the waveform change of the target signal at different offset speeds; Calculating a local feature of the target signal according to the amplitude to obtain a third weighting factor, wherein the local feature is a change in the amplitude of the target signal in a local area; A plant root profile function is calculated based on the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function, and the plant root profile function is used for spatial distribution imaging of the soil where the plant roots are located.

2. The method according to claim 1, characterized in that Determining the dip decomposition function at different migration velocities according to the wave field function includes: Performing inclination decomposition on the wave field function to decompose wave components with different inclination angles; The wave component is processed using different offset speeds to obtain the inclination decomposition function under different offset speeds, and the inclination decomposition function is I(x, v, p), where x is the spatial position of the plant root system, v is the offset speed, and p is the formation inclination.

3. The method according to claim 2, characterized in that Determining the amplitude of the target signal according to the inclination decomposition function, calculating the focusing degree of the target signal according to the amplitude, and obtaining a first weighting factor, including: Based on the inclination decomposition function, the amplitude A(x, v) of the plant root is obtained when the spatial position is x and the displacement speed is v; Determine the energy intensity of the target signal when the spatial position of the plant root is x and the offset speed is v according to the amplitude Among them, v i is the value of the i-th said offset speed among the different said offset speeds, and the said energy intensity is the total energy of the said target signal at the different said offset speeds; The focusing degree of the target signal is calculated according to the energy intensity of the target signal The focusing degree is determined as a first weighting factor.

4. The method according to claim 3, characterized in that Calculating the similarity of the target signal according to the amplitude to obtain a second weighting factor includes: A similarity measure is calculated based on the amplitude The similarity measure is used to evaluate the similarity of the amplitude of the target signal at different migration speeds; The variance of the offset velocity is calculated based on the similarity measure is the value of the i-th similarity measure among the different similarity measures; The similarity of the target signal is obtained according to the similarity measure and the variance calculation The similarity is determined as a second weighting factor.

5. The method according to claim 4, characterized in that Calculating the local feature of the target signal according to the amplitude to obtain a third weighting factor includes: The local feature of the target signal is obtained by calculating the amplitude and the Laplace operator Wherein, Δ is the Laplace operator, λ is the balance parameter, and the local feature is determined as the third weighting factor.

6. The method according to claim 5, characterized in that The plant root profile function is calculated according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function, including: The plant root profile function is obtained by calculating the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function. Among them, k is the nonlinear adjustment factor and m is the number of weighting factors.

7. A plant root distribution detection device, characterized in that: include: A first acquisition unit is used to acquire a wave field function of the soil where the plant root system is located, where the wave field function is a physical property of electromagnetic wave propagation in the soil where the plant root system is located, and the physical property includes a propagation speed, amplitude change, and phase change of the electromagnetic wave in the soil; A first determining unit is used to determine a dip decomposition function at different migration speeds according to the wave field function, wherein the dip decomposition function is the physical property of the electromagnetic wave propagation at different incident angles separated by the wave field function according to different migration speeds; a first calculation unit, configured to determine the amplitude of a target signal according to the inclination decomposition function, calculate the focusing degree of the target signal according to the amplitude, and obtain a first weighting factor, wherein the target signal is an electromagnetic wave signal reflected back after the electromagnetic wave interacts with the plant root system during propagation, the amplitude is a maximum value of the amplitude of the target signal, and the focusing degree is a clarity of the target signal; a second calculation unit, configured to calculate the similarity of the target signal according to the amplitude to obtain a second weighting factor, wherein the similarity is the consistency of the waveform change of the target signal at different offset speeds; a third calculation unit, configured to calculate a local feature of the target signal according to the amplitude to obtain a third weighting factor, wherein the local feature is a change in the amplitude of the target signal in a local area; The fourth calculation unit is used to calculate the plant root profile function according to the first weighting factor, the second weighting factor, the third weighting factor and the inclination decomposition function, wherein the plant root profile function is the spatial distribution characteristic of the soil where the plant roots are located.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium includes a stored program, wherein when the program is executed, the device where the computer-readable storage medium is located is controlled to execute the method according to any one of claims 1 to 6.

9. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the method according to any one of claims 1 to 6 is implemented.

10. A plant root distribution detection system, characterized in that: include: One or more processors, a memory, and one or more programs, wherein the one or more programs are stored in the memory and are configured to be executed by the one or more processors, and the one or more programs include methods for executing any one of claims 1 to 6.

Images