CN113793376A - Irrigation water body extraction method - Google Patents

Irrigation water body extraction method Download PDF

Info

Publication number
CN113793376A
CN113793376A CN202111076430.8A CN202111076430A CN113793376A CN 113793376 A CN113793376 A CN 113793376A CN 202111076430 A CN202111076430 A CN 202111076430A CN 113793376 A CN113793376 A CN 113793376A
Authority
CN
China
Prior art keywords
water body
index
calculating
spatial
moran
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
CN202111076430.8A
Other languages
Chinese (zh)
Other versions
CN113793376B (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.)
Anhui University of Science and Technology
Original Assignee
Anhui University of Science and Technology
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 Anhui University of Science and Technology filed Critical Anhui University of Science and Technology
Priority to CN202111076430.8A priority Critical patent/CN113793376B/en
Publication of CN113793376A publication Critical patent/CN113793376A/en
Application granted granted Critical
Publication of CN113793376B publication Critical patent/CN113793376B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/243Classification techniques relating to the number of classes
    • G06F18/24323Tree-organised classifiers
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10032Satellite or aerial image; Remote sensing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30181Earth observation
    • G06T2207/30184Infrastructure

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geometry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Image Processing (AREA)

Abstract

The invention discloses an irrigation water body extraction method, which comprises the following steps: (1) selecting remote sensing image data according to the requirements of irrigation time and geographic position; (2) calculating the sensitive spectral index of the water body to obtain the spatial distribution of the sensitive spectral index; (3) performing spatial autocorrelation analysis on the spatial information of the water body sensitive spectral index; (4) calculating a global Moran index (Moran's I), and judging whether aggregation or abnormal values occur in the research area according to the value range of the I value; (5) if an aggregate or outlier occurs, then calculate the local Moran's I index; (6) calculating a Getis-Ord Gi index to obtain the water space aggregation distribution; (7) and reasonable remote sensing images are utilized to obtain the planting structure of the research area through ground object classification, and further the irrigation water body area or the irrigation water body areas of different crop types are finally obtained.

Description

Irrigation water body extraction method
Technical Field
The invention relates to the fields of water resource optimization, remote sensing image processing technology, spatial data processing and the like, in particular to a monitoring method for efficiently and accurately acquiring the water body area.
Background
According to FAO reports, the global population reaches 97 billion by 2050, and 108 billion by 2080, and the food demand is expected to increase greatly under the push of population factors. In order to increase agricultural yields, irrigation is one of the effective measures in about one third of the world's farmlands, particularly dry land areas, to moderate or high deterioration. With the development of socio-economic, the contradiction between agricultural water and industrial urban water is growing, and the importance of irrigation in land management, as well as the impact on regional and global climate, is gradually recognized and brought to the scientific community's attention.
The most important thing to monitor the irrigation area is the extraction of the ground water. At present, the extraction method of regional water mainly utilizes a water sensitivity index to obtain the water area of a region by setting a threshold, but the extraction of surface water by utilizing the water sensitivity index faces two main problems: the first is that the results of surface water bodies obtained by using different water body indexes are inconsistent and unreliable; the second is that the threshold is not fixed and may vary from environment to environment and location to location. Therefore, how to improve the extraction accuracy of the surface water body is one of the hot spots of the current research.
In summary, in order to solve the above problems, the present invention provides an extraction method for irrigation water.
Disclosure of Invention
The invention provides an irrigation water body extraction method, which can better solve the problems of difficulty in threshold value determination and large difference of obtained water bodies due to different water body sensitivity indexes, further improve the irrigation area monitoring precision and have certain guiding significance for reasonable water quantity preparation and high-efficiency water use of cultivated land. The method comprises the following steps:
the invention adopts the following technical scheme for realizing the purpose:
an irrigation water body extraction method utilizes multi-temporal multi-source remote sensing images as data sources to calculate water body sensitivity indexes such as common normalized water difference index (NDWI), improved normalized water body index (MNDWI), automatic extraction index (AWEI) and the like, but not limited to the enumerated water body sensitivity indexes, after spatial distribution of water body sensitivity spectral indexes is calculated, spatial distribution obtained by the water body sensitivity indexes is processed by utilizing a spatial autocorrelation analysis method, water body gathering spatial distribution is extracted according to z score and p value in a Getis-Ord Gi index, the problems that threshold value determination is difficult, extraction of water body difference is large due to different water body sensitivity indexes are solved, and irrigation area extraction accuracy can be further improved.
An extraction method of irrigation water body comprises the following steps:
step 1, selecting remote sensing image data according to irrigation time and geographical position requirements;
step 2, calculating the water sensitive spectral index to obtain the spatial distribution of the water sensitive spectral index;
step 3, carrying out spatial autocorrelation analysis on the spatial information of the water body sensitive spectral index;
step 4, calculating a global Moran index (Moran's I), and judging whether aggregation or abnormal values occur in the research area according to the value range of the value I;
step 5, if the aggregation or abnormal value occurs, calculating a local Moran's I index;
step 6, calculating Getis-Ord GiIndexing to obtain water space aggregation distribution;
step 7, obtaining a planting structure of a research area by classifying ground objects by using a reasonable remote sensing image, and further finally obtaining the area of an irrigation water body or the areas of irrigation water bodies of different crop types;
further, the step 1 comprises:
step 8, the original image needs to be preprocessed before being applied, wherein preprocessing links comprise 3 links such as geometric correction, radiation correction and atmospheric correction, and selective preprocessing can be performed according to the type and the grade of the downloaded remote sensing image;
step 9, preprocessing the remote sensing image by utilizing an existing module on the ENVI software or a third-party open source program;
further, the step 2 comprises:
step 10, the water body sensitive spectral indexes are more, MNDWI or AWEI and other water body sensitive spectral indexes are selected, but the water body sensitive spectral indexes are not limited to use;
step 11, calculating an improved normalized water body index (MNDWI) according to the following formula:
Figure BDA0003262428080000021
wherein ρ is the reflectivity of the spectral band of the remote sensing image: BIUE is a blue band (0.45-0.51 μm), GRE is a green band (0.53-0.59 μm), NIR is a near infrared band (0.85-0.88 μm), SWIR1 is a short wave infrared band (1.57-1.65 μm), and SWIR2 is a short wave infrared band (2.21-2.29 μm).
Step 12, calculate the auto extraction index (AWEI), including AWEInshAnd AWEIshTwo types, the formula of which is as follows:
AWEInsh=4×(ρGREENSWIR1)-(0.25×ρNIR+2.75×ρSWIR1)
AWEIsh=ρBIUE+2.5×ρGREEN-1.5×(ρNIRSWIR1)-0.25×ρSWIR2)
wherein the meanings of the items are the same as the meanings of the items in the MNDWI calculation formula in the step 11;
further, the step 4 comprises:
step 13, before calculating the global Morland index, vectorizing the water sensitive spectrum index spatial information data, wherein the vectorizing can be processed in ArcGIS, and the file is converted into a file with a format of shp;
step 14, calculate the global Moran index (Moran's I), which is given by the formula:
Figure BDA0003262428080000031
Figure BDA0003262428080000032
where N is the number of data, WijIs a spatial weight, XiAnd XjThe property values of the spatial object at the ith and jth positions respectively,
Figure BDA0003262428080000033
is the average value of X; z is a normalized statistical quantity, E (I) is an expected value of the autocorrelation of the observed variable, var (I) is a squareAnd (4) poor.
The value range of the Moran's I index is [ -1,1], and I >0 represents that the attribute values present spatial positive correlation and tend to spatial aggregation characteristics; i <0 represents that the attribute value presents space negative correlation and tends to space dispersion characteristics; i-0 represents that the attribute values tend to be spatially randomly distributed.
Step 15, the global Morland index calculation can be processed in ArcGIS or GeoDa related modules;
further, the step 5 comprises:
step 16, judging according to the Moran's I index I value, if no abnormal value occurs, the method is invalid, no water body appears in a possibly researched area, and if aggregation or an abnormal value occurs, analyzing by using local space autocorrelation;
step 17, calculating a local Moran's I index, wherein the formula is as follows:
Figure BDA0003262428080000034
Figure BDA0003262428080000035
wherein each meaning is the same as each meaning in the calculation formula of the global Moran index in the step 4;
further, the step 6 comprises:
step 18, calculating Getis-Ord GiThe formula of the index is as follows:
Figure BDA0003262428080000041
wherein xjIs the attribute value of element j, wi,jIs the spatial weight between elements i and j, n is the total number of elements, and:
Figure BDA0003262428080000042
Figure BDA0003262428080000043
step 19, calculating Getis-Ord GiDuring indexing, the weight needs to be determined, and a fixed distance model, a surface adjacency model (shared edge and intersection), a surface adjacency model (shared edge, intersection and adjacent point), an inverse distance model, an inverse distance square model, a K nearest neighbor model and the like can be selected for construction according to actual conditions when a weight matrix file of a space is constructed;
step 20, calculating a z score and a p value by using relevant modules of ArcGIS or GeoDa software and the like;
step 21, selecting the spatial distribution with the Z score larger than 1.96 as the spatial aggregation distribution of the water body, wherein the confidence coefficient is 95%;
step 22, obtaining the spatial distribution of the water body by using a classification method such as a decision tree;
step 23, calculating the area of the water body of the region by using the spatial resolution of the pixels and the number of the pixels counted by the water body;
further, the step 7 includes:
step 24, selecting a remote sensing image with a proper time phase according to the characteristics of crops, calculating a normalized vegetation index (NDVI), classifying the research area according to the numerical value of the NDVI, and obtaining the spatial distribution of cultivated land and the spatial distribution of a crop planting structure of the research area;
step 25, registering spatial distribution of cultivated land and spatial distribution of water body in the research area, and obtaining spatial distribution of irrigation water body of cultivated land in the research area by using a mask method;
step 26, counting the number of pixels of the irrigation water body of the farmland, and calculating the area of the irrigation water body;
step 27, carrying out registration treatment on the spatial distribution of the crop planting structure and the spatial distribution of the water body in the research area, and obtaining the spatial distribution of the irrigation water body of different crop species by using a mask method;
and step 28, counting the pixel numbers of the farmland irrigation water bodies of different crop types, and calculating the irrigation water body area sizes of different crop types.
The embodiment of the invention has the following beneficial effects:
in the scheme, when the water body sensitivity index is used for extracting the water body in the research area, the uncertainty of threshold selection often causes larger difference of extraction areas, the water body extraction precision can be further improved by using the scheme, the subjectivity of threshold selection is avoided, and the working intensity of ground monitoring is reduced; in the scheme, the remote sensing image data of different time phases are utilized, the area and the spatial distribution of the irrigation water body in different time periods can be monitored, and a basis is provided for perfecting the irrigation system of the research area.
Drawings
FIG. 1 is a schematic view of the operational flow of an irrigation water extraction method of the present invention;
FIG. 2 is a schematic view of the monitoring results of an embodiment of the method for extracting irrigation water according to the present invention.
Detailed Description
The invention is further illustrated by the following specific examples. The following examples are intended to illustrate the invention but are not intended to limit the scope of the invention.
The embodiment of the invention provides an irrigation water body extraction method, which can quickly, effectively and accurately monitor the irrigation water body area of a research area and has certain guiding significance for further improving the regional water utilization efficiency and perfecting an irrigation system.
As shown in fig. 1, an embodiment of the present invention provides an irrigation water body extraction method, including the following main steps:
step 1, selecting remote sensing image data according to time and geographic position requirements;
step 2, calculating the water sensitive spectral index to obtain the spatial distribution of the water sensitive spectral index;
step 3, carrying out spatial autocorrelation analysis on the spatial information of the water body sensitive spectral index;
step 4, calculating a global Moran index (Moran's I), and judging whether aggregation or abnormal values occur in the research area according to the value range of the value I;
step 5, if aggregation or outliers occur, calculate the local Moran's I index:
step 6, calculating Getis-Ord GiIndexing to obtain water space aggregation distribution;
step 7, obtaining a planting structure of a research area by classifying ground objects by using a reasonable remote sensing image, and further finally obtaining the area of an irrigation water body or the areas of irrigation water bodies of different crop types;
in the invention, the selection of the water body sensitive spectral index suggests to select the MNDWI or AWEI and other sensitive spectral indexes of sensitive water bodies, and other sensitive spectral indexes of sensitive water bodies can also be selected.
In the invention, when the Getis-Ord Gi index is calculated, a fixed distance model, a surface adjacency model (shared edge and intersection), a surface adjacency model (shared edge, intersection and adjacency point), an inverse distance model, an inverse distance square model, a K nearest neighbor model and the like are selected according to actual conditions to construct a weight matrix.
The basic features, technical solutions and advantageous effects of the present invention have been described above. It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein.

Claims (7)

1. An extraction method for an irrigation water body, comprising:
step 1, selecting remote sensing image data according to irrigation time and geographical position requirements;
step 2, calculating the water sensitive spectral index to obtain the spatial distribution of the water sensitive spectral index;
step 3, carrying out spatial autocorrelation analysis on the spatial information of the water body sensitive spectral index;
step 4, calculating a global Moran index (Moran' sI), and judging whether an aggregation or an abnormal value occurs in a research area according to the value range of the I value;
step 5, if an aggregation or abnormal value occurs, calculating a local Moran' sI index;
step 6, calculating Getis-OrdGiIndexing to obtain water space aggregation distribution;
and 7, obtaining the planting structure of the research area by utilizing the reasonable remote sensing image through ground object classification, and further finally obtaining the irrigation water body area or the irrigation water body areas of different crop types.
2. The method of claim 1, wherein the step 1 comprises:
step 8, the original image needs to be preprocessed before being applied, wherein preprocessing links comprise 3 links such as geometric correction, radiation correction and atmospheric correction, and selective preprocessing can be performed according to the type and the grade of the downloaded remote sensing image;
and 9, preprocessing the remote sensing image by utilizing an existing module on the ENVI software or a third-party open source program.
3. The method of claim 1, wherein the step 2 comprises:
step 10, the water body sensitive spectral indexes are more, MNDWI or AWEI and other water body sensitive spectral indexes are selected, but the water body sensitive spectral indexes are not limited to use;
step 11, calculating an improved normalized water body index (MNDWI) according to the following formula:
Figure FDA0003262428070000011
wherein ρ is the reflectivity of the spectral band of the remote sensing image: BIUE is a blue band (0.45-0.51 μm), GRE is a green band (0.53-0.59 μm), NIR is a near infrared band (0.85-0.88 μm), SWIR1 is a short wave infrared band (1.57-1.65 μm), and SWIR2 is a short wave infrared band (2.21-2.29 μm).
Step 12, calculate the auto extraction index (AWEI), including AWEInshAnd AWEIshTwo types, the formula of which is as follows:
AWEInsh=4×(ρGREENSWIR1)-(0.25×ρNIR+2.75×ρSWIR1)
AWEIsh=ρBIUE+2.5×ρGREEN-1.5×(ρNIRSWIR1)-0.25×ρSWIR2)
the meaning of each item is the same as that of the calculation formula of MNDWI in step 11.
4. The method of claim 1, wherein the step 4 comprises:
step 13, before calculating the global Morland index, vectorizing the water sensitive spectrum index spatial information data, wherein the vectorizing can be processed in ArcGIS, and the file is converted into a file with a format of shp;
step 14, calculate the global Moran index (Moran' sI), which is expressed as follows:
Figure FDA0003262428070000021
Figure FDA0003262428070000022
where N is the number of data, WijIs a spatial weight, XiAnd XjThe property values of the spatial object at the ith and jth positions respectively,
Figure FDA0003262428070000023
is the average value of X; z is the normalized statistical value, E (I) is the expected value of the autocorrelation of the observed variable, and var (I) is the variance.
The value range of the Moran' sI index is [ -1,1], and I >0 represents that the attribute value presents spatial positive correlation and tends to spatial aggregation characteristics; i <0 represents that the attribute value presents space negative correlation and tends to space dispersion characteristics; i-0 represents that the attribute values tend to be spatially randomly distributed.
The global Morland index calculation can be processed in ArcGIS or GeoDa related modules, step 15.
5. The method of claim 1, wherein the step 5 comprises:
step 16, judging according to the Moran' sI index I value, if no abnormal value occurs, the method is invalid, no water body appears in a possibly researched area, and if aggregation or abnormal value occurs, analyzing by using local space autocorrelation;
step 17, calculating a local Moran' sI index, wherein the formula is as follows:
Figure FDA0003262428070000024
Figure FDA0003262428070000031
wherein the meanings of all the items are the same as the meanings of all the items in the calculation formula of the global Moran index in the step 4.
6. The method of claim 1, wherein the step 6 comprises:
step 18, calculating Getis-OrdGiThe formula of the index is as follows:
Figure FDA0003262428070000032
wherein xjIs the attribute value of element j, wi,jIs the spatial weight between elements i and j, n is the total number of elements, and:
Figure FDA0003262428070000033
Figure FDA0003262428070000034
step 19, in calculating Getis-OrdGiDuring indexing, the weight needs to be determined, and a fixed distance model, a surface adjacency model (shared edge and intersection), a surface adjacency model (shared edge, intersection and adjacent point), an inverse distance model, an inverse distance square model, a K nearest neighbor model and the like can be selected for construction according to actual conditions when a weight matrix file of a space is constructed;
step 20, calculating a z score and a p value by using relevant modules of ArcGIS or GeoDa software and the like;
step 21, selecting the spatial distribution with the Z score larger than 1.96 as the spatial aggregation distribution of the water body, wherein the confidence coefficient is 95%;
step 22, obtaining the spatial distribution of the water body by using a classification method such as a decision tree;
and step 23, calculating the area of the water body of the region by using the spatial resolution of the pixels and the number of the pixels counted by the water body.
7. The method of claim 1, wherein the step 7 comprises:
step 24, selecting a remote sensing image with a proper time phase according to the characteristics of crops, calculating a normalized vegetation index (NDVI), classifying the research area according to the numerical value of the NDVI, and obtaining the spatial distribution of cultivated land and the spatial distribution of a crop planting structure of the research area;
step 25, registering spatial distribution of cultivated land and spatial distribution of water body in the research area, and obtaining spatial distribution of irrigation water body of cultivated land in the research area by using a mask method;
step 26, counting the number of pixels of the irrigation water body of the farmland, and calculating the area of the irrigation water body;
step 27, carrying out registration treatment on the spatial distribution of the crop planting structure and the spatial distribution of the water body in the research area, and obtaining the spatial distribution of the irrigation water body of different crop species by using a mask method;
and step 28, counting the pixel numbers of the farmland irrigation water bodies of different crop types, and calculating the irrigation water body area sizes of different crop types.
CN202111076430.8A 2021-09-14 2021-09-14 Irrigation water body extraction method Active CN113793376B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111076430.8A CN113793376B (en) 2021-09-14 2021-09-14 Irrigation water body extraction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111076430.8A CN113793376B (en) 2021-09-14 2021-09-14 Irrigation water body extraction method

Publications (2)

Publication Number Publication Date
CN113793376A true CN113793376A (en) 2021-12-14
CN113793376B CN113793376B (en) 2023-12-05

Family

ID=78880317

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111076430.8A Active CN113793376B (en) 2021-09-14 2021-09-14 Irrigation water body extraction method

Country Status (1)

Country Link
CN (1) CN113793376B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2018101946A4 (en) * 2018-11-07 2019-01-17 Wuhan University Geographical multivariate flow data spatio-temporal autocorrelation analysis method based on cellular automaton
CN110288647A (en) * 2019-06-25 2019-09-27 中国水利水电科学研究院 One kind monitoring Irrigation Project Design area method based on high resolution satellite remote sensing image
CN110334880A (en) * 2019-07-16 2019-10-15 安徽理工大学 A kind of effective irrigation area monitoring method
CN111460003A (en) * 2020-04-04 2020-07-28 南京国准数据有限责任公司 Method for detecting coupling relation between land utilization and earth surface temperature based on urban grouping scale
CN112033937A (en) * 2020-09-08 2020-12-04 西安科技大学 Method for evaluating water body extraction precision
US20210199579A1 (en) * 2019-12-31 2021-07-01 Beijing Normal University Method and system for urban impervious surface extraction based on remote sensing

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2018101946A4 (en) * 2018-11-07 2019-01-17 Wuhan University Geographical multivariate flow data spatio-temporal autocorrelation analysis method based on cellular automaton
CN110288647A (en) * 2019-06-25 2019-09-27 中国水利水电科学研究院 One kind monitoring Irrigation Project Design area method based on high resolution satellite remote sensing image
CN110334880A (en) * 2019-07-16 2019-10-15 安徽理工大学 A kind of effective irrigation area monitoring method
US20210199579A1 (en) * 2019-12-31 2021-07-01 Beijing Normal University Method and system for urban impervious surface extraction based on remote sensing
CN111460003A (en) * 2020-04-04 2020-07-28 南京国准数据有限责任公司 Method for detecting coupling relation between land utilization and earth surface temperature based on urban grouping scale
CN112033937A (en) * 2020-09-08 2020-12-04 西安科技大学 Method for evaluating water body extraction precision

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
FATMA HAOUAS; ZOUHOUR BEN DHIAF; BASEL SOLAIMAN: "Fusion of spatial autocorrelation and spectral data for remote sensing image classification", 2016 2ND INTERNATIONAL CONFERENCE ON ADVANCED TECHNOLOGIES FOR SIGNAL AND IMAGE PROCESSING (ATSIP) *
时燕; 苏涛; 陆莉蓉: "长荡湖围网养殖区长时序时空演变遥感监测", 环境监控与预警 *
朱黎明: "面向灌溉事件的灌溉水量计算方法研究及其在区域上的应用", 中国博士学位论文全文数据库 (农业科技辑) *

Also Published As

Publication number Publication date
CN113793376B (en) 2023-12-05

Similar Documents

Publication Publication Date Title
US7058197B1 (en) Multi-variable model for identifying crop response zones in a field
CN109919083B (en) Early automatic winter wheat mapping method based on Sentinel-2 image data
GB2618896A (en) System and method for crop monitoring
CN112164062A (en) Wasteland information extraction method and device based on remote sensing time sequence analysis
CN104915674A (en) Landsat8 and MODIS fusion-construction high space-time resolution data identification autumn grain crop method
CN110796001A (en) Satellite image film-covering farmland identification and extraction method and system
CN111462223B (en) Sentinel-2 image-based planting area identification method for soybeans and corns in Jianghuai region
CN111345214A (en) Xinjiang cotton region identification method and system based on satellite image data
CN110334880B (en) Effective irrigation area monitoring method
Xu et al. Classification method of cultivated land based on UAV visible light remote sensing
CN114359745A (en) Information processing method, storage medium, and electronic device
She et al. Identification and mapping of soybean and maize crops based on Sentinel-2 data
CN113609899A (en) Backing-out land information positioning display system based on remote sensing time sequence analysis
CN117197668A (en) Crop lodging level prediction method and system based on deep learning
Xu et al. Research and application on corn crop identification and positioning method based on Machine vision
Lee et al. Single-plant broccoli growth monitoring using deep learning with UAV imagery
CN113793376B (en) Irrigation water body extraction method
CN114782835B (en) Crop lodging area proportion detection method and device
CN115880577A (en) Remote sensing technology-based dynamic monitoring method for autumn soil preparation progress of cultivated land
CN116124774A (en) Method for predicting nitrogen content of canopy based on unmanned aerial vehicle spectrum multi-source data
Zhou et al. Monitoring intra-annual and interannual variability in spatial distribution of plastic-mulched citrus in cloudy and rainy areas using multisource remote sensing data
Mashonganyika et al. Mapping of Winter Wheat Using Sentinel-2 NDVI Data
Shao et al. Large-scale mapping of new mixed rice cropping patterns in southern China with phenology-based algorithm and MODIS dataset
Zhang et al. Estimation of peanut seedling emergence rate of based on UAV visible light image
CN114119536B (en) Crop nitrogen content rapid monitoring and diagnosing method based on cloud platform technology

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