Next Article in Journal
Single Channel Circular SAR Moving Target Detection Based on Logarithm Background Subtraction Algorithm
Next Article in Special Issue
Evaluation of a Bayesian Algorithm to Detect Burned Areas in the Canary Islands’ Dry Woodlands and Forests Ecoregion Using MODIS Data
Previous Article in Journal
Assessment of Scattering Error Correction Techniques for AC-S Meter in a Tropical Eutrophic Reservoir
Previous Article in Special Issue
Burned Area Mapping of an Escaped Fire into Tropical Dry Forest in Western Madagascar Using Multi-Season Landsat OLI Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geostationary Sensor Based Forest Fire Detection and Monitoring: An Improved Version of the SFIDE Algorithm

by
Valeria Di Biase
1 and
Giovanni Laneve
2,*
1
Dipartimento di Ingegneria Astronautica, Elettrica e Energetica, Sapienza University of Rome, 00185 Roma, Italy
2
Scuola di Ingegneria Aerospaziale, Sapienza University of Rome, 00138 Roma, Italy
*
Author to whom correspondence should be addressed.
Remote Sens. 2018, 10(5), 741; https://doi.org/10.3390/rs10050741
Submission received: 18 April 2018 / Revised: 2 May 2018 / Accepted: 7 May 2018 / Published: 11 May 2018
(This article belongs to the Special Issue Remote Sensing of Wildfire)

Abstract

:
The paper aims to present the results obtained in the development of a system allowing for the detection and monitoring of forest fires and the continuous comparison of their intensity when several events occur simultaneously—a common occurrence in European Mediterranean countries during the summer season. The system, called SFIDE (Satellite FIre DEtection), exploits a geostationary satellite sensor (SEVIRI, Spinning Enhanced Visible and InfraRed Imager, on board of MSG, Meteosat Second Generation, satellite series). The algorithm was developed several years ago in the framework of a project (SIGRI) funded by the Italian Space Agency (ASI). This algorithm has been completely reviewed in order to enhance its efficiency by reducing false alarms rate preserving a high sensitivity. Due to the very low spatial resolution of SEVIRI images (4 × 4 km2 at Mediterranean latitude) the sensitivity of the algorithm should be very high to detect even small fires. The improvement of the algorithm has been obtained by: introducing the sun elevation angle in the computation of the preliminary thresholds to identify potential thermal anomalies (hot spots), introducing a contextual analysis in the detection of clouds and in the detection of night-time fires. The results of the algorithm have been validated in the Sardinia region by using ground true data provided by the regional Corpo Forestale e di Vigilanza Ambientale (CFVA). A significant reduction of the commission error (less than 10%) has been obtained with respect to the previous version of the algorithm and also with respect to fire-detection algorithms based on low earth orbit satellites.

Graphical Abstract

1. Introduction

Every year, annual reports highlight the extent of the forest fires phenomenon in the European Union, where more than half a million ha of forests burn in about 65,000 fires (Figure 1) [1,2]. This phenomenon is particularly relevant in the European Mediterranean region, where climatic factors produce a high level of vegetation stress during the summer months which causes a higher risk of inflammability.
In Italy, the amount of wildland and forest areas that burns each year is, on average, 50,000 ha, and this phenomenon is particularly relevant in regions such as Sardinia, where, from 1995 to 2009, an area of 16,600 ha per year burnt, and over 90% of these fires were human-induced [3].
For this reason, the adoption of a system capable of detecting as early as possible the trigger of new fires would considerably reduce environmental, material and social damage [4].
The use of satellite systems (NOAA/AVHRR, GOES, etc.), developed in the early 1970s by the US Forest Service [5], can be an extraordinary tool to detect hotspots and evaluate vegetation stress level.
A first classification of fire detection algorithms using remote data, distinguishes between those in which information is provided by a single channel (one-channel algorithms), and those in which data are provided from two or more channels (multi-channel algorithms). In both cases, it is essential to have information in the MWIR (Medium Wave InfraRed) region of the spectrum, and in particular in channels around 3.9 μm: the radiative emission peak for temperatures that characterize forest fires (800–1500 K) is in the range of 2–4 μm and, moreover, it is possible to work in atmospheric windows. In the case of multi-channel systems, it is common practice to choose a second channel in the TIR (Thermal InfraRed) region, at wavelengths between 8 μm and 12 μm, where the emission of bodies at temperatures close to those characterizing the Earth’s surface is maximum, as in Dozier [6].
Most of the current multi-channel algorithms use band at 3.9 μm in combination with band at 10.8 μm, available for the major satellite sensors used for fire detection. Remote data in these channels are often accompanied by the VIS (VISible) bands, used for the detection of pixels with cloud coverage and for the detection of false alarms.
A further distinction concerning fire detection systems is based on the characteristics of the orbit in which the satellite is located. For low orbit, the advantage of having images at a higher spatial resolution is accompanied by a low frequency of observation (approximately, two observations per day). Using instead weather satellites located in geostationary orbit, it is possible to have a continuous monitoring (images are usually acquired every 15’ or 30’), but with a significantly lower spatial resolution (from 3 to 5 km per pixel, depending on latitude). Despite the disadvantages of low spatial resolution as the minimum dimensions of detectable fires, this last solution is useful for the purpose of real-time monitoring.
MODIS (Moderate Resolution Imaging Spectroradiometer) and AVHRR (Advanced Very High Resolution Radiometer) are among the most commonly used low orbit satellite sensors for forest fire monitoring. In the last few years (2013), the new generation medium resolution imager, the VIIRS (Visible Infrared Imaging Radiometer Suite) satellite sensor, has been used to detect active fires, also.
MODIS radiometer allows the detection of data in 36 spectral bands from VIS to TIR with resolution of 250 m for channels 1–2, 500 m for channels 3–7, and 1 km for the remaining bands. MODIS sensors are located on NASA Aqua and Terra satellites. These two complementary satellites, allow the detection of images of the same region four times per day, recording two images during diurnal hours and two during nocturnal hours [7]. The works of Kaufman [8], Justice [9] and Giglio [10,11,12] are the most relevant for fire detection using MODIS data. AVHRR/3 is the third and last generation of the radiometer developed by NOAA (National Oceanic and Atmospheric Administration) in collaboration with NASA, used since 1978. It has 6 radiometric channels with 1 km resolution, and its orbit permits the acquisition of one or two images of the same area per day.
For fire detection using AVHRR data, the works of Robinson [13], Flasse and Ceccato [14] and Giglio [11,15] are of primary importance. The VIIRS sensor on the Suomi National Polar-orbiting Partnership (S-NPP) satellite incorporates fire-sensitive channels enabling active fire detection and characterization [16]. Concerning geostationary satellites: the works of Prins [17], Reid [18] and Xu [19] are most relevant for fire detection using GOES (Geostationary Operational Environmental Satellite) radiometer whereas the present work is focused on the images of the MSG satellite series, operated by EUMETSAT consortium. The SEVIRI (Spinning Enhanced Visible and Infrared Imager) sensor on board of MSG satellite can acquire images every 15 min, in twelve channels: three in Visible and Near Infra-Red (VNIR), and eight in Infra-Red (IR). These images have 3 km resolution at the sub-satellite point (the equator), with the exception of the High Resolution Visible (HRV) channel (12th), which provides 1 km resolution images [20].
For fire detection using MSG-SEVIRI data, works of Calle [21], Laneve [22], Roberts and Wooster [23] and Amraoui [24] are relevant. Calle proposes an algorithm capable of detecting fire with a minimum dimension of 0.7 ha on the Iberian land; Laneve shows a process for real-time coverage in Mediterranean area; Roberts and Wooster propose an algorithm in which false detection is less than 4% of observed fires; lastly, Amraoui proposes an algorithm for live coverage and combusted area rates on Africa.
Mediterranean areas are at the highest risk and suffer great losses of infrastructures, forested and agricultural land, human and livestock lives. The risk of such fires is expected to increase in forthcoming years under the impact of climate changes; vegetation becomes more inflammable (due to thermal stress and drought) and fire services are faced with difficulties when trying to suppress a fire due to increased inflammability and water shortage [4].
Previous studies demonstrate the capabilities of SEVIRI to detect fires, despite its low spatial resolution. In the Mediterranean areas thanks to its geostationary orbit, it can drastically reduce the reaction times compared to fire-detection based on low orbit systems. The main purpose of this paper is to show the utility of geostationary systems (MSG/SEVIRI) for fire-detection, focusing on fire occurrences on the Sardinian island, characterized by a non-homogeneous coverage (that is, the coexistence of urban settlements, infrastructure networks and vegetated areas (forest, agricultural and uncultivated areas) in a complex, dense and intricate patchwork) and high anthropization, issues neglected by current literature. In 2008, SIGRI (Sistema Integrato per la Gestione del Rischio Incendi) project was funded by ASI (Agenzia Spaziale Italiana) to face natural and human-induced disaster, and to demonstrate the abilities of supporting risk assessment, monitoring and management of fire by using remote-sensing satellite systems [25,26,27].
SEVIRI sensor is showing its ability to guarantee a live coverage by sending data every fifteen minutes. The CRPSM (Centro di Ricerca Progetto San Marco) has been studying for several years the possibility of using images acquired by SEVIRI: SFIDE (System for FIre DEtection) algorithm, presented in this paper, has been developed to exploit the images high update frequency by comparing temperature variation in subsequent images (change-detection), by using thresholds fit on the Sardinian island.
The paper is organized as follows: Section 2.1 describes exploited data and method for real time coverage, Section 2.2 describes SFIDE algorithm; Section 3 shows the results achieved by comparing hot spots detected by SFIDE, ground data provided by forest rangers CFVA (Corpo Forestale e di Vigilanza Ambientale) and hot spots detected by low orbit based fire-detection system.

2. Data and Methodology

The most common algorithms developed in the last decade, use remote data from low orbit satellite sensors (MODIS, AVHRR) or geostationary orbit satellite sensor (SEVIRI, GOES Imager) to apply absolute thresholds or contextual tests and obtain hot-spots as output results, together with the characterization of the fire pixels. In Section 2.1 a model of fire-detection algorithm and a description of input data is shown; in Section 2.2 how SFIDE algorithm uses absolute and contextual tests together with an innovative change-detection method to obtain a continuous monitoring of interested area from geostationary orbit data is illustrated.

2.1. Data & Methods

The European MSG satellite series are geostationary meteorological satellites operated by the EUMETSAT consortium. The most important sensor for the purpose of this paper is SEVIRI, comprising 11 spectral bands and a visible broadband (HRV) [Table 1].
Main applications of SEVIRI scanner are [28]: cloud detection through VIS 0.6 and VIS 0.8 channels; aerosol, soil humidity and vegetation index retrieval (IR 1.6, VIS 0.6, VIS 0.8); water vapor and wind determination (IR 6.2, IR 7.3); clouds and their temperature (IR 3.8, IR 8.7, IR 10.8, IR 12.0); high-atmosphere monitoring (IR 9.7); pressure at high altitude (IR 13.4).
SFIDE algorithm uses IR 3.8, IR 10.8, IR 12.0, VIS 0.6, VIS 0.8 for fire detection and estimation of surface temperature.
The common steps characterizing a fire-detection algorithm using remote sensing data are presented in Figure 2. The SEVIRI images are directly acquired at the University of Rome premises through the EUMETCAST service. After preliminary operations, consisting in the selection of interested area and the application of a sea-mask to exclude pixels in sea area, the brightness temperature (Tb) is evaluated whereas reflectance (rch) in VIS channels is retrieved knowing the values of radiance (R), Digital Number (DN), calibration coefficient (cf) and Offset (R0), data contained in the heading of SEVIRI images, archived in pgm (portable gray map) format.
One of the most common causes of false alarm in fire detection is the presence of clouds within the pixel, as clouds have similar features to fires: increase of Tb in MIR channel at 3.9 μm (Tb3.9) due to a major reflection of sunlight, decrease of Tb in TIR channel at 10.8 μm (Tb10.8), a consequent increase of brightness temperature difference (ΔT) between 3.9 μm and 10.8 μm channel. Various methods are adopted to detect cloud pixels. An algorithm similar to the one proposed by Saunders and Kriebel for AVHRR sensor [29] will be employed in this work. The algorithm applies fixed thresholds and distinguishes daily and nightly hours.
It is now possible to detect hot spots by applying specific thresholds. Using fixed thresholds for Tb3.9 and for ΔT, it is possible to distinguish two categories of pixels: non-fire pixel or potential-fire pixels. This distinction is generally made by using a multi-channel approach with fixed thresholds. A pixel is considered a potential hot-spot if it satisfies:
{ T b 3.9 > T b 3.9 M I N Δ T = T b 3.9 T b 10.8 > Δ T m i n
where Tb3.9MIN and ΔTmin are threshold values of Tb3.9 and ΔT, specifically selected considering vegetation coverage. Official reports by EUMETSAT [28] apply as fixed thresholds for a potential fire-pixel: 310 K and 5 K for daytime hours and 290 K and 0 K for nocturnal hours. In the last few years, a dependence on solar elevation, or solar zenith angle (SZA) has been shown. This condition allows the avoidance of omission errors in the hours in which SZA is higher.
Fixed thresholds expressed as linear function of SZA are:
{ T b 3.9 > T b 3.9 _ 0 c 3.9 × S Z A Δ T = T b 3.9 T b 10.8 > Δ T 0 c Δ T × S Z A
where Tb3.9_0 and ΔT0 are constants chosen for SZA = 0°; c3.9 and cΔT are coefficients evaluated as 0.3 [K/deg] and 0.0049 [K/deg] [23].
The algorithm usually continues by applying contextual thresholds, in some cases used for the validation of fire points [19,23], in other cases used as test for hot-spot detection [8,15]. Contextual algorithm is based on average values of T b 3.9 , Δ T and their standard deviations ( δ T b 3.9 e δ Δ T ) in n × n regions centered on the supposed fire point (from 3 × 3 to 5 × 5 pixels in case of GEO-Geostationary Earth Orbit-satellite based systems, to 21 × 21 pixels for LEO-Low Earth Orbit-satellite). Thresholds for the detection of fire points are:
{ T b 3.9 > T b 3.9 ¯ + k T b 3.9 δ T b 3.9 k 0 Δ T > Δ T ¯ + max ( k Δ T δ Δ T , Δ T M I N )
where T b 3.9 ¯ and Δ T ¯ are average values of T b 3.9 and Δ T in n × n region; k T b 3.9 and k Δ T are dimensionless coefficients with values ranging from 2 to 4; k 0 is a constant between 0 K and 3 K; Δ T M I N is a constant between 2.5 K and 6 K; max ( k Δ T δ Δ T , Δ T M I N ) is an operator that returns the maximum value between k Δ T δ Δ T and Δ T M I N . Giglio suggests for MODIS sensor: k T b 3.9 = 3 , k Δ T = 3.5 , k 0 = 0 K, Δ T M I N = 6 K [12]. Robert and Wooster suggest for SEVIRI sensor: k T b 3.9 = 2 , k Δ T = 2 , k 0 = 0 K, Δ T M I N = 2.5 K [23].
A further cause of false alarm is the sun-glint phenomenon, or the reflection of sunlight on water mirrors that causes an increase of T b 3.9 and Δ T . If there is a correspondence between the angle value ( ϑ g ) between the line connecting surface-satellite and the specular direction to the solar radiation affecting the water mirror (Figure 3), then high risk for coast line pixels occurs.
Even if water bodies are excluded from the analysis by applying a mask, the phenomenon could occur at water-land border pixels. This is particularly evident in GEO systems due to their low resolution, and it is given by ground and water co-presence in the same pixel, that lead to a contamination of the contextual algorithm.
More selective thresholds or further tests are therefore used for coast line pixels. Giglio [12] suggests for MODIS satellite:
{ ϑ g < 2 ° o r ϑ g < 10 °   and   r 0.6 > 0.1   a n d   r 0.8 > 0.2   a n d   r 2.1 > 0.12   o r ϑ g < 15 °   and   N w > 0  
where r 2.1 is the reflectance at 2.1 µm (channel 7); N w is the number of sea-pixels in 3 × 3 grid.
Hot-spots detected are now characterized in terms of dimensions (fire size in each pixel, ρ , and A f = ρ A P I X E L , where A P I X E L is the area of a pixel [m2]), fire temperature ( T f ), FRP (Fire Radiative Power) and confidence level of the fire point.
Dimension and temperature of fire can be obtained by the two-channel approach (at 3.9 µm and 10.8 µm) proposed by Dozier [6], based on the estimation of the radiant components of the pixel in absence of fire. Since the Dozier algorithm is based on a system of equations in two variables the active surface area can be estimated if the emissivity of the fire is known, otherwise we can determine their product.
FRP and its integration over time, FRE (Fire Radiative Energy), are further characterizations to estimate burned biomass of the pixel. In a first, two-channel approach, after solving equations proposed in [6], T f and ρ are obtained and it is possible to evaluate the Fire Radiative Power using Stefan–Boltzmann law [30]:
F R P S B = A f ε f σ ( T f 4 T b 4 )
where F R P S B is fire radiative power estimated by Stefan–Boltzmann in [W]; ε f is the emissivity of A f ; σ is Stefan–Boltzmann constant.
Another approach is offered by Wooster equation [31], and it does not require values of T f  and ρ :
F R P W = A P I X E L σ ε f a ε f , 3.9 L 3.9
where F R P W is fire radiative power estimated by Wooster in [W]; εf is the emissivity of the fire and εf,3.9 is spectral emissivity at about 3.9 μm, a is a typical constant dependent from used sensor (a = 3.06 × 10−9 for SEVIRI). In Wooster et al. [31] gray body behavior is assumed, that is (ɛf = ɛf,3.9).
After estimating FRP using described methods, it is possible to evaluate FRE by integrating FRP over time:
F R E = t f i r e F R P d t
where t f i r e is the fire duration.
Using the following equation, the estimation of burned biomass (BB) is obtained [23]:
B B C r F R E
where Cr is a conversion factor that links the FRP to combustion rate. In Wooster et al. [31] a universal value of C r = 0.368 [kg/MJ] is proposed, while Kaiser et al. [32] proposed specific values for different land cover classes. Combusted biomass density ( ρ B ) is [31]:
ρ B B B A f
Finally, it is possible to estimate a quality index of hot-spot which takes into account the amount of clouds in the area of interest and their position with respect to the hot spot. Current literature provides further heuristic methods able to estimate confidence level (C) of each hotspot, as the one proposed by Giglio for MODIS [10,12] or the ones proposed by Roberts and Wooster [23] and Xu et al. [19] for GOES.
The data used to validate the results of the algorithm herein described include:
  • Fire events reported by CFVA (Corpo Forestale e di Vigilanza Ambientale) of the Sardinia Region, that is the agency involved in the fire-fighting;
  • Hot spots detected by the MODIS and VIIRS satellite sensors.

2.2. The SFIDE Algorithm

In this section SFIDE algorithm and its innovative aspects are described. It is important to highlight that the algorithm has been developed to be efficient in Mediterranean areas and, in particular, validation tests have been carried out in Sardinia. The complexity of the area of interest, because of its non-homogeneous and anthropized coverage, led to the elaboration of different conditions compared to the ones used in current literature.
Because of its objective to perform real time monitoring using geostationary satellite, one of the innovation of SFIDE algorithm consists in fire detection based on change-detection thresholds, evaluated at 15 min and 30 min intervals. Anomalous variations of T b 3.9 and Δ T in these intervals are considered as a sign of new fire event.
Further innovation is the introduction of thresholds which aim at detecting fire pixels using statistical analysis of T b 3.9 and Δ T as function of SZA: false alarms due to fixed thresholds have been overcome, as well as the ones due to linearly dependent thresholds by SZA.
The SFIDE algorithm prospect for diurnal data analysis (SZA < 85°) is shown in Figure 4. Principal steps of the algorithm are the following:
  • Preliminary operations and applications of masks;
  • Fixed thresholds;
  • Potential hot-spot detection;
  • Hot-spot confirmation;
  • Hot-spot characterization.

2.2.1. Preliminary Operations

The first significant stage consists of applying typical preliminary operation: the selection of the area of interest, application of the sea-mask (Figure 5) and transformation of the DN (Digital Number) in radiances and then in reflectances (channels 1, 2) and brightness temperatures (channels 4, 9, 10, 11).
The next step is the application of a cloud-mask including new thresholds appropriately modified to be applied to Sardinia, as the introduction of thresholds on channel 2 and its relative reflectance ( r 0.8 ). Similarly to [29], a pixel is considered as cloud covered if it satisfies:
{ r 0.6 +   r 0.8 > 1 or T b 12.0 < 265   K   or r 0.6 + r 0.8 > 0.7   a n d   T b 12.0 < 285   K  
Other pixels, even not masked as cloudy could be excluded by the algorithm because of their high reflectance.
Reflectance in channel 2 is also used to exclude partly covered pixels because of their high reflectance, which could cause false alarms, and also pixels that satisfy:
r 0.8 > 0.35
are excluded.
Value in Equation (11) is more flexible than the one proposed in [15] to reduce omission errors due to the exclusion of elevated reflectance pixels.
The importance of excluding cloudy pixels will be shown again in the step of hotspot confirmation.

2.2.2. Fixed Thresholds

Omission errors can occur when both (10) and (11) are implemented; to avoid this, a fixed threshold has been imposed on :
T b 3.9 > 318 K
Pixels that satisfy (12) are classified as hot-spots. Temperature of 318 K has been chosen after analysis of T b 3.9 trend in the considered area. An investigation of seasonal trend of T b 3.9 in the study area showed that all pixels that satisfies (12), together with adequate values of FRP (see next), has been validated as fire-pixels. As shown in Figure 4, this represents the first test used to detect hot-spots.

2.2.3. Potential Hot-Spot Detection

Two kinds of thresholds are used to detect hot-spots: fixed thresholds and SZA dependent thresholds. Thresholds expressed as linear function of SZA are generally preferred to the fixed ones, but omission and commission errors are not avoided. To model with more accuracy the average values of T b 3.9 and Δ T , an analysis of their average trend as function of SZA on Sardinia has been conducted. Their trend and their standard deviation have been stored and described by a third-degree function, used as potential hot-spot detection threshold. Obtained functions are:
T b 3.9 = 6.2410 6 · S Z A 3 0.0027 · S Z A 2 ± 0.052 · S Z A + 305.43
Δ T = 4.7510 6 · S Z A 3 0.0011 · S Z A 2 ± 0.018 · S Z A + 3.69
where T b 3.9 and Δ T are respectively expected values of T b 3.9 and Δ T as a function of SZA. In the terms with odd exponent, the use of superior sign is for post-meridian hours, and inferior sign is for anti-meridian hours. A pixel is considered a hot-spot if it satisfies:
{ T b 3.9 > T b 3.9 Δ T = T b 3.9 T b 10.8 > Δ T

2.2.4. Hot-Spot Confirmation

At this stage, checks on potential hot-spots are carried out to confirm them. Checks follow two paths: change-detection and contextual analysis. For both methodologies, narrower thresholds have been used for high false alarm risk pixel due to its partial cloud coverage, sudden change in reflection over time, or Sun-glint.
Change-detection: This innovative approach is particularly efficient in detecting new fires earlier than contextual analysis. Statistical analysis has been conducted on average trend of Tb3.9 and ΔT as function by SZA at 15 and 30 min intervals. Data on their average variations have been interpolated in third degree function dependent on SZA and used as pre-thresholds. Lastly, these functions have been used to evaluate a kind of thresholds, named trigger thresholds, that allow one to confirm fires in growing phase because of the joint increase of Tb3.9 and ΔT above statistically evaluate values.
The following symbolism has been used:
  • Tb3.9_15’ and ΔT15’: Tb3.9 and ΔT values for image that precedes by 15 min the current one;
  • Tb3.9_30’ and ΔT30’: Tb3.9 and ΔT values for image that precedes by 30 min the current one;
  • ΔTb3.9_15’ = (Tb3.9 − Tb3.9_15’) and Δ(ΔT)15’ = (ΔT − ΔT15’): Tb3.9 and ΔT variations between current image and the one that precedes at 15 min;
  • ΔTb3.9_30’ = (Tb3.9 − Tb3.9_30’) and Δ(ΔT)30’ = (ΔT − ΔT30’): Tb3.9 and ΔT variations between current image and the one that precedes at 30 min.
For the thermal variations of 15 min, evaluated as ΔTb3.9_15’ and Δ(ΔT)15’ trend and their respective standard deviations are:
Δ T b 3.9 _ 15 = 2.91 E 7 SZA 3 1.75 E 5 SZA 2 ± 4.39 E 4 SZA + 0.49
σ Δ T b 3.9 _ 15 = ± 1.00 E 6 SZA 3 5.09 E 5 SZA 2 1.77 E 2 SZA + 0.21
Δ ( Δ T ) 15 = ± 5.03 E 7 SZA 3 1.21 E 6 SZA 2 6.84 E 3 SZA + 0.005
σ Δ ( Δ T ) 15 = 7.17 E 7 SZA 3 8.81 E 5 SZA 2 ± 1.75 E 3 SZA + 0.85
where Δ T b 3.9 _ 15 and Δ ( Δ T ) 15 are functions that approximate the average trend of ΔTb3.9_15’ and Δ(ΔT)15’, while σ Δ T b 3.9 _ 15 and σ Δ ( Δ T ) 15 are functions that approximate their respective average standard deviations. Similarly, an analysis considering 30 min intervals has been conducted:
Δ T b 3.9 _ 30 = ± 1.95 E 6 SZA 3 1.25 E 4 SZA 2 3.46 E 2 SZA + 0.48
σ Δ T b 3.9 _ 30 = 4.39 E 7 SZA 3 6.07 E 6 SZA 2 ± 1.21 E 3 SZA + 0.75
Δ ( Δ T ) 30 = ± 9.13 E 7 SZA 3 6.40 E 6 SZA 2 1.34 E 2 SZA + 0.026
σ Δ ( Δ T ) 30 = 1.18 E 6 SZA 3 1.09 E 4 SZA 2 ± 3.56 E 3 SZA + 1.16
where Δ T b 3.9 _ 30 and Δ ( Δ T ) 30 are functions that approximate the average trend of ΔTb3.9_30’ and Δ(ΔT)30’, while σ Δ T b 3.9 _ 30 and σ Δ ( Δ T ) 30 are functions that approximate their respective average standard deviations.
Functions described in Equations (16)–(23) are trigger thresholds imposed on potential hot-spots that satisfy Equation (15) to confirm fire in the pixel. Considering thermal variations in 15’, a potential hot-spot is confirmed if:
{ Δ T b 3.9 _ 15 = ( T b 3.9 T b 3.9 _ 15 ) > Δ T b 3.9 _ 15 + σ Δ T b 3.9 _ 15 Δ ( Δ T ) 15 = ( Δ T Δ T 15 ) > Δ ( Δ T ) 15 + σ Δ ( Δ T ) 15 + f 15 ( Δ r 0.6 ) T b 3.9 > T b 3.9 ¯ + 1.5 Δ T > Δ T ¯ + 0.5 N w + N c = 0
where T b 3.9 ¯ and Δ T ¯ are average values of Tb3.9 and ΔT in 3 × 3 grid centered on potential hot-spots; Nw and Nc are respectively number of water-pixels and cloudy-pixels in 3 × 3 grid centered on potential hot-spots; f 15 ( Δ r 0.6 ) = 100 ( r 0.6 r 0.6 _ 15 ) if r 0.6 > r 0.6 _ 15 , with r 0.6 _ 15 reflectance in channel 1 evaluated 15 min before current image. Function f15’ is necessary to correct ΔT increase due to r 0.6 increase in 15’ intervals.
Considering thermal variations in 30’, a potential hot-spot is confirmed if it satisfies:
{ Δ T b 3.9 _ 30 = ( T b 3.9 T b 3.9 _ 30 ) > Δ T b 3.9 _ 30 + σ Δ T b 3.9 _ 30 Δ ( Δ T ) 30 = ( Δ T Δ T 30 ) > Δ ( Δ T ) 30 + σ Δ ( Δ T ) 30 + f 30 ( Δ r 0.6 ) T b 3.9 > T b 3.9 ¯ + 1.5 Δ T > Δ T ¯ + 0.5 N w + N c = 0
where f 30 ( Δ r 0.6 ) = 100 ( r 0.6 r 0.6 _ 30 ) if r 0.6 > r 0.6 _ 30 , with r 0.6 _ 30 reflectance in channel 1 evaluated 30 min before current image.
Similar thresholds are applied to elevated risk of false alarm pixels:
| r 0.6 r 0.6 _ 15 | 0.03
| r 0.6 r 0.6 _ 30 | 0.03
( N w + N c ) > 0
( r 0.8 r 0.6 ) 0.1
Equations (26) and (27) are imposed to detect pixels with sudden variation of r 0.6 , and consequently ΔT increase, due to cloudy sky. Equation (27) allows the prevention of drawbacks due to coastal line (Nw > 0) or near cloudy pixels (Nc > 0) points.
If one of the conditions (26)–(29) is satisfied then conditions of Equation (24) are replaced with:
{ Δ T b 3.9 _ 15 = ( T b 3.9 T b 3.9 _ 15 ) > Δ T b 3.9 _ 15 + 2 σ Δ T b 3.9 _ 15 Δ ( Δ T ) 15 = ( Δ T Δ T 15 ) > Δ ( Δ T ) 15 + 2 σ Δ ( Δ T ) 15 + f 15 ( Δ r 0.6 ) T b 3.9 > T b 3.9 ¯ + 1.5 Δ T > Δ T ¯ + 0.5 N w + N c = 0
while limits of Equation (25) are replaced with:
{ Δ T b 3.9 _ 30 = ( T b 3.9 T b 3.9 _ 30 ) > Δ T b 3.9 _ 30 + 2 σ Δ T b 3.9 _ 30 Δ ( Δ T ) 30 = ( Δ T Δ T 30 ) > Δ ( Δ T ) 30 + 2 σ Δ ( Δ T ) 30 + f 30 ( Δ r 0.6 ) T b 3.9 > T b 3.9 ¯ + 1.5 Δ T > Δ T ¯ + 0.5 N w + N c = 0

2.2.5. Contextual Analysis

An alternative method for hot-spot confirmation is a contextual analysis that makes use of information extracted from the surrounding pixels.
Two confirmation conditions are imposed: low-probability and high-probability of fire point conditions. For low-probability condition point, the hot-spot is confirmed if:
{ T b 3.9 > T ¯ b 3.9   + max ( 1 , σ T b 3.9 3 ) Δ T > Δ T ¯ + max ( 1.25 , σ Δ T ) o r   Δ T > Δ T ¯ + max ( 2 , σ Δ T )   or   Δ T > 4.5
where σ T b 3.9 and σ Δ T are standard deviations of Tb3.9 and ΔT in the immediate neighborhood. If at least one condition between (26), (27), (29) or the following ones is satisfied:
r 0.6 > 0.15
r 0.6 > r 0.6 ¯ + σ r 0.6
r 0.6 ¯ < 0.1
min ( r 0.6 ) 3 x 3 < 0.08
the pixel is considered a high-probability fire. r 0.6 ¯ , σ r 0.6 and min ( r 0.6 ) 3 × 3  are respectively: average value, standard deviation and minimum value of  in 3 × 3 grid.
For high-probability points, a hot-spot is confirmed if:
{ T b 3.9 > T b 3.9 ¯ + max ( 2.5 , σ T b 3.9 3 ) Δ T > Δ T ¯ + min ( 4 , 2 σ Δ T )  

2.2.6. Hot-Spot Characterization: FRP and LST

For FRP evaluation, Wooster approximation (6) has been used [31]. FRPW allows one to obtain Tf and ρ values. FRP values are useful to evaluate fire dimension and power, and therefore to quantify resources needed for each fire. Furthermore, imposing a limit value for FRP:
F R P > 40   [ M W ]
false alarms can be avoided.
Land Surface Temperature (LST) has been introduced in fire characterization. In future developments, LST value could be used instead of Tb for fire-detection thresholds and FRP evaluations. For LST evaluation, double split window technique (DSWT) has been used, with Tb10.8 and Tb12 as input data [33]. Used algorithm is VZA (View Zenith Angle) dependent.

2.2.7. Nightly Hours Algorithm

For SZA < 85°, an alternative algorithm has been used due to the absence of reflected solar component. Principal differences with the daily one are:
  • Cloud-mask identified using a threshold on channel 10;
  • Hot-spot detection with fixed thresholds occurs with limit values of Tb3.9 and ΔT, adapted to nocturnal hours;
  • Fixed values of Tb3.9 and ΔT as hot-spot detection thresholds;
  • Contextual analysis carried out by using average value of Tb3.9 and ΔT on the whole of Sardinia;
  • No distinction between high and low probability fire pixels.
A model of the night-time algorithm is given in Figure 6.
A pixel is considered cloudy if:
T b 12.0 < 265 K
Fixed thresholds based on Tb3.9 and ΔT values [34] are:
{ T b 3.9 > 290   K Δ T > 1   K
Due to the high number of omission errors detected using fixed thresholds, also for nocturnal hours a detection method based on fixed pre-thresholds and confirmation has been used.
A pixel is considered as potential hot-spot if it satisfies the following revised thresholds:
{ T b 3.9 > 285   K Δ T > 2   K
Hot-spot confirmation is made with contextual-analysis. A pixel is confirmed as hot-spot if it satisfies:
{ T b 3.9 > T b 3.9 _ S A R + 1.5 σ T b 3.9 _ S A R Δ T > Δ T _ S A R + 1.5 σ Δ T _ S A R
where T b 3.9 _ S A R and σ T b 3.9 _ S A R are average values and standard deviation of Tb3.9 for the whole of Sardinia; Δ T _ S A R and σ Δ T _ S A R are average values of ΔT for the whole of Sardinia.
A final confirmation is given by FRP analysis. A pixel is confirmed as hot-spot only if (38) is satisfied.

3. Data Validation

Validation of SFIDE output has been made on Sardinian area using images acquired in 2014 (from 2 to 8 July and from 31 July to 6 August) and August 2016. These results have been compared with in situ data collected by CFVA (http://www.sardegnaambiente.it/corpoforestale/) and this comparison concerns fires detectable from geostationary orbit because of their dimension and coverage features. Another validation has been made by using the low orbit data of MODIS sensor (https://modis.gsfc.nasa.gov). It should be underlined that the fire season in the area of interest commonly goes from the 1 June to the 30 September.

3.1. Comparison with In Situ Measurements

Since CFVA measurements consist in total burned area evaluated at the end of the incendiary event, only fires larger- in total- than 5 ha have been considered. Instantaneous in situ measurements are not available, and minimum detectable fire size from geostationary orbit has been evaluated using [6] in 0.1–0.2 ha [35,36,37]. Fires in cloudiness and high reflectance area have been excluded from the comparison. The results (Table 2) show that 41 on 45 detectable fires were detected (91.1%), with less than 9% omission error.
Furthermore, a significant result is that not only fires that developed in forest area were detected: a certain percentage of detected hot-spots corresponds to fires that occurred in mixed or grassy covered areas. Even a considerable number of fire having burned less than 5 hectares have been detected, with an amount of 45 fires detected during the validation period.
SFIDE output consists in 464 hot-spots, of which only 32 (6.9%) are commission errors (false alarms). Trigger thresholds (Equations (24) and (25) or (30) and (31)) have been able to detect 31.25% of detectable fires. The utility of the change detection based approach is confirmed by the fact that the tests of Equations (23) and (24) (or Equations (30) and (31)) allow an early detection of the fires in 100% of the cases. It is necessary to highlight that these thresholds are decisive for the triggering phases and, moreover, trigger thresholds are able to detect fires only when they are in their increasing phase, that is essential for the purpose of real time monitoring of SFIDE.
Other hot-spots has also been detected from contextual thresholds (about 73.49%) and fixed thresholds (about 62.07%).
One of the principal causes of false alarms is the sudden variation of reflectance values in channel 2, together with a sudden variation of Tb3.9 and consequently ΔT. In nocturnal hours, when these phenomena do not appear, 100% of detected hot-spots has indeed been classified as real fire. An in-depth analysis of the hot-spots classified as false alarms, shows they mainly depends by fixed threshold.
By considering the purpose of real time monitoring, the optimal solution would be to use the trigger threshold for the detection of new fires, and the fixed threshold in the following phases to continue to consider pixel as burning and continue to estimate the FRP.
The tests carried out in 2016 shows how the purpose of SFIDE is achieved also using samples of data different from the ones used for the statistical analysis (as in 2014), and how the introduction of new thresholds is necessary for the purpose of real time monitoring.
Figure 7, processed with commercial software ArcMap, shows a comparison between SFIDE detected hot-spots (red square) and in situ measurements of CFVA (blue area) on 2, 3 and 8 of July and 5 of August 2014. Figure corresponding to 3rd July (Figure 7b) shows non-detected fires due to a local high cloudiness. Red squares near in situ reported fires represent a kind of detected hot-spot due to sensor saturation issues.

3.2. Comparison with MODIS Sensor

Comparison with MODIS sensor has been conducted with images acquired in 2014 (from 2 to 8 July and from 31 July to 6 August), MODIS algorithm results over evaluated period, consist in 3 fires compared to 20 fires detected by SFIDE algorithm and 22 reported by CFVA. Only one of the fires detected by MODIS was not-detected by SFIDE due to a meteorological condition (local high cloudiness). These results show:
the importance of a geostationary orbit sensor that allows real time monitoring due to a high update frequency: 96 daily images by SEVIRI and 2–4 by MODIS;
that, particularly relevant, fires larger than 100 [ha] have not been detected by MODIS sensor, and this result confirms the major advantages of SFIDE algorithm also regarding the number of detected fires.

3.3. Fire Monitoring Examples

Using the high refresh frequency of the SEVIRI images, in particular in the case of MSG 9 also known as RSS (Rapid Scanning Service) characterized by a refresh frequency of 5 min, it is, in general, possible to follow in a quasi-continuous way the evolution of the fire. This allows, possibly, the deduction of interesting information on the fire behavior and the effectiveness of the fire-fighting activity. Figure 8 allows one to explain this point better. The figure compares, through the FRP (measured in MW), the behavior of two fires that occurred in the Calabria region on 20 August 2017. The behavior of one of the two fires (the one for which the FRP is shown on the left plot) seems much more ‘natural’ or regular than the other. Further, this fire lasts also much less than the other, about 1.5 h (90 min) instead of 11 h (660 min). It is also less intense reaching a maximum FRP of 180 MW instead of 330 MW, as in the case of the other fire. This suggests that fire corresponding to left plot of Figure 8 extinguished itself for natural reasons not for the intervention of firefighters.

4. Conclusions

The SFIDE algorithm was developed by authors several years ago [22] in the framework of a project funded by the Italian Space Agency devoted to implement a forest fires early detection and monitoring system based on satellite images coming from the geostationary sensor SEVIRI on board of the MSG satellites series. This paper presents a new version of the algorithm aiming, in particular, at reducing the false alarm rate. This updated version of the algorithm takes into account other previously developed methods for fire detection in order to:
Improve the estimate of the reference temperature used to define a pixel as interested by a fire;
Improve the cloud mask accuracy thus reducing the risk of false alarms;
Exploit better the high refresh rate of the images to implement several tests for accurate detection of forest fires.
The results of the comparison of the events detected in Sardinia by the algorithm with those provided by local authorities and/or detected by other higher resolution satellite sensor (MODIS) confirm the improvement made. In fact, in the present version the algorithm provides an omission error of less than 9% and a commission error of less than 7%. At this point a further step is needed, which is, to reduce the approximation in the localization of the fire (the fire can occupy any position in the pixel of 4 × 4 km2). At present, the position of the fire corresponds to the coordinates of the center of the pixel. Possibly exploiting the simultaneous observation from 2 geostationary satellites like MSG 8 (RSS) and MSG 9 (IODC) and introducing a fuel map of the area of interest a more accurate localization of the fire position could be achieved. This will possibly be discussed in a future paper.

Author Contributions

G.L. conceived the SFIDE algorithm, supervised the research, provided suggestions and support in writing the paper and prepared the data sets. V.D.B. implemented the improvement to the SFIDE program, conducted the experiments, analyzed the results and wrote the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Schmuck, G.; San-Miguel-Ayanz, J.; Camia, A.; Durrant, T.; Santos De Oliveira, S.; Boca, R.; Whitmore, C.; Giovando, C.; Libertà, G.; Corti, P.; et al. Forest Fires in Europe; Publications Office of the European Union: Luxembourg, 2010. [Google Scholar]
  2. Laneve, G.; Fusilli, L.; Marzialetti, P.; De Bonis, R.; Bernini, G.; Tampellini, L. Development and Validation of Fire Damage-Severity Indices in the Framework of the PREFER Project. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 2806–2817. [Google Scholar] [CrossRef]
  3. Ager, A.A.; Preisler, H.K.; Arca, B.; Spano, D.; Salis, M. Wildfire risk estimation in the Mediterranean area. Environmetrics 2009, 25, 384–396. [Google Scholar] [CrossRef]
  4. Sifakis, N.I.; Iossifidis, C.; Kontoes, C.; Keramitsoglou, I. Wildfire Detection and Tracking over Greece Using MSG-SEVIRI Satellite Data. Remote Sens. 2011, 3, 524–538. [Google Scholar] [CrossRef]
  5. Martín, M.P.; Ceccato, P.; Flasse, S.; Downey, I. Fire detection and fire growth monitoring using satellite data. In Remote Sensing of Large Wildfires; Springer: Berlin/Heidelberg, Germany, 1999; pp. 101–122. [Google Scholar]
  6. Dozier, J. A method for satellite identification of surface temperature fields of subpixel resolution. Remote Sens. Environ. 1981, 11, 221–229. [Google Scholar] [CrossRef]
  7. Salomonson, V.V.; Barnes, W.L.; Maymon, P.W.; Montgomery, H.E.; Ostrow, H. MODIS: Advanced facility instrument for studies of the earth as a system. IEEE Trans. Geosci. Remote Sens. 1989, 27, 145–153. [Google Scholar] [CrossRef]
  8. Kaufman, Y.J.; Justice, C.O.; Flynn, L.P.; Kendall, J.D.; Prins, E.M.; Giglio, L.; Ward, D.E.; Menzel, W.P.; Setzer, A.W. Potential global fire monitoring from EOS-MODIS. J. Geophys. Res. Atmos. 1998, 103, 32215–32238. [Google Scholar] [CrossRef]
  9. Justice, C.O.; Giglio, L.; Korontzi, S.; Owens, J.; Morisette, J.T.; Roy, D.; Descloitres, J.; Alleaume, S.; Petitcolin, F.; Kaufman, Y. The MODIS fire products. Remote Sens. Environ. 2002, 83, 244–262. [Google Scholar] [CrossRef]
  10. Giglio, L.; Descloitres, J.; Justice, C.O.; Kaufman, Y.J. An enhanced contextual fire detection algorithm for MODIS. Remote Sens. Environ. 2003, 87, 273–282. [Google Scholar] [CrossRef]
  11. Giglio, L.; Csiszar, I.; Restas, A.; Morisette, J.T.; Schroeder, W.; Morton, D.; Justice, C.O. Active fire detection and characterization with the advanced spaceborne thermal emission and reflection radiometer (ASTER). Remote Sens. Environ. 2008, 112, 3055–3063. [Google Scholar] [CrossRef]
  12. Giglio, L.; Schroeder, W.; Justice, C.O. The collection 6 MODIS active fire detection algorithm and fire products. Remote Sens. Environ. 2016, 178, 31–41. [Google Scholar] [CrossRef]
  13. Robinson, J.M. Fire from space: Global fire evaluation using infrared remote sensing. Int. J. Remote Sens. 1991, 12, 3–24. [Google Scholar] [CrossRef]
  14. Flasse, S.P.; Ceccato, P. A contextual algorithm for AVHRR fire detection. Int. J. Remote Sens. 1996, 17, 419–424. [Google Scholar] [CrossRef]
  15. Giglio, L.; Kendall, J.D.; Justice, C.O. Evaluation of global fire detection algorithms using simulated AVHRR infrared data. Int. J. Remote Sens. 1999, 20, 1947–1985. [Google Scholar] [CrossRef]
  16. Schroeder, W.; Oliva, P.; Giglio, L.; Csiszar, I.A. The new VIIRS 375 m active fire detection data product: Algorithm description and initial assessment. Remote Sens. Environ. 2014, 143, 85–96. [Google Scholar] [CrossRef]
  17. Prins, E.M.; Feltz, J.M.; Menzel, W.P.; Ward, D.E. An overview of GOES-8 diurnal fire and smoke results for SCAR-B and 1995 fire season in South America. J. Geophys. Res. Atmos. 1998, 103, 31821–31835. [Google Scholar] [CrossRef]
  18. Reid, J.S.; Koppmann, R.; Eck, T.F.; Eleuterio, D.P. A review of biomass burning emissions part II: Intensive physical properties of biomass burning particles. Atmos. Chem. Phys. 2005, 5, 799–825. [Google Scholar] [CrossRef]
  19. Xu, W.; Wooster, M.J.; Roberts, G.; Freeborn, P. New GOES imager algorithms for cloud and active fire detection and fire radiative power assessment across North South and Central America. Remote Sens. Environ. 2010, 114, 1876–1895. [Google Scholar] [CrossRef]
  20. Aminou, D.M.A. MSG’s SEVIRI instrument. ESA Bull. 2002, 111, 15–17. [Google Scholar]
  21. Calle, A.; Casanova, J.L.; Romo, A. Fire detection and monitoring using MSG Spinning Enhanced Visible and Infrared Imager (SEVIRI) data. J. Geophys. Res. Biogeosci. 2006, 111, G04S06. [Google Scholar] [CrossRef]
  22. Laneve, G.; Castronuovo, M.M.; Cadau, E.G. Continuous monitoring of forest fires in the Mediterranean area using MSG. IEEE Trans. Geosci. Remote Sens. 2006, 44, 2761–2768. [Google Scholar] [CrossRef]
  23. Roberts, G.J.; Wooster, M.J. Fire detection and fire characterization over Africa using Meteosat SEVIRI. IEEE Trans. Geosci. Remote Sens. 2008, 46, 1200–1218. [Google Scholar] [CrossRef]
  24. Amraoui, M.; DaCamara, C.C.; Pereira, J.M.C. Detection and monitoring of African vegetation fires using MSG-SEVIRI imagery. Remote Sens. Environ. 2010, 114, 1038–1052. [Google Scholar] [CrossRef]
  25. Laneve, G.; Jahjah, M.; Ferrucci, F.; Battazza, F. SIGRI project: The development of the fire vulnerability index. In Proceedings of the 34th International Symposium on Remote Sensing of Environment, Sydney, Australia, 10–15 April 2011. [Google Scholar]
  26. Laneve, G.; Cadau, E.; Ferrucci, F.; Rongo, R.; Guarino, A.; Fortunato, G.; Loizzo, R. SIGRI-an Integrated System for Detecting, Monitoring, Characterizing Forest Fires and Assessing damage by LEO-GEO Data. Ital. J. Remote Sens. 2012, 44, 19–25. [Google Scholar] [CrossRef]
  27. Laneve, G.; Jahjah, M.; Ferrucci, F.; Hirn, B.; Battazza, F.; Fusilli, L.; De Bonis, R. SIGRI project: Products validation results. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 895–905. [Google Scholar] [CrossRef]
  28. Schmid, J. The SEVIRI instrument. In Proceedings of the 2000 EUMETSAT Meteorological Satellite Data User’s Conference, Bologna, Italy, 29 May–2 June 2000; Volume 29. [Google Scholar]
  29. Saunders, R.W.; Kriebel, K.T. An improved method for detecting clear sky and cloudy radiances from AVHRR data. Int. J. Remote Sens. 1988, 9, 123–150. [Google Scholar] [CrossRef]
  30. Zhukov, B.; Lorenz, E.; Oertel, D.; Wooster, M.; Roberts, G. Spaceborne detection and characterization of fires during the bi-spectral infrared detection (BIRD) experimental small satellite mission (2001–2004). Remote Sens. Environ. 2006, 100, 29–51. [Google Scholar] [CrossRef]
  31. Wooster, M.J.; Roberts, G.; Perry, G.L.W.; Kaufman, Y.J. Retrieval of biomass combustion rates and totals from fire radiative power observations: FRP derivation and calibration relationships between biomass consumption and fire radiative energy release. J. Geophys. Res. Atmos. 2005, 110. [Google Scholar] [CrossRef]
  32. Kaiser, J.W.; Heil, A.; Andreae, M.O.; Benedetti, A.; Chubarova, N.; Jones, L.; Morcrette, J.J.; Razinger, M.; Schultz, M.G.; Suttie, M.; et al. Biomass burning emissions estimated with a global fire assimilation system based on observed fire radiative power. Biogeosciences 2012, 9, 527–554. [Google Scholar] [CrossRef] [Green Version]
  33. Sobrino, J.A.; Romaguera, M. Land surface temperature retrieval from MSG1-SEVIRI data. Remote Sens. Environ. 2004, 92, 247–254. [Google Scholar] [CrossRef]
  34. Joro, S.; Samain, O.; Yildirim, A.; Van De Berg, L.; Lutz, H.J. Towards an Improved Active Fire Monitoring Product for MSG Satellites. EUMETSAT. 2008. Available online: http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.611.1456&rep=rep1&type=pdf (accessed on 9 May 2018).
  35. Giglio, L.; Schroeder, W. A global feasibility assessment of the bi-spectral fire temperature and area retrieval using MODIS data. Remote Sens. Environ. 2014, 152, 166–173. [Google Scholar] [CrossRef]
  36. Wooster, M.J.; Roberts, G.; Perry, G.L.W.; Kaufman, Y.J. Retrieval of biomass combustion rates and totals from fire radiative power observations: Application to southern Africa using geostationary SEVIRI imagery. J. Geophys. Res. 2005, 110, D21111. [Google Scholar]
  37. Laneve, G.; Castronuovo, M.; Cadau, E. Assessment of the fire detection limit using SEVIRI/MSG sensor. In Proceedings of the IGARSS, Denver, CO, USA, 31 July–4 August 2006; pp. 4157–4160. [Google Scholar]
Figure 1. Averaged annual burned area in the European countries most affected by fires between 1980 and 2013, compared with more recent years [2].
Figure 1. Averaged annual burned area in the European countries most affected by fires between 1980 and 2013, compared with more recent years [2].
Remotesensing 10 00741 g001
Figure 2. The common steps characterizing a fire-detection algorithm using remote sensing data (IR = InfraRed, VNIR = Visible & Near InfraRed, FRP = Fire Radiative Power).
Figure 2. The common steps characterizing a fire-detection algorithm using remote sensing data (IR = InfraRed, VNIR = Visible & Near InfraRed, FRP = Fire Radiative Power).
Remotesensing 10 00741 g002
Figure 4. Satellite FIre DEtection (SFIDE) algorithm prospect for diurnal data analysis (LST = land Surface Temperature).
Figure 4. Satellite FIre DEtection (SFIDE) algorithm prospect for diurnal data analysis (LST = land Surface Temperature).
Remotesensing 10 00741 g004
Figure 5. Sea-mask: black regions represent occulted pixels.
Figure 5. Sea-mask: black regions represent occulted pixels.
Remotesensing 10 00741 g005
Figure 6. SFIDE algorithm prospect for nocturnal data analysis.
Figure 6. SFIDE algorithm prospect for nocturnal data analysis.
Remotesensing 10 00741 g006
Figure 7. Comparison between SFIDE detected hot-spot and Corpo Forestale e di Vigilanza Ambientale (CFVA) measurements on 2014: 2 July (a); 3 July (b); 8 July (c) and 5 August (d).
Figure 7. Comparison between SFIDE detected hot-spot and Corpo Forestale e di Vigilanza Ambientale (CFVA) measurements on 2014: 2 July (a); 3 July (b); 8 July (c) and 5 August (d).
Remotesensing 10 00741 g007aRemotesensing 10 00741 g007b
Figure 8. Comparison of the behavior of two fires that occurred on 20 August 2017 through the temporal variation of the of Fire Radiative Power (FRP). The fire behavior in the two cases is very different. The evolution of the fire shown in the left plot seems much more “regular” than the one shown in the right plot.
Figure 8. Comparison of the behavior of two fires that occurred on 20 August 2017 through the temporal variation of the of Fire Radiative Power (FRP). The fire behavior in the two cases is very different. The evolution of the fire shown in the left plot seems much more “regular” than the one shown in the right plot.
Remotesensing 10 00741 g008
Table 1. Bands of SEVIRI radiometer (www.esa.int).
Table 1. Bands of SEVIRI radiometer (www.esa.int).
ChannelSpectral BandWavelength (μm)Central Wavelength (μm)Equator Resolution (km)
1VIS 0.60.56–0.710.6353
2VIS 0.80.74–0.880.813
3IR 1.61.50–1.781.643
4IR 3.93.48–4.363.923
5WV 6.25.35–7.156.253
6WV 7.36.85–7.857.353
7IR 8.78.30–9.108.703
8IR 9.79.38–9.949.663
9IR 10.89.80–11.8010.803
10IR 12.011.00–13.0012.003
11IR 13.412.40–14.4013.403
12HRV0.60–0.900.751
Table 2. SFIDE omission/commission matrix.
Table 2. SFIDE omission/commission matrix.
Detected Fire EventsOmitted Fire EventsDetected Hot-SpotFalse Alarms
454 (8.9%)46432 (6.9%)

Share and Cite

MDPI and ACS Style

Di Biase, V.; Laneve, G. Geostationary Sensor Based Forest Fire Detection and Monitoring: An Improved Version of the SFIDE Algorithm. Remote Sens. 2018, 10, 741. https://doi.org/10.3390/rs10050741

AMA Style

Di Biase V, Laneve G. Geostationary Sensor Based Forest Fire Detection and Monitoring: An Improved Version of the SFIDE Algorithm. Remote Sensing. 2018; 10(5):741. https://doi.org/10.3390/rs10050741

Chicago/Turabian Style

Di Biase, Valeria, and Giovanni Laneve. 2018. "Geostationary Sensor Based Forest Fire Detection and Monitoring: An Improved Version of the SFIDE Algorithm" Remote Sensing 10, no. 5: 741. https://doi.org/10.3390/rs10050741

APA Style

Di Biase, V., & Laneve, G. (2018). Geostationary Sensor Based Forest Fire Detection and Monitoring: An Improved Version of the SFIDE Algorithm. Remote Sensing, 10(5), 741. https://doi.org/10.3390/rs10050741

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop