CN107392503A - A kind of appraisal procedure of corn Climatic regionalization risk - Google Patents
A kind of appraisal procedure of corn Climatic regionalization risk Download PDFInfo
- Publication number
- CN107392503A CN107392503A CN201710712114.2A CN201710712114A CN107392503A CN 107392503 A CN107392503 A CN 107392503A CN 201710712114 A CN201710712114 A CN 201710712114A CN 107392503 A CN107392503 A CN 107392503A
- Authority
- CN
- China
- Prior art keywords
- corn
- temperature
- risk
- heat evil
- climatic regionalization
- 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
- 235000002017 Zea mays subsp mays Nutrition 0.000 title claims abstract description 127
- 240000008042 Zea mays Species 0.000 title claims abstract description 124
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 title claims abstract description 118
- 235000005822 corn Nutrition 0.000 title claims abstract description 118
- 230000007261 regionalization Effects 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000002159 abnormal effect Effects 0.000 claims abstract description 39
- 238000009826 distribution Methods 0.000 claims abstract description 24
- 241001269238 Data Species 0.000 claims abstract description 13
- 230000008569 process Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000011156 evaluation Methods 0.000 claims description 8
- 230000005856 abnormality Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 238000012732 spatial analysis Methods 0.000 claims description 2
- 230000002265 prevention Effects 0.000 abstract description 10
- 238000011160 research Methods 0.000 abstract description 4
- 238000012502 risk assessment Methods 0.000 abstract description 4
- 230000002123 temporal effect Effects 0.000 abstract description 4
- 241000196324 Embryophyta Species 0.000 description 6
- 235000016383 Zea mays subsp huehuetenangensis Nutrition 0.000 description 6
- 235000009973 maize Nutrition 0.000 description 6
- 230000006378 damage Effects 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 241000209094 Oryza Species 0.000 description 4
- 235000007164 Oryza sativa Nutrition 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 235000009566 rice Nutrition 0.000 description 4
- 241000482268 Zea mays subsp. mays Species 0.000 description 3
- 239000000843 powder Substances 0.000 description 3
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 210000005069 ears Anatomy 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000003337 fertilizer Substances 0.000 description 2
- 238000003306 harvesting Methods 0.000 description 2
- 230000003340 mental effect Effects 0.000 description 2
- 230000029553 photosynthesis Effects 0.000 description 2
- 238000010672 photosynthesis Methods 0.000 description 2
- 230000010152 pollination Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000013517 stratification Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 108091005804 Peptidases Proteins 0.000 description 1
- 239000004365 Protease Substances 0.000 description 1
- 102100037486 Reverse transcriptase/ribonuclease H Human genes 0.000 description 1
- XSQUKJJJFZCRTK-UHFFFAOYSA-N Urea Chemical compound NC(N)=O XSQUKJJJFZCRTK-UHFFFAOYSA-N 0.000 description 1
- 241000607479 Yersinia pestis Species 0.000 description 1
- 239000000443 aerosol Substances 0.000 description 1
- 238000005842 biochemical reaction Methods 0.000 description 1
- 238000009395 breeding Methods 0.000 description 1
- 230000001488 breeding effect Effects 0.000 description 1
- 239000004202 carbamide Substances 0.000 description 1
- 235000013339 cereals Nutrition 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 210000003763 chloroplast Anatomy 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000004720 fertilization Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000036512 infertility Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000021121 meiosis Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910000402 monopotassium phosphate Inorganic materials 0.000 description 1
- 235000019796 monopotassium phosphate Nutrition 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- PJNZPQUBCPKICU-UHFFFAOYSA-N phosphoric acid;potassium Chemical compound [K].OP(O)(O)=O PJNZPQUBCPKICU-UHFFFAOYSA-N 0.000 description 1
- 230000000243 photosynthetic effect Effects 0.000 description 1
- 230000003449 preventive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000017260 vegetative to reproductive phase transition of meristem 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
- 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/0635—Risk analysis of enterprise or organisation activities
-
- 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)
- Strategic Management (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Educational Administration (AREA)
- Animal Husbandry (AREA)
- General Health & Medical Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Marine Sciences & Fisheries (AREA)
- Agronomy & Crop Science (AREA)
- Development Economics (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
- Radiation Pyrometers (AREA)
Abstract
The present invention provides a kind of appraisal procedure of corn Climatic regionalization risk, and this method includes：MODIS LST remote sensing image datas are pre-processed, obtain surface temperature；According to surface temperature, by moving window algorithm, high temperature abnormal index is calculated, and temperature anomaly region is obtained according to high temperature abnormal index；According to surface temperature and corn florescence distribution situation, obtain in temperature anomaly region, the heat evil risk number of days and heat evil risk accumulated temperature in stage at corn florescence；And calculate accumulated temperature index.The present invention is using MODIS LST remote sensing image datas as foundation, and mutually the point based data of the meteorological site used in research, remote sensing image data can relatively accurately reflect the spatial and temporal distributions situation of high temperature, support provided for agriculture high temperature risk assessment than before；The Climatic regionalization risk class in the assessment stage at corn florescence of efficiently and accurately of the present invention, support, the efficiency and accuracy of the Climatic regionalization prevention of lifting corn are provided for the prevention of corn Climatic regionalization.
Description
Technical field
The present invention relates to agrometeorological hazard early warning field, more particularly to a kind of assessment side of corn Climatic regionalization risk
Method.
Background technology
During Maize Production, temperature is important meteorological factor, and temperature is too high to be caused to growing for corn
Have a strong impact on.As weather constantly warms in recent years, the frequency of occurrences of China corn producing region abnormal high temperature weather phenomenon is increasingly
Height, serious harm is brought to Maize Production.
High temperature influences the photosynthesis of corn.Under the high temperature conditions, the activity reduction of photosynthetic protease, chloroplast structure meet with
To destruction, stomata is caused to be closed, so that photosynthesis weakens；On the other hand, respiration strengthens under the high temperature conditions, consumption
Increase, dry-matter accumulation declines.The high temperature stress time is longer, and plant is aggrieved more serious, is more difficult to recover.High temperature forces corn to be given birth to
Educate various biochemical reactions in process to accelerate, each growing stage shortens.Such as the shortening of female fringe divergaence time, female fringe floret differentiation
Quantity is reduced, and fruit ear diminishes.Plant is become feeble and die too early in late growth stage high temperature, or terminate growing process in advance and enter into
Ripe phase, grouting time shorten, and dry-matter accumulation amount is reduced, and mass of 1000 kernel, unit weight, yield and quality reduce.
Maize male ears meiosis and loose powder phase are very sensitive to temperature.If running into dry weather in corn flowering period,
Abortive pollen is frequently formed, or makes pollen rapid dead and causes to be unable to normal fertilization solid.It is developed in maize male ears and blooms
During loose powder, high temperature is irrecoverable to the partial sterility process of some corn varieties, causes the setting percentage of corn significantly
Decline, tremendous influence is caused to production, serious meeting causes to have no harvest.Research shows that the high temperature stress duration is longer, plant by
Evil risk is more serious, more difficult to recover.
High temperature has highly important influence on corn growth, and corn growing season is different, and the critical value of high temperature is different therewith.It is beautiful
It is some stages in rice breeding time, such as corn florescence, sensitive to heat evil, if beautiful by Climatic regionalization, influence that can be serious
Rice grows, and corn sky bar, bald point occurs, lacks grain, lacks phenomena such as row, causes the decline of the yield and quality of corn, seriously
Influence peasant harvest yield and income.
Therefore, a kind of appraisal procedure of simple efficient corn high temperature risk how is proposed, turns into urgent problem to be solved.
The content of the invention
The present invention is the drawbacks described above for solving prior art, there is provided a kind of appraisal procedure of corn Climatic regionalization risk.
The present invention provides a kind of appraisal procedure of corn Climatic regionalization risk, including：
Step 1, according to the surface temperature for MODIS LST remote sensing image datas pre-process acquisition, Moving Window is passed through
Mental arithmetic method, high temperature abnormal index is calculated, and temperature anomaly region is obtained according to high temperature abnormal index；
Step 2, according to the surface temperature and corn florescence distribution situation, obtain in the temperature anomaly region, corn
The heat evil risk number of days and heat evil risk accumulated temperature in stage at florescence；According to the heat evil risk number of days and heat evil risk accumulated temperature, calculate
Accumulated temperature index.
Step 3, with reference to default corn Climatic regionalization grade scale, obtained according to the heat evil risk number of days and accumulated temperature index
In temperature anomaly region, the Climatic regionalization risk class in stage at corn florescence.
Wherein, in the step 1, the preprocessing process of the MODIS LST remote sensing image datas includes：
It is WGS-84 coordinate systems by remote sensing image data Coordinate Conversion, and pixel radiation brightness is converted into pixel earth's surface temperature
Degree.
Wherein, in the step 1, the high temperature abnormal index is：
The difference of MODIS LST remote sensing image datas center pel surface temperature values and pixel surface temperature average.
Wherein, in the step 1, described by moving window algorithm, calculating high temperature abnormal index includes：
Multiple window sizes are set, count the pixel surface temperature average value under different windows yardstick, and obtain in window
Imago member surface temperature value；
Calculate the difference of the window center pixel surface temperature value and the pixel surface temperature average value, the difference
For high temperature abnormal index.
Wherein, it is described to be included according to high temperature abnormal index acquisition temperature anomaly region in the step 1：
It is temperature anomaly region to choose stage at corn florescence high temperature abnormal index to continue larger pixel corresponding region.
Wherein, in the step 2, the heat evil risk number of days for obtaining the stage at corn florescence and heat evil risk accumulated temperature bag
Include：
Obtain in temperature anomaly region, in the stage at corn florescence, by the number of days of surface temperature >=default high temperature threshold value, definition
For heat evil risk number of days；
The heat evil risk accumulated temperature is the temperature build value in heat evil risk number of days, and calculation formula is：
In formula, Ts is the value of heat evil risk accumulated temperature；N is heat evil risk number of days；tiFor the heat evil risk temperature value of i-th day.
Wherein, the accumulated temperature index is the product of heat evil risk number of days and the harmful risk accumulated temperature of heat, and calculation formula is：
INt=Ts·n
In formula, INtFor the value of accumulated temperature index；Ts is the value of heat evil risk accumulated temperature；N is heat evil risk number of days.
Wherein, the default high temperature threshold value is 32~35 DEG C.
Wherein, in the step 2, the corn florescence distribution situation is entered by using the space interpolation instrument of spatial analysis
Row interpolation processing obtains.
Wherein, in the step 3, the default corn Climatic regionalization grade scale includes：
Using the accumulated temperature index in stage at corn florescence and heat evil risk number of days as the evaluation index of corn Climatic regionalization risk, lead to
Experiment and literature reading are crossed, corn Climatic regionalization risk class is divided into one-level~level Four by the order of severity from high to low.
The appraisal procedure of corn Climatic regionalization risk provided by the invention, based on MODIS LST remote sensing image datas, use
Moving window calculation window center pel and the difference of surrounding pixel surface temperature average, are defined as high temperature abnormal index.According to
The temperature anomaly region of high temperature abnormal index extraction Huang-Huai-Hai summer corn vitellarium, with the accumulated temperature index and heat in stage at corn florescence
Evil risk number of days is the evaluation index of corn Climatic regionalization risk, and Climatic regionalization risk class is divided into level Four, it is different to obtain temperature
In normal region, the Climatic regionalization risk class in stage at corn florescence.The present invention using MODIS LST remote sensing image datas as foundation,
Mutually than before research used in meteorological site point based data, remote sensing image data it is a technical advantage that can be more accurate
Ground reflects the spatial and temporal distributions situation of high temperature, accurately finds the high temperature abnormal area occurred in a wide range of, for the high warm air of agricultural
Danger, which is assessed, provides support；Climatic regionalization risk in stage at corn florescence etc. in the assessment temperature anomaly region of efficiently and accurately of the present invention
Level, support, the efficiency and accuracy of the Climatic regionalization prevention of lifting corn are provided for the prevention of corn Climatic regionalization.Avoid High Temperature Disaster
Caused by loss.
Brief description of the drawings
Fig. 1 is the appraisal procedure schematic flow sheet according to corn Climatic regionalization risk provided in an embodiment of the present invention；
Fig. 2 is according to MODIS LST invertings surface temperature data provided in an embodiment of the present invention；
Fig. 3 is according to high temperature abnormal index T-Tave under different scale in 1 day July in 2012 provided in an embodiment of the present invention
Distribution map；
Fig. 4 is according to Huang-Huai-Hai summer corn tasseling stage distribution map provided in an embodiment of the present invention.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
In accompanying drawing, the technical scheme in the embodiment of the present invention is explicitly described, it is clear that described embodiment be the present invention
Part of the embodiment, rather than whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art are not having
The every other embodiment obtained under the premise of creative work is made, belongs to the scope of protection of the invention.
Fig. 1 is according to the appraisal procedure schematic flow sheet of corn Climatic regionalization risk provided in an embodiment of the present invention, such as Fig. 1
Shown, this method includes:
Step S1, according to the surface temperature for MODIS LST remote sensing image datas pre-process acquisition, pass through Moving Window
Mental arithmetic method, high temperature abnormal index is calculated, and temperature anomaly region is obtained according to high temperature abnormal index；Step S2, according to described
Table temperature and corn florescence distribution situation, obtain in the temperature anomaly region, the stage at corn florescence heat evil risk number of days and
Heat evil risk accumulated temperature；According to the heat evil risk number of days and heat evil risk accumulated temperature, accumulated temperature index is calculated.Step S3, with reference to default
Corn Climatic regionalization grade scale, obtained according to the heat evil risk number of days and accumulated temperature index in temperature anomaly region, popcorn
The Climatic regionalization risk class in stage phase.
The appraisal procedure of corn Climatic regionalization risk provided in an embodiment of the present invention, based on MODISLST remote sensing image numbers
According to using the difference of moving window calculation window center pel and surrounding pixel surface temperature average, being defined as high temperature abnormality and refer to
Number.The temperature anomaly region of Huang-Huai-Hai summer corn vitellarium is extracted according to high temperature abnormal index, with the accumulated temperature in stage at corn florescence
Index and the evaluation index that heat evil risk number of days is corn Climatic regionalization risk, are divided into level Four by Climatic regionalization risk class, obtain
Take in temperature anomaly region, the Climatic regionalization risk class in stage at corn florescence.The embodiment of the present invention is with MODIS LST remote sensing shadows
As data be foundation, mutually study than before used in meteorological site point based data, the technical advantage of remote sensing image data
It is the spatial and temporal distributions situation that can relatively accurately reflect high temperature, accurately finds the high temperature exceptions area occurred in a wide range of
Domain, support is provided for agriculture high temperature risk assessment；Popcorn in the assessment temperature anomaly region of efficiently and accurately of the embodiment of the present invention
The Climatic regionalization risk class in stage phase, support, the effect of lifting corn Climatic regionalization prevention are provided for the prevention of corn Climatic regionalization
Rate and accuracy.Loss caused by avoiding High Temperature Disaster.
Wherein, in step S1, according to the surface temperature for MODIS LST remote sensing image datas pre-process acquisition, lead to
Moving window algorithm is crossed, calculates high temperature abnormal index, and temperature anomaly region is obtained according to high temperature abnormal index.
MODIS (moderate-resolution imaging spectroradiometer) full name be intermediate-resolution into
As spectrometer.LST (land surface temperature) surface temperature, the product serial number in NASA MODIS products
It is MOD11A2.
MODIS is an important sensor being mounted on terra and aqua satellites, will be uniquely seen in real time on satellite
Survey data directly to broadcast to the whole world by x wave bands, and can freely receive the spaceborne instrument of data and use without compensation, the whole world is permitted
More countries and regions are all receiving and used MODIS data.MODIS multi-wavelength data can provide reflection top simultaneously
Steam, aerosol, earth's surface temperature in situation, cloud border, cloud characteristic, ocean color, phytoplankton, biogeography, chemistry, air
The information of the features such as degree, cloud-top temperature, atmospheric temperature, ozone and cloud-top height.Multi-wavelength data can provide reaction land simultaneously
Steam, surface temperature, cloud top temperature in ground, cloud border, cloud characteristic, ocean color, phytoplankton, biogeography, chemistry, air
The information of the features such as degree, atmospheric temperature, ozone and cloud-top height, for land table, biosphere, Solid Earth, air and ocean
Carry out long-term global observation.
The preprocessing process of MODIS LST remote sensing image data includes, and MODIS LST remote sensing image datas coordinate is turned
WGS-84 coordinate systems are changed to, and pixel radiation brightness is converted into pixel surface temperature.Surface temperature distribution map in acquisition.WGS-
84 coordinate systems are a kind of geocentric coordinate systems used in the world.The origin of coordinates is earth centroid, its earth's core rectangular coordinate system in space
Z axis point to agreement earth pole (CTP) directions for defining of BIH (international time service organization) 1984.O, X-axis points to BIH
1984.0 zero meridian plane and the intersection point in CTP equator, Y-axis are vertically formed right-handed coordinate system, the referred to as world in 1984 with Z axis, X-axis
Geodetic coordinate system.
In the present embodiment, on July 1st, 2012 is obtained to the MODIS LST of the Yellow River and Huai He River sea region of August 31 remote sensing image
Data, and pre-processed, it is anti-to obtain the Yellow River and Huai He River sea region MODIS LST of on July 1st, 2012 to each day among August 31
Surface temperature distribution map is drilled, as shown in Figure 2.Fig. 2 is according to MODIS LST invertings surface temperature provided in an embodiment of the present invention point
Butut.
Further, the surface temperature obtained after being pre-processed according to MODIS LST remotely-sensed datas.Done a sum orally using Moving Window
Method, multiple window sizes are set, count the pixel surface temperature average value under different windows yardstick, and obtain window center pixel
Surface temperature value.Center pel surface temperature value and the difference of the pixel surface temperature average value are calculated, the difference is height
Temperature abnormality index (T-Tave).
Fig. 3 is according to high temperature abnormal index T-Tave under different scale in 1 day July in 2012 provided in an embodiment of the present invention
Distribution map, wherein, Fig. 3 (a) is high temperature abnormal index T-Tave distribution maps when yardstick is 3 × 3；During Fig. 3 (b) yardstick be 10 ×
High temperature abnormal index T-Tave distribution maps when 10；Fig. 3 (c) is high temperature abnormal index T-Tave distributions when yardstick is 50 × 50
Figure.It is temperature anomaly region to choose high temperature abnormal index to continue larger value of pixel corresponding region, and temperature anomaly region is entered
Row further evaluation.
Wherein, pixel is the minimum unit for forming digitized image.In remote sensing data acquiring, during such as scanning imagery, it is
Sensor is scanned the minimum unit of sampling to ground scenery；In Digital Image Processing, it is that analog image is swept
Retouch sampled point during digitlization.When resolution ratio is 1 km, a pixel represents the area of the km of the km of ground 1 × 1, i.e., 1 is flat
Square km；When resolution ratio is 30 meters, a pixel represents the area on 30 meters × 30 meters of ground；When resolution ratio is 1 meter, also
It is to say, area of on image the pixel equivalent to 1 meter of ground, 1 meter of x, i.e., 1 square metre.
Wherein, in step S2, according to the surface temperature and corn florescence distribution situation, the heat in acquisition stage at corn florescence
Evil risk number of days and heat evil risk accumulated temperature；According to the heat evil risk number of days and heat evil risk accumulated temperature, accumulated temperature index is calculated.
The present embodiment chooses Huang-Huai-Hai conduct, according to maize growth period data in Huang-Huai-Hai in table 1, uses sky
Between the space interpolation instrument analyzed carry out interpolation processing, the distribution map of corn tasseling stage is obtained, as shown in figure 4, Fig. 4 is according to this
The Huang-Huai-Hai summer corn tasseling stage distribution map that inventive embodiments provide.
The Huang-Huai-Hai maize growth period data few examples of table 1
Further, the heat evil risk number of days and heat evil risk accumulated temperature in stage at corn florescence are obtained, according to July 1 in 2012
Arrive the Yellow River and Huai He River sea region MODIS LST surface temperature data of each day among August 31, the heat evil in statistics stage at corn florescence day
Risk number of days, calculate heat evil risk accumulated temperature.And according to the heat evil risk number of days and heat evil risk accumulated temperature, calculate temperature anomaly area
The accumulated temperature index in domain.
Wherein, in step S3, with reference to default corn Climatic regionalization grade scale, according to the heat evil risk number of days and accumulated temperature
Index is obtained in temperature anomaly region, the Climatic regionalization risk class in stage at corn florescence.
Specifically, using assessment of the accumulated temperature index in stage at corn florescence with heat evil risk number of days as corn Climatic regionalization risk
Index, by experiment and literature reading, default Climatic regionalization grade scale is obtained, table 2 is high temperature provided in an embodiment of the present invention
Heat evil risk stratification standard.Corn Climatic regionalization risk class is divided into one to level Four by the order of severity from high to low, respectively with it is red,
Orange, yellow, Lan Si kinds color represent.
The accumulated temperature index of table 2 and high-temperature duration are the same as the Climatic regionalization risk class table of comparisons
Further, the accumulated temperature index and heat evil risk number of days in temperature anomaly region obtained according to calculating, and in table 2
Climatic regionalization risk stratification matches criteria.Obtain in temperature anomaly region, the Climatic regionalization risk class in stage at corn florescence.
Corn high temperature methods of risk assessment provided in an embodiment of the present invention, done harm to the accumulated temperature index and heat in stage at corn florescence
Risk number of days is the evaluation index of corn Climatic regionalization risk, corn Climatic regionalization risk class is divided into level Four, and obtain temperature
Spend in abnormal area, the Climatic regionalization risk class in stage at corn florescence.Support is provided for the prevention of corn Climatic regionalization, lifting is beautiful
The efficiency and accuracy of rice Climatic regionalization prevention.Loss caused by avoiding High Temperature Disaster.
On the basis of above-described embodiment, obtaining temperature anomaly region according to high temperature abnormal index includes：
It is temperature anomaly region to choose stage at corn florescence high temperature abnormal index to continue larger pixel corresponding region.
Specifically, the MODIS according to the 1 day July in 2012 that above-described embodiment obtains to the Yellow River and Huai He River sea region of August 31
LST remote sensing image data, stage at corn florescence daily high temperature abnormal index value being calculated, screening high temperature abnormal index is larger,
And occurring the more pixel of the larger frequency of high temperature abnormal index in the stage at corn florescence, the pixel corresponding region is temperature anomaly
Region.
On the basis of the various embodiments described above, according to the surface temperature and corn florescence distribution situation, popcorn is obtained
The heat evil risk number of days and heat evil risk accumulated temperature in stage phase include：
Obtain in temperature anomaly region, in the stage at corn florescence, the number of days of surface temperature >=default high temperature threshold value, be defined as
Heat evil risk number of days；
The heat evil risk accumulated temperature is the temperature build value in heat evil risk number of days, and calculation formula is：
In formula, Ts is the value of heat evil risk accumulated temperature；N is heat evil risk number of days；tiFor the heat evil risk temperature value of i-th day.
Wherein, it is 32~35 DEG C to preset high temperature threshold value, it is preferred that in the present embodiment, it is 34 DEG C to preset high temperature threshold value.
The accumulated temperature index is the product of heat evil risk number of days and the harmful risk accumulated temperature of heat, and calculation formula is：
INt=Ts·n
In formula, INtFor the value of accumulated temperature index；Ts is the value of heat evil risk accumulated temperature；N is heat evil risk number of days.
According to surface temperature and corn florescence distribution situation, obtain in temperature anomaly region, in the heat in stage at corn florescence
Evil risk number of days and accumulated temperature index, using the accumulated temperature index in stage at corn florescence and heat evil risk number of days as corn Climatic regionalization risk
Evaluation index, Climatic regionalization risk class is divided into level Four, and obtain in temperature anomaly region, the high temperature in stage at corn florescence
Heat evil risk class.
Fig. 4 is according to Huang-Huai-Hai summer corn tasseling stage distribution map provided in an embodiment of the present invention.Obtain temperature anomaly region
It is interior, can be targetedly for the higher area of Climatic regionalization risk class after the Climatic regionalization risk class in stage at corn florescence
Domain takes preventive measures.It it is the period most sensitive to high temperature when summer corn is in tasseling stage.It is beautiful to this part summer to mitigate high temperature
The harm of rice, it is conditional to take the measures such as cooling, artificial supplementary pollination, foliage application of pouring water.Field can be improved by pouring water
Miniclimate, temperature l~2 DEG C between strain are reduced, increase relative humidity, effectively mitigate high temperature damage.Artificial supplementary pollination, make to fall in post
Pollen amount increase on head, improves setting percentage.Repeatedly sprayed with urea, potassium dihydrogen phosphate or other high-quality foliar fertilizers etc.,
Increase plant fringe portion moisture, can decreasing temperature and increasing humidity, promote loose powder, while moisture and nutrient can be provided to blade, notice that solution is dense
Degree can not be too big.Impose and attack granulated fertilizer and can improve mass of 1000 kernel；Serious aggrieved plant is cut off, improves field microclimate, mitigates disease pest
Evil occurs, and increases individual plant production capacity by every possible means, mitigates production loss.
The appraisal procedure of corn Climatic regionalization risk provided in an embodiment of the present invention, based on MODISLST remote sensing image numbers
According to using the difference of moving window calculation window center pel and surrounding pixel surface temperature average, being defined as high temperature abnormality and refer to
Number.The temperature anomaly region of Huang-Huai-Hai summer corn vitellarium is extracted according to high temperature abnormal index, with the accumulated temperature in stage at corn florescence
Index and the evaluation index that heat evil risk number of days is corn Climatic regionalization risk, are divided into level Four by Climatic regionalization risk class, obtain
Take in temperature anomaly region, the Climatic regionalization risk class in stage at corn florescence.The present invention is with MODIS LST remote sensing image datas
For foundation, the mutually point based data of the meteorological site used in research than before, remote sensing image data it is a technical advantage that energy
Enough spatial and temporal distributions situations for relatively accurately reflecting high temperature, the high temperature abnormal area occurred in a wide range of is accurately found, is agriculture
The risk assessment of industry high temperature provides support；The high warm in stage at corn florescence in the assessment temperature anomaly region of efficiently and accurately of the present invention
Evil risk class, support, the efficiency and accuracy of the Climatic regionalization prevention of lifting corn are provided for the prevention of corn Climatic regionalization.Avoid
Loss caused by High Temperature Disaster.
Finally it should be noted that：The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although
The present invention is described in detail with reference to the foregoing embodiments, it will be understood by those within the art that：It still can
It is enough that technical scheme described in foregoing individual embodiments is modified, or which part technical characteristic is equally replaced
Change；And these modifications or replacement, the essence of appropriate technical solution is departed from each embodiment technical scheme of the present invention
Spirit and scope.
Claims (10)
- A kind of 1. appraisal procedure of corn Climatic regionalization risk, it is characterised in that including：Step 1, according to the surface temperature for MODIS LST remote sensing image datas pre-process acquisition, done a sum orally by Moving Window Method, high temperature abnormal index is calculated, and temperature anomaly region is obtained according to high temperature abnormal index；Step 2, according to the surface temperature and corn florescence distribution situation, obtain in the temperature anomaly region, the corn florescence The heat evil risk number of days and heat evil risk accumulated temperature in stage；According to the heat evil risk number of days and heat evil risk accumulated temperature, accumulated temperature is calculated Index；Step 3, with reference to default corn Climatic regionalization grade scale, temperature is obtained according to the heat evil risk number of days and accumulated temperature index In abnormal area, the Climatic regionalization risk class in stage at corn florescence.
- 2. the appraisal procedure of corn Climatic regionalization risk according to claim 1, it is characterised in that in the step 1, institute Stating the preprocessing process of MODIS LST remote sensing image datas includes：It is WGS-84 coordinate systems by remote sensing image data Coordinate Conversion, and pixel radiation brightness is converted into pixel surface temperature.
- 3. the appraisal procedure of corn Climatic regionalization risk according to claim 1, it is characterised in that in the step 1, institute Stating high temperature abnormal index is：The difference of MODIS LST remote sensing image datas center pel surface temperature values and pixel surface temperature average.
- 4. the appraisal procedure of corn Climatic regionalization risk according to claim 3, it is characterised in that in the step 1, institute State by moving window algorithm, calculating high temperature abnormal index includes：Multiple window sizes are set, count the pixel surface temperature average value under different windows yardstick, and obtain window center picture First surface temperature value；The difference of the window center pixel surface temperature value and the pixel surface temperature average value is calculated, the difference is height Temperature abnormality index.
- 5. the appraisal procedure of the corn Climatic regionalization risk according to claim 3 or 4, it is characterised in that the step 1 In, it is described to be included according to high temperature abnormal index acquisition temperature anomaly region：It is temperature anomaly region to choose stage at corn florescence high temperature abnormal index to continue larger pixel corresponding region.
- 6. the appraisal procedure of corn Climatic regionalization risk according to claim 1, it is characterised in that in the step 2, institute Stating the heat evil risk number of days for obtaining the stage at corn florescence and heat evil risk accumulated temperature includes：Obtain in temperature anomaly region, in the stage at corn florescence, by the number of days of surface temperature >=default high temperature threshold value, be defined as heat Evil risk number of days；The heat evil risk accumulated temperature is the temperature build value in heat evil risk number of days, and calculation formula is：<mrow> <msub> <mi>T</mi> <mi>s</mi> </msub> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>t</mi> <mi>i</mi> </msub> </mrow>In formula, Ts is the value of heat evil risk accumulated temperature；N is heat evil risk number of days；tiFor the heat evil risk temperature value of i-th day.
- 7. the appraisal procedure of corn Climatic regionalization risk according to claim 6, it is characterised in that the accumulated temperature index is The product of heat evil risk number of days and heat evil risk accumulated temperature, calculation formula are：INt=Ts·nIn formula, INtFor the value of accumulated temperature index；Ts is the value of heat evil risk accumulated temperature；N is heat evil risk number of days.
- 8. the appraisal procedure of corn Climatic regionalization risk according to claim 6, it is characterised in that the default high temperature threshold It is worth for 32~35 DEG C.
- 9. the appraisal procedure of corn Climatic regionalization risk according to claim 1, it is characterised in that in the step 2, institute State corn florescence distribution situation and interpolation processing acquisition is carried out by using the space interpolation instrument of spatial analysis.
- 10. the appraisal procedure of corn Climatic regionalization risk according to claim 1, it is characterised in that in the step 3, The default corn Climatic regionalization grade scale includes：Using the accumulated temperature index in stage at corn florescence and heat evil risk number of days as the evaluation index of corn Climatic regionalization risk, pass through reality Test and literature reading, corn Climatic regionalization risk class is divided into one-level~level Four by the order of severity from high to low.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710712114.2A CN107392503B (en) | 2017-08-18 | 2017-08-18 | Method for evaluating high-temperature heat damage risk of corn |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710712114.2A CN107392503B (en) | 2017-08-18 | 2017-08-18 | Method for evaluating high-temperature heat damage risk of corn |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107392503A true CN107392503A (en) | 2017-11-24 |
| CN107392503B CN107392503B (en) | 2020-12-25 |
Family
ID=60352955
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710712114.2A Active CN107392503B (en) | 2017-08-18 | 2017-08-18 | Method for evaluating high-temperature heat damage risk of corn |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107392503B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108764542A (en) * | 2018-05-16 | 2018-11-06 | 贾翔 | The remote sensing prediction method of phase and occurrence quantity occur for diversiform-leaved poplar spring looper |
| CN108876084A (en) * | 2018-04-08 | 2018-11-23 | 中国农业大学 | A kind of assessment prediction method of corn high temperature risk |
| CN108961089A (en) * | 2018-05-28 | 2018-12-07 | 浙江大学 | Based on star the In Middle And Lower Reaches of Changjiang River rice Climatic regionalization remote-sensing monitoring method of multi-source data |
| CN110163459A (en) * | 2018-03-23 | 2019-08-23 | 河南工业大学 | A method of building multiple index evaluation model is classified wheat quality |
| CN110956322A (en) * | 2019-11-28 | 2020-04-03 | 河南省气象科学研究所 | Summer corn flowering phase high-temperature disaster risk prediction method under climate warming trend |
| CN114568239A (en) * | 2022-03-01 | 2022-06-03 | 北京飞花科技有限公司 | Cotton high-temperature heat damage prediction method |
| CN116701859A (en) * | 2023-05-29 | 2023-09-05 | 河北省科学院地理科学研究所 | Plant activity accumulated temperature estimation method based on full remote sensing data |
| CN117892165A (en) * | 2024-03-14 | 2024-04-16 | 江苏省气象台 | Low-temperature disaster agricultural influence prediction method based on disaster analysis |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102116861A (en) * | 2011-02-01 | 2011-07-06 | 环境保护部卫星环境应用中心 | Method for extracting straw burning fire based on No. 1 environment satellite |
| CN103760619A (en) * | 2014-01-07 | 2014-04-30 | 中国神华能源股份有限公司 | Method and device for monitoring coal field fire zone |
| CN105557411A (en) * | 2015-12-15 | 2016-05-11 | 广西壮族自治区农业科学院水稻研究所 | Method for identifying heat-resistant seed setting performance of early rice at earing and flowering stages based on hour cumulative temperature |
| CN105760978A (en) * | 2015-07-22 | 2016-07-13 | 北京师范大学 | Agricultural drought grade monitoring method based on temperature vegetation drought index (TVDI) |
| CN103954973B (en) * | 2014-04-12 | 2016-08-31 | 中国农业科学院农业资源与农业区划研究所 | One method from VIIRS data inversion surface temperature |
| CN106171707A (en) * | 2016-07-21 | 2016-12-07 | 中国水稻研究所 | A kind of method for early warning of Oryza sativa L. Climatic regionalization |
| CN106778516A (en) * | 2016-11-24 | 2017-05-31 | 中国科学院地理科学与资源研究所 | A kind of method and device that surface temperature is determined using Chinese No. three satellite remote sensing dates of wind and cloud |
| CN106952171A (en) * | 2017-03-01 | 2017-07-14 | 中国农业大学 | A kind of corn high temperature risk class computational methods based on temperature record |
-
2017
- 2017-08-18 CN CN201710712114.2A patent/CN107392503B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102116861A (en) * | 2011-02-01 | 2011-07-06 | 环境保护部卫星环境应用中心 | Method for extracting straw burning fire based on No. 1 environment satellite |
| CN103760619A (en) * | 2014-01-07 | 2014-04-30 | 中国神华能源股份有限公司 | Method and device for monitoring coal field fire zone |
| CN103954973B (en) * | 2014-04-12 | 2016-08-31 | 中国农业科学院农业资源与农业区划研究所 | One method from VIIRS data inversion surface temperature |
| CN105760978A (en) * | 2015-07-22 | 2016-07-13 | 北京师范大学 | Agricultural drought grade monitoring method based on temperature vegetation drought index (TVDI) |
| CN105557411A (en) * | 2015-12-15 | 2016-05-11 | 广西壮族自治区农业科学院水稻研究所 | Method for identifying heat-resistant seed setting performance of early rice at earing and flowering stages based on hour cumulative temperature |
| CN106171707A (en) * | 2016-07-21 | 2016-12-07 | 中国水稻研究所 | A kind of method for early warning of Oryza sativa L. Climatic regionalization |
| CN106778516A (en) * | 2016-11-24 | 2017-05-31 | 中国科学院地理科学与资源研究所 | A kind of method and device that surface temperature is determined using Chinese No. three satellite remote sensing dates of wind and cloud |
| CN106952171A (en) * | 2017-03-01 | 2017-07-14 | 中国农业大学 | A kind of corn high temperature risk class computational methods based on temperature record |
Non-Patent Citations (3)
| Title |
|---|
| 刘哲 等: "基于MODIS 数据的黄淮海夏玉米高温风险空间分布", 《农 业 工 程 学 报》 * |
| 李德 等: "淮北平原夏玉米花期高温热害综合气候指标研究", 《中国生态农业学报》 * |
| 王越 等: "基于多窗口线性回归模型的浙北地区冬季气温估算", 《地理与地理信息科学》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110163459B (en) * | 2018-03-23 | 2023-04-21 | 河南工业大学 | Multi-index evaluation method for constructing wheat quality grading |
| CN110163459A (en) * | 2018-03-23 | 2019-08-23 | 河南工业大学 | A method of building multiple index evaluation model is classified wheat quality |
| CN108876084A (en) * | 2018-04-08 | 2018-11-23 | 中国农业大学 | A kind of assessment prediction method of corn high temperature risk |
| CN108764542A (en) * | 2018-05-16 | 2018-11-06 | 贾翔 | The remote sensing prediction method of phase and occurrence quantity occur for diversiform-leaved poplar spring looper |
| CN108961089A (en) * | 2018-05-28 | 2018-12-07 | 浙江大学 | Based on star the In Middle And Lower Reaches of Changjiang River rice Climatic regionalization remote-sensing monitoring method of multi-source data |
| CN110956322A (en) * | 2019-11-28 | 2020-04-03 | 河南省气象科学研究所 | Summer corn flowering phase high-temperature disaster risk prediction method under climate warming trend |
| CN110956322B (en) * | 2019-11-28 | 2023-06-20 | 河南省气象科学研究所 | Summer maize flowering phase high-temperature disaster risk prediction method under climate warming trend |
| CN114568239A (en) * | 2022-03-01 | 2022-06-03 | 北京飞花科技有限公司 | Cotton high-temperature heat damage prediction method |
| CN114568239B (en) * | 2022-03-01 | 2023-08-15 | 北京飞花科技有限公司 | Cotton high-temperature heat damage prediction method |
| CN116701859A (en) * | 2023-05-29 | 2023-09-05 | 河北省科学院地理科学研究所 | Plant activity accumulated temperature estimation method based on full remote sensing data |
| CN116701859B (en) * | 2023-05-29 | 2024-01-30 | 河北省科学院地理科学研究所 | Plant activity accumulated temperature estimation method based on full remote sensing data |
| CN117892165A (en) * | 2024-03-14 | 2024-04-16 | 江苏省气象台 | Low-temperature disaster agricultural influence prediction method based on disaster analysis |
| CN117892165B (en) * | 2024-03-14 | 2024-05-24 | 江苏省气象台 | Low-temperature disaster agricultural influence prediction method based on disaster analysis |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107392503B (en) | 2020-12-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107392503A (en) | A kind of appraisal procedure of corn Climatic regionalization risk | |
| 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 | |
| Benincasa et al. | Reliability of NDVI derived by high resolution satellite and UAV compared to in-field methods for the evaluation of early crop N status and grain yield in wheat | |
| Hall et al. | Optical remote sensing applications in viticulture‐a review | |
| Son et al. | Classification of multitemporal Sentinel-2 data for field-level monitoring of rice cropping practices in Taiwan | |
| CN106372592B (en) | A kind of winter wheat planting area calculation method based on winter wheat area index | |
| Duan et al. | Remote estimation of grain yield based on UAV data in different rice cultivars under contrasting climatic zone | |
| CN108985260A (en) | A kind of remote sensing and meteorological integrated rice yield estimation method | |
| Meng et al. | Corn yield forecasting in northeast china using remotely sensed spectral indices and crop phenology metrics | |
| CN105372672B (en) | Southern winter kind crops planting area extracting method based on time series data | |
| WO2015139356A1 (en) | Method and system for determining spatial scale for field investigation and field management | |
| Park et al. | Estimation of crop water stress in a nectarine orchard using high-resolution imagery from unmanned aerial vehicle (UAV) | |
| CN110954650A (en) | Satellite remote sensing monitoring method for cotton canopy nitrogen | |
| CN115950838A (en) | Summer corn drought unmanned aerial vehicle rapid monitoring and distinguishing method based on chlorophyll content | |
| CN116502901A (en) | Farmland drought spatial distribution evaluation method based on optical remote sensing data | |
| Zhang et al. | Evaluating maize evapotranspiration using high-resolution UAV-based imagery and FAO-56 dual crop coefficient approach | |
| Gusso et al. | Monitoring heat waves and their impacts on summer crop development in southern Brazil | |
| CN107491724B (en) | A kind of Spike Differentiation in Winter Wheat phase recognition methods and device | |
| Darmawan et al. | Seasonal analysis of precipitation, drought and vegetation index in Indonesian paddy field based on remote sensing data | |
| Johnson et al. | Vineyard canopy density mapping with IKONOS satellite imagery | |
| CN116011812A (en) | Weather disaster early warning system and method based on agricultural Internet of things | |
| Panda et al. | Blueberry crop growth analysis using climatologic factors and multi-temporal remotely sensed imagery | |
| CN115512233A (en) | Multi-temporal remote sensing extraction method for planting attributes of land scale cultivated land | |
| Bouskour et al. | The use of NDVI to improve cereals agriculture: A review | |
| Tatsumi | Cropping intensity and seasonality parameters across Asia extracted by multi-temporal SPOT vegetation data |
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 |