Multi-Method Tracking of Monsoon Floods Using Sentinel-1 Imagery

: Floods cause great losses in terms of human life and damages to settlements. Since the exposure is a proxy of the risk, it is essential to track ﬂood evolution. The increasing availability of Synthetic Aperture Radar (SAR) imagery extends ﬂood tracking capabilities because of its all-water and day / night acquisition. In this paper, in order to contribute to a better evaluation of the potential of Sentinel-1 SAR imagery to track ﬂoods, we analyzed a multi-pulse ﬂood caused by a typhoon in the Camarines Sur Province of Philippines between the end of 2018 and the beginning of 2019. Multiple simple classiﬁcation methods were used to track the spatial and temporal evolution of the ﬂooded area. Our analysis indicates that Valley Emphasis based manual threshold identiﬁcation, Otsu methodology, and K-Means Clustering have the potential to be used for tracking large and long-lasting ﬂoods, providing similar results. Because of its simplicity, the K-Means Clustering algorithm has the potential to be used in fully automated operational ﬂood monitoring, also because of its good performance in terms of computation time.


Introduction
Floods are among the most frequent and widespread natural hazards in the world. Being related to intense and/or extreme weather events, they cause great losses in terms of human life and damage to commercial and productive sites, infrastructures, and agriculture [1,2]. Particularly, it has been estimated that floods are responsible for approximately 40% of the total damage caused by natural hazards [3,4]. Exposure to flooding is considered a proxy of the risk [5], so that an evaluation of the extent of potentially inundated areas is crucial for hazard and risk assessment and represents the basis for land planning and policy decisions oriented towards flood mitigation (i.e., occupation restrictions, recommended uses, and flood insurance plan development). Flood hazard and risk evaluations can be completed using the results of deterministic hydrodynamic models that simulate water movement across the floodplain. In the presence of monitoring (fluvial stage or discharge) and topographic data, statistical models associated with GIS processing can provide a basis for such kinds of analyses [6][7][8]. Both deterministic and statistical models need to be calibrated and validated using available flood data in terms of spatial extent, persistence, and frequency.
Data derived from a number of satellite platforms can be used to image floods, providing a basis for a rapid and effective response to natural disasters. Satellite observations have the advantage of covering large areas with an increasingly short revisiting time, which makes them able to support continuous observation and operation monitoring of floods [9,10]. Among different satellite data, the central section of the Bicol Peninsula and is the largest province in the Bicol Region. The Bicol River is the main watercourse of the province and is surrounded by Mount Isarog (1966 m asl) and Mount Iriga (1196 m asl). The eastern portion of the province lies on the mountainous peninsula of Caramoan, which faces the island of Catanduanes to the east. The Bicol River drains the central and southern parts of the province toward the San Miguel Bay. Mount Iriga is surrounded by three lakes named Buhi, Bato, and Baao. The Albay Province is generally mountainous with scattered valleys. On the eastern part of the province is a line of volcanic mountains starting with the northernmost Malinao in Tiwi and followed by Mount Masaraga (1328 m asl) and the Mayon Volcano (2463 m asl).
valleys. On the eastern part of the province is a line of volcanic mountains starting with the northernmost Malinao in Tiwi and followed by Mount Masaraga (1328 m asl) and the Mayon Volcano (2463 m asl).
The Philippines archipelago has a tropical and strongly monsoonal climate characterized by relatively high temperature, oppressive humidity, and plenty of rainfall. Two main seasons alternate for the majority of islands, a wet season dominated by monsoon rainfall and a dry season. Some areas experience rainfall throughout the year, and clear season alternations are absent. Four climate regions are recognizable across the Philippines archipelago characterized by different distributions of rainfall [36,37]. The study area spans across three of these four climatic regions ( Figure 1): region II is characterized by the absence of a dry season with a very pronounced maximum rainfall during the months of November and December, region III is characterized by a slightly drier seasons between November and April and a wet season between May and October, and region IV is characterized by rainfall distributed throughout the entire year. From June to December, tropical cyclones (typhoons) often strike the Philippines. Most of these storms come from the southeast and are heaviest in Samar, Leyte, south-central Luzon, and the Batan Islands, and, when accompanied by floods or high winds, they may cause great loss of life and properties. In December 2018, as reported by the Department of Science and Technology (PAGASA) of the Republic of Philippines (http://bagong.pagasa.dost.gov.ph/tropical-cyclone), the study area was flooded as a consequence of a typhoon induced by a tropical depression (named "Usman"). This was the 21st and the last cyclone that intersected the Philippines in 2018. The "Usman" depression, developed from tropical disturbances, resulted in widespread heavy rainfall over large portions of Southern Luzon and Eastern Visayas. From 28 to 29 December 2018, prolonged and excessive rainfall over these areas resulted in multiple landslides and floods. Of the 139 stations with reliable data, the highest 2-day rainfall was estimated at 573.2 mm and was recorded in Daet (Camarines Norte), 50 km The Philippines archipelago has a tropical and strongly monsoonal climate characterized by relatively high temperature, oppressive humidity, and plenty of rainfall. Two main seasons alternate for the majority of islands, a wet season dominated by monsoon rainfall and a dry season. Some areas experience rainfall throughout the year, and clear season alternations are absent. Four climate regions are recognizable across the Philippines archipelago characterized by different distributions of rainfall [36,37]. The study area spans across three of these four climatic regions ( Figure 1): region II is characterized by the absence of a dry season with a very pronounced maximum rainfall during the months of November and December, region III is characterized by a slightly drier seasons between November and April and a wet season between May and October, and region IV is characterized by rainfall distributed throughout the entire year. From June to December, tropical cyclones (typhoons) often strike the Philippines. Most of these storms come from the southeast and are heaviest in Samar, Leyte, south-central Luzon, and the Batan Islands, and, when accompanied by floods or high winds, they may cause great loss of life and properties.
In December 2018, as reported by the Department of Science and Technology (PAGASA) of the Republic of Philippines (http://bagong.pagasa.dost.gov.ph/tropical-cyclone), the study area was flooded as a consequence of a typhoon induced by a tropical depression (named "Usman"). This was the 21st and the last cyclone that intersected the Philippines in 2018. The "Usman" depression, developed from tropical disturbances, resulted in widespread heavy rainfall over large portions of Southern Luzon and Eastern Visayas. From 28 to 29 December 2018, prolonged and excessive rainfall over these areas resulted in multiple landslides and floods. Of the 139 stations with reliable data, the highest 2-day rainfall was estimated at 573.2 mm and was recorded in Daet (Camarines Norte), 50 km away from the study area. In the week preceding the event, rainfall accumulations in excess of 100 mm were already observed by multiple rain gauges over the eastern section of Luzon and Visayas, especially over Bicol Region mainland, Samar island, Aurora, Quezon, and Rizal. The event caused hundreds of casualties and millions of dollars in damages to personal and public property. According to the report of the Department of Social Welfare and Development (DSWD), more than 680,000 people were affected by the storm. Nearly 55% were located primarily in the province of Camarines Sur in the Bicol Region. Over 1900 homes were destroyed, and more than 15,000 were damaged. Agricultural losses were estimated at PHP 2 billion (US $37 million), with over 56,000 farmers and fishers affected, according to the United Nations Food and Agriculture Organization.

Data
To quantify the extent of the flooded area, i.e., the area covered by open water, and track its spatial and temporal evolution, we used multiple Sentinel-1 A/B SAR images. The satellite data were downloaded from the Copernicus Open Data Access Hub (https://scihub.copernicus.eu/dhus/#/home). Sentinel-1 data products consist of imagery with a medium resolution (10-20 m) that represents the radar echo from the surface of the Earth of the signal emitted by the onboard antenna of the satellite. The Sentinel-1 satellite constellation is formed by two satellites (S-1 A and S-1 B) characterized by an acquisition frequency of six days. Each satellite is equipped with a 5.405 GHz C-band (λ ≈ 5.6 cm) imager payload (CSAR). The CSAR instrument supports operation in dual-polarization (HH + HV and VV + VH). The Sentinel-1 satellites can operate in three imaging modes for various observation approaches, spatial resolutions, and swath widths [38]. In this study, we used a total of 24 images with VV polarization. We chose this polarization because it characterized all of the selected data; only a limited number of images were acquired in double polarization VV + HH. VV polarization, in some cases, provides better results in detecting open water features [39,40]. The selected images had the same angle of incidence (38.8 • ) and footprint and were acquired in descending mode. The selected images covered the time period between 13 November 2018, prior to the typhoon, and 18 March 2019, after the typhoon. The first selected image was acquired 78 days before the main flooding event, occurring between 28 and 29 December 2018. This long pre-event time span allowed us to observe the normal fluctuations of the open water area before the flood.

Data Pre-Processing
Waterbody identification and extraction were completed considering the difference in the backscattering coefficient (Gamma γ • ) between open water and the land surface. This difference is related to the roughness of the target surface, which is regulated by the Rayleigh criterion. Since the open water theoretically has no roughness, the incident radiation is reflected away from the sensor; thus, the backscattering coefficient related to its flat surface is lower than the land surface [41][42][43]. In contrast, the land surface that has a variable roughness has a variable backscattering coefficient typically higher than that characterizing open water. SAR imagery, downloaded by Sentinel Data Hub, were pre-processed using the Sentinel Application Platform (SNAP) toolkit distributed by ESA (https://step.esa.int/main/toolboxes/snap/). The products were first cropped across the area of interest to allow for shorter processing times. Subsequently, radiometric calibration was applied, and terrain was corrected to the Gamma band (i.e., backscatter coefficient γ • ) [44]. To complete this last step, the STRM 1 arc-second digital elevation model (DEM) was used. Finally, in order to obtain high-quality SAR images and reduce the noise, we applied the Lee-Sigma filter [45,46] characterized by a combination of 5 × 5 and 7 × 7 windows.
As the last step of image pre-processing, calibrated and filtered SAR images were cropped over the land of the study area. This was done in order to remove (from images and relative histograms) the changing backscattering effects induced at the ocean surface by the difference in wind conditions that generating waves mimicked the surface roughness, which produced a double-bounce effect [47]. After pre-processing, the rearranged backscatter intensities were converted into dB (decibels). Figures 2 and 3 show a selection of 12 SAR images, out of 24 used, in which it is possible to visually identify the different flood traces of the events.      Figure 4 summarizes the workflow used to process images and quantify the total water area in the reference time period. In order to identify the open water features to determine the total area covered by water and its spatial and temporal evolution due to monsoon flood events, three different classification methods were applied that might have the potential to be used in operational monitoring.

Extraction of the Water Body Area
Our multi-method approach allowed us to compare specific results and quantify relative differences (e.g., [48,49]) in order to identify the best method to be used from an operational monitoring perspective. Particularly, we used two methods based on threshold identification and a third based on unsupervised K-Means Cluster analysis. Threshold-based methods are widely used because of their efficiency in identifying flooded areas in SAR images [50][51][52][53]. These methods consisted of defining a fixed value for the backscattering coefficient (i.e., intensity) that split the histogram in two subsets (clusters), one subset representing the land surface and the other representing the open water. Pixels with a low grayscale intensity typically corresponded to water bodies, while pixels with high values corresponded to the land (in our case, background). There are several techniques (supervised and unsupervised) to find the optimal threshold value [54][55][56][57]. Among these, for our comparative analysis, we selected visual identification of the threshold guided by the Valley Emphasis criterion, the Otsu method, and the unsupervised classification method based on K-Means Clustering. Figure 4 summarizes the workflow used to process images and quantify the total water area in the reference time period. In order to identify the open water features to determine the total area covered by water and its spatial and temporal evolution due to monsoon flood events, three different classification methods were applied that might have the potential to be used in operational monitoring. Our multi-method approach allowed us to compare specific results and quantify relative differences (e.g., [48,49]) in order to identify the best method to be used from an operational monitoring perspective. Particularly, we used two methods based on threshold identification and a third based on unsupervised K-Means Cluster analysis. Threshold-based methods are widely used because of their efficiency in identifying flooded areas in SAR images [50][51][52][53]. These methods consisted of defining a fixed value for the backscattering coefficient (i.e., intensity) that split the histogram in two subsets (clusters), one subset representing the land surface and the other representing the open water. Pixels with a low grayscale intensity typically corresponded to water bodies, while pixels with high values corresponded to the land (in our case, background). There are several techniques (supervised and unsupervised) to find the optimal threshold value [54][55][56][57]. Among these, for our comparative analysis, we selected visual identification of the threshold guided by the Valley Emphasis criterion, the Otsu method, and the unsupervised classification method based on K-Means Clustering. The first method consisted of visually inspecting a grayscale histogram of each SAR image and applying the Valley Emphasis criterion [58][59][60]. For the application of this criterion, we first analyzed the grayscale histogram in SNAP and, subsequently, chose the minimum backscatter intensity value that corresponded to the valley between the two peaks, in the case of bimodal distribution, or to the bottom rim of a single peak for a unimodal distribution. Once identified, we converted this intensity value from dB to a "grayscale" value between 0 to 255 (8-bit images). The identified threshold value was then used as a basis for manual image classification. Image classification was consistently The first method consisted of visually inspecting a grayscale histogram of each SAR image and applying the Valley Emphasis criterion [58][59][60]. For the application of this criterion, we first analyzed the grayscale histogram in SNAP and, subsequently, chose the minimum backscatter intensity value that corresponded to the valley between the two peaks, in the case of bimodal distribution, or to the bottom rim of a single peak for a unimodal distribution. Once identified, we converted this intensity value from dB to a "grayscale" value between 0 to 255 (8-bit images). The identified threshold value was then used as a basis for manual image classification. Image classification was consistently completed in the GIS environment (QGIS 3.7), assigning 0 to pixels with values higher than the threshold (i.e., land) and 1 to pixels with values lower than the threshold (i.e., water). After classifying the images, and considering that low backscattering intensity values can be also related to the presence of shadow zones (i.e., not only to the presence of water), especially in presence of complex topography, we used a further classification method based on topographic attributes to identify and exclude these areas from our analysis. Notably, we generated a slope map using the available STRM 1 arc-second DEM and used a slope threshold of 3 • to exclude the low backscattering areas with slopes higher than the threshold.

Extraction of the Water Body Area
The second method used to find the correct threshold value was the Otsu method. This method is widely used in applications where it is necessary to split grayscale images into two different classes. This method automatically chooses the threshold value from the grayscale histogram on the basis of the minimum within-class variance or the maximum between-class variance [61]. The Otsu method is commonly used because of its simplicity, and the best results are obtained if the image histogram is characterized by a bimodal or multimodal distribution. In the presence of a unimodal distribution, this method fails. The Otsu method for finding an optimal threshold value for image classification was applied in the Matlab™ environment. The Otsu threshold was then used to classify the images with the same procedure used as the previous method.
The last method was the K-Means Cluster analysis. This clustering method is a statistical technique widely used for dimensional reduction and has the potential to be used for unsupervised analysis of SAR images [62,63]. In our case, K-Means Clustering was used for dividing the image in two classes (i.e., water and land clusters). This method groups pixels on the basis of their grayscale level distribution and standard deviation. The purpose of this algorithm is to reduce variability within clusters, and the objective function is the sum of square distances between cluster centers and its assigned pixel value [64]. This clustering method, implemented in the SNAP toolbox, was applied considering two clusters, 30 interactions, and 32,000 random seeds. After cluster identification, the topographic filter and subsequent manual classification of the images were completed. Finally, we estimated the total water area in each image.
Once the spatial extent of open water features across the study area was evaluated, to further analyze the spatial and temporal distribution of flood intensity and persistence, and its change between the different methods, we estimated the submergence ratio across the whole of the study area and the reference period, and we compared this ratio in terms of arithmetic differences between products derived by our data (i.e., water area extent) obtained by different methods. A submergence ratio equal to 1 indicated the constant presence of water in the reference area, while a value of 0 indicated the constant absence of water in the reference area and across the reference period. A pixel value between 1 and 0 indicated that the reference area was only temporally submerged. This condition is consistent with the occurrence of the flood event, and the difference in the submergence ratio is representative of flood persistence and might be related to the water level. The submergence ratio map was constructed on the basis of the extent of open water features derived using the K-Means Clustering classification method. To determine the variation in the submergence ratio derived by the use of open water coverage derived by the three classification methods, we calculated the submerged ratio difference by iterative subtraction. Figures 5 and 6 show several examples of results in terms of open water identification from the different methods used to select the optimal threshold value of backscattering intensity. The first example reported in Figure 5a,c,e represents the water coverage at the beginning of the monitoring period before the occurrence of the rainfall events responsible for the flooding. The second example reported in Figure 5b,d,f represents the water coverage immediately after the first flood pulse induced by rainfall. Figure 6a,c,e represent the water coverage immediately after the second flood pulse induced by rainfall. The second example reported in Figure 6b,d,f represents the water coverage immediately after the third flood pulse induced by rainfall. For comparison purposes, in each map of Figures 5 and 6, polygons representing water bodies as depicted in the SWBD dataset (red polygon, e.g., [65]) were overlaid with areas classified as water bodies (blue area). The SWBD is the worldwide open water bodies boundary in vector format, which were generated by the National Geospatial Intelligence Agency [66]  (red polygon, e.g., [65]) were overlaid with areas classified as water bodies (blue area). The SWBD is the worldwide open water bodies boundary in vector format, which were generated by the National Geospatial Intelligence Agency [66] and represent the maximum extent of open water features. A comparison of the results in terms of water coverage indicated the consistency of all of the methods across most of the study area, with local relative underestimation of the Valley Emphasis method. An example is in the northeastern sector of the study area. No differences between the Otsu and K-Means Clustering were appreciable at the scale of Figures 5 and 6.    Figure 7 shows the evolution of the open water features during the period of interest, including the rainfall events responsible for multiple flood pulses, in terms of total extent for each considered method. Particularly, the red, blue, and green lines represent the total flooded area extracted by Valley Emphasis, Otsu, and K-Means Clustering methods, respectively. All three methods showed  Figure 7 shows the evolution of the open water features during the period of interest, including the rainfall events responsible for multiple flood pulses, in terms of total extent for each considered method. Particularly, the red, blue, and green lines represent the total flooded area extracted by Valley Emphasis, Otsu, and K-Means Clustering methods, respectively. All three methods showed the same pattern of water coverage over time with some variable differences in terms of total estimation. Before the first flood pulse, from 13 October to 24 December the graph indicated a slow increase of water coverage, from a minimum coverage area of about 52 km 2 (grey dashed line in the graph). After this first period, three individual flood pulses with decreasing intensities were identified in terms of total water area increase by our multi-method analysis. According to the PAGASA report, the first main flood pulse occurred between 30   Even if the three methods showed consistent trends in time, a number of differences in the total estimated water area with the different methods are depicted in the graph of Figure 7. The purple line of the graph provides an overview of this change in terms of the maximum percentage change. Indeed, the average difference of approximately 20% was observable. Because of the different estimations of the open water area at the beginning of the period of interest and the absence of an official estimation, our analysis was able to provide only the relative change of the open water area coverage estimated for each different flood pulse occurring in the reference period. Notably, our analysis indicated that the first (main) flood pulse produced an increase in total water area between 245 and 326 km 2 , the second pulse was responsible for an increase between 135 and 186 km 2 , and the third event induced an increase between 45 and 88 km 2 . As described above and observable from the graph of Figure 7, the minimum value of open water area coverage was consistently derived by the Valley Emphasis method, while the maximum value was consistently derived by the Otsu method. The K-Means Clustering method showed water coverage values between the minimum (i.e., Valley Emphasis method) and the maximum (i.e., Otsu method). The flood-induced increases in the open water area estimated with K-Means Clustering were equal to 300, 156, and 68 km 2 for the first, second, Even if the three methods showed consistent trends in time, a number of differences in the total estimated water area with the different methods are depicted in the graph of Figure 7. The purple line of the graph provides an overview of this change in terms of the maximum percentage change. Indeed, the average difference of approximately 20% was observable. Because of the different estimations of the open water area at the beginning of the period of interest and the absence of an official estimation, our analysis was able to provide only the relative change of the open water area coverage estimated for each different flood pulse occurring in the reference period. Notably, our analysis indicated that the first (main) flood pulse produced an increase in total water area between 245 and 326 km 2 , the second pulse was responsible for an increase between 135 and 186 km 2 , and the third event induced an increase between 45 and 88 km 2 . As described above and observable from the graph of Figure 7, the minimum value of open water area coverage was consistently derived by the Valley Emphasis method, while the maximum value was consistently derived by the Otsu method. The K-Means Clustering method showed water coverage values between the minimum (i.e., Valley Emphasis method) and the maximum (i.e., Otsu method). The flood-induced increases in the open water area estimated with K-Means Clustering were equal to 300, 156, and 68 km 2 for the first, second, and third flood pulses, respectively. The shift from the other two methods ranged, in relative terms and considering the entire reference period, between 4% and 31%.

Results and Discussion
The Otsu method consistently showed higher open water surface values across the reference period. The tendency to overestimate water coverage by the Otsu classification method in comparison with the Valley Emphasis method was observed also by Ba Duy [59], who confirmed that this method outperformed the Otsu method in terms of absolute identification capabilities. In addition, the author of the study underlines the limitation of the Valley Emphasis method for extracting water by SAR images, including inefficient identification of mixed water pixels, confusion of water bodies with background noise, and difficulty in choosing an optimal threshold value for very large areas. Bin Cui et al. [28] indicated that conversely from threshold-based methods like the Otsu and Valley Emphasis, the unsupervised methods, like the K-Means Clustering, have great potential in the context of image classification.
The good performance of the unsupervised method in comparison with threshold-based methods can be related to its ability in capturing the spatial variability of the backscattering coefficient that represents a basis for a self-adaptation to different environmental conditions captured by SAR images. Our application confirmed that despite the lack of supervision and parameter estimation (e.g., threshold), the K-Means Clustering performed at least as good as the other methods, and possibly better than the Otsu. In this way, our interpretation is that the unsupervised methods, like the K-Means Clustering, have the potential to support SAR imagery-based flood mapping, also in long-term operational frameworks. This is also related to better results in terms of computational time [28]. It is interesting to mention that additional methods for water body classification, like Bimodal Histogram and Local Adaptive Thresholding [67], have been introduced and have the potential to support this kind of analysis. However, their applicability and reliability need to be further analyzed. Figure 8a shows the persistence of the water in the reference period across the whole study area. In the submergence map of Figure 8a, permanent water bodies like the major lakes and the river are easily recognizable, as they are marked by a color corresponding to a value of the submergence ratio equal to 1. The Baao and Bato lakes expanded due to the flood events until they reached the boundary indicated in the SWBD dataset for a limited period of time (i.e., 20% of the reference time period). Most of the area hit by the flood stayed submerged between 5% and 25% of the time, with a number of spots in which the water persisted for approximately 40% of the reference time period. These spots were in the areas surrounding major lakes and next to the northern coastline of the study area. We expect that the temporal persistence of the flood was directly related to the water level and the presence of a draining channel. Figure 8b-d, report the difference in the submergence ratios calculated with products derived by image classification and the selected methods. Notably, the most important differences were between submergence ratios calculated with the products derived by K-Means Clustering and Valley Emphasis based threshold definitions and between Otsu and Valley Emphasis based threshold definitions. This is consistent with the results of our analysis summarized in Figure 7.

Conclusions
Our analysis indicates that Sentinel-1 data have the potential to be used in flood tracking and for monsoon events characterized by long-term rainfall associated with persistent cloud coverage, which excludes the possibility of using the other types of satellite products, like optical imagery. In this study, the multi-pulse flood caused by the Usman tropical depression in the Camarines Sur Province of the Philippines was analyzed with the aim to provide quantitative flooding data and test multiple methods from a relative perspective connected to the absence of accessible official data. All of the pulses characterizing the flood event occurred between 24 December 2018, and 23 January 2019, remembered for its dramatic effects over the Camarines Sur Province, and they were detected and characterized in terms of the open total flooded area. Since the discrimination between water features and the land surface might be a very challenging task in terms of both accuracy and computation time, in our study, the open water features extraction was completed with three different methods: the Valley Emphasis based manual threshold identification, the Otsu method, and the K-Means Clustering method. All of the methods identified a shape consistent with the flooded area, with slight differences in dimensions (i.e., total area) observed over the monitoring period. These differences seem to be not correlated to the size of the total area covered by open water. Although an official estimation of the flooded area is not available, these methods provide a total flooded area that Figure 8. The map of (a) shows the submergence ratio (i.e., the flood persistence) in the reference period. A pixel value equal to 1 corresponds to an always flooded area. A value of 0 corresponds to a never flooded area. The maps of (b-d) are the differences in submergence ratios calculated using open water features extracted by different methods.

Conclusions
Our analysis indicates that Sentinel-1 data have the potential to be used in flood tracking and for monsoon events characterized by long-term rainfall associated with persistent cloud coverage, which excludes the possibility of using the other types of satellite products, like optical imagery. In this study, the multi-pulse flood caused by the Usman tropical depression in the Camarines Sur Province of the Philippines was analyzed with the aim to provide quantitative flooding data and test multiple methods from a relative perspective connected to the absence of accessible official data. All of the pulses characterizing the flood event occurred between 24 December 2018, and 23 January 2019, remembered for its dramatic effects over the Camarines Sur Province, and they were detected and characterized in terms of the open total flooded area. Since the discrimination between water features and the land surface might be a very challenging task in terms of both accuracy and computation time, in our study, the open water features extraction was completed with three different methods: the Valley Emphasis based manual threshold identification, the Otsu method, and the K-Means Clustering method. All of the methods identified a shape consistent with the flooded area, with slight differences in dimensions (i.e., total area) observed over the monitoring period. These differences seem to be not correlated to the size of the total area covered by open water. Although an official estimation of the flooded area is not available, these methods provide a total flooded area that approximately averaged the results obtained with the Valley Emphasis based manual threshold identification and the Otsu method. The K-Means Clustering method is advantageous because (i) it is unsupervised, (ii) it is fully automatic, and (iii) it returned better results in terms of the computation time. A fast-computational time is a great advantage when it is necessary to have a quick overview of a flood situation. On this basis, this method might have the potential to be used in an operational flood monitoring perspective. In this framework, SAR images are an important means to study flood phenomena because of their day and night operational capabilities and their increasing frequency of acquisition related to the number of satellite constellations. In addition, these capabilities are enhanced by the free availability of data, like those derived by Sentinel-1 satellites.