Forest Fire Monitoring and Positioning Improvement at Subpixel Level: Application to Himawari-8 Fire Products

: Forest ﬁres are among the biggest threats to forest ecosystems and forest resources, and can lead to ecological disasters and social crises. Therefore, it is imperative to detect and extinguish forest ﬁres in time to reduce their negative impacts. Satellite remote sensing, especially meteorological satellites, has been a useful tool for forest-ﬁre detection and monitoring because of its high temporal resolution over large areas. Researchers monitor forest ﬁres directly at pixel level, which usually presents a mixture of forest and ﬁre, but the low spatial resolution of such mixed pixels cannot accurately locate the exact position of the ﬁre, and the optimal time window for ﬁre suppression can thus be missed. In order to improve the positioning accuracy of the origin of forest ﬁre (OriFF), we proposed a mixed-pixel unmixing integrated with pixel-swapping algorithm (MPU-PSA) model to monitor the OriFFs in time. We then applied the model to the Japanese Himawari-8 Geostationary Meteorological Satellite data to obtain forest-ﬁre products at subpixel level. In this study, the ground truth data were provided by the Department of Emergency Management of Hunan Province, China. To validate the positioning accuracy of MPU-PSA for OriFFs, we applied the model to the Himawari-8 satellite data and then compared the derived ﬁre results with ﬁfteen reference forest-ﬁre events that occurred in Hunan Province, China. The results show that the extracted forest-ﬁre locations using the proposed method, referred to as forest ﬁre locations at subpixel (FFLS) level, were far closer to the actual OriFFs than those from the modiﬁed Himawari-8 Wild Fire Product (M-HWFP). This improvement will help to reduce false ﬁre claims in the Himawari-8 Wild Fire Product (HWFP). We conducted a comparative study of M-HWFP and FFLS products using three accuracy-evaluation indexes, i.e., Euclidean distance, RMSE, and MAE. The mean distances between M-HWFP ﬁre locations and OriFFs and between FFLS ﬁre locations and OriFFs were 3362.21 m and 1294.00 m, respectively. The mean RMSEs of the M-HWFP and FFLS products are 1225.52 m and 474.93 m, respectively. The mean MAEs of the M-HWFP and FFLS products are 992.12 m and 387.13 m, respectively. We concluded that the newly proposed MPU-PSA method can extract forest-ﬁre locations at subpixel level, providing higher positioning accuracy of forest ﬁres for their suppression.


Introduction
Forest fires, when not suppressed in time, can spread over large areas and sometimes may last for years or decades. There were a total of 116,171 active forest-fire events in China from 2001 to 2015 [1]. Forest-associated biomass burning could be, in part, responsible for the increase in greenhouse gases (GHG) in the atmosphere [2,3]. Forest fires are among the biggest threats to plant ecological systems and resources [1,4], and can cause ecological disasters and social crises. For instance, in Liangshan Yi Autonomous Prefecture, Sichuan Province, China, serious forest fires occurred three times during 2019-2021 causing 49 casualties [5-7], and the 2019-2020 Australia megafires lasted for five months from September 2019 to February 2020, killing or displacing nearly 3 billion animals [8].
A series of precautions should be taken to mitigate adverse fire impacts and reduce the likelihood of forest fires [9]. Effective forest-fire detection and monitoring has always been a strong focus of research [10][11][12][13][14]. Meteorological satellites can acquire data continuously over large areas and in high temporal resolution. Thus, they have been widely used for forest-fire detection and monitoring [13,[15][16][17][18]. Over the past decades, based on the data from geostationary orbit (GEO) satellites or polar-orbiting satellites, researchers have proposed many forest-fire detection and monitoring algorithms. These can be divided into five categories [19]: (1) bi-spectral methods, such as that of Dozier [20], measuring the surface radiant temperature of subpixel spatial resolution; (2) modified bi-spectral methods, such as multiple endmember spectral mixture analysis (MESMA) [21], bi-spectral infrared detection (BIRD) [22], and the Giglio and Schroeder [23] method; (3) threshold methods, such as multi-channel threshold algorithms [24]; (4) spatial contextual methods, such as MODIS Collection 6 active fire algorithms [25]; and (5) multi-temporal fire detection methods, such as robust satellite techniques for FIRES detection and monitoring (RST-FIRES) algorithm [14]. However, due to the low spatial resolution of the satellites, researchers conduct forest-fire monitoring research directly at pixel scale, which can lead to significant fire-positioning error when the pixel shows a mixture of forest and fire [16,[26][27][28][29][30]. Atkinson [31,32] first proposed the subpixel concept to find specific spatial distribution information in a mixed pixel. Forest-fire spatial location must be specified in a mixed pixel at the subpixel level.
For satellite data of low spatial resolution, mixed pixels are a common phenomenon [33][34][35]. They can be an obstacle in remote sensing applications in forest-fire detection, causing the estimated location of a forest fire to vary widely from the actual position, depending on spatial resolution. To overcome this obstacle, researchers proposed the traditional exponential model, which is often found to produce inconsistent results, since a mixed pixel includes both objective and non-objective materials [36]. A series of unmixing methods have been proposed, such as maximum-margin criterion and derivative weights (MD-WSU) [37], minimum-volume transform (MVT) [38], vertex component analysis (VCA) [39], simplex growing algorithm (SGA) [40], minimum-volume-constrained nonnegative matrix factorization (MVC-NMF) [41], successive projection algorithm (SPA) [42], etc. He et al. [36] showed that spectral unmixing methods have advantages over the traditional exponential model due to their ability to decompose the mixed pixel into several fractional abundances. Shao and Lan [37] used the MDWSU technique to obtain more accurate endmembers and abundance estimates. Nascimento and Dias [39] performed a series of experiments using simulated and real data, and found that the VCA algorithm performs better than the pixel purity index (PPI) method, and better than or similar to the N-FINDR algorithms [43]. Miao and Qi [41] proposed the MVC-NMF method to extract unsupervised endmembers from highly mixed image data. Zhang et al. [42] proposed SPA, which can provide a general guideline to constrain the total number of endmembers.
However, spectral unmixing methods only solve the problem of endmember types and abundances, without determining endmember location in a mixed pixel. To solve this problem, Atkinson [31,32] proposed the subpixel concept for identifying precise spatial locations of endmembers in a mixed pixel based on unmixing analysis results [44]. On the basis of Atkinson's study, many algorithms have been developed [45], such as spatial attraction model (SAM), double-calculated spatial attraction model (DSAM) [34], subpixel learning algorithms [46], subpixel edge-detection method [47], and so on. These subpixel techniques have found many applications [48][49][50][51]. For instance, the random forests and Spatial Attraction Model (RFSAM) was applied to remote sensing images to improve the accuracy of subpixel mapping of wetland flooding [52]. Li et al. [53] proposed the spatiotemporal subpixel land cover mapping (STSPM) method, demonstrating that it can predict land cover maps accurately. Ling et al. [54] proposed a new approach aimed at efficiently and accurately monitoring reservoir surface water area variations using daily Remote Sens. 2022, 14, 2460 3 of 19 moderate resolution imaging spectroradiometer (MODIS) images to explore subpixel scale information. Deng and Zhu [55] proposed the continuous subpixel monitoring (CSM) method and successfully applied it to mapping urban impervious surface area (ISA %) at subpixel scale and characterizing its dynamics in Broome County, New York.
However, there are very few applications of these subpixel algorithms to forest-fire location detection and fire monitoring [21]. In order to rapidly detect and monitor forestfire locations to reduce damage and save lives and properties, ascertaining the exact locations of the origins of fires is important in order to gain quick access for firefighters. Subpixel algorithms can help locate fires with subpixel accuracy. However, current forestfire monitoring methods using remote sensing data can locate forest fires only at pixel-size level. Location accuracy can be improved with subpixel-level detection.
The objective of this study is to develop a new algorithm to improve the detection of fire locations at subpixel level. To this end, we proposed a mixed pixel unmixing integrated with pixel-swapping algorithm (MPU-PSA) model, in which spectral unmixing analysis is performed using the Newton's Method, then the unmixing analysis results are processed with the pixel-swapping algorithm to obtain OriFFs of forest fires. We then applied the new algorithm to the Himawari-8 L1 gridded data (HLGD) and Himawari-8 Wild Fire Product (HWFP), taking advantage of the high temporal resolution of the Himawari-8 satellite data to achieve more accurate forest-fire locations. The forest-fire information provided by the Department of Emergency Management of Hunan Province, China, was used as ground truth reference data. The result showed that the algorithm improved the location accuracy of the origin of forest fire (OriFF), saving time for forest-fire extinguishment.

Materials and Methods
The present research used the following data: (a) satellite data, including Himawari-8, Sentinel 2, and Landsat 8 satellite data (Section 2.2); (b) ground truth data (Section 2.4). First, we ordered the Himawari-8 (https://www.eorc.jaxa.jp/ptree/, accessed on 10 December 2021) HLGD and HWFP data for the dates closest to the starting time of a forest fire. Then, band 7 and band 14 of the HLGD data were selected and fused to produce dataset B7_14 in Matlab R2016a (v 9.0, Mathworks, Natick, MA, USA). Wildfires detected in the HWFP data that did not match the ground truth forest-fire information were eliminated to generate the modified Himawari-8 Wild Fire Product (M-HWFP) dataset, but this remained at pixel level. Finally, the datasets B7_14 and M-HWFP were used as inputs to our proposed model to obtain the fire dataset forest-fire locations at subpixel (FFLS) level. Figure 1 shows the flowchart of data processing and analysis, which are explained in detail over the following sub-sections.

Study Area
Hunan Province is located in the middle Yangtze River, central China, within an area enclosed by 108 • 47 ~114 • 15 longitude and 24 • 38 ~30 • 08 latitude (Figure 2). Over a year, the temperature in Hunan Province typically varies from 7.78 • C to 33.33 • C. The rainiest month is June, with an average rainfall of 177.8 mm. The driest month is December, with an average rainfall of 35.56 mm. The relatively warm and humid weather results in lush vegetation and an evergreen landscape. Hunan province has a subtropical evergreen broad-leaved forest area with a forest coverage rate of 59.96%. Forest and wetland ecosystems are two major natural ecosystems. The ecological service function of the forest ecosystems in Hunan Province is worth 1.01 trillion Chinese Yuan. There are about 200 significant ecological zones around the world. Hunan province has two forest and wetland ecosystems (subtropical evergreen broad-leaved forest ecological zones in the Wuling Xuefeng Mountains and the Nanling Luoxiao Mountains), and they are known as the world's most valuable ecological regions within their latitude zone. Due to the high forest coverage rate in Hunan Province, China, forest fires are significant potential hazards. Therefore, monitoring and rapidly suppressing forest fires is imperative for maintaining ecological integrity.

Study Area
Hunan Province is located in the middle Yangtze River, central China, within an are enclosed by 108°47′∼114°15′ longitude and 24°38′∼30°08′ latitude (Figure 2). Over a year the temperature in Hunan Province typically varies from 7.78 °C to 33.33 °C. The rainies month is June, with an average rainfall of 177.8 mm. The driest month is December, wit an average rainfall of 35.56 mm. The relatively warm and humid weather results in lus vegetation and an evergreen landscape. Hunan province has a subtropical evergree broad-leaved forest area with a forest coverage rate of 59.96%. Forest and wetland ecosys tems are two major natural ecosystems. The ecological service function of the forest eco systems in Hunan Province is worth 1.01 trillion Chinese Yuan. There are about 200 sig Xuefeng Mountains and the Nanling Luoxiao Mountains), and they are known as the world's most valuable ecological regions within their latitude zone. Due to the high forest coverage rate in Hunan Province, China, forest fires are significant potential hazards. Therefore, monitoring and rapidly suppressing forest fires is imperative for maintaining ecological integrity.

Satellite Data Processing
Himawari-8, Sentinel 2, and Landsat 8 satellite data were reprojected into a common Universal Transverse Mercator (UTM) projection with WGS84 as the datum so that all datasets could be overlaid within the same coordinate system. After each satellite dataset was clipped, the cloud cover was less than 10%. The Himawari-8 satellite data were used for detecting and monitoring the OriFF attribute at high temporal resolution (10 min). The Sentinel 2 and Landsat 8 images with higher spatial resolution were used as base maps to verify whether the hotspots detected in the Himawari-8 data were within the forest.
In this study, the HLGD data in Network Common Data Form (NetCDF) format and the HWFP data in comma-separated values (CSV) format were downloaded via File Transfer Protocol (FTP) in FileZilla Client 3.3.2 software. The Advanced Himawari Imager (AHI) carried by the Himawari-8 satellite scans the Pacific Ocean hemispheric region every 10 min, resulting in the production of 142 images over a specific scene per day [30,56]. The high temporal resolution of the Himawari-8 data makes it possible to semicontinuously monitor and detect forest fires, even though its spatial resolution is low. The JAXA Himawari Monitor User's Guide is available on website https://www.eorc.jaxa.jp/ptree/userguide.html (accessed on 12 December 2021). Table 1 shows the characteristics of band 7 (3.85 μm) and band 14 (11.20 μm) selected for this study. These two bands are the most important bands used for forest-fire monitoring and detection [13,16,30]. This study acquired 2130 HLGD images and 2130 HWFP images for the fifteen forest-fire events used as reference fires. To reduce the influence of clouds, we filtered the clouds for each pixel when calculating the average background brightness

Satellite Data Processing
Himawari-8, Sentinel 2, and Landsat 8 satellite data were reprojected into a common Universal Transverse Mercator (UTM) projection with WGS84 as the datum so that all datasets could be overlaid within the same coordinate system. After each satellite dataset was clipped, the cloud cover was less than 10%. The Himawari-8 satellite data were used for detecting and monitoring the OriFF attribute at high temporal resolution (10 min). The Sentinel 2 and Landsat 8 images with higher spatial resolution were used as base maps to verify whether the hotspots detected in the Himawari-8 data were within the forest.
In this study, the HLGD data in Network Common Data Form (NetCDF) format and the HWFP data in comma-separated values (CSV) format were downloaded via File Transfer Protocol (FTP) in FileZilla Client 3.3.2 software. The Advanced Himawari Imager (AHI) carried by the Himawari-8 satellite scans the Pacific Ocean hemispheric region every 10 min, resulting in the production of 142 images over a specific scene per day [30,56]. The high temporal resolution of the Himawari-8 data makes it possible to semi-continuously monitor and detect forest fires, even though its spatial resolution is low. The JAXA Himawari Monitor User's Guide is available on website https://www.eorc.jaxa. jp/ptree/userguide.html (accessed on 12 December 2021). Table 1 shows the characteristics of band 7 (3.85 µm) and band 14 (11.20 µm) selected for this study. These two bands are the most important bands used for forest-fire monitoring and detection [13,16,30]. This study acquired 2130 HLGD images and 2130 HWFP images for the fifteen forest-fire events used as reference fires. To reduce the influence of clouds, we filtered the clouds for each pixel when calculating the average background brightness temperature. Furthermore, to improve the running speed, we used the subset data from shapefile batch of ENVI 5.3 to clip HLGD. 0.50~0. 68 15 The Sentinel 2 satellite data were downloaded from https://earthexplorer.usgs.gov/ (accessed on 15 December 2021) as Level 1C product. This product was orthorectified to topof-atmosphere (TOA) reflectance. The Sentinel 2 satellite data were packaged into images according to the spatial resolution of 10 m (band 2, 3, 4, and 8), 20 m (bands 5, 6, 7, 8b, 11, and 12), and 60 m (bands 1, 9, and 10), respectively. The Sentinel 2 TOA reflectance data were processed to bottom-of-atmosphere (BOA) reflectance after radiometric calibration and atmospheric correction using Sen2cor (V2.8). Table 1 shows the characteristics of band 2 (Blue), band 3 (Green), band 4 (Red), and band 8 (NIR) selected for this study.

Input Data Preparation
Using Visual Studio 2012 software based on C# programming (https://github.com/ HZXu1/MPU_PSA.git, accessed on 10 May 2022), we batch-processed the 142 HWFP daily images and extracted and saved hotspot information in a CSV file. Furthermore, all hotspot information of the 142 HWFP images, especially the earliest UTC time of every forest-fire event, was obtained using ArcGIS 10.4 software through the tabulated collected fire information (TCFI). For example, taking the first wildfire event (i.e., 1HeCo, which is the serial number of forest-fire event) (Figure 3), the red dot inside the yellow circle in the HWFP image gives the location of a hotspot in the IR images detected 64 times that day, as indicated in the third red circle (Figure 3b). False fires in the HWFP images were eliminated by comparing them to the TCFI data and burned area base map, and a Modified Himawari-8 Wild Fire Product (M-HWFP) was obtained. Combining the base map and TCFI data, we can see that the earliest monitoring time of the forest-fire event was 11:32 am (UTC time). Then, the M-HWFP data and the preprocessed HLGD data obtained at the timepoint closest to a fire's starting time in a target area were used as input to the Visual Studio 2012 software based on C# programming.

Ground Truth Data
Tabulated collected fire information (TCFI) as provided by the Department of Emergency Management of Hunan Province was used as ground truth data. Information about each fire event includes location, weather, starting and extinguishing times, damaged area, and cause. The Department of Emergency Management of Hunan Province confirms and reports forest-fire events when a forest fire occurs (Figure 4). Field measurements, including those pertaining to OriFF and damaged areas from a forest fire, were performed by well-trained personnel with global navigation satellite system survey equipment and unmanned aerial vehicles (UAV. The type of platform was Phantom 4 Pro V2.0, and the sensor was 1-inch CMOS with 20 million effective pixels. After the route was designed, the data were automatically collected and stored in the fuselage memory card). From the TCFI data between 2018 and 2021, the fifteen forest-fire events that resulted in the most environmental damage were selected as reference data.
14, x FOR PEER REVIEW 7 of 20

Ground Truth Data
Tabulated collected fire information (TCFI) as provided by the Department of Emergency Management of Hunan Province was used as ground truth data. Information about each fire event includes location, weather, starting and extinguishing times, damaged area, and cause. The Department of Emergency Management of Hunan Province confirms and reports forest-fire events when a forest fire occurs (Figure 4). Field measurements, including those pertaining to OriFF and damaged areas from a forest fire, were performed by well-trained personnel with global navigation satellite system survey equipment and unmanned aerial vehicles (UAV. The type of platform was Phantom 4 Pro V2.0, and the sensor was 1-inch CMOS with 20 million effective pixels. After the route was designed, the data were automatically collected and stored in the fuselage memory card). From the TCFI data between 2018 and 2021, the fifteen forest-fire events that resulted in the most environmental damage were selected as reference data.

Ground Truth Data
Tabulated collected fire information (TCFI) as provided by the Department of Emergency Management of Hunan Province was used as ground truth data. Information about each fire event includes location, weather, starting and extinguishing times, damaged area, and cause. The Department of Emergency Management of Hunan Province confirms and reports forest-fire events when a forest fire occurs (Figure 4). Field measurements, including those pertaining to OriFF and damaged areas from a forest fire, were performed by well-trained personnel with global navigation satellite system survey equipment and unmanned aerial vehicles (UAV. The type of platform was Phantom 4 Pro V2.0, and the sensor was 1-inch CMOS with 20 million effective pixels. After the route was designed, the data were automatically collected and stored in the fuselage memory card). From the TCFI data between 2018 and 2021, the fifteen forest-fire events that resulted in the most environmental damage were selected as reference data.  information using pixel-Swapping algorithm (PSA). We used the MPU-PSA method to improve determination of the spatial locations of forest fires at a subpixel scale. The M-HWFP and B7_14 datasets were processed using the MPU-PSA method to obtain the final result. Based on the Himawari-8 satellite data, each pixel was evenly divided into 5 × 5 subpixels. The spatial resolution of each subpixel was thus 400 m. The final result was forest-fire spatial distribution information in a mixed pixel.

Mixed Pixel-Unmixing (MPU) Analysis
The Planck function for emittance at wavelength λ of a blackbody at kinetic temperature T is shown in Equation (1): The method for satellite identification of surface temperature fields at subpixel resolution was proposed by Dozier [20]. Assuming a mixed pixel is composed of active fire with high temperature in proportion P (where 0 ≤ P ≤ 1) and a background with normal temperature in the remaining proportion (1 − P), radiance of a mixed pixel can be expressed by the following Equation (2): where subscript i indicates the i-th band; L i is the radiance at i-th band (W·m −2 ·µm −1 ·sr −1 ); L i, fire the target forest-fire radiance (W·m −2 ·µm −1 ·sr −1 ); L i, bg the background (bg) radiance (W·m −2 ·µm −1 ·sr −1 ); C 1 and C 2 are the first (3.74 × 10 −16 W·m 2 ) and second (1.44 × 10 −2 m·K) Planck constants, respectively; λ i is the wavelength (µm) of i-th band; and T i, fire and T i,bg are the temperature (K) of the target forest fire and the background, respectively. We chose band 7 and band 14 in this study, and the radiances of a mixed pixel for band 7 and band 14 were given as below according to Equation (2).

Pixel-Swapping Algorithm (PSA)
Atkinson [44] proposed the PSA method, which was designed to process an image of land cover proportions in K = 2 classes. In this study, PSA needed to be carried out on the basis of forest-fire abundance information. Thus, the PSA method was applied based on the unmixing results of mixed pixels (Section 2.5.1). If 10 × 10 (=100) subpixels are to Remote Sens. 2022, 14, 2460 9 of 19 be mapped within each mixed pixel, a land cover class (with or without forest fire) with a proportion of 57 percent would mean that 57 subpixels were allocated to that class.
PSA is comprised of three basic steps. Firstly, for every subpixel, the attractiveness A(p c i ) of a subpixel p i in c-th class is predicted as a distance-weighted function of its j = 1, 2, . . . , J neighbors: where z(x c j ) is the value of j-th subpixel belonging to c-th class, and λ ij is a distance dependent weight predicted by Equation (6): where h ij is the distance between two subpixels, p i and p j , and a is the non-linear parameter of the exponential model. Secondly, once the attractiveness of each subpixel has been predicted based on the current arrangement of subpixel classes, the subpixel algorithm ranks the values on a pixel-by-pixel basis. For each pixel, the least attractive subpixel currently allocated to a "1" (i.e., a "1" surrounded mainly by "0"s) is stored (shown in Equation (7)): The most attractive subpixel currently allocated to a "0" (i.e., a "0" surrounded mainly by "1"s) is also stored (shown in Equation (8)): Lastly, classes of subpixels are swapped as follows: if the attractiveness of the least attractive subpixel is less than that of the most attractive subpixel, the classes are swapped for the subpixels in question (shown in Equation (9)): If it is more attractive, no change is made.

Accuracy Assessment
In forestry firefighting, firefighters must arrive at the OriFF location as soon as possible to extinguish the fire. Therefore, it is important to know the distance between the estimated OriFF location and the actual OriFF location. It is necessary to compare the FFLS and M-HWFP datasets to decide which derived OriFF is closer to the actual OriFF. To this end, we added Euclidean distance as one of the accuracy-evaluation indexes. Moreover, using distance was able show more intuitively whether the M-HWFP data or the FFLS data had better OriFF estimation; and using distance could also reflect the positioning accuracy more intuitively. To evaluate the performance of the MPU-PSA for OriFF detection, fifteen forest-fire events with TCFI data were selected as references, in order to assess: (1) the accuracy of distance comparison analysis (distance from the M-HWFP-estimated OriFF or from the FFLS-estimated OriFF to the actual OriFF), and (2) accuracy-comparison analysis (forest-fire positioning accuracy comparison between M-HWFP and FFLS). Three widely used indexes were calculated for the evaluation: Euclidean distance, root mean square error (RMSE), and mean absolute error (MAE). The two metrics, RMSE and MAE, are measures of the difference between the values predicted by the model and those of the actual phenomenon [57]. The lower the RMSE and MAE, the better the precision and accuracy of the model [58]. The RMSE is significantly sensitive to large values and outliers, and the MAE is suitable to describe uniformly distributed errors [59,60]. Therefore, it is necessary to combine RMSE with MAE to evaluate the variation in model errors [61,62]. They can be formulated as follows: To illustrate the better positioning accuracy of d FFLS-OriFF , we introduced positioningprecision rate as given below: where d M-HWFP-OriFF is the Euclidean distance between the M-HWFP-estimated OriFF and the actual OriFF; d FFLS-OriFF is the Euclidean distance between the FFLS-estimated OriFF and the actual OriFF; (x i , y i ) and (x j , y j ) are coordinates of samples in two-dimensional space. The second and third indices are RMSE and MAE, as given below: where ∧ d i represents the estimated fire locations in FFLS of fire event i that were obtained using the MPU-PSA method; d i is the actual OriFF of fire event i from the TCFI data; and m (=15) is the number of all forest-fire events evaluated; RMSE and MAE were calculated for both M-HWFP and FFLS relative to the actual OriFFs.

Forest-Fire Detection
For convenience, we used green triangles to represent the actual OriFFs, or the true forest-fire locations. We used purple dots to represent the M-HWFP-estimated OriFFs, and red rectangles to represent the FFLS-estimated OriFFs, which have subpixels of 400 m by 400 m. As shown in Figure 5, among the fifteen forest-fire events, red rectangles (i.e., locations of FFLS-estimated OriFFs) and green triangles (i.e., locations of actual OriFFs) of fourteen forest-fire events were surrounded by purple dots (i.e., locations of M-HWFPestimated OriFFs). It can be clearly seen that the red rectangles are closer than the purple dots to the green triangles. Eight green triangles are closest to or even overlap with the subpixel results. The red line represents the county boundary. We can see that among the fourteen forest-fire events, nine occurred at the junction of counties. Combining this knowledge with geographic information system (GIS) analysis, we can infer that these places have denser forests or more forest fuel, factors which need further verification. In Figure 6, green triangles again represent the actual OriFF locations, purple dots represent forest-fire locations detected by Himawari-8 satellite, and red rectangles represent the locations of the FFLS-estimated OriFFs, as in Figure 5. Sentinel 2 was the first choice for the base map, but if there were no suitable Sentinel 2 data for the forest-fire event area, we used Landsat 8 instead. From Figure 6, it can be seen that among the forest vegetation, the area burnt by the fire can be easily recognized by visual interpretation. Remote Sens. 2022, 14, x FOR PEER REVIEW 12 of 20   Therefore, we can conclude that fire locations determined at subpixel level in FFLS were closer to the actual OriFF locations, while fire locations determined at pixel level in M-HWFP were further away from the actual OriFF locations. These results demonstrate that the MPU-PSA technique has the advantage of detecting fire locations at subpixel level, thus performs more accurately than traditional forest-fire detection method at pixel level.

Comparison of M-HWFP and FFLS in Positioning Accuracy
Comparative analysis results from d M-HWFP-OriFF and d FFLS-OriFF indices data are shown in Table 2 to illustrate the better performance of FFLS in comparison with M-HWFP in locating forest fires. M-HWFP provides the pixel-level fire locations, which were retrieved from the infrared band (IR) images of the Himawari-8 satellite using the method developed by JAXA/EORC [63,64]. Using the proposed method, we attempted to obtain more precise spatial locations of OriFFs than those obtained using traditional forest-fire detection methods, such as M-HWFP. Knowing the exact locations of forest fires is imperative in effective forestry firefighting, especially in virgin and dense forests among high mountains with bumpy or even no roads. As shown in Table 2, the coordinates of fire locations are projected coordinates with easting and northing in meters. Taking the first forest-fire event (i.e., 1HeCo) as an example, TCFI data showed that the forest fire started at 12:00 UTC on 24 September 2019, and the closest timepoint to forest-fire starting time recorded by Himawari-8 satellite was 11:30 UTC on 24 September 2019. In Table 2, d M-HWFP-OriFF and d FFLS-OriFF are 6294.90 m and 1573.00 m, respectively. Therefore, for the first forest-fire event (i.e., 1HeCo), the FFLS-estimated OriFF was closer to the actual OriFF than the M-HWFP-estimated OriFF, with a positioning-precision rate of 75.01%.
As shown in Figure 7 and Table 2, the positioning-precision rates of the fifteen forestfire events all exceeded 50%. The mean d M-HWFP-OriFF , d FFLS-OriFF, and positioning-precision rate were 3362.21 m, 1294.00 m, and 60.99%, respectively. The largest positioning-precision rate was 79.58%, and the lowest positioning-precision rate was 51.28%. The positioningprecision rates of two forest-fire events exceeded 70.00%, and the positioning-precision rates of five forest-fire events exceeded 60.00%. The MPU-PSA method has an advantage in monitoring OriFFs over the traditional forest-fire detection method.

Comparison of M-HWFP and FFLS in Positioning Accuracy
Comparative analysis results from dM-HWFP-OriFF and dFFLS-OriFF indices data are shown in Table 2 to illustrate the better performance of FFLS in comparison with M-HWFP in locating forest fires. M-HWFP provides the pixel-level fire locations, which were retrieved from the infrared band (IR) images of the Himawari-8 satellite using the method developed by JAXA/EORC [63,64]. Using the proposed method, we attempted to obtain more precise spatial locations of OriFFs than those obtained using traditional forest-fire detection methods, such as M-HWFP. Knowing the exact locations of forest fires is imperative in effective forestry firefighting, especially in virgin and dense forests among high mountains with bumpy or even no roads. As shown in Table 2, the coordinates of fire locations are projected coordinates with easting and northing in meters. Taking the first forest-fire event (i.e., 1HeCo) as an example, TCFI data showed that the forest fire started at 12:00 UTC on 24 September 2019, and the closest timepoint to forest-fire starting time recorded by Himawari-8 satellite was 11:30 UTC on 24 September 2019. In Table 2, dM-HWFP-OriFF and dFFLS-OriFF are 6294.90 m and 1573.00 m, respectively. Therefore, for the first forest-fire event (i.e., 1HeCo), the FFLS-estimated OriFF was closer to the actual OriFF than the M-HWFPestimated OriFF, with a positioning-precision rate of 75.01%.
As shown in Figure 7 and Table 2, the positioning-precision rates of the fifteen forestfire events all exceeded 50%. The mean dM-HWFP-OriFF, dFFLS-OriFF, and positioning-precision rate were 3362.21 m, 1294.00 m, and 60.99%, respectively. The largest positioning-precision rate was 79.58%, and the lowest positioning-precision rate was 51.28%. The positioning-precision rates of two forest-fire events exceeded 70.00%, and the positioning-precision rates of five forest-fire events exceeded 60.00%. The MPU-PSA method has an advantage in monitoring OriFFs over the traditional forest-fire detection method.

Accuracy Comparison between M-HWFP and FFLS
We chose two commonly used accuracy-evaluation indexes, RMSE and MAE, to evaluate the accuracy of M-HWFP and FFLS. Compared with M-HWFP, FFLS showed more accurate spatial location of OriFFs.  Table 3, the largest decreases in RMSE and MAE were 80.35% and 79.58%, respectively; and the smallest decreases in RMSE and MAE were 50.50% and 51.28%, respectively. For the 15 reference forest-fire events, the decrease in RMSE and MAE exceeded 50.00%. There was one forest-fire event with an RMSE decrease exceeding 80.00%, and there were six forest-fire events with RMSE decrease exceeding 60.00%. There were two forest-fire events with an MAE decrease exceeding 70.00%, and seven forest-fire events with an MAE decrease exceeding 60.00%. Furthermore, the average RMSE and MAE of M-HWFP were 1225.52 m and 992.12 m, respectively, owing to the low spatial resolution of the Himawari-8 satellite. On the other hand, the average RMSE and MAE of FFLS were 474.93 m and 387.13 m, respectively. FFLS generally showed a more refined detection of spatial location of forest fires than M-HWFP.

Discussion
Traditionally, forest-fire detection has mainly relied on artificial ground patrols and watchtowers [65]. Slow acquisition of fire information, such as location and starting time, using the conventional methods often lead to delayed action and ultimately to forestfire disasters. The advent of remote-sensing techniques, such as satellite-based [66,67], aerial-based [68], and ground-surveillance-camera-based [69] techniques, offers researchers alternative or even better methods for forest-fire detection and monitoring. Meteorological satellites are widely used in forest-fire detection and monitoring because of their high temporal resolution. For example, Jang et al. [13] detected most forest fires successfully based on the Himawari-8 satellite data. In this study, we could monitor and detect forest fires every 10 min using the Himawari-8 satellite data. The high temporal resolution makes it possible to detect forest fires in almost real time.
However, forest-fire monitoring at pixel scale can cause large errors in fire locations. The ability of subpixels to specify the distribution location of mixed pixel composition could somewhat improve this situation [31,32]. Burnt-area estimation at subpixel level using Advanced Very-High-Resolution Radiometer (AVHRR) data has a mean error of 6.5% [70]. In this study, we proposed the MPU-PSA method to detect OriFFs. Compared with M-HWFP, OriFF detection by the new method (MPU-PSA) has higher positioning accuracy. On the one hand, the subpixel result (i.e., FFLS) maintains the high temporal resolution of the Himawari-8 satellite but provides more accurate positions of forest fires. Identification of fire locations at subpixel level enhanced forest-fire monitoring. It is of great importance that forest fires are extinguished in a timely manner at their OriFF locations. The MPU-PSA method was successfully applied to fifteen forest fires in Hunan province, China. However, in future research, it should be a priority to carry out research to validate and conduct the MPU-PSA in other regions. Furthermore, it is necessary to further improve forest-fire positioning accuracy at subpixel scale by improving the algorithm.

Conclusions
In this study, we proposed a forest-fire detection model for forest-fire detection at subpixel level (MPU-PSA), and applied the model to the Himawari-8 satellite data for fire detection. The proposed model involves two main procedures: (1) pixel-unmixing analysis to obtain forest-fire distribution, and (2) adaptation of the pixel-swapping algorithm to generate the FFLS dataset for OriFF estimation. The MPU-PSA method was applied to fifteen forest-fire events occurring in Hunan Province, China. Compared with the M-HWFP product, the FFLS dataset has better positioning accuracy, which was reflected in three accuracy-evaluation indices: Euclidean Distance, RMSE, and MAE. Positioning error was significantly improved, decreasing from 3362.21 m to 1294.00 m (with the average positioning-precision rate being 60.99%). The mean RMSE decreased from 1225.52 m to 474.93 m (decrease rate = 60.96%), and the mean MAE decreased from 992.12 m to 387.13 m (decrease rate = 60.99%). All these results show that (1) the MPU-PSA method developed in this study can estimate forest-fire locations at subpixel level, which is more accurate than the conventional pixel-level method, especially for geostationary satellite products, such as the Himawari-8 Wild Fire Product, which has low spatial resolution but high temporal resolution, and (2) the FFLS product with geographic information can provide more accurate OriFF locations than the M-HWFP product. This new forest-fire detection technique can assist in forest-fire extinguishment by quickly and accurately detecting OriFF locations.