1. Introduction
Satellites equipped with a synthetic aperture radar (SAR) are being used in various applications utilizing an analysis of signal backscatter intensity and eventually its change in time. An indisputable advantage of this active microwave technology is the ability to work at any time of the day and under any meteorological conditions. Satellite-based flood detection is one of the applications which can enormously benefit from this advantage as optical imagery often suffers from high cloud coverage during flood events.
The flood detection from SAR imagery is based on the fact that the SAR signal backscattered by water bodies has a very low intensity. Therefore, a SAR intensity threshold allowing the pixels representing the water surface to be distinguished from others (mainly corresponding to land) is sought in most of the references (such as [
1,
2,
3]). Hereafter, we call this threshold a water–land threshold. According to the cited references, the threshold, which depends on many factors, including meteorological conditions (wind, temperature), differs for each image, and has to be estimated from the evaluated data. The estimation of this water–land threshold is, therefore, a crucial part of the algorithm.
To estimate this water–land threshold, authors in [
2,
3] firstly used a method based on a bimodal histogram. To handle cases where only a small part of the AOI is flooded, they divided the AOI into tiles and searched for the threshold in each tile individually. We suppose this tiling method is suitable for very high-resolution imagery (such as TerraSAR-X for which it was developed) and the evaluation of solely flood images: in the case of a small flooded area, bimodal histograms are not guaranteed.
Several algorithms or even complete processing chains for automatic flood detection from SAR imagery have been presented to date. In 2009, Martinis et al. in [
1] reported the successful development of a near real-time (NRT) flood detection solution from high-resolution TerraSAR-X data. Split-based histogram thresholding and segmentation-based classification formed the basis of this approach. It was more recently substantially refined and extended as described in [
2]. Very high-resolution X-band SAR images were also used to test a solution developed by utilizing fuzzy logic and auxiliary data in the classification step [
4].
The deployment of the Sentinel-1 mission, delivering a reasonable and uniform revisit acquisition time increased effort in this area. Firstly, Twele et al. in [
3] adapted the solution presented in [
2] for this satellite system. The authors in [
5] tried to improve histogram thresholding using bilateral filtering as a smooth labelling method and tested the approach on three European events. The authors in [
6] used the Otsu iterative thresholding algorithm applied to Sentinel-1 GDR data with single VV (vertical transmit-vertical receive) polarisation. Another approach is mentioned in [
7]. It is based on a normalized difference ratio comparing pre-event and post-event images, using both VV and VH (vertical transmit and horizontal receive) polarisations. Furthermore, more recently, automatic flood detection algorithms were developed from Sentinel-1 imagery based on neural networks [
8,
9], kernel algorithms and machine learning [
10], or alternating decision trees [
11].
Although radar satellite remote sensing (RS) is mainly used to detect open-water flooding in rural environments, some authors dealt with possibilities to use it in urban areas [
12,
13,
14,
15,
16,
17,
18]. Moreover, authors in [
19,
20,
21,
22] worked on the detection of flooded vegetation.
Some authors were applying a fusion of SAR imagery with other types of data, such as optical imagery [
7,
17,
23,
24], LiDAR [
9], WorldDEM [
18], land cover [
25], the shuttle radar topography mission (SRTM) digital elevation model [
25], or flood hazard maps [
14,
26].
For more detailed information about the state-of-the-art of using SAR imagery for flood detection and mapping, the reader is referred to [
27].
Most of the studies mentioned above were performed only to evaluate a specific event or a limited set of events. Only study of Yang et al. was broader, covering the continental part of the USA for three and a half years [
10]. However, even this study was tied to US- specific conditions. The developed system was closely linked to other existing structures in the USA, such as the United States Geological Survey, a network of gauge stations, etc. Some other solutions, e.g., [
6,
7,
8,
9,
10,
11,
13,
14,
17,
18,
24,
26], relied on the use of computational intelligence tools, indexes and thresholding, or on specific auxiliary data sets. The computational intelligence tools used, such as neural networks, must always be trained on sample data first and only afterwards can be used to evaluate the current situation. Their disadvantage is that the training data set binds them to a specific geographical area and spatial or spatiotemporal context. However, they may not work properly in another spatial context. Other methods (e.g., indexes and thresholding) have been mostly tuned to work best in the test area.
The auxiliary data sets used often do not have global coverage, and thus the use of the procedures developed over them is also geographically limited. It is impossible to use them in any place and at any time.
None of the studies mentioned above offer an assessment of the performance of the developed system over a long time period, especially with respect to so-called false positives, i.e., cases when flooding was identified by the system but did not happen in reality. These false positives can happen, e.g., due to agricultural works or seasonal changes.
The solution proposed in this paper is an attempt to overcome these barriers. Our goal was to develop a method that can be applied anytime and anywhere, not dependent on training, using only a minimum auxiliary data which will also have a global validity.
We developed an automated processing chain for NRT open-water flood detection based on Sentinel-1 imagery which can be run in an on-demand or permanent monitoring mode. When the solution is implemented, automatic notifications for a given area of interest (AOI) can be provided to a user without any need for an intervention or manual steps from his side. It can be beneficial mainly for stakeholders who operate infrastructure over vast or distant areas and want to be warned about potential threats or damage to their facilities.
The developed solution is based on a combination of two separate algorithms utilizing Sentinel-1 images to suppress false positives and to enhance the quality of the flood detection results. The first algorithm addresses the detection of areas where a statistically significant change in backscatter intensity happened compared to the previous state. The second algorithm, inspired by [
1,
2,
3], applies thresholding to the backscatter intensity to find open-water flooded pixels and includes fuzzy and region growing approaches in the classification process. Although multi-temporal change detection and intensity thresholding was previously combined by [
28], this is the first time it was performed in an automatic processing context. After a detailed description of the developed solution, the results from processing long periods, including several flood events in Austria, Germany, the Czech Republic and Italy, are demonstrated, and the potential limits of the developed solution are discussed.
4. Experimental Results
Validation results from four European countries, including five flood events, are presented separately in the following sub-sections. As the processing chain was developed for an automatic operation without user intervention in any given AOI, identical settings of flood detection parameters were used for all evaluated events and areas. No fine-tuning of these parameters for specific areas or flood events was therefore applied.
4.1. Germany, Lower Saxony, Summer 2017
An extensive flood event on 25 July–30 July 2017 happened in southern Lower Saxony close to the cities of Hannover, Hildesheim and Braunschweig. This territory is mainly constituted of lowlands. Hydrological report of the flood event is provided in [
30]. Based on the available sources, nine separate sub-areas with a total surface area of 97 km
2 (see
Figure 5) were selected for validation. The localities were selected in order to test the developed solution in a diverse set of conditions. Some of them are constituted mainly by a rural land cover (agricultural fields, meadows); others mainly by urban areas in which the developed solution was not expected to work. The officially reported level of flooding in the selected sub-areas correspond to a return period from 20 to 200 years depending on the locality [
30], but was unfortunately not specified for some of them (n. 3, 4, 5, 8).
Sentinel-1 data from one ascending (A117) and one descending track (D66) were processed with the developed solution from the beginning of March until the end of November 2017. The evaluated flood event was correctly identified in five sub-areas with rural landscape, all from both processed tracks. The flooding was not identified in four urban sub-areas (n. 3, 4, 5, and 6). In terms of false positives, there were only three events in total within all nine sub-areas during the whole processed period. An example of false positive flood detection in sub-area two is shown in
Figure 6. Similar to the other two false-positive events, the majority of the areas wrongly detected as flooded were located on agricultural fields.
The outputs of the developed solution were compared with a reference data set to assess its performance in detecting the extent of real flooding. The reference dataset was based on Sentinel-2 images acquired on days of flooding or shortly after. In sub-areas 3, 4, 5, and 6, no flooding was identified from the Sentinel-2 images acquired on July 29 or on July 31. These are the sub-areas in which the flooding was also not detected by the developed solution.
Figure 7 shows the extent of flooding from both sources for the selected sub-areas 1 and 7. It is apparent from the figure that the flooded polygons detected by the developed solution (in red) correspond well to the flooded polygons obtained from the reference Sentinel-2 imagery (in blue). A similar situation was found in the case of all other sub-areas with detected flooding.
Table 1 contains statistical evaluation between the results of the developed solution and the flooding extents derived from the reference dataset. Overall accuracy (OA) represents a percentage of all pixels classified correctly, producer accuracy (PA) for open water is the probability that an area identified as flooded in reference dataset was flooded in reality. Similarly, user accuracy (UA) for open water represents the probability that an area identified as flooded by the developed solution was flooded in reality. Mean OA reached 97.3%, mean UA (PA) for open water was lower with 77.8% (66.0%) and varied for individual sub-areas. Lower values of PA compared to UA indicate that flooding was detected on a higher number of pixels from reference Sentinel-2 imagery than by the developed solution. The reasons for this situation are discussed in
Section 4.
4.2. Austria, Southern Carinthia, Autumn 2018
From 28 October to 31 October 2018, flooding took place in the southern Carinthia region of Austria, west of Villach city. The floods mainly affected several valleys in the local Alps and consequently areas in lower parts of affected water streams. A general overview of the event is available in the hydrological report [
31]. For validation, the mountain valley lying on the Gail river was chosen for which officially reported flood levels ranged between the 45 and 115 years return periods. Meadows are the prevailing landscape in the selected 54 km
2 area, together with various villages and forests.
Data from one ascending and two descending Sentinel-1 tracks were processed for the period from the beginning of March until the end of November, 2018. The flood event was correctly identified in the selected area on October 29 and November 4 from track A44, on November 2 from track D95 and on November 3 from track D22. A high level of flooding allowed flood detection in some places also a few days after the water level culmination. No false positives occurred during the whole processed period in any processed track.
Similar to the above-presented results for the flood event in Germany, flooded polygons were manually digitized from the Sentinel-2 image acquired on October 31, 2018. The quality of the used Sentinel-2 image was unfortunately not ideal, with overall low contrast and some places covered with low foggy cloudiness. This situation could reduce the quality of flood extent data derived from the image.
The map showing the extents of flooding from both sources is presented in
Figure 8, while a statistical evaluation is found in
Table 2. The UA for open water reached 94.7 %, meaning that almost all pixels identified as flooded by the developed solution were also identified as flooded from the Sentinel-2 image. On the other hand, the PA for open water was only 51.8%, indicating that large areas identified as flooded from the Sentinel-2 image were not detected as flooded from the Sentinel-1 data processing. In this case, the reason for the situation can be at least partly found in the time of acquisition of used images. While Sentinel-2 images were acquired approximately 36 h after the water level culmination (see
Figure 9), one Sentinel-1 image was available shortly after the beginning of water level rising and the second one was from approximately 72 h after the water level culmination. In that time, the extent of open water flooding was already significantly reduced.
4.3. The Czech Republic, Central Moravia, Spring 2020
Severe torrential storms that hit the area north of Uničov city during the night from 7 June to 8 June 2020 resulted in a devastating flash flood. It caused the loss of life of two people and significant property damage (see, e.g.,
https://floodlist.com/europe/czech-republic-flash-floods-olomouc-june-2020, accessed on 1 September 2021). The officially reported level of flooding was equal to a 100-year return period. The course of flood event is provided in
Figure 10, where the water level measured on the Oslava river at the hydrological station located in Dlouhá Loučka is plotted. In the first hours of June 8, a rapid increase of water level was visible, exceeding the third level of flooding activity (marked as SPA in the figure), which corresponds to a level of emergency. Shortly after culmination, the water level started to decrease. Due to the character of the flood event, the water passed quickly through the area and only several temporarily flooded lagoons were created, mostly on agricultural fields around lower parts of water streams in the selected AOI. Own photo documentation was acquired on June 8. The temporarily flooded lagoons disappeared quickly over the next several days.
For the validation of the developed flood detection solution, an area of 21 km
2 was selected. Data from the beginning of February until the end of August 2020 were processed from all four available tracks (A73, A175, D22, D124). In terms of false positives, three events were identified over the whole processed period. All of them occurred in track A175 during March and April. The extent of two of them are visualized in
Figure 11.
While it was not expected to detect the flooding in urban areas lying in the higher parts of water streams, some temporary water lagoons (see
Figure 12) were supposed to be detected. Unfortunately, the flood was not detected in any track. Reasons leading to the failure in the detection of flooded lagoons were therefore analysed and are well understandable from the outputs provided in
Figure 13, where the processing results of the Sentinel-1 image from track D124 acquired on 8 June 2020 are presented. Open water detection in both polarizations (VV, VH) provided correct results, as all four lagoons were detected as flooded. On the other hand, the lagoons were not detected by the change detection algorithm as areas where the backscatter intensity significantly changed compared to the state before the flood event. In this regard, an increased sensitivity of change detection in both polarizations was tested. When applied, lagoons n. 1, 2 and 4 were identified as changed in the VV polarization output, but none of them were identified as changed in the case of the VH output. Areas of all these four lagoons are agricultural fields where vegetation started to grow during the springtime, but a significant part of their surface was still bare soil when the flooding happened. Therefore, the change in backscatter intensity between this type of land cover and the flooded water was not high enough in VH to identify the area as changed, which prevented the overall flood detection.
Although taking only the outputs of the open water detection and neglecting the change detection results would bring good flood detection results in the described case, it cannot be applied generally, or at least not in the case of using the developed solution for autonomous flood monitoring. During validation, many situations were found when some area was wrongly detected as open water flooded, and only the negative result of change detection prevented the overall false positive. On the other hand, if the developed solution is run on demand for specific dates over an area where some flooding occurred, the outputs of open water detection themselves can be used for flood delineation if the change detection blocks the overall detection.
4.4. The Czech Republic, Northeastern Moravia, Autumn 2020
Regional floods hit several river basins in north-eastern and eastern Moravia during 14–16 October 2020 [
33]. They were caused by continuous rainfalls starting on October 13 when previous rains had already saturated the soils. Although floods occurred in several territories, only two of them were selected for the validation because these were the only cases where reference optical imagery without strong cloud coverage was found. The selected two sub-areas are located east of Opava city, in the basin of the Opava river. The landscape is flat there, with prevailing meadows, agricultural fields, wetlands and various permanent water areas. The total surface area of both localities is 62 km
2.
Data from four Sentinel-1 tracks from the beginning of August to the middle of November, 2020 were processed for validation. No false positives occurred in the selected sub-areas during the processed period. For images acquired on October 14 (track A73) and 15 (track A175), floods were successfully detected in both sub-areas. In sub-area 2, floods with a much smaller extent were also detected from images collected on October 18 (track D124) and 19 (track D51). This was mainly due to the character of the selected sub-areas where sub-area 2 contains a complex of protected natural zones where flooding tends to remain longer and could be therefore detected from several consequent radar images.
A multispectral optical image acquired on 15 October 2020 by a PlanetScope satellite was used to create a reference data set. The eastern part of the image was entirely covered with cloudiness and was therefore excluded from the validation process.
Figure 14 shows results from both sources, while
Table 3 contains their statistical comparison. Values of UA for the open water reached 92% in both sub-areas, while values of PA for the open water were 63% and 73%. Therefore, the situation was similar to the one described in previous sections, with a weaker flood detection from the developed solution compared to the reference data set.
Figure 15 and
Figure 16 provide the original Sentinel-1 imagery acquired over both processed sub-areas in track A73 on the day of flooding (14 October 2020) and several days before the flooding (8 October 2020). Dark areas in the images correspond to places with low backscatter intensity, typically represented by open water. The extent of flooding identified by the developed solution from the data acquired in this satellite track on 14 October 2020 are shown in yellow colour. It is well visible that most of dark-coloured areas in the Sentinel-1 imagery from 14 October 2020 were successfully identified as flooded by the developed solution.
4.5. Italy, Piedmont and Lombardy Regions, Autumn 2020
On 2–4 October 2020, extensive floods hit the Piedmont, Lombardy and Liguria regions. Storm Alex brought heavy rains reaching recording rainfalls; values exceeding 500 mm in 24 h were reported at some stations. For validation of the developed solution, a flood event at the river Sesia close to Vercelli city was selected, since it was also processed by the Copernicus Emergency Management Service (EMS), see
https://emergency.copernicus.eu/mapping/list-of-components/EMSR468 (accessed on 31 August 2021). Two casualties and significant losses on the property were reported for this flood event.
Data from a single lowland area of 292 km2 were processed, covering the period from the middle of September until the middle of October in three different Sentinel-1 tracks. Floods were successfully detected in the whole AOI from Sentinel-1 images acquired on October 3 (track A88) and October 4 (track A15). A few small polygons in the southern part of the processed area were also detected as flooded from the images acquired between October 8 and 15 (tracks D66, A15, A88). No false positives occurred during the processed period.
The extent of flood detection from the developed solution was compared with the flood delineation product provided by the Copernicus EMS. Its output was based on a RADARSAT image with a ground sample distance (GSD) of 3.0 m acquired on 6 October 2020 at 17:21 UTC. For a visual evaluation, maps showing the extent of flooding detected by the developed solution and Copernicus EMS are presented in
Figure 17. Since the Sentinel-1 images used by the developed solution were acquired during the flood event while the RADARSAT image used by Copernicus EMS approximately three days after it, significantly larger areas were detected as flooded from the first source. Significantly more extensive flooding in the area than provided by the Copernicus EMS output are also visible in the Sentinel-2 image from October 3 (not shown).
On the other hand, the Copernicus EMS product also provided some flooded polygons in areas that were not detected as flooded by the developed solution. Usually, they were of small size and were therefore problematically detectable from Sentinel-1 imagery due to its coarser spatial resolution compared to RADARSAT. Still, a few larger polygons detected as flooded only by the Copernicus EMS product could be found, e.g., on urban surfaces close to Vercelli city or natural surfaces in the southern part of the processed area.
Table 4 contains the statistical evaluation of the comparison between the outputs of the developed solution and the Copernicus EMS. Values of PA and particularly of UA for open water were very low compared to results presented in previous sections for other flood events. This outcome can be attributed to reasons described in the previous paragraphs, mainly to the time difference in the acquisition of the Sentinel-1 and RADARSAT satellite images. That is why we do not include this locality in the summary of statistics provided in the Conclusion Section.
4.6. False Positives Testing in the South Moravia, 2019
Data from a 27 km
2 rural locality in south Moravia were processed to thoroughly test the developed solution in terms of false positive occurrence. The selected area contains mainly extensively used agricultural fields and some urban areas or vineyards. Sentinel-1 images from three individual tracks spanning from the beginning of February until the end of December 2019 were tested. There were no false-positive alarms found in tracks D22 and D124. In track A73, three epochs with false-positive alarms occurred. All of them happened in April 2019; two of them are presented in the
Figure 18. Since 55 epochs were processed for the A73 track, false positives occurred in 5.5% of epochs. Similar to the results presented in the previous sections, most areas wrongly detected as flooded were located on agricultural fields.
4.7. Flooding Detection from Individual Data Polarizations
As mentioned in
Section 3.1.4, the developed solution combines processing runs in both the VV and VH polarizations for a final result to suppress occurrence of false-positive detections.
Figure 19 provides an overview of the frequency of false positive detections from processing based on individual polarizations and their combinations. From top to bottom, outputs for long-term processing over areas introduced in previous sections (Germany, Austria, central Moravia, south Moravia) are shown. It is evident that the number of final false-positives (red crosses in the upper row far away from the black line) is significantly lower than the number of false detections in either polarization (blue and green crosses in the upper row). In the VH polarization, a larger number of false-positive detections occurred compared to the VV polarization. The frequency of false-positive detections for the Austrian locality was much smaller compared to all three other localities. This is given by the prevailing land use. While the alpine valley in Austria consists mainly of meadows, the other three localities are intensively used for growing agricultural crops. Bare soils and various vegetation cover therefore change throughout the year, impacting the backscatter intensity of the radar signal, leading to false positives.
In
Figure 19, there are several cases where the floods were detected in both polarizations, but not in their combination. It can be caused by two factors: (a) the floods were detected in another processed sub-area (
Figure 19 displays aggregate results), or (b) more frequently, the floods were detected in the same sub-area, but on a different part of the sub-area. It is the most significant benefit of the use of both polarizations. In addition to the instances mentioned above, dates can be found where both the “FLOODS” and “NO FLOODS” states are present for a combination of both polarizations. This means that in some processed sub-areas, the flood was detected, but in some others, it was not.
5. Discussion
As confirmed by the validation, the developed solution can detect only open water pixels and only in non-urban areas due to Sentinel-1 resolution and SAR signal backscatter behaviour in the urban environment. In this regard, processing other SAR satellite data with the developed solution would be straightforward, except that with commercial data with a smaller coverage, it could be more difficult to find the needed permanent water bodies in the same image. We recommend using a permanent water bodies layer for the estimation of water–land thresholds. If the approach based on histogram bimodality is used instead, the reliability of the algorithm gets lower and can even lead to no detection of flooding in areas where the floods were correctly identified while using permanent water bodies in the input.
Generally, the quality of flood detection from SAR depends strongly on the acquisition time of the satellite images with respect to the flood culmination. As shown in
Section 4.2 or
Section 4.5, obtaining an image a few days after water culmination typically leads to a much weaker or even unsuccessful detection. A similar situation applies to flash floods, where the open water does not last long enough to be captured by the satellite. While using the developed solution, areas hit by a flash flood are often visible in the change-detection maps. However, without external information about flood occurrence in the area, they are not distinguishable from other changes such as those caused by common agricultural works. On the other hand, if the information on flood occurrence exists and the AOI covering potentially affected areas is provided, the method could be used in the on-demand mode while using only its change detection results for approximate delimitation of floods extent. A similar approach was described in [
10] where a unique triggering mechanism was developed based on gauge stations and satellite precipitation estimations, identifying potentially flooded zones. Similarly, alerts issued by national meteorological services could be used.
Apart from the acquisition time of SAR imagery, differences between flood detection from the developed solution and reference data sets are supposed to be caused by two reasons.
Firstly, it can be caused by imperfections in open water manual identification from optical satellite imagery due to the varying quality of images and possibilities to mark large polygons as open water flooded, although some vegetation could have been above the water level, changing the radar’s reflection signal and limiting the detection.
Secondly, differences can be caused by weakening open water detection from the developed solution due to using universal parameter settings in automatic processing, as they were tuned to minimize the occurrence of false-positive detections.
In this regard, the results achieved by the developed solution can be compared to the results achieved by automatic SAR flood detectors presented by [
2,
3]. They both reported, for evaluated flood events, UAs for open water ranging from 82.4% to 99.3% and PAs for open water ranging from 83.6% to 98.5%. In [
3], the kappa coefficient reached up to 0.91. While the UA for open water achieved by the developed solution typically exceeded 85%, The PA mostly was only between 60% and 70% and the maximum reached kappa coefficient was 0.77. On the other hand, [
2,
3] did not mention using automated flood detection solutions to process data over a long-term period to evaluate potential false positives. They focused only on the selected flood events themselves. Other authors, such as [
10], do not operate their system permanently, but only in areas and at times when their triggering mechanism reported potential floods. We consider the presented long-term tests to be an essential part of the quality assessment of our developed solution.
The above-mentioned search for universal parameters was necessary as the developed solution is supposed to be run primarily as an automatic continuous monitoring over an area of client’s interest and not on-demand over areas with already proven flooding. As found by performing a set of tests, increasing flood detection sensitivity led to an increasing extent of flood detection, but also to an increasing number of false positives. Increased sensitivity is mainly valid for the change detection algorithm where a statistically significant change in backscatter intensity must happen between the pre-event stack of images and the image evaluated for flooding to allow for flood detection.
6. Conclusions
In this work, we presented a developed automatic flood detection solution based on Sentinel-1 multi-temporal imagery. Although it can be run in on-demand mode, it was prepared mainly for autonomous continuous monitoring of potential flooding in a given area. The solution combines change detection and open-water detection algorithms; both run in the VH and VV polarizations. The only auxiliary data which should be provided by the user are a layer defining AOI and a layer defining permanent water bodies within or near the AOI.
The results of validating the developed solution on five flood events in Europe and an additional false positive test in a rural landscape were provided. Data from long-term periods were processed for at least two satellite tracks and the delivered outputs were compared with prepared reference datasets based on optical or SAR imagery. Floods were successfully detected in four riverine flood events; the method failed in the case of the flash flood event. The classification metrics were computed to statistically compare the extent of floods identified from the developed solution and the reference datasets. The overall accuracy values mostly exceeded 95%. With an exception for the Italian flood event, user accuracy values for the open-water class ranged between 55% and 95%, with a mean of 83%. A high percentage of pixels identified as flooded by the developed solution were also flooded according to the used reference dataset. On the other hand, producer accuracy values for the open-water class were on average only 65%. The developed solution usually provided flooded polygons of a smaller size or did not detect a flood in some places where the reference dataset did. This influenced kappa coefficient values which ranged between 0.59 and 0.77. According to [
34], values between 0.61 and 0.80 correspond a substantial level of agreement. Except for the quality of the prepared reference datasets (influenced by differences in real flooding extent happening between times of acquisition of the Sentinel-1 and reference dataset imagery and imperfections during the manual digitization of reference imagery) the described situation is mainly attributed to the applied tuning of sensitivity parameters to limit false positives. Increasing the sensitivity of flood detection leads to increased flood detection, but also to a higher number of false positives. The flash flood at the beginning of June, 2020 near Uničov city in the Czech Republic was not detected with the developed solution since temporarily flooded water lagoons at agricultural fields were not identified with the change detection algorithm to differ significantly from their pre-event state.
Zero or minimum false positives occurred in all processed events and included sub-areas with various land cover and land use characteristics. All of the false positives were identified during the spring season (from late March till the end of May) and the majority of areas wrongly detected as flooded lay at agricultural land. Therefore, the origin of the false positives is related to changes in agricultural fields during the beginning of the vegetation season and at times of crop harvesting or extensive agricultural work. The developed solution currently does not use any auxiliary data to suppress false positives. In this regard, implementing, e.g., the topo-hydrological factor HAND (height above nearest drainage, [
35]), which is also used in SAR-based flood detection solutions by [
3] or [
5], could be helpful.