Topographic Correction to Landsat Imagery through Slope Classiﬁcation by Applying the SCS + C Method in Mountainous Forest Areas

: The aim of the topographic normalization of remotely sensed imagery is to reduce reﬂectance variability caused by steep terrain and thus improve further processing of images. A process of topographic correction was applied to Landsat imagery in a mountainous forest area in the south of Mexico. The method used was the Sun Canopy Sensor + C correction (SCS + C) where the C parameter was differently determined according to a classiﬁcation of the topographic slopes of the studied area in nine classes for each band, instead of using a single C parameter for each band. A comparative, visual, and numerical analysis of the normalized reﬂectance was performed based on the corrected images. The results showed that the correction by slope classiﬁcation improves the elimination of the effect of shadows and relief, especially in steep slope areas, modifying the normalized reﬂectance values according to the combination of slope, aspect, and solar geometry, obtaining reﬂectance values more suitable than the correction by non-slope classiﬁcation. The application of the proposed method can be generalized, improving its performance in forest mountainous areas.


Background
The determination of the reflectance value of the surface by removing topographic effects is an important task in remote sensing studies [1]. However, the ability to provide reliable information through satellite imagery is constrained by the effects that changes in slope and terrain orientation angles (aspect) in combination to the solar geometry (zenithal and azimuthal angles) at the time of the image acquisition cause on the spectral direct and diffuse irradiance [2].
This occurs especially in mountainous areas where slope surfaces directly oriented to sun rays receive more light and consequently appear brighter in images than those surfaces that are not receiving the sunlight directly [3]. C parameter may not produce an accurate topographic correction for the whole scene because of slope variation.
The aim of this study is to assess, qualitatively and quantitatively, the results of topographic correction of Landsat images through the SCS + C method, by non-slope classification (N-SC) and by slope classification (SC). To achieve the proposed aim, a comparative analysis between reflectance images normalized by the standard SCS + C topographic correction method (with a single global C parameter) and reflectance images normalized by the same approach but where nine C parameters were determined by applying a classification based on the slopes of the studied zone.
For this work, Dinamica EGO software [30][31][32] was used to develop the described models and processes.

Study Area and Dataset
The studied area consists of a group of basins in the central region of the Guerrero state in Mexico, around the largest gold mine in Latin America, "Filos-Bermejal" [33], that covers nearly 4000 km 2 of mountain terrain with altitudes from 580 to 2840 m above the sea level, where slopes greater than 40 • are present ( Figure 1). The aim of this study is to assess, qualitatively and quantitatively, the results of topographic correction of Landsat images through the SCS + C method, by non-slope classification (N-SC) and by slope classification (SC). To achieve the proposed aim, a comparative analysis between reflectance images normalized by the standard SCS + C topographic correction method (with a single global C parameter) and reflectance images normalized by the same approach but where nine C parameters were determined by applying a classification based on the slopes of the studied zone.
For this work, Dinamica EGO software [30][31][32] was used to develop the described models and processes.

Study Area and Dataset
The studied area consists of a group of basins in the central region of the Guerrero state in Mexico, around the largest gold mine in Latin America, "Filos-Bermejal" [33], that covers nearly 4000 km 2 of mountain terrain with altitudes from 580 to 2840 m above the sea level, where slopes greater than 40° are present ( Figure 1). The land cover distribution of the zone is 45% of oak, pine, and mesophyll forest, 36% of deciduous forest, and 19% of grassland, agricultural, bare soil, water, and urban area [34,35].
The vegetation types refer to the different groups that are presented in the zone and that are part of the large of vegetation ecological information, floristics, and physiognomy. The types of vegetation included as land cover are as follows [36]: Oak, pine, and mesophyll forest, including secondary bushy vegetation that corresponds to areas where the particular ecological original conditions of the community have been altered and where each plant community has a group of minor species that cover the altered space.
The oak forest refers to arboreal, sub-arboreal, or occasional shrub communities integrated by multiple species of the genus Quercus that, in Mexico, except in very arid conditions, are located at altitudes from 300 to 2800 m. It is closely related to pine forests, forming mixed forests with species of both genera.
The pine forest includes arboreal, sub-arboreal, or occasionally shrub communities of northern (Holarctic) origin, mainly of milder and middle-cold weather regions with different degrees of humidity. Usually with minimal variation of species of conifers. Its development is a consequence of the climate and the soil of a region, which have not significantly influenced other factors for its establishment. The land cover distribution of the zone is 45% of oak, pine, and mesophyll forest, 36% of deciduous forest, and 19% of grassland, agricultural, bare soil, water, and urban area [34,35].
The vegetation types refer to the different groups that are presented in the zone and that are part of the large of vegetation ecological information, floristics, and physiognomy. The types of vegetation included as land cover are as follows [36]: Oak, pine, and mesophyll forest, including secondary bushy vegetation that corresponds to areas where the particular ecological original conditions of the community have been altered and where each plant community has a group of minor species that cover the altered space.
The oak forest refers to arboreal, sub-arboreal, or occasional shrub communities integrated by multiple species of the genus Quercus that, in Mexico, except in very arid conditions, are located at altitudes from 300 to 2800 m. It is closely related to pine forests, forming mixed forests with species of both genera.
The pine forest includes arboreal, sub-arboreal, or occasionally shrub communities of northern (Holarctic) origin, mainly of milder and middle-cold weather regions with different degrees of humidity. Usually with minimal variation of species of conifers. Its development is a consequence of the climate and the soil of a region, which have not significantly influenced other factors for its establishment.
The mesophyll forest refers to plant communities characterized by the presence of dense to very dense tree vegetation, with epiphytes and ferns. These are located mainly in mountains, ravines, and sites that present favorable conditions of humidity and frequent fogs. In Mexico, these are located at altitudes between 600 m and 2800 m.
The deciduous forest refers to arboreal or sub-arboreal communities of tropical origin that grow in places with seasonal precipitation and where more than 75% of their plant components lose their leaves during the dry season of the year.
The SR information is produced by the specialized software "Landsat Ecosystem Disturbance Adaptive Processing System" for Landsat 5-TM and L8SR for Landsat 8-OLI, developed by the Making Earth System Data Records for Use in Research Environments of the Goddard Space Flight Center of NASA and the University of Maryland [40].
The software applies the MODIS atmospheric correction routines. In addition to the Landsat data, the water vapor, ozone, geopotential height, aerosol optical thickness, and digital elevation data are presented to the radiative transfer model "Second Simulation of a Satellite Signal in the Solar Spectrum" (6S), to generate the top of atmosphere reflectance, SR, brightness temperature and cloud masks, cloud shadows, and adjacent clouds, land, and water.
The land cover and vegetation map "Series I to IV" scale 1:250,000 were used, and a 30 m grid cell size DEM, which was generated from the vectorial data (contours) of the topographic maps, courtesy of the National Institute of Statistics and Geography of Mexico (INEGI), with a scale of 1:50,000 and 20 m for vertical accuracy.

Overview
The topographic corrections by SCS + C method were performed in two stages: first with a global C parameter determined by N-SC considering all pixels, and second with a set of nine C parameters determined by SC considering all those pixels with a slope greater than 0 • (Figure 2). The mesophyll forest refers to plant communities characterized by the presence of dense to very dense tree vegetation, with epiphytes and ferns. These are located mainly in mountains, ravines, and sites that present favorable conditions of humidity and frequent fogs. In Mexico, these are located at altitudes between 600 m and 2800 m.
The deciduous forest refers to arboreal or sub-arboreal communities of tropical origin that grow in places with seasonal precipitation and where more than 75% of their plant components lose their leaves during the dry season of the year.
The SR information is produced by the specialized software "Landsat Ecosystem Disturbance Adaptive Processing System" for Landsat 5-TM and L8SR for Landsat 8-OLI, developed by the Making Earth System Data Records for Use in Research Environments of the Goddard Space Flight Center of NASA and the University of Maryland [40].
The software applies the MODIS atmospheric correction routines. In addition to the Landsat data, the water vapor, ozone, geopotential height, aerosol optical thickness, and digital elevation data are presented to the radiative transfer model "Second Simulation of a Satellite Signal in the Solar Spectrum" (6S), to generate the top of atmosphere reflectance, SR, brightness temperature and cloud masks, cloud shadows, and adjacent clouds, land, and water.
The land cover and vegetation map "Series I to IV" scale 1:250,000 were used, and a 30 m grid cell size DEM, which was generated from the vectorial data (contours) of the topographic maps, courtesy of the National Institute of Statistics and Geography of Mexico (INEGI), with a scale of 1:50,000 and 20 m for vertical accuracy.

Overview
The topographic corrections by SCS + C method were performed in two stages: first with a global C parameter determined by N-SC considering all pixels, and second with a set of nine C parameters determined by SC considering all those pixels with a slope greater than 0° ( Figure 2).  The resulting corrected images were visually (qualitatively) explored and numerically (quantitatively) validated based on the normalized reflectance values obtained.
In concordance with previous works showing a close connection between the red band and the chlorophyll content and the near-infrared band linked to the density of vegetation [41,42], the resulting C parameters are presented and discussed for the red and near-infrared bands. Nevertheless, the reflectance values resulting from the topographic corrections are shown for each date, band, and slope class.

Correction C Parameter Calculation
The C parameter is determined by C = b/m, where b is the intercept and m is the gradient of the line obtained from the linear regression between the cosine of the solar incidence angle i (illumination) as the independent variable, and spectral data (reflectance) as the dependent variable.
The cosine of the incident solar angle (cos i), referred to as illumination [3,10,43], was calculated by where θp is the terrain slope; θi is the solar zenith; Φa is the solar azimuth; Φ0 is the aspect. The slope and aspect data were derived from the DEM, and the solar zenith and azimuth angles were obtained from the metadata files.
In the first stage, a global C(i) parameter was calculated (where i corresponds to the band), through a linear regression (illumination-reflectance) and considering all pixels, regardless of the slope variations.
In the second stage, due to the wide range of terrain slopes in the studied area, the C(i,j) parameter was defined for each band (i), considering a slope classification grouped into nine classes (j), with ranges of 5 • each, for those pixels greater than 0 • . In this case, for each slope class, a linear regression was performed just for the pixels considered in the range of the slope classes and excluding the rest.
Thus, for each band (i) and date image, first a single C(i) parameter value by N-SC was generated, and nine C(i,j) parameter values (one per class) by SC were later generated.

Topographic Correction
Soenen et al. [6] used the SCS + C method incorporating the C parameter to the SCS method, to characterize the diffuse atmospheric irradiance in a better way. The following equation expresses the model: Ln = L (cos θp cos θi + C)/(cos i + C) (2) where Ln is the equivalent flat surface reflectance (normalized), L is the measured radiance in the remotely sensed data, and C is the semi-empirical parameter.
With the C(i,j) parameters obtained in each case (SC and N-SC), the SCS + C topographic correction method was applied, and normalized-reflectance images have been achieved for each band. In the case of topographic correction by SC, this process was performed for each class, and the resulting nine normalized-reflectance individual images (one per class) were later combined in a global final single image for each band.

Qualitative Validation
The visual exploration of images in remote sensing studies is a common method of assessing the quality of results [13,44]. The quantitative validation essentially consists of developing a visual comparative analysis between the reflectance images before topographic corrections, and images after the application of topographic correction either by N-SC or by SC; however, this is a subjective evaluation that can be inconvenient [45].

Quantitative Validation
To overcome the qualitative validation limitation, additional quantitative validations were performed: First, the effect of corrections over "oak, pine, and mesophyll forest," "deciduous forest," and "bare soil" land cover samples were analyzed, all of which had a broad range of slopes, enabling the identification of all slope classes used in the topographic correction process, and a wide variety of terrain orientations (aspect) were considered to have a broad range of shades as well.
For each date of the study, selected land cover, band, and slope class; the mean reflectance values before and after applying topographical corrections were obtained, and an analysis of standard deviation was performed.
In a second validation, it was considered that the topographic correction is based on the terrain slope, aspect, and geometry of the sun at the time of image acquisition; therefore, if the slope is null (flat terrain), the corrected reflectance values obtained by the SCS + C model equations are equal to the originally measured reflectance values.
Under this consideration for each date, land cover, band, and slope class, the absolute and relative differences between the mean normalized reflectance values of slope classes greater than 0 • and mean normalized reflectance values of those pixels with a slope equal to 0 • were calculated, considering these last pixels as topographically corrected, and an analysis of these differences was then performed.
To assess the relationship between the illumination (cos i) and reflectance values, an analysis through scatterplot graphs was conducted to the forest cover sample, before and after applying the topographic corrections (SC and N-SC), for each slope class.

C(i,j) Parameters
The variation of C(i,j) parameters obtained from the red and near-infrared bands considering SC can be seen (Figure 3), even becoming negative in the lower slope classes, unlike the global parameter determined by N-SC, the results of which are always positive. the identification of all slope classes used in the topographic correction process, and a wide variety of terrain orientations (aspect) were considered to have a broad range of shades as well.
For each date of the study, selected land cover, band, and slope class; the mean reflectance values before and after applying topographical corrections were obtained, and an analysis of standard deviation was performed.
In a second validation, it was considered that the topographic correction is based on the terrain slope, aspect, and geometry of the sun at the time of image acquisition; therefore, if the slope is null (flat terrain), the corrected reflectance values obtained by the SCS + C model equations are equal to the originally measured reflectance values.
Under this consideration for each date, land cover, band, and slope class, the absolute and relative differences between the mean normalized reflectance values of slope classes greater than 0° and mean normalized reflectance values of those pixels with a slope equal to 0° were calculated, considering these last pixels as topographically corrected, and an analysis of these differences was then performed.
To assess the relationship between the illumination (cos i) and reflectance values, an analysis through scatterplot graphs was conducted to the forest cover sample, before and after applying the topographic corrections (SC and N-SC), for each slope class.

C(i,j) Parameters
The variation of C(i,j) parameters obtained from the red and near-infrared bands considering SC can be seen (Figure 3), even becoming negative in the lower slope classes, unlike the global parameter determined by N-SC, the results of which are always positive. It can be noted that the C(i,j) parameters by SC increase, according to the upper slopes, and critical values can be identified for those slopes where the values of both (N-SC and SC) C parameters are equal (at crossing lines).
It can also be noted that the patterns differ with the graphs obtained from the red band (Figure 3ac) versus the graphs obtained for the infrared band (Figure 3d-f). For infrared graphs, the global C(i) parameters calculated by N-SC have less variation compared to those calculated by SC, which follows that the resulting reflectance values after N-SC and SC corrections are more similar to each other. It can be noted that the C(i,j) parameters by SC increase, according to the upper slopes, and critical values can be identified for those slopes where the values of both (N-SC and SC) C parameters are equal (at crossing lines).
It can also be noted that the patterns differ with the graphs obtained from the red band (Figure 3a-c) versus the graphs obtained for the infrared band (Figure 3d-f). For infrared graphs, the global C(i) parameters calculated by N-SC have less variation compared to those calculated by SC, which follows that the resulting reflectance values after N-SC and SC corrections are more similar to each other.
This difference between C parameters calculated by SC or N-SC undoubtedly produces different results in the normalized reflectance obtained by the topographic corrections, as shown below.
All C parameters (N-SC) and by slope classes (SC) resulting for each band, slope class, and date of study are included in Supplementary Materials Table S1.

Qualitative Validation
The elimination of the shadows caused by relief and sun geometry at the time of image acquisition from the original image (Figure 4a) can be confirmed after the application of topographic correction processes (Figure 4c,d). This difference between C parameters calculated by SC or N-SC undoubtedly produces different results in the normalized reflectance obtained by the topographic corrections, as shown below.
All C parameters (N-SC) and by slope classes (SC) resulting for each band, slope class, and date of study are included in Supplementary Materials Table S1.

Qualitative Validation
The elimination of the shadows caused by relief and sun geometry at the time of image acquisition from the original image (Figure 4a) can be confirmed after the application of topographic correction processes (Figure 4c,d). Visual results of both corrections are similar for areas with slopes that are not very steep; however, for areas that combine higher values of shadows and slope (greater than 40° as rounded in Figure 4c,d), it can be noted that the corrected image by SC (Figure 4d) shows a visual improvement compared to the corrected image by N-SC (Figure 4c). Table 1 shows the standard deviation values of each band and slope class before and after the topographic corrections, including by each band, the global mean deviation standard value. Visual results of both corrections are similar for areas with slopes that are not very steep; however, for areas that combine higher values of shadows and slope (greater than 40 • as rounded in Figure 4c,d), it can be noted that the corrected image by SC (Figure 4d) shows a visual improvement compared to the corrected image by N-SC (Figure 4c). Table 1 shows the standard deviation values of each band and slope class before and after the topographic corrections, including by each band, the global mean deviation standard value. According to Hantson and Chuvieco [46], it is expected that the standard deviation values decrease after successful shade removal, meaning that the impact of illumination has been reduced. Thus, based on the lower standard deviation, it can be seen that the values obtained by SC correction are smaller than the values of N-SC correction, indicating an improvement.

Quantitative Validation
In Supplementary Materials Table S2 and S3, the mean standard deviation values before and after the topographic correction for each band, slope class, and date of study have been included for the land covers of forest (Supplementary Materials Table S2) and deciduous forest (Supplementary  Materials Table S3).
In the corresponding standard deviation graph (Figure 5a), this improvement can be depicted for oak, pine, and mesophyll forest areas. However, if the same task is performed for deciduous forest and bare soil areas, different patterns of standard deviations can be obtained (Figure 5b,c).
Comparing the patterns before and after corrections (red versus green and black lines) on different land covers can be seen the improvement of the application of the topographic correction. Additionally, the patterns depicted show an improvement in SC correction in vegetation areas (Figure 5a,d,g and 5b,e,h) since, in all bands, the standard deviation values are lower (black line).
On the other hand, for bare soil areas (Figure 5c,f,h), it can be seen that this improvement is not so obvious, as even the correction by no-slope classification seems slightly better in the Landsat OLI image from 7 March 2015 (Figure 5i). This fact may seem evident since, according to previous works, the SCS + C topographic correction model is the appropriated framework for forested scenes [6,15,16]. Regarding this fact, Gu and Gillespie [17] suggested the performance of separate corrections using a sun-canopy-sensor model for forest areas and the C-correction model for non-forest areas [2].
It can be seen in the graphs that bare soil land covers similar patterns (Figure 5c,f,i). However, for the vegetation land cover areas (forest and deciduous forest), it can be seen that the depicted patterns show some differences between the applied topographic corrections (N-SC: green line; SC: black line). According to Hantson and Chuvieco [46], it is expected that the standard deviation values decrease after successful shade removal, meaning that the impact of illumination has been reduced. Thus, based on the lower standard deviation, it can be seen that the values obtained by SC correction are smaller than the values of N-SC correction, indicating an improvement.
In Supplementary Materials Table S2 and S3, the mean standard deviation values before and after the topographic correction for each band, slope class, and date of study have been included for the land covers of forest (Supplementary Materials Table S2) and deciduous forest (Supplementary  Materials Table S3).
In the corresponding standard deviation graph (Figure 5a), this improvement can be depicted for oak, pine, and mesophyll forest areas. However, if the same task is performed for deciduous forest and bare soil areas, different patterns of standard deviations can be obtained (Figure 5b,c).
Comparing the patterns before and after corrections (red versus green and black lines) on different land covers can be seen the improvement of the application of the topographic correction. Additionally, the patterns depicted show an improvement in SC correction in vegetation areas (Figure 5a,d,g and 5b,e,h) since, in all bands, the standard deviation values are lower (black line).
On the other hand, for bare soil areas (Figure 5c,f,h), it can be seen that this improvement is not so obvious, as even the correction by no-slope classification seems slightly better in the Landsat OLI image from 7 March 2015 (Figure 5i). This fact may seem evident since, according to previous works, the SCS + C topographic correction model is the appropriated framework for forested scenes [6,15,16]. Regarding this fact, Gu and Gillespie [17] suggested the performance of separate corrections using a sun-canopy-sensor model for forest areas and the C-correction model for non-forest areas [2].
It can be seen in the graphs that bare soil land covers similar patterns (Figure 5c,f,i). However, for the vegetation land cover areas (forest and deciduous forest), it can be seen that the depicted patterns show some differences between the applied topographic corrections (N-SC: green line; SC: black line).  For the deciduous forest land cover, although similar global patterns are observed, some individual differences between the study dates (sensors) can be seen, while for the forest land cover, the difference can be viewed in the Landsat TM image from 27 November 1995.
These differences may be due to the angular values of solar zenith (higher for the TM image than the other study dates and sensors). Additionally, it can be seen that the differences are observed in the shortwave infrared bands (SWIR-1 and SWIR-2). Thus, it is possible that this difference is due to sensor calibration.
The scatterplots (Figure 6) before and after corrections for each class confirm, first, the improvement in normalization of reflectance after the performance of the topographic correction.
Although the scatterplot graphs corresponds to a forest sample, the general distribution of the slopes in the study area can be seen, the majority being concentrated within a range of 20-35 • .
In the first column ( Figure 6) corresponding to the scatterplots before the topographic correction, it can be seen that the correlation between the illumination and reflectance is increased as the slopes of the terrain are greater than in the second and third columns corresponding to the scatterplots after the topographic correction (N-SC and SC).
In the scatterplot graphs, we can see more clearly the moderating effect of the C(i,j) parameter determined by SC (third column, Figure 6). It is possible to observe that the values of normalized reflectance of those pixels corresponding to lower illumination values increase with the topographic correction by SC in comparison to the N-SC correction. On the other hand, it is also possible to see that the values of the normalized reflectance of those pixels corresponding to higher illumination values decrease with the topographic correction by SC compared to the N-SC correction. Finally, the values of normalized reflectance of those pixels that correspond to mid-range illumination values, show a similar result, either by applying one or another correction (SC or N-SC). The moderating effect of the C(i,j) parameter increases according to the slope, i.e., the effect is much more noticeable on slopes greater than 40 • (Figure 6).
In agreement with previous works [2], it can be seen that the relationship between reflectance and illumination (cos i) is broken by the effect of the topographic corrections; being better for SC correction as the slopes increases, this fact is more perceptible in the change of the coefficient of determination (R 2 ) ( Table 2). Another important issue to note is that the number of analyzed pixels has to do directly with the determination of the C parameter value and its moderating effect. This means that, in the class of slopes with the highest number of pixels (e.g., 25-30 • ), the normalized reflectance for both corrections (SC or N-SC) is very similar, and the difference between the effects of the moderating C(i,j) parameter is less noticeable. However, it is still possible to see the difference ( Figure 6) in the trend line and the value of R 2 , whose value for N-SC correction is 0.0173 versus 0.0168 for SC correction.
This fact is also relevant since, in our work, unlike other studies, we have used all the pixels of the study area organized into classes to determine the values of the C(i,j) parameter. Thus, the individual C(i,j) parameter that most resembles the global C(i) parameter will be the one corresponding to the class that has the highest number of pixels evaluated (Figure 3b). This fact also explains the pattern of scatterplots for classes with a significant number of pixels. Table 3 shows the absolute and relative differences between mean reflectance values of pixels with slope classes greater than 0 • and the average reflectance values of pixels with slopes equal to 0 • (flat), before and after topographic corrections, including an average percentage of reflectance variation and its corresponding standard deviation for oak, pine, and mesophyll forest land cover. In agreement with previous works [2], it can be seen that the relationship between reflectance and illumination (cos i) is broken by the effect of the topographic corrections; being better for SC correction as the slopes increases, this fact is more perceptible in the change of the coefficient of determination (R 2 ) ( Table 2).  According to the mean values, the major differences are seen before the topographic corrections and specifically on the highest slopes. It can be seen how these differences are smaller after the topographic corrections, which indicates a homogenization of the original reflectance values and, thus, an improvement due to the correction. Additionally, it can be seen that the average values obtained by SC correction, compared with the values obtained by N-SC correction, show differences closer to 0 • .
In Supplementary Material Table S4, the absolute and relative reflectance differences, with respect to flat terrain for oak, pine, and mesophyll forest for all dates of study, are shown.
This pattern described for oak, pine, and mesophyll forest land cover is depicted in the corresponding graph (Figure 7a). Similar patterns can be noted (differences closer to 0 • for correction by SC), in graphs corresponding to the same land cover performed on additional dates of the study (Figure 7b,c). Table S4, the absolute and relative reflectance differences, with respect to flat terrain for oak, pine, and mesophyll forest for all dates of study, are shown.

In Supplementary Material
This pattern described for oak, pine, and mesophyll forest land cover is depicted in the corresponding graph (Figure 7a). Similar patterns can be noted (differences closer to 0° for correction by SC), in graphs corresponding to the same land cover performed on additional dates of the study (Figure 7b,c). This same task was performed for each date of the study for the bare soil land cover, and the results were not the same since, in this case, the differences closer to 0° are now corresponding to reflectance values before applying the topographic correction (Figure 7d-f). This fact, once again, suggests that the SCS + C topographic correction model is not the most suitable approach for bare soil land cover, which is consistent with previous studies [2,17].
The analysis of the differences before and after corrections in different land covers seems to confirm that the correction by SC is better, and that the SCS + C topographic correction model is the appropriate framework for forested scenes, but not for non-forest areas.

Conclusions
The incorporation of semi-empirical C parameter to the topographic correction methods reduces the over-correction that previous models had. However, the definition of different C(i,j) parameters, each generated from various linear regressions based on the slopes arranged in classes, seems to be a more accurate technique in a topographic correction because the final reflectance values are normalized in a way that is better than applying a unique C parameter determined by a single linear regression without slope differentiation.
This fact seems more evident in more shady areas with high steep slopes and can represent significant benefits in subsequent imagery processes.
The SC can be applied to any topographic correction method that includes a semi-empirical moderator determined by a linear regression in their model equations. The ranks of the slope classes can be defined according to the needs of the study and can be performed by automated unsupervised processes. The final user should choose their application depending on the specific characteristics and needs of the study, considering that processing efforts may be required. Relative reflectance and differences with respect to flat terrain, for oak, pine, and mesophyll forest and bare soil for the dates of study. This same task was performed for each date of the study for the bare soil land cover, and the results were not the same since, in this case, the differences closer to 0 • are now corresponding to reflectance values before applying the topographic correction (Figure 7d-f). This fact, once again, suggests that the SCS + C topographic correction model is not the most suitable approach for bare soil land cover, which is consistent with previous studies [2,17].
The analysis of the differences before and after corrections in different land covers seems to confirm that the correction by SC is better, and that the SCS + C topographic correction model is the appropriate framework for forested scenes, but not for non-forest areas.

Conclusions
The incorporation of semi-empirical C parameter to the topographic correction methods reduces the over-correction that previous models had. However, the definition of different C(i,j) parameters, each generated from various linear regressions based on the slopes arranged in classes, seems to be a more accurate technique in a topographic correction because the final reflectance values are normalized in a way that is better than applying a unique C parameter determined by a single linear regression without slope differentiation.
This fact seems more evident in more shady areas with high steep slopes and can represent significant benefits in subsequent imagery processes.
The SC can be applied to any topographic correction method that includes a semi-empirical moderator determined by a linear regression in their model equations. The ranks of the slope classes can be defined according to the needs of the study and can be performed by automated unsupervised processes. The final user should choose their application depending on the specific characteristics and needs of the study, considering that processing efforts may be required.

Future Research Lines
Through the development of the present work, it has been possible to improve the results of the conventional application of the topographic correction SCS + C method. However, we have also seen that it is possible to enrich our research, looking for other ways to calculate the empirical C parameter and compare its performance in the topographic correction. Thus, we believe that it is possible to suggest, as future works, the calculation of parameters based on the land covers or the performance of classification tests using other variables, such as illumination (cos i) or aspect.
In the same line of our research, and based on previous works [43], we believe that it is possible to evolve our approach if, instead of developing the correction by slope ranges, we could apply the topographic correction using continuous values of C parameters through the construction of C parameter continuous image pixel-based data.
The data found and tests implemented in the present study are not enough to conclude if the observed variations in the results obtained from the applied topographic corrections are caused by the influence of the sun geometry according to the date of image acquisition. Thus, particular research should also be developed with much more remote sensing data from specific dates.