CN113984212B - Agricultural irrigation area extraction method and system - Google Patents

Abstract

The invention provides an agricultural irrigation area extraction method and system, wherein in the extraction process, the temperature vegetation drought index of a current year crop planting area is determined by performing downscaling treatment on thermal infrared remote sensing image data of a current year of a historical crop planting area and according to a characteristic space constructed by refined thermal infrared temperature data and vegetation coverage data, so that the aim of accurately extracting the agricultural irrigation area can be fulfilled.

Description

Method and system for extracting agricultural irrigation areas

technical field

The invention relates to the technical field of information extraction, in particular to an extraction method and system for agricultural irrigation areas.

Background technique

Irrigated areas refer to artificially irrigated areas. The area and distribution information of agricultural irrigated areas located in arid and semi-arid areas is very important for crop growth and yield estimation. In addition, due to the relative scarcity of water resources (including precipitation, groundwater, river and lake water resources, etc.) important.

At present, the extraction methods of agricultural irrigation area information in arid and semi-arid areas usually include manual statistical methods and remote sensing methods. The artificial statistical method mainly obtains the information of agricultural irrigation areas through artificial field statistics and summary reporting. The problem of this method is that it is not objective enough and it is difficult to obtain the spatial and temporal dynamic distribution information of agricultural irrigation areas. The remote sensing method is the most commonly used method at home and abroad, with the advantages of objectivity, rapidity, and low cost, including land use analysis methods, growth monitoring methods, and soil moisture analysis methods. The land use analysis method is to obtain cultivated land information through remote sensing technology and estimate it combined with experience. This method can effectively obtain cultivated land distribution information, but it needs a lot of experience to identify. The growth monitoring method mainly uses remote sensing information observed at different times, such as canopy moisture analysis method, vegetation index time series analysis method, etc. The canopy moisture analysis method quickly and accurately identifies the canopy moisture information and compares the canopy information in other areas to determine whether it is an agricultural irrigation area. However, this method requires relatively high observation time points and spectral resolution of remote sensing data; vegetation The index time series analysis method mainly analyzes and identifies the growth difference between irrigated and non-irrigated areas through the time series vegetation index. The soil moisture analysis method mainly uses microwave data or thermal infrared data to analyze soil moisture. The mechanism of this method is clear, but limited by the spatial resolution, it is difficult to apply to fragmented or small areas.

For this reason, there is an urgent need to provide a method for extracting agricultural irrigation areas.

Contents of the invention

The invention provides a method and system for extracting agricultural irrigation areas to solve the defects in the prior art.

The invention provides a method for extracting agricultural irrigation areas, comprising:

Obtain the historical crop planting area in the agricultural planting area monitoring area, and determine the current year crop planting area in the agricultural planting area monitoring area based on the thermal infrared remote sensing image data of the historical crop planting area in the current year;

Based on the thermal infrared remote sensing image data, determine thermal infrared temperature data and vegetation coverage data of time series within the monitoring period, and perform downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data;

Based on the feature space constructed by the refined thermal infrared temperature data and the vegetation coverage data, the temperature vegetation drought index of the crop planting area of the current year is determined, and based on the temperature vegetation drought index, the current year crop The agricultural irrigation area within the planting area is used for extraction.

According to a method for extracting agricultural irrigation areas provided by the present invention, the feature space constructed based on the refined thermal infrared temperature data and the vegetation coverage data is used to determine the temperature vegetation drought index of the crop planting area in the current year, specifically include:

For any pixel of any scene in the thermal infrared remote sensing image data of the current year in the current year thermal infrared remote sensing image, based on the thermal infrared temperature data at any pixel, the vegetation coverage data and the features Wet edge temperature in the space, determine the temperature vegetation drought index at any pixel.

According to a method for extracting agricultural irrigation areas provided by the present invention, the extraction of agricultural irrigation areas in the crop planting area of the current year based on the temperature vegetation drought index specifically includes:

For any pixel position in the thermal infrared remote sensing image data of the current year, based on the temperature and vegetation of the pixels of the thermal infrared remote sensing image of each scene in the current year thermal infrared remote sensing image at the arbitrary pixel position in the thermal infrared remote sensing image data of the current year Drought index, to determine the irrigated area identification index;

Based on the irrigated area identification index, it is judged whether any pixel position in the thermal infrared remote sensing image data of the current year is an irrigated area pixel position, and based on the judgment result, the agricultural irrigation area in the crop planting area of the current year is determined. extract.

According to a method for extracting agricultural irrigation areas provided by the present invention, the acquisition of historical crop planting areas in the monitoring area of agricultural planting areas specifically includes:

Obtaining the 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;

Determining the first-type normalized difference vegetation index at each pixel in each scene remote sensing image in the remote sensing image data of each historical year, and based on each pixel in each scene remote sensing image in the remote sensing image data of each historical year The first type of normalized difference vegetation index at , to determine the historical crop planting area.

According to a method for extracting agricultural irrigation areas provided by the present invention, the first-type normalized difference vegetation index at each pixel in each scene remote sensing image in the remote sensing image data of each historical year is used to determine the planting of the historical crops. areas, including:

For the remote sensing image data of any historical year, determine the average value of the first normalized difference vegetation index at each pixel position in the remote sensing image data of any historical year;

For any remote sensing image of any scene in the remote sensing image data of any historical year, based on the first type of normalized normalized vegetation index at each pixel in the remote sensing image of any scene and the first type of normalized vegetation The index mean value determines the crop identification index at each pixel in the remote sensing image of any scene;

Based on the crop recognition index at each pixel in the remote sensing image of any scene, judge whether each pixel in the remote sensing image of any scene is a crop planting pixel, and based on the judgment result, determine the remote sensing image of any scene The corresponding crop growing area;

Based on the crop planting area corresponding to each remote sensing image in the remote sensing image data of any historical year, determine the crop planting area corresponding to the remote sensing image data of any historical year;

The historical crop planting area is determined based on the crop planting area corresponding to the remote sensing image data of each historical year.

According to a method for extracting agricultural irrigation areas provided by the present invention, the determination of the crop planting area of the current year 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 includes:

Determining the second normalized difference vegetation index at the start point and end point of each pixel position in the thermal infrared remote sensing image data at the monitoring time period;

Based on the second normalized normalized vegetation index of each pixel position in the thermal infrared remote sensing image data at the start point and the end point, the crop planting area of the current year in the monitoring area of the agricultural planting area is determined.

According to a method for extracting agricultural irrigation areas provided by the present invention, the second type of normalized normalized vegetation index based on the position of each pixel in the thermal infrared remote sensing image data at the starting point and the ending point is determined to determine the The crop planting area of the current year within the agricultural planting area monitoring area, including:

For any pixel position in the thermal infrared remote sensing image data, based on the length of the monitoring time period, the second normalized vegetation of the any pixel position at the start point and the end point Index, determine the crop planting index of the current year at any pixel position;

The crop planting area of the current year is determined based on the crop planting index of the current year at each pixel position in the historical crop planting area.

According to a method for extracting agricultural irrigation areas provided by the present invention, the agricultural irrigation areas in the crop planting area of the current year are extracted based on the temperature vegetation drought index, and then further include:

Carry out vectorization processing on the crop planting area of the current year to obtain several farmland patches, and determine the agricultural irrigation area corresponding to each farmland patch;

For any patch of farmland, if the relative area deviation between the patch of any farmland and the agricultural irrigation area corresponding to the patch of any farmland is less than or equal to a deviation threshold, then the patch of any farmland is regarded as an irrigation area of agriculture.

According to a method for extracting agricultural irrigation areas provided by the present invention, the acquisition of historical crop planting areas in the monitoring area of agricultural planting areas also includes:

Obtain the precipitation information of the target area in the monitoring time period of many years including the current year;

Based on the precipitation information, the monitoring area of the agricultural planting area in the target area is determined.

The present invention also provides an extraction system for agricultural irrigation areas, comprising:

The first acquisition module is used to acquire the historical crop planting area in the agricultural planting area monitoring area, and determine 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 in the current year crop growing areas;

The second acquisition module is used to determine the time-series thermal infrared temperature data and vegetation coverage data within the 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;

The extraction module is used to determine 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 based on the temperature vegetation drought index, the The agricultural irrigation area within the crop planting area of the current year is extracted.

The present invention also provides an electronic device, including a memory, a processor, and a computer program stored on the memory and operable on the processor. When the processor executes the program, the extraction of agricultural irrigation areas as described in any of the above-mentioned method steps.

The present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of any method for extracting agricultural irrigation areas described above are realized.

The present invention also provides a computer program product, including a computer program. When the computer program is executed by a processor, the steps of any method for extracting agricultural irrigation areas described above are realized.

The method and system for extracting agricultural irrigation areas provided by the present invention first obtain the historical crop planting areas in the agricultural planting area monitoring area, and determine the current crop planting area in the agricultural planting area monitoring area based on the thermal infrared remote sensing image data of the historical crop planting area in the current year. The annual crop planting area; then based on the thermal infrared remote sensing image data, determine the thermal infrared temperature data and vegetation coverage data of the time series within the monitoring period, and downscale the thermal infrared temperature data to obtain refined thermal infrared temperature data; Finally, based on the feature space constructed by refined thermal infrared temperature data and vegetation coverage data, the temperature vegetation drought index of the crop planting area in the current year is determined, and based on the temperature vegetation drought index, the agricultural irrigation areas in the crop planting area of the current year are extracted. In the extraction process, by downscaling the thermal infrared remote sensing image data of the historical crop planting area in the current year, and according to the feature space constructed by the refined thermal infrared temperature data and vegetation coverage data, the crop planting area of the current year is determined. The temperature vegetation drought index can achieve the purpose of accurately extracting agricultural irrigation areas.

Description of drawings

In order to more clearly illustrate the present invention or the technical solutions in the prior art, the accompanying drawings that need to be used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the accompanying drawings in the following description are the present invention. For some embodiments of the invention, those skilled in the art can also obtain other drawings based on these drawings without creative effort.

Fig. 1 is the schematic flow chart of the extraction method of agricultural irrigation area provided by the present invention;

Fig. 2 is the structural representation of the agricultural irrigation area extraction system provided by the present invention;

Fig. 3 is a schematic structural diagram of an electronic device provided by the present invention.

Detailed ways

In order to make the purpose, technical solutions and advantages of the present invention clearer, the technical solutions in the present invention will be clearly and completely described below in conjunction with the accompanying drawings in the present invention. Obviously, the described embodiments are part of the embodiments of the present invention , but not all examples. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without creative efforts fall within the protection scope of the present invention.

Fig. 1 is a schematic flow chart of a method for extracting agricultural irrigation areas provided in an embodiment of the present invention. As shown in Fig. 1, the method includes:

S1. Obtain the historical crop planting area in the monitoring area of the agricultural planting area, and determine the crop planting area of the current year 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;

S2. Based on the thermal infrared remote sensing image data, determine thermal infrared temperature data and vegetation coverage data of time series within the monitoring period, and perform downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data;

S3, based on the feature space constructed by the refined thermal infrared temperature data and the vegetation coverage data, determine the temperature vegetation drought index of the crop planting area in the current year, and based on the temperature vegetation drought index, calculate the current The agricultural irrigation area within the annual crop planting area is extracted.

Specifically, the method for extracting agricultural irrigation areas provided in the embodiments of the present invention is executed by a server. The server may be a local server or a cloud server. The local server may specifically be a computer, a tablet computer, or a smart phone. This is not specifically limited in the embodiments.

Firstly, step S1 is executed to obtain the historical crop planting area in the agricultural planting area monitoring area. The agricultural planting area monitoring area refers to the area in the target area that needs to be monitored whether it is an agricultural irrigation area. The historical crop planting area can be the union of crop planting areas in many years of history, can be determined by land use information, or can be obtained by extracting remote sensing image data for many years in history, which is not specifically limited in the embodiments of the present invention. The number of years involved in the historical years can be selected as required, for example, 10 years, 9 years, etc. can be selected.

Then the thermal infrared remote sensing image data of the historical crop planting area in the current year can be obtained. The thermal infrared remote sensing image data of the current year can include multiple thermal infrared remote sensing images, and each thermal infrared remote sensing image can include reflectance data at each pixel. and thermal infrared temperature data.

Furthermore, according to the thermal infrared remote sensing image data of the current year, the crop planting area of the current year in the monitoring area of the agricultural planting area can be determined. In the embodiment of the present invention, the area change of the crop planting area in the monitoring area of the agricultural planting 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 monitoring area of the agricultural planting area can be determined. The crop planting area of the current year can also be determined from the historical crop planting area based on the thermal infrared remote sensing image data of the current year.

Then step S2 is executed to determine the time series thermal infrared temperature data and vegetation coverage data within the monitoring period according to the thermal infrared remote sensing image data of the current year. The monitoring period can be determined by the following methods:

Obtain the phenological curve of the monitoring area of the agricultural planting area; wherein, the phenological curve is determined according to the crops planted in the monitoring area of the agricultural planting area, and different crops correspond to different phenological curves.

The monitoring time period of the current year is determined according to the phenological curve; the monitoring time period of the irrigated area in the current year may be the active time period of the crop, for example, it may be the entire time period from the beginning of crop seedling growth to the end of crop maturity.

Determine the thermal infrared temperature data and vegetation coverage data of the time series within the monitoring period, that is, the thermal infrared temperature data and vegetation coverage data of the time series at each pixel position in the thermal infrared remote sensing image data of the current year can be obtained.

Since the obtained thermal infrared temperature data usually has a low resolution, in order to improve its resolution, the thermal infrared temperature data can be downscaled, that is, downsampled, and the processing method can be realized by the following formula:

LST _high ＝F(PRS _high )+LST _low -F(PRS _low )

Among them, LST _high and LST _low are the simulated high-resolution land surface temperature and low-resolution land surface temperature, respectively; LST _sd and LST _sw are the highest temperature on the bare soil surface when the vegetation coverage rate is 0 and the vegetation coverage under full vegetation coverage, respectively. The maximum temperature on the surface; k _dry and k _wet are the slopes of the simulated surface temperature and the dry line and wet line in the feature space determined by the vegetation coverage data; β is the Bowen ratio; FVC is the vegetation coverage, FVC _high and FVC _low are simulated high-resolution vegetation coverage and low-resolution vegetation coverage respectively; LST is observed thermal infrared data; ε is error; F(PRS _high ) and F(PRS _low ) are high-resolution and low-resolution The simulation function of low-resolution thermal infrared remote sensing images, PRS _high and PRS _low are simulated high-resolution and low-resolution thermal infrared remote sensing images respectively, and PRS includes PRS _high and PRS _low .

Let the value of LST _low in the above formula be the thermal infrared temperature data, and the obtained value of LST _high is the refined thermal infrared temperature data.

Finally, step S3 is executed to construct a feature space according to the obtained refined thermal infrared temperature data and vegetation coverage data. That is, in a two-dimensional coordinate system with the refined thermal infrared temperature data as the vertical axis and the vegetation coverage data as the horizontal axis, draw the relationship between the refined thermal infrared temperature data and the vegetation coverage data to obtain the feature space. In this feature space, there are wet lines and dry lines, the wet lines are parallel to the horizontal axis, the end points of the dry lines coincide with the end points of the wet lines, and the slope of the dry lines is negative.

According to the constructed feature space, combined with thermal infrared temperature data, the temperature vegetation drought index of the crop planting area in the current year can be determined. The temperature vegetation drought index can be used to characterize whether there is vegetation drought in the crop planting area of the current year. Each pixel of the current year thermal infrared remote sensing image data in the current year thermal infrared remote sensing image corresponds to a temperature vegetation drought index. The temperature vegetation drought index at each pixel can be determined through thermal infrared temperature data and vegetation coverage data at each pixel, which is not specifically limited in this embodiment 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 of the current year can be extracted. When extracting agricultural irrigation areas, the temperature and vegetation drought index at each pixel of each scene in the current year thermal infrared remote sensing image data can be used to determine the temperature and vegetation drought index of each image in the current year thermal infrared remote sensing image data. According to the irrigated area identification index at the pixel position, it can be determined whether each pixel position is the irrigated area pixel position, and then the agricultural irrigation area in the crop planting area of the current year can be extracted.

The agricultural irrigation area extraction method provided in the embodiment of the present invention first obtains the historical crop planting area in the monitoring area of the agricultural planting area, and determines the 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. The crop planting area in the current year; then based on the thermal infrared remote sensing image data, determine the thermal infrared temperature data and vegetation coverage data of the time series within the monitoring period, and downscale the thermal infrared temperature data to obtain refined thermal infrared temperature data ; Finally, based on the feature space constructed by refined thermal infrared temperature data and vegetation coverage data, the temperature vegetation drought index of the crop planting area in the current year is determined, and based on the temperature vegetation drought index, the agricultural irrigation area in the crop planting area of the current year is extracted . In the extraction process, by downscaling the thermal infrared remote sensing image data of the historical crop planting area in the current year, and according to the feature space constructed by the refined thermal infrared temperature data and vegetation coverage data, the crop planting area of the current year is determined. The temperature vegetation drought index can achieve the purpose of accurately extracting agricultural irrigation areas.

On the basis of the above-mentioned embodiments, in the method for extracting agricultural irrigation areas provided in the embodiments of the present invention, the feature space constructed based on the refined thermal infrared temperature data and the vegetation coverage data is used to determine the current year crop planting The temperature and vegetation drought index of the region, including:

For any pixel of any scene in the thermal infrared remote sensing image data of the current year in the current year thermal infrared remote sensing image, based on the thermal infrared temperature data at any pixel, the vegetation coverage data and the features Wet edge temperature in the space, determine the temperature vegetation drought index at any pixel.

Specifically, in the embodiment of the present invention, the temperature vegetation drought index at a certain pixel can be determined by the following formula:

In the formula, TDVI is the temperature vegetation drought index at a certain pixel, a and b are the constants and slopes of the fitting equation of the main line in the feature space, respectively, LST is the thermal infrared temperature data at the pixel, and LST _min is The thermal infrared temperature data corresponding to the wet line in the feature space, and FVC is the vegetation coverage data at this pixel.

In the embodiment of the present invention, by refining the feature space constructed by the thermal infrared temperature data and the vegetation coverage data, the accuracy of the temperature vegetation drought index can be improved.

On the basis of the above-mentioned embodiments, the agricultural irrigation area extraction method provided in the embodiment of the present invention is to extract the agricultural irrigation area in the crop planting area of the current year based on the temperature vegetation drought index, which specifically includes:

For any pixel position in the thermal infrared remote sensing image data of the current year, based on the temperature and vegetation of the pixels of the thermal infrared remote sensing image of each scene in the current year thermal infrared remote sensing image at the arbitrary pixel position in the thermal infrared remote sensing image data of the current year Drought index, to determine the irrigated area identification index;

Based on the irrigated area identification index, it is judged whether any pixel position in the thermal infrared remote sensing image data of the current year is an irrigated area pixel position, and based on the judgment result, the agricultural irrigation area in the crop planting area of the current year is determined. extract.

Specifically, in the embodiment of the present invention, when using the temperature vegetation drought index to extract the agricultural irrigation area in the crop planting area of the current year, for any pixel position in the thermal infrared remote sensing image data of the current year, it can first be based on the current year The temperature and vegetation drought index of the pixel of thermal infrared remote sensing image in the current year of each scene at any pixel position in the thermal infrared remote sensing image data is used to determine the irrigated area identification index. For example, the irrigated area identification index can be determined by the following formula:

Among them, IAI is the irrigated area identification index, TVDI _i is the temperature and vegetation drought index of the pixel of the i-th scene at any pixel position in the current year thermal infrared remote sensing image, n1 is the number of scenes in the current year thermal infrared remote sensing image data, C ₀ is the preset threshold.

Then, according to the irrigated area identification index, it is judged whether any pixel position in the thermal infrared remote sensing image data of the current year is the irrigated area pixel position, that is, it can be judged whether the irrigated area identification index is greater than or equal to the first threshold, and the first threshold can be set as required. Definitely, it is not specifically limited in this embodiment of the present invention.

Finally, according to the judgment result, the agricultural irrigation area in the crop planting area of the current year is extracted. That is, if the judgment result is that the identification index of the irrigation area is greater than or equal to the first threshold, then any pixel position can be determined to be the pixel position of the irrigation area, and the area formed by all the pixel positions of the irrigation area in the thermal infrared remote sensing image data of the current year can be determined as For agricultural irrigation. Otherwise, if the judging result is that the irrigated area identification index is less than the first threshold, it can be determined that any pixel position is not the irrigated area pixel position.

In the embodiment of the present invention, in the process of extracting the agricultural irrigation areas in the crop planting area of the current year, the irrigation area identification index is introduced to make the extraction results more reliable.

On the basis of the above-mentioned embodiments, the agricultural irrigation area extraction method provided in the embodiment of the present invention, the acquisition of the historical crop planting area in the monitoring area of the agricultural planting area specifically includes:

Obtaining the 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;

Determining the first-type normalized difference vegetation index at each pixel in each scene remote sensing image in the remote sensing image data of each historical year, and based on each pixel in each scene remote sensing image in the remote sensing image data of each historical year The first type of normalized difference vegetation index at , to determine the historical crop planting area.

Specifically, in the embodiment of the present invention, when obtaining the historical crop planting area in the monitoring area of the agricultural planting area, the historical remote sensing image data of many years close to the current year can be obtained first, that is, the remote sensing image data of multiple historical years before the current year . Furthermore, the first type of normalized difference vegetation index at each pixel in each scene remote sensing image in the historical remote sensing image data can be determined. The remote sensing image data of each historical year may include n1 scenes of remote sensing images, and each scene of remote sensing images may include n2 pixels. In each remote sensing image, the first-type normalized difference vegetation index at each pixel can be determined through the reflectance data at the pixel, which is not specifically limited in this embodiment of the present invention.

Then, according to the first-class normalized difference vegetation index at each pixel in the remote sensing image data of each historical year, the crop planting area corresponding to the remote sensing image data of each historical year can be determined first, and then Based on the crop planting area corresponding to the remote sensing image data of each historical year, the historical crop planting area in the monitoring area of the agricultural planting area can be determined. The historical crop planting area is the crop planting area before the current year, and may be a union of the crop planting areas corresponding to the remote sensing image data of each historical year.

In the embodiment of the present invention, historical crop planting areas are obtained by using historical remote sensing image data close to the current year for many years, which can make the result more accurate.

On the basis of the above-mentioned embodiments, the method for extracting agricultural irrigation areas provided in the embodiments of the present invention is based on the first type of normalized vegetation at each pixel in each scene remote sensing image in the remote sensing image data of each historical year Index, identifying the historical crop growing areas, including:

For the remote sensing image data of any historical year, determine the average value of the first normalized difference vegetation index at each pixel position in the remote sensing image data of any historical year;

For any remote sensing image of any scene in the remote sensing image data of any historical year, based on the first type of normalized normalized vegetation index at each pixel in the remote sensing image of any scene and the first type of normalized vegetation The index mean value determines the crop identification index at each pixel in the remote sensing image of any scene;

Based on the crop recognition index at each pixel in the remote sensing image of any scene, judge whether each pixel in the remote sensing image of any scene is a crop planting pixel, and based on the judgment result, determine the remote sensing image of any scene The corresponding crop growing area;

Based on the crop planting area corresponding to each remote sensing image in the remote sensing image data of any historical year, determine the crop planting area corresponding to the remote sensing image data of any historical year;

The historical crop planting area is determined 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, determine the first type of normalized vegetation at each pixel position in the remote sensing image data of any historical year Exponential mean, that is, for any pixel position in the remote sensing image data of any historical year, the normalization of the pixels of each remote sensing image in the remote sensing image data of any historical year at any pixel position The average value of the vegetation index can be expressed as

确定第j景遥感影像中第k个像元处的作物识别指数。即有：For the j-th remote sensing image in the remote sensing image data of the i-th historical year, based on the first-type normalized difference vegetation index NDVI _i and the first-type normalized difference vegetation index at the k-th pixel in the j-th remote sensing image average

Determine the crop identification index at the kth pixel in the jth remote sensing image. That is:

Among them, CRI _k is the crop identification index at the kth pixel in the remote sensing image of the jth scene.

Then, according to the crop recognition index at each pixel in the remote sensing image of scene j, it is judged whether each pixel in the remote sensing image of scene j is a crop planting pixel. That is, it is judged whether the CRI _k is greater than or equal to a second threshold, and the second threshold can be set according to requirements, which is not specifically limited in this embodiment of the present invention.

If the CRI _k is greater than or equal to the second threshold, it can be determined that the kth pixel is a crop planting pixel, otherwise it is indicated.

Finally, according to the judgment result, the crop planting area corresponding to the remote sensing image of scene j can be determined, that is, the number of crop planting pixels in the remote sensing image of scene j can be determined, and the crop planting area can be represented by this number.

According to the crop planting area corresponding to each remote sensing image in the remote sensing image data of the i-th historical year, the crop planting area corresponding to the remote sensing image data of the i-th historical year can be determined, that is, all the remote sensing image data in the i-th historical year The union of crop planting areas corresponding to landscape remote sensing images.

According to the union of the crop planting areas corresponding to the remote sensing image data of each historical year, the historical crop planting area can be determined. That is:

ZWI＝A ₁ ∪A ₂ ∪...∪A _n3

Among them, ZWI is the historical crop planting area; A ₁ , A ₂ , ..., A _n3 are the crop planting areas corresponding to the remote sensing image data of each historical year, represented by the number of crop planting pixels, and n3 is the number of historical years .

In the embodiment of the present invention, when determining the historical crop planting area in the monitoring area of the agricultural planting area, the crop planting area corresponding to the remote sensing image data of each historical year is used to determine the historical crop planting area, which can ensure the feasibility of the scheme.

On the basis of the above-mentioned embodiments, the agricultural irrigation area extraction method provided in the embodiment of the present invention, based on the thermal infrared remote sensing image data of the historical crop planting area in the current year, determines the current crop area in the monitoring area of the agricultural planting area Annual crop planting area, including:

Determining the second normalized difference vegetation index at the start point and end point of each pixel position in the thermal infrared remote sensing image data at the monitoring time period;

Based on the second normalized normalized vegetation index of each pixel position in the thermal infrared remote sensing image data at the start point and the end point, the crop planting area of the current year in the monitoring area of the agricultural planting area is determined.

Specifically, in the embodiment of the present invention, when determining the crop planting area in the current year, the thermal infrared remote sensing image data of the current year may be determined first, that is, the thermal infrared remote sensing image data of the historical crop planting area in the current year. Then determine the second normalized difference vegetation index at the start point and end point of each pixel position in the current year's thermal infrared remote sensing image data at the monitoring time period. Similarly, the thermal infrared remote sensing image data of the current year includes n1 scenes of thermal infrared remote sensing images, and each scene of thermal infrared remote sensing images may include n2 pixels.

According to the second normalized difference vegetation index of each pixel position in the thermal infrared remote sensing image data at the start point of the monitoring time period and the end point of the monitoring time period, the crop planting area of the current year in the agricultural planting area monitoring area is determined, That is, the cropping area of the current year within the historical cropping area. In the embodiment of the present invention, the crop planting area of the current year can be determined by the second-type normalized normalized vegetation index at the start point of the monitoring time period and the end point of the monitoring time period, combined with the length of the monitoring time period. This is not specifically limited in the example.

In the embodiment of the present invention, when determining the crop planting area of the current year, the historical crop planting area in the monitoring area of the agricultural planting area is determined first, and then the crop planting area of the current year in the historical crop planting area is determined. This can improve the efficiency of determining the crop planting area for the current year.

On the basis of the above-mentioned embodiments, the agricultural irrigation area extraction method provided in the embodiment of the present invention, the second classification based on the position of each pixel in the thermal infrared remote sensing image data at the starting point and the ending point A vegetation index is used to determine the crop planting area of the current year in the monitoring area of the agricultural planting area, specifically including:

For any pixel position in the thermal infrared remote sensing image data, based on the length of the monitoring time period, the second normalized vegetation of the any pixel position at the start point and the end point Index, determine the crop planting index of the current year at any pixel position;

The crop planting area of the current year is determined based on the crop planting index of the current year 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, it can first be based on the duration of the monitoring period of the irrigation area, the location of any pixel position The second normalized difference vegetation index at the start point of the irrigation area monitoring period and the end point of the irrigation area monitoring period determines the crop planting index of the current year at any pixel position. That is:

Among them, DZWI is the crop planting index of the current year at any pixel position, NDVI _tb is the second normalized difference vegetation index at the end point tb of the irrigation area monitoring period after crop planting at any pixel position, and NDVI _t0 is The second normalized difference vegetation index at the starting point t0 of the irrigation area monitoring period after crop planting at any pixel position. t _b -t ₀ is the monitoring period of the irrigation area.

Then, according to the crop planting index of the current year at each pixel position in the historical crop planting area, the crop planting area of the current year is determined. For any pixel position, it can be judged whether the crop planting index of the current year at any pixel position is greater than or equal to the third threshold, and the third threshold can be set according to needs, which is not specifically limited in the embodiments of the present invention .

If the crop planting index of the current year at any pixel position is greater than or equal to the third threshold, it can be determined that any pixel position is a crop pixel position, and all crop pixel positions constitute the crop planting area of the current year.

In the embodiment of the present invention, the crop planting index of the current year is determined through the duration of the monitoring time period, the second normalized vegetation index of each pixel position at the starting point of the monitoring time period and the end point of the monitoring time period, and then determined The crop planting area of the current year can ensure the accuracy of the crop planting area of the current year, and only consider the second normalized difference vegetation index at the starting point and the end point of the monitoring time period of each pixel position, which can be Improve determination efficiency.

On the basis of the above-mentioned embodiments, the method for extracting agricultural irrigation areas provided in the embodiments of the present invention is to extract the agricultural irrigation areas in the crop planting area of the current year based on the temperature vegetation drought index, and then further include:

Carry out vectorization processing on the crop planting area of the current year to obtain several farmland patches, and determine the agricultural irrigation area corresponding to each farmland patch;

For any patch of farmland, if the relative area deviation between the patch of any farmland and the agricultural irrigation area corresponding to the patch of any farmland is less than or equal to a deviation threshold, then the patch of any farmland is regarded as an irrigation area of agriculture.

Specifically, in the embodiment of the present invention, after extracting the agricultural irrigation area in the crop planting area of the current year, the crop planting area of the current year can be further vectorized. Since the crop planting area of the current year is a farmland area, the vector Several farmland patches can be obtained after chemical treatment. The grid size of each farmland plate is a×a, where a is the grid side length of the farmland plate. Combined with the extracted agricultural tank area, the agricultural irrigation area corresponding to each farmland patch can be determined.

Then, for any farmland patch, if the relative area deviation between the any farmland patch and its corresponding agricultural irrigation area is less than or equal to the deviation threshold, the any farmland patch is regarded as an agricultural irrigation area. The relative area deviation can be expressed as the ratio of the difference to the area of the farmland patch, that is:

Among them, GPI is the relative area deviation, S _0j is the area of the jth farmland patch, which can be obtained through statistics, and S _j is the area of the agricultural tank farm corresponding to the jth farmland patch.

If the relative area deviation GPI between the jth farmland patch and its corresponding agricultural irrigation area is less than or equal to the deviation threshold, the jth farmland patch can be considered as an agricultural irrigation area, that is, the agricultural tank area will be expanded into a farmland patch, which can expand the agricultural area. The area of the tank farm.

On the basis of the above-mentioned embodiments, the agricultural irrigation area extraction method provided in the embodiment of the present invention, the acquisition of the historical crop planting area in the monitoring area of the agricultural planting area also includes:

Obtain the precipitation information of each sub-region in the target area in the monitoring time period of many years including the current year;

Based on the precipitation information, the monitoring area of the agricultural planting area in the target area is determined.

Specifically, in the embodiment of the present invention, when determining the monitoring area of the agricultural planting area, the target area is first obtained, and then the target area is divided into multiple sub-areas, and the historical multi-year monitoring time of each sub-area including the current year is determined. The precipitation information of each sub-region is the precipitation information of the monitoring time period in each historical year of each sub-region.

Then, according to the precipitation information in each historical year of each sub-region, the monitoring area of the agricultural planting area in the target area is determined. For any sub-region, the precipitation anomaly percentage of any sub-region in each historical year can be determined according to the precipitation information of any sub-region in each historical year of the monitoring period, that is:

为任一子区域在所有历史年份的监测时间段的降水量信息的平均值。Among them, PJI is the precipitation anomaly percentage of any sub-region in each historical year, P is the precipitation information of any sub-region in the monitoring period of each historical year,

is the average value of the precipitation information of any sub-region in the monitoring period of all historical years.

Then, according to the precipitation anomaly percentage of any sub-region in each historical year, and the precipitation information of any sub-region during the monitoring time period in each historical year, the monitoring area of the agricultural planting area in the target area is determined. That is, for any sub-region, judge whether the precipitation anomaly percentage of any sub-region in each year is less than or equal to the fourth threshold, and whether the precipitation information of any sub-region in the monitoring time period in each historical year is less than the fourth threshold Five thresholds. Wherein, the fourth threshold may be set as required, which is not specifically limited in this embodiment of the present invention. The fifth threshold can be set as required, for example, it can be set to 800mm.

If there is a target year, the precipitation anomaly percentage of any sub-region in the target year is less than or equal to the fourth threshold, and the precipitation information of any sub-region in the monitoring period of the irrigated area in the target year is less than the fifth threshold, then determine Any of the sub-regions is the monitoring area of the candidate agricultural planting area, and the union of all the monitoring areas of the candidate agricultural planting area is the monitoring area of the agricultural planting area.

In the embodiment of the present invention, combined with meteorological information such as precipitation information, the determination of the monitoring area of the agricultural planting area can be made more quickly and accurately, and the extraction efficiency of the agricultural irrigation area can be improved.

As shown in Figure 2, on the basis of the foregoing embodiments, an agricultural irrigation area extraction system is provided in an embodiment of the present invention, including:

The first acquisition module 21 is used to acquire the historical crop planting area in the agricultural planting area monitoring area, and determine the current crop planting area in the agricultural planting area monitoring area based on the thermal infrared remote sensing image data of the historical crop planting area in the current year. annual crop planting area;

The second acquisition module 22 is used to determine the time-series thermal infrared temperature data and vegetation coverage data within the monitoring period based on the thermal infrared remote sensing image data, and perform downscaling processing on the thermal infrared temperature data to obtain detailed Thermal infrared temperature data;

The extraction module 23 is used to determine 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 based on the temperature vegetation drought index, The agricultural irrigation areas in the crop planting area of the current year are extracted.

On the basis of the above-mentioned embodiments, in the agricultural irrigation area extraction system provided in the embodiments of the present invention, the extraction module is specifically used for:

For any pixel of any scene in the thermal infrared remote sensing image data of the current year in the current year thermal infrared remote sensing image, based on the thermal infrared temperature data at any pixel, the vegetation coverage data and the features Wet edge temperature in the space, determine the temperature vegetation drought index at any pixel.

On the basis of the above-mentioned embodiments, the agricultural irrigation area extraction system provided in the embodiments of the present invention, the extraction module is also used for:

For any pixel position in the thermal infrared remote sensing image data of the current year, based on the temperature and vegetation of the pixels of the thermal infrared remote sensing image of each scene in the current year thermal infrared remote sensing image at the arbitrary pixel position in the thermal infrared remote sensing image data of the current year Drought index, to determine the irrigated area identification index;

Based on the irrigated area identification index, it is judged whether any pixel position in the thermal infrared remote sensing image data of the current year is an irrigated area pixel position, and based on the judgment result, the agricultural irrigation area in the crop planting area of the current year is determined. extract.

On the basis of the above embodiments, the agricultural irrigation area extraction system provided in the embodiments of the present invention, the first acquisition module is specifically used for:

Obtaining the 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;

Determining the first-type normalized difference vegetation index at each pixel in each scene remote sensing image in the remote sensing image data of each historical year, and based on each pixel in each scene remote sensing image in the remote sensing image data of each historical year The first type of normalized difference vegetation index at , to determine the historical crop planting area.

On the basis of the above embodiments, the agricultural irrigation area extraction system provided in the embodiments of the present invention, the first acquisition module is also specifically used for:

For the remote sensing image data of any historical year, determine the average value of the first normalized difference vegetation index at each pixel position in the remote sensing image data of any historical year;

For any remote sensing image of any scene in the remote sensing image data of any historical year, based on the first type of normalized normalized vegetation index at each pixel in the remote sensing image of any scene and the first type of normalized vegetation The index mean value determines the crop identification index at each pixel in the remote sensing image of any scene;

Based on the crop recognition index at each pixel in the remote sensing image of any scene, judge whether each pixel in the remote sensing image of any scene is a crop planting pixel, and based on the judgment result, determine the remote sensing image of any scene The corresponding crop growing area;

Based on the crop planting area corresponding to each remote sensing image in the remote sensing image data of any historical year, determine the crop planting area corresponding to the remote sensing image data of any historical year;

The historical crop planting area is determined based on the crop planting area corresponding to the remote sensing image data of each historical year.

On the basis of the above embodiments, the agricultural irrigation area extraction system provided in the embodiments of the present invention, the first acquisition module is also specifically used for:

For the remote sensing image data of any historical year, determine the average value of the first normalized difference vegetation index at each pixel position in the remote sensing image data of any historical year;

For any remote sensing image of any scene in the remote sensing image data of any historical year, based on the first type of normalized normalized vegetation index at each pixel in the remote sensing image of any scene and the first type of normalized vegetation The index mean value determines the crop identification index at each pixel in the remote sensing image of any scene;

Based on the crop recognition index at each pixel in the remote sensing image of any scene, judge whether each pixel in the remote sensing image of any scene is a crop planting pixel, and based on the judgment result, determine the remote sensing image of any scene The corresponding crop growing area;

Based on the crop planting area corresponding to each remote sensing image in the remote sensing image data of any historical year, determine the crop planting area corresponding to the remote sensing image data of any historical year;

The historical crop planting area is determined based on the crop planting area corresponding to the remote sensing image data of each historical year.

On the basis of the above embodiments, the agricultural irrigation area extraction system provided in the embodiments of the present invention, the first acquisition module is also specifically used for:

Determining the second normalized difference vegetation index at the start point and end point of each pixel position in the thermal infrared remote sensing image data at the monitoring time period;

Based on the second normalized normalized vegetation index of each pixel position in the thermal infrared remote sensing image data at the start point and the end point, the crop planting area of the current year in the monitoring area of the agricultural planting area is determined.

On the basis of the above embodiments, the agricultural irrigation area extraction system provided in the embodiments of the present invention, the first acquisition module is also specifically used for:

For any pixel position in the thermal infrared remote sensing image data, based on the length of the monitoring time period, the second normalized vegetation of the any pixel position at the start point and the end point Index, determine the crop planting index of the current year at any pixel position;

The crop planting area of the current year is determined based on the crop planting index of the current year at each pixel position in the historical crop planting area.

On the basis of the above embodiments, the agricultural irrigation area extraction system provided in the embodiments of the present invention also includes a correction module for:

Carry out vectorization processing on the crop planting area of the current year to obtain several farmland patches, and determine the agricultural irrigation area corresponding to each farmland patch;

For any patch of farmland, if the relative area deviation between the patch of any farmland and the agricultural irrigation area corresponding to the patch of any farmland is less than or equal to a deviation threshold, then the patch of any farmland is regarded as an irrigation area of agriculture.

On the basis of the above-mentioned embodiments, the agricultural irrigation area extraction system provided in the embodiments of the present invention also includes a third acquisition module, which is used for:

Obtain the precipitation information of each sub-region in the target area in the monitoring time period of many years including the current year;

Based on the precipitation information, the monitoring area of the agricultural planting area in the target area is determined.

Specifically, the function of each module in the agricultural irrigation area extraction system provided in the embodiment of the present invention is in one-to-one correspondence with the operation process of each step in the above-mentioned method embodiment, and the achieved effect is also consistent. Please refer to the above-mentioned embodiment for details. This will not be described in detail in the embodiments of the present invention.

FIG. 3 illustrates a schematic diagram of the physical structure of an electronic device. As shown in FIG. 3 , the electronic device may include: a processor (processor) 310, a communication interface (Communications 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 through the communication bus 340 . The processor 310 can call the logic instructions in the memory 330 to execute the agricultural irrigation area extraction method provided in the above-mentioned embodiments, the method includes: obtaining the historical crop planting area in the monitoring area of the agricultural planting area, and based on the historical crop planting area Based on the thermal infrared remote sensing image data of the region in the current year, determine the crop planting area of the current year in the agricultural planting area monitoring area; based on the thermal infrared remote sensing image data, determine the thermal infrared temperature data and vegetation of the time series within the monitoring period Coverage data, and perform downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data; based on the feature space constructed by the refined thermal infrared temperature data and the vegetation coverage data, determine the current The temperature vegetation drought index of the crop planting area of the year, and based on the temperature vegetation drought index, the agricultural irrigation area in the crop planting area of the current year is extracted.

In addition, the above-mentioned logic instructions in the memory 330 may be implemented in the form of software functional units and may be stored in a computer-readable storage medium when sold or used as an independent product. Based on this understanding, the essence of the technical solution of the present invention or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including Several instructions are used to make a computer device (which may be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disk or optical disk and other media that can store program codes. .

On the other hand, the present invention also provides a computer program product. The computer program product includes a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can Executing the method for extracting agricultural irrigation areas provided in each of the above embodiments, the method includes: obtaining historical crop planting areas in the monitoring area of agricultural planting areas, and determining the selected area based on the thermal infrared remote sensing image data of the historical crop planting areas in the current year. The crop planting area of the current year in the agricultural planting area monitoring area; based on the thermal infrared remote sensing image data, determine the thermal infrared temperature data and vegetation coverage data of the time series within the monitoring period, and perform the thermal infrared temperature data Downscaling processing to obtain refined thermal infrared temperature data; based on the feature space constructed by the refined thermal infrared temperature data and the vegetation coverage data, determine the temperature vegetation drought index of the crop planting area in the current year, and based on the The temperature vegetation drought index is used to extract the agricultural irrigation area in the crop planting area of the current year.

In yet another aspect, the present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, it is implemented to perform the agricultural irrigation area extraction method provided in the above-mentioned embodiments, the The method includes: obtaining the historical crop planting area in the monitoring area of the agricultural planting area, and determining the crop planting area of the current year 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; Based on the thermal infrared remote sensing image data, determine the thermal infrared temperature data and vegetation coverage data of the time series within the monitoring period, and perform downscaling processing on the thermal infrared temperature data to obtain refined thermal infrared temperature data; According to the feature space constructed by the refined thermal infrared temperature data and the vegetation coverage data, the temperature vegetation drought index of the crop planting area of the current year is determined, and based on the temperature vegetation drought index, the crop planting area of the current year is calculated Extraction in agricultural irrigation areas within.

The device embodiments described above are only illustrative, and the units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in One place, or it can be distributed to multiple network elements. Part or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. It can be understood and implemented by those skilled in the art without any creative effort.

Through the above description of the implementations, those skilled in the art can clearly understand that each implementation can be implemented by means of software plus a necessary general-purpose hardware platform, and of course also by hardware. Based on this understanding, the essence of the above technical solution or the part that contributes to the prior art can be embodied in the form of software products, and the computer software products can be stored in computer-readable storage media, such as ROM/RAM, magnetic discs, optical discs, etc., including several instructions to make a computer device (which may be a personal computer, server, or network device, etc.) execute the methods described in various embodiments or some parts of the embodiments.

Finally, it should be noted that: the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: it can still be Modifications are made to the technical solutions described in the foregoing embodiments, or equivalent replacements are made to some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the various embodiments of the present invention.

Claims (10)

1. An agricultural irrigation area extraction method, characterized by comprising the following steps:

acquiring a historical crop planting area in an agricultural planting area monitoring area, and determining a current year crop planting area in the agricultural planting area monitoring area based on thermal infrared remote sensing image data of the historical crop planting area in the current year;

based on the thermal infrared remote sensing image data, determining thermal infrared temperature data and vegetation coverage data of a time sequence in a monitoring time period, and performing downscaling treatment 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 characteristic space constructed by the refined thermal infrared temperature data and the vegetation coverage data, and extracting an agricultural irrigation area in the current year crop planting area based on the temperature vegetation drought index.

2. The agricultural irrigation area extraction method 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 the thermal infrared remote sensing image of any scene in the thermal infrared remote sensing image data of the current year, 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 feature space.

3. The method for extracting an agricultural irrigation area according to claim 1, wherein the extracting the agricultural irrigation area in the current year crop planting area based on the temperature vegetation drought index specifically comprises:

determining an irrigation area identification index for any pixel position in the current year thermal infrared remote sensing image data based on a temperature vegetation drought index of a pixel of each scene of the current year thermal infrared remote sensing image at any pixel position in the current year thermal infrared remote sensing image data;

and judging whether any pixel position in the current year thermal infrared remote sensing image data is a irrigation area pixel position based on the irrigation area identification index, and extracting an agricultural irrigation area in the current year crop planting area based on a judgment result.

4. The method for extracting an agricultural irrigation area according to claim 1, wherein the step of obtaining a historical crop planting area in an agricultural planting area monitoring area comprises the following steps:

Acquiring historical years of remote sensing image data adjacent to the current year, and determining the remote sensing image data of each historical year in the historical years of remote sensing image data;

and determining a first type of normalized vegetation index at each pixel in each scene of the remote sensing image data of each historical year, and determining the historical crop planting area based on the first type of normalized vegetation index at each pixel in each scene of the remote sensing image data of each historical year.

5. The method according to claim 4, wherein the determining the historical crop planting area based on the first type of 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 one scene remote sensing image in the remote sensing image data of any historical year, determining a crop identification index at each pixel in the any one scene remote sensing image based on a first type normalized vegetation index at each pixel in the any one scene remote sensing image and the first type normalized vegetation index mean;

Judging whether each pixel in any one scene remote sensing image is a crop planting pixel or not based on a crop identification index at each pixel in the any one scene remote sensing image, and determining a crop planting area corresponding to the any one scene remote sensing image based on a judging 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 remote sensing image in the remote sensing image data of any historical year;

and determining the historical crop planting areas based on the crop planting areas corresponding to the remote sensing image data of each historical year.

6. The agricultural irrigation area extraction method according to claim 1, wherein the determining the current year crop planting area in the agricultural planting area monitoring 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 the starting point and the ending point of each pixel position in the thermal infrared remote sensing image data in the monitoring time period;

and determining a current year crop planting area in the agricultural planting area monitoring area based on second-class normalized vegetation indexes 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 an agricultural irrigation area according to claim 6, wherein the determining the current year crop planting area in the agricultural planting area monitoring area based on the second type of normalized vegetation indexes 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 current year crop planting index at the any pixel position based on the duration of the monitoring time period and a second type of normalized vegetation index of the any pixel position at the starting point and the ending point;

and determining the current year crop planting area based on the year crop planting index at each pixel position in the historical crop planting area.

8. The agricultural irrigation area extraction method according to any one of claims 1-7, wherein the extracting the agricultural irrigation area within the current year crop planting area based on the temperature vegetation drought index further comprises:

vectorizing the current year crop planting area to obtain a plurality of farmland patches, and determining an agricultural irrigation area corresponding to each farmland patch;

And 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 smaller than or equal to a deviation threshold value, taking the any farmland patch as the agricultural irrigation area.

9. The agricultural irrigation area extraction method according to any one of claims 1 to 7, wherein the acquiring a historical crop planting area within the agricultural planting area monitoring area further comprises, before:

acquiring precipitation information of each subarea in the target area in the monitoring time period for a plurality of years including the current year;

and determining the agricultural planting area monitoring area in the target area based on the precipitation amount information.

10. An agricultural irrigation area extraction system, comprising:

the first acquisition module is used for acquiring a historical crop planting area in the agricultural planting area monitoring area and determining a current year crop planting area in the agricultural planting area monitoring 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;

The extraction module is used for determining a temperature vegetation drought index of the current year crop planting area based on the characteristic space constructed by the refined thermal infrared temperature data and the vegetation coverage data, and extracting an agricultural irrigation area in the current year crop planting area based on the temperature vegetation drought index.

Images