CN107292755A - A kind of Analysis on Selecting method and device in corn planting environment Typical Representative area - Google Patents
A kind of Analysis on Selecting method and device in corn planting environment Typical Representative area Download PDFInfo
- Publication number
- CN107292755A CN107292755A CN201610221767.6A CN201610221767A CN107292755A CN 107292755 A CN107292755 A CN 107292755A CN 201610221767 A CN201610221767 A CN 201610221767A CN 107292755 A CN107292755 A CN 107292755A
- Authority
- CN
- China
- Prior art keywords
- zoning
- area
- minimum
- typical
- cluster
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 240000008042 Zea mays Species 0.000 title claims abstract description 38
- 235000002017 Zea mays subsp mays Nutrition 0.000 title claims abstract description 38
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 title claims abstract description 37
- 235000005822 corn Nutrition 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000004458 analytical method Methods 0.000 title claims abstract description 18
- 238000013316 zoning Methods 0.000 claims abstract description 157
- 230000010354 integration Effects 0.000 claims abstract description 17
- 238000009395 breeding Methods 0.000 claims abstract description 10
- 230000001488 breeding effect Effects 0.000 claims abstract description 10
- 230000007613 environmental effect Effects 0.000 claims description 23
- 230000001413 cellular effect Effects 0.000 claims description 18
- 230000035558 fertility Effects 0.000 claims description 9
- 238000005192 partition Methods 0.000 abstract description 5
- 238000009825 accumulation Methods 0.000 description 11
- 230000007261 regionalization Effects 0.000 description 7
- 238000010899 nucleation Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 5
- 238000009331 sowing Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000010606 normalization Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000003139 buffering effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 230000032696 parturition Effects 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 238000000638 solvent extraction Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 238000012952 Resampling Methods 0.000 description 1
- 235000016383 Zea mays subsp huehuetenangensis Nutrition 0.000 description 1
- 238000012271 agricultural production Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 235000009973 maize Nutrition 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of Analysis on Selecting method and device in corn planting environment Typical Representative area, method includes：According to minimum zoning dividing elements and area data, aggregate-value and year average in each division index annual breeding time in each minimum zoning unit are calculated；Space attribute integration is carried out to minimum zoning unit according to the year average of each division index to cluster, and obtains the integrated environment zoning of crop-planting environment；According to the weight of each division index and year average, it is determined that the feature description of each integrated environment zoning；All minimum zoning units of all years are clustered, the fluctuation situation of zoning after each cluster is calculated；Zoning after cluster is analyzed according to feature description and fluctuation situation, the Typical Representative area of corn planting environment is chosen.The present invention is clustered by the year average of division index to zoning unit, and partition boundaries are clearly fine；And described by the feature of different integrated environment zonings and fluctuation situation, choose corn planting environment Typical Representative area.
Description
Technical field
The present invention relates to IT application to agriculture technical field, and in particular to a kind of corn planting environment allusion quotation
The Analysis on Selecting method and device of type Representative Region.
Background technology
Agricultural planting zoning will be ground according to principle of similarity according to agricultural production conditions and feature
Study carefully region and carry out subregion scribing so that weather conditions, cropping system, soil in same subregion
The characteristics such as condition have similar feature, and have larger otherness between different subregions.I
Therefore the condition differences such as heat, moisture, illumination, the soil of state various regions significantly, also determine me
The otherness of the inequality and variety selection, planting system of state's crop varieties plantation distribution etc., be
Optimize crop varieties planting structure distribution, planting scale is realized, according to different ecological condition Yin Shiyin
Selection kind, the mode of production on ground, according to local environment situation selection, setting breed breeding and
The demand of the test station of popularization etc. so that carry out crop varieties planting regionalization tool from different perspectives
It is significant.
Research for the zoning methods of proportion of crop planting environment is more, and main zoning methods are big
Cause is divided to two classes, the first kind be according to expertise, using corn planting environment fitness as foundation,
Zoning environment attribute value is divided into some suitable levels, and the suitability progress of multiple attributes is comprehensive
Overlap-add procedure is closed, suitability commentary is carried out to obtained subregion, the attribute suitability degree of the party is interval
Value is difficult to determine, gained Zoning subjectivity is strong, and zoning border is difficult to determine.Equations of The Second Kind is
The attribute of selected planting regionalization, is clustered, the method such as disaggregated model carries out region with space attribute
Divide, this method determines territorial classification in the way of quantification, relative to first way result
It is more objective but explanatory not strong.Therefore by data and method are limited, existing crop-planting
Often obscure boundary, subregion yardstick are larger for zoning, it is difficult to express in each biome different subprovinces,
And the difference between niche.Therefore, existing crop-planting environmental regionalization method can not
The attributive character of the different zonings of crops is recognized simultaneously, clearly fine partition boundaries, so as to carry out
Fluction analysis.
The content of the invention
Because existing crop-planting environmental regionalization method can not recognize crops not same district simultaneously
The attributive character drawn, clearly fine partition boundaries are of the invention so as to the problem of carrying out fluction analysis
Propose a kind of Analysis on Selecting method and device in corn planting environment Typical Representative area.
In a first aspect, the present invention proposes a kind of Analysis on Selecting in corn planting environment Typical Representative area
Method, including：
According to minimum zoning dividing elements and area data, calculate every in each minimum zoning unit
Aggregate-value in individual division index annual breeding time, each division index is calculated according to aggregate-value
Year average；
Space attribute one is carried out to minimum zoning unit according to the year average of each division index
Change cluster, obtain the integrated environment zoning of crop-planting environment；
According to each division index in the weight of each division index and each minimum zoning unit
Year average, it is determined that the feature description of each integrated environment zoning；
Space attribute integration cluster was carried out to all minimum zoning units of all years, calculates every
The fluctuation situation of zoning after individual cluster；
Zoning after cluster is analyzed according to feature description and fluctuation situation, corn planting is chosen
The Typical Representative area of environment.
Preferably, the Typical Representative area includes：Typical average Representative Volume Element area, Dian Xingjun
It is worth special cellular zone, Typical stability Representative Volume Element area, the special cellular zone of Typical stability.
Preferably, it is described according to minimum zoning dividing elements and area data, calculate each minimum
Aggregate-value in zoning unit in each division index annual breeding time, calculates every according to aggregate-value
The year of individual division index is before average, in addition to：
Obtain the area data in region to be analyzed.
Preferably, described pair of all minimum zoning units of all years carry out space attribute integration
Cluster, is calculated after each cluster after the fluctuation situation of zoning, in addition to：
According to the attributive character of each integrated environment zoning, the class of each integrated environment zoning is calculated
Not Bo Dong the frequency and ownership probability, and the frequency and the ownership probability are fluctuated according to the classification,
The fluctuation situation of zoning after being clustered.
Preferably, the area data includes：Geo-spatial data, crop are in region to be analyzed
For many years fertility issue according to this and region to be analyzed environmental data for many years.
Second aspect, the present invention also proposes a kind of selection in corn planting environment Typical Representative area point
Analysis apparatus, including：
Year average computing module, for according to minimum zoning dividing elements and area data, calculating
Aggregate-value in each minimum zoning unit in each division index annual breeding time, according to accumulative
Value calculates the year average of each division index；
Integrated environment zoning acquisition module, for according to each division index year average to minimum
Zoning unit carries out space attribute integration cluster, obtains the integrated environment area of crop-planting environment
Draw；
Feature describes determining module, for the weight according to each division index and each smallest region
The year average of each division index in unit is drawn, it is determined that the feature of each integrated environment zoning is retouched
State；
Fluctuation situation computing module, space is carried out for all minimum zoning units to all years
Attribute integration cluster, calculates the fluctuation situation of zoning after each cluster；
Typical Representative area chooses module, for situation to be described and fluctuated according to feature to cluster back zone
Draw and analyzed, choose the Typical Representative area of corn planting environment.
Preferably, the Typical Representative area includes：Typical average Representative Volume Element area, Dian Xingjun
It is worth special cellular zone, Typical stability Representative Volume Element area, the special cellular zone of Typical stability.
Preferably, in addition to：
Region query module, the area data for obtaining region to be analyzed.
Preferably, in addition to：
Computing module is fluctuated, for the attributive character according to each integrated environment zoning, calculates every
The classification fluctuation frequency and ownership probability of individual integrated environment zoning, and frequency is fluctuated according to the classification
Secondary and described ownership probability, the fluctuation situation of zoning after being clustered.
Preferably, the area data includes：Geo-spatial data, crop are in region to be analyzed
For many years fertility issue according to this and region to be analyzed environmental data for many years.
As shown from the above technical solution, the present invention by division index year average to zoning unit
Clustered, partition boundaries are clearly fine；And described by the feature of different integrated environment zonings
With fluctuation situation, corn planting environment Typical Representative area is chosen.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below
The accompanying drawing used required in embodiment or description of the prior art will be briefly described, show and
Easy insight, drawings in the following description are only some embodiments of the present invention, for this area
, on the premise of not paying creative work, can also be according to these for those of ordinary skill
Figure obtains other accompanying drawings.
A kind of corn planting environment Typical Representative area that Fig. 1 provides for one embodiment of the invention
The schematic flow sheet of Analysis on Selecting method；
A kind of Three Eastern Provinces corn planting environment comprehensive ring that Fig. 2 provides for one embodiment of the invention
Border Zoning and specifically distinguish Butut；
A kind of Three Eastern Provinces corn planting environment year border ripple that Fig. 3 provides for one embodiment of the invention
Dynamic zoning belongs to probability analysis figure with pixel；
A kind of 20 year year of Three Eastern Provinces corn planting environment that Fig. 4 provides for one embodiment of the invention
Fluctuate the frequency in border；
A kind of Three Eastern Provinces corn planting environment comprehensive ring that Fig. 5 provides for one embodiment of the invention
Border zoning border and year border fluctuation zoning stack result；
A kind of corn planting environment Typical Representative area that Fig. 6 provides for one embodiment of the invention
The structural representation of Analysis on Selecting device.
Embodiment
Below in conjunction with the accompanying drawings, the embodiment to invention is further described.Implement below
Example is only used for clearly illustrating technical scheme, and can not limit this hair with this
Bright protection domain.
Fig. 1 shows a kind of corn planting environment Typical Representative that one embodiment of the invention is provided
The Analysis on Selecting method in area, including：
S1, according to minimum zoning dividing elements and area data, calculate each minimum zoning unit
In aggregate-value in each division index annual breeding time, each zoning is calculated according to aggregate-value and referred to
Target year average；
S2, according to each division index year average to minimum zoning unit carry out space attribute one
Bodyization is clustered, and obtains the integrated environment zoning of crop-planting environment；
Each zoning refers in S3, the weight according to each division index and each minimum zoning unit
Target year average, it is determined that the feature description of each integrated environment zoning；
S4, all minimum zoning units of all years are carried out with space attributes integration cluster, meter
Calculate the fluctuation situation of zoning after each cluster；
S5, according to feature description and fluctuation situation zoning after cluster is analyzed, selection corn
The Typical Representative area of planting environment.
The present embodiment is clustered by the year average of division index to zoning unit, partition boundaries
It is clearly fine；And described by the feature of different integrated environment zonings and fluctuation situation, choose beautiful
Rice planting environment Typical Representative area.
Alternatively, the Typical Representative area includes：Typical average Representative Volume Element area, Dian Xingjun
It is worth special cellular zone, Typical stability Representative Volume Element area, the special cellular zone of Typical stability.
The Typical Representative area different by choosing, can be promoted and testing station reconnaissance for crop varieties
Decision support is provided.
Further, before S1, in addition to：
S0, the area data for obtaining region to be analyzed.
Further, in order to which the fluctuation situation for obtaining crop-planting environment is distributed and environmental characteristic
Further to be analyzed, after S4, in addition to：
S45, the attributive character according to each integrated environment zoning, calculate each integrated environment area
The classification fluctuation frequency and ownership probability drawn, and the frequency and the ownership are fluctuated according to the classification
Probability, the fluctuation situation of zoning after being clustered.
Specifically, the area data includes：Geo-spatial data, crop are in region to be analyzed
For many years fertility issue according to this and region to be analyzed environmental data for many years.
, can by collecting geo-spatial data, for many years fertility issue evidence and for many years environmental data
The characteristics of reflecting crop-planting environment more fully hereinafter.
In order to which the corn planting environment Typical Representative area that the present embodiment is provided is described in more detail
Specific steps are illustrated by Analysis on Selecting method below：
Step 1：The area data in region to be analyzed is obtained, the area data includes：Basis
The issue of fertility for many years of geodata, crop in the region to be analyzed is according to this and described to be analyzed
The environmental data for many years in region；
Step 2：It is determined that minimum zoning unit and zoning attribute, calculate each zoning minimum unit
In annual sowing time to maturity period planting environment zoning attribute accumulation value and for many years sowing time to giving birth to
Educate the average of attribute accumulation value in the phase.
Step 3:Using attribute average for many years as cluster attribute, all zoning units in research area are carried out
Space attribute integration cluster, obtains crop-planting environment comprehensive environmental regionalization；Identification is special single
Member；In each subregion, it is determined that the attribute larger to division result otherness contribution rate, and respectively
Attribute interval value, the final agronomy description information for obtaining Zoning；
To in the clustering of zoning minimum unit, first against many attributes of zoning minimum unit
Hierarchical cluster attribute is carried out, wherein cluster attribute is the attribute average of minimum zoning unit, according to each category
Property to crop-planting zoning influence size setting different attribute different weights and carry out normalizing
Change processing and eliminate dimension, be secondly the spatial continuity for ensureing Zoning, can be by minimum zoning
The X of unit, Y-coordinate and set weight as the attribute of hierarchical cluster attribute.Simultaneously should be according to sky
Between Continuous Adjustment criterion carry out the adjustment of unit in small, broken bits, obtain space attribute integration cluster knot
Really.
Calculate each zoning cell attribute value and its affiliated all zoning list of subregion in all subregions
The distance between attribute average of member, because these distance values meet Normal Distribution Characteristics, in setting
Significance a after, the zoning unit of region of rejection is turned into special area.
To each subregion, variance yields is calculated after all each attribute weight of unit normalization (wherein
Weight is consistent when setting with cluster), obtain minimum and maximum to each subregion otherness contribution rate
Attribute.
The statistics such as extreme value, variance size before each attribute weight normalization of each subregion are calculated to obtain
The agronomy description information of scoring area result.
Step 4：, will be all using the annual property value of all minimum zoning units as cluster attribute
All zoning minimum units in time carry out space attribute integration cluster, calculate each zoning list
Member for many years between the classification fluctuation frequency and ownership probability, obtain crop-planting environmental fluctuating situation point
Cloth, and recognize strong, the weak variation zone of environmental characteristic for many years；
According to the same Attributions selection of step 3 and weight setting method, by owning for all times
The mixing of zoning minimum unit is carried out in clustering, obtained cluster result, on geographical space
Same unit may belong to different classifications between not the same year.
Each minimum zoning unit probability for belonging to some class interior for many years on geographical space is calculated,
According to the environment of the Probit Analysis unit year border fluctuation situation, using the classification of maximum probability as
The final belonging kinds of the zoning unit, classification belongs to probability to the unit in generation research area for many years
Figure, year border most probable value be less than threshold value the fluctuation of unit year border it is strong, be defined as typical single
Member.
The distance between any each two class value is calculated, setup unit fluctuation frequency initial value is 0,
By each minimum zoning unit on geographical space, in annual belonging kinds value, sequence is arranged per year
Row, when unit generic changes between year border, fluctuation frequency value plus change class it
Between distance value, the fluctuation frequency accumulated value between each unit completes all years borders, generation
The unit year border classification fluctuation frequency diagram in area is studied, the year of unit synthesized attribute between border is used as
One metric of fluctuation situation.
Step 5：Integrated environment average is realized to planting environment with the result of fluctuation situation zoning
Minimum zoning unit environments are representational to be evaluated, and typical environment unit is selected, including typical case
Average Representative Volume Element area, the typical special cellular zone of average, Typical stability Representative Volume Element area,
The special cellular zone of Typical stability, promotes for crop varieties and provides decision-making branch with testing station reconnaissance
Hold.
The class center of each subregion is used as the subregion kind in computing environment integrated environment Zoning
Environment long-run average representative degree highest mesh region is planted, the representative grid of typical average is used as；
The special area of integrated environment zoning that is obtained by step 3 is selected as the typical case of each subregion
The special grid of average.
Be superimposed environment integrated environment zoning zoning border and unit for many years classification ownership probability graph,
Unit year border classification fluctuation frequency diagram, all rings to finding presence in each Partitioning boundary
Situation divisional type is fluctuated in border, and selection class makes a variation small grid as Typical stability and represents personality
Net, the larger grid of simultaneous selection class variation is used as the special grid of Typical stability.
Mean apparent and Nian Jibo of the method that the present embodiment is provided from crop-planting environmental aspect
The dynamic aspect of performance two carries out planting environment zoning, the identification of region subenvironment, stress risk assessment,
Using hierarchical cluster attribute, zoning border is specified, the agronomy description of division result is realized, for instructing
Kind is tested and popularization is significant, while being conducive to improving to proportion of crop planting ring
The cognition of border spatial distribution, determine region main environment feature and recognize special subenvironment region with
Feature, promotes kind to be accurately positioned advantage and promotes area, improve kind testing efficiency.
For example, Liaoning, Jilin, Heilungkiang San Sheng are the main of eastern North China Spring maize seeding area
Planting area, the present invention is using the Three Eastern Provinces as region to be analyzed, using corn as zoning crop, including
Following steps：
Step 1：The area data in region to be analyzed is obtained, the area data includes：Basis
The issue of fertility for many years of geodata, crop in the region to be analyzed is according to this and described to be analyzed
The environmental data for many years in region；
Obtaining the region base geodata by region to be analyzed includes administrative stroke provincial, at county level
Divided data, DEM altitude datas.The agriculture of 20 years in acquisition Three Eastern Provinces area buffering area 100km
Industry meteorological site data, the agricultural weather station data mainly includes the annual corn seeding date
With the ripe date.The meteorological site data of 20 years in acquisition Three Eastern Provinces area buffering area 100km,
The meteorological site data include mean daily temperature, max. daily temperature, Daily minimum temperature, day drop
Rainfall, day sunshine time etc., wherein the meteorological site and agricultural weather website are in the region
It is evenly distributed, and rejects height above sea level and the excessive meteorological site of zone leveling height above sea level difference.
Calculate respectively agricultural weather website it is annual from January 1 to the corn seeding phase with to corn into
The number of days of ripe phase, using radial base interpolation method obtain full survey region 20 years from January 1
Day is converted to specific date date to corn seeding and ripe number of days, by two dates
Assignment is in meteorological site, obtaining the annual sowing of each meteorological site and ripe date.
Step 2：It is determined that minimum zoning unit and zoning attribute, calculate each zoning minimum unit
In annual sowing time to maturity period planting environment zoning attribute accumulation value and for many years sowing time to giving birth to
Educate the average of attribute accumulation value in the phase.
The present embodiment is divided into 10kmx10km grid by zoning is studied, using grid as minimum zoning
Unit, to study the accumulation accumulated temperature in area's corn seeding phase to maturity period, accumulation rainfall, accumulation day
According to when number, DEM height values be zoning attribute.
Calculate the accumulation mean daily temperature in meteorological site annual corn seeding phase to maturity period, tire out
Product rainfall, accumulation sunshine time, utilize each grid of Kriging regression method acquisition tiring out every year
Product accumulated temperature, accumulation rainfall, accumulation sunshine value, each grid is obtained using method for resampling
DEM height values.
In addition, calculating the average value of each attribute in each grid 20 years.
Step 3:Many attribute space attributes are carried out to the equal value attribute for studying all zoning units in area
Integration cluster, obtains crop-planting environment average zoning；Recognize special unit；At each point
In area, similitude and the larger attribute of otherness, and each attribute interval value are determined, it is final to obtain
The agronomy description information of Zoning；
In the present embodiment, to 20 years of all grids of survey region each attribute mean normalization
To eliminating dimension between 0-100 and assigning different weights, hierarchical cluster attribute, this case study are carried out
Middle that accumulative accumulated temperature, accumulative rainfall, accumulative sunlight weight are set into 0.25, elevation weight is set to
0.15, X, Y-coordinate weight are entered as 0.05 respectively where grid element center, using k-means
The close and distant degree calculating of clustering, wherein sample uses euclidean distance method, according to the R partially of R side, half
The size of square statistic determines final clusters number (preferable clustering number mesh is 8 in the present embodiment),
Wherein X, the addition of Y-coordinate ensure that the preferable spatial continuity of Zoning.
Special area choosing method is as follows：Calculate in all subregions each grid property value with belonging to it
Euclidean distance between each attribute average of all grids of subregion, to different subregions, special area
Selection should be based on identical standard, because of grid and class center where it of all different subregions
Distance value meet Normal Distribution Characteristics, after the significance a of setting, by region of rejection
Zoning unit turns to special area, and special area is the outlier for clustering zoning, is in practical application
It should be specifically noted that, it is as shown in Figure 2 that zoning chooses result with special area.
Variance yields after the normalization of each attribute weight of all grids is calculated in each subregion (wherein to weigh
Reset consistent when putting with cluster), draw minimum and maximum to each subregion otherness contribution rate
Attribute, it is as shown in table 1 below.
The maximum category of Three Eastern Provinces corn planting environment comprehensive environmental regionalization result difference contribution of table 1
Property
| VAR (accumulated temperature) | VAR (rainfall) | VAR (elevation) | VAR (sunshine) | |
| 1 | 18.00398 | 16.54014 | 6.591056 | 15.70887 |
| 2 | 8.154028 | 34.62356 | 12.80553 | 8.623617 |
| 3 | 17.2428 | 12.71221 | 2.072623 | 12.53993 |
| 4 | 6.777262 | 9.356085 | 9.795064 | 5.142063 |
| 5 | 18.91218 | 18.15093 | 4.197132 | 15.55342 |
| 6 | 24.29578 | 7.252329 | 3.867112 | 17.12033 |
| 7 | 10.52315 | 4.678987 | 12.83869 | 4.939667 |
| 8 | 10.56928 | 18.94914 | 13.40226 | 11.96777 |
The true standard difference and extreme value of all each attributes of grid in each subregion are calculated, this point is obtained
The attributive character in area, obtains the agronomy description of the subregion.
Step 4：, will be all using the annual property value of all minimum zoning units as cluster attribute
All zoning minimum units in time carry out space attribute integration cluster, calculate each zoning list
Member for many years between the classification fluctuation frequency and ownership probability, obtain crop-planting environmental fluctuating situation point
Cloth, and recognize strong, the weak variation zone of environmental characteristic for many years；
The grid divided on survey region geographical space is referred to as " geographic grid ", this research is real
7832 geographic grids are had in example, each geographic grid arranges annual right in temporal sequence
A grid for containing multiple attribute value attributes is answered, is referred to as " value grid ", in this case study
In have the data of 20 years, therefore each geographic grid has 20 value grids, all geography networks
The value grid of lattice has 7832*20.All these value grids are put together cluster, wherein
Attribute and attribute weight setting, the determination mode of clusters number, clustering method etc. with step 3
Just the same (preferable clustering number mesh is 7 in this example).Each geographic grid pair is known more than
The multiple value grids answered must have identical dem height values and X, Y scale value.What is drawn is poly-
The value grid in not the same year of same geographic grid may belong to different classifications in class result.
The probability that each geographic grid belongs to some class is calculated, circular is to calculate ground
Reason grid all values grid belongs to the frequency of some class, with the frequency divided by year, that is, obtains ground
Reason grid belongs to such probability, if unnecessary one of the classification of value grid, the geographic grid
The class belonged to there is a variety of possibility, and this case study is regard the maximum class of rate of imputation as this
The final ownership class of geographic grid, classification belongs to probability graph to generation research area's grid for many years, in agriculture
Learn and said in meaning, if to belong to some class probability very big for some geographic grid, illustrated comparatively
The environmental modification very little of the geographic grid, if the probability of the maximum class of some geographic grid degree of membership
Value is less than threshold value (being set to 30% here), then is set to typical wave grid.It is as follows
In legend in shown in Fig. 3, black grid is typical wave grid, and remaining mesh color contains
Justice is：Hundreds are the class number that geographic grid finally belongs to, ten and a bit identification ownership
In such percent probability, such as 295, which represent the geographic grid, belongs to the probability of the 2nd class and is
95%, it is 100% that 700, which represent the geographic grid to belong to the probability of the 7th class,.
Each geographic grid environment category fluctuation frequency is calculated, specific method is：It is every more than obtaining
The cluster centre of individual class, calculates the Euclidean distance of all cluster centres, is used as the distance between class
Metric.For the setting of each geographic grid " the value grid fluctuation frequency " attribute, initial value is set
For 0, by all values grid of each geographic grid, sequence is arranged per year, the value net in relatively more adjacent year
Lattice generic, when two class label differences, the frequency is fluctuated by the value grid of each geographic grid
Plus the distance between change class value, until each geographical grid completes the value net between all years borders
The classification of lattice compares to be superimposed with the value grid fluctuation frequency.The border classification fluctuation of generation research area's grid year
Frequency diagram, as shown in figure 4, regarding value grid fluctuation frequency value as geographical grid year border fluctuation feelings
The measurement attribute of condition.
Step 5：Integrated environment average is realized to planting environment with the result of fluctuation situation zoning
Minimum zoning unit environments are representational to be evaluated, and typical environment unit is selected, including typical case
Average Representative Volume Element area, the typical special cellular zone of average, Typical stability Representative Volume Element area,
The special cellular zone of Typical stability, promotes for crop varieties and provides decision-making branch with testing station reconnaissance
Hold.
The class center of each subregion is used as the subregion kind in computing environment integrated environment Zoning
Environment long-run average representative degree highest mesh region is planted, as the representative grid of typical average,
The special area of integrated environment zoning obtained by step 3 is selected as the special grid of typical average.It is folded
Plus zoning border and the environmental fluctuating situation Zoning (such as Fig. 5) of environment comprehensive environmental regionalization,
All environmental fluctuating situation divisional types to finding presence in each Partitioning boundary, selection
Class makes a variation small grid as the representative grid of Typical stability, and the variation of simultaneous selection class is larger
Grid is used as the special grid of Typical stability.
Fig. 6 shows a kind of choosing in corn planting environment Typical Representative area that the present embodiment is provided
Analytical equipment is taken, including：
Year average computing module 11, for according to minimum zoning dividing elements and area data,
The aggregate-value in each division index annual breeding time in each minimum zoning unit is calculated, according to
Aggregate-value calculates the year average of each division index；
Integrated environment zoning acquisition module 12, for the year average pair according to each division index
Minimum zoning unit carries out space attribute integration cluster, obtains the synthesis ring of crop-planting environment
Border zoning；
Feature describes determining module 13, for the weight according to each division index and it is each most
The year average of each division index in small zoning unit, it is determined that the feature of each integrated environment zoning
Description；
Fluctuation situation computing module 14, is carried out for all minimum zoning units to all years
Space attribute integration cluster, calculates the fluctuation situation of zoning after each cluster；
Typical Representative area chooses module 15, for situation to be described and fluctuated according to feature to cluster
Zoning is analyzed afterwards, chooses the Typical Representative area of corn planting environment.
Alternatively, the Typical Representative area includes：Typical average Representative Volume Element area, Dian Xingjun
It is worth special cellular zone, Typical stability Representative Volume Element area, the special cellular zone of Typical stability.
Further, in addition to：
Region query module, the area data for obtaining region to be analyzed.
Further, in addition to：
Computing module is fluctuated, for the attributive character according to each integrated environment zoning, calculates every
The classification fluctuation frequency and ownership probability of individual integrated environment zoning, and frequency is fluctuated according to the classification
Secondary and described ownership probability, the fluctuation situation of zoning after being clustered.
Further, the area data includes：Geo-spatial data, crop are to be analyzed
Region for many years fertility issue according to this and region to be analyzed environmental data for many years.
A kind of Analysis on Selecting device in corn planting environment Typical Representative area described in the present embodiment
It can be used for performing above method embodiment, its principle is similar with technique effect, no longer go to live in the household of one's in-laws on getting married herein
State.
In the specification of the present invention, numerous specific details are set forth.It is to be appreciated, however, that this
The embodiment of invention can be put into practice in the case of these no details.In some examples
In, known method, structure and technology is not been shown in detail, so as not to fuzzy to this specification
Understanding.
Claims (10)
1. a kind of Analysis on Selecting method in corn planting environment Typical Representative area, it is characterised in that
Including：
According to minimum zoning dividing elements and area data, calculate every in each minimum zoning unit
Aggregate-value in individual division index annual breeding time, each division index is calculated according to aggregate-value
Year average；
Space attribute one is carried out to minimum zoning unit according to the year average of each division index
Change cluster, obtain the integrated environment zoning of crop-planting environment；
According to each division index in the weight of each division index and each minimum zoning unit
Year average, it is determined that the feature description of each integrated environment zoning；
Space attribute integration cluster was carried out to all minimum zoning units of all years, calculates every
The fluctuation situation of zoning after individual cluster；
Zoning after cluster is analyzed according to feature description and fluctuation situation, corn planting is chosen
The Typical Representative area of environment.
2. according to the method described in claim 1, it is characterised in that the Typical Representative area
Including：Typical average Representative Volume Element area, the special cellular zone of typical average, Typical stability generation
The special cellular zone of table cellular zone, Typical stability.
3. method according to claim 2, it is characterised in that described according to smallest region
Dividing elements and area data are drawn, each division index in each minimum zoning unit is calculated annual
Aggregate-value in breeding time, calculates year of each division index before average, also according to aggregate-value
Including：
Obtain the area data in region to be analyzed.
4. method according to claim 3, it is characterised in that described pair all years
All minimum zoning units carry out space attribute integration cluster, calculate zoning after each cluster
After fluctuation situation, in addition to：
According to the attributive character of each integrated environment zoning, the class of each integrated environment zoning is calculated
Not Bo Dong the frequency and ownership probability, and the frequency and the ownership probability are fluctuated according to the classification,
The fluctuation situation of zoning after being clustered.
5. method according to claim 4, it is characterised in that the area data bag
Include：The issue of fertility for many years of geo-spatial data, crop in region to be analyzed is according to this and to be analyzed
The environmental data for many years in region.
6. a kind of Analysis on Selecting device in corn planting environment Typical Representative area, it is characterised in that
Including：
Year average computing module, for according to minimum zoning dividing elements and area data, calculating
Aggregate-value in each minimum zoning unit in each division index annual breeding time, according to accumulative
Value calculates the year average of each division index；
Integrated environment zoning acquisition module, for according to each division index year average to minimum
Zoning unit carries out space attribute integration cluster, obtains the integrated environment area of crop-planting environment
Draw；
Feature describes determining module, for the weight according to each division index and each smallest region
The year average of each division index in unit is drawn, it is determined that the feature of each integrated environment zoning is retouched
State；
Fluctuation situation computing module, space is carried out for all minimum zoning units to all years
Attribute integration cluster, calculates the fluctuation situation of zoning after each cluster；
Typical Representative area chooses module, for situation to be described and fluctuated according to feature to cluster back zone
Draw and analyzed, choose the Typical Representative area of corn planting environment.
7. device according to claim 6, it is characterised in that the Typical Representative area
Including：Typical average Representative Volume Element area, the special cellular zone of typical average, Typical stability generation
The special cellular zone of table cellular zone, Typical stability.
8. device according to claim 7, it is characterised in that also include：
Region query module, the area data for obtaining region to be analyzed.
9. device according to claim 8, it is characterised in that also include：
Computing module is fluctuated, for the attributive character according to each integrated environment zoning, calculates every
The classification fluctuation frequency and ownership probability of individual integrated environment zoning, and frequency is fluctuated according to the classification
Secondary and described ownership probability, the fluctuation situation of zoning after being clustered.
10. device according to claim 9, it is characterised in that the area data bag
Include：The issue of fertility for many years of geo-spatial data, crop in region to be analyzed is according to this and to be analyzed
The environmental data for many years in region.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610221767.6A CN107292755B (en) | 2016-04-11 | 2016-04-11 | Selection and analysis method and device for typical representative area of corn planting environment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610221767.6A CN107292755B (en) | 2016-04-11 | 2016-04-11 | Selection and analysis method and device for typical representative area of corn planting environment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107292755A true CN107292755A (en) | 2017-10-24 |
| CN107292755B CN107292755B (en) | 2020-08-04 |
Family
ID=60095778
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610221767.6A Active CN107292755B (en) | 2016-04-11 | 2016-04-11 | Selection and analysis method and device for typical representative area of corn planting environment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107292755B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102682202A (en) * | 2012-04-28 | 2012-09-19 | 清华大学 | Green space system community planning method for urban and rural ecological planning |
| CN102742462A (en) * | 2011-04-22 | 2012-10-24 | 黑龙江省黑土大森林山产品开发有限公司 | High-cold area winter wheat life limit division and north shift and freezing injury check planting method |
| CN104899786A (en) * | 2015-05-13 | 2015-09-09 | 中国农业大学 | Corn variety planting suitability fine dividing method and system thereof |
| CN104956907A (en) * | 2015-07-26 | 2015-10-07 | 黑龙江超前农业科技开发有限公司 | Beidou navigation global winter wheat satellite zoning expansion and green snow wheat fire belt |
-
2016
- 2016-04-11 CN CN201610221767.6A patent/CN107292755B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102742462A (en) * | 2011-04-22 | 2012-10-24 | 黑龙江省黑土大森林山产品开发有限公司 | High-cold area winter wheat life limit division and north shift and freezing injury check planting method |
| CN102682202A (en) * | 2012-04-28 | 2012-09-19 | 清华大学 | Green space system community planning method for urban and rural ecological planning |
| CN104899786A (en) * | 2015-05-13 | 2015-09-09 | 中国农业大学 | Corn variety planting suitability fine dividing method and system thereof |
| CN104956907A (en) * | 2015-07-26 | 2015-10-07 | 黑龙江超前农业科技开发有限公司 | Beidou navigation global winter wheat satellite zoning expansion and green snow wheat fire belt |
Non-Patent Citations (1)
| Title |
|---|
| 刘哲等: "黄淮海地区玉米品种适宜性精细区划研究", 《农业机械学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107292755B (en) | 2020-08-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Rahman et al. | Remote sensing, spatial multi criteria evaluation (SMCE) and analytical hierarchy process (AHP) in optimal cropping pattern planning for a flood prone area | |
| CN109409261B (en) | Crop classification method and system | |
| Wang et al. | Integrating remotely sensed leaf area index and leaf nitrogen accumulation with RiceGrow model based on particle swarm optimization algorithm for rice grain yield assessment | |
| CN107103378B (en) | Corn planting environment test site layout method and system | |
| CN104899786A (en) | Corn variety planting suitability fine dividing method and system thereof | |
| de Freitas et al. | Land-use and land-cover change processes in the Upper Uruguay Basin: linking environmental and socioeconomic variables | |
| CN109145885B (en) | Remote sensing classification method and system for large-scale crops | |
| Rodriguez et al. | Using entropy theory to improve the definition of homogeneous regions in the semi-arid region of Brazil | |
| CN115618021B (en) | Method and device for recommending planting area suitable for crop variety | |
| CN113592132B (en) | Rainfall objective forecasting method based on numerical weather forecast and artificial intelligence | |
| CN109858180B (en) | Regional crop water consumption space pattern optimization design method | |
| CN102880752B (en) | Optimized design method of regional crop evapotranspiration spatial-temporal pattern | |
| CN113139717B (en) | Crop seedling condition grading remote sensing monitoring method and device | |
| CN117611993B (en) | Method for estimating vegetation classification based on remote sensing actual evapotranspiration | |
| CN107609687B (en) | Crop variety testing station layout method and device | |
| CN115661670B (en) | Method and device for identifying transition region on urban and rural natural landscape gradient | |
| Mirhashemi et al. | Using combined clustering algorithms and association rules for better management of the amount of water delivered to the irrigation network of Abyek Plain, Iran | |
| CN115310719B (en) | Farmland soil sampling scheme design method based on three-stage k-means | |
| CN107292755A (en) | A kind of Analysis on Selecting method and device in corn planting environment Typical Representative area | |
| CN115719453A (en) | Rice planting structure remote sensing extraction method based on deep learning | |
| CN111476434B (en) | GIS-based soil heavy metal fractal dimension spatial variation analysis method | |
| Cadenas et al. | Development of an application to make knowledge available to the farmer: Detection of the most suitable crops for a more sustainable agriculture | |
| Manikandan et al. | Concepts on land evaluation | |
| Saranya et al. | Multi-model ensemble depth adaptive deep neural network for crop yield prediction | |
| Zhao et al. | Risk assessment for comprehensive meteorological disaster of winter wheat based on multisource gridded data: A case study in the upper Huaihe River basin, China |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |