CN116433408A

CN116433408A - Method and device for determining phase change critical point of tilling and utilizing intensity

Info

Publication number: CN116433408A
Application number: CN202310680381.1A
Authority: CN
Inventors: 宋长青; 叶思菁; 蒋嘉益; 高培超; 穆望舒
Original assignee: Beijing Normal University
Current assignee: Beijing Normal University
Priority date: 2023-06-09
Filing date: 2023-06-09
Publication date: 2023-07-14

Abstract

The invention provides a method and a device for determining a phase change critical point of farmland utilization intensity, and relates to the field of farmland intensive utilization research, wherein the method for determining the phase change critical point of farmland utilization intensity comprises the following steps: acquiring farmland input intensity of a plurality of farmland in a designated input dimension and farmland output intensity in a designated output dimension; calculating a correlation coefficient and a significance coefficient corresponding to each cultivated land input intensity range according to the cultivated land input intensity and cultivated land output intensity of each cultivated land in the cultivated land input intensity range aiming at each preset cultivated land input intensity range; and determining the phase change critical point of the tilling and utilizing intensity under the appointed input dimension and the appointed output dimension according to the correlation coefficient and the significance coefficient corresponding to the input intensity range of each tilling area. The dynamic expression of the intensity phase change process and the phase change critical point of the cultivated land is realized, the basis can be provided for diagnosing the phase change inflection point of the input-output coupling effect of the cultivated land, and the method has the advantages of clear expression and easiness in popularization and application.

Description

Method and device for determining phase change critical point of tilling and utilizing intensity

Technical Field

The invention relates to the field of intensive utilization research of cultivated land, in particular to a method and a device for determining a phase change critical point of cultivated land utilization intensity.

Background

Today, global food security faces a great challenge. The phase change critical point of the farmland utilization intensity is evaluated on the national scale and the global scale, and guidance can be provided for adjusting the farmland utilization intensity and the farmland utilization structure so as to better cope with the conflict among grain requirements, economic development and ecological protection.

In most tilling use strength assessment frameworks, the relationship of input strength and output strength is considered an important dimension. However, in a specific calculation, this relationship is expressed in terms of an efficiency index and is typically calculated by the ratio of the output intensity to the input intensity. This ignores the complex nonlinear relationship between input and output. In addition, the higher the tilling strength is, the better, and the higher the tilling strength is, the burden of the farmland ecosystem is increased, and the sustainability of the tilling is threatened. Therefore, how to quantitatively describe the impact of input intensity on output intensity, and thus estimate the phase transition critical point of the tilling intensity remains a challenge.

Disclosure of Invention

Aiming at the problems existing in the prior art, the embodiment of the invention provides a method and a device for determining the phase change critical point of the tillage utilization intensity.

The invention provides a method for determining a phase change critical point of a tilling strength, which comprises the following steps:

acquiring farmland input intensity of a plurality of farmland in a designated input dimension and farmland output intensity in a designated output dimension, wherein the designated input dimension is any input dimension, and the designated output dimension is any output dimension;

calculating a correlation coefficient and a significance coefficient corresponding to each preset cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range;

and determining the phase change critical point of the tilling and utilizing intensity in the appointed input dimension and the appointed output dimension according to the correlation coefficient and the significance coefficient corresponding to each tilling and throwing intensity range.

According to the method for determining the phase transition critical point of the cultivated land utilization intensity provided by the invention, before calculating the correlation coefficient and the significance coefficient corresponding to each cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range, the method further comprises:

Acquiring environmental characteristics of a plurality of cultivated lands, wherein the environmental characteristics comprise at least one of climate conditions, topography characteristics, soil properties, cultivated land facility management and cultivated land crushing degree;

dividing the plurality of cultivated lands into a plurality of cultivated land types according to the environmental characteristics;

correspondingly, for each preset cultivated land input intensity range, calculating a correlation coefficient and a significance coefficient corresponding to the cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range, including:

calculating a correlation coefficient and a significance coefficient corresponding to each preset cultivated land input intensity range according to the cultivated land input intensity and the cultivated land generation intensity of each cultivated land of a designated cultivated land type in the cultivated land input intensity range, wherein the designated cultivated land type is any cultivated land type;

determining the phase change critical point of the tilling strength in the specified input dimension and the specified output dimension according to the correlation coefficient and the significance coefficient corresponding to each tilling input strength range comprises the following steps:

And determining the phase change critical point of the tilling utilization intensity of the designated tilling type under the designated input dimension and the designated output dimension according to the correlation coefficient and the significance coefficient corresponding to each tilling input intensity range.

According to the method for determining the phase transition critical point of the cultivated land utilization intensity provided by the invention, for each preset cultivated land input intensity range, according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range, a correlation coefficient and a significance coefficient corresponding to the cultivated land input intensity range are calculated, and the method comprises the following steps:

constructing a two-dimensional coordinate graph according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land, and setting the window width and the sliding step length of a sliding window, wherein the transverse direction of the two-dimensional coordinate graph is the cultivated land input intensity, the longitudinal direction of the two-dimensional coordinate graph is the cultivated land output intensity, and the sliding window slides along the longitudinal direction;

and moving the sliding window backwards from the longitudinal starting position according to the sliding step length, calculating a correlation coefficient and a significance coefficient corresponding to the window width for the cultivated land input intensity and the cultivated land output intensity of each cultivated land in each window width, wherein each window width corresponds to one cultivated land input intensity range.

According to the method for determining the phase change critical point of the tilling strength provided by the invention, the determining the phase change critical point of the tilling strength in the appointed input dimension and the appointed output dimension according to the correlation coefficient and the significance coefficient corresponding to each tilling input strength range comprises the following steps:

for each cultivated land input intensity range, drawing points in a two-dimensional coordinate system according to a specified display mode by taking the average value of the cultivated land input intensity ranges as an abscissa and the correlation coefficient corresponding to the cultivated land input intensity ranges as an ordinate, wherein the specified display mode is a display method corresponding to a specified significance region, and the specified significance region is a significance region where the significance coefficient corresponding to the cultivated land input intensity ranges is located;

and determining the phase change critical points of the tilling and utilizing intensity under the specified input dimension and the specified output dimension according to each point in the two-dimensional coordinate system.

According to the method for determining the phase change critical point of the tilling strength provided by the invention, the phase change critical point of the tilling strength under the specified input dimension and the specified output dimension is determined according to each point in the two-dimensional coordinate system, and the method comprises the following steps:

Identifying a positive correlation stage, a non-correlation stage and a negative correlation stage between the cultivated land input intensity and the cultivated land output intensity according to each point in the two-dimensional coordinate system;

determining a first boundary cultivated land input intensity connecting the positive correlation stage and the non-correlation stage, and taking a first initial value of a cultivated land input intensity range where the first boundary cultivated land input intensity is located as a first phase change critical point of cultivated land utilization intensity in the specified input dimension and the specified output dimension;

determining a second boundary farmland input intensity connecting the non-correlation stage and the negative correlation stage, and taking a second initial value of a farmland input intensity range where the second boundary farmland input intensity is located as a second phase change critical point of the farmland utilization intensity in the specified input dimension and the specified output dimension.

According to the method for determining the phase change critical point of the tilling and utilizing intensity provided by the invention, after acquiring the tilling input intensity of a plurality of tilling in a designated input dimension and the tilling output intensity in a designated output dimension, the method further comprises the following steps:

and normalizing the cultivated land input intensity and/or the cultivated land generation intensity.

According to the method for determining the phase change critical point of the tilling and utilizing intensity, the correlation coefficient is a partial correlation coefficient;

correspondingly, for each preset cultivated land input intensity range, calculating a correlation coefficient corresponding to the cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range, including:

and aiming at each preset cultivated land input intensity range, taking the cultivated land input intensity in a non-designated input dimension as a control variable, and calculating a correlation coefficient corresponding to the cultivated land input intensity range according to the cultivated land input intensity of each cultivated land in the designated input dimension and the cultivated land output intensity in the designated output dimension in the cultivated land input intensity range.

The invention also provides a device for determining the phase change critical point of the tillage utilization intensity, which comprises:

the system comprises an acquisition module, a control module and a control module, wherein the acquisition module is configured to acquire farmland input intensity of a plurality of farmland in a designated input dimension and farmland output intensity in a designated output dimension, wherein the designated input dimension is any input dimension, and the designated output dimension is any output dimension;

A calculation module configured to calculate, for each preset cultivated land input intensity range, a correlation coefficient and a significance coefficient corresponding to the cultivated land input intensity range from the cultivated land input intensity and the cultivated land output intensity of each cultivated land within the cultivated land input intensity range;

and the determining module is configured to determine a phase change critical point of the tilling utilization intensity in the appointed input dimension and the appointed output dimension according to the correlation coefficient and the significance coefficient corresponding to each tilling input intensity range.

The invention also provides an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing a method of determining a phase change critical point of a tilling strength as described in any of the above when executing the program.

The present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements a method of determining a phase change critical point of a tilling use intensity as described in any of the above.

The invention also provides a computer program product comprising a computer program which when executed by a processor implements a method of determining a phase transition critical point of a tilling strength as described in any of the above.

According to the method and the device for determining the phase change critical point of the tilling and utilizing strength, the tilling input strength of a plurality of tilling areas in the appointed input dimension and the tilling output strength of the tilling areas in the appointed output dimension are obtained, wherein the appointed input dimension is any input dimension, and the appointed output dimension is any output dimension; calculating a correlation coefficient and a significance coefficient corresponding to each preset cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range; and determining the phase change critical point of the tilling and utilizing intensity in the appointed input dimension and the appointed output dimension according to the correlation coefficient and the significance coefficient corresponding to each tilling and throwing intensity range. By calculating the promotion or inhibition effect of the farmland input intensity change on the farmland output intensity in different farmland input intensity ranges, the dynamic expression of the farmland utilization intensity phase change process and the phase change critical point is realized, the basis can be provided for diagnosing the phase change inflection point of the farmland input-output coupling effect, and the method has the advantages of clear expression and easiness in popularization and application.

Drawings

In order to more clearly illustrate the invention or the technical solutions of the prior art, the following description will briefly explain the drawings used in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are some embodiments of the invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.

FIG. 1 is a flow chart of a method for determining a phase transition critical point of a tilling strength according to the present invention;

FIG. 2 is one of the two-dimensional graphs provided by the present invention;

FIG. 3 is a second two-dimensional graph provided by the present invention;

FIG. 4 is a schematic flow chart of calculating correlation coefficients and significance coefficients according to the present invention;

FIG. 5 is a schematic diagram showing the determination of the phase change critical point according to the present invention

FIG. 6 is a second schematic diagram of determining a phase transition threshold according to the present invention;

FIG. 7 is a schematic structural view of an apparatus for determining the phase transition critical point of the tilling strength according to the present invention;

fig. 8 is a schematic structural diagram of an electronic device provided by the present invention.

Detailed Description

For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.

In order to facilitate a clearer understanding of various embodiments of the present invention, some relevant background knowledge is first presented as follows.

Today, global food security faces a great challenge. Efforts to develop agriculture to meet food demand have also required balancing the health and stability of the ecosystem, a fact that has made food safety challenges more complex. The sustainable agriculture development road is explored, on the premise of meeting the increasing grain demands, the harm of agriculture to the environment is reduced, and the fortune of all human beings is related. Intensive utilization of cultivated land is widely regarded as an important approach for relieving contradiction between grain gaps and ecological system health, and its contribution to limiting cultivated land expansion and protecting biodiversity has been confirmed by many studies. The phase change critical point of the farmland utilization intensity is evaluated on the national scale and the global scale, and guidance can be provided for adjusting the farmland utilization intensity and the farmland utilization structure so as to better cope with the conflict among grain requirements, economic development and ecological protection.

In most tilling use strength assessment frameworks, the relationship of input strength and output strength is considered an important dimension. However, in a specific calculation, the relationship is expressed by an efficiency index, and is generally calculated by a ratio of output intensity to input intensity, for example, some methods calculate economic benefits of cultivated land based on a ratio of output value of grain to input cost of cultivation, and other methods calculate nitrogen utilization rate based on a ratio of output nitrogen to input nitrogen, etc. This ignores the complex nonlinear relationship between input and output. For example, a highly efficient cultivation under extensive management may be obtained due to the extremely low investment. In addition, the higher the tilling strength is, the better, and the higher the tilling strength is, the burden of the farmland ecosystem is increased, and the sustainability of the tilling is threatened.

Therefore, to quantitatively describe the impact of input intensity on output intensity, the phase transition critical point of the tilling intensity is estimated. The invention provides a method and a device for determining a phase change critical point of tilling and utilizing strength, which are characterized in that tilling input strength of a plurality of tilling in a designated input dimension and tilling output strength of the tilling in a designated output dimension are obtained, wherein the designated input dimension is any input dimension, and the designated output dimension is any output dimension; calculating a correlation coefficient and a significance coefficient corresponding to each preset cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range; and determining the phase change critical point of the tilling and utilizing intensity in the appointed input dimension and the appointed output dimension according to the correlation coefficient and the significance coefficient corresponding to each tilling and throwing intensity range. By calculating the promotion or inhibition effect of the farmland input intensity change on the farmland output intensity in different farmland input intensity ranges, the dynamic expression of the farmland utilization intensity phase change process and the phase change critical point is realized, the basis can be provided for diagnosing the phase change inflection point of the farmland input-output coupling effect, and the method has the advantages of clear expression and easiness in popularization and application.

The method and apparatus for determining the phase transition critical point of the tilling strength of the present invention are described below with reference to fig. 1-7.

Fig. 1 is a schematic flow chart of a method for determining a phase transition critical point of a tilling strength according to the present invention, referring to fig. 1, including steps 101 to 103, wherein:

step 101: the method comprises the steps of obtaining farmland input intensity of a plurality of farmland in a designated input dimension and farmland output intensity in a designated output dimension, wherein the designated input dimension is any input dimension, and the designated output dimension is any output dimension.

It should be noted that the execution body of the present invention may be any electronic device for determining the phase change critical point, for example, any one of a smart phone, a smart watch, a desktop computer, a laptop computer, and the like.

Specifically, the input dimension, namely the type of the cultivated land input intensity, comprises a substance type input dimension and a yield increase type input dimension, wherein the substance type input dimension comprises a chemical fertilizer input dimension, a pesticide input dimension, an agricultural film input dimension and the like, and the yield increase type input dimension comprises a labor input dimension, a diesel input dimension, an agricultural machine input dimension and the like. The cultivated land input intensity can be expressed by means of the energy value input per unit area, the amount of matter input per unit area, and the like, for example, the cultivated land input intensity in the fertilizer input dimension, the pesticide input dimension, the labor input dimension, the diesel input dimension, and the agricultural machine input dimension are the chemical input amount per unit area (fertilizer input intensity), the pesticide input amount per unit area (pesticide input intensity), the labor input amount per unit area (labor input intensity), the diesel input amount per unit area (diesel input intensity), and the total power per unit area (agricultural machine input intensity), respectively.

The yield dimension, namely the type of cultivated yield strength, can be divided into yield dimensions of different types of crops, such as corn yield dimension, wheat yield dimension, rice yield dimension and the like, and also can be divided into main grain crop yield dimension, oil crop yield dimension, sugar crop yield dimension, vegetable yield dimension and the like. The cultivated land yield strength can be expressed by means of energy yield per unit area, value yield per unit area, quality of yield per unit area, and the like, for example, the cultivated land yield strength in the corn yield dimension, the wheat yield dimension, and the rice yield dimension are corn yield per unit area (corn yield strength), wheat yield per unit area (wheat yield strength), and rice yield per unit area (rice yield strength), respectively; for example, the yield of the main grain crop, the yield of the oil crop, the yield of the sugar crop and the yield of the vegetable in the yield of the vegetable are respectively the yield of the main grain crop (yield of the main grain crop), the yield of the oil crop (yield of the oil crop), the yield of the sugar crop (yield of the sugar crop) and the yield of the vegetable (yield of the vegetable). Different types of cultivated land input intensity dimensions are inconsistent and need to be converted into uniform dimensions. Different types of cultivated land yield strength have different crop types and need to be converted into uniform energy yield per unit area or standard grain yield per unit area.

In practical application, after a plurality of cultivated lands of a research area are determined, cultivated land input intensities of the cultivated lands in each input dimension and cultivated land output intensities of the cultivated lands in each output dimension are obtained. There are various methods for acquiring the input intensities and the output intensities of the cultivated lands, for example, a user uploads the input intensities of the cultivated lands and the output intensities of the cultivated lands in each input dimension respectively through an uploading page, and accordingly, the executing body receives the input intensities of the cultivated lands and the output intensities of the cultivated lands in each input dimension respectively; for another example, the execution body receives the acquisition instruction or the phase change critical point determination instruction, and accordingly, the execution body acquires the cultivated land input intensity of the plurality of cultivated lands in each input dimension and the cultivated land output intensity in each output dimension from the storage area to which the instruction is directed. The invention is not limited in this regard.

Alternatively, the cultivated land input intensity at each input dimension and the cultivated land output intensity at each output dimension of the plurality of cultivated lands may be collected from the internet.

Alternatively, the cultivated land input amount of the plurality of cultivated lands in each input dimension and the cultivated land output amount of the plurality of cultivated lands in each output dimension may be obtained, and then calculation is performed to obtain cultivated land input intensities of the plurality of cultivated lands in each input dimension and cultivated land output intensities of the plurality of cultivated lands in each output dimension. Namely, acquiring the cultivated land input intensity of a plurality of cultivated lands in a designated input dimension and the cultivated land output intensity in a designated output dimension, comprising:

Acquiring a sowing area of any cultivated land, an input amount of the cultivated land in a specified input dimension and an output amount of the cultivated land in a specified output dimension;

and calculating the cultivated land input intensity of the cultivated land in the appointed input dimension according to the sowing area and the cultivated land input intensity, and calculating the cultivated land output intensity of the cultivated land in the appointed output dimension according to the sowing area and the cultivated land output quantity.

Specifically, the sowing area may be the actual area of the cultivated land, or the area of the cultivated land in which different crops are sown in the year. The cultivated land input amount is the amount of the cultivated land input, such as fertilizer input amount. The cultivated land output refers to the number of crops produced by the cultivated land.

In practical application, the cultivated land input intensity and the cultivated land output intensity are calculated by acquiring the cultivated land input intensity and the cultivated land output intensity, so that the cultivated land input intensity of each input dimension and the cultivated land output intensity of each output dimension can be ensured to be acquired, and the problem that the determination of the phase change critical point cannot be normally performed due to the fact that the cultivated land input intensity and the cultivated land output intensity cannot be directly acquired is avoided.

Illustratively, calculating the cultivated land input intensity (unit: solar joules/hectare) of provincial unit years from the chemical fertilizer input dimension, the pesticide input dimension, the labor input dimension and the agricultural machinery input dimension by 4 dimensions, respectively; the cultivated land yield strength (units: joules/hectare) for the provincial unit year is calculated from the staple grain crop yield dimension.

Taking province M as an example, the chemical fertilizer input strength (ALUI) _{in_fer} ) Namely the farmland input intensity of the fertilizer input dimension, the calculation formula is shown in the formula (1); pesticide input intensity (ALUI) for that year _{in_pes} ) Namely the farmland throwing strength in the pesticide throwing dimension, and the calculation formula is shown as the formula (2); agricultural machine input intensity (ALUI) for that year _{in_am} ) Namely the farmland input intensity of the agricultural machine input dimension, and the calculation formula is shown as a formula (3); labor input intensity for that year (ALUI) _{in_lab} ) Tilling in a labor input dimensionThe intensity is shown in a formula (4); the yield strength (ALUI) of the main grain crop of that year _out ) Namely the cultivated land strength under the yield dimension of the staple food crops, and the calculation formula is shown as formula (5).

(1)

(2)

(3)

(4)

(5)

In the formulas (1) to (5), SET _FN 、SET _FP 、SET _FK And SET _FC Solar value conversion coefficients (unit: solar joules/ton) of nitrogen fertilizer, potassium fertilizer, phosphate fertilizer and compound fertilizer are respectively represented; q (Q) _FN 、Q _FP 、Q _FK And Q _FC The total input amount (unit: ton) of nitrogen fertilizer, potash fertilizer, phosphate fertilizer and compound fertilizer in the year of province M is respectively represented; s is S _i Is the sowing area (unit: hectare) of the crop i; n represents n crops; i represents an ith crop; q (Q) _pes 、Q _am And Q _lab The total input of pesticide in the year (unit: ton), the total kinetic energy of agricultural machinery (unit: kilowatt-hour) and the agricultural labor population (unit: person) are respectively the province M; SET (SET) _pes 、SET _am And SET _lab Solar value conversion coefficients respectively representing pesticides, agricultural machinery and labor force; y is Y _i Represents the annual total yield (in kg) of the ith crop, NCC _i Indicating the corresponding nutritional energy transfer of the ith cropAnd (5) coefficient conversion.

Step 102: and calculating a correlation coefficient and a significance coefficient corresponding to each cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range aiming at each preset cultivated land input intensity range.

Specifically, the correlation coefficient may be a Pearson (Pearson) correlation coefficient, a partial correlation coefficient, or other types of correlation coefficients. The significance coefficient, i.e., the significance difference (Statistical Significance), i.e., the P-value, is a statistical evaluation of the data differences.

In practical application, on the basis of determining the cultivated land input intensity and the cultivated land output intensity, determining a target cultivated land of which the cultivated land input intensity belongs to a preset cultivated land input intensity range for each cultivated land input intensity range, and calculating a correlation coefficient and a significance coefficient between a designated cultivated land input dimension and a designated cultivated land output dimension in the cultivated land input intensity range according to a set calculation strategy according to the cultivated land input intensity and the cultivated land output intensity of all target cultivated lands.

Step 103: and determining the phase change critical point of the tilling and utilizing intensity in the appointed input dimension and the appointed output dimension according to the correlation coefficient and the significance coefficient corresponding to each tilling and throwing intensity range.

Specifically, the phase transition refers to a phase transition process in which the correlation between the input intensity of cultivated land and the intensity of cultivated land undergoes a significant positive correlation-uncorrelation-significant negative correlation. The phase transition critical point refers to the turning point in the phase transition process.

In practical application, on the basis of obtaining the correlation coefficient and the significance coefficient corresponding to each cultivated land input intensity range, further, the correlation coefficient and the significance coefficient corresponding to each cultivated land input intensity range determine that the correlation between cultivated land input intensity and cultivated land output intensity is subjected to a phase change process of significant positive correlation, non-correlation and significant negative correlation, and a phase change critical point is extracted from the phase change process.

According to the method for determining the phase change critical point of the tilling and utilizing strength, the tilling input strength of a plurality of tilling areas in the appointed input dimension and the tilling output strength of the tilling areas in the appointed output dimension are obtained, wherein the appointed input dimension is any input dimension, and the appointed output dimension is any output dimension; calculating a correlation coefficient and a significance coefficient corresponding to each preset cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range; and determining the phase change critical point of the tilling and utilizing intensity in the appointed input dimension and the appointed output dimension according to the correlation coefficient and the significance coefficient corresponding to each tilling and throwing intensity range. By calculating the promotion or inhibition effect of the farmland input intensity change on the farmland output intensity in different farmland input intensity ranges, the dynamic expression of the farmland utilization intensity phase change process and the phase change critical point is realized, the basis can be provided for diagnosing the phase change inflection point of the farmland input-output coupling effect, and the method has the advantages of clear expression and easiness in popularization and application.

In one or more optional embodiments of the present invention, before calculating the correlation coefficient and the significance coefficient corresponding to the cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range for each preset cultivated land input intensity range, the method further includes:

In practical application, because of the large difference of the phase change critical points of the tilling utilization intensities of the tilling of different tilling types, in order to improve the accuracy and the reliability of the determined phase change critical points, a plurality of tilling can be divided into n groups according to climate conditions, terrain features, soil properties, tilling crushing degree features, tilling crushing degree and the like, each group corresponds to one tilling type, and n is a positive integer. Then, for each cultivated land type cultivated land, the correlation coefficient and the saliency coefficient are calculated, and the phase change critical point is determined based on the correlation coefficient and the saliency coefficient. Therefore, the corresponding relation between the cultivated land type and the phase change critical point is ensured, the phase change critical point is determined in a finer granularity, and the accuracy and the reliability of the phase change critical point are improved.

For research areas, such as a province, a country and the like, an index system (an environment characteristic system) for grouping or dividing the types of cultivated lands is designed, index data (data corresponding to environment characteristics) are collected, and weather condition indexes such as average precipitation, average air temperature, minimum air temperature, maximum air temperature and the like are included in a month by month of a certain year; taking the average elevation and the average gradient of the cultivated land area as terrain characteristic indexes; soil organic carbon density, soil volume weight, soil thickness, clay/powder/sand ratio, cation exchange capacity and the like are taken as soil property indexes; the average plaque area, the cultivated land density and the area weighted shape index are taken as cultivated land crushing degree indexes. Based on the data corresponding to the indexes, the plurality of cultivated lands are divided into a plurality of cultivated land types, namely n groups.

The multiple cultivated lands may be grouped by a subjective grouping method, for example, by an average grouping according to environmental characteristics, or by a clustering algorithm.

Preferably, the K-means clustering algorithm is selected for grouping. And taking the environmental characteristics as input, and applying a K-means clustering algorithm to divide all the cultivated lands into n groups. The specific value of n can be set manually or can be determined through variance analysis, for example, a plurality of values are set, for each value, the environmental characteristics are taken as input, a K-means clustering algorithm is applied to divide a plurality of cultivated lands into a plurality of groups of the values, intra-group variances among cultivated lands in each group are calculated, and the intra-group variances are added to obtain the total variance corresponding to the value. The value with the smallest total variance is determined as the specific value of n.

In one or more optional embodiments of the present invention, for each preset cultivated land input intensity range, according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range, a correlation coefficient and a significance coefficient corresponding to the cultivated land input intensity range are calculated, and the specific implementation process may be as follows:

Specifically, the window width and the sliding step size of the sliding window may be adjusted according to the actual data characteristics.

In practical application, on the basis of determining the input intensity of cultivated land and the output intensity of cultivated land, as shown in fig. 2 and 3, fig. 2 is one of two-dimensional coordinates provided by the present invention, and fig. 3 is the second two-dimensional coordinates provided by the present invention: a two-dimensional graph is established with the cultivated land input intensity as a horizontal axis (x-axis: cultivated land input intensity) and the cultivated land output intensity as a vertical axis (y-axis: cultivated land output intensity), and then each cultivated land is drawn in the two-dimensional graph according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land, and each point in fig. 2 and 3 represents one cultivated land. Then, a sliding window is set, i.e., the window width and the sliding step size of the sliding window are set. Starting from the initial position of the two-dimensional coordinate graph, calculating a correlation coefficient and a significance coefficient according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the current sliding window, wherein the hollow points are each cultivated land in the current sliding window, and the black points are each cultivated land outside the current sliding window; and then the sliding window is moved along the x-axis according to the sliding step length, and the correlation coefficient and the saliency system are calculated according to the cultivated land input intensity and the cultivated land generation intensity of each cultivated land in the current sliding window, and the like until the cultivated lands are not contained in the sliding window. Therefore, the sliding window is used for calculating the promotion or inhibition effect of the farmland input intensity change on the farmland output intensity in different farmland input intensity ranges, so that the dynamic expression of the farmland utilization intensity phase change process and the phase change critical point is realized, and the method has the advantages of clear expression and easiness in popularization and application.

The calculation of the correlation coefficient and the significance coefficient will be described with reference to fig. 2, 3 and 4, and fig. 4 is a schematic flow chart of the calculation of the correlation coefficient and the significance coefficient provided by the invention:

step 401: each cultivated land is arranged in ascending order according to the cultivated land input intensity under the appointed input dimension.

Step 402: setting the window size (window width) of the sliding window as W and the sliding step length as L, respectively extracting each cultivated land in the cultivated land input intensity range [0, W ] (the sliding window is positioned at the initial position, namely at the leftmost side of the x axis), and calculating the correlation coefficient f1 and the significance coefficient p1 of the cultivated land input intensity in the designated input dimension and the cultivated land output intensity in the designated output dimension.

Step 403: and moving the sliding window to a farmland input intensity range [ L, W+L ] and extracting each farmland, and calculating a correlation coefficient f2 and a saliency coefficient p2 of farmland input intensity in a designated input dimension and farmland output intensity in a designated output dimension.

Step 404, moving a sliding window, extracting each cultivated land in the cultivated land input intensity ranges of [2×l, w+2×l ], [3×l, w+3×l ], …, and [ k×l, w+k×l ], and calculating a correlation coefficient and a significance coefficient until w+k×l is greater than or equal to the maximum value of cultivated land input intensity in the specified input dimension.

Illustratively, in fig. 2, the cultivated land input intensity ranges are [32.5, 72.51], the window width w=40, the sliding step l=2.5, the correlation coefficient (partial correlation coefficient) f1= 0.1406, and the saliency coefficient p1<0.01. In fig. 3, the cultivated land input intensity ranges are [37.5, 77.5], the window width w=40, the sliding step length l=2.5, and the significance coefficient p1 > 0.05.

In one or more optional embodiments of the present invention, the determining, according to the correlation coefficient and the significance coefficient corresponding to each of the cultivated land input intensity ranges, a phase transition critical point of cultivated land utilization intensity in the specified input dimension and the specified output dimension includes:

Specifically, the display mode may include at least one of a shape, a color, a line shape, and a size of a dot. The two-dimensional coordinate system may be a cartesian coordinate system.

In practical application, on the basis of determining the correlation coefficient and the significance coefficient, further, taking the average value of the cultivated land input intensity range as an abscissa, taking the correlation coefficient corresponding to the cultivated land input intensity range as an ordinate, and drawing points in a two-dimensional coordinate system according to a specified display mode. Until all the cultivated land input intensity ranges are traversed. After the drawing is completed, determining a designated input dimension and a phase change critical point of the tilling utilization intensity under the designated output dimension according to each point in the two-dimensional coordinate system. Therefore, the correlation coefficient and the significance coefficient of different farmland input intensity ranges are displayed in a coordinate system mode and different display modes, the phase change process can be intuitively and rapidly determined, and further the efficiency of determining the phase change critical point is improved.

Illustratively, each cultivated land is extracted in a cultivated land input intensity range [0, W ] (the sliding window is at the starting position, i.e., at the leftmost side of the x-axis), a correlation coefficient f1 and a saliency coefficient p1 of the cultivated land input intensity in a specified input dimension and the cultivated land output intensity in a specified output dimension are calculated, a value corresponding to the abscissa is an average value W/2 of the cultivated land input intensity range [0, W ], a point (W/2, f 1) is plotted in a cartesian coordinate system with the correlation coefficient f1 as a value corresponding to the ordinate, and the plotted point is displayed according to a display mode corresponding to the saliency coefficient p 1. Moving the sliding window to a cultivated land input intensity range [ L, W+L ] and extracting each cultivated land, calculating a correlation coefficient f2 and a saliency coefficient p2 of cultivated land input intensity in a designated input dimension and cultivated land output intensity in a designated output dimension, using an average value L+W/2 of the cultivated land input intensity range [ L, W+L ] as a numerical value corresponding to an abscissa, using the correlation coefficient f2 as a numerical value corresponding to an ordinate, drawing a point (L+W/2, f 2) in a Cartesian coordinate system, and displaying the drawing point according to a display mode corresponding to the saliency coefficient p 2. And the like, until W+k is equal to or greater than the maximum value of the farmland throwing strength in the appointed throwing dimension, wherein k is any positive integer.

Optionally, the determining the phase change critical point of the tilling strength under the specified input dimension and the specified output dimension according to each point in the two-dimensional coordinate system may be implemented as follows:

taking the first boundary cultivated land input intensity connecting the positive correlation stage and the non-correlation stage as a first phase transition critical point of cultivated land utilization intensity in the appointed input dimension and the appointed output dimension; and taking the second boundary farmland input intensity connecting the uncorrelated stage and the negative correlated stage as a second phase change critical point of the farmland utilization intensity in the designated input dimension and the designated output dimension.

Specifically, the positive correlation stage refers to a stage in which the cultivated land yield strength increases with the increase of the cultivated land input strength; the uncorrelated phase, i.e., the uncorrelated phase, refers to a phase in which the change in the input intensity of cultivated land and the change in the produced intensity are not significantly related. The negative correlation phase is a phase in which the cultivated land input strength increases and the cultivated land yield strength decreases.

In practical application, as the input intensity of the cultivated land increases, the correlation between the input intensity of the cultivated land and the output intensity of the cultivated land undergoes a phase change process of obvious positive correlation, irrelevant, obvious negative correlation, and the positive correlation stage, the irrelevant stage and the negative correlation stage can be identified by analyzing each point in a two-dimensional coordinate system. Then, the farmland input intensity corresponding to the connection point of the positive correlation stage and the non-correlation stage, namely the first boundary farmland input intensity, is used as a first phase change critical point, namely a positive correlation and non-correlation phase change critical point; and taking the farmland input intensity corresponding to the connecting point of the uncorrelated stage and the negative correlated stage, namely the second boundary farmland input intensity, as a second phase change critical point, namely the phase change critical point of uncorrelated and negative correlation. Namely, the value is the lower limit of the numerical range of the sliding window corresponding to the phase change. Thus, the determination efficiency and accuracy of the phase transition critical point can be improved.

In practical application, as the input intensity of the cultivated land increases, the correlation between the input intensity of the cultivated land and the output intensity of the cultivated land undergoes a phase change process of obvious positive correlation, irrelevant, obvious negative correlation, and the positive correlation stage, the irrelevant stage and the negative correlation stage can be identified by analyzing each point in a two-dimensional coordinate system. Then determining the farmland input intensity corresponding to the connection point of the positive correlation stage and the non-correlation stage, namely the first boundary farmland input intensity, identifying the farmland input intensity range corresponding to the first boundary farmland input intensity, and taking the initial value of the farmland input intensity range, namely the first initial value, as a first phase change critical point, namely the positive correlation and non-correlation phase change critical point; determining the farmland input intensity corresponding to the connection point of the uncorrelated stage and the negative correlated stage, namely the second boundary farmland input intensity, identifying the farmland input intensity range corresponding to the second boundary farmland input intensity, taking the starting value of the farmland input intensity range, namely the second starting value, as a first phase change critical point, and taking the starting value of the farmland input intensity range as a second phase change critical point, namely the uncorrelated and negative correlated phase change critical point. Namely, the value is the lower limit of the numerical range of the sliding window corresponding to the phase change. Thus, the determination efficiency and accuracy of the phase transition critical point can be improved.

It should be noted that, because the input intensities of the cultivated lands and the output intensities of the cultivated lands have different dimensions, in order to ensure that the input intensities of the cultivated lands in different input dimensions have the same dimensions, and/or the output intensities of the cultivated lands in different output dimensions have the same dimensions, the method is convenient for data processing, and further improves the determination efficiency of the phase change critical point, and after obtaining the input intensities of the cultivated lands in the specified input dimensions and the output intensities of the cultivated lands in the specified output dimensions, the method further comprises: and normalizing the cultivated land input intensity and/or the cultivated land generation intensity. That is, the input intensity and/or the output intensity of different types of cultivated lands are normalized.

In addition, because of different cultivated land investment dimensions, in order to ensure that the accuracy and the reliability of the critical point are improved, the correlation coefficient is preferably a partial correlation coefficient; correspondingly, for each preset cultivated land input intensity range, calculating a correlation coefficient corresponding to the cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range, wherein the implementation process is as follows: and aiming at each preset cultivated land input intensity range, taking the cultivated land input intensity in a non-designated input dimension as a control variable, and calculating a correlation coefficient corresponding to the cultivated land input intensity range according to the cultivated land input intensity of each cultivated land in the designated input dimension and the cultivated land output intensity in the designated output dimension in the cultivated land input intensity range.

Since there are a plurality of types of input intensities (tilling input intensities in different tilling input dimensions), the output intensity (tilling input intensity) is comprehensively affected. The correlation between the tilling input intensity and the tilling input intensity in each tilling input dimension is represented by a bias correlation coefficient so as to reduce the influence of the tilling input intensity in other tilling input dimensions. The partial correlation coefficient is calculated by the formula (6).

(6)

In formula (6), p _xy|z X is the farmland input intensity in the designated farmland input dimension, y is the farmland input intensity in the other farmland input dimension, z is the farmland input intensity,

residual error, which is a multiple linear regression established between X and Z,>

is the residual of the multiple linear regression established between Y and Z. />

Is->

And->

Is used for the correlation coefficient of the (c) for the (c),

is->

And->

Covariance of->

Is->

Variance of->

Is->

Is a variance of (c).

Due to residual error

The linear correlation between X and Z is eliminated, and (2)>

The linear correlation between Y and Z is eliminated. Thus, calculate +.>

And->

The correlation coefficient between them gives the bias correlation coefficient.

Referring to fig. 5, fig. 5 is one of the schematic diagrams for determining the phase change critical point according to the present invention: the display icons of the cultivated lands with the cultivated land type A are set to be round, the display icons of the cultivated lands with the cultivated land type B are set to be triangular, and the display icons of the cultivated lands with the cultivated land type C are set to be square. For the significance factor, i.e., the p-value belongs to (- ≡0.001), the display mode is: the lines are single solid lines; the p value is [0.001,0.01) is displayed in the following manner: the lines are broken lines; the p-value belonging to [0.01, 0.05) is displayed in the following manner: black is filled in the interior; the value of p is of the order of 0.05, ++ infinity) is displayed in the following manner: the lines are double solid lines. Correspondingly, the p value belongs to (- ≡0.001) and the cultivated land of cultivated land type a is represented by a circle of single solid line, the p value belongs to (-fact0.001) and the cultivated land type B is represented by a triangle of single solid line, and the p value belongs to (-fact0.001) and the cultivated land type C is represented by a square of single solid line; the p value is [0.001,0.01) and the cultivated land of cultivated land type A is represented by a dotted circle, the p value is [0.001,0.01) and the cultivated land of type B is represented by a dotted triangle, and the p value is [0.001,0.01) and the cultivated land of type C is represented by a dotted square; a p value of [0.01, 0.05) and a cultivated land of cultivated land type a is represented by an internally filled black circle, a p value of [0.01, 0.05) a cultivated land of cultivated land type B is represented by an internally filled black triangle, and a p value of [0.01, 0.05) a cultivated land of cultivated land type C is represented by an internally filled black square; the value of p is of the order of 0.05, +++) and the cultivated land of type a is represented by a double solid circle, the value of p is of the order of 0.05, the +++) type of tilling is B tilling and is indicated by the triangle with double solid lines, the value of p is of the order of 0.05, ++ infinity) the cultivated land type C is represented by a square with double solid lines.

Referring next to fig. 6, fig. 6 is a schematic diagram illustrating a second embodiment of determining a phase transition critical point according to the present invention: according to each cultivated land input intensity range and the corresponding partial correlation coefficient and significance coefficient, taking the average value of the cultivated land input intensity range (normalized fertilizer input intensity) as an abscissa, taking the partial correlation coefficient corresponding to the cultivated land input intensity range as an ordinate, and drawing points in a two-dimensional coordinate system according to the display mode and the display icon set in fig. 5. Until all the cultivated land input intensity ranges are traversed.

For cultivated lands of different cultivated land types, calculating the phase transition critical point of the cultivated land utilization intensity according to the influence characteristics of the cultivated land input intensity on the cultivated land output intensity in different cultivated land input intensity ranges, as shown in fig. 5 and 6. Taking the bias correlation coefficient of the chemical fertilizer input intensity to the yield intensity of the main grain crops as an example, the cultivated land type A, the cultivated land type B and the cultivated land type C are respectively expressed by symbols with different shapes, and different display modes are set according to the p value. Each point in the coordinate system corresponds to a cultivated land input intensity range of the chemical fertilizer input intensity, the abscissa of the point is the average value of the cultivated land input intensity ranges, and the ordinate of the point is the offset correlation coefficient of the chemical fertilizer input intensity of the cultivated land and the yield intensity of the main grain crops in the cultivated land input intensity range. For different cultivated land types, as the normalized fertilizer input intensity increases, the correlation between the normalized fertilizer input intensity and the yield intensity of the main grain crops is subjected to a phase change process of obvious positive correlation, irrelevant correlation and obvious negative correlation. Taking cultivated land type A as an example, when the normalized fertilizer input intensity is less than 50, the normalized fertilizer input intensity is obviously and positively correlated with the yield intensity of the main grain crops, the p value is less than 0.01, and the normalized fertilizer input intensity is represented as a circle with a single solid line or a circle with a broken line in a coordinate system. As the normalized fertilizer input strength increases, its impact on the yield strength of the staple food crop is transformed from significant forward drive to insignificant correlation. Based on this, the start value k1 of the fertilizer input intensity range [ k1, k2] corresponding to the circle s1 can be used as a critical point (first phase transition critical point) at which the fertilizer utilization intensity is converted from the positive correlation ordered phase state to the disordered phase state. As the normalized fertilizer input intensity continues to increase, its impact on the yield intensity of the staple food crop is shifted from insignificant correlation to significant negative drive. Based on this, the circle s2 can be used as the critical point for converting the unordered phase state into the negative correlation ordered phase state of the fertilizer utilization intensity corresponding to the start value k3 of the fertilizer input intensity range [ k3, k4 ]. Then, the chemical fertilizer input intensity x1 of each cultivated land in the cultivated land type A is equal to or less than k1, which indicates that the chemical fertilizer input intensity is proper and is in a positive correlation ordered phase state; x1 is greater than k1 and less than or equal to k3, which indicates that the chemical fertilizer input strength is critically unsuitable and is in a disordered phase state; x1> k3 indicates that fertilizer input strength is unsuitable, in a negatively-correlated ordered phase.

In addition, all the steps of the method are realized by using a programming technology, and the suitability of the tilling strength can be determined by adopting the method for determining the phase change critical point of the tilling strength, so that the influence process of the tilling input strength change of different tilling input dimensions on the tilling output strength can be clearly and intuitively displayed, and a basis can be provided for diagnosing the phase change inflection point of the tilling input-output coupling effect. Meanwhile, the method is suitable for calculation of different scales such as point location scale, farm scale, regional scale and national scale.

The device for determining the phase change critical point of the tilling strength provided by the invention is described below, and the device for determining the phase change critical point of the tilling strength and the method for determining the phase change critical point of the tilling strength described above can be referred to correspondingly.

Fig. 7 is a schematic structural diagram of an apparatus for determining a phase transition critical point of a tilling strength according to the present invention, and as shown in fig. 7, the apparatus 700 for determining a phase transition critical point of a tilling strength includes: an acquisition module 701, a calculation module 702 and a determination module 703, wherein:

an acquisition module 701 configured to acquire a cultivated land input intensity of a plurality of cultivated lands in a specified input dimension and a cultivated land output intensity in a specified output dimension, the specified input dimension being any input dimension, the specified output dimension being any output dimension;

A calculation module 702 configured to calculate, for each preset cultivated land input intensity range, a correlation coefficient and a significance coefficient corresponding to the cultivated land input intensity range from the cultivated land input intensity and the cultivated land output intensity of each cultivated land within the cultivated land input intensity range;

a determining module 703 configured to determine a phase transition critical point of the tilling strength in the specified input dimension and the specified output dimension according to the correlation coefficient and the significance coefficient corresponding to each of the tilling input strength ranges.

According to the device for determining the phase change critical point of the tilling and utilizing strength, provided by the invention, the tilling input strength of a plurality of tilling areas in the appointed input dimension and the tilling output strength of the tilling areas in the appointed output dimension are obtained, wherein the appointed input dimension is any input dimension, and the appointed output dimension is any output dimension; calculating a correlation coefficient and a significance coefficient corresponding to each preset cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range; and determining the phase change critical point of the tilling and utilizing intensity in the appointed input dimension and the appointed output dimension according to the correlation coefficient and the significance coefficient corresponding to each tilling and throwing intensity range. By calculating the promotion or inhibition effect of the farmland input intensity change on the farmland output intensity in different farmland input intensity ranges, the dynamic expression of the farmland utilization intensity phase change process and the phase change critical point is realized, the basis can be provided for diagnosing the phase change inflection point of the farmland input-output coupling effect, and the method has the advantages of clear expression and easiness in popularization and application.

Optionally, the apparatus 700 for determining the phase transition critical point of the tilling strength further includes a partitioning module configured to:

accordingly, the computing module 702 is further configured to:

the determining module 703 is further configured to:

Optionally, the computing module 702 is further configured to:

Optionally, the determining module 703 is further configured to:

Optionally, the apparatus 700 for determining a phase transition critical point of a tilling strength further includes a normalization module configured to:

Optionally, the correlation coefficient is a partial correlation coefficient;

accordingly, the computing module 702 is further configured to:

Fig. 8 illustrates a physical structure diagram of an electronic device, as shown in fig. 8, which may include: processor 810, communication interface (Communications Interface) 820, memory 830, and communication bus 840, wherein processor 810, communication interface 820, memory 830 accomplish communication with each other through communication bus 840. The processor 810 may invoke logic instructions in the memory 830 to perform a method of determining a phase change critical point for a tilling strength, the method comprising: acquiring farmland input intensity of a plurality of farmland in a designated input dimension and farmland output intensity in a designated output dimension, wherein the designated input dimension is any input dimension, and the designated output dimension is any output dimension; calculating a correlation coefficient and a significance coefficient corresponding to each preset cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range; and determining the phase change critical point of the tilling and utilizing intensity in the appointed input dimension and the appointed output dimension according to the correlation coefficient and the significance coefficient corresponding to each tilling and throwing intensity range.

Further, the logic instructions in the memory 830 described above may be implemented in the form of software functional units and may be stored in a computer-readable storage medium when sold or used as a stand-alone product. Based on this understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server, a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.

In another aspect, the present invention also provides a computer program product comprising a computer program, the computer program being storable on a non-transitory computer readable storage medium, the computer program, when executed by a processor, being capable of performing the method of determining a phase change critical point of a tilling use intensity provided by the methods described above, the method comprising: acquiring farmland input intensity of a plurality of farmland in a designated input dimension and farmland output intensity in a designated output dimension, wherein the designated input dimension is any input dimension, and the designated output dimension is any output dimension; calculating a correlation coefficient and a significance coefficient corresponding to each preset cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range; and determining the phase change critical point of the tilling and utilizing intensity in the appointed input dimension and the appointed output dimension according to the correlation coefficient and the significance coefficient corresponding to each tilling and throwing intensity range.

In yet another aspect, the present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, is implemented to perform the method of determining a phase change critical point of a tilling strength provided by the methods described above, the method comprising: acquiring farmland input intensity of a plurality of farmland in a designated input dimension and farmland output intensity in a designated output dimension, wherein the designated input dimension is any input dimension, and the designated output dimension is any output dimension; calculating a correlation coefficient and a significance coefficient corresponding to each preset cultivated land input intensity range according to the cultivated land input intensity and the cultivated land output intensity of each cultivated land in the cultivated land input intensity range; and determining the phase change critical point of the tilling and utilizing intensity in the appointed input dimension and the appointed output dimension according to the correlation coefficient and the significance coefficient corresponding to each tilling and throwing intensity range.

The apparatus embodiments described above are merely illustrative, wherein the elements illustrated as separate elements may or may not be physically separate, and the elements shown as elements may or may not be physical elements, may be located in one place, or may be distributed over a plurality of network elements. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment. Those of ordinary skill in the art will understand and implement the present invention without undue burden.

From the above description of the embodiments, it will be apparent to those skilled in the art that the embodiments may be implemented by means of software plus necessary general hardware platforms, or of course may be implemented by means of hardware. Based on this understanding, the foregoing technical solution may be embodied essentially or in a part contributing to the prior art in the form of a software product, which may be stored in a computer readable storage medium, such as ROM/RAM, a magnetic disk, an optical disk, etc., including several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the method described in the respective embodiments or some parts of the embodiments.

Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and are not limiting; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method of determining a phase transition critical point of a tilling strength, comprising:

2. The method for determining a phase transition critical point of a tilling strength according to claim 1, wherein before calculating a correlation coefficient and a significance coefficient corresponding to each tilling input strength range from the tilling input strength and the tilling output strength of each tilling within the tilling input strength range for each preset tilling input strength range, further comprising:

3. The method for determining a phase transition critical point of a tilling use intensity according to claim 1, wherein the calculating a correlation coefficient and a significance coefficient corresponding to each tilling input intensity range according to the tilling input intensity and the tilling output intensity of each tilling within the tilling input intensity range for each preset tilling input intensity range includes:

4. The method of determining a phase transition critical point of a tillable strength according to any one of claims 1-3, wherein the determining a phase transition critical point of a tillable strength in the specified input dimension and the specified output dimension from the correlation coefficient and the saliency coefficient corresponding to each of the tillable input strength ranges comprises:

5. The method of determining a phase transition critical point of a tillable strength according to claim 4, wherein determining a phase transition critical point of a tillable strength in the specified input dimension and the specified output dimension from each point in the two-dimensional coordinate system comprises:

6. The method of determining a phase transition critical point of a tilling strength according to claim 1, wherein after the acquiring of the tilling input strength of the plurality of tilling areas in the designated input dimension and the tilling output strength in the designated output dimension, further comprising:

7. The method of determining a phase transition critical point of a tilling strength of claim 1 wherein the correlation coefficient is a partial correlation coefficient;

8. An apparatus for determining a phase transition critical point of a tilling strength, comprising:

9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of determining the phase transition critical point of a tilling strength as claimed in any one of claims 1 to 7 when the program is executed.

10. A non-transitory computer readable storage medium having stored thereon a computer program, wherein the computer program when executed by a processor implements a method of determining a phase transition critical point of a tilling strength as claimed in any one of claims 1 to 7.