CC BY
37
УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА. СЕРИЯ ЕСТЕСТВЕННЫЕ НАУКИ
2020, Т. 162, кн. 2 С.302-313
ISSN 2542-064X (Print) ISSN 2500-218X (Online)
UDC 528.8+528.94
doi: 10.26907/2542-064X.2020.2.302-313
MAPPING CROPLANDS WITH A LONG HISTORY OF CROP CULTIVATION USING TIME SERIES OF MODIS VEGETATION INDICES
M.A. Ivanov, S.S. Mukharamova, O.P. Yermolaev, B. Essuman-Quainoo
Kazan Federal University, Kazan, 420008 Russia
In the study, a method of automated detection of croplands and natural grasslands using multi-year time series of the vegetation indices was developed and implemented. The methodology is based on a new recognition feature enabling identification of all lands cultivated over a long period. The new recognition feature was calculated from time series of the MODIS vegetation indices. It takes into account both intraseasonal and interannual characteristics of the vegetation period for arable lands with spring and winter crops, as well as for natural grasslands. The methodology determines: the procedure of obtaining time series of the NDVI and EVI vegetation indices for the period from March to November in the given years; generalization of the time series by land cover categories and agro-climatic zones; dates for calculating the recognition feature in each zone; calculation of the recognition feature; adaptation of the decision rules to the regional environmental conditions (threshold values for the agro-climatic zones); procedure of classification into croplands and grasslands. Croplands within the territory of the European part of Russia were mapped for the period from 2000 to 2016. The assessment of recognition accuracy and comparison of the results with the existing products confirm that the developed recognition feature can be successfully used for cropland detection and can improve the accuracy of data interpretation when applied along with other features.
Keywords: MODIS, NDVI, EVI, cropland, European part of Russia
During the recent years, monitoring of agricultural land using remote sensing data has become an increasingly urgent issue, from both practical and research perspective. Of major interest are automated mapping of cultivated and abandoned croplands [1, 2], crop recognition [3, 4], and cropland-use intensity assessment [5]. On a regional scale, this requires data with large spatial coverage, high temporal resolution, sufficient spatial resolution, and minimal cost. Data collected by the MODIS instrument installed either on the TERRA and AQUA satellites appear to meet all these criteria. Most standard techniques for cropland mapping are based on the time-series analysis of vegetation indices using parametric (maximum likelihood classifier method [6-8]) and non-parametric (decision tree classifier [8-10], random forest [11-13], support vector machine [1, 14], artificial neural network [5, 8, 15], etc.) classifiers. The most complete list of methods for cropland mapping has been compiled by Chen et al. [16].
This study proposes a new recognition feature for cropland mapping based on the analysis of long-term time series of the MODIS vegetation indices. A method for
Abstract
Introduction
automated detection of crop- and grasslands using the proposed feature was developed and realized for the European part of Russia.
1. Materials and Methods
1.1. Source data. The analysis was carried out on 16-day composites of the NDVI and EVI vegetation indices based on the MODIS data with a spatial resolution of 250 m (MOD13Q1 and MYD13Q1 products), which were obtained for the observation period of 2000-2016 over the entire territory of European Russia covered by seven MODIS granules (h19v02, h19v03, h20v02, h20v03, h20v04, h21v03, and h21v04). We also used spectral data (MOD/MYD09Q1 products) and the MODIS land cover model (MCD12Q1), as well as suitable (cloudless) data from the Landsat satellites for the observation period of 2000-2015 (path 172-174, row 15-29) and the Landsat TM image mosaics 5 by 6 degrees size (Tri-Decadal TM Mosaics) from the US Geological Survey (USGS) [17].
All calculations were carried out with software written in R-language using the following packages: MBA [18], mgcv [19], rgdal [20], and raster [21].
1.2. Preparation of training sites. At the initial stage of the study, we prepared training sites for the land cover of three types: fields planted with spring and winter crops, as well as areas with natural grass vegetation. The Landsat data were used as the reference geodata to differentiate between two land cover types: cropland and land with meadow grass vegetation. Then, we analyzed the multi-season images (the MOD/MYD09Q1 data synthesized in the band1-band2-band1 combination), which made it possible to assign the fields to either spring or winter crops of certain years. In the spring images, spring crops were pink or purple and most often had a uniform texture. In June-July, depending on the region and the farming system, a peak in the vegetation of spring crops was observed; the fields turned green (of various shades, depending on the crop type). After the harvesting period (August, September), the fields with open soil began to take on shades from pink to brown again. Winter crops were clearly visible in the spring images, as well as in the images taken in late autumn; their color was bright green and they had a homogeneous structure. Due to the spatial resolution (250 m) of the MODIS data, the size of the training sites was at least 1 per 1 km, i.e., 4 by 4 pixels. When choosing trainings, the maximum possible regularity of their spatial distribution was ensured. As a result, 4020 training sites were selected over the territory of European Russia for different years: 618 for 2000, 988 for 2008, 782 for 2011, 734 for 2012, and 898 for 2015.
Since the time of meadow grass vegetation and the schedule of crop sowing and harvesting differ between the regions, time series of the vegetation indices must be analyzed separately for different agro-climatic zones. For this reason, each agro-climatic zone was assigned an individual code based on the corresponding map [22]. Moreover, each MODIS granule was processed separately, because their combined mosaic is too large.
As a result, an attributive database was created for the selected trainings. The database includes the following information: Year - the year of the image, in which the training was allocated; Type - the land cover type, to which the training site belongs (1 - spring crops, 2 - winter crops, 3 - natural grass vegetation); Zone - belonging to
particular agro-climatic zones digitized from the map of agro-climatic resources of the USSR [22]; Scene - the nomenclature of the granule.
1.3. Analysis of time series of the vegetation indices. For each training, using 16-day composites of the MODIS vegetation indices, the NDVI and EVI values were assigned to all selected dates of the corresponding year. The index value for a certain date was considered acceptable if the pixel value on reliability raster equaled 0 (good data) or 1 (marginal data). Otherwise, the value of the vegetation index attribute field was considered 0 and this record was excluded from further analysis.
With the help of the obtained database, the NDVI and EVI time series were analyzed for each land cover type in all agro-climatic zones and during all years under consideration. The agro-climatic zones with fewer than 20 trainings of each type were excluded from the analysis. These are the focal agriculture zone and the northward zones.
The next step was to analyze of the seasonal dynamics of vegetation indices in training pixels. For each training, according to MOD/MYD13Q1, the NDVI and EVI time series were plotted for the period from March to November of the corresponding year (taking into account the quality band values). A total 8040 time series were plotted (2 indices for 4020 trainings). The time series were generalized for each land cover type and agro-climatic zone using the average, median, and smoothing splines included in the method of generalized additive models [19]. The seasonal dynamics of the vegetation indices in 2008 for the training sites of all three land cover types in one of the studied zones is shown in Fig. 1. Here, the dynamics of the average, median, and model approximation of each index can also be seen.
Fig. 1. Seasonal dynamics of the NDVI and EVI indices obtained using the MODIS data in 2008 for the fields with spring and winter crops, as well as for the natural grasslands (the zone with the accumulated temperatures above 10 °C from 2800 °C to 3400 °C)
100 150 200 250 300
day of the year
150 200 250 300
day of the year
- Spring crops - Winter crops-Natural grasslands
100 150 200 250 300 100 150 200 250 300
day of the year day of the year
- Spring crops-Winter crops
Fig. 2. Seasonal dynamics of the NDVI and EVI vegetation indices generalized by three types of land cover (a) and difference between the indices of spring and winter crops (b) in 2008 for the zone with the accumulated temperatures above 10 °C from 2800 °C to 3400 °C
The behavior of the generalized time series for all three land cover types in the same axis and the dynamics of difference between the indices values of spring and winter crops are presented in Fig. 2. The plots clearly demonstrate the seasonal variations in the vegetation indices of spring and winter crops, which are also typical for other years. The following recognition feature, which can separate croplands from grasslands was proposed based on the analysis of the dynamics of the indices. On the plots (Fig. 2, b), we can see that the absolute value of difference between the indices for spring and winter crops has two maximums: the first one occurs in May, while the other takes place in July. Due to the rotation of spring and winter crops in the arable area, the difference between the vegetation indices calculated on these dates will have greater variability from year to year than in the area with natural grass vegetation.
As a measure of variability, a variance or standard deviation of the difference between the values of the vegetation index on these dates as estimated over 17 years (from 2000 to 2016) can be assessed. The dispersion for croplands should be greater than for grasslands.
1.4. Calculation of the recognition feature and classification. To calculate the recognition feature, we determined two dates for each agro-climatic zone of European Russia that are the minimum and maximum values of the difference between the indices of spring and winter crops (based on the analysis of the NDVI and EVI generalized time series in training pixels), respectively. Dates are based on the MOD/MYD13Q1 time resolution (Table 1).
The calculation of the recognition feature - the standard deviation (SD) of the difference between the values of the vegetation index on these dates was carried out both for NDVI and EVI. At the same time, the study area was initially masked to exclude lands
Table 1
Dates corresponding to the minimum and maximum values of the difference between the indices of spring and winter crops
Agro-climatic zone code NDVI EVI
Date (day of the year) of minimum difference Date (day of the year) of maximum difference Date (day of the year) of minimum difference Date (day of the year) of maximum difference
3 May 17 (137) July 28 (209) May 17 (137) July 28 (209)
4 May 9 (129) Aug 5 (217) May 17 (137) July 28 (209)
5 May 9 (129) July 28 (209) May 17 (137) July 20 (201)
6 May 9 (129) July 20 (201) May 9 (129) July 20 (201)
7 May 1 (121) July 20 (201) May 1 (121) July 20 (201)
8 May 1 (121) Aug 13 (225) May 1 (121) Aug 5 (217)
0 0.15 0.2 0.4 0.6
Fig. 3. Standard deviation of the difference of EVI values: territory of the Melekessky district of the Ulyanovsk region (on the left) and the Alekseevsky district of the Volgograd region (on the right)
that obviously did not belong to the category of croplands or grasslands using the MODIS land cover model (MCD12Q1). For each pixel of the unmasked territory, the standard deviations sdNDVI and sdEVI were calculated, for which: 1) a particular agro-climatic zone and two corresponding dates were defined for each pixel; 2) in each year from 2000 to 2016, differences between the index values on these two dates were calculated; 3) the standard deviation of the difference was calculated.
Some fragments of the result for sdEVI are shown in Fig. 3. In the image, typical regular spatial structures with high sdEVI values (brown) corresponding to cultivated croplands are clearly visible. Low values of sdEVI, which are displayed in green, correspond to grasslands and are mostly allocated along the river valleys.
The assignment of pixels to croplands and grasslands was carried out using the threshold classification. For a more accurate adaptation of the decision rules to the regional environments of European Russia, a detailed zoning of the territory was carried out. The accumulated temperatures above 10 °C were used for zoning. To calculate this indicator, open data from the RosHydromet weather stations (long-term daily observations of the air temperature and precipitation at weather stations in Russia and the former USSR) were processed [23]. Spatial interpolation to the study area was performed using the multilevel B-spline approximation method [18, 24]. Following that, zones with a step of 250 °C were allocated.
Natural grasslands ■ Croplands
Fig. 4. Estimates of the density functions of conditional probability: left column -P(sdNDVj | cropland) and P(sdNDVl | grassland); right column - P(sdEVl | cropland) and P(sdEV1 | grassland); top row - estimate of the entire sample of training pixels; middle row -estimate by subsample within the zone with the accumulated temperatures above 0 °C from 2500 °C to 2750 °C; bottom row - estimate by subsample within the zone with the accumulated temperatures above 0 °C from 3000 °C to 3250 °C
The sample of pixels belonging to the training sites allowed us to explore the distribution of values of the sdNDVI and sdEVI recognition features for the classes of croplands and grasslands, to assess the degree of their separability, and to determine threshold values.
In Fig. 4, the estimates are given for the conditional probability density functions (| cropland), | grassland), P(sdEm | cropland), P(sdEm | grassland))
which were obtained both for the entire sample and for subsamples within different zones. In the parts of this figure, the red histograms correspond to the sample of training pixels of croplands, while the green histograms indicate the sample of training pixels of grasslands.
The Hellinger distance was calculated as a separation measure of distributions [25]:
H = ( p, P2) =
(
1 - J yjp (u), p (u) du
Y
where p (u), P2 (u) are conditional probabilities of the classification feature.
The H value can vary from 0 (complete coincidence of histograms) to V2 (complete separability of histograms). In our case, the Hellinger distance for both features (sdNDVI and sdEVI) varies from 0.7 to 1.3 depending on the zone. The separability of classes is satisfying (the lower threshold of the Hellinger distance is normally 0.4). Then, an independent classification was carried out for each of the classification features (sdNDVI and sdEVI). The results of the independent classifications were integrated: only the pixels assigned to the corresponding class by both recognition features were classified as croplands and grasslands.
2. Results and Accuracy Assessment
As a result of the classification, a raster layer of croplands and grasslands for the European territory of Russia was obtained (Fig. 5).
20°E 30°E 40°E 50°E 60°E 70°E
Fig. 5. Cropland recognition results
To assess the accuracy of mapping, the results of manual digitizing of croplands from the multi-seasonal Landsat 8 images for 2014-2016 at eight control sites were used: 1 - a part of the Izh River basin (Udmurt Republic); 2 - the Mesha River basin (Republic of Tatarstan); 3 - a part of the Sviyaga River basin (Ulyanovsk region); 4 - the Ulema River basin (Republic of Tatarstan); 5 - the Veduga River basin (Voronezh region); 6 - a part of the Medveditsa River basin (Saratov region); 7 - a part of the Samara River basin (Orenburg region); 8 - a part of the Kuma River basin (Stavropol territory) [26].
To assess the accuracy of cropland mapping, the boundaries obtained by automated classification were compared with the boundaries of croplands obtained by the manual vectorization of Landsat 8. The sites of agreement and disagreement were identified. Their areas were calculated. An error matrix was created. The accuracy and omission and commission errors were calculated (Table 2).
Table 2
Accuracy and errors
River basin Izh Mesha Sviyaga Ulema Veduga Medveditsa Samara Kuma
User's accuracy, % 57.76 74.44 75.59 74.77 74.31 80.44 79.94 77.04
Producer's accuracy,% 59.79 83.69 78.84 90.42 85.67 89.09 78.17 84.62
Commission error, % 42.24 25.56 24.41 25.23 25.69 19.56 20.06 22.96
Omission error, % 40.21 16.31 21.16 9.58 14.33 10.91 21.83 15.38
Table 3
Accuracy of cropland recognition on the TerraNorte RLC model and errors
River basin Izh Mesha Sviyaga Ulema Veduga Medveditsa Samara Kuma
User's accuracy, % 77.48 85.88 81.93 84.63 83.67 88.74 90.00 83.68
Producer's accuracy,% 43.79 78.75 84.26 79.48 79.38 82.22 70.53 77.22
Commission error, % 22.52 14.12 18.07 15.37 16.33 11.26 10.00 16.32
Omission error, % 56.21 21.25 15.74 20.52 20.62 17.78 29.47 22.78
The mapping accuracy in almost all cases is higher than 75%. Thus, the proposed recognition feature that separates croplands from grasslands can be used among a number of other features.
In addition, the mapping results were compared with the boundaries of croplands obtained from the TerraNorte RLC model (the land cover map of Russia) [27]. These boundaries were also compared with the data obtained by the Landsat 8 vectorization. An error matrix was created, and the accuracy and errors were estimated (Table 3).
As evident from Tables 2 and 3, for all control sites, the user's accuracy for the obtained classification results is lower than the same accuracy for the data from the TerraNorte RLC map, while the producer's accuracy is higher (except for the Svi-yaga River basin). This means that the results of the classification carried out according to the proposed methodology have fewer omission errors, i.e., fewer pixels belonging to croplands were skipped. Moreover, a greater number of pixels, which do not belong to croplands according to the control data, were classified as croplands.
The differences in accuracy can be explained by the different duration of time series of the used input and verification data. Since the MODIS data obtained in 2000-2016 were used in this study for classification, the result is the boundaries of all croplands cultivated during this period. The borders of the cultivated croplands on the TerraNorte RLC model are relevant for 2014, so its area here is obviously smaller. Furthermore, verification data were obtained from the 2014-2016 images, i.e., they almost coincide in date with the TerraNorte model. Taken together, these two points led to higher values of commission errors in the results obtained in this work.
Acknowledgments. The work was supported by the Russian Science Foundation (project no. 19-17-00064).
References
1. Alcantara C., Kuemmerle T., Prishchepov A.V., Radeloff V.C. Mapping abandoned agriculture with multi-temporal MODIS satellite data. Remote Sens. Environ., 2012, vol. 124, pp. 334-347. doi: 10.1016/j.rse.2012.05.019.
2. Bartalev S.A., Plotnikov D.E., Loupian E.A. Mapping of arable land in Russia using multi-year time series of MODIS data and the LAGMA classification technique. Remote Sens. Lett., 2016, vol. 7, no. 3, pp. 269-278. doi: 10.1080/2150704X.2015.1130874.
3. Wardlow B.D., Egbert S.L., Kastens J.H. Analysis of time-series MODIS 250 m vegetation index data for crop classification in the U.S. Central Great Plains. Remote Sens. Environ.,
2007, vol. 108, no. 3, pp. 290-310. doi: 10.1016/j.rse.2006.11.021.
4. Wardlow B.D., Egbert S.L. Large-area crop mapping using time-series MODIS 250 m NDVI data: An assessment for the U.S. Central Great Plains. Remote Sens. Environ.,
2008, vol. 112, no. 3, pp. 1096-1116. doi: 10.1016/j.rse.2007.07.019.
5. Estel S., Kuemmerle T., Levers Chr., Baumann M., Hostert P. Mapping cropland-use intensity across Europe using MODIS NDVI time series. Environ. Res. Lett., 2016, vol. 11, no. 2, art. 024015, pp. 1-10. doi: 10.1088/1748-9326/11/2/024015.
6. Arvor D., Jonathan M., Meirelles M.S.P., Dubreuil V., Durieux L. Classification of MODIS EVI time series for crop mapping in the state of Mato Grosso, Brazil. Int. J. Remote Sens., 2011, vol. 32, no. 22, pp. 7847-7871. doi: 10.1080/01431161.2010.531783.
7. Bartalev S.A., Egorov V.A., Lupyan E.A., Plotnikov D.E., Uvarov I.A. Arable land mapping based on long-term satellite data from a MODIS spectroradiometer and locally adaptive classification. Komp'yut. Opt., 2011, vol. 35, no. 1, pp. 103-116. (In Russian)
8. Shao Y., Lunetta R.S. The use of MODIS-NDVI data for mapping cropland across the Great Lakes Basin, USA. Proc. IGARSS 2008 - 2008 IEEE Int. Geosci. Remote Sens. Symp. Boston, Mass., 2008, pp. 196-199. doi: 10.1109/IGARSS.2008.4780061.
9. Wardlow B.D., Egbert S.L. A comparison of MODIS 250-m EVI and NDVI data for crop mapping: A case study for southwest Kansas. Int. J. Remote Sens., 2010, vol. 31, no. 3, pp. 805-830. doi: 10.1080/01431160902897858.
10. Brown J.C., Kastens J.H., Coutinho A.C., Victoria D. de C., Bishop C.R. Classifying mul-tiyear agricultural land use data from Mato Grosso using time-series MODIS vegetation index data. Remote Sens. Environ., 2013, vol. 130, pp. 39-50. doi: 10.1016/j.rse.2012.11.009.
11. Estel S., Kuemmerle T., Alcántara C., Levers Chr., Prishchepov A., Hostert P. Mapping farmland abandonment and recultivation across Europe using MODIS NDVI time series. Remote Sens. Environ., 2015, vol. 163, pp. 312-325. doi: 10.1016/j.rse.2015.03.028.
12. Hao P., Zhan Y., Wang L., Niu Z., Shakir M. Feature selection of time series MODIS data for early crop classification using random forest: A case study in Kansas, USA. Remote Sens., 2015, vol. 7, no. 5, pp. 5347-5369. doi: 10.3390/rs70505347.
13. Zhong L., Yu L., Li X., Hu L., Gong P. Rapid corn and soybean mapping in US Corn Belt and neighboring areas. Sci. Rep., 2016, vol. 6, art. 36240, pp. 1-14. doi: 10.1038/srep36240.
14. Alcantara C., Kuemmerle T., Baumann M., Bragina E.V., Griffiths P., Hostert P., Knorn J., Müller D., Prishchepov A.V., Schierhorn F., Sieber A., Radeloff V.C. Mapping the extent of abandoned farmland in Central and Eastern Europe using MODIS time series satellite data. Environ. Res. Lett., 2013, vol. 8, no. 3, art. 035035, pp. 1-9. doi: 10.1088/17489326/8/3/035035.
15. Lunetta R.S., Shao Y., Ediriwickrema J., Lyon J.G. Monitoring agricultural cropping patterns across the Laurentian Great Lakes Basin using MODIS-NDVI data. Int. J. Appl. Earth Obs. Geoinf., 2010, vol. 12, no. 2, pp. 81-88. doi: 10.1016/j.jag.2009.11.005.
16. Chen Y., Lu D., Moran E., Batistella M., Dutra L.V., Sanches I.D., da Silva R.F.B., Huang J., Luiz A.J.B., de Oliveira M.A.F. Mapping croplands, cropping patterns, and
crop types using MODIS time-series data. Int. J. Appl. Earth Obs. Geoinf., 2018, vol. 69, pp. 133-147. doi: 10.1016/j.jag.2018.03.005.
17. U.S. Geological Survey. Available at: https://www.usgs.gov/core-science-systems/nli/ landsat/data-tools.
18. Finley A., Baneijee S., Hjelle O., Bivand R. MBA: Multilevel B-spline approximation. R package version 0.0-9, 2017. Available at: https://cran.r-project.org/web/packages/ MBA/index.html.
19. Wood S.N. Generalized Additive Models: An Introduction with R. CRC Press, 2006. 410 p.
20. Bivand R., Keitt T., Rowlingson B., Pebesma E., Sumner M., Hijmans R., Rouault E., War-merdam F., Ooms J., Rundel C. rgdal: Bindings for the geospatial data abstraction library. R package version 1.2-18, 2018. Available at: https://cran.r-project.org/web/ packages/rgdal/index.html.
21. Hijmans R.J., van Etten J., Cheng J., Mattiuzzi M., Sumner M., Greenberg J.A., Lami-gueiro O.P., Bevan A., Racine E.B., Shortridge A., Ghosh A. Raster: Geographic Data Analysis and Modeling. R package version 2.6-7, 2017. Available at: https://cran.r-project.org/web/packages/raster/index.html.
22. Geograficheskii atlas dlya uchitelei srednei shkoly [Geographic Atlas for Secondary School Teachers]. Moscow, Gl. Upr. Geod. Kartogr. Sov. Minist. SSSR, 1982. 238 p. (In Russian)
23. Bulygina O.N., Razuvaev V.N., Aleksandrova T.M. Description of the dataset of daily air temperature and precipitation for weather stations of Russia and the former USSR (TTTR). Database State Registration Certificate no. 2014620942. 2014. (In Russian)
24. Lee S., Wolberg G., Shin S.Y. Scattered data interpolation with multilevel B-splines.
IEEE Trans. Visualization Comput. Graphics, 1997, vol. 3, no. 3, pp. 228-244. doi: 10.1109/2945.620490.
25. Upton G., Cook I. Oxford Dictionary of Statistics. Oxford, Oxford Univ. Press, 2002. 420 p.
26. Ivanov M.A., Prishchepov A.V., Golosov V.N., Zalyaliev R.R., Efimov K.V., Kond-rat'eva A.A., Kinyashova A.D., Ionova Yu.K. Methods for mapping croplands dynamics in river basins of the European part of Russia for 1985-2015. Sovrem. Probl. Distantsionnogo Zondirovaniya Zemli Kosmosa, 2017, vol. 14, no. 5, pp. 161-171. doi: 10.21046/20707401-2017-14-5-161-171. (In Russian)
27. Bartalev S.A., Egorov V.A., Loupian E.A., Khvostikov S.A. A new locally-adaptive classification method LAGMA for large-scale land cover mapping using remote-sensing data. Remote Sens. Lett., 2014, vol. 5, no. 1, pp. 55-64. doi: 10.1080/2150704X.2013.870675.
Recieved April 7, 2020
Ivanov Maxim Andreevich, PhD in Geography, Associate Professor of Landscape Ecology Department
Kazan Federal University
ul. Kremlyovskaya 18, Kazan, 420008 Russia E-mail: maximko-87@mail.ru
Mukharamova Svetlana Sayasovna, PhD in Biology, Associate Professor of Department of Ecological
Systems Modeling
Kazan Federal University
ul. Kremlyovskaya 18, Kazan, 420008 Russia E-mail: smukhara@gmail.com
Yermolaev Oleg Petrovich, Doctor of Geography, Professor of Landscape Ecology Department
Kazan Federal University
ul. Kremlyovskaya 18, Kazan, 420008 Russia E-mail: oyermol@gmail.com
Essuman-Quainoo Benedict, Graduate Student of Landscape Ecology Department
Kazan Federal University
ul. Kremlyovskaya 18, Kazan, 420008 Russia E-mail: benessquainoo@gmail.com
УДК 528.8+528.94 doi: 10.26907/2542-064X.2020.2.302-313
Распознавание пахотных угодий длительного пользования по данным временных рядов вегетационных индексов MODIS
М.А. Иванов, С.С. Мухарамова, О.П. Ермолаев, Б. Эссуман-Куэйну Казанский (Приволжский) федеральный университет, г. Казань, 420008, Россия
Аннотация
В статье предложена методика автоматизированного распознавания пахотных угодий, используемых в многолетний период, по временным рядам вегетационных индексов. Разработан новый дешифровочный признак, который одновременно учитывает внутрисезонные и межгодовые особенности вегетации, с целью отделения пахотных угодий от естественного травостоя. Для этого с помощью методов геостатистики и геоинформационного пространственного анализа были решены следующие задачи: формирование архива данных ДЗЗ, включая вегетационные индексы NDVI и EVI MODIS за период 2000-2016 гг.; создание обучающих выборок для яровых и озимых культур и естественной луговой растительности; построение и анализ временных рядов вегетационных индексов с учетом агроклиматического зонирования; расчет дешифровочного признака; распознавание пахотных земель и оценка точности. На основе разработанного признака проведено распознавание пахотных угодий на территории Европейской части России. Точность результатов дешифрирования оценивалась путем сопоставления с результатами ручной векторизации пахотных угодий по данным Landsat на восьми контрольных участках, расположенных в различных ландшафтных зонах. Оценка проводилась с помощью матрицы ошибок. Кроме того, результаты сравнивались с моделью земного покрова TerraNorte RLC, полученной другим методом. Рассчитанные значения показали высокую точность и сопоставимость данных. Таким образом, разработанный дешифровочный признак может использоваться как самостоятельно, так и совместно с другими признаками для повышения точности дешифрирования. Научно-практической ценностью обладают и полученные пространственные данные о пахотных угодьях, используемых на протяжении большого периода времени.
Ключевые слова: MODIS, NDVI, EVI, пахотные угодья, Европейская часть России
Благодарность. Исследование выполнено в рамках научного проекта Российского научного фонда № 19-17-00064.
Поступила в редакцию 07.04.2020
Иванов Максим Андреевич, кандидат географических наук, доцент кафедры ландшафтной экологии
Казанский (Приволжский) федеральный университет ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: maximko-87@mail.ru
Мухарамова Светлана Саясовна, кандидат биологических наук, доцент кафедры моделирования экологических систем
Казанский (Приволжский) федеральный университет ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: smukhara@gmail.com
Ермолаев Олег Петрович, доктор географических наук, профессор кафедры ландшафтной экологии
Казанский (Приволжский) федеральный университет ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: oyermol@gmail.com
Эссуман-Куэйну Бенедикт, аспирант кафедры ландшафтной экологии
Казанский (Приволжский) федеральный университет ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: benessquainoo@gmail.com
<For citation: Ivanov M.A., Mukharamova S.S., Yermolaev O.P., Essuman-Quainoo B. Mapping croplands with a long history of crop cultivation using time series of MODIS vegetation indices. Uchenye Zapiski Kazanskogo Universiteta. Seriya Estestvennye Nauki, 2020, vol. 162, no. 2, pp. 302-313. doi: 10.26907/2542-064X.2020.2.302-313.
<Для цитирования: Ivanov M.A., Mukharamova S.S., Yermolaev O.P., Essuman-Quainoo B Mapping croplands with a long history of crop cultivation using time series of MODIS vegetation indices // Учен. зап. Казан. ун-та. Сер. Естеств. науки. - 2020. - Т. 162, кн. 2. -С. 302-313. - doi: 10.26907/2542-064X.2020.2.302-313.
