Multispectral Image Super-Resolution Burned-Area Mapping Based on Space-Temperature Information

: Multispectral imaging (MI) provides important information for burned-area mapping. Due to the severe conditions of burned areas and the limitations of sensors, the resolution of collected multispectral images is sometimes very rough, hindering the accurate determination of burned areas. Super-resolution mapping (SRM) has been proposed for mapping burned areas in rough images to solve this problem, allowing super-resolution burned-area mapping (SRBAM). However, the existing SRBAM methods do not use su ﬃ ciently accurate space information and detailed temperature information. To improve the mapping accuracy of burned areas, an improved SRBAM method utilizing space–temperature information (STI) is proposed here. STI contains two elements, a space element and a temperature element. We utilized the random-walker algorithm (RWA) to characterize the space element, which encompassed accurate object space information, while the temperature element with rich temperature information was derived by calculating the normalized burn ratio (NBR). The two elements were then merged to produce an objective function with space–temperature information. The particle swarm optimization algorithm (PSOA) was employed to handle the objective function and derive the burned-area mapping results. The dataset of the Landsat-8 Operational Land Imager (OLI) from Denali National Park, Alaska, was used for testing and showed that the STI method is superior to the traditional SRBAM method.


Introduction
A challenging problem of the earth's ecosystem is wildland fires, which affect the balance of greenhouse gases, plant distribution, and inhabitant safety. The distribution of burned areas is fundamental for the study of wildland fires [1], which can be achieved using multispectral imaging (MI). However, due to the severe conditions of burned areas and the limitations of sensors, the resolution of collected multispectral images is sometimes very rough, presenting many mixed pixels, which hinders the accurate determination of burned areas [2]. Since there is more than one land-cover class in one mixed pixel, traditional classification technology that assigns one land-cover class to one pixel often cannot effectively deal with mixed pixels. Therefore, burned-area mapping using the classification results of coarse images is usually not ideal [3]. To address this issue, the super-resolution mapping (SRM) technology has been proposed to handle mixed pixels and obtain burned-area mapping, resulting in the method named super-resolution burned-area mapping (SRBAM). information with the temperature information. According to this objective function, the particle swarm optimization algorithm (PSOA) is employed to obtain the final burned-area mapping.
The innovations of STI are twofold: (1) We utilize RWA to characterize the space element with object space information. Since this takes into account the space information among and within objects, it is more comprehensive and accurate than the subpixel space information [3]. (2) The temperature information is utilized in the proposed STI by calculating the NBR. The experiments here reported using the Landsat-8 Operational Land Imager (OLI) dataset show the superiority of STI over other state-of-the-art methods.

Dataset
Lightning ignited the Castle Rocks fire in the deep backcountry of Denali National Park, Alaska, in July 2013. More than 12,900 acres were burned in two months. The fire caused great losses to the local ecosystem and economy. It is very important to obtain a fine space distribution of burned areas for firefighting and disaster relief. The experimental dataset used was an image of this area obtained by Landsat-8 OLI on 26 August 2013, which can be downloaded from the US Geological Survey (USGS) website: http://earthexplorer.usgs.gov/. The image has a size of 2968 × 2052 pixels, 30 m space resolution, and is centered at 64 • 31 N, 152 • 52 W. As shown in Figure 1a, the five visible main burned areas are marked in red due to the false color of the image. For quantitative evaluation, we needed a reference image derived from Figure 1a by a classification algorithm based on a least-squares support vector [42]. There are two class labels (burned area and background) in the reference image shown in Figure 1b. To highlight the burned area, the label is marked in red, and the background label is marked in black.
Remote Sens. 2019, 11, x FOR PEER REVIEW 3 of 20 temperature information is utilized in the proposed STI by calculating the NBR. The experiments here reported using the Landsat-8 Operational Land Imager (OLI) dataset show the superiority of STI over other state-of-the-art methods.

Dataset
Lightning ignited the Castle Rocks fire in the deep backcountry of Denali National Park, Alaska, in July 2013. More than 12,900 acres were burned in two months. The fire caused great losses to the local ecosystem and economy. It is very important to obtain a fine space distribution of burned areas for firefighting and disaster relief. The experimental dataset used was an image of this area obtained by Landsat-8 OLI on 26 August 2013, which can be downloaded from the US Geological Survey (USGS) website: http://earthexplorer.usgs.gov/. The image has a size of 2968 × 2052 pixels, 30 m space resolution, and is centered at 64°31'N, 152°52'W. As shown in Figure 1a, the five visible main burned areas are marked in red due to the false color of the image. For quantitative evaluation, we needed a reference image derived from Figure 1a by a classification algorithm based on a least-squares support vector [42]. There are two class labels (burned area and background) in the reference image shown in Figure 1b. To highlight the burned area, the label is marked in red, and the background label is marked in black.

Space Element
Here, we introduce the space element spa T with the object space information to obtain more accurate space information. Figure 2 shows the process of production of the space element. The rough multispectral image was upsampled by bicubic interpolation. A fractional image with burned-area classes of subpixel proportions was obtained by unmixing the upsampled image. The first principal component (PC) was then extracted from the upsampled image through principal component analysis (PCA). Because there was a lot of space information contained in this first PC, we segmented it to obtain objects using a multiresolution segmentation method [43]. Q is defined as the segmentation scale parameter, which determines the object size and the condition of merger termination. The segmentation method is given by

Space Element
Here, we introduce the space element T spa with the object space information to obtain more accurate space information. Figure 2 shows the process of production of the space element. The rough multispectral image was upsampled by bicubic interpolation. A fractional image with burned-area classes of subpixel proportions was obtained by unmixing the upsampled image.
where b represents a spectral band ( b = 1, 2, …, B ; with B being the total band number), b D is the th b band spectral value standard deviation in the object region, and spectral b λ is the free parameter here.
Among adjacent object regions, we merged two objects with the minimum difference. When H was larger than Q , we terminated the merging process and extracted the final objects.  The first principal component (PC) was then extracted from the upsampled image through principal component analysis (PCA). Because there was a lot of space information contained in this first PC, we segmented it to obtain objects using a multiresolution segmentation method [43]. Q is defined as the segmentation scale parameter, which determines the object size and the condition of merger termination. The segmentation method is given by where H represents regional differences, and λ is a free parameter to balance the shape difference H shape and the spectral difference H spectral . The shape difference H shape is calculated by where A is the actual frontier length of the object region, R is the rectangular frontier length of the object region, N is the subpixel number in the object region, A/ √ N and A/R represent the smoothness and compactness of the object region, respectively, and λ shape is a free parameter.
The spectral difference H spectral is defined by where b represents a spectral band (b = 1, 2, . . . , B; with B being the total band number), D b is the b th band spectral value standard deviation in the object region, and λ spectral b is the free parameter here. Among adjacent object regions, we merged two objects with the minimum difference. When H was larger than Q, we terminated the merging process and extracted the final objects.
Third, the space element T spa with the object space information was derived by RWA.
The space part T(i) spa corresponding to the i th subpixel was then obtained by RWA, according to Equation (5): where T among (G) is the space information among objects, T within (G) represents the space information within each object, ] is a column vector, and β is set to 0.5. T among (G) is given by where L is the Laplacian matrix: . The free parameter ε was set to 0.6, andv m was the m th object O m spectral value.v where v i is the spectral value of the i th subpixel in object O m . T within (G) is defined as where the diagonal values in the diagonal matrix Λ are the background class proportions, and the diagonal values in the diagonal matrix Λ represent the burned area class proportion. We minimized the space element T spa for all land-cover classes. The minimized formula is x(p i ) = 1, if subpixel p i belongs to burnt-area class 0, otherwise (10)

Temperature Element
A new temperature part T tem is here proposed to fully utilize the temperature information. T tem aims to minimize the difference in spectrum between the observed NBR value (NBR obe ) and the simulated NBR value (NBR sim ). The near-infrared (NIR) band and short-wave infrared band 1 (SWIR1) were used to calculate NBR obe here [41]: where the observed reflectance of both NIR band ρ obe NIR and SWIR1 band ρ obe SWIR1 are obtained directly from the original MI, and K is the number of mixed pixels.
Suppose r bur NIR and r bur SWIR1 are the reflectance of the burned area in the NIR and SWIR1 bands, and r non NIR and r non SWIR1 are the corresponding reflectance of the background. For each mixed pixel in these bands, the ratio of burned area subpixel number to total subpixel number is the proportion of burned area a bur NIR or a bur SWIR1 . The proportions of background in the 2 bands are then 1−a bur NIR and 1−a bur SWIR1 , respectively. We built a linear mixture including all subpixel spectra to consider the reflectance of each mixed pixel. Then, each mixed pixel-simulated reflectance in the NIR band ρ sim NIR and SWIR1 band ρ sim SWIR1 were calculated using Equations (13) and (14), respectively: NBR sim is given by: The temperature part T tem was then obtained by minimizing the difference between NBR obe and NBR sim :

Implementation of STI
To improve the burned area mapping result, STI is proposed as shown in Figure 3. It includes the following 3 steps: Step 1. Bicubic interpolation, segmentation, and RWA were utilized to obtain the space element T spa with more accurate space information. At the same time, the temperature element T tem , which contains rich temperature information, was obtained by calculating the NBR.
Step 2. We merged the space element T spa and the temperature element T tem through a trade-off parameter θ to produce the objective function T with space-temperature information. The aim of the proposed STI is to minimize T. In STI, we consider a weighted sum of the space and temperature elements of STI, because this information-fusion method has a simple physical meaning and is easy to implement. Of course, we can also use other more effective information-fusion techniques, such as multiobjective optimization [44], alpha integration [45], and so on.
Step 3. To optimize the objective function, PSOA was employed. First, we randomly assigned a burned area or background label to all subpixels. Second, the labels of these subpixels were iteratively changed until the minimum value of T was derived. During each iteration, the burned-area label was changed to the background label and vice versa. If T increased, the change was rejected, otherwise it was accepted. When less than 0.1% of labels were changed, the PSOA was terminated.
The temperature part tem T was then obtained by minimizing the difference between obe NBR and sim NBR :

Implementation of STI
To improve the burned area mapping result, STI is proposed as shown in Figure 3. It includes the following 3 steps: Step 1. Bicubic interpolation, segmentation, and RWA were utilized to obtain the space element spa T with more accurate space information. At the same time, the temperature element tem T , which contains rich temperature information, was obtained by calculating the NBR.
Step 2. We merged the space element spa T and the temperature element tem T through a tradeoff parameter θ to produce the objective function T with space-temperature information. The aim of the proposed STI is to minimize T . In STI, we consider a weighted sum of the space and temperature elements of STI, because this information-fusion method has a simple physical meaning and is easy to implement. Of course, we can also use other more effective information-fusion techniques, such as multiobjective optimization [44], alpha integration [45], and so on.
Step 3. To optimize the objective function, PSOA was employed. First, we randomly assigned a burned area or background label to all subpixels. Second, the labels of these subpixels were iteratively changed until the minimum value of T was derived. During each iteration, the burned-area label was changed to the background label and vice versa. If T increased, the change was rejected, otherwise it was accepted. When less than 0.1% of labels were changed, the PSOA was terminated.

Experimental Settings
Fine images of five visible main burned areas from the experimental dataset are shown in Figure 4. The test sizes of the five burned areas were 720 × 720 pixels, 300 × 300 pixels, 720 × 720 pixels, 400 × 400 pixels, and 500 × 500 pixels. A flowchart of the experimental process is shown in Figure 5. We used the most commonly used experimental process of SRM to conduct the experiments. The five visible main burned areas were downsampled via an S × S mean filter to produce a rough multispectral image. Here, scale S was set to 8, namely 8 × 8 pixels in the original fine image were merged into one mixed pixel in the simulated rough image. In this case, we could directly evaluate the impact of the error of image registration on SRM. In addition, a quantitative evaluation could be carried out more reasonably in this way; a reference image could be derived from the classification result of the fine image, which was compared with the SRM result from the simulated rough image. Rough images of the five burned areas are shown in Figure 6. Although the false color rough image could highlight the burned area, it was difficult to obtain more accurate distribution and boundary information of the burned area due to the rough resolution. For example, because of the presence of many mixed pixels at the edge of the burned area, there is an obvious vertical line at the edge of the burned area in Figure 6d compared with Figure 4d. In addition, it was difficult for the classification technology to handle the mixed pixels, because one mixed pixel contained more than one land-cover class. To solve this problem, SRM was utilized to handle the mixed pixels and produce accurate burned-area mapping. The least-squares linear mixture model (LSLMM) [46] was applied to the rough images to derive fractional images as inputs. In the segmentation method, the selected λ, λ shape , and λ

Experimental Settings
Fine images of five visible main burned areas from the experimental dataset are shown in Figure  4. The test sizes of the five burned areas were 720 × 720 pixels, 300 × 300 pixels, 720 × 720 pixels, 400 × 400 pixels, and 500 × 500 pixels. A flowchart of the experimental process is shown in Figure 5. We used the most commonly used experimental process of SRM to conduct the experiments. The five visible main burned areas were downsampled via an S S × mean filter to produce a rough multispectral image. Here, scale S was set to 8, namely 8 8 × pixels in the original fine image were merged into one mixed pixel in the simulated rough image. In this case, we could directly evaluate the impact of the error of image registration on SRM. In addition, a quantitative evaluation could be carried out more reasonably in this way; a reference image could be derived from the classification result of the fine image, which was compared with the SRM result from the simulated rough image. Rough images of the five burned areas are shown in Figure 6. Although the false color rough image could highlight the burned area, it was difficult to obtain more accurate distribution and boundary information of the burned area due to the rough resolution. For example, because of the presence of many mixed pixels at the edge of the burned area, there is an obvious vertical line at the edge of the burned area in Figure 6d compared with Figure 4d. In addition, it was difficult for the classification technology to handle the mixed pixels, because one mixed pixel contained more than one land-cover class. To solve this problem, SRM was utilized to handle the mixed pixels and produce accurate burned-area mapping. The least-squares linear mixture model (LSLMM) [46] was applied to the rough images to derive fractional images as inputs. In the segmentation method, the selected λ , We tested four SRM methods: the hybrid spatial attraction model (HSAM) [15], the object-scale spatial SRM (OSRM) [39], the SRBAM [3], and the proposed STI. The ratio between the number of correct mapping subpixels belonging to burned areas derived from each SRM result and the total number of subpixels belonging to burned areas derived from the reference image was defined as burned area (%), and the ratio between the number of correct mapping subpixels derived from each SRM result belonging to background and the total number of subpixels belonging to background derived from the reference image was defined as background (%). The four methods were evaluated on the basis of the determination accuracy of each class (burned area (%) and background (%)), overall accuracy (OA (%)), and kappa coefficient (Kappa) [5].

Results Analysis
First, a visual comparison was performed. The burned-area mapping results of the SRM method in the five test areas are given in Figures 7-11. The detailed area is marked in a rectangular white frame. When we compared the reference images with the four experimental results, we found that STI outperformed the other three SRM methods, and the results from STI were more similar to the reference images. For burned areas with complex distribution, such as Areas 1, 3, 4, and 5, there were many disconnected patches, and some small areas disappeared in the results of HSAM, OSRM, and SRBAM. For some areas with simple distribution, such as Area 2, there were many obvious burrs at the boundaries between burned area in the results of HSAM, OSRM, and SRBAM. There are two reasons for these phenomena: First, the space information was not accurate enough. HSAM and We tested four SRM methods: the hybrid spatial attraction model (HSAM) [15], the object-scale spatial SRM (OSRM) [39], the SRBAM [3], and the proposed STI. The ratio between the number of correct mapping subpixels belonging to burned areas derived from each SRM result and the total number of subpixels belonging to burned areas derived from the reference image was defined as burned area (%), and the ratio between the number of correct mapping subpixels derived from each Remote Sens. 2019, 11, 2695 9 of 18 SRM result belonging to background and the total number of subpixels belonging to background derived from the reference image was defined as background (%). The four methods were evaluated on the basis of the determination accuracy of each class (burned area (%) and background (%)), overall accuracy (OA (%)), and kappa coefficient (Kappa) [5].

Results Analysis
First, a visual comparison was performed. The burned-area mapping results of the SRM method in the five test areas are given in Figures 7-11. The detailed area is marked in a rectangular white frame. When we compared the reference images with the four experimental results, we found that STI outperformed the other three SRM methods, and the results from STI were more similar to the reference images. For burned areas with complex distribution, such as Areas 1, 3, 4, and 5, there were many disconnected patches, and some small areas disappeared in the results of HSAM, OSRM, and SRBAM. For some areas with simple distribution, such as Area 2, there were many obvious burrs at the boundaries between burned area in the results of HSAM, OSRM, and SRBAM. There are two reasons for these phenomena: First, the space information was not accurate enough. HSAM and SRBAM only consider pixel-level space information, which is less detailed than object-level space information. Although OSRM utilizes object-level space information, it only calculates space information among object regions and does not consider space information within object regions. Since the proposed STI utilizes object-level space information among and within object regions through RWA, the space information is more accurate in STI than in the other three methods. In addition, STI is better able to make full use of temperature information than the other methods, so it obtains better burned-area mapping results.

Results Analysis
First, a visual comparison was performed. The burned-area mapping results of the SRM method in the five test areas are given in Figures 7-11. The detailed area is marked in a rectangular white frame. When we compared the reference images with the four experimental results, we found that STI outperformed the other three SRM methods, and the results from STI were more similar to the reference images. For burned areas with complex distribution, such as Areas 1, 3, 4, and 5, there were many disconnected patches, and some small areas disappeared in the results of HSAM, OSRM, and SRBAM. For some areas with simple distribution, such as Area 2, there were many obvious burrs at the boundaries between burned area in the results of HSAM, OSRM, and SRBAM. There are two reasons for these phenomena: First, the space information was not accurate enough. HSAM and SRBAM only consider pixel-level space information, which is less detailed than object-level space information. Although OSRM utilizes object-level space information, it only calculates space information among object regions and does not consider space information within object regions. Since the proposed STI utilizes object-level space information among and within object regions through RWA, the space information is more accurate in STI than in the other three methods. In addition, STI is better able to make full use of temperature information than the other methods, so it obtains better burned-area mapping results.    Second, we analyzed the accuracy evaluation index. The performance of the four SRM methods was evaluated on the basis of burned area (%), background (%), OA (%), and Kappa. Checking the evaluating indicators in Table 1, the burned area (%) measured by STI was higher than that measured by the other three methods. Compared with SRBAM, the burned area (%) of STI increased by 3.29%, 4.73%, 2.85%, 0.63% and 3.55% in the five test areas. With the aid of space-temperature information, the proposed STI produced the highest OA (%) and Kappa. Second, we analyzed the accuracy evaluation index. The performance of the four SRM methods was evaluated on the basis of burned area (%), background (%), OA (%), and Kappa. Checking the evaluating indicators in Table 1, the burned area (%) measured by STI was higher than that measured by the other three methods. Compared with SRBAM, the burned area (%) of STI increased by 3.29%, 4.73%, 2.85%, 0.63% and 3.55% in the five test areas. With the aid of space-temperature information, the proposed STI produced the highest OA (%) and Kappa.  Second, we analyzed the accuracy evaluation index. The performance of the four SRM methods was evaluated on the basis of burned area (%), background (%), OA (%), and Kappa. Checking the evaluating indicators in Table 1, the burned area (%) measured by STI was higher than that measured by the other three methods. Compared with SRBAM, the burned area (%) of STI increased by 3.29%, 4.73%, 2.85%, 0.63% and 3.55% in the five test areas. With the aid of space-temperature information, the proposed STI produced the highest OA (%) and Kappa. Third, we tested the performance of SRM by different scales S, which represent the simulated rough images with different resolution as inputs. The scales S set to different values confirmed that the STI still had the best performance for inputs with different resolutions. HSAM, OSRM, SRBAM, and STI were tested using the other two scales (5 and 10) in the five test areas. The burned area (%) of these methods in relation to S = 5 and S = 10 are shown in Figure 12. We found that, as the value of S increased, the burned area (%) determined by the four methods decreased. This is because as S became larger, the input image became rougher, producing a greater challenge to SRM. The experimental results showed that STI still determined the highest burned area (%) when using different scales S.
Fourth, the influence of the selected parameter θ on the proposed method was studied. Five test areas (S = 8) were rerun for 10 combinations of θ from 0 to 0.9, at an interval of 0.1. The results are shown in Figure 13. There was no contribution of the temperature element T tem when θ = 0. At this time, only the space element T spa was working, therefore the value of burned area (%) was low. As θ increased, the burned area (%) increased. This is because the use of temperature information from the temperature element T tem increased as θ increased. When θ = 0.4, θ = 0.5, θ = 0.4, θ = 0.6, and θ = 0.4 in the five test areas, the burned area (%) reached its highest value. At this time, the contributions of the space term T spa and the temperature term T tem reached a state of balance. However, when θ increased, the space term T spa reduced its contribution to Equation (17). The burned-area mapping accuracy was affected as a consequence of the decreased space information from the space term T spa . The parameter θ in STI required several experiments to be determines, that is, adaptability is not ideal in STI.  Figure 13. There was no contribution of the temperature element tem T when =0 θ . At this time, only the space element spa T was working, therefore the value of burned area (%) was low.
As θ increased, the burned area (%) increased. This is because the use of temperature information from the temperature element tem T increased as θ increased. When =0.4 θ , =0.5 θ , =0.4 θ , =0.6 θ , and =0.4 θ in the five test areas, the burned area (%) reached its highest value. At this time, the contributions of the space term spa T and the temperature term tem T reached a state of balance.
However, when θ increased, the space term spa T reduced its contribution to Equation (17). The burned-area mapping accuracy was affected as a consequence of the decreased space information from the space term spa T . The parameter θ in STI required several experiments to be determines, that is, adaptability is not ideal in STI. Fifth, the impact of the segmentation scale parameter Q on the proposed method was studied.
In the proposed STI, the space element spa T was obtained by calculating the class proportion of objects through RWA. Therefore, the step of segmentation that produces the object is very important for STI, and the quality of the objects is determined by the segmentation scale parameter Q in the Fifth, the impact of the segmentation scale parameter Q on the proposed method was studied. In the proposed STI, the space element T spa was obtained by calculating the class proportion of objects through RWA. Therefore, the step of segmentation that produces the object is very important for STI, and the quality of the objects is determined by the segmentation scale parameter Q in the segmentation method. Because Q determines the object size and the condition of merger termination, we studied the optimal selection of Q in this experiment. Ten Q values from 5 to 50 at an interval of 5 were tested in the five test areas (S = 8). As shown in Figure 14, the selection of Q had an impact on the final mapping accuracy. When the value of Q was not properly selected, the burned area (%) was low. This is because an inappropriate Q resulted in low-quality objects, which affected the accuracy of space information in the space element T spa . After many experiments, it was noted that the best Q values of the five test areas were 15, 10, 15, 20, and 15. The segmentation scale parameter Q also required many experiments to be determined; therefore, this also proves that the adaptability of STI is not ideal.  Finally, we analyzed the operation time (s). Figure 15 shows the operation time (s) of the four SRM methods in the five test areas (S = 8). The results showed that STI required the longest time. This is because the proposed STI is characterized by a more complex processing. Although STI requires more computation time than the other SRM methods, it shows improved performance. Therefore, the long computation and running times are disadvantages of the proposed method.

Conclusion
Often, traditional classification technology cannot effectively deal with mixed pixels in rough multispectral images. In order to obtain accurate distribution results of land-cover classes in rough images, the SRM technology has been proposed. SRM can produce better mapping results than traditional classification methods when processing rough images. In this paper, STI is proposed to improve burned-area mapping by fully utilizing the space-temperature information of burned areas. The space element and the temperature element are used in STI. The RWA is used to compute the

Conclusions
Often, traditional classification technology cannot effectively deal with mixed pixels in rough multispectral images. In order to obtain accurate distribution results of land-cover classes in rough images, the SRM technology has been proposed. SRM can produce better mapping results than traditional classification methods when processing rough images. In this paper, STI is proposed to improve burned-area mapping by fully utilizing the space-temperature information of burned areas. The space element and the temperature element are used in STI. The RWA is used to compute the segmented objects to obtain the space element with accurate and comprehensive space information. At the same time, the temperature element with full temperature information is obtained by calculating the difference between NBR obe and NBR sim . An objective element with the space-temperature information is derived by integrating the space element and the temperature element. Finally, the PSOA is utilized to optimize the objective element to produce burned-area mapping results. Thanks to the space-temperature information, the proposed STI obtains better burned-area mapping results than the existing SRBAM. Experiments on Landsat-8 OLI images of burned areas in Denali National Park, Alaska, showed that STI produced the highest OA (%), achieving 92.39%, 94.21%, 94.42%, 99.01%, and 95.48% OA in the five tested areas. Although it seems that STI performed only 1-2% better in OA (%) than the other SRM methods, in fact it successfully corrected thousands of pixels. For example, the OA (%) of STI was around 1.5% greater than that of SRBAM in Area 1. According to the definition of OA (%), since Area 1 had 720 × 720 pixels, the number of corrected pixels obtained by STI included about 7776 pixels more than that that obtained by SRBAM. Therefore, the gain in accuracy of the proposed method is obvious.
The appropriate values of the parameters θ and Q were selected by multiple tests when using STI. To improve the final mapping results, it is worth studying an adaptive method for selecting the most appropriate value of the parameters θ and Q in future work. In addition, it is worth further studying how to simplify the structure of the proposed STI and improve the running speed. Finally, the use of new artificial intelligence technology and a large number of auxiliary data to improve STI and obtain OA values closer to 100% is also worth studying in the future.