[go: up one dir, main page]

CN113984212A - Agricultural irrigated area extraction method and system - Google Patents

Agricultural irrigated area extraction method and system Download PDF

Info

Publication number
CN113984212A
CN113984212A CN202111257130.XA CN202111257130A CN113984212A CN 113984212 A CN113984212 A CN 113984212A CN 202111257130 A CN202111257130 A CN 202111257130A CN 113984212 A CN113984212 A CN 113984212A
Authority
CN
China
Prior art keywords
area
remote sensing
sensing image
crop planting
planting area
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
Application number
CN202111257130.XA
Other languages
Chinese (zh)
Other versions
CN113984212B (en
Inventor
周广胜
王树东
汲玉河
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chinese Academy of Meteorological Sciences CAMS
Aerospace Information Research Institute of CAS
Original Assignee
Chinese Academy of Meteorological Sciences CAMS
Aerospace Information Research Institute of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chinese Academy of Meteorological Sciences CAMS, Aerospace Information Research Institute of CAS filed Critical Chinese Academy of Meteorological Sciences CAMS
Priority to CN202111257130.XA priority Critical patent/CN113984212B/en
Publication of CN113984212A publication Critical patent/CN113984212A/en
Application granted granted Critical
Publication of CN113984212B publication Critical patent/CN113984212B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J2005/0077—Imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N2021/1793—Remote sensing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N2021/1793—Remote sensing
    • G01N2021/1797—Remote sensing in landscape, e.g. crops

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Image Processing (AREA)

Abstract

本发明提供一种农业灌区提取方法及系统，在提取过程中，通过对历史作物种植区域在当前年份的热红外遥感影像数据进行降尺度处理，并根据精细化热红外温度数据以及植被覆盖度数据构建的特征空间，确定当前年份作物种植区域的温度植被干旱指数，可以达到精准提取农业灌区的目的。

Figure 202111257130

The invention provides an agricultural irrigation area extraction method and system. In the extraction process, the thermal infrared remote sensing image data of the historical crop planting area in the current year is downscaled, and the refined thermal infrared temperature data and vegetation coverage data are processed according to the downscaling process. The constructed feature space can determine the temperature and vegetation drought index of the crop planting area in the current year, which can achieve the purpose of accurately extracting agricultural irrigation areas.

Figure 202111257130

Description

Agricultural irrigated area extraction method and system
Technical Field
The invention relates to the technical field of information extraction, in particular to an agricultural irrigation area extraction method and system.
Background
The irrigation area refers to an area for artificial irrigation, and the area and distribution information of the agricultural irrigation area in arid and semi-arid areas are very important for estimating the growth vigor and the yield of crops. In addition, as water resources (including rainfall, groundwater, river and lake water resources and the like) in the arid and semi-arid regions are relatively short, accurate agricultural irrigation area distribution and area information are of great importance to scientific allocation and management of the water resources in the arid and semi-arid regions.
At present, the extraction method of agricultural irrigation area information in arid and semi-arid regions generally comprises a manual statistical method and a remote sensing method. The manual statistical method is mainly used for obtaining information of the agricultural irrigation area through manual field statistics and summary reporting, and the method has the problems that the method is not objective enough and space-time dynamic agricultural irrigation area distribution information is difficult to obtain. The remote sensing method is the most common method applied at home and abroad, has the advantages of objectivity, rapidness, low cost and the like, and comprises a land utilization analysis method, a growth monitoring method, a soil moisture analysis method and the like. The land utilization analysis method is to estimate the farmland information obtained by a remote sensing technology method and combine experience, and the method can effectively obtain the farmland distribution information, but can be identified by a large amount of experience. The growth monitoring method mainly comprises the steps of observing remote sensing information at different times, such as a canopy water analysis method, a vegetation index time sequence analysis method and the like. The canopy water analysis method judges whether the irrigation area is an agricultural irrigation area or not by quickly and accurately identifying canopy water information and comparing canopy information of other areas, but the method has higher requirements on observation time points and spectral resolution of remote sensing data; the vegetation index time series analysis method mainly analyzes the growth difference between an irrigated area and a non-irrigated area through the vegetation index of a time series and identifies the growth difference. The soil moisture analysis method mainly analyzes the soil moisture condition through microwave data or thermal infrared data, has a clear mechanism, is limited by spatial resolution, and is difficult to popularize and apply to broken or small areas of a land parcel.
Therefore, an agricultural irrigation area extraction method is urgently needed.
Disclosure of Invention
The invention provides an agricultural irrigated area extraction method and system, which are used for overcoming the defects in the prior art.
The invention provides an extraction method of an agricultural irrigated area, which comprises the following steps:
acquiring a historical crop planting area in a monitoring area of an agricultural planting area, and determining a crop planting area of a current year in the monitoring area of the agricultural planting area based on thermal infrared remote sensing image data of the historical crop planting area in the current year;
determining thermal infrared temperature data of a time sequence and vegetation coverage data in a monitoring time period based on the thermal infrared remote sensing image data, and performing downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data;
and determining a temperature vegetation drought index of the current year crop planting area based on the refined thermal infrared temperature data and the feature space constructed by the vegetation coverage data, and extracting the agricultural irrigation area in the current year crop planting area based on the temperature vegetation drought index.
According to the agricultural irrigated area extraction method provided by the invention, the determining of the temperature vegetation drought index of the crop planting area in the current year based on the feature space constructed by the refined thermal infrared temperature data and the vegetation coverage data specifically comprises the following steps:
and for any pixel of any scene current year thermal infrared remote sensing image in the current year thermal infrared remote sensing image data, determining a temperature vegetation drought index at any pixel based on the thermal infrared temperature data at any pixel, the vegetation coverage data and the wet edge temperature in the characteristic space.
According to the agricultural irrigated area extraction method provided by the invention, the extraction of the agricultural irrigated area in the crop planting area of the current year is carried out based on the temperature vegetation drought index, and the method specifically comprises the following steps:
for any pixel position in the current year thermal infrared remote sensing image data, determining an irrigated area identification index based on the temperature vegetation drought index of the pixel of each scene current year thermal infrared remote sensing image at the any pixel position in the current year thermal infrared remote sensing image data;
and judging whether any pixel position in the thermal infrared remote sensing image data of the current year is an irrigation area pixel position or not based on the irrigation area identification index, and extracting the agricultural irrigation area in the crop planting area of the current year based on a judgment result.
According to the agricultural irrigation area extraction method provided by the invention, the step of obtaining the historical crop planting area in the agricultural planting area monitoring area specifically comprises the following steps:
acquiring historical multi-year remote sensing image data close to the current year, and determining the remote sensing image data of each historical year in the historical multi-year remote sensing image data;
and determining a first-class normalized vegetation index at each pixel position in each scene remote sensing image in the remote sensing image data of each historical year, and determining the historical crop planting area based on the first-class normalized vegetation index at each pixel position in each scene remote sensing image in the remote sensing image data of each historical year.
According to the agricultural irrigated area extraction method provided by the invention, the historical crop planting area is determined based on the first normalized vegetation index of each pixel in each scene remote sensing image in the remote sensing image data of each historical year, and the method specifically comprises the following steps:
for remote sensing image data of any historical year, determining a first type normalized vegetation index mean value at each pixel position in the remote sensing image data of any historical year;
for any scene remote sensing image in the remote sensing image data of any historical year, determining a crop identification index at each pixel in any scene remote sensing image based on the first-class normalized vegetation index and the first-class normalized vegetation index mean value at each pixel in any scene remote sensing image;
judging whether each pixel in any scene remote sensing image is a crop planting pixel or not based on the crop identification index at each pixel in any scene remote sensing image, and determining a crop planting area corresponding to any scene remote sensing image based on a judgment result;
determining a crop planting area corresponding to the remote sensing image data of any historical year based on the crop planting area corresponding to each scene of remote sensing image in the remote sensing image data of any historical year;
and determining the historical crop planting area based on the crop planting area corresponding to the remote sensing image data of each historical year.
According to the agricultural irrigated area extraction method provided by the invention, the determination of the crop planting area of the current year in the agricultural planting area monitoring area based on the thermal infrared remote sensing image data of the historical crop planting area of the current year specifically comprises the following steps:
determining a second type of normalized vegetation index of each pixel position in the thermal infrared remote sensing image data at the starting point and the ending point of the monitoring time period;
and determining a crop planting area of the current year in the agricultural planting area monitoring area based on the second type normalized vegetation index of each pixel position in the thermal infrared remote sensing image data at the starting point and the ending point.
According to the agricultural irrigation area extraction method provided by the invention, the determining of the current year crop planting area in the agricultural planting area monitoring area based on the second normalized vegetation indexes of the positions of the pixels in the thermal infrared remote sensing image data at the starting point and the ending point specifically comprises the following steps:
for any pixel position in the thermal infrared remote sensing image data, determining a crop planting index of the current year at any pixel position based on the duration of the monitoring time period and the second normalized vegetation indexes of the pixel position at the starting point and the ending point;
and determining the current year crop planting area based on the current year crop planting index at each pixel position in the historical crop planting area.
According to the agricultural irrigated area extraction method provided by the invention, the extraction of the agricultural irrigated area in the crop planting area of the current year is carried out based on the temperature vegetation drought index, and then the method further comprises the following steps:
vectorizing the crop planting area in the current year to obtain a plurality of farmland patches, and determining agricultural irrigation areas corresponding to the farmland patches;
and for any farmland patch, if the relative area deviation of the farmland patch and the agricultural irrigation area corresponding to the farmland patch is less than or equal to a deviation threshold value, taking the farmland patch as the agricultural irrigation area.
According to the agricultural irrigated area extraction method provided by the invention, the step of obtaining the historical crop planting area in the agricultural planting area monitoring area also comprises the following steps:
acquiring rainfall information of a target area in the monitoring time period of years including the current year;
and determining the agricultural planting area monitoring area in the target area based on the precipitation information.
The invention also provides an agricultural irrigated area extraction system, which comprises:
the system comprises a first acquisition module, a second acquisition module and a third acquisition module, wherein the first acquisition module is used for acquiring a historical crop planting area in a monitoring area of an agricultural planting area and determining a current year crop planting area in the monitoring area of the agricultural planting area based on thermal infrared remote sensing image data of the historical crop planting area in the current year;
the second acquisition module is used for determining thermal infrared temperature data and vegetation coverage data of a time sequence in a monitoring time period based on the thermal infrared remote sensing image data, and performing downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data;
and the extraction module is used for determining the temperature vegetation drought index of the current year crop planting area based on the refined thermal infrared temperature data and the feature space constructed by the vegetation coverage data, and extracting the agricultural irrigation area in the current year crop planting area based on the temperature vegetation drought index.
The invention also provides electronic equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor executes the program to realize the steps of any one of the agricultural irrigation area extraction methods.
The present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of any of the agricultural irrigated area extraction methods described above.
The invention also provides a computer program product comprising a computer program which, when executed by a processor, carries out the steps of any of the methods of extraction of agricultural irrigated areas as described above.
The invention provides an extraction method and system for an agricultural irrigated area, which comprises the steps of firstly obtaining a historical crop planting area in a monitoring area of an agricultural planting area, and determining a current year crop planting area in the monitoring area of the agricultural planting area based on thermal infrared remote sensing image data of the historical crop planting area in the current year; then determining thermal infrared temperature data of a time sequence and vegetation coverage data in a monitoring time period based on the thermal infrared remote sensing image data, and carrying out downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data; and finally, determining the temperature vegetation drought index of the crop planting area in the current year based on the feature space constructed by the refined thermal infrared temperature data and the vegetation coverage data, and extracting the agricultural irrigation area in the crop planting area in the current year based on the temperature vegetation drought index. In the extraction process, the thermal infrared remote sensing image data of the historical crop planting area in the current year is subjected to downscaling processing, and the temperature vegetation drought index of the crop planting area in the current year is determined according to the feature space constructed by the refined thermal infrared temperature data and the vegetation coverage data, so that the purpose of accurately extracting the agricultural irrigation area can be achieved.
Drawings
In order to more clearly illustrate the technical solutions of the present invention or the prior art, the drawings needed for the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
FIG. 1 is a schematic flow chart of the agricultural irrigation area extraction method provided by the invention;
FIG. 2 is a schematic structural diagram of an agricultural irrigation area extraction system provided by the present invention;
fig. 3 is a schematic structural diagram of an electronic device provided in the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is obvious that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Fig. 1 is a schematic flow chart of an agricultural irrigation area extraction method provided in an embodiment of the present invention, and as shown in fig. 1, the method includes:
s1, acquiring a historical crop planting area in a monitoring area of an agricultural planting area, and determining a current year crop planting area in the monitoring area of the agricultural planting area based on thermal infrared remote sensing image data of the historical crop planting area in the current year;
s2, determining thermal infrared temperature data and vegetation coverage data of a time sequence in a monitoring time period based on the thermal infrared remote sensing image data, and carrying out downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data;
s3, determining the temperature vegetation drought index of the current year crop planting area based on the refined thermal infrared temperature data and the feature space constructed by the vegetation coverage data, and extracting the agricultural irrigation area in the current year crop planting area based on the temperature vegetation drought index.
Specifically, in the agricultural irrigation area extraction method provided in the embodiment of the present invention, an execution main body is a server, the server may be a local server or a cloud server, and the local server may specifically be a computer, a tablet computer, a smart phone, and the like, which is not specifically limited in the embodiment of the present invention.
Step S1 is executed first, and a historical crop planting area in a monitored area of an agricultural planting area is obtained, where the monitored area of the agricultural planting area is an area in the target area that needs to be monitored whether the area is an agricultural irrigation area. The historical crop planting area can be a union of the crop planting areas in the historical years, can be determined through land utilization information, and can also be obtained through extraction of historical multi-year remote sensing image data. The number of years involved in the history of years can be selected according to needs, for example, 10 years, 9 years and the like can be selected.
And then acquiring thermal infrared remote sensing image data of the historical crop planting area in the current year, wherein the thermal infrared remote sensing image data of the current year can comprise a plurality of thermal infrared remote sensing images, and each thermal infrared remote sensing image can comprise reflectivity data and thermal infrared temperature data of each pixel.
Furthermore, the crop planting area of the current year in the agricultural planting area monitoring area can be determined according to the thermal infrared remote sensing image data of the current year. According to the embodiment of the invention, the area change of the crop planting area in the agricultural planting area monitoring area can be determined according to the thermal infrared remote sensing image data, and then the crop planting area of the current year in the agricultural planting area monitoring area is determined. And determining the crop planting area of the current year from the historical crop planting area according to the thermal infrared remote sensing image data of the current year.
And then executing a step S2, and determining thermal infrared temperature data and vegetation coverage data of a time sequence in the monitoring time period according to the thermal infrared remote sensing image data of the current year. The monitoring period may be determined by:
acquiring a phenological curve of a monitoring area of an agricultural planting area; the phenological curve is determined according to crops planted in a monitoring area of the agricultural planting area, and different crops correspond to different phenological curves.
Determining the monitoring time period of the current year according to the phenological curve; the irrigation area monitoring time period of the current year may be an active time period of the crop, for example, the whole time period from the beginning of the crop seedling to the end of the crop before the crop is mature.
And determining thermal infrared temperature data and vegetation coverage data of the time sequence in the monitoring time period, namely obtaining the thermal infrared temperature data and vegetation coverage data of the time sequence at each pixel position in the thermal infrared remote sensing image data of the current year.
Because the obtained thermal infrared temperature data is generally low in resolution, in order to improve the resolution, downscaling processing, namely down-sampling processing, can be performed on the thermal infrared temperature data, and the processing mode can be realized by adopting the following formula:
LSThigh＝F(PRShigh)+LSTlow-F(PRSlow)
Figure BDA0003324448860000091
Figure BDA0003324448860000092
Figure BDA0003324448860000093
Figure BDA0003324448860000094
wherein, LSThighAnd LSTlowRespectively simulating a high-resolution earth surface temperature and a low-resolution earth surface temperature; LSTsdAnd LSTswRespectively the highest temperature of the bare soil surface when the vegetation coverage rate is 0 and the highest temperature of the vegetation surface when the vegetation is covered completely; k is a radical ofdryAnd kwetDetermining slopes of a trunk line and a wet line in a characteristic space for the simulated surface temperature and vegetation coverage data respectively; beta is the Bowen ratio; FVC is vegetation coverage, FVChighAnd FVClowRespectively simulating vegetation coverage with high resolution and vegetation coverage with low resolution; LST is observed thermal infrared data; epsilon is the error; f (PRS)high) And F (PRS)low) Simulation functions, PRS, of thermal infrared remote sensing images of high and low resolution, respectivelyhighAnd PRSlowRespectively simulated high-resolution and low-resolution thermal infrared remote sensing images, the PRS comprising PRShighAnd PRSlow。
Let LST in the above formulalowThe value of (A) is thermal infrared temperature data to obtain LSThighThe value is the refined thermal infrared temperature data.
And finally, executing a step S3, and constructing a feature space according to the obtained refined thermal infrared temperature data and the vegetation coverage data. The method comprises the steps of drawing a relation between refined thermal infrared temperature data and vegetation coverage data in a two-dimensional coordinate system with the refined thermal infrared temperature data as a vertical axis and the vegetation coverage data as a horizontal axis to obtain a feature space. The feature space includes a wet line and a main line, the wet line is parallel to the horizontal axis, the end point of the main line coincides with the end point of the wet line, and the slope of the main line is negative.
According to the constructed characteristic space, the temperature vegetation drought index of the crop planting area in the current year can be determined by combining thermal infrared temperature data. The temperature vegetation drought index can be used for representing whether vegetation drought exists in a crop planting area of the current year, and each pixel of each scene of the current year thermal infrared remote sensing image in the current year thermal infrared remote sensing image data corresponds to one temperature vegetation drought index. The temperature vegetation drought index at each pixel can be determined by thermal infrared temperature data and vegetation coverage data at each pixel, which are not specifically limited in the embodiments of the present invention.
According to the temperature vegetation drought index of the crop planting area in the current year, the agricultural irrigation area in the crop planting area in the current year can be extracted. When the agricultural irrigated area is extracted, the temperature vegetation drought index of each pixel of each scene of the current year thermal infrared remote sensing image in the current year thermal infrared remote sensing image data can be used for determining the irrigated area identification index of each pixel position in the current year thermal infrared remote sensing image data, whether each pixel position is the irrigated area pixel position can be determined according to the irrigated area identification index, and then the agricultural irrigated area in the current year crop planting area can be extracted.
The agricultural irrigated area extraction method provided by the embodiment of the invention comprises the steps of firstly, obtaining a historical crop planting area in a monitoring area of an agricultural planting area, and determining a current year crop planting area in the monitoring area of the agricultural planting area based on thermal infrared remote sensing image data of the historical crop planting area in the current year; then determining thermal infrared temperature data of a time sequence and vegetation coverage data in a monitoring time period based on the thermal infrared remote sensing image data, and carrying out downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data; and finally, determining the temperature vegetation drought index of the crop planting area in the current year based on the feature space constructed by the refined thermal infrared temperature data and the vegetation coverage data, and extracting the agricultural irrigation area in the crop planting area in the current year based on the temperature vegetation drought index. In the extraction process, the thermal infrared remote sensing image data of the historical crop planting area in the current year is subjected to downscaling processing, and the temperature vegetation drought index of the crop planting area in the current year is determined according to the feature space constructed by the refined thermal infrared temperature data and the vegetation coverage data, so that the purpose of accurately extracting the agricultural irrigation area can be achieved.
On the basis of the above embodiment, the method for extracting an agricultural irrigated area provided by the embodiment of the invention determines the temperature vegetation drought index of the crop planting area in the current year based on the feature space constructed by the refined thermal infrared temperature data and the vegetation coverage data, and specifically includes:
and for any pixel of any scene current year thermal infrared remote sensing image in the current year thermal infrared remote sensing image data, determining a temperature vegetation drought index at any pixel based on the thermal infrared temperature data at any pixel, the vegetation coverage data and the wet edge temperature in the characteristic space.
Specifically, in the embodiment of the present invention, the drought index of temperature vegetation at a certain pixel can be determined by the following formula:
Figure BDA0003324448860000111
in the formula, TDVI is a temperature vegetation drought index at a certain pixel, a and b are respectively a constant and a slope of a fitting equation of a main line in a characteristic space, LST is thermal infrared temperature data at the pixel, and LST isminThe method is characterized in that thermal infrared temperature data corresponding to a wet line in a characteristic space are provided, and the FVC is vegetation coverage data at the pixel.
In the embodiment of the invention, the accuracy of the temperature vegetation drought index can be improved by refining the characteristic space constructed by the thermal infrared temperature data and the vegetation coverage data.
On the basis of the above embodiment, the method for extracting an agricultural irrigated area provided by the embodiment of the invention, based on the temperature vegetation drought index, extracts the agricultural irrigated area in the crop planting area of the current year, and specifically includes:
for any pixel position in the current year thermal infrared remote sensing image data, determining an irrigated area identification index based on the temperature vegetation drought index of the pixel of each scene current year thermal infrared remote sensing image at the any pixel position in the current year thermal infrared remote sensing image data;
and judging whether any pixel position in the thermal infrared remote sensing image data of the current year is an irrigation area pixel position or not based on the irrigation area identification index, and extracting the agricultural irrigation area in the crop planting area of the current year based on a judgment result.
Specifically, in the embodiment of the present invention, when the agricultural irrigation area in the crop planting area of the current year is extracted through the temperature vegetation drought index, for any pixel position in the thermal infrared remote sensing image data of the current year, the irrigation area identification index may be determined according to the temperature vegetation drought index of the pixel of each scene of the thermal infrared remote sensing image of the current year at any pixel position in the thermal infrared remote sensing image data of the current year. For example, the irrigated area identification index may be determined by the following formula:
Figure BDA0003324448860000121
wherein, IAI is irrigation area identification index, TVDIiIs the temperature vegetation drought index of the pixel of the ith scene current year thermal infrared remote sensing image at any pixel position, n1 is the scene number in the current year thermal infrared remote sensing image data, C0Is a preset threshold.
Then, according to the irrigation area identification index, whether any pixel position in the thermal infrared remote sensing image data of the current year is the irrigation area pixel position is judged, that is, whether the irrigation area identification index is larger than or equal to a first threshold value can be judged, and the first threshold value can be set according to needs.
And finally, extracting the agricultural irrigation areas in the crop planting areas of the current year according to the judgment result. Namely, if the irrigation area identification index is larger than or equal to the first threshold value according to the judgment result, the position of any pixel can be determined to be the irrigation area pixel position, and the area formed by all irrigation area pixel positions in the thermal infrared remote sensing image data of the current year can be determined to be the agricultural irrigation area. Otherwise, if the judging result is that the irrigation area identification index is smaller than the first threshold, determining that any pixel position is not the irrigation area pixel position.
In the embodiment of the invention, the irrigated area identification index is introduced in the process of extracting the agricultural irrigated area in the crop planting area of the current year, so that the extraction result is more reliable.
On the basis of the above embodiment, the method for extracting an agricultural irrigation area provided in the embodiment of the present invention for obtaining a historical crop planting area in a monitored area of an agricultural planting area specifically includes:
acquiring historical multi-year remote sensing image data close to the current year, and determining the remote sensing image data of each historical year in the historical multi-year remote sensing image data;
and determining a first-class normalized vegetation index at each pixel position in each scene remote sensing image in the remote sensing image data of each historical year, and determining the historical crop planting area based on the first-class normalized vegetation index at each pixel position in each scene remote sensing image in the remote sensing image data of each historical year.
Specifically, in the embodiment of the invention, when the historical crop planting area in the agricultural planting area monitoring area is obtained, the historical multi-year remote sensing image data close to the current year, that is, the remote sensing image data of a plurality of historical years before the current year, can be obtained first. Further, the first type of normalized vegetation index at each pixel in each scene remote sensing image in the historical multi-year remote sensing image data can be determined. The remote sensing image data of each historical year can comprise n1 remote sensing images, and each remote sensing image can comprise n2 image elements. In each remote sensing image, the first normalized vegetation index at each pixel can be determined by the reflectivity data at the pixel, which is not specifically limited in the embodiment of the present invention.
Then, according to the first-class normalized vegetation index of each pixel in each scene remote sensing image in the remote sensing image data of each historical year, a crop planting area corresponding to the remote sensing image data of each historical year can be determined, and then based on the crop planting area corresponding to the remote sensing image data of each historical year, a historical crop planting area in the agricultural planting area monitoring area can be determined. The historical crop planting area is a crop planting area before the current year, and can be a union of crop planting areas corresponding to remote sensing image data of all historical years.
In the embodiment of the invention, the historical multi-year remote sensing image data close to the current year is adopted to obtain the historical crop planting area, so that the result is more accurate.
On the basis of the above embodiment, the method for extracting an agricultural irrigation area provided in the embodiment of the present invention is a method for determining a historical crop planting area based on a first normalized vegetation index at each pixel in each scene remote sensing image in the remote sensing image data of each historical year, and specifically includes:
for remote sensing image data of any historical year, determining a first type normalized vegetation index mean value at each pixel position in the remote sensing image data of any historical year;
for any scene remote sensing image in the remote sensing image data of any historical year, determining a crop identification index at each pixel in any scene remote sensing image based on the first-class normalized vegetation index and the first-class normalized vegetation index mean value at each pixel in any scene remote sensing image;
judging whether each pixel in any scene remote sensing image is a crop planting pixel or not based on the crop identification index at each pixel in any scene remote sensing image, and determining a crop planting area corresponding to any scene remote sensing image based on a judgment result;
determining a crop planting area corresponding to the remote sensing image data of any historical year based on the crop planting area corresponding to each scene of remote sensing image in the remote sensing image data of any historical year;
and determining the historical crop planting area based on the crop planting area corresponding to the remote sensing image data of each historical year.
Specifically, in the embodiment of the present invention, when determining the historical crop planting area, for the remote sensing image data of any historical year, the first normalized vegetation index average value at each pixel position in the remote sensing image data of any historical year is determined, that is, for any pixel position in the remote sensing image data of any historical year, the average value of the normalized vegetation index of the pixel at any pixel position of each scene remote sensing image in the remote sensing image data of any historical year can be represented as
Figure BDA0003324448860000141
For the jth scene remote sensing image in the remote sensing image data of the ith historical year, based on the first-class normalized vegetation index NDVI at the kth pixel position in the jth scene remote sensing imageiAnd normalized mean vegetation index of the first kind
Figure BDA0003324448860000142
Determining the crop at the kth pixel position in the jth scene remote sensing imageAn index is identified. Namely, the method comprises the following steps:
Figure BDA0003324448860000143
wherein the CRIkAnd identifying the index of the crop at the kth pixel in the jth scene remote sensing image.
And then, judging whether each pixel in the jth scene remote sensing image is a crop planting pixel or not according to the crop identification index at each pixel in the jth scene remote sensing image. Namely, CRI judgmentkWhether the threshold value is greater than or equal to the second threshold value may be set according to needs, which is not particularly limited in the embodiment of the present invention.
If CRIkIf the k-th pixel element is larger than or equal to the second threshold value, the k-th pixel element can be determined as a crop planting pixel element, otherwise, the k-th pixel element is represented.
And finally, according to the judgment result, determining a crop planting area corresponding to the jth scene remote sensing image, namely determining the number of crop planting pixels in the jth scene remote sensing image, and representing the crop planting area through the number.
According to the crop planting areas corresponding to the remote sensing images of the scenes in the remote sensing image data of the ith historical year, the crop planting areas corresponding to the remote sensing image data of the ith historical year can be determined, namely, a union of the crop planting areas corresponding to the remote sensing images of all scenes in the remote sensing image data of the ith historical year.
And determining the historical crop planting areas according to the union set of the crop planting areas corresponding to the remote sensing image data of each historical year. Namely, the method comprises the following steps:
ZWI＝A1∪A2∪...∪An3
wherein, ZWI is a historical crop planting area; a. the1、A2、…、An3The crop planting areas corresponding to the remote sensing image data of each historical year are respectively represented by the number of crop planting pixels, and n3 is the number of the historical years.
In the embodiment of the invention, when the historical crop planting area in the agricultural planting area monitoring area is determined, the historical crop planting area is determined through the crop planting area corresponding to the remote sensing image data of each historical year, so that the feasibility of the scheme can be ensured.
On the basis of the above embodiment, the method for extracting an agricultural irrigation area provided in the embodiment of the present invention, which determines a crop planting area of a current year in a monitoring area of an agricultural planting area based on thermal infrared remote sensing image data of the historical crop planting area of the current year, specifically includes:
determining a second type of normalized vegetation index of each pixel position in the thermal infrared remote sensing image data at the starting point and the ending point of the monitoring time period;
and determining a crop planting area of the current year in the agricultural planting area monitoring area based on the second type normalized vegetation index of each pixel position in the thermal infrared remote sensing image data at the starting point and the ending point.
Specifically, in the embodiment of the present invention, when determining the crop planting area of the current year, the thermal infrared remote sensing image data of the current year, that is, the thermal infrared remote sensing image data of the historical crop planting area in the current year, may be determined first. And then determining a second type normalized vegetation index of each pixel position in the thermal infrared remote sensing image data of the current year at the starting point and the ending point of the monitoring time period. Similarly, the thermal infrared remote sensing image data of the current year includes n1 scenes, and each scene may include n2 pixels.
And determining a current annual crop planting area in the monitoring area of the agricultural planting area, namely the current annual crop planting area in the historical crop planting area according to the second normalized vegetation index of each pixel position in the thermal infrared remote sensing image data at the starting point of the monitoring time period and the ending point of the monitoring time period. In the embodiment of the present invention, the crop planting area of the current year may be determined by the second-type normalized vegetation index at the starting point of the monitoring time period and the ending point of the monitoring time period, and by combining the duration of the monitoring time period, which is not specifically limited in the embodiment of the present invention.
In the embodiment of the invention, when the current year crop planting area is determined, the historical crop planting area in the agricultural planting area monitoring area is determined, and then the current year crop planting area in the historical crop planting area is determined. Therefore, the determining efficiency of the crop planting area in the current year can be improved.
On the basis of the above embodiment, the method for extracting an agricultural irrigation area provided in the embodiment of the present invention is to determine a current annual crop planting area in a monitored area of the agricultural planting area based on the second normalized vegetation indexes of the positions of the pixels in the thermal infrared remote sensing image data at the starting point and the ending point, and specifically includes:
for any pixel position in the thermal infrared remote sensing image data, determining a crop planting index of the current year at any pixel position based on the duration of the monitoring time period and the second normalized vegetation indexes of the pixel position at the starting point and the ending point;
and determining the current year crop planting area based on the current year crop planting index at each pixel position in the historical crop planting area.
Specifically, in the embodiment of the present invention, when determining the crop planting area of the current year, for any pixel position in the thermal infrared remote sensing image data of the current year, the crop planting index of the current year at any pixel position may be determined according to the duration of the irrigation area monitoring time period, the second-type normalized vegetation index of any pixel position at the starting point of the irrigation area monitoring time period and the ending point of the irrigation area monitoring time period. Namely, the method comprises the following steps:
Figure BDA0003324448860000171
wherein DZWI is the current year crop planting index, NDVI, at any pixel positiontbA second type normalized vegetation index, NDVI, of the end point tb of the post-planting irrigated area monitoring period of the crop at any pixel positiont0Post-planting irrigation area monitoring for crops at any pixel positionA second type of normalized vegetation index at the start point t0 of the time period. t is tb-t0And monitoring the time period for the irrigation area.
And then determining the crop planting area of the current year according to the crop planting index of the current year at each pixel position in the historical crop planting area. For any pixel position, whether the annual crop planting index at any pixel position is greater than or equal to a third threshold value can be judged, and the third threshold value can be set according to needs, which is not specifically limited in the embodiment of the invention.
If the annual crop planting index at any pixel position is greater than or equal to the third threshold value, it can be determined that any pixel position is a crop pixel position, and all crop pixel positions form the current annual crop planting area.
In the embodiment of the invention, the crop planting index of the current year is determined by monitoring the duration of the time period, the second normalized vegetation index of each pixel position at the starting point of the monitoring time period and the second normalized vegetation index of each pixel position at the ending point of the monitoring time period, and then the crop planting area of the current year is determined, so that the accuracy of the crop planting area of the current year can be ensured, and the determining efficiency can be improved only by considering the second normalized vegetation indexes of each pixel position at the starting point of the monitoring time period and the ending point of the monitoring time period.
On the basis of the above embodiment, the method for extracting an agricultural irrigated area provided by the embodiment of the invention extracts the agricultural irrigated area in the crop planting area of the current year based on the temperature vegetation drought index, and then further includes:
vectorizing the crop planting area in the current year to obtain a plurality of farmland patches, and determining agricultural irrigation areas corresponding to the farmland patches;
and for any farmland patch, if the relative area deviation of the farmland patch and the agricultural irrigation area corresponding to the farmland patch is less than or equal to a deviation threshold value, taking the farmland patch as the agricultural irrigation area.
Specifically, in the embodiment of the present invention, after the agricultural irrigation areas in the crop planting area of the current year are extracted, vectorization processing may be further performed on the crop planting area of the current year, and because the crop planting area of the current year is a farmland area, a plurality of farmland patches may be obtained after vectorization processing. The grid size of each farmland plate is a multiplied by a, and a is the grid side length of the farmland plate. And determining the agricultural irrigation areas corresponding to the patches of the farmlands by combining the extracted agricultural tank areas.
And then regarding any farmland patch, if the relative area deviation of the any farmland patch and the agricultural irrigation area corresponding to the any farmland patch is less than or equal to a deviation threshold value, taking the any farmland patch as the agricultural irrigation area. The relative area deviation can be expressed as the ratio of the difference to the area of the field patch, i.e.:
Figure BDA0003324448860000181
wherein GPI is relative area deviation, S0jThe area of the jth farmland patch can be obtained by statistics, SjThe area of the agricultural tank field corresponding to the jth farmland patch.
If the relative area deviation GPI of the jth farmland patch and the corresponding agricultural irrigation area is less than or equal to the deviation threshold value, the jth farmland patch can be regarded as the agricultural irrigation area, namely, the agricultural tank area is expanded into the farmland patch, and thus the area of the agricultural tank area can be enlarged.
On the basis of the above embodiment, the method for extracting an agricultural irrigation area provided in the embodiment of the present invention for obtaining a historical crop planting area in a monitored area of an agricultural planting area further includes:
acquiring precipitation information of each sub-area in the target area in the monitoring time period of years including the current year;
and determining the agricultural planting area monitoring area in the target area based on the precipitation information.
Specifically, in the embodiment of the invention, when the agricultural planting area monitoring area is determined, the target area is firstly obtained, then the target area is divided into a plurality of sub-areas, and the precipitation information of each sub-area in the historical multi-year monitoring time period including the current year is determined, namely the precipitation information of each sub-area in each historical year in the monitoring time period.
And then determining an agricultural planting area monitoring area in the target area according to the precipitation information of the monitoring time period in each historical year of each sub-area. For any sub-region, the precipitation distance flat percentage of any sub-region in each historical year can be determined according to the precipitation amount information of the monitoring time period of any sub-region in each historical year, namely:
Figure BDA0003324448860000191
wherein PJI is the precipitation flat percentage of any sub-region in each historical year, P is the precipitation information of any sub-region in each historical year monitoring time period,
Figure BDA0003324448860000192
is the average value of the precipitation information of any subregion in the monitoring time period of all historical years.
And then determining an agricultural planting area monitoring area in the target area according to the precipitation range flat percentage of any sub-area in each historical year and the precipitation information of any sub-area in the monitoring time period in each historical year. That is, for any sub-region, whether the precipitation distance flat percentage of the sub-region in each year is smaller than or equal to a fourth threshold value or not is judged, and whether the precipitation information of the sub-region in each historical year is smaller than a fifth threshold value or not is judged. The fourth threshold may be set as needed, which is not specifically limited in the embodiment of the present invention. The fifth threshold may be set as desired, for example, 800 mm.
And if the target year exists, the precipitation distance flat percentage of any sub-area in the target year is smaller than or equal to a fourth threshold, and the precipitation information of any sub-area in the irrigation area monitoring time period in the target year is smaller than a fifth threshold, determining that any sub-area is a candidate agricultural planting area monitoring area, wherein the union set formed by all the candidate agricultural planting area monitoring areas is the agricultural planting area monitoring area.
In the embodiment of the invention, the meteorological information of the precipitation information is combined, so that the monitoring area of the agricultural planting area can be determined more quickly and accurately, and the extraction efficiency of the agricultural irrigation area can be improved.
As shown in fig. 2, on the basis of the above embodiment, an agricultural irrigation area extraction system is provided in an embodiment of the present invention, including:
the first acquisition module 21 is configured to acquire a historical crop planting area in a monitoring area of an agricultural planting area, and determine a current year crop planting area in the monitoring area of the agricultural planting area based on thermal infrared remote sensing image data of the historical crop planting area in the current year;
the second acquisition module 22 is configured to determine thermal infrared temperature data and vegetation coverage data of a time sequence within a monitoring time period based on the thermal infrared remote sensing image data, and perform downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data;
and the extraction module 23 is configured to determine a temperature vegetation drought index of the current year crop planting area based on the refined thermal infrared temperature data and the feature space constructed by the vegetation coverage data, and extract an agricultural irrigation area in the current year crop planting area based on the temperature vegetation drought index.
On the basis of the above embodiment, in the agricultural irrigation area extraction system provided in the embodiment of the present invention, the extraction module is specifically configured to:
and for any pixel of any scene current year thermal infrared remote sensing image in the current year thermal infrared remote sensing image data, determining a temperature vegetation drought index at any pixel based on the thermal infrared temperature data at any pixel, the vegetation coverage data and the wet edge temperature in the characteristic space.
On the basis of the above embodiment, in the agricultural irrigation area extraction system provided in the embodiment of the present invention, the extraction module is further configured to:
for any pixel position in the current year thermal infrared remote sensing image data, determining an irrigated area identification index based on the temperature vegetation drought index of the pixel of each scene current year thermal infrared remote sensing image at the any pixel position in the current year thermal infrared remote sensing image data;
and judging whether any pixel position in the thermal infrared remote sensing image data of the current year is an irrigation area pixel position or not based on the irrigation area identification index, and extracting the agricultural irrigation area in the crop planting area of the current year based on a judgment result.
On the basis of the above embodiment, in the agricultural irrigation area extraction system provided in the embodiment of the present invention, the first obtaining module is specifically configured to:
acquiring historical multi-year remote sensing image data close to the current year, and determining the remote sensing image data of each historical year in the historical multi-year remote sensing image data;
and determining a first-class normalized vegetation index at each pixel position in each scene remote sensing image in the remote sensing image data of each historical year, and determining the historical crop planting area based on the first-class normalized vegetation index at each pixel position in each scene remote sensing image in the remote sensing image data of each historical year.
On the basis of the above embodiment, in the agricultural irrigation area extraction system provided in the embodiment of the present invention, the first obtaining module is further specifically configured to:
for remote sensing image data of any historical year, determining a first type normalized vegetation index mean value at each pixel position in the remote sensing image data of any historical year;
for any scene remote sensing image in the remote sensing image data of any historical year, determining a crop identification index at each pixel in any scene remote sensing image based on the first-class normalized vegetation index and the first-class normalized vegetation index mean value at each pixel in any scene remote sensing image;
judging whether each pixel in any scene remote sensing image is a crop planting pixel or not based on the crop identification index at each pixel in any scene remote sensing image, and determining a crop planting area corresponding to any scene remote sensing image based on a judgment result;
determining a crop planting area corresponding to the remote sensing image data of any historical year based on the crop planting area corresponding to each scene of remote sensing image in the remote sensing image data of any historical year;
and determining the historical crop planting area based on the crop planting area corresponding to the remote sensing image data of each historical year.
On the basis of the above embodiment, in the agricultural irrigation area extraction system provided in the embodiment of the present invention, the first obtaining module is further specifically configured to:
for remote sensing image data of any historical year, determining a first type normalized vegetation index mean value at each pixel position in the remote sensing image data of any historical year;
for any scene remote sensing image in the remote sensing image data of any historical year, determining a crop identification index at each pixel in any scene remote sensing image based on the first-class normalized vegetation index and the first-class normalized vegetation index mean value at each pixel in any scene remote sensing image;
judging whether each pixel in any scene remote sensing image is a crop planting pixel or not based on the crop identification index at each pixel in any scene remote sensing image, and determining a crop planting area corresponding to any scene remote sensing image based on a judgment result;
determining a crop planting area corresponding to the remote sensing image data of any historical year based on the crop planting area corresponding to each scene of remote sensing image in the remote sensing image data of any historical year;
and determining the historical crop planting area based on the crop planting area corresponding to the remote sensing image data of each historical year.
On the basis of the above embodiment, in the agricultural irrigation area extraction system provided in the embodiment of the present invention, the first obtaining module is further specifically configured to:
determining a second type of normalized vegetation index of each pixel position in the thermal infrared remote sensing image data at the starting point and the ending point of the monitoring time period;
and determining a crop planting area of the current year in the agricultural planting area monitoring area based on the second type normalized vegetation index of each pixel position in the thermal infrared remote sensing image data at the starting point and the ending point.
On the basis of the above embodiment, in the agricultural irrigation area extraction system provided in the embodiment of the present invention, the first obtaining module is further specifically configured to:
for any pixel position in the thermal infrared remote sensing image data, determining a crop planting index of the current year at any pixel position based on the duration of the monitoring time period and the second normalized vegetation indexes of the pixel position at the starting point and the ending point;
and determining the current year crop planting area based on the current year crop planting index at each pixel position in the historical crop planting area.
On the basis of the above embodiment, the agricultural irrigation area extraction system provided in the embodiment of the present invention further includes a correction module, configured to:
vectorizing the crop planting area in the current year to obtain a plurality of farmland patches, and determining agricultural irrigation areas corresponding to the farmland patches;
and for any farmland patch, if the relative area deviation of the farmland patch and the agricultural irrigation area corresponding to the farmland patch is less than or equal to a deviation threshold value, taking the farmland patch as the agricultural irrigation area.
On the basis of the above embodiment, the agricultural irrigation area extraction system provided in the embodiment of the present invention further includes a third obtaining module, configured to:
acquiring precipitation information of each sub-area in the target area in the monitoring time period of years including the current year;
and determining the agricultural planting area monitoring area in the target area based on the precipitation information.
Specifically, the functions of the modules in the agricultural irrigation area extraction system provided in the embodiment of the present invention correspond to the operation flows of the steps in the method embodiments one to one, and the implementation effects are also consistent, for which, reference is specifically made to the embodiments described above, which is not described again in the embodiments of the present invention.
Fig. 3 illustrates a physical structure diagram of an electronic device, which may include, as shown in fig. 3: a processor (processor)310, a communication Interface (communication Interface)320, a memory (memory)330 and a communication bus 340, wherein the processor 310, the communication Interface 320 and the memory 330 communicate with each other via the communication bus 340. The processor 310 may call the logic instructions in the memory 330 to execute the agricultural irrigation area extraction method provided in the above embodiments, the method includes: acquiring a historical crop planting area in a monitoring area of an agricultural planting area, and determining a crop planting area of a current year in the monitoring area of the agricultural planting area based on thermal infrared remote sensing image data of the historical crop planting area in the current year; determining thermal infrared temperature data of a time sequence and vegetation coverage data in a monitoring time period based on the thermal infrared remote sensing image data, and performing downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data; and determining a temperature vegetation drought index of the current year crop planting area based on the refined thermal infrared temperature data and the feature space constructed by the vegetation coverage data, and extracting the agricultural irrigation area in the current year crop planting area based on the temperature vegetation drought index.
In addition, the logic instructions in the memory 330 may be implemented in the form of software functional units and stored in a computer readable storage medium when the software functional units are sold or used as independent products. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In another aspect, the present invention also provides a computer program product comprising a computer program, the computer program being stored on a non-transitory computer readable storage medium, the computer program, when executed by a processor, being capable of executing the method for extracting agricultural irrigation areas provided in the above embodiments, the method comprising: acquiring a historical crop planting area in a monitoring area of an agricultural planting area, and determining a crop planting area of a current year in the monitoring area of the agricultural planting area based on thermal infrared remote sensing image data of the historical crop planting area in the current year; determining thermal infrared temperature data of a time sequence and vegetation coverage data in a monitoring time period based on the thermal infrared remote sensing image data, and performing downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data; and determining a temperature vegetation drought index of the current year crop planting area based on the refined thermal infrared temperature data and the feature space constructed by the vegetation coverage data, and extracting the agricultural irrigation area in the current year crop planting area based on the temperature vegetation drought index.
In yet another aspect, the present invention also provides a non-transitory computer readable storage medium having stored thereon a computer program, which when executed by a processor, is implemented to perform the agricultural irrigated area extraction method provided in the above embodiments, the method comprising: acquiring a historical crop planting area in a monitoring area of an agricultural planting area, and determining a crop planting area of a current year in the monitoring area of the agricultural planting area based on thermal infrared remote sensing image data of the historical crop planting area in the current year; determining thermal infrared temperature data of a time sequence and vegetation coverage data in a monitoring time period based on the thermal infrared remote sensing image data, and performing downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data; and determining a temperature vegetation drought index of the current year crop planting area based on the refined thermal infrared temperature data and the feature space constructed by the vegetation coverage data, and extracting the agricultural irrigation area in the current year crop planting area based on the temperature vegetation drought index.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
Through the above description of the embodiments, those skilled in the art will clearly understand that each embodiment can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware. With this understanding in mind, the above-described technical solutions may be embodied in the form of a software product, which can be stored in a computer-readable storage medium such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the methods described in the embodiments or some parts of the embodiments.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (10)

1. An agricultural irrigated area extraction method is characterized by comprising the following steps:
acquiring a historical crop planting area in a monitoring area of an agricultural planting area, and determining a crop planting area of a current year in the monitoring area of the agricultural planting area based on thermal infrared remote sensing image data of the historical crop planting area in the current year;
determining thermal infrared temperature data of a time sequence and vegetation coverage data in a monitoring time period based on the thermal infrared remote sensing image data, and performing downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data;
and determining a temperature vegetation drought index of the current year crop planting area based on the refined thermal infrared temperature data and the feature space constructed by the vegetation coverage data, and extracting the agricultural irrigation area in the current year crop planting area based on the temperature vegetation drought index.
2. The method for extracting agricultural irrigated area according to claim 1, wherein the determining the temperature vegetation drought index of the current year crop planting area based on the feature space constructed by the refined thermal infrared temperature data and the vegetation coverage data specifically comprises:
and for any pixel of any scene current year thermal infrared remote sensing image in the current year thermal infrared remote sensing image data, determining a temperature vegetation drought index at any pixel based on the thermal infrared temperature data at any pixel, the vegetation coverage data and the wet edge temperature in the characteristic space.
3. The method for extracting agricultural irrigated area according to claim 1, wherein the extracting the agricultural irrigated area in the crop planting area of the current year based on the temperature vegetation drought index specifically comprises:
for any pixel position in the current year thermal infrared remote sensing image data, determining an irrigated area identification index based on the temperature vegetation drought index of the pixel of each scene current year thermal infrared remote sensing image at the any pixel position in the current year thermal infrared remote sensing image data;
and judging whether any pixel position in the thermal infrared remote sensing image data of the current year is an irrigation area pixel position or not based on the irrigation area identification index, and extracting the agricultural irrigation area in the crop planting area of the current year based on a judgment result.
4. The agricultural irrigated area extraction method of claim 1, wherein the obtaining of the historic crop planting area within the agricultural planting area monitoring area specifically comprises:
acquiring historical multi-year remote sensing image data close to the current year, and determining the remote sensing image data of each historical year in the historical multi-year remote sensing image data;
and determining a first-class normalized vegetation index at each pixel position in each scene remote sensing image in the remote sensing image data of each historical year, and determining the historical crop planting area based on the first-class normalized vegetation index at each pixel position in each scene remote sensing image in the remote sensing image data of each historical year.
5. The agricultural irrigated area extraction method of claim 4, wherein the determining the historical crop planting area based on the first normalized vegetation index at each pixel in each scene remote sensing image in the remote sensing image data of each historical year specifically comprises:
for remote sensing image data of any historical year, determining a first type normalized vegetation index mean value at each pixel position in the remote sensing image data of any historical year;
for any scene remote sensing image in the remote sensing image data of any historical year, determining a crop identification index at each pixel in any scene remote sensing image based on the first-class normalized vegetation index and the first-class normalized vegetation index mean value at each pixel in any scene remote sensing image;
judging whether each pixel in any scene remote sensing image is a crop planting pixel or not based on the crop identification index at each pixel in any scene remote sensing image, and determining a crop planting area corresponding to any scene remote sensing image based on a judgment result;
determining a crop planting area corresponding to the remote sensing image data of any historical year based on the crop planting area corresponding to each scene of remote sensing image in the remote sensing image data of any historical year;
and determining the historical crop planting area based on the crop planting area corresponding to the remote sensing image data of each historical year.
6. The method for extracting agricultural irrigated area according to claim 1, wherein the determining the current year crop planting area in the monitoring area of the agricultural planting area based on the thermal infrared remote sensing image data of the historical crop planting area in the current year specifically comprises:
determining a second type of normalized vegetation index of each pixel position in the thermal infrared remote sensing image data at the starting point and the ending point of the monitoring time period;
and determining a crop planting area of the current year in the agricultural planting area monitoring area based on the second type normalized vegetation index of each pixel position in the thermal infrared remote sensing image data at the starting point and the ending point.
7. The method for extracting agricultural irrigated areas according to claim 6, wherein the determining the crop planting area of the current year in the monitoring area of the agricultural planting area based on the second normalized vegetation index of each pixel position in the thermal infrared remote sensing image data at the starting point and the ending point specifically comprises:
for any pixel position in the thermal infrared remote sensing image data, determining a crop planting index of the current year at any pixel position based on the duration of the monitoring time period and the second normalized vegetation indexes of the pixel position at the starting point and the ending point;
and determining the current year crop planting area based on the current year crop planting index at each pixel position in the historical crop planting area.
8. The method of extracting agricultural irrigated area according to any one of claims 1-7, wherein the extracting the agricultural irrigated area within the current year crop planting area based on the temperature vegetation drought index further comprises:
vectorizing the crop planting area in the current year to obtain a plurality of farmland patches, and determining agricultural irrigation areas corresponding to the farmland patches;
and for any farmland patch, if the relative area deviation of the farmland patch and the agricultural irrigation area corresponding to the farmland patch is less than or equal to a deviation threshold value, taking the farmland patch as the agricultural irrigation area.
9. The agricultural irrigated area extraction method of any one of claims 1 to 7, wherein the obtaining of the historic crop planting area within the monitored area of the agricultural planting area further comprises:
acquiring precipitation information of each sub-area in the target area in the monitoring time period of years including the current year;
and determining the agricultural planting area monitoring area in the target area based on the precipitation information.
10. An agricultural irrigated area extraction system, comprising:
the system comprises a first acquisition module, a second acquisition module and a third acquisition module, wherein the first acquisition module is used for acquiring a historical crop planting area in a monitoring area of an agricultural planting area and determining a current year crop planting area in the monitoring area of the agricultural planting area based on thermal infrared remote sensing image data of the historical crop planting area in the current year;
the second acquisition module is used for determining thermal infrared temperature data and vegetation coverage data of a time sequence in a monitoring time period based on the thermal infrared remote sensing image data, and performing downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data;
and the extraction module is used for determining the temperature vegetation drought index of the current year crop planting area based on the refined thermal infrared temperature data and the feature space constructed by the vegetation coverage data, and extracting the agricultural irrigation area in the current year crop planting area based on the temperature vegetation drought index.
CN202111257130.XA 2021-10-27 2021-10-27 Agricultural irrigation area extraction method and system Active CN113984212B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111257130.XA CN113984212B (en) 2021-10-27 2021-10-27 Agricultural irrigation area extraction method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111257130.XA CN113984212B (en) 2021-10-27 2021-10-27 Agricultural irrigation area extraction method and system

Publications (2)

Publication Number Publication Date
CN113984212A true CN113984212A (en) 2022-01-28
CN113984212B CN113984212B (en) 2023-06-27

Family

ID=79742724

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111257130.XA Active CN113984212B (en) 2021-10-27 2021-10-27 Agricultural irrigation area extraction method and system

Country Status (1)

Country Link
CN (1) CN113984212B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114724029A (en) * 2022-04-02 2022-07-08 中科海慧（天津）科技有限公司 Remote sensing extraction method and system of mining goaf's influence on the surface
CN114858219A (en) * 2022-06-09 2022-08-05 西安微媒软件有限公司 Agricultural Internet of things management and analysis platform for greenhouse
CN115063328A (en) * 2022-05-05 2022-09-16 广东省深圳生态环境监测中心站(广东省东江流域生态环境监测中心) Thermal infrared image reconstruction method based on feature space interpolation
CN119196561A (en) * 2024-09-14 2024-12-27 铜陵有色金属集团股份有限公司动力厂 A method for monitoring leakage of overhead thermal pipelines based on high-resolution inversion of ground surface temperature

Citations (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2722251A1 (en) * 1994-07-05 1996-01-12 Transenergie Solar powered water pumping station for use in remote locations such as semi-arid zones
US5764819A (en) * 1991-10-18 1998-06-09 Dekalb Genetics Corporation Methods for classifying plants for evaluation and breeding programs by use of remote sensing and image analysis technology
US7058197B1 (en) * 1999-11-04 2006-06-06 Board Of Trustees Of The University Of Illinois Multi-variable model for identifying crop response zones in a field
US20080238423A1 (en) * 2007-03-29 2008-10-02 Guo Bin Li Method and apparatus for accelerating magnetic resonance temperature imaging
CN102176002A (en) * 2010-12-30 2011-09-07 中国科学院地理科学与资源研究所 Surface water heat flux remote sensing inversion-based drought monitoring method and system
KR101570297B1 (en) * 2014-08-01 2015-11-18 성균관대학교산학협력단 System, method and computer program for drought mapping using satellite imageries in northeast asia
CN105469195A (en) * 2015-11-18 2016-04-06 国家电网公司 Power transmission line corridor environment fire danger class evaluation method
CN105760978A (en) * 2015-07-22 2016-07-13 北京师范大学 Agricultural drought grade monitoring method based on temperature vegetation drought index (TVDI)
US20160290918A1 (en) * 2014-09-12 2016-10-06 The Climate Corporation Forecasting national crop yield during the growing season
CN106226260A (en) * 2016-08-10 2016-12-14 武汉大学 A kind of combination microwave and the Soil Moisture Inversion method of infrared remote sensing image
CN107782701A (en) * 2017-09-20 2018-03-09 北京师范大学 A kind of agricultural arid monitoring method of multi- source Remote Sensing Data data
US20180189564A1 (en) * 2016-12-30 2018-07-05 International Business Machines Corporation Method and system for crop type identification using satellite observation and weather data
CN108955620A (en) * 2018-02-13 2018-12-07 中国科学院遥感与数字地球研究所 A kind of method and system of farmland irrigated area area Remotely sensed acquisition
US20190230875A1 (en) * 2015-08-05 2019-08-01 Clearag, Inc. Customized land surface modeling in a soil-crop system using satellite data to detect irrigation and precipitation events for decision support in precision agriculture
CN110175932A (en) * 2019-05-16 2019-08-27 中国农业科学院农业资源与农业区划研究所 The Agriculture Drought monitoring method coupled with more ten days precipitation is extracted based on remote sensing information
CN111199185A (en) * 2019-11-26 2020-05-26 广州地理研究所 Surface temperature downscaling method, system and device based on XGBoost learning algorithm
WO2020112013A1 (en) * 2018-11-27 2020-06-04 Chawade Aakash Methods, models and systems for predicting yellow rust in wheat crops
CN111582554A (en) * 2020-04-17 2020-08-25 北京农业信息技术研究中心 Crop growth prediction method and system
CN111610159A (en) * 2020-05-11 2020-09-01 中科海慧（天津）科技有限公司 Earth surface temperature downscaling estimation method and vegetation water stress monitoring method
CN112232230A (en) * 2020-10-20 2021-01-15 中国气象科学研究院 Method and device for remote sensing identification of abandoned farmland information
CN112418050A (en) * 2020-11-18 2021-02-26 中国科学院空天信息创新研究院 Method and device for remote sensing identification of abandoned farmland information
CN112858632A (en) * 2021-01-14 2021-05-28 中国科学院空天信息创新研究院 Grassland drought monitoring method comprehensively considering temperature and water stress
CN113010994A (en) * 2021-01-21 2021-06-22 成都信息工程大学 Regional evapotranspiration monitoring method based on surface evapotranspiration remote sensing
CN113189015A (en) * 2021-04-23 2021-07-30 宿州学院 Hyperspectral remote sensing-based insect pest monitoring system and method
CN113421273A (en) * 2021-06-30 2021-09-21 中国气象科学研究院 Remote sensing extraction method and device for forest and grass collocation information

Patent Citations (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5764819A (en) * 1991-10-18 1998-06-09 Dekalb Genetics Corporation Methods for classifying plants for evaluation and breeding programs by use of remote sensing and image analysis technology
FR2722251A1 (en) * 1994-07-05 1996-01-12 Transenergie Solar powered water pumping station for use in remote locations such as semi-arid zones
US7058197B1 (en) * 1999-11-04 2006-06-06 Board Of Trustees Of The University Of Illinois Multi-variable model for identifying crop response zones in a field
US20080238423A1 (en) * 2007-03-29 2008-10-02 Guo Bin Li Method and apparatus for accelerating magnetic resonance temperature imaging
CN102176002A (en) * 2010-12-30 2011-09-07 中国科学院地理科学与资源研究所 Surface water heat flux remote sensing inversion-based drought monitoring method and system
KR101570297B1 (en) * 2014-08-01 2015-11-18 성균관대학교산학협력단 System, method and computer program for drought mapping using satellite imageries in northeast asia
US20160290918A1 (en) * 2014-09-12 2016-10-06 The Climate Corporation Forecasting national crop yield during the growing season
CN105760978A (en) * 2015-07-22 2016-07-13 北京师范大学 Agricultural drought grade monitoring method based on temperature vegetation drought index (TVDI)
US20190230875A1 (en) * 2015-08-05 2019-08-01 Clearag, Inc. Customized land surface modeling in a soil-crop system using satellite data to detect irrigation and precipitation events for decision support in precision agriculture
CN105469195A (en) * 2015-11-18 2016-04-06 国家电网公司 Power transmission line corridor environment fire danger class evaluation method
CN106226260A (en) * 2016-08-10 2016-12-14 武汉大学 A kind of combination microwave and the Soil Moisture Inversion method of infrared remote sensing image
US20180189564A1 (en) * 2016-12-30 2018-07-05 International Business Machines Corporation Method and system for crop type identification using satellite observation and weather data
CN107782701A (en) * 2017-09-20 2018-03-09 北京师范大学 A kind of agricultural arid monitoring method of multi- source Remote Sensing Data data
CN108955620A (en) * 2018-02-13 2018-12-07 中国科学院遥感与数字地球研究所 A kind of method and system of farmland irrigated area area Remotely sensed acquisition
WO2020112013A1 (en) * 2018-11-27 2020-06-04 Chawade Aakash Methods, models and systems for predicting yellow rust in wheat crops
CN110175932A (en) * 2019-05-16 2019-08-27 中国农业科学院农业资源与农业区划研究所 The Agriculture Drought monitoring method coupled with more ten days precipitation is extracted based on remote sensing information
CN111199185A (en) * 2019-11-26 2020-05-26 广州地理研究所 Surface temperature downscaling method, system and device based on XGBoost learning algorithm
CN111582554A (en) * 2020-04-17 2020-08-25 北京农业信息技术研究中心 Crop growth prediction method and system
CN111610159A (en) * 2020-05-11 2020-09-01 中科海慧（天津）科技有限公司 Earth surface temperature downscaling estimation method and vegetation water stress monitoring method
CN112232230A (en) * 2020-10-20 2021-01-15 中国气象科学研究院 Method and device for remote sensing identification of abandoned farmland information
CN112418050A (en) * 2020-11-18 2021-02-26 中国科学院空天信息创新研究院 Method and device for remote sensing identification of abandoned farmland information
CN112858632A (en) * 2021-01-14 2021-05-28 中国科学院空天信息创新研究院 Grassland drought monitoring method comprehensively considering temperature and water stress
CN113010994A (en) * 2021-01-21 2021-06-22 成都信息工程大学 Regional evapotranspiration monitoring method based on surface evapotranspiration remote sensing
CN113189015A (en) * 2021-04-23 2021-07-30 宿州学院 Hyperspectral remote sensing-based insect pest monitoring system and method
CN113421273A (en) * 2021-06-30 2021-09-21 中国气象科学研究院 Remote sensing extraction method and device for forest and grass collocation information

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ARAYA Z.GHEBREAMLAK ET,AL: "Satellite-Based Mapping of Cultivated Area in Gash Delta Spate Irrigation System,Sudan", 《REMOTE SENSING》, vol. 10, no. 2, pages 1 - 14 *
SELVER SENTURK ET,AL: "Application of Remote Sensing Techniques in Locating Dry and Irrigated Farmland Parcels", 《IEEE》, pages 1 - 4 *
张智韬 等: "泾惠渠灌区作物种植结构变化对灌溉需水量的影像", 《农业机械学报》, vol. 47, no. 10, pages 122 - 130 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114724029A (en) * 2022-04-02 2022-07-08 中科海慧（天津）科技有限公司 Remote sensing extraction method and system of mining goaf's influence on the surface
CN115063328A (en) * 2022-05-05 2022-09-16 广东省深圳生态环境监测中心站(广东省东江流域生态环境监测中心) Thermal infrared image reconstruction method based on feature space interpolation
CN114858219A (en) * 2022-06-09 2022-08-05 西安微媒软件有限公司 Agricultural Internet of things management and analysis platform for greenhouse
CN119196561A (en) * 2024-09-14 2024-12-27 铜陵有色金属集团股份有限公司动力厂 A method for monitoring leakage of overhead thermal pipelines based on high-resolution inversion of ground surface temperature

Also Published As

Publication number Publication date
CN113984212B (en) 2023-06-27

Similar Documents

Publication Publication Date Title
CN113984212A (en) Agricultural irrigated area extraction method and system
KR102540762B1 (en) Reservoir monitoring method using satellite informations
Vélez et al. Mapping the spatial variability of Botrytis bunch rot risk in vineyards using UAV multispectral imagery
US9430499B2 (en) Automated feature extraction from imagery
CN109635731B (en) A method and device for identifying effective cultivated land, storage medium and processor
Han et al. The RapeseedMap10 database: annual maps of rapeseed at a spatial resolution of 10 m based on multi-source data
Bai et al. A fast and robust method for plant count in sunflower and maize at different seedling stages using high-resolution UAV RGB imagery
CN110288647A (en) A Method for Monitoring Irrigated Area of Irrigated Areas Based on High Resolution Satellite Data
CN110942039A (en) Remote sensing monitoring and forecasting system and method for high-temperature disasters of main crops
CN113223040B (en) Banana estimated yield method and device based on remote sensing, electronic equipment and storage medium
CN109522899B (en) Detection method and device for ripe coffee fruits and electronic equipment
CN111610159B (en) Earth surface temperature downscaling estimation method and vegetation water stress monitoring method
CN114419367B (en) A method and system for high-precision crop mapping
CN107680098A (en) A kind of recognition methods of sugarcane sugarcane section feature
Carmona et al. Grazing management or physiography? Factors controlling vegetation recovery in Mediterranean grasslands
CN107437262B (en) Crop Planting Area Early Warning Method and System
CN112585505A (en) Determining location-specific weather information for agronomic decision support
CN112418016B (en) SAR-based irrigation information extraction method and device
CN112001809A (en) Method for acquiring farmland returning information of agriculture and forestry area
CN116129258A (en) Citrus land parcel identification method and system based on canopy density and multisource remote sensing image
CN113570273B (en) Spatialization method and system for irrigation farmland statistical data
CN113421273B (en) Remote sensing extraction method and device for forest and grass collocation information
CN102621075A (en) Method for automatically detecting rice heading stage
CN116912702B (en) Weed coverage determination method, system and device and electronic equipment
CN111583050A (en) Crop pest and disease early warning method and system fusing multi-scale remote sensing images

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