A soil profile sample analysis system and method

Abstract

The present invention discloses a soil profile sample analysis system and method, which relates to the technical field of soil analysis data management. By identifying the hyperspectral images of the soil, a certain sampling record in the historical sampling records is recorded as the target sampling. The change functions corresponding to all types of soil components in the target slice are obtained, the number of target changes in all the change functions is calculated, the corresponding values of the soil component contents in all the unit sampling layers of each profile sampling slice are obtained, the change ranges of the contents of various soil components are calculated respectively, the current sampling record of the soil is obtained, the weight coefficient is calculated by comparing the differences in the change ranges of the soil component contents between the current sampling record and the historical sampling record, the differences in the target changes between the current sampling record and the historical sampling record are compared, and the difference evaluation value is obtained by measuring the differences with the weight coefficient.

Description

Technical Field

[0001] The present invention relates to the technical field of soil analysis data management, and specifically to a soil profile sample analysis system and method. Background Art

[0002] By analyzing the soil in the garden, the obtained soil data is helpful for the operation decision-making of the garden development. In order to reduce the impact on the garden landform, soil information in the garden is obtained by means of profile sampling of the soil. In actual application scenarios, random sampling is usually used to obtain soil samples in the garden. Since the number of soil samples obtained each time is limited, the range of garden land covered by the soil samples is limited, resulting in incomplete sampling data, which in turn affects the judgment of soil data by relevant management personnel. Summary of the Invention

[0003] The purpose of the present invention is to provide a soil profile sample analysis system and method to solve the problems raised in the prior art.

[0004] To achieve the above purpose, the present invention provides the following technical solution: A soil profile sample analysis method, the method includes:

[0005] Step S100: Obtain a profile sampling template of the soil, collect sampling records through the profile sampling template to obtain a profile sampling slice of the soil, and identify the hyperspectral image of the soil to obtain a change function of the corresponding values of the soil component content on the profile sampling slice;

[0006] Step S200: Obtain the historical sampling records of the soil, record a certain sampling record in the historical sampling records as the target sampling, obtain all the profile sampling slices of the target sampling, take a certain profile sampling slice among all the profile sampling slices as the target slice, for the change functions corresponding to all types of soil components in the target slice, take the change rate greater than the threshold in the change function as the target change, and calculate the number of target changes in all the change functions;

[0007] Step S300: Obtain all the profile sampling slices of several land sample collections, calculate the mean value of the number of all target changes, and take the mean value as the sampling determination threshold;

[0008] Step S400: Obtain all the profile sampling slices in the target sampling, obtain the corresponding values of the soil component content in all the unit sampling layers of each profile sampling slice, and calculate the change ranges of various soil component contents respectively;

[0009] Step S500: Obtain the current sampling record of the soil. Calculate the weight coefficient by comparing the difference in the change range of soil component contents between the current sampling record and the historical sampling record. Compare the difference in the target change between the current sampling record and the historical sampling record. Measure the difference through the weight coefficient to obtain a difference evaluation value. When the difference evaluation value is less than the sampling determination threshold, an information prompt is given.

[0010] Further, step S100 includes:

[0011] Step S101: Obtain a profile sampling template with a depth of W1 and a width of W2, satisfying the condition W1 > W2. The profile sampling template includes n unit sampling layers. Save the sampling samples of the soil through the profile sampling template to obtain a profile sampling slice of the soil. Save all the profile sampling slices in the same direction, and the direction is from the ground to the underground or from the underground to the ground;

[0012] Divide the profile sampling template into a short side and a long side. In this solution, analyze the change trend of the soil component content in the long side direction of the profile sampling slice;

[0013] Since the physical or chemical properties of different soil layers are different, the diffusion methods of soil components in the soil are different. By capturing the rate of diffusion change, obtain the soil layer information of component diffusion, and further obtain the range of the soil layer covered by the sampling sample;

[0014] Step S102: Identify k soil components of the profile sampling slice respectively to obtain a soil component change function. The change function of the i-th soil component among the k soil components is denoted as Gi(λ), where λ is the distance variable along the direction of the profile sampling slice.

[0015] Further, step S200 includes:

[0016] Step S201: Obtain the change function Gi ’ i(λ) of the i-th soil component in the target slice, and obtain the derivative function g ’ i(λ) of Gi(λ) in the direction of the sampling slice; ’ i(λ);

[0017] Step S202: Set a change rate threshold f, and record a continuous part greater than f in gi ’ i(λ) as a target change, and obtain the number of target changes in gi ’ i(λ);

[0018] Step S203: Obtain the soil component change functions corresponding to the k soil components in the target slice respectively, obtain the target changes corresponding to the soil component change functions, and obtain the total number of all target changes in the target slice, denoted as m.

[0019] Further, step S300 includes:

[0020] Step S301: Obtain d1 land sample collection records. Each land sample collection record includes d2 profile sampling slices, and obtain the total number M of all target changes of q profile sampling slices of the p-th land sampling record, pq where 1 ≤ p ≤ d1 and 1 ≤ q ≤ d2;

[0021] Step S302: Calculate the sampling decision threshold γ. ;

[0022] γ represents the average number of target changes in each profile sampling slice in historical data, and represents the average number of sampled soil layers in historical records, serving as a reference standard for the number of sampled soil layers in the current sampling process.

[0023] Further, step S400 includes:

[0024] Step S401: Obtain the j-th unit sampling layer L in the target slice j , and obtain the corresponding values of the contents of various soil components in L j . Denote the corresponding value of the i-th soil component in L j as V ij ;

[0025] Step S402: Obtain the corresponding values of the i-th soil component of all the j-th unit sampling layers in all profile sampling slices, and gather the corresponding values into the sampling set E ij , and obtain the maximum value V ij and the minimum value V max in E, and calculate the coefficient of variation ω min corresponding to the i-th soil component of the j-th unit sampling layer, ij ,

[0026] ω ij =lg(V max / V min ), where lg represents the logarithmic function with base 10;

[0027] The coefficient of variation represents the sampling coverage range in the direction parallel to the ground. When the distance between profile sampling slices is relatively close, the values of V max and V min are relatively close, and the value of V max / V min tends to 1. When the coverage range is relatively wide, there are obvious differences between V max and V min . Through the coefficient of variation, the coverage degree of profile sampling slices in the horizontal direction is measured;

[0028] Step S403: Obtain h difference coefficients from several soil sample collections, and take the minimum value among the h difference coefficients as the difference reference value, denoted as B min 。

[0029] Further, step S500 includes:

[0030] Step S501: Obtain the current sampling of the garden soil, denoted as the current sampling. Denote all the profile sampling slices of the current sampling as the current sampling slices. Obtain all the soil component change functions of the current sampling slices, and obtain the number of target changes in the current sampling slices;

[0031] Step S502: Obtain the number of target changes in all the profile sampling slices of the current sampling, calculate the total number of target changes T, and calculate the average value R of the target changes of the current sampling. R = T / c, where c represents the number of all the profile sampling slices in the current sampling;

[0032] Step S503: Obtain the difference coefficients corresponding to the unit sampling layers of each profile sampling slice in the current sampling, calculate the average value of the difference coefficients, and denote the average value as a;

[0033] Step S504: Calculate the weight coefficient α, α = a / B min , when α×R<γ, give an information prompt to the relevant sampling personnel.

[0034] To better implement the above method, a soil profile sample analysis system is also proposed. The system includes:

[0035] A hyperspectral image recognition module, a historical record management module, a sampling decision threshold management module, a change range management module, and a sampling evaluation module. Among them, the hyperspectral image recognition module is used to collect the hyperspectral images of the profile sampling slices of the soil, and identify the component content in the soil. The historical record management module is used to obtain the historical sampling records of the soil and manage the historical data in the historical sampling records. The sampling decision threshold management module is used to collect the historical sampling records and calculate the sampling decision threshold. The change range management module is used to manage the changes in the soil components in the same layer of different profile sampling slices in the historical sampling records. The sampling evaluation module is used to obtain the current sampling records of the soil, evaluate the sampling records, and give an information reminder according to the current sampling records of the soil that meet the conditions;

[0036] Further, the hyperspectral image recognition module includes: a collection template management unit, a component recognition unit, and a change function management unit. Among them, the collection template management unit is used to manage the profile collection template. The component recognition unit is used to identify the content values of various soil components in the profile sampling slices. The change function management unit is used to manage the soil component change functions in the profile sampling slices;

[0037] Furthermore, the historical record management module includes: a threshold management unit, a target change recognition unit, and a change statistics unit. Among them, the threshold management unit is used to manage the change rate threshold of target changes, the target change recognition unit is used to identify function changes greater than the change rate threshold, and the change statistics unit is used to record the quantity of target changes;

[0038] Furthermore, the sampling determination threshold management module includes: a data aggregation unit and a sampling determination threshold calculation unit. Among them, the data aggregation unit is used to aggregate the quantities of all target changes in the historical records, and the sampling determination threshold calculation unit is used to calculate the sampling determination threshold;

[0039] Furthermore, the change range management module includes: a sampling set management unit, a difference coefficient calculation unit, and a difference reference value management unit. Among them, the sampling set management unit is used to aggregate the data of soil components in the unit sampling layer into the sampling set, the difference coefficient calculation unit is used to calculate the difference coefficient, and the difference reference value management unit is used to obtain the difference reference value;

[0040] Furthermore, the sampling evaluation module includes: a sampling management unit, a target change calculation unit, a weight coefficient calculation unit, and an information reminder unit. Among them, the sampling management unit is used to manage the current sampling records of the soil, the target change calculation unit is used to calculate the total quantity of target changes in the sampling records, the weight coefficient calculation unit is used to calculate the weight coefficient, and the information reminder unit is used to give an information reminder for the difference evaluation values that meet the judgment conditions.

[0041] Compared with the prior art, the beneficial effects of the present invention are as follows: By analyzing the characteristics in the vertical and horizontal directions of the soil sampling section, classifying and sorting the sampling data of the soil, a standard for comparing the current sampling data is established. By comparing the sampling data of the current soil with the historical data, analyzing from two dimensions of the vertical and horizontal directions, the weight value and the difference value are obtained respectively. By comparing the comparison model with the reference threshold, it is verified whether the sampling data of the current soil fully covers the garden. When the sampling data is insufficient to cover the garden, relevant personnel are reminded. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] Figure 1 It is a schematic structural diagram of a soil profile sample analysis system of the present invention;

[0043] Figure 2 It is a schematic flowchart of a soil profile sample analysis method of the present invention;

[0044] Figure 3 It is a schematic example diagram of a soil profile sample analysis method of the present invention. DETAILED DESCRIPTION OF THE INVENTION

[0045] The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all of the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0046] Embodiment: As Figure 1 、 Figure 2 and Figure 3 shown, the present invention provides a technical solution, a method for analyzing soil profile samples:

[0047] Step S100: Obtain a profile sampling template of the soil, collect sampling records through the profile sampling template to obtain a profile sampling slice of the soil, and identify the hyperspectral image of the soil to obtain a variation function of the corresponding values of the soil component content on the profile sampling slice;

[0048] The hyperspectral image is obtained by imaging the target area by a hyperspectral sensor on multiple consecutive and subdivided spectral bands. These bands cover the ultraviolet, visible, near-infrared, and mid-infrared regions of the electromagnetic spectrum and can capture the spectral characteristics of the soil sample at different bands. These characteristics contain information about soil composition, structure, water content, etc.;

[0049] Among them, step S100 includes:

[0050] Step S101: Obtain a profile sampling template with a depth of W1 and a width of W2, satisfying the condition W1 > W2. The profile sampling template includes n unit sampling layers, save the sampling samples of the soil through the profile sampling template to obtain a profile sampling slice of the soil, and save all the profile sampling slices in the same direction, and the direction is from the ground to the underground or from the underground to the ground;

[0051] Step S102: Identify k soil components of the profile sampling slice respectively to obtain a soil component variation function. The variation function of the i-th soil component among the k soil components is denoted as Gi(λ), where λ is the distance variable along the direction of the profile sampling slice;

[0052] In the embodiment, the sampling template represents a standardized soil sampling slice. During the implementation process, a sampling template with the same depth is set, and the depth of the sampling template is greater than the sampling depth. The non-soil layer part is removed from the sampling template to obtain a sampling slice;

[0053] Select the depth of a certain sampling template penetrating into the soil as the reference depth. All sampling templates penetrate to the horizontal plane where the reference depth is located, and the sampling slices are aligned at the bottom to reduce the influence of different terrains on the sampling slices.

[0054] Step S200: Obtain the historical sampling records of the soil. Denote a certain sampling record in the historical sampling records as the target sampling. Obtain all the profile sampling slices of the target sampling. Select a certain profile sampling slice among all the profile sampling slices as the target slice. For the variation functions corresponding to all types of soil components in the target slice, take the variation rate greater than the threshold in the variation functions as the target variation, and calculate the number of target variations in all the variation functions.

[0055] In the schematic diagram as Figure 3 shown, the depth is W1’, the width is W2’, and the profile sampling slice includes several unit sampling layers.

[0056] Among them, step S200 includes:

[0057] Step S201: Obtain the variation function Gi(λ) of the i-th soil component in the target slice, and obtain the derivative function gi(λ) of Gi(λ) in the direction of the sampling slice. ’ i(λ), and obtain Gi ’ i(λ)’s derivative function gi ’ i(λ) in the direction of the sampling slice.

[0058] Step S202: Set the variation rate threshold. Denote a continuous part greater than f in gi(λ) as a target variation, and obtain the number of target variations in gi(λ). ’ i(λ) as a target variation, and obtain the number of target variations in gi ’ i(λ).

[0059] Step S203: Respectively obtain the soil component variation functions corresponding to k soil components in the target slice, obtain the target variations corresponding to the soil component variation functions, and record the total number of all target variations in the target slice as m.

[0060] Step S300: Obtain all the profile sampling slices of several land sample collections, calculate the mean value of the number of all target variations, and take the mean value as the sampling determination threshold.

[0061] Among them, step S300 includes:

[0062] Step S301: Obtain d1 land sample collection records. Among them, each land sample collection record includes d2 profile sampling slices. Obtain the total number M of all target variations of the q profile sampling slices of the p-th land sampling record pq where 1 ≤ p ≤ d1 and 1 ≤ q ≤ d2;

[0063] Step S302: Calculate the sampling determination threshold γ. 。

[0064] Step S400: Obtain all profile sampling slices in the target sampling, obtain the corresponding values of the soil component contents in all unit sampling layers of each profile sampling slice, and calculate the change ranges of various soil component contents respectively;

[0065] Among them, step S400 includes:

[0066] Step S401: Obtain the j-th unit sampling layer L in the target slice j , obtain the corresponding values of various soil component contents in L j , and denote the corresponding value of the i-th soil component in L j as V ij ;

[0067] In this solution, W2 is the short side of the profile sampling slice. In the short side direction, in each profile sampling slice, the corresponding values of the component contents change little and are approximately taken as a fixed value. In the unit sampling layer, the average value of the corresponding values of the component contents in the unit sampling layer is taken as the corresponding value in the corresponding unit sampling layer;

[0068] Step S402: Obtain the corresponding values of the i-th soil component of all the j-th unit sampling layers in all profile sampling slices, collect the corresponding values into the sampling set E ij , obtain the maximum value V ij and the minimum value V max in E min , and calculate the difference coefficient ω ij corresponding to the i-th soil component of the j-th unit sampling layer,

[0069] ω ij = lg(V max / V min ), where lg represents the logarithmic function with base 10;

[0070] Step S403: Obtain h difference coefficients from several land sample collections, take the minimum value among the h difference coefficients as the difference reference value, and denote the difference reference value as B min ;

[0071] In the embodiment, in the same sampling record, the corresponding values of the soil component contents in each unit sampling layer are respectively collected, and the difference coefficients of the corresponding values of each soil component content in the same layer are calculated;

[0072] Collect the difference coefficients of each layer of each sampling record. At this time, the total number of difference coefficients is h, and the minimum value among them is selected to obtain B min ;

[0073] Alternatively, randomly sample h coefficient of variation values and take the minimum value among them to obtain B min .

[0074] Step S500: Obtain the current sampling record of the soil. Calculate the weight coefficient by comparing the difference in the change range of soil component content between the current sampling record and the historical sampling record. Compare the difference in the target change between the current sampling record and the historical sampling record. Measure the difference through the weight coefficient to obtain a difference evaluation value. When the difference evaluation value is less than the sampling determination threshold, give an information prompt;

[0075] Among them, step S500 includes:

[0076] Step S501: Obtain the current sampling of the garden land as the current sampling. Denote all the profile sampling slices of the current sampling as the current sampling slices. Obtain all the soil component change functions of the current sampling slices. Obtain the number of target changes in the current sampling slices;

[0077] Step S502: Obtain the number of target changes in all the profile sampling slices in the current sampling. Calculate the total number of target changes T. Calculate the average value R of the target changes in the current sampling. R = T / c, where c represents the number of all the profile sampling slices in the current sampling;

[0078] Step S503: Obtain the coefficient of variation corresponding to the unit sampling layer of each profile sampling slice in the current sampling. Calculate the average value of the coefficient of variation and denote the average value as a;

[0079] Step S504: Calculate the weight coefficient α, α = a / B min , when α×R<γ, give an information prompt to the relevant sampling personnel.

[0080] A soil profile sample analysis system, the system includes: a hyperspectral image recognition module, a historical record management module, a sampling determination threshold management module, a change range management module, and a sampling evaluation module;

[0081] Among them, the hyperspectral image recognition module is used to collect hyperspectral images of the profile sampling slices of the soil and identify the component content in the soil. Among them, the hyperspectral image recognition module includes: a collection template management unit, a component recognition unit, and a change function management unit. Among them, the collection template management unit is used to manage the profile collection template, the component recognition unit is used to identify the content values of various components in the soil of the profile sampling slices, and the change function management unit is used to manage the soil component change functions in the profile sampling slices.

[0082] Among them, the historical record management module is used to obtain the historical sampling records of the soil and manage the historical data in the historical sampling records. Among them, the historical record management module includes: a threshold management unit, a target change identification unit, and a change statistics unit. Among them, the threshold management unit is used to manage the change rate threshold of the target change, the target change identification unit is used to identify the function changes greater than the change rate threshold, and the change statistics unit is used to record the number of target changes;

[0083] Among them, the sampling determination threshold management module is used to collect the historical sampling records and calculate the sampling determination threshold. Among them, the sampling determination threshold management module includes: a data collection unit and a sampling determination threshold calculation unit. Among them, the data collection unit is used to collect the number of all target changes in the historical records, and the sampling determination threshold calculation unit is used to calculate the sampling determination threshold;

[0084] Among them, the change range management module is used to manage the changes in the soil components in the same layer of different profile sampling slices in the historical sampling records. Among them, the change range management module includes: a sampling set management unit, a difference coefficient calculation unit, and a difference reference value management unit. Among them, the sampling set management unit is used to collect the data of the soil components in the unit sampling layer and record them in the sampling set, the difference coefficient calculation unit is used to calculate the difference coefficient, and the difference reference value management unit is used to obtain the difference reference value;

[0085] Among them, the sampling evaluation module is used to obtain the current sampling records of the soil, evaluate the sampling records, and give an information reminder according to the current sampling records of the soil that meet the conditions. Among them, the sampling evaluation module includes: a sampling management unit, a target change calculation unit, a weight coefficient calculation unit, and an information reminder unit. Among them, the sampling management unit is used to manage the current sampling records of the soil, the target change calculation unit is used to calculate the total number of target changes in the sampling records, the weight coefficient calculation unit is used to calculate the weight coefficient, and the information reminder unit is used to give an information reminder for the difference evaluation value that meets the judgment conditions.

[0086] For those skilled in the art, it is obvious that the present invention is not limited to the details of the above exemplary embodiments, and without departing from the spirit or basic characteristics of the present invention, the present invention can be implemented in other specific forms. Therefore, from any point of view, the embodiments should be regarded as exemplary and non-limiting. The scope of the present invention is defined by the appended claims rather than the above description. Therefore, all changes falling within the meaning and scope of the equivalent elements of the claims are intended to be included in the present invention. Any reference signs in the claims should not be regarded as limiting the claimed rights.

Claims

1. A method for analyzing soil profile samples, characterized in that: The method includes the following steps: Step S100: Obtain a profile sampling template of the soil, collect sampling records through the profile sampling template to obtain a profile sampling slice of the soil, and identify the hyperspectral image of the soil to obtain the variation function of the corresponding values of the soil component content on the profile sampling slice; Step S200: Obtain the historical sampling records of the soil, record a certain sampling record in the historical sampling records as the target sampling, obtain all the profile sampling slices of the target sampling, take a certain profile sampling slice among all the profile sampling slices as the target slice, for the variation functions corresponding to all types of soil components in the target slice, take the variation rate greater than the threshold in the variation function as the target variation, and calculate the number of target variations in all the variation functions; Step S300: Obtain all the profile sampling slices of several land sample collections, calculate the mean value of the number of all target variations, and take the mean value as the sampling determination threshold; Step S300 includes: Step S301: Obtain the d1 land sample collection records. Each land sample collection record includes d2 profile sampling slices, and obtain the total number M of all target changes of the q profile sampling slices of the p-th land sampling record pq , where 1 ≤ p ≤ d1 and 1 ≤ q ≤ d2; Step S302: Calculate the sampling decision threshold γ, ; Step S400: Obtain all the profile sampling slices in the target sampling, obtain the corresponding values of the soil component content in all the unit sampling layers of each profile sampling slice, and calculate the variation range of each soil component content respectively; Step S400 includes: Step S401: Obtain the j-th unit sampling layer L in the target slice j , and obtain L j The corresponding values of the contents of various soil components in it. Denote the corresponding value of the i-th soil component in L j as V ij ; Step S402: Obtain the corresponding values of the i-th soil component in all the j-th unit sampling layers among all the profile sampling slices, and collect the corresponding values into the sampling set E ij , obtain E ij the maximum value V max and the minimum value V min , calculate the coefficient of variation ω corresponding to the i-th soil component in the j-th unit sampling layer ij , ω ij =lg(V max / V min ), where lg represents the logarithmic function with base 10; Step S403: Obtain h difference coefficients from several soil sample collections, and take the minimum value among the h difference coefficients as the difference reference value, denoted as B min ; Step S500: Obtain the current sampling record of the soil, calculate the weight coefficient by comparing the difference in the variation range of the soil component content between the current sampling record and the historical sampling record, compare the difference in the target variation between the current sampling record and the historical sampling record, measure the difference through the weight coefficient to obtain a difference evaluation value, and perform an information prompt when the difference evaluation value is less than the sampling determination threshold; Step S500 includes: Step S501: Obtain the current sampling of the garden land as the current sampling, record all the profile sampling slices of the current sampling as the current sampling slices, obtain all the soil component variation functions of the current sampling slices, and obtain the number of target variations in the current sampling slices; Step S502: Obtain the number of target variations in all the profile sampling slices in the current sampling, calculate the total number of target variations T, calculate the average value R of the target variations of the current sampling, R = T / c, where c represents the number of all the profile sampling slices in the current sampling; Step S503: Obtain the difference coefficient corresponding to the unit sampling layer of each profile sampling slice in the current sampling, calculate the average value of the difference coefficients, and record the average value as a; Step S504: Calculate the weight coefficient α, where α = a / B min , and when α × R < γ, give an information prompt to the relevant sampling personnel.

2. The soil profile sample analysis method according to claim 1, characterized in that: Step S100 includes: Step S101: Obtain a profile sampling template with a depth of W1 and a width of W2, satisfying the condition W1 > W2. The profile sampling template includes n unit sampling layers. Save the sampling samples of the soil through the profile sampling template to obtain a profile sampling slice of the soil, and save all the profile sampling slices in the same direction, and the direction is from the ground to the underground or from the underground to the ground; Step S102: Identify k soil components of the profile sampling slice respectively to obtain the soil component variation function. The variation function of the i-th soil component among the k soil components is denoted as Gi(λ), where λ is the distance variable along the direction of the profile sampling slice.

3. The soil profile sample analysis method according to claim 2, wherein: Step S200 includes: Step S201: Obtain the change function G ’ i(λ) of the i-th soil component in the target slice, and obtain G ’ The derivative function g ’ i(λ) of i(λ) in the direction of the sampling slice; Step S202: Set the change rate threshold f, and record a part of g ’ where one continuous part greater than f in i(λ) as a target change, and obtain the number of target changes in g ’ in i(λ); Step S203: Obtain the soil component change functions corresponding to k types of soil components in the target slice respectively, obtain the target changes corresponding to each soil component change function, and record the total number of all target changes of the target slice as m.

4. A soil profile sample analysis system for performing the soil profile sample analysis method according to any one of claims 1-3, characterized in that: The system includes: A hyperspectral image recognition module, a historical record management module, a sampling decision threshold management module, a change range management module, and a sampling evaluation module. Among them, the hyperspectral image recognition module is used to collect hyperspectral images of the profile sampling slices of the soil and identify the component contents in the soil. The historical record management module is used to obtain the historical sampling records of the soil and manage the historical data in the historical sampling records. The sampling decision threshold management module is used to collect the historical sampling records and calculate the sampling decision threshold. The change range management module is used to manage the changes of soil components in the same layer of different profile sampling slices in the historical sampling records. The sampling evaluation module is used to obtain the current sampling records of the soil, evaluate the sampling records, and give an information reminder according to the current sampling records of the soil that meet the conditions.

5. A soil profile sample analysis system according to claim 4, characterized in that: The hyperspectral image recognition module includes: a collection template management unit, a component recognition unit, and a change function management unit. Among them, the collection template management unit is used to manage the profile collection template. The component recognition unit is used to identify the content values of various soil components in the profile sampling slice. The change function management unit is used to manage the soil component change functions in the profile sampling slice.

6. The soil profile sample analysis system according to claim 4, characterized in that: The historical record management module includes: a threshold management unit, a target change recognition unit, and a change statistics unit. Among them, the threshold management unit is used to manage the change rate threshold of the target change. The target change recognition unit is used to identify the function changes greater than the change rate threshold. The change statistics unit is used to record the number of target changes. The sampling decision threshold management module includes: a data collection unit and a sampling decision threshold calculation unit. Among them, the data collection unit is used to collect the number of all target changes in the historical records. The sampling decision threshold calculation unit is used to calculate the sampling decision threshold. The change range management module includes: a sampling set management unit, a difference coefficient calculation unit, and a difference reference value management unit. Among them, the sampling set management unit is used to collect the data of soil components in the unit sampling layer and record them in the sampling set. The difference coefficient calculation unit is used to calculate the difference coefficient. The difference reference value management unit is used to obtain the difference reference value.

7. A soil profile sample analysis system according to claim 4, characterized in that: The sampling evaluation module includes: a sampling management unit, a target change calculation unit, a weight coefficient calculation unit, and an information reminder unit. Among them, the sampling management unit is used to manage the current sampling records of the soil. The target change calculation unit is used to calculate the total number of target changes in the sampling records. The weight coefficient calculation unit is used to calculate the weight coefficient. The information reminder unit is used to give an information reminder for the difference evaluation value that meets the judgment conditions.

Images