CN107491724B - A kind of Spike Differentiation in Winter Wheat phase recognition methods and device - Google Patents
A kind of Spike Differentiation in Winter Wheat phase recognition methods and device Download PDFInfo
- Publication number
- CN107491724B CN107491724B CN201710515269.7A CN201710515269A CN107491724B CN 107491724 B CN107491724 B CN 107491724B CN 201710515269 A CN201710515269 A CN 201710515269A CN 107491724 B CN107491724 B CN 107491724B
- Authority
- CN
- China
- Prior art keywords
- identified
- winter wheat
- time series
- temperature
- temperature record
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 76
- 230000004069 differentiation Effects 0.000 title claims abstract description 41
- 241000209140 Triticum Species 0.000 claims abstract description 226
- 235000021307 Triticum Nutrition 0.000 claims abstract description 226
- 238000004088 simulation Methods 0.000 claims abstract description 46
- 230000004807 localization Effects 0.000 claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims description 32
- 238000005259 measurement Methods 0.000 claims description 22
- 238000009395 breeding Methods 0.000 claims description 11
- 230000001488 breeding effect Effects 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 10
- 238000013179 statistical model Methods 0.000 claims description 9
- 238000003860 storage Methods 0.000 claims description 9
- 238000004891 communication Methods 0.000 claims description 5
- 230000001052 transient effect Effects 0.000 claims description 5
- 241000208340 Araliaceae Species 0.000 claims description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 3
- 235000008434 ginseng Nutrition 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 25
- 238000010586 diagram Methods 0.000 description 19
- 238000011161 development Methods 0.000 description 8
- 230000008901 benefit Effects 0.000 description 6
- 238000012544 monitoring process Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000007774 longterm Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000000611 regression analysis Methods 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 235000013339 cereals Nutrition 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000012271 agricultural production Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 235000016709 nutrition Nutrition 0.000 description 1
- 230000035764 nutrition Effects 0.000 description 1
- 230000029553 photosynthesis Effects 0.000 description 1
- 238000010672 photosynthesis Methods 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 238000009331 sowing Methods 0.000 description 1
- 230000005068 transpiration Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/188—Vegetation
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- General Health & Medical Sciences (AREA)
- Marketing (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Animal Husbandry (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Primary Health Care (AREA)
- Agronomy & Crop Science (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Multimedia (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the invention provides a kind of Spike Differentiation in Winter Wheat phase recognition methods and devices, inventive method includes: according to the parameter in the winter wheat growing area to be identified being collected into, the model parameter of Wheat Simulation Model is obtained, localization Wheat Simulation Model is obtained；According to winter wheat growing area agricultural meteorological station data to be identified, the first temperature record time series of winter wheat to be identified is obtained, according to the multi-temporal remote sensing image in winter wheat growing area to be identified, obtains the second temperature record time series of winter wheat to be identified；Using data assimilation method, the first temperature record time series and the second temperature record time series are subjected to data assimilation, obtain winter wheat assimilation temperature record time series to be identified；According to assimilation temperature record time series and localization Wheat Simulation Model, the Jointing stage recognition result in winter wheat growing area to be identified is obtained.The embodiment of the present invention realizes the identification of the Spike Differentiation in Winter Wheat phase on regional scale.
Description
Technical field
The present embodiments relate to agricultural technology fields, and in particular to a kind of Spike Differentiation in Winter Wheat phase recognition methods and dress
It sets.
Background technique
A kind of index of the Crop growing stage as reflection crop growth process, for determining that it is important that potential production has
Effect.Useful information and section can be provided for the field management in each stage breeding time work by grasping Crop development process in real time
Foundation is learned, for realizing that the stable yields high receipts of grain are of great significance to.In recent years, with to plant physiological ecology basis
Theoretical further investigation, the crop growth model based on the mechanism process such as crop photosynthesis, breathing, transpiration, nutrition gradually develop
Come, crop growth model by it physiology course and kinetic mechanism, can be compared with accurate simulation crop object with " day "
For the growth and development state of crop on the single-point scale of time step, and the temperature record collection of long-term sequence is driving crop mould
One of the most important input parameter of type operation, directly determines the simulation precision of crop growth situation.
In recent years, global climate and Changes in weather are violent, abnormal and extreme climate takes place frequently, so as to cause more meteorology
Disaster occurs, and has brought tremendous economic losses to agricultural production.Spring frost is exactly Main Agricultural in China's Production of Winter Wheat
One of meteorological disaster, and the Spike Differentiation in Winter Wheat phase is the tricky time that Spring frost occurs, main cause evening is frost damage
When, young fringe subject to damage directly affects " yield forming three elements " (i.e. unit area spike number, number of grain per ear and mass of 1000 kernel).Cause
This, is recognized accurately Tillering of Winter Wheat each period, instructs work for the defence of Spring frost and realizes winter wheat
High and stable yields have great importance.
In the prior art, to the most methods using traditional field experiment of the identification of Spike Differentiation in Winter Wheat phase, have larger
Limitation.Mechanistic crop modeling WheatGrow based on winter wheat physiological and biochemical procedure can be realized with ' day ' as the time
Wheat Jointing stage simulation on the single-point scale of unit.For a long time, the temperature record of driving crop modeling operation mainly leads to
It crosses the long-term observation weather station fixed point being distributed in research area to obtain, there are following both sides defects for this method: on the one hand,
Due to the limitation of economic technology condition, studies the long-term observation weather station limited amount in area and spatial distribution is uneven, see
The temperature record collection measured is only capable of reflecting a certain range of temperature change situation near the website, can not express non-viewing
The temperature change feature of point other than website.Although can be evaluated whether the temperature value of unknown point using the method for space interpolation, by
In by landform, sample point density and distribution and interpolation method etc. influenced, the error of interpolation result is larger.Another party
Face, for crops, its growth and development and yield composition are directly influenced by agricultural microclimate locating for it, and weather station is usual
Be arranged on the spacious vacant lot in urban fringe or suburb, by urban heat land effect and underlying surface influenced, weather station temperature
There are larger differences with farmland temperature.
Therefore, how to propose a kind of scheme, can be realized the identification of Regional Fall Wheat Jointing stage, become urgently to be resolved
Problem.
Summary of the invention
For the defects in the prior art, the embodiment of the invention provides a kind of Spike Differentiation in Winter Wheat phase recognition methods and dresses
It sets.
On the one hand, the embodiment of the invention provides a kind of Spike Differentiation in Winter Wheat phase recognition methods, comprising:
According to the work of geographic factor, meteorologic parameter and winter wheat to be identified in the winter wheat growing area to be identified being collected into
Object parameter obtains the model parameter of Wheat Simulation Model, obtains localization Wheat Simulation Model；
According to the agricultural meteorological station data in the winter wheat growing area preset range to be identified, the winter to be identified is obtained
The first temperature record time series in wheat growth stage, according to the multi-temporal remote sensing shadow in the winter wheat growing area to be identified
Picture obtains the second temperature record time series in the During Growing Period of Winter Wheat to be identified；
Using data assimilation method, by the first temperature record time series and the second temperature record time series
Data assimilation is carried out, the assimilation temperature record time series in the During Growing Period of Winter Wheat to be identified is obtained；
According to the assimilation temperature record time series and the localization Wheat Simulation Model, obtain described wait know
Jointing stage recognition result in other winter wheat growing area.
Further, the agricultural meteorological station data according in the winter wheat growing area preset range to be identified, are obtained
Take the first temperature record time series in the During Growing Period of Winter Wheat to be identified, comprising:
According to the agricultural meteorological station data in the winter wheat growing area preset range to be identified, GIS-Geographic Information System is utilized
Spatial Interpolation Method carries out temperature record interpolation, obtains the first temperature record time series.
Further, described to utilize GIS-Geographic Information System Spatial Interpolation Method, temperature record interpolation is carried out, obtains described first
Temperature record time series, comprising:
Using following formula (1), temperature record interpolation is carried out:
In formula: λiFor weight；N is agricultural weather website number；diFor the reality between interpolation point and i-th of agricultural meteorological station
Ranging from；Z is the temperature estimated value of interpolation point；ZiFor the temperature measured value of i-th (i=1,2,3 ... n) a agricultural meteorological station.
Further, the multi-temporal remote sensing image according in the winter wheat growing area to be identified, obtain it is described to
Identify the second temperature record time series in During Growing Period of Winter Wheat, comprising:
The multi-temporal remote sensing image in the winter wheat growing area to be identified is obtained, according to the multi-temporal remote sensing image,
Obtain the sample Remote Sensing temperature-time sequence data of the winter wheat to be identified within a preset time at default sample point；
Obtain actual measurement temperature time of the winter wheat to be identified at the default sample point in the preset time
Sequence data；
According to the sample Remote Sensing temperature-time sequence data and the actual measurement temperature time series data, institute is obtained
State the second temperature record time series in During Growing Period of Winter Wheat to be identified.
Further, described according to the sample Remote Sensing temperature-time sequence data and the actual measurement temperature time sequence
Column data obtains the second temperature record time series in the During Growing Period of Winter Wheat to be identified, comprising:
According to the sample Remote Sensing temperature-time sequence data and the actual measurement temperature time series data, institute is established
State the Time Series Regression statistical model of winter wheat to be identified；
According to the multi-temporal remote sensing image in the winter wheat growing area to be identified, it is entire to obtain the winter wheat to be identified
Calculating Remote Sensing temperature-time sequence data in breeding time；
According to the calculating Remote Sensing temperature-time sequence data and the Time Series Regression statistical model, institute is obtained
State the second temperature record time series.
Further, described to utilize data assimilation method, by the first temperature record time series and second gas
Warm data time series carry out data assimilation, obtain the assimilation temperature record time sequence in the During Growing Period of Winter Wheat to be identified
Column, comprising:
Using Kalman filtering algorithm, by the first temperature record time series and the second temperature record time sequence
Column carry out sequence assimilation obtains the assimilation temperature record time series.
Further, described to utilize Kalman filtering algorithm, the first temperature record time series and described the are stated by described
Two temperature record time series carry out sequence assimilations obtain the assimilation temperature record time series, comprising:
First temperature record time series and the second temperature record time series are stated for described using following formula (2)
Carry out sequence assimilation obtains the assimilation temperature record time series:
In formula: Yt fFor the predicted value of temperature, XtFor predictor, i.e., the described first temperature record time series, Bt-1It is slotting
Value coefficient vector, RtFor extrapolated value BtError covariance matrix, Ct-1For Bt-1Error covariance matrix, W be dynamic noise error variance
Battle array；σtFor prediction error variance matrix,For the transposed matrix of predictor Xt, V is the error covariance matrix of observation noise；AtTo increase
Beneficial matrix,For σtInverse matrix；Yt oFor the second temperature record time series.
Further, the basis is collected into the geographic factor in winter wheat growing area to be identified, meteorologic parameter and to
The crop parameter for identifying winter wheat, obtains the model parameter of Wheat Simulation Model, comprising:
According to geographic factor, meteorologic parameter and the crop of the winter wheat to be identified ginseng in the winter wheat growing area to be identified
Number, seeks excellent method using loop iteration, obtains the model parameter of the Wheat Simulation Model.
On the other hand, the embodiment of the present invention provides a kind of device for the identification of Spike Differentiation in Winter Wheat phase, comprising:
At least one processor；And
At least one processor being connect with the processor communication, in which:
The memory is stored with the program instruction that can be executed by the processor, and the processor calls described program to refer to
Order is able to carry out above-mentioned Spike Differentiation in Winter Wheat phase recognition methods.
Another aspect, the embodiment of the present invention provide a kind of non-transient computer readable storage medium, which is characterized in that described
Non-transient computer readable storage medium stores computer instruction, and it is small that the computer instruction makes the computer execute the above-mentioned winter
The recognition methods of wheat head idiophase.
Spike Differentiation in Winter Wheat phase recognition methods provided in an embodiment of the present invention is obtained to be identified by satellite remote sensing technology
Winter wheat area remote sensing images obtain the remote sensing temperature-time sequence of winter wheat to be identified, i.e. the second temperature record time series, knot
Two groups of data are utilized data assimilation by the first temperature record time series for closing the winter wheat to be identified that agricultural meteorological station obtains,
Obtain assimilation temperature record time series.Recycle localization Wheat Simulation Model, it can complete winter wheat to be identified
The identification of Spike development.Satellite remote sensing technology is utilized has that the period is short, range is wide in the variation of reflecting regional crop growth
And advantage at low cost, combined data assimilation technique carry out the data and remotely-sensed data that obtain according to agricultural meteorological station
Data assimilation realizes the identification of the Spike Differentiation in Winter Wheat phase on regional scale.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is this hair
Bright some embodiments for those of ordinary skill in the art without creative efforts, can be with root
Other attached drawings are obtained according to these attached drawings.
Fig. 1 is the flow diagram of Spike Differentiation in Winter Wheat phase recognition methods in the embodiment of the present invention；
Fig. 2 is the flow diagram of another Spike Differentiation in Winter Wheat phase recognition methods in the embodiment of the present invention；
Fig. 3 is the highest temperature data processed result contrast schematic diagram that sample point is preset in the embodiment of the present invention；
Fig. 4 is the lowest temperature data processed result contrast schematic diagram that sample point is preset in the embodiment of the present invention；
Fig. 5 (a)-Fig. 5 (g) is Spike Differentiation in Winter Wheat phase recognition result schematic diagram in the embodiment of the present invention；
Fig. 6 is structural schematic diagram of one of the inventive embodiments for the device of Spike Differentiation in Winter Wheat phase identification.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
Every other embodiment obtained without creative efforts, shall fall within the protection scope of the present invention.
Fig. 1 is the flow diagram of Spike Differentiation in Winter Wheat phase recognition methods in the embodiment of the present invention, as shown in Figure 1, this hair
The Spike Differentiation in Winter Wheat phase recognition methods that bright embodiment provides includes:
Geographic factor, meteorologic parameter and the winter wheat to be identified in winter wheat growing area to be identified that S1, basis are collected into
Crop parameter, obtain the model parameter of Wheat Simulation Model, obtain localization Wheat Simulation Model；
Specifically, it is small to use Wheat Simulation Model, that is, WheatGrow crop modeling progress winter for the embodiment of the present invention
The identification model of wheat head idiophase, but Wheat Simulation Model in use, need the model is demarcated, that is, fitted
Close the localization Wheat Simulation Model of winter wheat to be identified.The embodiment of the present invention collects winter wheat growing area to be identified in advance
Interior geographic factor, meteorologic parameter and the crop parameter of winter wheat to be identified carries out wheat growth according to the parameters of acquisition
The calibration of simulation model obtains the localization Wheat Simulation Model for being suitble to winter wheat to be identified.Wherein, geographic factor includes
Longitude and latitude data；Meteorologic parameter mainly includes temperature record, including the highest temperature and lowest temperature data day by day；Crop parameter packet
Include winter wheat variety characterisitic parameter and winter wheat phenological period characteristic etc..Certainly, other needs can also be acquired as needed
Parameter information, the embodiment of the present invention is not specifically limited.
S2, according to the agricultural meteorological station data in the winter wheat growing area preset range to be identified, obtain described wait know
The first temperature record time series in other During Growing Period of Winter Wheat, it is distant according to the multidate in the winter wheat growing area to be identified
Feel image, obtains the second temperature record time series in the During Growing Period of Winter Wheat to be identified；
Specifically, it obtains after being suitble to the localization Wheat Simulation Model of winter wheat to be identified, needs to obtain to be identified
The growth data of winter wheat, to carry out the identification of Spike Differentiation in Winter Wheat phase using the localization Wheat Simulation Model.This
Inventive embodiments obtain the agricultural meteorological station data in winter wheat growing area preset range to be identified, and according to the agricultural got
The first temperature record time series in weather station data acquisition During Growing Period of Winter Wheat to be identified.Wherein, agricultural meteorological station data
The main temperature record day by day including in Winter Wheat-Growing, winter wheat growing area preset range to be identified can be according to practical need
It is configured, the agricultural meteorological station of all agricultural meteorological stations where available winter wheat growing area to be identified in municipality directly under the Central Government
Data, also can according to need the size that setting preset range makes, and the embodiment of the present invention is not specifically limited.Due to general agricultural
The website of weather station is arranged in city suburbs, and spatial distribution is uneven, and agricultural meteorological station data can not accurately represent to be identified
The temperature record of winter wheat.In order to keep the temperature record for getting winter wheat to be identified more accurate, the embodiment of the present invention is utilized
Satellite remote sensing technology obtains more in winter wheat growing area to be identified according to the latitude and longitude information of winter wheat growing area to be identified
Phase remote sensing image obtains the second temperature record time sequence in During Growing Period of Winter Wheat to be identified according to multi-temporal remote sensing image
Column.
Wherein temperature-time sequence refers to that the temperature by the winter wheat area to be identified of acquisition is temporally ranked up, obtain with
Temperature data on the basis of time.
S3, using data assimilation method, by the first temperature record time series and the second temperature record time
Sequence carries out data assimilation, obtains the assimilation temperature record time series in the During Growing Period of Winter Wheat to be identified；
Specifically, the first temperature record time series and the second temperature record time series of winter wheat to be identified are got
Afterwards, the first temperature record time series that will acquire using data assimilation method and the second temperature record time series are carried out
Data assimilation obtains the assimilation temperature record time series in During Growing Period of Winter Wheat to be identified.
S4, according to the assimilation temperature record time series and the localization Wheat Simulation Model, described in acquisition
Jointing stage recognition result in winter wheat growing area to be identified.
Specifically, it after getting the assimilation temperature record time series in During Growing Period of Winter Wheat to be identified, will acquire
Assimilation temperature record time series in During Growing Period of Winter Wheat to be identified, substitute into localization Wheat Simulation Model, obtain to
It identifies the Jointing stage recognition result in winter wheat growing area, that is, obtains the growth and development situation of winter wheat to be identified.
Wherein, the first temperature record time series and the second temperature record time series can be winter wheat plantation to be identified
The daily highest temperature and lowest temperature data in area.
Spike Differentiation in Winter Wheat phase recognition methods provided in an embodiment of the present invention is obtained to be identified by satellite remote sensing technology
Winter wheat area remote sensing images obtain the remote sensing temperature-time sequence of winter wheat to be identified, i.e. the second temperature record time series, knot
Two groups of data are utilized data assimilation by the first temperature record time series for closing the winter wheat to be identified that agricultural meteorological station obtains,
Obtain assimilation temperature record time series.Recycle localization Wheat Simulation Model, it can complete winter wheat to be identified
The identification of Spike development.Satellite remote sensing technology is utilized has that the period is short, range is wide in the variation of reflecting regional crop growth
And advantage at low cost, combined data assimilation technique carry out the data and remotely-sensed data that obtain according to agricultural meteorological station
Data assimilation realizes the identification of the Spike Differentiation in Winter Wheat phase on regional scale.
On the basis of the above embodiments, the agriculture gas according in the winter wheat growing area preset range to be identified
As data of standing, the first temperature record time series in the During Growing Period of Winter Wheat to be identified is obtained, comprising:
According to the agricultural meteorological station data in the winter wheat growing area preset range to be identified, GIS-Geographic Information System is utilized
Spatial Interpolation Method carries out temperature record interpolation, obtains the first temperature record time series.
Specifically, the agricultural meteorological station number in winter wheat growing area preset range to be identified is got in the embodiment of the present invention
According to wherein agricultural meteorological station data include temperature record.Obtain the agriculture gas in winter wheat growing area preset range to be identified
As collected temperature record of standing carries out temperature record difference, i.e., by agricultural weather using GIS-Geographic Information System Spatial Interpolation Method
Acquisition temperature record of standing is converted into the corresponding temperature data in winter wheat area to be identified, obtains the first temperature number of winter wheat to be identified
According to time series.
On the basis of the above embodiments, described to utilize GIS-Geographic Information System Spatial Interpolation Method, temperature record interpolation is carried out,
Obtain the first temperature record time series, comprising:
Using following formula (1), temperature record interpolation is carried out:
In formula: λiFor weight；N is agricultural weather website number；diFor the reality between interpolation point and i-th of agricultural meteorological station
Ranging from；Z is the temperature estimated value of interpolation point；ZiFor the temperature measured value of i-th (i=1,2,3 ... n) a agricultural meteorological station.
Pair specifically, GIS-Geographic Information System Spatial Interpolation Method is used in the embodiment of the present invention, carries out temperature record difference, i.e.,
Each collected temperature record of agricultural meteorological station generates the temperature record time day by day using inverse distance weighted interpolation method (IDW) interpolation
Sequence, specific method can carry out temperature record difference using above-mentioned formula (1).Firstly, according to interpolation point and each agriculture gas
As the measured distance between station obtains a weighted value λi, further according to the weighted value λ of acquisitioniWith the temperature of each agricultural meteorological station
Measured value successively estimates temperature estimated value Z of each agricultural meteorological station relative to interpolation point, that is, obtains the of winter wheat to be identified
One temperature record time series.Wherein, interpolation point is the position where winter wheat area to be identified, it can is existed as needed
Interpolation point is arranged in winter wheat area to be identified, and the number of interpolation point, which also can according to need, to be configured, and the embodiment of the present invention is not done
Specific to limit, the measured distance between interpolation point and each agricultural meteorological station can be obtained according to latitude and longitude coordinates.
Spike Differentiation in Winter Wheat phase recognition methods provided in an embodiment of the present invention, since general agricultural meteorological station is all disposed within city
The suburb in city, the temperature record that agricultural meteorological station obtains cannot indicate the temperature record in winter wheat area to be identified.The present invention is real
Example is applied by obtaining the agricultural meteorological station data in winter wheat preset range to be identified, using GIS-Geographic Information System space interpolation
Method carries out temperature record difference for the collected temperature record of agricultural meteorological station and is converted into the temperature in winter wheat area to be identified
Data obtain the first temperature record time series of winter wheat to be identified.Improve Spike Differentiation in Winter Wheat phase identification data used
Accuracy, further improve the accuracy of Spike Differentiation in Winter Wheat phase recognition result.
On the basis of the above embodiments, the multi-temporal remote sensing shadow according in the winter wheat growing area to be identified
Picture obtains the second temperature record time series in the During Growing Period of Winter Wheat to be identified, comprising:
The multi-temporal remote sensing image in the winter wheat growing area to be identified is obtained, according to the multi-temporal remote sensing image,
Obtain the sample Remote Sensing temperature-time sequence data of the winter wheat to be identified within a preset time at default sample point；
Obtain actual measurement temperature time of the winter wheat to be identified at the default sample point in the preset time
Sequence data；
According to the sample Remote Sensing temperature-time sequence data and the actual measurement temperature time series data, institute is obtained
State the second temperature record time series in During Growing Period of Winter Wheat to be identified.
Specifically, since agricultural meteorological station data can not accurately indicate the temperature record in winter wheat area to be identified, this hair
Bright embodiment obtains the remote sensing image data in winter wheat growing area to be identified, according to the remote sensing shadow of acquisition by remote sensing technology
As data, During Growing Period of Winter Wheat to be identified within a preset time at the default sample point in winter wheat growing area to be identified is obtained
Interior sample Remote Sensing temperature-time sequence data, that is, MODIS LST time series data.MODIS is that Terra and Aqua are defended
One of main sensors carried on star, LST are a kind of surface temperature products on MODIS, by satellite remote sensing technology, specifically
The MODIS Aqua surface temperature product MYD11A1 that multidate can be used obtains MODIS LST time series data i.e. sample
Remote Sensing temperature-time sequence data.Due to being the earth's surface temperature of winter wheat growing area to be identified by satellite remote sensing technology acquisition
Degree evidence needs to convert surface temperature data to the temperature record that can indicate winter wheat to be identified.The embodiment of the present invention is logical
It crosses and obtains the actual measurement temperature time series data of the winter wheat to be identified at default sample point within a preset time, it specifically can be with
The temperature above winter wheat to be identified is surveyed by ground, as: it can be by winter wheat growing area to be identified presets sample point
Air temperature sensor is set at top preset height such as 2m and observes acquisition actual measurement temperature time series data.Further according to acquisition
Sample Remote Sensing temperature-time sequence data and actual measurement temperature time series data, obtain in During Growing Period of Winter Wheat to be identified
Second temperature record time series.
Such as: multiple default sample points are set in winter wheat growing area to be identified in advance, are obtained by satellite remote sensing technology
The satellite remote sensing images for obtaining winter wheat kind to be identified are obtained in winter wheat growing area to be identified according to the satellite remote sensing images and are preset
The MODIS LST time series data of winter wheat to be identified within a preset time at sample point, such as obtain default sample point
When the MODIS LST of winter wheat to be identified daily 13:30 on March 12nd, 2016 to May 31 at place or so and evening 01:30 or so
Between sequence data, certain preset time, which can according to need, to be set, and the embodiment of the present invention is not especially limited.Again wait know
At default sample point in other winter wheat growing area be arranged air temperature sensor, the air temperature sensor can be set away from
Ground level is to acquire default sample point within a preset time such as by air temperature sensor: in March, 2016 at the position of 2m
Actual measurement temperature time series data to daily 13:30 on May 31 or so and evening 01:30 or so on the 12.According to what is got
MODIS LST time series data and actual measurement temperature time series data, obtain the second gas in During Growing Period of Winter Wheat to be identified
Warm data time series.
On the basis of the above embodiments, described according to the sample Remote Sensing temperature-time sequence data and the reality
Temperature time series data is surveyed, the second temperature record time series in the During Growing Period of Winter Wheat to be identified is obtained, comprising:
According to the sample Remote Sensing temperature-time sequence data and the actual measurement temperature time series data, institute is established
State the Time Series Regression statistical model of winter wheat to be identified；
According to the multi-temporal remote sensing image in the winter wheat growing area to be identified, it is entire to obtain the winter wheat to be identified
Calculating Remote Sensing temperature-time sequence data in breeding time；
According to the calculating Remote Sensing temperature-time sequence data and the Time Series Regression statistical model, institute is obtained
State the second temperature record time series.
Specifically, the sample Remote Sensing temperature-time sequence preset at sample point in planting winter wheat area to be identified is being got
After column data and actual measurement temperature time series data, to the sample Remote Sensing temperature-time sequence data and actual measurement temperature time
Sequence data carries out regression analysis, establishes the Time Series Regression statistical model of winter wheat to be identified.According to winter wheat kind to be identified
Multi-temporal remote sensing image in growing area obtains the calculating Remote Sensing temperature-time sequence in winter wheat to be identified entire breeding time
Data, by the calculating Remote Sensing temperature-time sequence data input time serial regression statistical model, it can will be to be identified
The Remote Sensing temperature inversion of winter wheat is the air themperature of winter wheat to be identified, obtains the second temperature number of winter wheat to be identified
According to time series.
Spike Differentiation in Winter Wheat phase recognition methods provided in an embodiment of the present invention obtains the winter to be identified using satellite remote sensing technology
The MODIS LST time series data of wheat, then by regression analysis, convert MODIS LST time series data to wait know
The second temperature record time series of other winter wheat.Satellite remote sensing technology has the characteristics that macroscopic view, quick and dynamic, in reflection area
There is the advantage that the period is short, range is wide and at low cost in the variation of domain crop growth.The embodiment of the present invention utilizes satellite remote sensing
Technology obtains the second temperature record time series of winter wheat to be identified, provides for the subsequent identification for carrying out the Spike Differentiation in Winter Wheat phase
Better data basis improves the accuracy of Spike Differentiation in Winter Wheat phase identification data and recognition result.
On the basis of the above embodiments, described to utilize data assimilation method, by the first temperature record time series
Data assimilation is carried out with the second temperature record time series, obtains the assimilation temperature in the During Growing Period of Winter Wheat to be identified
Data time series, comprising:
Using Kalman filtering algorithm, by the first temperature record time series and the second temperature record time sequence
Column carry out sequence assimilation obtains the assimilation temperature record time series.
Specifically, in the first temperature record time series for getting winter wheat to be identified by agricultural meteorological station, and it is logical
It crosses after satellite remote sensing technology gets the second temperature record time sequence of winter wheat to be identified, the embodiment of the present invention utilizes Kalman
Filtering algorithm, the first temperature record time series that will acquire and the sequence carry out sequence assimilation of the second temperature record time.Press
According to the corresponding relationship of time, the first temperature record time series and the second temperature record time sequence are assimilated, assimilated
Temperature record time series.
On the basis of the above embodiments, described to utilize Kalman filtering algorithm, stated for the first temperature record time for described
Sequence and the second temperature record time series carry out sequence assimilation obtain the assimilation temperature record time series, comprising:
First temperature record time series and the second temperature record time series are stated for described using following formula (2)
Carry out sequence assimilation obtains the assimilation temperature record time series:
In formula: Yt fFor the predicted value of temperature, XtFor predictor, i.e., the described first temperature record time series, Bt-1It is slotting
Value coefficient vector, RtFor extrapolated value BtError covariance matrix, Ct-1For Bt-1Error covariance matrix, W be dynamic noise error variance
Battle array；σtFor prediction error variance matrix,For predictor XtTransposed matrix, V be observation noise error covariance matrix；AtTo increase
Beneficial matrix,For σtInverse matrix；Yt oFor the second temperature record time series.
Specifically, the embodiment of the present invention is getting the first temperature record time series and the second temperature record time sequence
Afterwards, Kalman filtering (Kalman Fileter, KF) algorithm, principle is divided into prediction and updates two steps, specific to utilize card
Kalman Filtering algorithm carries out data assimilation using above-mentioned formula (2).Yt f=XtBt-1For prognostic equation, Yt fFor the predicted value of temperature,
XtFor predictor, i.e., meteorological site temperature value in the During Growing Period of Winter Wheat period to be identified, that is, when the first temperature record
Between sequence；Bt-1For interpolation coefficient vector, i.e., the weighted value λ of each meteorological site in the anti-distance weighting of temperature (IDW) interpolation algorithm；
RtFor extrapolated value BtError covariance matrix, Ct-1For Bt-1Error covariance matrix, W be dynamic noise error covariance matrix；σtIt is missed for forecast
Poor variance matrix,For predictor XtTransposed matrix, V be observation noise error covariance matrix；AtFor gain matrix,For
σtInverse matrix；Yt oFor the second temperature record time series, that is, utilize MODIS surface temperature product (LST) through regression equation meter
Temperature value above the entire breeding time wheatland of obtained winter wheat to be identified.
By the above method, after the assimilation temperature record time series for obtaining winter wheat to be identified, temperature record will be assimilated
Time series is input to localization Wheat Simulation Model, realizes the real-time monitoring of winter wheat to be identified, obtains the winter to be identified
The recognition result of wheat Jointing stage.
Satellite remote sensing technology has the characteristics that macroscopic view, quick and dynamic, has in the variation of reflecting regional crop growth
Have that the period is short, range is wide and advantage at low cost, but it is due to the limitation of the factors such as spatial resolution, is also difficult to really reflect and make
The relationship of object growth and development and meteorological factor.Data assimilation passes through coupling remote sensing technology and GIS-Geographic Information System space interpolation
Technology rationally can effectively estimate temperature record above winter wheat growing area to be identified, to drive crop in conjunction with the two advantage
The real-time monitoring of model realization crop growth situation is provided for correct cultivation step of formulating with the high and stable yields for realizing crop
Decision support.
Spike Differentiation in Winter Wheat phase recognition methods provided in an embodiment of the present invention, by acquiring winter wheat growing area pair to be identified
Two groups of data are carried out data assimilation, base using kalman filter method by the agricultural meteorological station data and satellite remote sensing date answered
Accurately identifying for the Spike Differentiation in Winter Wheat phase on regional scale is realized in the method for remote sensing and data assimilation, with traditional Spike development
Phase recognition methods is compared, easy to operate, while having saved time and manpower and material resources, is winter wheat growing way and meteorological disaster monitoring
Reference is provided with administrative decision, while providing the new thinking solved the problems, such as to improve crop modeling simulation precision.
The geographic factor in winter wheat growing area to be identified that on the basis of the above embodiments, the basis is collected into,
The crop parameter of meteorologic parameter and winter wheat to be identified, obtains the model parameter of Wheat Simulation Model, comprising:
According to geographic factor, meteorologic parameter and the crop of the winter wheat to be identified ginseng in the winter wheat growing area to be identified
Number, seeks excellent method using loop iteration, obtains the model parameter of the Wheat Simulation Model.
Specifically, in order to make Wheat Simulation Model i.e. WheatGrow crop modeling be more suitable winter wheat to be identified
Monitor environment, the embodiment of the present invention according in winter wheat growing area to be identified geographic factor, meteorologic parameter and the winter to be identified it is small
The crop parameter of wheat asks the calibration of excellent method progress model using loop iteration.It specifically can be in winter wheat growing area to be identified
The continuous 10 years highest temperatures day by day and lowest temperature data that National primary standard weather station obtains are that driving data combines the region
The phenological period data of winter wheat completes the localization calibration of model parameter, obtains the model parameter of Wheat Simulation Model, builds
It is vertical to be suitble to local localization Wheat Simulation Model.
Fig. 2 is the flow diagram of another Spike Differentiation in Winter Wheat phase recognition methods in the embodiment of the present invention, as shown in Fig. 2,
The embodiment of the present invention mainly passes through the day highest temperature data and day minimum gas that agricultural meteorological station obtains winter wheat growing area to be identified
Warm data obtain interpolation lowest temperature time series (on face) and interpolation highest temperature time sequence by inverse distance weighted interpolation method
It arranges in (on face), that is, obtains the first temperature record time series.Simultaneously.It is obtained by satellite remote sensing technology MODIS LST to be identified
The remote sensing images of winter wheat growing area, and remote sensing images are pre-processed, obtain the surface temperature time of winter wheat to be identified
Sequence data (MYD11A1 time series data) is divided into surface temperature time series on daytime (MYD11A1-day time series number
According to) and nighttime surface temperature-time sequence (MYD11A1-night time series data).Again by regression analysis, earth's surface is obtained
The regression equation of temperature-time sequence and actual measurement temperature time series data, is divided into highest temperature surface temperature time series and actual measurement
The regression equation and lowest temperature surface temperature time series of temperature time series data and returning for actual measurement temperature time series data
Return equation.The the second temperature record time series for further obtaining winter wheat to be identified is divided into the lowest temperature time above wheatland
Highest temperature time series above sequence (on face) and wheatland (on face).Using Kalman filtering algorithm, data assimilation is carried out, is obtained
The lowest temperature time series after highest temperature time series (on face) and assimilation after must assimilating (on face).By what is be obtained ahead of time
The crop parameter of winter wheat to be identified such as winter wheat variety genetic parameter, assimilation after highest temperature time series (on face) with
And the lowest temperature time series (on face) after assimilation is input to localization Wheat Simulation Model WheatGrow crop modeling
In, the Spike Differentiation in Winter Wheat recognition result on acquisition face.
Carry out the technology of the present invention is further explained embodiment with the identification of Henan Province's Shangqiu City Spike Differentiation in Winter Wheat phase below
Scheme:
Selecting ripe main breed Wenmai 6 in Shangqiu City winter wheat growing area semi-winterness is standard crop, is based on
WheatGrow crop modeling asks excellent method to solve calibration Cultivar parameter value using loop iteration.Specific steps are as follows: (1) according to Shangqiu
Each annual winter wheat phenological observation data of base station 2005-2006 to 2014-2015, determines that winter wheat universal sowing date is
October 15, date of generally emerging are October 21, and heading stage is April 17, and it is 178d that emergence to heading, which averagely lasts number of days,
In this, as emergence to the standard number of days of heading；(2) according to related literatures, temperature sensitivity TS, the life of Wenmai 6 are determined
Manage vernalization time PVT, the value range of photoperiod sensitivity PS and basic prematureness IE be respectively [1.4,1.5], [20,25],
[0.004,0.005] and [0.80,0.85], and the circulation step-length that 4 parameters are respectively set is 0.01,1,0.0001 and 0.01；
(3) based on each annual diurnal meterorological data of the above parameter and Shangqiu base station 2005-2006 to 2014-2015 and (1) step
Determining emergence date, Cultivar parameter value range and step-length, loop iteration run WheatGrow model, extract each run
Calculation average time is sought after the emergence of acquisition to the number of days of heading, and the number of days is compared with standard number of days, by difference minimum
When Cultivar parameter combination be determined as the combination of optimal Cultivar parameter, obtained the value difference of tetra- parameters of TS, PVT, PS and IE
It is 1.5,25,0.0048 and 0.80, that is, completes the single-point scale calibration of WheatGrow crop modeling, obtain localization wheat
Growth simulation model.
Based on 2015-2016 year Shangqiu City is local and periphery in totally 21 agricultural weather website Winter Wheat-Growings by
Daily temperature data obtain region daily maximum temperature and daily minimal tcmperature time sequence using inverse distance weighted interpolation method (IDW) interpolation
Column data obtains the first temperature record time series of winter wheat growing area, can specifically be carried out using above-mentioned formula (1) slotting
Value calculates.
Sample point is set up in the winter wheat contiguous plant area of the town Shangqiu City Shuan Ba, obtains wheatland by setting field sensors
It surveys temperature record and surveys temperature time series data in top.It extracts in observation period (12 days~May 31 March in 2016)
MODIS LST product daytime 13:30 or so and night 01:30 or so time series data collection, obtained respectively with field sensing
The daily maximum temperature of 2m, the lowest temperature carry out linear regression analysis and settling time serial regression statistics mould above the wheatland taken
Type.MODIS LST product MYD11A1 13:30 on daytime of winter wheat entire breeding time or so and night 01 are extracted by pixel unit:
30 or so time series data collection is calculated in winter wheat growing area using the Time Series Regression statistical model of foundation
2m daily maximum temperature and daily minimal tcmperature data above wheatland, that is, obtain the second temperature record time sequence of winter wheat growing area
Column.
Agricultural meteorological station data will be passed through) interpolation obtains region daily maximum temperature and daily minimal tcmperature time series data
(the first temperature record time series), and the daily maximum temperature and daily minimal tcmperature data (that are obtained by satellite remote sensing technology
Two temperature record time serieses), using Kalman filtering carry out sequence assimilation, can specifically be counted using above-mentioned formula (2)
According to assimilation, the assimilation temperature record time series in During Growing Period of Winter Wheat to be identified is obtained.Fig. 3 is to preset in the embodiment of the present invention
The highest temperature data processed result contrast schematic diagram of sample point, Fig. 4 are the minimum gas that sample point is preset in the embodiment of the present invention
Warm data processed result contrast schematic diagram has the data obtained by interpolation, by default sample as shown in Figure 3 and Figure 4, in figure
This is put the data that air temperature sensor is arranged at high 2m and obtains, is converted by the surface temperature data that satellite remote sensing technology obtains
Obtained by the data that obtain of temperature data, that is, LST conversion data and data assimilation.As shown in Figure 3 and Figure 4, the present invention is implemented
The assimilation temperature record time series finally chosen in example, comprehensively considered the temperature record obtained by agricultural meteorological station and
Surface temperature data are obtained by satellite remote sensing, improve the accuracy of data acquisition.
The assimilation temperature record time series of acquisition is inputted into WheatGrow crop modeling and localization wheat growth simulation
The Spike Differentiation in Winter Wheat phase recognition result of Shangqiu City can be obtained in model, moving model.Fig. 5 (a)-Fig. 5 (g) is that the present invention is implemented
Spike Differentiation in Winter Wheat phase recognition result schematic diagram in example, Jointing stage Start Date indicate that Fig. 5 (a) is with number of days after emergence
The schematic diagram of wheat list rib phase, 65d and 75d is respectively indicated after emergence 65 days and 75 days in figure, and Fig. 5 (b) is the two rib phase of wheat
Schematic diagram, 85d and 110d is respectively indicated after emergence 85 days and 110 days in figure, Fig. 5 (c) for the wheat small florescence schematic diagram, in figure
127d and 137 is respectively indicated after emergence 127 days and 137 days, and Fig. 5 (d) is the schematic diagram of wheat Pistil And Stamen idiophase, 133d in figure
It is respectively indicated with 143d after emerging 133 days and 143 days, Fig. 5 (e) is schematic diagram of the wheat medical every the phase, 137d and 151d points in figure
137 days and 151 days after Biao Shi not emerging, Fig. 5 (f) is the schematic diagram of four body idiophase of wheat, and 149d and 160d distinguishes table in figure
It shows after seedling 149 days and 160 days, Fig. 5 (g) is the schematic diagram of wheat heading stage, after 160d and 170d respectively indicates emergence in figure
160 days and 170 days.As shown in Fig. 5 (a)-Fig. 5 (g), method through the embodiment of the present invention can obtain winter wheat kind to be identified
Wheat in growing area different zones changes with time, entire breeding time variation, that is, obtain the fringe of winter wheat to be identified
Break up recognition result.
Spike Differentiation in Winter Wheat phase recognition methods provided in an embodiment of the present invention, by remote sensing technology, GIS-Geographic Information System space
Interpolation technique and crop modeling are merged by data assimilation method, and it is simple to realize the Spike Differentiation in Winter Wheat phase on regional scale
Quickly identification, improves the accuracy of Spike Differentiation in Winter Wheat phase recognition result, while having saved time and manpower and material resources, is the winter
Wheat growing way and meteorological disaster monitoring and administrative decision provide reference, while providing for raising crop modeling simulation precision new
The thinking solved the problems, such as.
Fig. 6 is structural schematic diagram of one of the inventive embodiments for the device of Spike Differentiation in Winter Wheat phase identification, such as Fig. 6
It is shown, the apparatus may include: processor (processor) 61, memory (memory) 62 and communication bus 63, wherein
Processor 61, memory 62 complete mutual communication by communication bus 63.Processor 61 can call in memory 62
Logical order, to execute following method: according in the winter wheat growing area to be identified being collected into geographic factor, meteorologic parameter and
The crop parameter of winter wheat to be identified obtains the model parameter of Wheat Simulation Model, obtains localization wheat growth simulation
Model；According to the agricultural meteorological station data in the winter wheat growing area preset range to be identified, it is small to obtain the winter to be identified
The first temperature record time series in wheat breeding time, according to the multi-temporal remote sensing shadow in the winter wheat growing area to be identified
Picture obtains the second temperature record time series in the During Growing Period of Winter Wheat to be identified；It, will be described using data assimilation method
First temperature record time series and the second temperature record time series carry out data assimilation, and it is small to obtain the winter to be identified
Assimilation temperature record time series in wheat breeding time；According to the assimilation temperature record time series and the localization wheat
Growth simulation model obtains the Jointing stage recognition result in the winter wheat growing area to be identified.
In addition, the logical order in above-mentioned memory 62 can be realized and as only by way of SFU software functional unit
Vertical product when selling or using, can store in a computer readable storage medium.Based on this understanding, this hair
Substantially the part of the part that contributes to existing technology or the technical solution can be with soft in other words for bright technical solution
The form of part product embodies, which is stored in a storage medium, including some instructions are to make
It obtains a computer equipment (can be personal computer, server or the network equipment etc.) and executes each embodiment of the present invention
The all or part of the steps of the method.And storage medium above-mentioned include: USB flash disk, mobile hard disk, read-only memory (ROM,
Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic or disk etc. it is various
It can store the medium of program code.
The embodiment of the present invention provides a kind of non-transient computer readable storage medium, the non-transient computer readable storage
Medium storing computer instruction, the computer instruction make the computer execute side provided by above-mentioned each method embodiment
Method, for example, according to geographic factor, meteorologic parameter and the winter wheat to be identified in the winter wheat growing area to be identified being collected into
Crop parameter, obtain the model parameter of Wheat Simulation Model, obtain localization Wheat Simulation Model；According to described
Agricultural meteorological station data in winter wheat growing area preset range to be identified obtain the in the During Growing Period of Winter Wheat to be identified
One temperature record time series obtains described wait know according to the multi-temporal remote sensing image in the winter wheat growing area to be identified
The second temperature record time series in other During Growing Period of Winter Wheat；Using data assimilation method, when by first temperature record
Between sequence and the second temperature record time series carry out data assimilation, obtain same in the During Growing Period of Winter Wheat to be identified
Change temperature record time series；According to the assimilation temperature record time series and the localization Wheat Simulation Model,
Obtain the Jointing stage recognition result in the winter wheat growing area to be identified.
The above examples are only used to illustrate the technical scheme of the present invention, rather than its limitations；Although with reference to the foregoing embodiments
Invention is explained in detail, those skilled in the art should understand that: it still can be to aforementioned each implementation
Technical solution documented by example is modified or equivalent replacement of some of the technical features；And these are modified or replace
It changes, the spirit and scope for technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution.
Claims (5)
1. a kind of Spike Differentiation in Winter Wheat phase recognition methods characterized by comprising
According to geographic factor, meteorologic parameter and the crop of the winter wheat to be identified ginseng in the winter wheat growing area to be identified being collected into
Number obtains the model parameter of Wheat Simulation Model, obtains localization Wheat Simulation Model；
According to the agricultural meteorological station data in winter wheat growing area preset range to be identified, obtain in During Growing Period of Winter Wheat to be identified
The first temperature record time series the winter to be identified is obtained according to the multi-temporal remote sensing image in winter wheat growing area to be identified
The second temperature record time series in wheat growth stage；
Using data assimilation method, the first temperature record time series and the second temperature record time series are carried out
Data assimilation obtains the assimilation temperature record time series in During Growing Period of Winter Wheat to be identified；
According to the assimilation temperature record time series and the localization Wheat Simulation Model, winter wheat to be identified is obtained
Jointing stage recognition result in growing area；
Wherein, the agricultural meteorological station data according in winter wheat growing area preset range to be identified, it is small to obtain the winter to be identified
The first temperature record time series in wheat breeding time, comprising:
According to the agricultural meteorological station data in winter wheat growing area preset range to be identified, GIS-Geographic Information System space interpolation is utilized
Method carries out temperature record interpolation, obtains the first temperature record time series,
It is described to utilize GIS-Geographic Information System Spatial Interpolation Method, temperature record interpolation is carried out, the first temperature record time is obtained
Sequence, comprising:
Using following formula (1), temperature record interpolation is carried out:
In formula: λiFor weight；N is agricultural weather website number；diActual measurement between interpolation point and i-th of agricultural meteorological station away from
From；Z is the temperature estimated value of interpolation point；ZiFor the temperature measured value of i-th (i=1,2,3 ... n) a agricultural meteorological station；
The multi-temporal remote sensing image according in winter wheat growing area to be identified obtains in During Growing Period of Winter Wheat to be identified
Two temperature record time serieses, specifically include:
According to remote sensing image data when more in winter wheat growing area to be identified, obtain default in winter wheat growing area to be identified
Sample Remote Sensing temperature-time sequence data in During Growing Period of Winter Wheat to be identified at sample point within a preset time, then lead to
It crosses ground actual measurement and obtains winter wheat actual measurement temperature time series data to be identified, according to the sample Remote Sensing temperature-time sequence
Column data and the winter wheat to be identified survey temperature time series data, obtain the in the During Growing Period of Winter Wheat to be identified
Two temperature record time serieses；
It is described according to the sample Remote Sensing temperature-time sequence data and the actual measurement temperature time series data, obtain institute
State the second temperature record time series in During Growing Period of Winter Wheat to be identified, comprising:
According to the sample Remote Sensing temperature-time sequence data and the actual measurement temperature time series data, establish to be identified
The Time Series Regression statistical model of winter wheat；
According to the multi-temporal remote sensing image in the winter wheat growing area to be identified, obtain in winter wheat to be identified entire breeding time
Calculating Remote Sensing temperature-time sequence data；
According to the calculating Remote Sensing temperature-time sequence data and the Time Series Regression statistical model, described the is obtained
Two temperature record time serieses；
It is described to utilize data assimilation method, by the first temperature record time series and the second temperature record time series
Data assimilation is carried out, is specifically included:
First temperature record time series and the second temperature record time series progress are stated for described using following formula (2)
Sequence is assimilated, and the assimilation temperature record time series is obtained:
In formula: Yt fFor the predicted value of temperature, XtFor predictor, i.e., the described first temperature record time series, Bt-1For interpolation system
Number vector, Ct-1For Bt-1Error covariance matrix, W be dynamic noise error covariance matrix, RtFor extrapolated value BtError covariance matrix and dynamic
The sum of state noise error variance matrix；σtFor prediction error variance matrix,For predictor XtTransposed matrix, V is observation noise
Error covariance matrix；AtFor gain matrix,For σtInverse matrix；Yt oFor the second temperature record time series；CtFor Bt
Error covariance matrix.
2. the method according to claim 1, wherein described utilize data assimilation method, by first temperature
Data time series and the second temperature record time series carry out data assimilation, obtain the During Growing Period of Winter Wheat to be identified
Interior assimilation temperature record time series, comprising:
Using Kalman filtering algorithm, by the first temperature record time series and the second temperature record time series into
Row sequence is assimilated, and the assimilation temperature record time series is obtained.
3. -2 described in any item methods according to claim 1, which is characterized in that the winter wheat to be identified that the basis is collected into
The crop parameter of geographic factor, meteorologic parameter and winter wheat to be identified in growing area, obtains the mould of Wheat Simulation Model
Shape parameter, comprising:
According to geographic factor, meteorologic parameter and the crop parameter of winter wheat to be identified in the winter wheat growing area to be identified,
Excellent method is asked using loop iteration, obtains the Wheat Simulation Model.
4. a kind of device for the identification of Spike Differentiation in Winter Wheat phase characterized by comprising
At least one processor；And
At least one processor being connect with the processor communication, in which:
The memory is stored with the program instruction that can be executed by the processor, and the processor calls described program to instruct energy
Enough execute method as described in any one of claims 1 to 3.
5. a kind of non-transient computer readable storage medium, which is characterized in that the non-transient computer readable storage medium is deposited
Computer instruction is stored up, the computer instruction makes the computer execute method as described in any one of claims 1 to 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710515269.7A CN107491724B (en) | 2017-06-29 | 2017-06-29 | A kind of Spike Differentiation in Winter Wheat phase recognition methods and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710515269.7A CN107491724B (en) | 2017-06-29 | 2017-06-29 | A kind of Spike Differentiation in Winter Wheat phase recognition methods and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107491724A CN107491724A (en) | 2017-12-19 |
| CN107491724B true CN107491724B (en) | 2019-11-26 |
Family
ID=60643471
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710515269.7A Active CN107491724B (en) | 2017-06-29 | 2017-06-29 | A kind of Spike Differentiation in Winter Wheat phase recognition methods and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107491724B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108280410B (en) * | 2018-01-10 | 2020-10-27 | 北京农业信息技术研究中心 | Crop identification method and system based on binary coding |
| CN110751094B (en) * | 2019-10-21 | 2020-08-28 | 北京师范大学 | Crop yield estimation method based on GEE comprehensive remote sensing image and deep learning method |
| CN111296201B (en) * | 2020-03-05 | 2023-04-07 | 山西省农业科学院棉花研究所 | Method for field identification and evaluation of late frost freezing resistance of wheat |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105809148A (en) * | 2016-03-29 | 2016-07-27 | 中国科学院遥感与数字地球研究所 | Crop drought recognition and risk evaluation method based on remote sensing time-space-spectrum fusion |
-
2017
- 2017-06-29 CN CN201710515269.7A patent/CN107491724B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105809148A (en) * | 2016-03-29 | 2016-07-27 | 中国科学院遥感与数字地球研究所 | Crop drought recognition and risk evaluation method based on remote sensing time-space-spectrum fusion |
Non-Patent Citations (3)
| Title |
|---|
| Estimating net radiation flux in the Tibetan Plateau by assimilating MODIS LST products with an ensemble Kalman filter and particle filter;Zhen Li 等;《International Journal of Applied Earth Observation and Geoinformation》;20121031;第19卷;第1-11页 * |
| 基于时间序列LAI和ET同化的冬小麦遥感估产方法比较;黄健熙 等;《农业工程学报》;20150228;第31卷(第4期);第197-203页 * |
| 基于集合Kalman滤波的作物生长模型与遥感数据的同化研究;范丽颖;《中国优秀硕士学位论文全文数据库农业科技辑》;20121015(第10期);第1-12,27-30页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107491724A (en) | 2017-12-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109829234B (en) | A kind of across scale Dynamic High-accuracy crop condition monitoring and yield estimation method based on high-definition remote sensing data and crop modeling | |
| CN103034910B (en) | Based on the regional scale disease and pest Forecasting Methodology of multi-source information | |
| CN111368736B (en) | Rice refined estimation method based on SAR and optical remote sensing data | |
| Zhang et al. | Study on factors affecting corn yield based on the Cobb-Douglas production function | |
| WO2016040662A1 (en) | Methods and systems for managing crop harvesting activities | |
| CN101858971A (en) | Rice yield remote sensing estimation method based on MODIS data | |
| JP7176980B2 (en) | Crop management system and crop management method | |
| Bannayan et al. | A stochastic modelling approach for real-time forecasting of winter wheat yield | |
| Wang et al. | Increased uncertainty in simulated maize phenology with more frequent supra-optimal temperature under climate warming | |
| US11716985B2 (en) | Method for remediating developmentally delayed plants | |
| US20150185196A1 (en) | Characterization of field sites for utility in agronomic stress trials | |
| CN111316802B (en) | Crop water and fertilizer monitoring method, device, equipment and medium based on large-scale sprinkler | |
| Zhu et al. | Optimization of multi-source UAV RS agro-monitoring schemes designed for field-scale crop phenotyping | |
| CN107491724B (en) | A kind of Spike Differentiation in Winter Wheat phase recognition methods and device | |
| Amankulova et al. | Time-series analysis of Sentinel-2 satellite images for sunflower yield estimation | |
| CN113223040A (en) | Remote sensing-based banana yield estimation method and device, electronic equipment and storage medium | |
| Yang et al. | A VI-based phenology adaptation approach for rice crop monitoring using UAV multispectral images | |
| CN112819227A (en) | County scale winter wheat yield per unit prediction method and system | |
| Belal et al. | Precision farming technologies to increase soil and crop productivity | |
| Amankulova et al. | Sunflower crop yield prediction by advanced statistical modeling using satellite-derived vegetation indices and crop phenology | |
| JP2022082636A (en) | Information processing device | |
| Plaščak et al. | An overview of precision irrigation systems used in agriculture | |
| US20130096832A1 (en) | Method and system for partitioning farming field | |
| Peng et al. | Improving the spatial and temporal estimation of maize daytime net ecosystem carbon exchange variation based on unmanned aerial vehicle multispectral remote sensing | |
| Kushwaha et al. | Survey of impact of technology on effective implementation of precision farming in India |
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 |