# A Fuzzy Inference System Using Gaussian Distribution Curves for Forest Fire Risk Estimation

Lazaros Iliadis, Stergios Skopianos, Stavros Tachos, Stefanos Spartalis

► To cite this version:

Lazaros Iliadis, Stergios Skopianos, Stavros Tachos, Stefanos Spartalis. A Fuzzy Inference System Using Gaussian Distribution Curves for Forest Fire Risk Estimation. 6th IFIP WG 12.5 International Conference on Artificial Intelligence Applications and Innovations (AIAI), Oct 2010, Larnaca, Cyprus. pp.376-386, 10.1007/978-3-642-16239-8_49. hal-01060639

HAL Id: hal-01060639

https://inria.hal.science/hal-01060639

Submitted on 17 Nov 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A Fuzzy Inference System using Gaussian distribution curves for forest fire risk estimation

Lazaros Iliadis¹, Stergios Skopianos¹, Stavros Tachos², Stefanos Spartalis³

¹Democritus University of Thrace, Pandazidou 193 str., Orestiada, Greece

²Aristotle University of Thessaloniki, Department of Informatics

³Democritus University of Thrace, Xanthi, Greece

Lazaros Iliadis: liliadis@fmenr.duth.gr

**Abstract.** This paper describes the development of a fuzzy inference system under the MATLAB platform. The system uses three distinct Gaussian distribution fuzzy membership functions in order to estimate the partial and the overall risk indices due to wild fires in the southern part of Greece. The behavior of each curve has been investigated in order to determine which one fits better for the specific problem and for the specific areas. Regardless the characteristics of each function, the risky areas have been spotted from 1984 till 2007. The results have shown a reliable performance over time and they encourage its wider use in the near future.

**Keywords.** Fuzzy algebra, Gaussian distribution curves, Fuzzy Inference System, wildfire risk

## 1 Introduction

Every year, forest fires destroy important forest areas, with enormous consequences. The problem of forest fires occurs mainly in the Mediterranean countries, in Australia, in Canada and the United States [6]. Especially in Greece due to the lack of a proper cadastral survey, the problem of forest fires is related to political elections and also to the cycles of drought [9][10]. More detailed analysis on the problem of forest fires in Greece can be found in [14][16][17].

Fuzzy logic systems are quite flexible and they are used for several purposes such as modeling, forecasting and classification [8]. This research effort aims in the development of a fuzzy inference system, that estimates an integrated fuzzy forest fire risk index offering an overall view of the problem. This is achieved by employing three distinct fuzzy membership functions (of the same family) namely the Gaussian, the Sigmoid and the S-shaped one. The application of fuzzy logic enhances the ability to model proper linguistics that are used in the consequent part of the developed rule-based system. The application area concerns the Greek forest departments of Attica and Peloponnesus which are located in the southern part of the country. The whole research has been performed for the period 1984 to 2007 and the estimation of the integrated risk index has been achieved by considering the annual number of forest fire breakouts and the annual burned forest area for each forest department.

### 1.1 Literature Review

In order for the appropriate authorities to adopt efficient fire prevention and protection measures, the assistance of modern technology is required towards the evaluation of forest fire risk [6][8][13]. Deeming et al [4], developed the national fire danger rating system (NFDRS), which is still used in the USA after some improvements. It is a system that produces the daily ignition component index and also the man caused risk index. There have been many studies in the United States, in Australia, in New Zealand and in Europe for the long term evaluation of fire risk. A system that evaluates daily as well as long term fire danger is EFFIS (European Forest Fire Information System). It has been developed by the Institute for Environmental and Sustainability Land Management which resides in Italy [1]. In reality EFFIS is an attempt towards the creation of a common European indicator for fire danger. In Greece, the lab of forest informatics of the Democritus University of Thrace has developed decision support systems that are employing heuristic reasoning and fuzzy logic (fuzzy expected intervals) in order to estimate wild fire risk for the following year. The heuristic reasoning is based on the effect of political elections and drought cycles (following a periodicity of 3 to 5 years) in the number of forest fires in Greece [7].
# 2 Intelligent Information Systems

## 2.1 Fuzzy Algebra

Fuzzy logic [18],[2],[11], allows the definition of real world concepts, such as “area with high vulnerability to forest fires” in the form of fuzzy sets. Every member of the universe of discourse belongs to a fuzzy set with a different degree of membership, ranging from 0 to 1. The degree of membership can be estimated by the use of proper membership functions (MF). Fuzzy sets correspond to real life linguistics. Fuzzy logic is an human knowledge embodying tool through operational algorithms [6][11]. In fuzzy algebra there are many kinds of membership functions but the most characteristic are the triangular, the trapezoidal ones and the family of the sigmoid functions. Sigmoid MF have been used in several research efforts in order to develop fuzzy inference systems [3]. This research effort has been implemented by applying three distinct fuzzy functions that all belong to the sigmoid family.

### 2.1.1 The functions of the Gaussian distribution curve family

The fuzzy toolbox of the Matlab platform has been used in order to obtain the membership degrees for each forest department, given by three MF that belong to the Gaussian distribution curve, namely the Gaussian, the Sigmoid, and the S-shaped one. The Gaussian fuzzy membership function (GFMF) is symmetrical and it depends on two parameters c and σ which represent the center and the width of a fuzzy set [15]. More details can be seen in the following figure 2. Gaussian fuzzy membership functions are quite popular in the fuzzy logic literature, as they are the basis for the connection between fuzzy systems and radial basis function neural networks. A typical form of a GFMF can be seen in the following equation 1 [7].

$$f(x, \sigma, c) = e^{\frac{-(x-c)^2}{2\sigma^2}} \qquad (1)$$

The Sigmoid fuzzy membership function (SFMF) is given by the following equation 2. It is a continuous monotonic mapping of the input into a value between 0.0 and 1.0.

$$f(x;a,c) = \frac{1}{1 + e^{-a(x-c)}} \qquad (2)$$

In the case of the SFMF the parameters a, c define the schema and the position of the curve [7]. Depending on the sign of the parameter a, it is inherently open to the right or to the left, and thus it is appropriate for representing concepts such as "very large" or "very negative". The S-shaped membership function (SSMF) is a spline-based curve and it is a special version of the Sigmoid one. The following table 1, presents the commands that were employed in MATLAB in order to perform the membership functions.

<table><caption>**Table 1.** Sample T-Norms</caption><thead><tr><td>Fuzzy membership function</td><td>MATLAB's implementation command</td></tr></thead><tbody><tr><td>GFMF</td><td>y = Gaussmmf (x, [sig c])</td></tr><tr><td>SFMF</td><td>y = Sigmf (x, [a c])</td></tr><tr><td>SSMF</td><td>y = Smf (x, [a b])</td></tr></tbody></table>

It should be specified that the parameters a and b locate the extremes of the sloped portion of the curve. The following figure 1 shows the actual shape of the three sigmoid membership function that have been applied.
<figure><img src="image_2.png" /><figcaption>Fig. 1. Graphics of the three employed membership functions</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>Fig. 2. The center and the width of a Gaussian membership function</figcaption></figure>

2.1.2 Fuzzy Aggregation Operators

According to fuzzy algebra, every element of the universe of discourse can belong to a fuzzy set $\tilde{A}$ and simultaneously to its complement $\tilde{A}C$ with a degree of membership complementary to 1. T-Norms are the fuzzy conjunction operators between fuzzy sets. Fuzzy algebra contains several T-Norms, each one of them offering a different approach of the fuzzy AND, operation [7],[11]. The following table 2 below contains equations 3 and 4 and it presents a small sample of the most common T-Norms. They have a wide range of applications including the evaluation of Artificial Neural Networks [5].

<table><caption>Table 2. Sample T-Norms</caption><thead><tr><th>T-Norm name</th><th>T-Norm equation</th></tr></thead><tbody><tr><td>Min (3)</td><td>$\tilde{A} \wedge \tilde{B} = \{ (x, \mu_{\tilde{A} \wedge \tilde{B}}(x)) | \mu_{\tilde{A} \wedge \tilde{B}}(x) = \mu_{\tilde{A}}(x) \wedge \mu_{\tilde{B}}(x) = \min(\mu_{\tilde{A}}(x), \mu_{\tilde{B}}(x)) \}</td></tr><tr><td>Drastic Product (4)</td><td>$\tilde{A} \wedge \tilde{B} = \min\{\mu_{\tilde{A}}(x), \mu_{\tilde{B}}(x)\}$ if $\max\{\mu_{\tilde{A}}(x), \mu_{\tilde{B}}(x)\}=1$ else $\tilde{A} \wedge \tilde{B} = 0$</td></tr></tbody></table>
The choice of the appropriate norm depends on the perspective under which the problem is examined. For example, if the developer needs an optimistic norm that will assign the least potential overall degree of risk, then the min T-norm should be applied, whereas if the developer wishes to distinguish the cases that are characterized by extreme values for one or more risk parameters, the drastic product is the proper norm to be employed.

## 3 Developing the fuzzy inference system

The risk estimation is based on historical data that refer to the annual forest fire frequency (FFF) and to the annual burned area (ABA) in each area under study. The fuzzy inference system has been developed in the MatLab integrated environment and more specifically by writing MatLab code in the development environment and also by the use of the fuzzy toolbox for the construction of the fuzzy membership functions and for the Ruleset construction. The primitive data that have been used as input have been initially stored in Microsoft Excel sheets and they are input to the fuzzy inference system by employing

commands of the type: $\vec{C} = xlsread('data1.xls','1990','B2:B89');$ This sample command imports data concerning burned areas for 1990 into a MATLAB vector named C.

Several commands of this type have been written in order to store each column containing the annual FFF and the ABA in proper one dimensional Matrices (vectors). All of these commands have been stored in a MatLab programming file called *Input_Data.m*. The benefit from the creation of this file is that it runs just by typing the name of the file in the command line. For each one of the two input variables three fuzzy sets have been formed corresponding to proper Linguistics. The following six fuzzy sets have been formed, namely:

$$ \tilde{A}_1 = \left\{ \left( \vec{X}_i \mu_i \right) / i = 1..n..\text{"low risky due to FFF"} \right\} \quad \tilde{B}_1 = \left\{ \left( \vec{X}_i \mu_i \right) / i = 1..n..\text{"risky due to FFF"} \right\} $$

$$ \tilde{C}_1 = \left\{ \left( \vec{X}_i \mu_i \right) / i = 1..n..\text{"extremely risky due to FFF"} \right\} \quad \tilde{A}_2 = \left\{ \left( \vec{X}_i \mu_i \right) / i = 1..n..\text{"low risky due to ABA"} \right\} $$

$$ \tilde{B}_2 = \left\{ \left( \vec{X}_i \mu_i \right) / i = 1..n..\text{"risky due to ABA"} \right\} \quad \tilde{C}_2 = \left\{ \left( \vec{X}_i \mu_i \right) / i = 1..n..\text{"extremely risky due to ABA"} \right\} $$

where n is the number of case under study. Each forest department belongs to all six fuzzy sets with a distinct degree of membership. Of course the characteristic linguistic that describes a case is the one with the highest degree of membership. This approach not only allows the proper definition of linguistics that show clearly the status of an area due to one factor, but it also enables the distinction between two cases of the same linguistic. For example if two forest departments are both high risky, their degree of membership determines which one is more risky than the other.

During the design of the fuzzy inference system, the developer has to input the maximum and the minimum annual values ever recorded in Greece for the forest fire frequency and for the total burned area. Then the fuzzy toolbox determines the values of the σ and c parameters of the Gaussian function and the parameter values α, c, used by the Sigmoid and the S-shaped membership functions. For the output of the system also three fuzzy sets corresponding to three linguistics have been defined namely:

$$ \tilde{A}_3 = \left\{ \left( \vec{X}_i \mu_i \right) / i = 1..n..\text{"low risky due to integrated risk"} \right\} \quad \tilde{B}_3 = \left\{ \left( \vec{X}_i \mu_i \right) / i = 1..n..\text{"risky due to integrated risk"} \right\} $$

$$ \tilde{C}_3 = \left\{ \left( \vec{X}_i \mu_i \right) / i = 1..n..\text{"extremely risky due to integrated risk"} \right\} $$

After the storage of the data in vectors and the definitions of the fuzzy membership functions the Ruleset has been constructed. Figure 3, presents the design of the Ruleset. Several approaches have been proposed lately in order to design decision tables and decision trees [12]. In this case, due to the small number of parameters a rather heuristic approach has been followed in order to minimize the number of required rules as much as possible and to produce a quite rational system.
<table><caption>Fig. 3. A sample of the Ruleset of the fuzzy inference system.</caption><thead><tr><th>forest fire frequency</th><th>LR</th><th>R</th><th>LR</th><th>R</th><th>ER</th><th>R</th><th>ER</th><th>ER</th><th>LR</th></tr></thead><tbody><tr><th>annual burned area</th><td>LR</td><td>LR</td><td>R</td><td>R</td><td>R</td><td>ER</td><td>ER</td><td>LR</td><td>ER</td></tr><tr><th>OUTPUT</th><td>LR</td><td>LR</td><td>R</td><td>R</td><td>R</td><td>ER</td><td>ER</td><td>R</td><td>R</td></tr></tbody><tfoot><tr><td colspan="10" style="text-align: center;">The RuleSet</td></tr></tfoot></table>

The unification of the partial risk degrees due to the number of forest fire frequencies and to the total burned area to an integrated risk index has been performed by the system by the use of the minimum T-Norm.

Three distinct systems have been developed, each one using a specific sigmoid membership function and they have been saved as a *ffriskGaus.fis*, *ffriskSigm.fis* and *ffriskSm.fis* using the Gaussian, the Sigmoid and the S-shaped membership functions respectively. Each one of the developed fuzzy inference systems can be executed from the Matlab command line by using the commands *readfis* and *evalfis*.

For example the commands: *ffrisk = readfis('ffriskGaus.fis');* and *intergrated_risk = evalfis(ffriskGaus, fffrisk)*; have been used to execute the Gaussian system and they store the output (unified degrees of membership taking values from 0 to 1) into a proper matrix named *integrated_risk*. The results have been exported to sheets of Microsoft Excel again by using the MatLab command *xlswrite*.

## 4 Results - Discussion

The system has been tested and evaluated with actual data related to forest fire incidents in Greece for the period 1984-2007 which is a quite long time interval.

<table><caption>Table 3. PRI based on the FFF using the Gaussian and S shaped functions</caption><thead><tr><th>Year</th><th>forest department</th><th>GFMF Extremely risky index FFF</th><th>SFMF Extremely risky index FFF</th><th>Year</th><th>forest department</th><th>GFMF Extremely risky index FFF</th><th>SFMF Extremely risky index FFF</th></tr></thead><tbody><tr><td rowspan="3">1985</td><td>Pyrgos</td><td>0.9737</td><td>1</td><td rowspan="3">1999</td><td>Hleia</td><td>1</td><td>1</td></tr><tr><td>Olympia</td><td>0.9538</td><td>0.966</td><td>Messinia</td><td>1</td><td>1</td></tr><tr><td>Patra</td><td>0.8275</td><td></td><td>Argolida</td><td>0.9967</td><td>1</td></tr><tr><td rowspan="4">1986</td><td>Olympia</td><td>0.9802</td><td>1</td><td rowspan="4">2000</td><td>Messinia</td><td>1</td><td>1</td></tr><tr><td>Pyrgos</td><td>0.995</td><td>1</td><td>Korinthos</td><td>1</td><td>1</td></tr><tr><td>Patra</td><td>0.956</td><td>0.9938</td><td>Hleia</td><td>0.9994</td><td>1</td></tr><tr><td>Korinthos</td><td>0.9231</td><td>0.9753</td><td>Lakonia</td><td>0.9508</td><td>1</td></tr><tr><td rowspan="2">1987</td><td>Pyrgos</td><td>1</td><td>1</td><td rowspan="2">2001</td><td>Messinia</td><td>1</td><td>1</td></tr><tr><td>Olympia</td><td>0.9692</td><td>1</td><td>Korinthos</td><td>1</td><td>1</td></tr><tr><td rowspan="2">1988</td><td>Olympia</td><td>0.9523</td><td>1</td><td rowspan="2">2003</td><td>Arkadia</td><td>1</td><td>1</td></tr><tr><td>Argolida</td><td>0.9692</td><td>0.962</td><td>Lakonia</td><td>0.9569</td><td>1</td></tr><tr><td rowspan="4">1989</td><td>Pyrgos</td><td>0.9802</td><td>1</td><td rowspan="4">2004</td><td>Messinia</td><td>1</td><td>1</td></tr><tr><td>Olympia</td><td>0.9802</td><td>1</td><td>Hleia</td><td>1</td><td>1</td></tr><tr><td>Tripolh</td><td>0.9802</td><td>0.9988</td><td>Arkadia</td><td>0.9999</td><td>1</td></tr><tr><td>Patra</td><td>0.956</td><td>1</td><td>Korinthos</td><td>0.9928</td><td>1</td></tr><tr><td rowspan="2">1990</td><td>Pyrgos</td><td>0.9756</td><td>1</td><td rowspan="2">2004</td><td>Messinia</td><td>1</td><td>1</td></tr><tr><td>Olympia</td><td>0.9756</td><td>0.9882</td><td>Hleia</td><td>1</td><td>1</td></tr><tr><td rowspan="3">1991</td><td>Argolida</td><td>0.9692</td><td>0.9644</td><td rowspan="3">2004</td><td>Korinthos</td><td>0.9889</td><td>0.9828</td></tr><tr><td>Korinthos</td><td>0.9692</td><td>0.9644</td><td>Arkadia</td><td>0.9595</td><td>0.946</td></tr><tr><td>Olympia</td><td>0.9321</td><td>1</td><td>Lakonia</td><td>0.9406</td><td>0.92</td></tr><tr><td>1992</td><td>Korinthos</td><td>0.9978</td><td>1</td><td></td><td>Argolida</td><td>0.9231</td><td>1</td></tr></tbody></table>
<table><tbody><tr><td></td><td></td><td>Olympia</td><td>0.9802</td><td>1</td><td>2005</td><td>Messinia</td><td>1</td><td>1</td></tr><tr><td>1993</td><td>Kalamata</td><td>0.9396</td><td>1</td><td></td><td></td><td>Hleia</td><td>1</td><td>1</td></tr><tr><td>1994</td><td>Argolida</td><td>0.9287</td><td>1</td><td></td><td></td><td>Korinthos</td><td>0.9996</td><td>0.9929</td></tr><tr><td></td><td>Kalamata</td><td>0.9737</td><td>0.9444</td><td>2006</td><td></td><td>Hleia</td><td>1</td><td>1</td></tr><tr><td>1995</td><td>Kalamata</td><td>1</td><td>1</td><td></td><td></td><td>Messinia</td><td>1</td><td>1</td></tr><tr><td></td><td>Amaliada</td><td>0.9577</td><td>1</td><td></td><td></td><td>Arkadia</td><td>0.9928</td><td>1</td></tr><tr><td>1996</td><td>Kalamata</td><td>1</td><td>1</td><td></td><td></td><td>Argolida</td><td>0.9851</td><td>0.9812</td></tr><tr><td></td><td>Pyrgos</td><td>0.96</td><td>0.933</td><td>2007</td><td></td><td>Hleia</td><td>1</td><td>1</td></tr><tr><td></td><td>Pendeli</td><td></td><td>1</td><td></td><td></td><td>Korinthos</td><td>1</td><td>1</td></tr><tr><td></td><td>Parnitha</td><td></td><td>1</td><td></td><td></td><td>Patra</td><td>0.96</td><td>1</td></tr><tr><td>1997</td><td>Kalamata</td><td>1</td><td>1</td><td></td><td></td><td>Messinia</td><td>0.9675</td><td></td></tr><tr><td></td><td>Olympia</td><td>1</td><td>1</td><td></td><td></td><td>Arkadia</td><td>0.8</td><td>1</td></tr><tr><td></td><td>Pyrgos</td><td>0.9596</td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

The above Table 3 presents clearly the obtained partial risk indices based on the number of forest fire incidents' frequencies, with the use of both Gaussian and S-shaped membership functions only for the cases of the Extremely risky forest departments.

<table><caption>**Table 4.** PRI for the ABA and overall risk indices for the combination of GFMF and Min (1985-2007)</caption><thead><tr><th>Year</th><th>Forest department</th><th>GFMF Extremely risky Partial index ABA</th><th>SSMF F Extremely risky Partial index ABA</th><th>OVE RALL risk index GFMF and Min Norm</th><th>Year</th><th>Forest department</th><th>GFMF Extremely risky Partial index ABA</th><th>SSMF Extremely risky Partial index ABA</th><th>OVERALL risk index GFMF and Min Norm</th></tr></thead><tbody><tr><td rowspan="3">1984</td><td>Kapandriti</td><td>1</td><td>1</td><td></td><td rowspan="3">1994</td><td>Patra</td><td>1</td><td>1</td><td>XR 0.181</td></tr><tr><td>Aigaleo</td><td>0.9817</td><td>1</td><td></td><td>Pendelh</td><td>1</td><td>1</td><td>XR 0.606</td></tr><tr><td>Korinthos</td><td>0.9247</td><td>1</td><td>R</td><td>Gytheio</td><td>1</td><td>1</td><td>R 0.323</td></tr><tr><td rowspan="6">1985</td><td>Aigaleo</td><td>1</td><td>1</td><td>0.150</td><td rowspan="6">1995</td><td>Patra</td><td>1</td><td>1</td><td>R 0.969</td></tr><tr><td>Megara</td><td>1</td><td>1</td><td>0.110</td><td>Gytheio</td><td>1</td><td>1</td><td>XR 0.369</td></tr><tr><td>Patra</td><td>0.9999</td><td>1</td><td>0.827</td><td>Molaoi</td><td>0.9584</td><td>1</td><td>0.310</td></tr><tr><td>Lavrio</td><td>0.9903</td><td>1</td><td>0.582</td><td>Kalama ta</td><td>0.9525</td><td>0.9191</td><td>XR 0.952</td></tr><tr><td>Amaliada</td><td>0.9522</td><td>1</td><td>0.457</td><td>Molaoi</td><td>1</td><td>1</td><td>XR 0.230</td></tr><tr><td>Amaliada</td><td>0.9522</td><td>1</td><td>0.923</td><td>Patra</td><td>0.9718</td><td>1</td><td>0.234</td></tr><tr><td rowspan="4">1986</td><td>Korinthos</td><td>1</td><td>1</td><td>R</td><td rowspan="4">1997</td><td>Korinth os</td><td>0.9969</td><td>0.9904</td><td>XR 0.774</td></tr><tr><td>Molaoi</td><td>0.9867</td><td>1</td><td>0.946</td><td>Olympi a</td><td>1</td><td>1</td><td>XR 1</td></tr><tr><td>Olymbia</td><td>0.9356</td><td>0.9249</td><td>0.935</td><td>Gytheio</td><td>0.9886</td><td>1</td><td>R 0.946</td></tr><tr><td>Pyrgos</td><td>0.9097</td><td>0.8943</td><td>0.909</td><td>Korinth os</td><td>1</td><td>1</td><td>XR 0.700</td></tr><tr><td rowspan="5">1987</td><td>Patra</td><td>1</td><td>1</td><td>0.681</td><td rowspan="5">1998</td><td>Pendeli</td><td>1</td><td>1</td><td>XR 0.290</td></tr><tr><td>Kalamata</td><td>0.9868</td><td>1</td><td>0.754</td><td>Kapand riti</td><td>0.95</td><td>1</td><td>R 0.950</td></tr><tr><td>Moaloi</td><td>1</td><td>1</td><td>0.457</td><td></td><td></td><td></td><td>XR 0.910</td></tr><tr><td>Argolida</td><td>0.9712</td><td>0.9675</td><td>0.339</td><td>Lakonia</td><td>1</td><td>1</td><td>XR 0.967</td></tr><tr><td>Kalamata</td><td>1</td><td>1</td><td>0.216</td><td>Messhn ia</td><td>0.9679</td><td>1</td><td></td></tr></tbody></table>
<table><tbody><tr><td></td><td>Tripolh</td><td>0.9024</td><td>1</td><td>XR<br>0.153<br>XR</td><td>Argolid<br>a</td><td>0.9998</td><td>1</td><td>XR<br>0.996</td></tr><tr><td rowspan="4">1989</td><td>Pyrgos</td><td>1</td><td>1</td><td>0.980<br>XR</td><td>2000</td><td>Korinth<br>os</td><td>1</td><td>1</td><td>XR 1<br>XR</td></tr><tr><td>Xylokastro</td><td>0.9201</td><td>1</td><td>0.920<br>XR</td><td>Xylokat<br>sro</td><td>1</td><td>1</td><td>0.433<br>XR</td></tr><tr><td>Patra</td><td>1</td><td>1</td><td>0.956<br>XR</td><td>Arkadia</td><td>0.8479</td><td>1</td><td>0.847<br>XR</td></tr><tr><td>Sparti</td><td>0.9978</td><td>0.9835</td><td>0.882<br>XR</td><td>Kapand<br>riti</td><td>1</td><td>1</td><td>0.244<br>XR</td></tr><tr><td rowspan="3">1990</td><td>Kapandriti</td><td>1</td><td>1</td><td>0.657<br>R</td><td>Lakonia</td><td>1</td><td>1</td><td>0.965<br>XR</td></tr><tr><td>Korinthos</td><td>1</td><td>1</td><td>0.324<br>XR</td><td>Arkadia</td><td>0.9708</td><td>1</td><td>0.970<br>XR</td></tr><tr><td>Molaoi</td><td>0.9737</td><td>1</td><td>0.673<br>XR</td><td>Arkadia</td><td>1</td><td>1</td><td>0.999</td></tr><tr><td rowspan="5">1991</td><td>Kapandriti</td><td>1</td><td>1</td><td>0.135<br>XR</td><td>Hleia</td><td>1</td><td>1</td><td>XR 1<br>XR</td></tr><tr><td>Korinthos</td><td>0.8667</td><td>1</td><td>0.866<br>XR</td><td>Argolid<br>a</td><td>0.9166</td><td>1</td><td>0.928<br>XR</td></tr><tr><td>Argolida</td><td>0.9704</td><td>1</td><td>0.969<br>XR</td><td>Lakonia</td><td>1</td><td>1</td><td>0.940<br>XR</td></tr><tr><td>Olympia</td><td>0.9584</td><td>0.9349</td><td>0.932<br>XR</td><td>Lavrio</td><td>1</td><td>1</td><td>0.740<br>XR</td></tr><tr><td>Pendeli</td><td>1</td><td>1</td><td>0.766<br>R</td><td>Hleia</td><td>0.9136</td><td>1</td><td>0.913<br>XR</td></tr><tr><td rowspan="6">1992</td><td>Kapandriti</td><td>1</td><td>1</td><td>0.1543<br>R</td><td>2005</td><td>Pendeli</td><td>1</td><td>1</td><td>0.432<br>XR</td></tr><tr><td>Sparti</td><td>0.9065</td><td>1</td><td>0.2226<br>XR</td><td>Lakonia</td><td>1</td><td>1</td><td>0.8795<br>XR</td></tr><tr><td>Lavrio</td><td>0.9794</td><td>1</td><td>0.657<br>XR</td><td>Messhn<br>ia</td><td>0.9393</td><td>1</td><td>0.939<br>XR</td></tr><tr><td>Kalamata</td><td>0.9850</td><td>0.9771</td><td>0.882<br>XR</td><td>Arkadia</td><td>0.9761</td><td>1</td><td>0.841<br>XR</td></tr><tr><td>Argolida</td><td>0.9845</td><td>0.9764</td><td>0.748<br>XR</td><td>Lakonia</td><td>1</td><td>1</td><td>0.718<br>XR</td></tr><tr><td>Kalamata</td><td>1</td><td>1</td><td>0.939<br>R</td><td>Xylokat<br>sro</td><td>0.8276</td><td>1</td><td>0.102<br>XR</td></tr><tr><td rowspan="4">1993</td><td>Tripolh</td><td>1</td><td>1</td><td>0.923<br>XR</td><td>Arkadia</td><td>0.9874</td><td>0.9892</td><td>0.992</td></tr><tr><td>Molaoi</td><td>0.8604</td><td>1</td><td>0.439<br>XR</td><td>Hleia</td><td>1</td><td>1</td><td>XR 1<br>XR</td></tr><tr><td>Xylokastro</td><td>0.9258</td><td>1</td><td>0.354</td><td>Arkadia</td><td>1</td><td>1</td><td>0.798<br>XR</td></tr><tr><td></td><td></td><td></td><td></td><td>Korinth<br>os</td><td>0.9337</td><td>1</td><td>0.933</td></tr></tbody></table>

The determination of the most risky areas based on partial risk indices (due to forest fire frequencies or due to the annual burned area) is similar regardless the use of the Gaussian, the Sigmoid or the S-shaped membership functions. However the S-shaped function assigns the highest risk value more often whereas the Gaussian one differentiates the forest departments more clearly by assigning them with a wide range of risk indices.

Table 4 presents the partial risk indices for the extremely risky areas according to the annual burned area, in the cases of the Gaussian and of the S-shaped membership functions. The estimated overall risk indices' vector (obtained as the output of the fuzzy inference system) for the Gaussian membership function and for the Min T-Norm is also presented clearly in the above Table 4. It contains the most risky forest departments based on partial risk indices and also on the overall risk index (considering risk due to both FFF and ABA). It estimates their corresponding linguistics (risky R, or extremely risky XR) only for the Gaussian function due to the fact that it offers a more rational production of risk indices, whereas the other two functions assign the highest index with the value 1 for too many cases and thus there is no proper distinction between the cases. The determination of the unified overall risk indices has been performed by the fuzzy inference system based on the Ruleset which is presented in figure 3.
# 5 Conclusions

The annual ‘extremely risky’ PRIs and the overall integrated ‘extremely risky’ indices obtained by the system for each year were compared to the actually most risky areas based on the number of FFF and on the ABA for the following year in order to determine the validity of the system. As far as the Gaussmf curve is concerned, for the parameter FFF the percentages of successful risk characterization reaches as high as 70%. The same function for the parameter of the annual burned area offers a percentage of success equal to 48%. The Sigmoid function (Sigmf), for the parameter FFF has also a percentage of successful risk characterizations as high as 70%, but for the burned area the percentage is 51%. Finally Smf has a percentage of success as high as 74% for the number of forest fires and 46% for the burned area. During the period 1985-2007 the overall integrated risk index obtained by employing the min T-Norm offers an average accuracy of 65%, 54% and 52% for the cases of the Gaussian, Sigmoid and S-shaped membership functions respectively. The Gaussian function achieves the highest percentages of success towards overall risk estimation, whereas the other two functions give also satisfying results. The system will be extended by applying weights in the unification of the risk indices according to the importance of each factor.

# References

1. Ayanz, J.S.M., Barbosa P.M., Schmuck G., Liberta G., 2003 The European Forest Fire Information System (EFFIS). European Commission. JIR Institute for Environment and sustainability. URL: http://natural-hazards.jrc.it/documents/fires/2003-publications/EFFIS-earsel-paper.pdf

2. Cox E., 2005 Fuzzy Modeling and Genetic Algorithms for Data Mining and Exploration Elsevier Inc., USA

3. Duprey B., Taheri S., 2010 A Fuzzy Based Stability Index Using a Right Sigmoid Membership Function SAE International Journal of Commercial Vehicles March 2010 vol. 2 no. 2

4. Deeming, J., Burgan R., Cohen J., 1978 The National Fire Danger Rating System USDA

5. Iliadis L, S. Spartalis, S. Tachos, 2008 "Application of fuzzy T-norms towards Artificial Neural Networks' evaluation: A case from wood industry" Journal Information Sciences (Informatics and Computer Science Intelligent Systems Applications) Volume 178, Issue 20, pp. 3828-3839

6. Iliadis L., 2005 A decision support system applying an integrated Fuzzy model for long - term forest fire risk estimation Environmental Modeling and Software Vol. 20, Issue 5, pp. 613-621

7. Iliadis L., 2007 Intelligent Information Systems in Risk Estimation Stamoulis Editors (Greek)

8. Iliadis L., Papastavrou A., Lefakis P., 2002 A computer-system that classifies the prefectures of Greece in forest fire risk zones using fuzzy sets. Journal Forest policy and Economics, vol 4, Issue 1, pp. 43-54

9. Kailidis, D. 1990. Forest fires. 3rd edition. Giahoudi-Giapouli editions. Thessaloniki, Greece. 510 p. (in Greek).

10. Kailidis, D. & Xanthopoulos, G. 1991. The forest fire problem in Greece. Aristotelian University of Thessaloniki, Greece, Department of Forestry and Natural Environment, Forest Protection Laboratory. No. 310 pp.

11. Kecman V., 2001 Learning and Soft Computing MIT Press

12. Kuncheva L.I and J.J. Rodriguez, An experimental study on Rotation Forest ensembles, Proc 7th International Workshop on Multiple Classifier Systems, MCS'07, Prague, Czech Republic, 2007, LNCS 4472, 459-468

13. Lin C., 2002 The Development, Systems, and Evaluation of Forest Fire Danger Rating: A Review. Taiwan Forestry Research Institute. URL: http://fire.tfri.gov.tw./DocRefl4.pdf

14. Markalas, S. & Pantelis, D. 1996. Forest fires in Greece in 1993. Aristotelian Univ. of Thessaloniki, Department of Forestry and the Natural Environment, Forest Protection Laboratory. No. 3, 40 pp. (in Greek).

15. Pieczynski A., Obuchowicz A., 2004 Application of the General Gaussian Membership Function for the Fuzzy Model Parameters Tunning. pp. 350-355 Artificial Intelligence and Soft Computing. Lecture Notes in Computer Science. Springer Berlin/Heidelberg

16. Xanthopoulos, G. 1988. Greek forest fires and property damage: A brief history. p. 199-200. In: Proceedings - Symposium and Workshop on "Protecting People and Homes from Wildfire in the Interior West", October 6-8, 1987, Missoula, Montana, USA. USDA For. Serv. Gen. Tech. Rep. INT-251, 213 pp.

17. Xanthopoulos, G. 1998. Forest fires in Greece: Past, present and future. Epikentra (Foundation for Political Research and Advanced Education). 6: 62-71. (Special issue titled "Greek forests: myths and reality - Proposals for the present and for the future). (in Greek).

18. Zadeh, L., A. 1968 Inform Control vol12 p.p. 94-102