INTRODUCTION
Monitoring land cover changes and particularly water bodies represents an important contemporary issue of our society. It means that society is interested in new findings, which can help to understand how the Earth and its components (biophysical spheres, e.g. hydrosphere and biosphere) behave. Understanding functions, relationships, and rules can, beside others, support landscape management and further sustainable development. Studies focusing on land cover/land use change are often realized at national or international level. The local level of the monitoring of land cover/land use change and water bodies is important too but it is not so frequently covered by the studies because it may require a higher spatial resolution that satellite imagery provides.
Remotely sensed images are very often used for observation and monitoring changes of landscape and terrain. There are several reasons why remote sensing is more reasonable than in situ observations, sampling and measurements, and land surveying. The costs are usually lower (namely in the case when reasonably large areas are monitored or the area is monitored periodically), measurement does not require a direct access to the monitored area (e.g. useful in the case of swamplands), accuracy and spatial resolution is adequate for monitoring, and data measured in various parts of the electromagnetic spectrum are available (in the case of proper sensors). Multispectral and hyperspectral imagery is important for land cover monitoring because of possibility of calculation of various spectral indices, which are suitable for land cover identification. Satellite or aerial imagery is suitable for monitoring larger areas. UAVs are increasingly used to monitor smaller areas, e.g. small water bodies (ponds), mainly because of faster data collection, data collection on demand, higher spatial resolution, and lower costs in comparison to aerial imagery collection (Iizuka et al., 2018; Pásler, Komárková and Sedlák, 2015; Salamí, Barrado and Pastor, 2014).
The aim of the paper is to propose a suitable way of monitoring small water bodies and land cover changes in the surroundings, which is based on data collected by a middle-class UAV and visible spectral indices (based on red, green, and blue bands). It extends the study of Cermakova, Komarkova and Sedlak (2019). At first, a suitable way of utilization of a UAV for data collection is described. Next, several visible sspectral indices (CIVE, ExG, ExR, ExG – ExR, GLI, and VARI) are calculated from the collected data. Evaluation of indices is the following part of the paper. Only a part of the selected pond and its surroundings are used to demonstrate the proposed procedure.
STATE OF THE ART
Water Bodies and Land Cover Changes
Identification of land use/land cover changes and their quality monitoring have been based on remotely sensed data since the 70s along with Landsat existence (Bukata, Harris and Bruton, 1974; Lin et al., 2014; Maleki et al., 2018; Pechanec et al., 2014; Saharan et al., 2018). Landsat data was successfully used for water bodies observation and identification of their shorelines many times, e.g. (Bhagat and Sonawane, 2011; Jiang et al., 2014; Jones et al., 2017; Nandi et al., 2016).
Many articles were particularly focused on various water bodies and their shorelines. In some cases, the observation takes many years, e.g. (Gallop et al., 2015). This study observed changes of shoreline in the period of 34 years in southwestern Australia. Aerial photographs were used as a data source. It was realized that over 96 % of beaches straightened during the time period. In this case, there were large changes in shoreline position that did not correspond to the indices of annual wave conditions. In future, researches want to use fuzzy-based computer algorithms. (Nandi et al., 2016) used the Digital Shoreline Analysis System (DSAS) for the detection of changes of shoreline in a monitored area. DSAS is a special extension of United States Geological Survey, which can calculate the rate of shoreline change.
Small water inland bodies identification from remotely sensed data was emphasized as a specific issue by several authors, e.g. (Cermakova, Komarkova and Sedlak, 2019; Haas, Bartholomé and Combal, 2009; Jones et al., 2017; Li et al., 2015; Ogilvie et al., 2018; Pásler, Komárková and Sedlák, 2015). In the Czech Republic, e.g. water erosion of the river Trkmanka shoreline (a small stream in South Moravia) was studied in (Pechanec et al., 2015). They focused on erosion susceptibility and shoreline changes because of the destructive influences of water erosion on the soil surface. This method provides the maximum volume of soil erosion caused by water.
Unmanned Aerial Vehicles for Monitoring
UAVs can provide very low elevation and very high-resolution data on demand, of course with respect to limitations given by law regulations and weather conditions (Iizuka et al., 2018; Lejot et al., 2007; Pásler, Komárková and Sedlák, 2015).
UAVs are increasingly used for monitoring land cover/land use changes (Feng, Liu and Gong, 2015; Iizuka et al., 2018; Salamí, Barrado and Pastor, 2014) and natural disasters, e.g. floods or earthquakes (Abdelkader et al., 2013; Casado et al., 2018). It is a reasonably cheap and fast possibility of how to determine the situation in the observed area. It can be used as an early warning before e.g. landslide or pond failure, or to fast monitor impacts of a disaster because the changes of monitored area and situation are observable immediately. So, utilization of UAVs for prevention and monitoring actual situation can save money and peoples’ lives. Using an UAV for floods warning is also possible according to (Srikudkao et al., 2015). The study is focused on creating a basis for floods warning systems and prevention systems. (Abdelkader et al., 2013) presented a novel flood sensing architecture to monitor large scale desert hydrological basins surrounding metropolitan areas, based on unmanned air vehicles. The reason for their proposal was that the satellite-based weather forecast was not enough reliable. On the opposite, a combination of UAV and mobile microsensors using many parameters provided more accurate results.
Spectral visible indices based on red (R), green (G), and blue (B) bands are increasingly used to process data from low-cost and middle-class UAVs. Many indices have been proposed for an agriculture and phenological purposes – to observe vegetation, its growth, and its health status, e.g. CIVE – Colour Index of Vegetation Extraction (Kataoka et al., 2003), ExG – Excess Green (Woebbecke et al., 1995), ExR – Excess Red (Meyer, Hindman and Laksmi, 1999), ExG - ExR difference (Meyer et al., 2004), GLI – Green Leaf Index (Louhaichi, Borman and Johnson, 2001), NExG – Normalized Excess Green (Woebbecke et al., 1995), NGRDI – Normalized Red-Green Difference Index (Tucker, 1979), RGBVI – Red-Green-Blue Vegetation Index (Bendig et al., 2015), VARI – Visible Atmospherically Resistant Index (Gitelson et al., 2002), and many others.
Shorelines can be used to autonomously navigate UAVs as proposed by (Baker and Kamgar-Parsi, 2010). A near-infrared sensor allows the vision system to distinguish water from land. The vision system then issues commands to the autopilot to follow the coastline or the riverbank.
CASE STUDY
Methods and Software
The following visible spectral indices are calculated in the paper: CIVE, ExG, ExR, ExG – ExR, GLI, and VARI. All used indices, their formulas, and sources describes Table 1. GLI is focused on using reflectance of green vegetation cover to highlight it in the images (Louhaichi, Borman and Johnson, 2001). VARI is used to highlight vegetation in images as well (Gitelson et al., 2002). ExG uses the green part of the electromagnetic spectrum to highlight the green vegetation cover (Woebbecke et al., 1995).
DJI Ground Station on iOS was used for flights planning and DJI GO 4 for flights controlling. ENVI OneButton 5.0 was used for creating mosaics from particular images. QGIS Desktop (version 3.8.1) was used for calculation of vegetation indices. ArcGIS for Desktop (version 10.5.1) was used for next processing and cartographic visualisations.
Table 1. Used spectral indices and their formulas
Index
|
Formula
|
|
CIVE
|
0.441*R – 0.81*G + 0.385*B + 18.78745
|
(Kataoka et al., 2003)
|
ExG
|
2*g – r – b
components are normalized:
g = G/(R + G + B)
r = R/(R + G + B)
b = B/(R + G + B)
|
(Woebbecke et al., 1995)
|
ExR
|
1.3*R – G
|
(Meyer, Hindman and Laksmi, 1999)
|
ExG – ExR (ExGR)
|
ExG – ExR
|
(Meyer et al., 2004)
|
GLI
|
(2*G – R – B)/(2*G + R + B)
|
(Louhaichi, Borman and Johnson, 2001)
|
VARI
|
(G – R)/(G + R – B)
|
(Gitelson et al., 2002)
|
|
Area of Interest
A western part of the pond Skrin, which is located close to city Lázně Bohdaneč city (the Pardubicky region, the Czech Republic), is used as an area of interest (see Figure 1). The total area of the pond is 0.27 km2 and the pond is used for fish farming (Czech Fishing Union, 2018). The area of interest is a flat area, which lies approx. 220 m a.s.l. It covers clear water surface, seasonally flooded greenery, various vegetation including tree crowns, dry reeds (including dry grass), and dry trees.
The appearance of the area during the monitored period was changing (see Figure 2). During the first flight (February 9, 2018), a part of the pond Skříň was frozen. Simultaneously, some parts were covered by snow and there was no green vegetation. The whole area of the pond was frozen on March 8, 2018. Start of vegetation growth is visible in the last image – acquired on April 20, 2018.
Data Collection and Pre-Processing
Phantom 3 Professional with built-in 2,7K Ultra HD DJI camera was used for data collection. Camera has 1/2.3” CMOS sensor with 12.4 M of effective pixels. Images from the DJI camera fixed on the Phantom 3 contain the visible bands (red, green, and blue) only. Phantom 3 has the following characteristics: weight 1 216 g, 4 motors, max. speed 16 m.s-1, max. range 1 000 m, max. flight time 25 minutes. (DJI, 2018). The drone itself (including the DJI camera) costs approx. 700 – 800 Eur (DJI, 2018).
Images were acquired once per a month. Particular dates were influenced by weather conditions and vegetation changes, which are fast in April. The Pond Skříň was observed in the following dates: February 9, 2018; March 8, 2018; and April 20, 2018. In all cases, images were acquired between 9 – 10 am. The route of the recording was planned in advance in DJI Ground Station and it was the same in all cases. All images were acquired during one flight lasting approx. 15 min. The total length of the flight was 1 545 m, size of the covered area was 0.0285 m2. The flight was planned in 7 lines with both front and side overlap 60 % and at 39.6 m of altitude. In total, 64 images were acquired during each flight. Spatial resolution is 1.7 cm per pixel.
The atmospheric correction was not calculated because of the low altitude, a short time span, and good weather conditions. Information about coordinates is recorded during the flight by the software. Verification of the location was conducted by a comparison with aerial imagery provided by the State Administration of Land Surveying and Cadastre in a form of WMS at the first stage. ENVI OneButton is used as a tool to correct particular images and mosaic them to create an output orthorectified mosaic with WGS 84 – UTM zone 33 as the coordinate system.
Data Analysis – Visible Spectral Indices
The following indices were calculated: CIVE, ExG, ExR, ExG - ExR, GLI, and VARI for all three mosaics. Next, low-pass filtering on a 3×3 window (ArcGIS tool Filter) was used as a post-processing tool. Visualization of the results was the last step. Visualization in a form of classifies raster was chosen and ‘Natural breaks – Jenks’ was used as a classification method in ArcGIS. ‘Pink to YellowGreen Diverging, Dark’ scheme was used for visualization purposes. To simplify interpretation, we preferred to use pink/violet colour for all water forms and green for all vegetation forms.
Because of the different appearance of the landscape in the three involved months and limited possibilities of visible spectral indices, each index provides the most information value for another time period and land cover. Figures 3 – 5 show results.
DISCUSSION
Cermakova, Komarkova and Sedlak (2019) shows calculations of the land cover changes in the area of interest in the observed period. It is clear that the changes are significant. In this article we extend utilization of visible spectral indices by involving other indices, by low-pass filtering and by a fast but yet useful visualization.
Difference ExG – ExR together with VARI allows a highlighting of frozen water surface. Almost all indices support green vegetation identification as it can be seen in Figure 5. CIVE, ExG, and GLI well highlight a green vegetation. CIVE and VARI are able to highlight a clear water surface. All of the indices allow distinguishing a seasonally flooded vegetation from the surrounding land covers.
The results were compared with visual interpretation and field observation. A small size of the area of interest and a very-high spatial resolution of the images allows utilization of this method. After comparison of both methods, it is clear that land cover changes can be observed in the monitored time period.
The key question is if regional state administration bodies will be able to use UAVs as a tool to quickly and cheaply identify land cover changes, not only of the small water bodies. Or if the state administration bodies will still use traditional data sources, which provide lover spatial resolution, or which may be more expensive and less frequently available.
CONCLUSIONS
Land use and land cover changes monitoring is increasingly important with respect to a sustainable development. Climate changes and natural disasters, e.g. droughts in the last years, can be used as example of monitoring drivers. Water and water sources belong to the key contemporary sources and issues. Researches are still mainly focused on a world or national level of the water bodies monitoring, but the regional level of this issue is important as well, mainly for the citizens of the region.
This is one of the reasons why the article is focused on a small water body – a pond Skříň and its surroundings, which is located close to Lázně Bohdaneč city in the Czech Republic. A middle-class UAV (Phantom 3 Pro) was used for the data acquisition because it represents a quick and cheap way how to collect the spatial data with a very-high spatial resolution. In this case, the images were acquired from the camera sensing in the visible bands only.
The monitoring of the part of the pond and its surrounding was conducted 3 times: February 9, 2018; March 8, 2018; and April 20, 2018. For the study of changes, visible spectral indices were used. The indices are based on different reflectivity of various land covers and mostly highlight the green vegetation. All of them require the availability of red, green, and blue bands only.
In the article, there are calculated various visible spectral indices, which use only visible bands. The calculated indices are the following ones: CIVE, ExG, ExR, ExG - ExR, GLI, and VARI.
Visual comparison was used in this case because of the small size of the area and high resolution of the images. It is clearly visible that the water area changed during the monitored time period.
The article shows that it is possible to identify land cover types and to monitor land cover changes on demand and using a low-cost middle-class UAV. UAV allows data acquisition on demand, of course with a respect to legal and weather limitations. Spatial resolution of acquired data is very high (centimetres level). Middle-class UAVs are usually equipped with a standard visible (RGB) camera. The article demonstrates the ease of utilization of this type of UAV for land cover identification and its changes monitoring.
ACKNOWLEDGEMENT
Authors thank the University of Pardubice, Students grant competition, SGS_2019_17 project for the support.