Difference between revisions of "ALGORITHM2"

From MyDewetra World
Jump to: navigation, search
(Created page with "caption [Home] - [Algorithms] ---- '''HASARD''' is a probabilistic flood mapping algorithm developed by [https:/...")
 
 
(30 intermediate revisions by the same user not shown)
Line 3: Line 3:
 
----
 
----
  
'''HASARD''' is a probabilistic flood mapping algorithm developed by [https://www.list.lu LIST] based on [https://en.wikipedia.org/wiki/Bayesian_inference Bayesian inference]. <br>
+
[https://www.dlr.de/ DLR]'s flood detection algorithm (provisionally, '''Algorithm2''') is a [https://it.wikipedia.org/wiki/Fuzzy_Logic fuzzy logic]-based water class membership assignment, described in Martinis et al. (2015)<sup>[1]</sup> and Twele et al. (2016)<sup>[2]</sup>. <br>
  
Following the procedure introduced in Giustarini et al. (2016) <sup>[1]</sup> the probability of each pixel in any newly acquired SAR image being flooded given its measured backscatter value is estimated in near-real time.<br>
+
It aims to exclude water-lookalikes and to reduce underestimations from initial classification by constructing a fuzzy set that consists of:
The approach is based on the two distribution functions of backscatter values attributed to flooded and non-flooded terrain in the area covered by the image. <br>
+
*(a) the backscatter level
The procedure to calibrate the four parameters of the [https://en.wikipedia.org/wiki/Gaussian_function#:~:text=The%20graph%20of%20a%20Gaussian,symmetric%20%22bell%20curve%22%20shape.&text=Gaussian%20functions%20are%20often%20used,variance%20%CF%83%202%20%3D%20c%202. Gaussian probability distribution functions (PDFs)] based on the histograms of every newly acquired SAR image is explained in detail in Chini et al. (2017)<sup>[2]</sup>. <br>
+
*(b) the elevation of an image element in comparison to the mean elevation of the initially derived water areas
 +
*(c) topographic slope information
 +
*(d) the size of an individual flood object; degree of an element’s membership to the class water is determined by standard S and Z membership functions.
  
When the parameters of the PDFs are known, the conditional probability of a pixel being flooded given its backscatter value can be derived using [https://en.wikipedia.org/wiki/Bayes%27_theorem Bayes’ theorem].<br>
 
Its computation, thus, requires both the probability distributions of backscatter values of flooded and non-flooded pixels and the prior probabilities of a pixel being flooded and non-flooded. <br>
 
In the Bayesian framework, the prior probability ''p(F)'' of a pixel being flooded corresponds to the probability of it being flooded prior to the backscatter value being measured. <br>
 
In principle, for a given pixel, any information other than the backscatter value of the pixel itself, ''σ0'', can be used to estimate these prior probabilities. <br>
 
However, when no such ancillary data is available, the simplest choice is to use a non-informative prior, where ''p(F) = 0.5'', and where flooded and non-flooded pixels are equally likely a priori. <br>
 
  
Findings in Giustarini et al. (2016)<sup>[1]</sup> showed that the non-informative prior leads to acceptable performance in terms of reliability and little is to be gained by attempting to estimate prior probabilities from other sources of data.<br>
+
The average of the individual membership degrees is computed for each pixel in order to combine all fuzzy elements into a single composite fuzzy set (see figure below).
The non-informative priori is, thus, considered an acceptable choice, especially in the absence of any other exploitable data. It also worth mentioning that [https://en.wikipedia.org/wiki/Bayesian_inference Bayesian inference] is applied independently to each pixel, regardless of its spatial location.<br>
+
 
Preliminary masking with the exclusion layer and the reference water extent layer is undertaken so that no probabilities are computed for those areas. <br>
+
[[File:alg2.jpg|centre|800px|border]]
Probability map values range from 0 and 1.
+
''The automatic tile-based thresholding procedure is applied to a given SAR image to select for representative tiles therein. A sub-histogram is generated for a selected tile, based on the backscatter values, one of the parameters of the fuzzy set (left); the corresponding fuzzy standard Z function for backscatter is shown (middle), along with the resulting fuzzy mask and classification maps (right; adapted from Martinis et al. (2015)<sup>[1]</sup>).''
 +
 
 +
The membership degree strongly depends on the position of the crossover point xc (i.e. the half of the width of the fuzzy curve), which is defined by the fuzzy thresholds ''x1'' and ''x2''.The fuzzy threshold values for each element are either determined according to statistical computations or are set empirically. The membership degree of the subsequently, the flood mask is derived through a threshold defuzzification step, which transforms each image element with a membership degree > ''0.6'' into a discrete thematic class. Furthermore, the region-growing tolerance criterion is defined by a relaxed fuzzy threshold of > ''0.45'' based on the composite fuzzy set.
 +
Uncertainty information is contained in the fuzzy mask generated over classified water pixels, which takes on values ranging from ''0'' and ''1''.
 +
 
 +
 
 +
====<u>''Workflow of the algorithm''</u>====
 +
[[File:alg2_wf.jpg|centre|800px|border]]
 +
 
 +
 
 +
====<u>''In a nutshell''</u>====
 +
The key strength of '''Algorithm 2''' is the automatic identification of flooded areas in the SAR data using '''hierarchical tile-based thresholding''' and the optimization of the classification by combining various information sources using '''fuzzy-logic theory''' and '''region growing'''.
  
  
Line 24: Line 32:
  
  
----
 
 
====<u>''References''</u>====
 
====<u>''References''</u>====
<sup>[1]</sup> Giustarini, L., Hostache, R., Kavetski, D., Chini, M., Corato, G., Schlaffer, S., and Matgen P. (2016). Probabilistic flood mapping using Synthetic Aperture Radar data, IEEE Transactions on Geoscience and Remote Sensing, 54 (12), 6958-6969 <br>
+
<sup>[1]</sup> [https://www.sciencedirect.com/science/article/pii/S0924271614001981 <u>Martinis, S., Twele, A., Kersten, J., (2015). A fully automated TerraSAR-X based flood service. ISPRS Journal of Photogrammetry and Remote Sensing, 104, 203-212.</u>] <br>
[https://ieeexplore.ieee.org/abstract/document/7548327 <u>download Giustarini et al. (2016)</u>]
 
  
  
<sup>[2]</sup> Chini M., Hostache R., Giustarini L., Matgen P., A Hierarchical Split-Based Approach for parametric thresholding of SAR images: flood inundation as a test case, IEEE Transactions on Geoscience and Remote Sensing, 55 (12), 6975-6988, 2017.<br>
+
<sup>[2]</sup> [https://www.tandfonline.com/doi/abs/10.1080/01431161.2016.1192304 <u>Twele, A., Cao, W., Plank, S., Martinis, S., 2016. Sentinel-1 based flood mapping: a fully automated processing chain. International Journal of Remote Sensing, 37 (13), 2990-3004.</u>]
[https://ieeexplore.ieee.org/document/8017436 <u>download Chini et al. (2017)</u>]
 
  
  
 
----
 
----
 
[[GFMS | [Home]]] - [[ALGORITHMS | [Algorithms]]]
 
[[GFMS | [Home]]] - [[ALGORITHMS | [Algorithms]]]

Latest revision as of 11:11, 10 March 2021

caption

[Home] - [Algorithms]


DLR's flood detection algorithm (provisionally, Algorithm2) is a fuzzy logic-based water class membership assignment, described in Martinis et al. (2015)[1] and Twele et al. (2016)[2].

It aims to exclude water-lookalikes and to reduce underestimations from initial classification by constructing a fuzzy set that consists of:

  • (a) the backscatter level
  • (b) the elevation of an image element in comparison to the mean elevation of the initially derived water areas
  • (c) topographic slope information
  • (d) the size of an individual flood object; degree of an element’s membership to the class water is determined by standard S and Z membership functions.


The average of the individual membership degrees is computed for each pixel in order to combine all fuzzy elements into a single composite fuzzy set (see figure below).

Alg2.jpg
The automatic tile-based thresholding procedure is applied to a given SAR image to select for representative tiles therein. A sub-histogram is generated for a selected tile, based on the backscatter values, one of the parameters of the fuzzy set (left); the corresponding fuzzy standard Z function for backscatter is shown (middle), along with the resulting fuzzy mask and classification maps (right; adapted from Martinis et al. (2015)[1]).

The membership degree strongly depends on the position of the crossover point xc (i.e. the half of the width of the fuzzy curve), which is defined by the fuzzy thresholds x1 and x2.The fuzzy threshold values for each element are either determined according to statistical computations or are set empirically. The membership degree of the subsequently, the flood mask is derived through a threshold defuzzification step, which transforms each image element with a membership degree > 0.6 into a discrete thematic class. Furthermore, the region-growing tolerance criterion is defined by a relaxed fuzzy threshold of > 0.45 based on the composite fuzzy set. Uncertainty information is contained in the fuzzy mask generated over classified water pixels, which takes on values ranging from 0 and 1.


Workflow of the algorithm

Alg2 wf.jpg


In a nutshell

The key strength of Algorithm 2 is the automatic identification of flooded areas in the SAR data using hierarchical tile-based thresholding and the optimization of the classification by combining various information sources using fuzzy-logic theory and region growing.


For further details the Reader is referred to the dedicated section of the Product Description Document: https://www.gfm_pdd.org/Algorithm2


References

[1] Martinis, S., Twele, A., Kersten, J., (2015). A fully automated TerraSAR-X based flood service. ISPRS Journal of Photogrammetry and Remote Sensing, 104, 203-212.


[2] Twele, A., Cao, W., Plank, S., Martinis, S., 2016. Sentinel-1 based flood mapping: a fully automated processing chain. International Journal of Remote Sensing, 37 (13), 2990-3004.



[Home] - [Algorithms]