Burnt-Area Quick Mapping Method with Synthetic Aperture Radar Data

: Forest and ﬁeld ﬁres have become a frequent phenomenon in recent years caused by human activities in Indonesia, affecting all forms of forest and ﬁeld cover. Forest ﬁre-degraded land is more prone to burn again, due to the nature of peatland in Kalimantan. Rapid mapping technology for burnt areas affected by forest ﬁres is needed to obtain this information. The use of remote sensing technology, in the form of synthetic aperture radar (SAR) images, combined with cloud computing processing speeds up data processing and is not affected by the existing cloud cover. The Quick-Mapping employed in this research provides faster mapping time, compared to the currently employed method, based on ﬁeld report data, to enable a better and more efﬁcient ﬁreﬁghting effort. The data processing is carried out using cloud computing, enabling the processing of large amounts of data. The processing is carried out starting with importing the data, preprocessing to classiﬁcation running, simultaneously, using the JavaScript programming language. The research classiﬁes the burnt area from backscatter patterns before and after the event in two measurements, namely the radar burn ratio (RBR) and the radar burn difference (RBD). The RBR is deﬁned as the average backscatter ratio at a certain polarization, while RBD is the difference between the average scattering conditions. The composite image for the classiﬁcation utilizes images from the RBR and RBD with co-polarized (VV) and cross-polarized (VH) backscatter. The burnt area difference is − 1.9 for VH and − 1.7 for VV, which indicates a lower backscatter, due to forest ﬁre. The classiﬁcation of the burnt area yields the best overall accuracy of 88.26% with a support vector machine and processing time of 1 h, compared to the currently 12 h needed to provide burnt area maps from ﬁeld observation data.


Introduction
Forest and land fires have become a frequent and human-caused phenomenon in Indonesia, affecting the biodiversity and economic loss [1,2]. For reference, during the great 2015 forest and land fires in Indonesia, around 2.6 million trees were burnt, costing the country around USD 16.1 billion or the equivalent of 1.9% of Indonesia's gross domestic product (GDP), in 2015 [2]. Most of the fire ignition is coming from the peat land, compared to the mineral soils and the estimated the total emissions were 1.2 billion CO2 equivalent, during the Indonesian forest fire event in 2015 [3]. The haze from the fire ignition caused acute respiratory infections for more than half a million people and also impacted the surrounding countries, such as Malaysia and Singapore [2,4]. The unsustainable conversion of land use practices, such as the use of fire in converting forest and peat to agricultural areas and the widespread use of fire in cultivation, are considered to be the main driving factors [5,6]. In 2019, another episode of forest and land fires hit in Indonesia [7]. Although the 2019, the Indonesian forest and land fires were not as devastating as in 2015, but at least about 296.942 hectares were burnt in 2019. Some previous studies stated that forest and land fire episodes always hit Indonesia every year, although the magnitude of the forest and land fire ignition is closely related to the El Niño event for example in 1997-1998 and 2015 [8][9][10].
Various impacts arise after forest and land fires, the most concerning issue about post-forest and land fire events is land degradation, due to the burning of the existing vegetation cover and fire impacts on the soil environment [11]. The effect of fire is not only directly associated with the destruction of vegetation cover but also alters the hydrological condition of the soil, i.e., reducing the resistance of soil to erosion and increasing the overland flow during storm events [12][13][14]. As the consequences of the post-fire effect to the soil environment is really harmful, a rapid assessment is needed to map the distribution of the post-fire land degradation area, widely known as a burnt scar area [15].
Remote sensing in particular, due the variety of its sensors, is known as a powerful tool to assess the burnt scars rapidly [15][16][17]. A few studies have explored the application of remotely-sensed data for mapping the burnt scars in Indonesia, and found that the radar sensor is most suitable, rather than other optical sensors [17][18][19]. The combination of using radar sensors with the state-of-the-art computational techniques, using cloud computing, can improve the speed of image processing and avoid the bias of cloud cover [20][21][22].
The usability of the Google Earth Engine as one of the cloud computing platforms for satellite imagery processing, has grown over the last decade, since 2011 [23]. GEE is host to many satellite imageries for over 40 years, from such satellites as Landsat, MODIS, the National Oceanographic and Atmospheric Administration Advanced-Very-High Resolution Radiometer (NOAA-AVHRR), Sentinel 1,2,3, and 5-P, and the Advanced Land Observing Satellite (ALOS) [20]. The Google Earth Engine (GEE) not only offers the parallelized processing of remote sensing data on a global scale, using the Google Cloud, but also provides several machine learning algorithms for improving the classification techniques of the imageries [20,24]. Another advantage using GEE for image processing, GEE was first created and optimized for satellite imagery processing that can process petabytes of imagery data, both on large geographical scales and in long temporal coverages, hence it is a great tool for analyzing regional, national, continental, and global-scale applications [20,25].
A machine learning algorithm, as a classifier feature in GEE, offers a huge application for mapping any earth surface, including forest environments [20]. The most applicable machine learning algorithm for mapping the forest environment, i.e., random forest, classification and regression trees (CART), and the support vector machine have been widely used [26][27][28][29][30]. Random forest formulates to construct assemblage-based decision trees which are independently sampled during the training process and improve the classification results, rather than a single decision tree model [27,31].
Finding the best algorithm for GEE to detect burnt areas in Indonesia has never been achieved before so that in this study, we develop a rapid assessment of the burnt area detection method using GEE, and evaluate the performance of the available machine learning algorithm in GEE, for detecting the post-fire burnt area in Indonesia, using radar sensors. This research aims to quickly map the burnt areas caused by forest fires, using cloud computing of SAR image data with SVM, RF, and CART classification methods being considered, and the best result is selected to provide the forest fire impact map for the effective disaster mitigation efforts. The proposed method provides a faster mapping time, compared to the currently employed method, based on field report data, to enable a better and more efficient firefighting effort.

Location and Data
The location of this research is Kalimantan island, with the focused area of interest (AoI) marked with the red square on the regencies of Kapuas, Pulang Pisau, East Barito, South Barito, Barito Kuala, Hulu Sungai Utara, Hulu Sungai Tengah, Hulu Sungai Selatan, and Tapin, as shown in Figure 1.
The location of this research is Kalimantan island, with the focused area of interest (AoI) marked with the red square on the regencies of Kapuas, Pulang Pisau, East Barito, South Barito, Barito Kuala, Hulu Sungai Utara, Hulu Sungai Tengah, Hulu Sungai Selatan, and Tapin, as shown in Figure 1. The radar image is obtained from reflected microwave emissions originating from the radar platform. As the radar sensor moves along a trajectory, the area highlighted by the radar also moves along with the platform, known as the footprint, to form an image. The digital radar image consists of pixel dots representing the backscatter of a surface. The radar system generally uses wavelengths that are undisturbed by the interference from particles and vapor in the air (i.e., clouds and rain), and since the system is independent of the illumination (irradiation) from the sun or other sources, it can function in all weather conditions. The synthetic aperture radar (SAR) works by detecting the phase-change of the reflected signals caused by the movement of the platform, to obtain the surface image with a good resolution (i.e., visually discernible) [22,32]. Example of the difference between optical and SAR sensor is shown in Figure 2.
The data used in this study is Sentinel-1a (SAR) data, which has the advantage of not depending on weather conditions. This data processing uses cloud computing data processing, making it possible to process large amounts of data. The processing process, which consists of importing data, and preprocessing to classification, is executed in cloud computing, simultaneously, with commands using the JavaScript programming language and Python. The advantage of using this cloud computing method, compared to processing using a desktop, besides being faster in processing, the processor does not need high software specifications. When using desktop applications, raw data and output files from the processing results often fill the hard disk capacity, while using cloud computing, all files are stored in the cloud and can be downloaded, if needed. The cloud computing application that will be used in this research is the Google Earth Engine [20]. The data used in this research are The radar image is obtained from reflected microwave emissions originating from the radar platform. As the radar sensor moves along a trajectory, the area highlighted by the radar also moves along with the platform, known as the footprint, to form an image. The digital radar image consists of pixel dots representing the backscatter of a surface. The radar system generally uses wavelengths that are undisturbed by the interference from particles and vapor in the air (i.e., clouds and rain), and since the system is independent of the illumination (irradiation) from the sun or other sources, it can function in all weather conditions. The synthetic aperture radar (SAR) works by detecting the phase-change of the reflected signals caused by the movement of the platform, to obtain the surface image with a good resolution (i.e., visually discernible) [22,32]. Example of the difference between optical and SAR sensor is shown in Figure 2.
The data used in this study is Sentinel-1a (SAR) data, which has the advantage of not depending on weather conditions. This data processing uses cloud computing data processing, making it possible to process large amounts of data. The processing process, which consists of importing data, and preprocessing to classification, is executed in cloud computing, simultaneously, with commands using the JavaScript programming language and Python. The advantage of using this cloud computing method, compared to processing using a desktop, besides being faster in processing, the processor does not need high software specifications. When using desktop applications, raw data and output files from the processing results often fill the hard disk capacity, while using cloud computing, all files are stored in the cloud and can be downloaded, if needed. The cloud computing application that will be used in this research is the Google Earth Engine [20]. The data used in this research are the Sentinel 2A radar and Hotspot data, which comes from the SIPONGI system, used by the Ministry of Environment and Forestry [33]. The Sentinel data is used to detect the area affected, while the hotspot data is used to verify the burnt land. The data processing in this research uses cloud computing technology, by utilizing the Earth Engine. The processing is carried out by inputting a script in JavaScript language, to give commands to the server. The initial stage of the process is to form a composite image, as the basis for the classification process.
the Sentinel 2A radar and Hotspot data, which comes from the SIPONGI system, used by the Ministry of Environment and Forestry [33]. The Sentinel data is used to detect the area affected, while the hotspot data is used to verify the burnt land. The data processing in this research uses cloud computing technology, by utilizing the Earth Engine. The processing is carried out by inputting a script in JavaScript language, to give commands to the server. The initial stage of the process is to form a composite image, as the basis for the classification process. Figure 2. The advantages of the SAR data, compared to the optical sensors for the detection of the burn scar areas where visually, it is not covered by clouds and smoke.

Related Method
The support vector machine (SVM) algorithm allocates pixels into classes by means of maximizing the class separation from the training data, and labels the pixels, according to the nearest class [34]. The SVM principle is to search for the best classifier which gives the largest margin of accuracy [35]. For a training set = , , … , ∈ × , which consists of two classes that are denoted as = −1, +1 , the SVM aims to maximize the hyperplane margin while minimizing the error. The optimization is solved as follows: where αi is Lagrange multipliers and K(x, xi) is the kernel function, the decision function is written as follows: The SVM maps the original data into a high dimensional feature space using kernel functions and finds the optimal hyperplane that separates the true and false classes [36]. The SVM's most interesting property is its high capacity for generalization with relatively small numbers of training data points. Recent studies have suggested that the support vector machine (SVM) can provide good results for the hyperspectral remote sensing

Related Method
The support vector machine (SVM) algorithm allocates pixels into classes by means of maximizing the class separation from the training data, and labels the pixels, according to the nearest class [34]. The SVM principle is to search for the best classifier which gives the largest margin of accuracy [35]. For a training set x = {x 0 , x 1 , . . . , x N } ∈ R B×N , which consists of two classes that are denoted as L = {−1, +1}, the SVM aims to maximize the hyperplane margin while minimizing the error. The optimization is solved as follows: where α i is Lagrange multipliers and K(x, x i ) is the kernel function, the decision function is written as follows: The SVM maps the original data into a high dimensional feature space using kernel functions and finds the optimal hyperplane that separates the true and false classes [36]. The SVM's most interesting property is its high capacity for generalization with relatively small numbers of training data points. Recent studies have suggested that the support vector machine (SVM) can provide good results for the hyperspectral remote sensing classification and superior results have been reported, compared to the traditional remote sensing classification algorithms, such as the maximum likelihood (ML), k-nearest neighbor, and neural networks (NNs) [37]. In addition to the SVM-based categorical classification, there is also growing interest in the SVM regression for estimating the subpixel land cover proportions. The SVM with the radial basis function (RBF) kernel can achieve an excellent performance in classification/prediction tasks, due to the following facts: First, the RBF kernel can map a sample to a higher dimensional space, and the linear kernel function is essentially a special case of the RBF. Meanwhile, the kernel of the RBF and the sigmoid have a similar performance with certain parameters. Second, compared with the polynomial kernel function, only a few parameters need to be modulated by the RBF and the number of kernel function parameters directly influences the complexity of this function. Finally, the RBF kernel has fewer numerical difficulties [36].
The random forest method generates an ensemble of regression trees through the aggregated sampling of the training data. At each decision tree's node, a random selection parameter of the predictor variables is evaluated for their ability to split the training data into respective target classes, with the variable leading to the most homogeneous classification being selected [38]. Previous studies revealed that the random forest algorithm was acceptable to map the burnt scars area, by showing an overall accuracy >0.90 [39]. The decision tree-based algorithm and support vector machine was are also powerful to identify the burnt areas with an accuracy 0.85-0.91 and >0.99, respectively [40,41].
In line with the development of the field of artificial intelligence, image processing methods also develop by making use of artificial intelligence functions. Several artificial intelligence methods that are widely used in image processing are artificial neural networks or ANNs. The method that has recently begun to be applied in the studies of mapping flooding potential and vulnerability, is to use machine learning [22]. Neural network classification algorithms have long been used for remote sensing image classification. Many have suggested that these types of models are superior to traditional statistical classification approaches (i.e., maximum-likelihood classification), because they do not make assumptions about the nature of the data distribution, and the function is simply learned from training samples. Several neural network models are commonly applied.
Classification and regression trees (CART), are known as decision tree classifiers, which apply a multi-stage binary system to classify imageries [42]. The CART principle is based on the tree framework that has been widely used in remote sensing applications. The key principle is to determine the structure of a decision tree to select an input feature and threshold value at each splitting point. The pixels then are divided, according to the binary classification rule. Groups of pixels are divided, based on the tree growing and pruning parameters, until the optimal classification is found [43]. However, CART easily suffers from overfitting, if allowed to grow uncontrolled, to fit all of the training data. Therefore, the limiting process called pruning, can be employed to reduce the tree levels and increase the generalization ability of CART [37]. The large number of input features requires a significant amount of training pixels to generate a robust classification result.

Proposed Method
Forest fire mapping employed by Indonesian agencies mostly use optical imagery, such as Landsat-8 and Sentinel-2 [44][45][46], but extensive and persistent cloud cover in Indonesia requires other sensors which are unaffected by clouds [22]. The burnt area, during the 2019 Indonesia land and forest events, was generated from the Sentinel-1a synthetic aperture radar (SAR) satellite imageries and processed using random forest, classification and regression trees (CART), and the support vector machine.
One of the task of forest fire mapping is to map the fire perimeters and to predict the areas of potential fire hazards from the areas of vegetation regrowth after fires. The burn area index (BAI) highlights the burnt land between two bands of the spectrum, by emphasizing the difference of reflectance of the burnt vegetation in post-fire images. The index is computed from the spectral distance from each pixel to a reference spectral point where recently burnt areas converge, in which brighter pixels indicate burnt areas. The BAI is calculated as in Equation (3) [47].
6 of 14 The normalized burn ratio index (NBR), one of the most widely used image enhancements for mapping wildfires, is used to determine the burnt area and was improved by the machine learning classifier algorithm [48] and the basic formula, is expressed in Equation (4).
The method is used to detect burnt tracks is the burn ratio index technique and the difference ratio index technique. This index uses a temporal indication, which is the difference in the conditions before and after a fire for each polarization. The recording after a fire is selected through the October recording period, while the pre-fire condition data is the recording for the period from March to May. The backscatter value generated from each condition (before and after fire) is the average of several recordings, because it aims to minimize the backscatter anomaly or bias.
The radar burn ratio (RBR) in Equation (5) is modified from the NBR and is defined as the average backscatter ratio before and after a fire occurs at a certain polarization [49], meanwhile, the radar burn difference (RBD) in Equation (6) is the difference between the average scattering conditions before and after the fire [50].

RBR xy =
Postfire average backscatter xy Prefire average backscatter xy RBD xy = Postfire average backscatter xy − Prefire average backscatter xy The results of the two indexes with the two polarizations, produced four variables, namely the VV polarization RBR, VH polarization RBR, VV variation RBD, and the VH polarized RBD. The four layers are used as variables to determine the pixels of the burnt land object. The supervised classification is used to obtain the pixels on the burnt land objects by referring to a sample of the selected area. To obtain the best accuracy, three types of classification algorithms were selected, namely the support vector machine, random forest, and CART. The three classification results from these different algorithms are then compared with the field data on the burnt land, as a reference image.
The sample locations of this study are the peatlands in the Kalimantan Island region. The peatlands in the area often experience forest and land fires, related to the characteristics of the deep peat depth, which makes the organic content in the area very high, combined with unsustainable agricultural and plantation land conversion practices, and a dropping groundwater table, especially during the dry season [48].
The calculation of the accuracy test is carried out by taking a reference image on the burnt area, since it is in vector polygon format, the classification model was first converted into vectors to obtain a high calculation accuracy. The accuracy test uses desktop-based calculations for convenience in terms of a summarized database. The sample area was selected, based on random sampling, specifically for the areas with large burnt areas. Figure 3 is a co-polarized (VV) and cross-polarized (VH) SAR image from Sentinel-1a, along with the difference values and the Hotspot distribution data. The range of the ratio values for the RBR variable, ranges from 1 to 1.5 dB, while the RBD variable has a pixel value range of −2.5 to 1 dB. The hotspot distribution pattern on the RBR variable has the same pattern, which is in bright pixels or high ratio values in both the RBR VV and RBR VH variables. Meanwhile, in the RBD variable, the hotspot distribution pattern can be identified, but in reverse where the hotspot points are in dark pixels or low difference values. From the visualization above, it can be seen that the distribution of the hotspots has the same distribution pattern for the four variables where changes in land cover, due to land fires, cause changes in the backscatter value during the time series observations (before-after). the same pattern, which is in bright pixels or high ratio values in both the RBR VV and RBR VH variables. Meanwhile, in the RBD variable, the hotspot distribution pattern can be identified, but in reverse where the hotspot points are in dark pixels or low difference values. From the visualization above, it can be seen that the distribution of the hotspots has the same distribution pattern for the four variables where changes in land cover, due to land fires, cause changes in the backscatter value during the time series observations (before-after). The results of the experiment, is shown in Figure 4, in which, the radar burn difference shows significant burnt areas in South East and South West Kalimantan Island, as a dark area. While the radar burn ratio shows a difference between co-polarized (VV) and cross-polarized (VH), where there are more significant bright areas of the RBR VV that correlate with burnt area of the RBD. The RBR bright spots and the dark spots of the RBD correlates to the red fire location from the hotspot data The results of the experiment, is shown in Figure 4, in which, the radar burn difference shows significant burnt areas in South East and South West Kalimantan Island, as a dark area. While the radar burn ratio shows a difference between co-polarized (VV) and crosspolarized (VH), where there are more significant bright areas of the RBR VV that correlate with burnt area of the RBD. The RBR bright spots and the dark spots of the RBD correlates to the red fire location from the hotspot data.

Radar Burn Ratio and the Difference for the Burnt Area Detection
i. 2022, 12, x FOR PEER REVIEW 7 of 15 the same pattern, which is in bright pixels or high ratio values in both the RBR VV and RBR VH variables. Meanwhile, in the RBD variable, the hotspot distribution pattern can be identified, but in reverse where the hotspot points are in dark pixels or low difference values. From the visualization above, it can be seen that the distribution of the hotspots has the same distribution pattern for the four variables where changes in land cover, due to land fires, cause changes in the backscatter value during the time series observations (before-after). The results of the experiment, is shown in Figure 4, in which, the radar burn difference shows significant burnt areas in South East and South West Kalimantan Island, as a dark area. While the radar burn ratio shows a difference between co-polarized (VV) and cross-polarized (VH), where there are more significant bright areas of the RBR VV that correlate with burnt area of the RBD. The RBR bright spots and the dark spots of the RBD correlates to the red fire location from the hotspot data The resulting images for RBR and RBD for oth polarization can be seen on Figure 5 with corresponding RBR bright spots and RBD dark spots. The average RBD value for both the VH and VV polarization in the burnt area is very low, ranging from −1.5 to −2 dB, in contrast to the unburnt area, which has an average value of 0.1 to 0.3 dB. These values indicated that the burnt area can be detected while the RBR value of both the burnt and unburnt areas did not show significantly different values with 1.1 dB for the cross-polarized RBR and 1.3 dB for the co-polarized RBR, the unburnt area is at 1.0 dB as the baseline value for both components. From the values represented in Figure 6, we selected three components that have the largest margin between the burnt and unburnt areas to create a composite image, namely: RBD VH and VV and RBR VV. The composite image is then quantified with CART, random forest, and SVM, and compared with the Hotspot data.  The resulting images for RBR and RBD for oth polarization can be seen on Figure 5 with corresponding RBR bright spots and RBD dark spots. The average RBD value for both the VH and VV polarization in the burnt area is very low, ranging from −1.5 to −2 dB, in contrast to the unburnt area, which has an average value of 0.1 to 0.3 dB. These values indicated that the burnt area can be detected while the RBR value of both the burnt and unburnt areas did not show significantly different values with 1.1 dB for the crosspolarized RBR and 1.3 dB for the co-polarized RBR, the unburnt area is at 1.0 dB as the baseline value for both components. From the values represented in Figure 6, we selected three components that have the largest margin between the burnt and unburnt areas to create a composite image, namely: RBD VH and VV and RBR VV. The composite image is then quantified with CART, random forest, and SVM, and compared with the Hotspot data.

Burn Scar Mapping Machine-Learning Algorithms
The process of determining the burnt area is carried out using a machine learning algorithm on the composite data from four variables. The burnt area classification process in this research uses three types of classifiers, namely: SVM, random forest, and CART

Burn Scar Mapping Machine-Learning Algorithms
The process of determining the burnt area is carried out using a machine learning algorithm on the composite data from four variables. The burnt area classification process in this research uses three types of classifiers, namely: SVM, random forest, and CART with the parameters defined in Table 1. For the SVM classifier, the radial basis function (RBF) or the Gaussian kernel, is the most widely used kernel because of its high accuracy value. Usually used for datasets that are not linearly separated [51].  Figure 7 shows the differences in the classification results of the three classification algorithms and the distribution of the hotspots. The results of the classification using the CART classifier still show a lot of noise that appears as a burn area. While the SVM and random forest classifiers show the maximum noise reduction results and the burn area distribution pattern matches the hotspot distribution. Based on the significant burn area detected, we selected three locations as the area of interest, as indicated in Figure 7, to measure the model accuracy. Comparative data used to measure the accuracy of the model was the field data on the burnt land in the same time period. This data was obtained from integration of the Hotspot data, the fire extinguishing Based on the significant burn area detected, we selected three locations as the area of interest, as indicated in Figure 7, to measure the model accuracy. Comparative data used to measure the accuracy of the model was the field data on the burnt land in the same time period. This data was obtained from integration of the Hotspot data, the fire extinguishing points, and the manual delineation of the Landsat-8 and Sentinel-2 optical images.

Discussion
The composite formed in Figure 8 is a combination of RBR_VH, RBR_VV, RBD_VH and RBD_VV. The results of the experiment are shown in Figure 7, by taking samples in the burnt and unburnt areas the showed that the RBD variable was better at distinguishing between burnt and non-burnt objects, than the RBR variable.  Figure 9 shows the comparison of the Sentinel-2A satellite imagery, three classification results, and reference image. It can be seen that the model results from the three classifiers show the same pattern as the appearance of the burnt land, depicted in the Sentinel-2A image, as well as the reference image. The model generated by the random forest and the SVM algorithms has a high degree of similarity, compared to the model generated by the CART algorithm, where a lot of noise was found. In the CART classification, there are many pixel points that should not be part of the burnt area, so this noise causes the total burnt area to be larger than the random forest and SVM models.  Figure 9 shows the comparison of the Sentinel-2A satellite imagery, three classification results, and reference image. It can be seen that the model results from the three classifiers show the same pattern as the appearance of the burnt land, depicted in the Sentinel-2A image, as well as the reference image. The model generated by the random forest and the SVM algorithms has a high degree of similarity, compared to the model generated by the CART algorithm, where a lot of noise was found. In the CART classification, there are many pixel points that should not be part of the burnt area, so this noise causes the total burnt area to be larger than the random forest and SVM models. The results of the accuracy test of the three classifiers used in this research are shown in Table 2. Where the overall accuracy to detect the burnt and unburnt areas are above 75%. The model generated by the SVM classifier have the highest accuracy of 88.26%, compared to the field data, while CART classification has a much lower accuracy at 78.40%. The examination on the confusion matrix shows that CART classifies the unburnt area as burnt with a significantly wider area, this means that the classifier is more prone to mis-classification on the edges of the burnt areas. Based on these results, the SVM yields the best accuracy for classifying the forest area as burnt. As our research aims for a quick mapping method, the second parameter to test is the processing time for each method. The burnt area mapping, using the SVM and RF classification methods each require the processing time of 1 h, while CART takes 2 h. This processing time represents the temporal resolution of the map produced or, in other The results of the accuracy test of the three classifiers used in this research are shown in Table 2. Where the overall accuracy to detect the burnt and unburnt areas are above 75%. The model generated by the SVM classifier have the highest accuracy of 88.26%, compared to the field data, while CART classification has a much lower accuracy at 78.40%. The examination on the confusion matrix shows that CART classifies the unburnt area as burnt with a significantly wider area, this means that the classifier is more prone to mis-classification on the edges of the burnt areas. Based on these results, the SVM yields the best accuracy for classifying the forest area as burnt. As our research aims for a quick mapping method, the second parameter to test is the processing time for each method. The burnt area mapping, using the SVM and RF classification methods each require the processing time of 1 h, while CART takes 2 h. This processing time represents the temporal resolution of the map produced or, in other words, the period of updating the map. When compared to the minimum updating period of the forest fire mapping process that is currently valid and which is 12 h [33], the quick mapping method can speed up the process so that this information can more quickly reach the fire mitigation team in the field.

Conclusions
The characteristics of SAR data that are able to identify objects, based on their physical form, are an advantage in identifying burn scar areas. The characteristics of burnt land in the form of open land will have a low backscatter value, compared to unburnt forest areas, which have a coarser texture, so that they have a higher backscatter. In addition, another advantage, when compared to the optical imagery is that it is able to provide visual objects without being constrained by clouds or weather and smoke that usually covers objects, especially when forest fires occur.
The identification of the burnt areas using SAR backscatter differences data, in the conditions before and after a fire, is able to produce s good accuracy. This can be seen from the results of the accuracy test of the three classifiers used in this research, showing an overall accuracy of over 88% with a processing time reduced from 12 h to 1 h. Based on the results, we select the support vector machine method as the best to be applied to support the forest fire mitigation system. Author Contributions: R. is the principal researcher, R.A. provided with field data and machine learning. I.R. also contributed to the machine learning and paper writing. A., S.I., and L.M. contributed to the GIS mapping and analysis of the data. All authors have read and agreed to the published version of the manuscript.