Skip to content

Commit

Permalink
Update bayts.rst
Browse files Browse the repository at this point in the history
Proofreading complete.
  • Loading branch information
alexryangregor authored Nov 17, 2023
1 parent c07e119 commit 45f53e7
Showing 1 changed file with 24 additions and 8 deletions.
32 changes: 24 additions & 8 deletions docs/source/workflows/bayts.rst
Original file line number Diff line number Diff line change
@@ -1,27 +1,43 @@
BayTS-based near real-time (NRT) forest disturbance alerts
==========================================================
BayTS NRT-FDM
=============

Background
----------

Near real-time forest disturbance monitoring (NRT-FDM) entails remote sensing techniques that are based on dense time series, targeting the generation of recent disturbance events in forested areas. The main objective is to detect those changes as early as possible. This prioritization usually comes at the cost of accuracy, as the methods are optimized for speed and timeliness. NRT-FDM is a rapidly evolving field of research and many different approaches have been proposed. The most common ones are based on the detection of abrupt changes in the time series, which are then classified as disturbance events.
Near real-time forest disturbance monitoring (NRT-FDM) entails remote sensing techniques that are based on dense time series, targeting the generation of recent disturbance events in forested areas.

This page provides background information and a detailed how-to guide for the implementation of the BayTS algorithm on the SEPAL platform that is optimized for the generation of NRT alerts using Sentinel-1 radar data. The advantage of using radar data is its independence from cloud cover and sensitivity to structural elements of the vegetation cover.
The main objective is to detect those changes as early as possible. This prioritization usually comes at the cost of accuracy, as the methods are optimized for speed and timeliness.

NRT-FDM is a rapidly evolving field of research and many different approaches have been proposed. The most common ones are based on the detection of abrupt changes in the time series, which are then classified as disturbance events.

This page provides background information and a detailed how-to guide for the implementation of the BayTS algorithm on the SEPAL platform that is optimized for the generation of NRT alerts using Sentinel-1 radar data. The advantage of using radar data is its independence from cloud cover and sensitivity to structural elements of vegetation cover.

Methodology
-----------
To detect forest cover loss in NRT in dense Sentinel-1 time series, a pixel-based approach is applied. First, a historic reference is created by calculating the mean and standard deviation of the backscatter values for each pixel. The reference is then used to calculate the probability of a pixel being in a disturbed state using the Bayes Theorem, a statistical approach that allows for the calculation of the probability of an event based on prior knowledge of conditions that might be related to the event. In this case, the event is a forest disturbance and the conditions are the backscatter values of the pixel.
To detect forest cover loss in NRT in dense Sentinel-1 time series, a pixel-based approach is applied.

First, a historic reference is created by calculating the mean and standard deviation of the backscatter values for each pixel.

Then, the reference is used to calculate the probability of a pixel being in a disturbed state using the Bayes Theorem, a statistical approach that allows for the calculation of the probability of an event based on prior knowledge of conditions that might be related to the event. In this case, the event is a forest disturbance and the conditions are the backscatter values of the pixel.

The method works without any training data, as a probability of being forest is derived from its historic reference expressed by the mean and standard deviation, from which a probability density function is derived. The non-forest state is assumed to have the same probability density function centered at a 4 decibel (dB) backscatter below the forested one. Once the probability of being non-forest exceeds a user-defined threshold (0.6 by default), Bayesian Updating is applied until the alert either gets confirmed or rejected. Therefore, the updating is repeated until the probability of being a change is above a high-confidence threshold (0.975 by default). If the alert is not confirmed within a certain time range (90 days by default), the alert is rejected.
The method works without any training data, as a probability of being forest is derived from its historic reference expressed by the mean and standard deviation, from which a probability density function is derived.

The non-forest state is assumed to have the same probability density function centered at a 4 decibel (dB) backscatter below the forested one. Once the probability of being non-forest exceeds a user-defined threshold (0.6 by default), Bayesian Updating is applied until the alert either gets confirmed or rejected. Therefore, updating is repeated until the probability of being a change is above a high-confidence threshold (0.975 by default). If the alert is not confirmed within a certain time range (90 days by default), the alert is rejected.

The method capitalizes on both VV and VH polarized channels, whereas the highest non-forest probability is considered for updating. Initial alerts from one orbit can only be confirmed by another orbit, as the image geometry changes, so that small-scale changes might be invisible in one orbit but visible in the other.

.. figure:: ../_images/workflows/bayts/bayts_pdfs.png
:alt: A single-pixel time series with apparent change, and the forest and non-forest probability to the right.
:alt: A single-pixel time series with apparent change, and the forest and non-forest probability to the right
:width: 800
:align: center

SEPAL allows for the displaying of different stages of confidence, which is the initial one (where an alert is triggered by the probability of being non-forest above 0.6), low-confidence alerts (where the probability of change being > 0.85) and high-confidence alerts (where the probability of change being > 0.975). In addition, the thresholds can be changed.
SEPAL allows for the displaying of different stages of confidence:

- the initial one (where an alert is triggered by the probability of being non-forest above 0.6);
- low-confidence alerts (where the probability of change being > 0.85); and
- high-confidence alerts (where the probability of change being > 0.975).

In addition, the thresholds can be changed.

The main sensitivity parameter affects the difference of the forest and non-forest probability functions. By default, it is 4 dB; a low-sensitivy multiplies the 4 dB by a factor of 1.2, while a higher sensitivity is achieved by decreasing the 4 dB range by a factor of 0.8, resulting in a 4.8 dB or 3.2 dB difference, respectively. This will affect both detectability of change and timeliness.

Expand Down

0 comments on commit 45f53e7

Please sign in to comment.