Correction for the Impact of the Surface Characteristics on the Estimation of the Effective Emissivity at Fine Resolution in Urban Areas

Most of the methods used to retrieve land surface temperature (LST) from thermal infrared (TIR) satellite data in urban areas do not take into account the complexity of the surface. Cities are characterized by high surface roughness and one of the main constraints to estimate LST over those areas is the difficulty to define an effective emissivity for a given pixel at a given scale. When working with mixed pixels, the emissivity used to estimate the LST is an effective emissivity composed of the emissivities of each basic element constituting the pixel. In urban areas, the surface geometry has a strong impact on this effective emissivity. Its estimation from TIR satellite data must be carried out considering multiple surface reflections and diffusions within the urban canopy in order to retrieve accurate LST values. The objective of this study is then to evaluate the impact of the surface geometry within the pixel on effective emissivity estimation and to propose a method to derive an effective emissivity corrected for those effects. Emissivity can be derived at 90 m of spatial resolution from the TIR data acquired by ASTER. To evaluate the impact of the geometry at the scale of an ASTER pixel, several urban canyon configurations are designed to develop and test the correction method. The basic principle behind the method is to accurately estimate the downwelling TIR radiation received by a pixel integrating contributions from both the atmosphere and the scene inside this pixel and then derive the corrected effective emissivity from ASTER data using the TES (temperature emissivity separation) algorithm. First, the total downwelling TIR radiation is estimated from the geometric characteristics of the scene, using morphological indicators and integrating the non-isothermal behavior of the pixel thanks to 3D thermo-radiative model simulations. The validation of those estimations for each canyon configuration provides a maximum RMSE (Root Mean Square Error) value of 2.2 W·m−2. The validation performed over a district extracted from the 3D numerical model of Strasbourg (France) shows a RMSE of 2.5 W·m−2. Once the method to estimate the total downwelling TIR radiation is validated, LSE and LST maps are retrieved from an ASTER image over three districts of Strasbourg, showing that accounting for the surface geometry highlights thermal behavior differences inside districts, and that the impact of the geometry seems more influenced by building height than street width or building density.


Introduction
The growing trend towards urbanization raises a number of issues such as deteriorating air quality or rising temperatures [1,2].The adaptation of cities to climate change has become necessary and many actions are currently carried out in order to understand and improve the comfort and quality of life in urban areas.However, the effectiveness of the implemented or planned actions must be quantified to evaluate their relevance and performance.To that purpose, several studies highlighted the potential of satellite remote sensing, each image providing instantaneous snapshots of an entire urban area [3][4][5][6][7][8].Those studies mainly deal with the urban heat island phenomenon and strongly rely on the estimation of land surface temperature (LST) and emissivity (LSE) from different thermal infrared (TIR) space borne sensors.
The finest spatial resolution offered by the current TIR sensors, such as ASTER (asterweb.jpl.nasa.gov) or Landsat (landsat.usgs.gov), is about 100 m, which is considered as sufficient to represent most intraurban temperature variations for urban climate studies [9].At this scale, each pixel is a mixture of different urban elements such as buildings, streets, courtyards or squares, and its surface is then characterized by a strong roughness.This surface complexity inside the pixels is one of the main obstacles to the accurate retrieval of the LST [10][11][12][13][14][15] and is often neglected when working over urban areas.When dealing with mixed pixels, the estimation of the LST requires the use of an effective emissivity, depending on the emissivity of each component of the pixel.In cities, this effective emissivity is strongly influenced by the surface geometry and its estimation must be carried out by integrating the emission and multiple reflections of the urban canopy.It is crucial to derive a correct effective emissivity to accurately estimate LST from TIR satellite data over urban areas.Currently, LST and LSE data products are operationally derived at 90 m spatial resolution from ASTER data using the TES (temperature emissivity separation) algorithm [16].Although this method provides satisfying results in most cases [17,18], it does not take into account the complexity of the surface over urban areas while the surface geometry was identified as one of the factors of error [18].The objective of the study is then to address the issue of geometry within the pixel by investigating the impact of this geometry on the surface emissivity and temperature retrieved using TES and to propose a method to correct for those geometric effects.
The TES requires two inputs: the surface-leaving radiance measured in the five TIR bands of ASTER and the corresponding downwelling sky irradiance for each band, currently simulated as the top-of-canopy TIR radiation coming from the atmosphere using MODTRAN [19].The idea is to improve the estimations of the TES by providing better estimates of its inputs, in this case the downwelling TIR radiation taking into account the geometry within the pixel.To do so, the impact of the surface geometry is integrated into the estimation of the downwelling TIR radiation for each pixel by taking into account all contributions, namely radiation coming from the atmosphere, radiation emitted by the surrounding elements, and radiation from multiple reflections in the scene inside the pixels.This total downwelling TIR radiation can be estimated precisely using 3D thermo-radiative models, such as DART [20], TUF-3D [21], LASER/F [22], SOLENE-Microclimat [23], ENVI-met [24], etc., which simulate the energy transfers at the soil-atmosphere interface for urban environment.Those models are applied and evaluated in several studies [25][26][27][28][29][30].However, most of those models provide broadband estimations, which are then difficult to convert to narrowband estimations corresponding to the ASTER spectral bands.For those allowing narrowband simulations, they require complex parameterization datasets.Setting out from these facts, this study proposes a method to estimate the total downwelling TIR radiation narrowband based on the geometric characteristics of the surface derived from a digital surface model (DSM) combined with metric spatial resolution LST simulated with a simple 3D radiative model parameterization.This hybrid approach has the advantage to overcome the spectral dimension problem by using the LST, which can be simulated using rather simple parameterization as compared to spectral radiances, while still accounting for the pixel non-isothermal behavior.Once validated, the total downwelling TIR radiation can be provided as input to the TES to produce LSE and LST maps corrected for the impact of surface geometry.Those maps can then be compared to equivalent non-corrected maps and to the ASTER Global Emissivity Dataset v3.The results of the study are divided into two parts:

•
Validation of the estimation of the total downwelling TIR radiation based on the geometric characteristics of the surface and the 3D model simulations.

•
Production of effective LSE and LST maps using the TES accounting for the surface geometry and comparison with equivalent maps computed with the original method.
The method developed to estimate the total downwelling TIR radiation is described in Section 2, the study area is presented in Section 3, while the data and the 3D thermo-radiative models used in this study are detailed in Section 4.Then, Section 5 provides results concerning the validation of the estimated total downwelling TIR radiation along with the comparison of three LSE and LST maps derived from an ASTER image produced with and without geometric correction.Finally, Section 6 presents the discussion and conclusions.

Method
The TES algorithm simultaneously retrieves the LSE and LST from multispectral TIR satellite data based on the relation between surface-leaving radiance measured by the sensor (compensated for atmospheric absorption and path radiance) and the downwelling TIR radiation received by the surface [16].The ASTER emissivity product provided by the NASA's Jet Propulsion Laboratory is derived from the surface-leaving radiance (W•m −2 •sr −1 •µm −1 ) and the downwelling sky irradiance (W•m −2 •µm −1 ) for the five ASTER TIR bands at 90 m of spatial resolution from the AST09 product [31].The downwelling sky irradiance from AST09 is simulated using MODTRAN and corresponds to the amount of TIR radiation received at the top of the urban canopy.In urban areas, the total downwelling TIR radiation received for a given pixel is not equal to the top-of-canopy downwelling TIR radiation due to the roughness of the surface.As shown in Figure 1, several TIR radiative intra-pixel processes contribute to the total downwelling TIR radiation: (1) radiation coming from the atmosphere (R ↓ at ); (2) radiation coming from the scene, divided in two components: the emission (R ↓ em ) and multiple reflections (R ↓ mr ) from the surrounding elements.The method developed to estimate the total downwelling TIR radiation is described in Section 2, the study area is presented in Section 3, while the data and the 3D thermo-radiative models used in this study are detailed in Section 4.Then, Section 5 provides results concerning the validation of the estimated total downwelling TIR radiation along with the comparison of three LSE and LST maps derived from an ASTER image produced with and without geometric correction.Finally, section 6 presents the discussion and conclusions.

Method
The TES algorithm simultaneously retrieves the LSE and LST from multispectral TIR satellite data based on the relation between surface-leaving radiance measured by the sensor (compensated for atmospheric absorption and path radiance) and the downwelling TIR radiation received by the surface [16].The ASTER emissivity product provided by the NASA's Jet Propulsion Laboratory is derived from the surface-leaving radiance (W•m −2 •sr −1 •μm −1 ) and the downwelling sky irradiance (W•m −2 •μm −1 ) for the five ASTER TIR bands at 90 m of spatial resolution from the AST09 product [31].The downwelling sky irradiance from AST09 is simulated using MODTRAN and corresponds to the amount of TIR radiation received at the top of the urban canopy.In urban areas, the total downwelling TIR radiation received for a given pixel is not equal to the top-of-canopy downwelling TIR radiation due to the roughness of the surface.As shown in Figure 1, several TIR radiative intrapixel processes contribute to the total downwelling TIR radiation: (1) radiation coming from the atmosphere ( ↓ ); (2) radiation coming from the scene, divided in two components: the emission ( ↓ ) and multiple reflections ( ↓ ) from the surrounding elements.To integrate the impact of the surface geometry into the TES, the idea is to estimate and use the total downwelling TIR radiation ( , ↓ ) by pixel for each ASTER spectral band instead of the downwelling sky irradiance (W•m −2 •μm −1 ) from the AST09 product (Figure 2).At pixel scale, the total downwelling TIR radiation should take into account the three incident contributions as follow: To integrate the impact of the surface geometry into the TES, the idea is to estimate and use the total downwelling TIR radiation (R ↓ T,λ ) by pixel for each ASTER spectral band instead of the downwelling sky irradiance (W•m −2 •µm −1 ) from the AST09 product (Figure 2).At pixel scale, the total downwelling TIR radiation should take into account the three incident contributions as follow: Remote Sens. 2018, 10, x FOR PEER REVIEW 3 of 20 The method developed to estimate the total downwelling TIR radiation is described in Section 2, the study area is presented in Section 3, while the data and the 3D thermo-radiative models used in this study are detailed in Section 4.Then, Section 5 provides results concerning the validation of the estimated total downwelling TIR radiation along with the comparison of three LSE and LST maps derived from an ASTER image produced with and without geometric correction.Finally, section 6 presents the discussion and conclusions.

Method
The TES algorithm simultaneously retrieves the LSE and LST from multispectral TIR satellite data based on the relation between surface-leaving radiance measured by the sensor (compensated for atmospheric absorption and path radiance) and the downwelling TIR radiation received by the surface [16].The ASTER emissivity product provided by the NASA's Jet Propulsion Laboratory is derived from the surface-leaving radiance (W•m −2 •sr −1 •μm −1 ) and the downwelling sky irradiance (W•m −2 •μm −1 ) for the five ASTER TIR bands at 90 m of spatial resolution from the AST09 product [31].The downwelling sky irradiance from AST09 is simulated using MODTRAN and corresponds to the amount of TIR radiation received at the top of the urban canopy.In urban areas, the total downwelling TIR radiation received for a given pixel is not equal to the top-of-canopy downwelling TIR radiation due to the roughness of the surface.As shown in Figure 1, several TIR radiative intrapixel processes contribute to the total downwelling TIR radiation: (1) radiation coming from the atmosphere ( ↓ ); (2) radiation coming from the scene, divided in two components: the emission ( ↓ ) and multiple reflections ( ↓ ) from the surrounding elements.To integrate the impact of the surface geometry into the TES, the idea is to estimate and use the total downwelling TIR radiation ( , ↓ ) by pixel for each ASTER spectral band instead of the downwelling sky irradiance (W•m −2 •μm −1 ) from the AST09 product (Figure 2).At pixel scale, the total downwelling TIR radiation should take into account the three incident contributions as follow:  This paper proposes a method to derive R ↓ T,λ for the 5 ASTER TIR bands based on the geometric characteristics of the considered scene and simulated LST.

Atmospheric TIR Radiation
When considering a given location x in an urban canyon (e.g., the grey dot in Figure 1), the downwelling TIR radiation coming from the atmosphere in a given band (R ↓ at,λ,x ) is derived from the downwelling sky radiation at the top of the urban canopy in that band (R ↓ at,λ,H ) weighted by the portion of sky visible from that location, namely its sky view factor (SVF), and is expressed as: To estimate R ↓ at,λ at pixel level, the SVF should be representative of the considered pixel.Groleau et al. and Bernabé et al. [32,33] showed that the effective SVF (SVF T ) for a given urban area varies linearly with facade density (D F ) and established the following relationship: with: S T being the total land surface within the considered area, in this case a pixel.It corresponds to the sum of the total surface of roof elements (S R ), the total surface of facade elements (S F ) and the total surface of ground (S G ).Then, the atmospheric contribution is estimated at pixel scale for a given ASTER band as:

Scene Emitted Radiation
The contribution of the scene emission to the energy travelling below the rooftops at location x (R ↓ em,λ,x ) is computed by weighting the scene emitted TIR radiation (R em,λ,s ) by the terrain view factor, the portion of scene elements visible from location x.In this study, the terrain view factor is estimated as (1 − SVF).Furthermore, considering only the radiative exchanges below the rooftops, roofs emitted radiation is assumed as negligible as most of it is sent back to the atmosphere and very few contributes to the amount of TIR radiations remaining within the canopy.Thus, only the emission from the facades and the ground are considered.The contribution coming from the scene emission is then equal to: Using SVF T (Equation ( 3)), the scene emitted TIR radiation contribution for a given pixel is estimated as: R ↓ em,λ,s , the TIR radiation emitted by all facade and ground surfaces comprised in a pixel, is derived from the LST simulated by the 3D thermo-radiative model.The latter provides LST values for each mesh of a 3D scene fragmented at metric resolution.Then, for each mesh i, the emitted TIR radiation R ↓ em,λ,i is computed from the simulated LST and the spectral emissivity set for each element of the scene, respectively R ↓ Ts,i and ε λ,i , as follow: Then, R ↓ em,λ,s is computed as the surface weighted average of R ↓ em,λ,i emitted by each facade and ground meshes within the pixel, S i being the surface the meshes:

Multiple Reflections
The amount of TIR radiation received at location x from the multiple reflections (R ↓ mr,λ,x ) is estimated based on the capacity of the scene to reflect the TIR radiation coming from the atmosphere and the scene emission.The reflectivity of the scene is expressed as (1 − ε s,λ ) and the portion of reflected TIR radiation remaining in the urban canopy is estimated based on the terrain view factor as follow: As for the emitted scene radiation, the emission from the roofs is neglected as we only consider the radiation trapped into the urban canopy.Then, the first iteration to estimate the amount of reflected radiation for a given pixel is computed as follow: Using ), and after the several reflections: Then, the total amount of reflected TIR radiation remaining in the urban canopy can be written as: With a = (1 When n tends to infinity, a n−1 tends to 0, then:

Study Area
The method is developed and tested over the urban area of Strasbourg.Located in the east of France, this agglomeration covers 78.26 km 2 and counts slightly more than 276,000 inhabitants.The effective LSE and the LST are estimated for three districts presenting different urban configurations, the objective being to investigate the impact of the geometric correction and to quantify the errors caused by neglecting the surface geometry.The three districts, namely the historical center, the university district and the Neudorf district, are presented in Figure 3 and their morphological characteristics are summarized in Table 1.

3D Surface Geometry
The 3D geometric information used to compute the morphological characteristics of the surface by pixel is derived from the IGN (French National Institute for Geographical Information) dataset BD Topo (Figure 4) [34].The BD Topo is a simplified 3D vectorial description of the objects composing the surface with a metric resolution and scales ranging from 1:5000 to 1:50,000.This dataset covers all the geographical and administrative entities of France.It is used to derive the facade density within each pixel required for the computation of the effective SVF as well as the digital surface model (DSM) used by the 3D thermo-radiative models.

ASTER Products
The ASTER products used in this study are distributed at the NASA Land Processes Distributed Active Archive Center (LP DAAC) at the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS).

3D Surface Geometry
The 3D geometric information used to compute the morphological characteristics of the surface by pixel is derived from the IGN (French National Institute for Geographical Information) dataset BD Topo (Figure 4) [34].The BD Topo is a simplified 3D vectorial description of the objects composing the surface with a metric resolution and scales ranging from 1:5000 to 1:50,000.This dataset covers all the geographical and administrative entities of France.It is used to derive the facade density within each pixel required for the computation of the effective SVF as well as the digital surface model (DSM) used by the 3D thermo-radiative models.

3D Surface Geometry
The 3D geometric information used to compute the morphological characteristics of the surface by pixel is derived from the IGN (French National Institute for Geographical Information) dataset BD Topo (Figure 4) [34].The BD Topo is a simplified 3D vectorial description of the objects composing the surface with a metric resolution and scales ranging from 1:5000 to 1:50,000.This dataset covers all the geographical and administrative entities of France.It is used to derive the facade density within each pixel required for the computation of the effective SVF as well as the digital surface model (DSM) used by the 3D thermo-radiative models.

ASTER Products
The ASTER products used in this study are distributed at the NASA Land Processes Distributed Active Archive Center (LP DAAC) at the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS).

ASTER Products
The ASTER products used in this study are distributed at the NASA Land Processes Distributed Active Archive Center (LP DAAC) at the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS).

TES Inputs
The TES inputs consist in the surface-leaving radiance measured in the 5 TIR bands of ASTER and the corresponding downwelling sky irradiance for each band.In this study, the TIR at surface-leaving radiance provided by the AST09 product (bands 10 to 14) is used.This product was derived from an ASTER image acquired on the 28 May 2013 and also provide the sky irradiance at the top of the urban canopy, which will be used to compare to the proposed method.The validation of those products showed that, under ideal conditions and for a large homogeneous scene, the temperature and emissivity standard products are predicted to be within +1.5 K and +0.015 of correct values, respectively [17].

ASTER Global Emissivity Dataset
The ASTER GEDv3 propose an average emissivity at global scale computed using all ASTER scenes from 2000 to 2008.The product is available at 100 m, 1 km and 5 km spatial resolution.The emissivity was retrieved using the TES algorithm combined with a water vapor scaling (WVS) atmospheric correction approach [35].The validation of the product shows an agreement to within 1.5% [36].In this study, the dataset at 100 m is used to compare with the LSE maps retrieved after applying the proposed geometric correction.

3D Thermo-Radiative Models
To consider the non-isothermal behavior of the pixel, the LST used to estimate the TIR radiation emitted by the scene is simulated at mesh level using a 3D thermo-radiative model.In this study, two models are used: LASER/F [22] and SOLENE-Microclimate [23].Both models simulate 3D radiative and energy balance for each mesh of a discretized urban area based on thermodynamics principles.Those models have been validated and used in several studies [29,30,37,38].As compared to other models based on the street canyon assumption [39][40][41], the distribution of the radiative energy is performed integrating the complexity and variety of the 3D geometric configurations.The urban canopy is reproduced using 3D shapes of all the objects constituting the urban environment, as provided by the DSM; considering their geometrical and thermophysical properties.Here, the DSM is derived from the BD Topo.Only buildings are considered-no vegetation is included in the 3D surface model.Prior to the computations, the scene is fragmented in small elementary meshes.The simulations are performed for each mesh, providing results at a resolution of a few square meters.The granularity of the simulation depends of the fragmentation size, the latter influencing the computation time as well.LASER/F is used to quantify the impact of the surface geometry on the estimation of downwelling TIR radiation and to validate the correction method, both at canyon level.To extend the application of the method to the district level, simulations are performed using SOLENE-Microclimate, whose implementation allows to compute significantly larger 3D scenes as compared to LASER/F.
Both models require ground meteorological forcing data, such as air humidity and temperature, incoming radiative fluxes (top of canopy short and longwave radiations), precipitation, wind speed and orientation.These data are acquired on a regular basis in Strasbourg by fixed meteorological stations located in the city center [42].

Results
This section is divided in two parts.In the first part, the method proposed to estimate the total downwelling TIR radiation is validated at canyon level based on LASER/F simulations.In the second part, the total downwelling TIR radiation is estimated at district scale with the validated method and SOLENE-microclimat and used in the TES to derive LSE and LST maps from an ASTER image.Those maps are compared to equivalent maps produced using TOC-downwelling TIR radiation to evaluate the differences in LSE and LST obtained with the TES with and without geometric correction and to the ASTER Global Emissivity Dataset v3.

Validation of R ↓
T at Canyon Scale The validation presented in this section compares the total, atmospheric and scene downwelling TIR radiations estimated using the proposed method (hereafter mentioned as 'estimated') to the corresponding fine-grained simulated values, computed broadband by the reference 3D model LASER/F (hereafter mentioned as 'reference').First, the validation is performed for 27 different canyon configurations with three different building heights (10, 30, 50 m) and street width (5, 25, 50 m) representative of Strasbourg urban area and three different orientations (0, 45, 90 • ) (Figure 5).Four periods of 3 to 6 clear days taken in March, June, September and December are extracted from a set of meteorological data measured at the ground to force the model.For the fragmentation, the mesh minimum size is set to 15 m 2 .The albedo, emissivity, and thermal properties of the materials of the scene are detailed in Table 2.Those settings are defined according to the outputs of several experiments carried out with the model over Strasbourg over the past decade [38,43].In this case, the emissivity is kept constant for all surfaces of the scene, which is not the case in reality.However, this assumption is acceptable as a sensitivity analyses performed with LASER/F showed that varying the emissivity values from 0.85 to 0.95 lead to variations in surface temperature lower than 0.24 • C [38].The reference values are simulated for each mesh of the scene and then averaged at canyon level.The overall computation time for each period and each configuration considering the described scenes is provided in Appendix A to this paper (Table A1).

Validation of ↓ at Canyon Scale
The validation presented in this section compares the total, atmospheric and scene downwelling TIR radiations estimated using the proposed method (hereafter mentioned as 'estimated') to the corresponding fine-grained simulated values, computed broadband by the reference 3D model LASER/F (hereafter mentioned as 'reference').First, the validation is performed for 27 different canyon configurations with three different building heights (10, 30, 50 m) and street width (5, 25, 50 m) representative of Strasbourg urban area and three different orientations (0, 45, 90°) (Figure 5).Four periods of 3 to 6 clear days taken in March, June, September and December are extracted from a set of meteorological data measured at the ground to force the model.For the fragmentation, the mesh minimum size is set to 15 m 2 .The albedo, emissivity, and thermal properties of the materials of the scene are detailed in Table 2.Those settings are defined according to the outputs of several experiments carried out with the model over Strasbourg over the past decade [38,43].In this case, the emissivity is kept constant for all surfaces of the scene, which is not the case in reality.However, this assumption is acceptable as a sensitivity analyses performed with LASER/F showed that varying the emissivity values from 0.85 to 0.95 lead to variations in surface temperature lower than 0.24 °C [38].The reference values are simulated for each mesh of the scene and then averaged at canyon level.The overall computation time for each period and each configuration considering the described scenes is provided in appendix to this paper (Table A1).  Figure 6 presents a set of nine validation graphs for canyon configurations with different building and street dimensions but the same orientation.The latter was found to have no significant impact on the results.The validation of ↓ is divided into atmospheric and scene contributions, the latter comprising the scene emitted radiation and the multiple reflections.Each graph shows ↓ represented with black circles, ↓ with dark grey dots and ↓ with light grey dots along with the RMSE computed for ↓ only.A detailed validation of atmospheric and scene contributions is  Figure 6 presents a set of nine validation graphs for canyon configurations with different building and street dimensions but the same orientation.The latter was found to have no significant impact on the results.The validation of R ↓ T is divided into atmospheric and scene contributions, the latter comprising the scene emitted radiation and the multiple reflections.Each graph shows R ↓ T represented with black circles, R ↓ at with dark grey dots and R ↓ s with light grey dots along with the RMSE computed for R ↓ T only.A detailed validation of atmospheric and scene contributions is provided further.Figure 6 highlights that very accurate values of R ↓ T can be retrieved using the proposed method based on the geometry characteristics of the canyon.The largest RMSE values (1.8, 2.1 and 2.2 W•m −2 ) occur for canyons at 50-m streetwide, no matter what the height of the buildings.Both R ↓ at and R ↓ s are also correctly estimated, even if a slight overestimation of the scene contribution and underestimation of the atmospheric contribution can be observed for the canyons with buildings higher than 30 m and streets wider than 25 m.This figure shows that the atmospheric contribution decreases when street width decreases and when building height increases.With building height below 30 m, the atmospheric radiation represents the major contribution to the total downwelling TIR radiation.On the contrary, the scene contribution increases together with building height and decreases when street width increases.Nevertheless, both contributions seem to be more influenced by the height of the building than the street width.The same behavior was observed for the four time periods.

Layer composition (c (J•kg
Remote Sens. 2018, 10, x FOR PEER REVIEW 9 of 20 Both ↓ and ↓ are also correctly estimated, even if a slight overestimation of the scene contribution and underestimation of the atmospheric contribution can be observed for the canyons with buildings higher than 30 m and streets wider than 25 m.This figure shows that the atmospheric contribution decreases when street width decreases and when building height increases.With building height below 30 m, the atmospheric radiation represents the major contribution to the total downwelling TIR radiation.On the contrary, the scene contribution increases together with building height and decreases when street width increases.Nevertheless, both contributions seem to be more influenced by the height of the building than the street width.The same behavior was observed for the four time periods.Based on the different canyon configurations and the four time periods spread over the year, the impact of the geometry on the total downwelling TIR radiation is quantified.Figure 7 shows the difference between the non-corrected downwelling TIR radiation at the top of the canopy (TOC) and the computed ↓ integrating the impact of the surface geometry for each time period and configuration.Differences up to 68 W•m −2 (25% of the TOC value) are reached after correcting for the geometry.Thus, neglecting the impact of geometry could lead to errors of about 0.05 concerning the surface emissivity and 4 °C for the surface temperature.If expressed as the percentage of TOCdownwelling TIR radiation, the largest differences are obtained in March and December, when the amount of TOC-downwelling TIR radiation is the lowest.For those periods with lower atmospheric contribution, accounting for the scene contribution has a large impact on the total downwelling TIR Based on the different canyon configurations and the four time periods spread over the year, the impact of the geometry on the total downwelling TIR radiation is quantified.Figure 7 shows the difference between the non-corrected downwelling TIR radiation at the top of the canopy (TOC) and the computed R ↓ T integrating the impact of the surface geometry for each time period and configuration.Differences up to 68 W•m −2 (25% of the TOC value) are reached after correcting for the geometry.Thus, neglecting the impact of geometry could lead to errors of about 0.05 concerning the surface emissivity and 4 • C for the surface temperature.If expressed as the percentage of TOC-downwelling TIR radiation, the largest differences are obtained in March and December, when the amount of TOC-downwelling TIR radiation is the lowest.For those periods with lower atmospheric contribution, accounting for the scene contribution has a large impact on the total downwelling TIR radiation.In all cases, the largest differences are observed for canyons with tall buildings and narrow streets.As previously, those differences appear to be more linked to the building height than to the street width.No significant correlation between the height/width ratio of the canyons and those differences could be highlighted.Based on the validation presented in Figure 6, the RMSE of the estimation of ↓ is provided in Table 3.The largest RMSE (11.6 W•m −2 ) is observed for the configuration with the largest street and the highest building.The main parameter required to estimate the contribution of the atmosphere is the effective SVF ( ) for a given canyon.Errors in ↓ estimation are directly linked to the estimation of .In this method, is derived from the surface geometry based on its direct relation with the facade density as expressed in Equation ( 3).The validation of the SVF estimated from the facade density ( ) according to [32,33] using the SVF computed by LASER/F as reference, provides very accurate results (RMSE = 0.017), and shows that the proposed method is relevant to estimate for a given canyon.Further analysis pointed out that, when both building height and street width increase, the tends to be underestimated.This explains the slight underestimation of the atmospheric contribution observed in Figure 6 for canyons with buildings higher than 30 m and streets wider than Based on the validation presented in Figure 6, the RMSE of the estimation of R ↓ at is provided in Table 3.The largest RMSE (11.6 W•m −2 ) is observed for the configuration with the largest street and the highest building.The main parameter required to estimate the contribution of the atmosphere is the effective SVF (SVF T ) for a given canyon.Errors in R ↓ at estimation are directly linked to the estimation of SVF T .In this method, SVF T is derived from the surface geometry based on its direct relation with the facade density as expressed in Equation ( 3).The validation of the SVF estimated from the facade density (SVF T ) according to [32,33] using the SVF computed by LASER/F as reference, provides very accurate results (RMSE = 0.017), and shows that the proposed method is relevant to estimate SVF T for a given canyon.Further analysis pointed out that, when both building height and street width increase, the SVF T tends to be underestimated.This explains the slight underestimation of the atmospheric contribution observed in Figure 6 for canyons with buildings higher than 30 m and streets wider than 25 m.

Scene-Emitted Radiation and Multiple Reflections
Based on the validation presented in Figure 6, the RMSE of the estimation of R ↓ s is provided in Table 4.As for the atmospheric contribution, the largest RMSE (13.6 W•m −2 ) is observed for the configurations with building height above 30 m and street width larger than 25 m.As for the atmospheric contribution, the accuracy of the scene contribution estimation is directly linked to the SVF T .However, it also depends on the LST of the different elements of the canyon and the assumptions made to estimate the emitted radiation.This LST is estimated at mesh level using the 3D thermo-radiative model LASER/F in order to consider the non-isothermal behavior of the pixel.Very large temperature differences can be observed in an urban canyon, especially around solar noon.A sensitivity analysis is carried out to quantify the impact of a variation of the LST on the final R ↓ T estimation.To do so, R ↓ T is estimated using the LST simulated by the model increased or decreased by 1, 2 or 5 • C. The differences are computed as percentage of the value derive from the original LST as shown in Figure 8.As for the atmospheric contribution, the accuracy of the scene contribution estimation is directly linked to the .However, it also depends on the LST of the different elements of the canyon and the assumptions made to estimate the emitted radiation.This LST is estimated at mesh level using the 3D thermo-radiative model LASER/F in order to consider the non-isothermal behavior of the pixel.Very large temperature differences can be observed in an urban canyon, especially around solar noon.A sensitivity analysis is carried out to quantify the impact of a variation of the LST on the final ↓ estimation.To do so, ↓ is estimated using the LST simulated by the model increased or decreased by 1, 2 or 5 °C.The differences are computed as percentage of the value derive from the original LST as shown in Figure 8.The results presented in Figure 8 highlight the fact that it is very important to have realistic surface temperature for the considered scene or pixel, since a variation of 5 °C already leads to an error larger than 5% of the reference value.In a canyon, during the daytime, temperature differences up to several tens of degrees can be observed between lighted and shadowed areas.Furthermore, if well parameterized, the thermo-radiative models such as Laser/F simulate LST with an accuracy of about 2 °C [38], thus using models outputs would not introduce more than 2% of error in the final estimation.
To illustrate the necessity to consider a pixel as non-isothermal, Figure 9 shows the variability of the amount of radiation emitted by the scene according to the LST value used to estimate it (Equation [9]).The dots represent ↓ estimated according to the proposed method and the error areas are estimated taking the minimum and maximum LST values observed over the scene for each canyon configuration for each time step.The figure presents hourly estimates for four clear days The results presented in Figure 8 highlight the fact that it is very important to have realistic surface temperature for the considered scene or pixel, since a variation of 5 • C already leads to an error larger than 5% of the reference value.In a canyon, during the daytime, temperature differences up to several tens of degrees can be observed between lighted and shadowed areas.Furthermore, if well parameterized, the thermo-radiative models such as Laser/F simulate LST with an accuracy of about 2 • C [38], thus using models outputs would not introduce more than 2% of error in the final estimation.
To illustrate the necessity to consider a pixel as non-isothermal, Figure 9 shows the variability of the amount of radiation emitted by the scene according to the LST value used to estimate it (Equation [9]).The dots represent R ↓ em estimated according to the proposed method and the error areas are estimated taking the minimum and maximum LST values observed over the scene for each canyon configuration for each time step.The figure presents hourly estimates for four clear days selected in March (blue), June (green), September (red) and December (yellow).The largest errors areas are observed in June and September, when the range of LST values observed over each scene is the wider.In most of the cases, the isothermal assumption would lead to an overestimation of the contribution of the scene (maximum +196 W•m −2 ), except for canyon with low buildings and large streets.A second test two different values of hourly emitted scene radiation for each time step, both estimated using Equation ( 9), but based on two different model simulations.The first value comes from a simulation using the 3D geometry derived from the BD Topo, while the second value is simulated using the same thermodynamic characteristics but over a flat surface.For the latter, the same proportion of street, building and roof materials are kept, only the geometry is removed.The distributions of the differences observed between the two values are presented in Figure 10.As for Figure 9, the results by months (Figure 10a) confirm that larger errors are observed in June and September, when the LST range over the scene is the wider.In Figure 10b, the results also show larger discrepancies for the configuration with tall buildings and narrow streets.It also shows that considering the surface as flat would mostly lead to an overestimation of the contribution of the scene, except for canyon with low buildings and large streets.A second test compares two different values of hourly emitted scene radiation for each time step, both estimated using Equation ( 9), but based on two different model simulations.The first value comes from a simulation using the 3D geometry derived from the BD Topo, while the second value is simulated using the same thermodynamic characteristics but over a flat surface.For the latter, the same proportion of street, building and roof materials are kept, only the geometry is removed.The distributions of the differences observed between the two values are presented in Figure 10.As for Figure 9, the results by months (Figure 10a) confirm that larger errors are observed in June and September, when the LST range over the scene is the wider.In Figure 10b, the results also show larger discrepancies for the configuration with tall buildings and narrow streets.It also shows that considering the surface as flat would mostly lead to an overestimation of the contribution of the scene, except for canyon with low buildings and large streets.
distributions of the differences observed between the two values are presented in Figure 10.As for Figure 9, the results by months (Figure 10a) confirm that larger errors are observed in June and September, when the LST range over the scene is the wider.In Figure 10b, the results also show larger discrepancies for the configuration with tall buildings and narrow streets.It also shows that considering the surface as flat would mostly lead to an overestimation of the contribution of the scene, except for canyon with low buildings and large streets.

Validation Over a District
After validating the proposed method to derive R ↓ T integrating all the contributions over different canyon configurations, the method is applied and tested over a larger and more complex scene.This case deals with a real scene, extracted from the BD Topo covering the Strasbourg urban area (Figure 11).The average building height of the scene is 21.5 m and the average street width is 13.5 m.

Validation Over a District
After validating the proposed method to derive ↓ integrating all the contributions over different canyon configurations, the method is applied and tested over a larger and more complex scene.This case deals with a real scene, extracted from the BD Topo covering the Strasbourg urban area (Figure 11).The average building height of the scene is 21.5 m and the average street width is 13.5 m.As for the canyons, the estimation of over the entire scene is validated.The value computed by the 3D model LASER/F is 0.40 and the one derived from the facade density (Equation [3]) is 0.39.This estimation of the SVF based on the facade density within the pixel also provides satisfying results at district scale and can then be applied over pixels covering larger and more complex areas than just a canyon.However, the estimation of the SVF is sensitive to the size of ground area used at the border of the scene.In fact, the estimation of the facade density is directly linked to the proportion of facade, roof and ground elements.Here, the average street width is used to define the area of ground at the border of the scene.
As previously, the estimation of the total downwelling TIR radiation ( ↓ ) using the proposed method is compared to the one simulated by the 3D model for the month of September.The model setup is the same as for the validation over the canyon configurations.Figure 12 presents the validation graph for the entire scene: for ↓ (black circles), ↓ (dark dots) and ↓ (light dots) along with their respective RMSE.This validation shows that very accurate values of ↓ (RMSE = 2.5 W•m −2 ) can be retrieved using the proposed method.This graph also points out an overestimation of the atmospheric contribution and underestimation of the scene contribution, both compensating each other leading to a better estimation of the total downwelling TIR radiation.This trend is opposed to what is observed in Figure 6 for basic canyon configurations, where the scene contribution is As for the canyons, the estimation of over the entire scene is validated.The value computed by the 3D model LASER/F is 0.40 and the one derived from the facade density (Equation [3]) is 0.39.This estimation of the SVF based on the facade density within the pixel also provides satisfying results at district scale and can then be applied over pixels covering larger and more complex areas than just a canyon.However, the estimation of the SVF is sensitive to the size of ground area used at the border of the scene.In fact, the estimation of the facade density is directly linked to the proportion of facade, roof and ground elements.Here, the average street width is used to define the area of ground at the border of the scene.
As previously, the estimation of the total downwelling TIR radiation (R ↓ T ) using the proposed method is compared to the one simulated by the 3D model for the month of September.The model setup is the same as for the validation over the canyon configurations.Figure 12  respective RMSE.This validation shows that accurate values of R ↓ T (RMSE = 2.5 W•m −2 ) can be retrieved using the proposed method.This graph also points out an overestimation of the atmospheric contribution and underestimation of the scene contribution, both compensating each other leading to a better estimation of the total downwelling TIR radiation.This trend is opposed to what is observed in Figure 6 for basic canyon configurations, where the scene contribution is overestimated while the atmospheric contribution is underestimated.The observed differences may be caused by the fact that the model simulations were performed using a larger scene, to avoid boundary effects.The values corresponding to the validation area were extracted afterward, but still accounting for the effect of the surrounding buildings.Even if the value of interest to proceed with the estimation of the effective emissivity is the total radiation, this should be further investigated.Remote Sens. 2018, 10, x FOR PEER REVIEW 14 of 20

Land Surface Emissivity
Figure 13 provides two LSE maps derived using the TES algorithm applied over the ASTER image for the three selected districts.The first row of the figure shows a LSE map derived using the total downwelling TIR radiation ( ↓ ) derived from the proposed method for each ASTER pixel.The second row shows the LSE map derived using the TOC downwelling TIR, without geometric correction radiation (TOC ↓ ).The last row shows the difference between those two maps.The results are displayed for the band 12 (9.1 μm), as it is the band showing the largest differences.The statistics derived for each ASTER band are provided in Table 5. Figure 13 along with Table 5 show that integrating the impact of the surface geometry when deriving LSE from ASTER data highlights the heterogeneity inside the three districts.It also points out that, without correction, the effective emissivity is likely to be overestimated.The largest differences between corrected and non-corrected LSE is observed for the historical center, this district presenting the highest average building height (Table 1).The university district, which has a similar average building height but lower building density, is less affected by the correction.However, the Neudorf district, which presents the highest building density and the lowest average building height, is the district with the lower difference.This confirms that the building height is a key morphological parameter to account for when mapping LSE.

Land Surface Emissivity
Figure 13 provides two LSE maps derived using the TES algorithm applied over the ASTER image for the three selected districts.The first row of the figure shows a LSE map derived using the total downwelling TIR radiation (R ↓ T ) derived from the proposed method for each ASTER pixel.The second row shows the LSE map derived using the TOC downwelling TIR, without geometric correction radiation (TOC R ↓ ).The last row shows the difference between those two maps.The results are displayed for the band 12 (9.1 µm), as it is the band showing the largest differences.The statistics derived for each ASTER band are provided in Table 5. Figure 13 along with Table 5 show that integrating the impact of the surface geometry when deriving LSE from ASTER data highlights the heterogeneity inside the three districts.It also points out that, without correction, the effective emissivity is likely to be overestimated.The largest differences between corrected and non-corrected LSE is observed for the historical center, this district presenting the highest average building height (Table 1).The university district, which has a similar average building height but lower building density, is less affected by the correction.However, the Neudorf district, which presents the highest building density and the lowest average building height, is the district with the lower difference.This confirms that the building height is a key morphological parameter to account for when mapping LSE.
emissivity is likely to be overestimated.The largest differences between corrected and non-corrected LSE is observed for the historical center, this district presenting the highest average building height (Table 1).The university district, which has a similar average building height but lower building density, is less affected by the correction.However, the Neudorf district, which presents the highest building density and the lowest average building height, is the district with the lower difference.This confirms that the building height is a key morphological parameter to account for when mapping LSE.

Land Surface Temperature
As in Figure 13, Figure 14 provides two LST maps derived using the TES algorithm over the same ASTER image for the three districts of interest.First, a LST map derived using the total downwelling TIR radiation (R ↓ T ), then a LST map derived using the TOC-downwelling TIR radiation with no geometric correction (TOC R ↓ ).The third row shows the difference between those two maps and the last line displays the building footprints, extracted from the BD Topo, colored according to the building height and overlaid on the difference map.This last line aims at highlighting the spatial correlation between larger temperature differences-with and without geometric correction-and building height.Indeed, the largest differences (darker blue) are mostly located close to tall buildings (dark red).The map of the university district shows very low difference values in the southern part (bottom), which is something expected since there is nearly no building.Actually, this area is an urban park with tall trees and grass, and as such should be processed as a vegetated zone.However, as stated in Section 4, urban vegetation coverage processing is out of the scope of this study but should be considered in the future as it can lead to large discrepancies.
As for the LSE, accounting for the impact of the surface geometry highlights the LST heterogeneity inside the districts.The temperature difference statistics for each district are summarized in Table 6.The maximum observed LST difference is −1.3 • C.However, this image was taken in May and, as shows in the previous section, larger differences could be reached for warmer days and other urban configurations.
be considered in the future as it can lead to large discrepancies.
As for the LSE, accounting for the impact of the surface geometry highlights the LST heterogeneity inside the districts.The temperature difference statistics for each district are summarized in Table 6.The maximum observed LST difference is −1.3 °C.However, this image was taken in May and, as shows in the previous section, larger differences could be reached for warmer days and other urban configurations.
Figure 14.LST maps derived from ASTER for each district using total downwelling TIR radiation ( ↓ , 1st row), top-of-canopy-downwelling TIR radiation (TOC ↓ , 2nd row), the difference between both (3rd row) and difference between both + building footprints (4th row).LST maps derived from ASTER for each district using total downwelling TIR radiation (R ↓ T , 1st row), top-of-canopy-downwelling TIR radiation (TOC R ↓ , 2nd row), the difference between both (3rd row) and difference between both + building footprints (4th row).Table 6.Maximum, mean and standard deviation of the difference between LST estimated using total downwelling TIR radiation or top-of-canopy-downwelling TIR radiation for each district.

University District
Neudorf District Historical Center

Comparison with ASTER Global Emissivity Dataset v3
Finally, the LSE maps produced using the proposed method have been compared with the ASTER GED v3.This dataset is also generated using the TES algorithm but with an improved version combined with a water vapor scaling (WVS) atmospheric correction approach.However, the impact of the surface geometry is not considered during the process.Figure 15 shows the difference maps between the two datasets for band 12 for each district.The statistics derived for each ASTER band are provided in Table 7.As for Figure 13, this confirms the probability of overestimation of the emissivity over urban areas when the surface geometry is not taken into account during the estimation of the LSE.The larger mean differences in all bands are still observed for the historical center (Table 5) where the buildings are the tallest.The differences observed between the results obtained with the proposed method and the ASTER GED v3 are globally larger than the ones presented in Section 5.2.1.Those larger differences may be explained by the different version of the TES used by both methods, especially knowing the sensitivity of this algorithm to the atmospheric correction performed upstream [18,44].It would then be interesting to test the geometric correction proposed in this study on the improved TES used to generate the ASTER GED v3.

Discussion and Conclusions
Multispectral TIR satellite radiometers such as ASTER allow for simultaneous retrieval of LSE and LST over urban areas using the TES algorithm.This algorithm is operationally run to generate LSE and LST products from the surface-leaving radiance measured by ASTER and TOC-downwelling TIR radiation simulated using MODTRAN.Because of the urban surface geometrical heterogeneity, the actual amount of TIR radiation received by a pixel is not equal to the TOC-downwelling TIR radiation, which can lead to large discrepancies in LSE and LST estimates over cities.The proposed method accounts for this heterogeneity by integrating the impact of the surface geometry on LSE and LST through estimates of the total downwelling TIR radiation, including atmospheric and scene contributions, and providing it as input of the TES algorithm, in place of the TOC-downwelling TIR radiation.The total downwelling TIR radiation (R ↓ T ) is derived from surface morphological indicators, such as facade density, and 3D thermo-radiative model simulations.While the surface morphological indicators are used to estimate the SVF, the 3D model simulations are used to account for the non-isothermal behavior of the pixel when estimating the scene emitted TIR radiation.To evaluate the benefits of the proposed method, the estimations are compared to radiances simulated by a reliable 3D model, LASER/F.Various canyon configurations are presented, showing that the proposed method accurately estimates R ↓ T , accounting for the contribution from the atmosphere as well as the scene emission and multiple reflections, with a maximum RMSE of 2.2 W•m −2 .Slight overestimation of the scene contribution and underestimation of the atmospheric contribution are identified for canyons with high buildings (50 m) and a large street (50 m).Concerning the estimation of the scene contribution, results highlight the need to account for the non-isothermal behavior of the scene.The different tests show that an isothermal assumption can lead to an overestimation of the contribution of the scene, through overestimation of the scene emission up to 196 W•m -2 when taking the maximum LST observed over the scene and of 56 W•m -2 when assuming the surface as flat.The analyses of the difference between TOC-downwelling TIR radiation and the computed R ↓ T integrating the impact of the surface geometry show variations up to 68 W•m −2 (25% of the TOC value), which can lead to errors of about 0.05 concerning the surface emissivity and 4 • C for the surface temperature.The largest differences are observed when the scene contribution significantly impacts the total downwelling TIR radiation, like the test periods in March or December, or for canyons with tall buildings and narrow streets.
The same validation of R ↓ T performed over a real, larger and more complex scene, shows satisfactory results with a maximum error of 2.5 W•m −2 .In this case, slight overestimation of the atmospheric contribution and underestimation of the scene contribution are observed.Even if those errors should be further investigated, the proposed method produces very good R ↓ T estimates that can be used as inputs of the TES algorithm to retrieve effective emissivity and LST corrected for the surface geometry.Furthermore, even if no strong correlation with morphological indicators could be highlighted, the differences between TOC and R ↓ T seem to be mostly linked to the height of the building, more than street width.An analysis over three districts of the city of Strasbourg is presented to quantify the impact of this method for a real urban environment based on ASTER images.Even if this analysis does not provide information about the estimation accuracy, it illustrates that the proposed method better highlights district heterogeneity as compared to maps produced by the TES without correction.As for the canyon configurations cases, the largest differences between non-corrected and corrected values are observed in areas with tall buildings while the building density does not seem to significantly affect the results.A comparison with the ASTER GED v3 also confirms those conclusions and leads to the possibility to test the geometric correction proposed in this study on the improved TES used to generate this dataset.
The proposed method, based on geometric characteristics and 3D simulations, accurately estimates the total downwelling TIR radiation.Even if using R ↓ T to feed the TES algorithm and derive LSE and LST surely improves the estimation, it is necessary to further validate those estimations to confirm this assumption and to evaluate the real benefit of the proposed correction method.This validation requires the use of very high-resolution data, such as airborne data, to quantify the improvement brought by the correction method.Furthermore, as the effect of urban geometry seems to be mostly driven by building height, it would be interesting to extend the tests and apply the method on cities with higher average building height than Strasbourg.However, the correlation between the impact of the surface geometry on the estimations and the surface characteristics should be further investigated.It should also be noted that the proposed method should later account for vegetation structure and properties to better characterize large urban parks.Considering the importance of accounting for the actual geometrical heterogeneity, as demonstrated in this study, and the availability of 3D datasets for more and more cities worldwide, the proposed method can be implemented to provide higher quality LSE and LST products with other existing or forthcoming multispectral TIR sensors.

Figure 1 .
Figure 1.Scheme of the different radiative fluxes contributing to the total downwelling TIR radiation: atmospheric radiation (thick arrow), scene radiations divided between emission (thin arrow) and multiple reflections (dash arrow).

Figure 2 .Figure 1 .
Figure 2. Processing steps to estimate corrected and non-corrected land surface temperature ( ) and emissivity ( ) from ASTER data (TOC = top of canopy, DSM = digital surface model)

Figure 1 .
Figure 1.Scheme of the different radiative fluxes contributing to the total downwelling TIR radiation: atmospheric radiation (thick arrow), scene radiations divided between emission (thin arrow) and multiple reflections (dash arrow).

Figure 2 .
Figure 2. Processing steps to estimate corrected and non-corrected land surface temperature ( ) and emissivity ( ) from ASTER data (TOC = top of canopy, DSM = digital surface model)

Figure 2 .
Figure 2. Processing steps to estimate corrected and non-corrected land surface temperature (T s ) and emissivity (ε) from ASTER data (TOC = top of canopy, DSM = digital surface model).

Figure 3 .
Figure 3. Strasbourg urban area (top) and the three districts used for land surface emissivity and temperature estimation (bottom).

Figure 4 .
Figure 4. 3D datasets of the three districts derived from the BD Topo (IGN).

Figure 3 .
Figure 3. Strasbourg urban area (top) and the three districts used for land surface emissivity and temperature estimation (bottom).

20 Figure 3 .
Figure 3. Strasbourg urban area (top) and the three districts used for land surface emissivity and temperature estimation (bottom).

Figure 4 .
Figure 4. 3D datasets of the three districts derived from the BD Topo (IGN).

Figure 4 .
Figure 4. 3D datasets of the three districts derived from the BD Topo (IGN).

Figure 5 .Table 2 .
Figure 5.The different canyon configurations used to set and validate ↓ estimation.

Figure 6 .
Figure 6.Comparison of estimated and reference total (black circles), atmospheric (dark grey dots) and scene (light grey dots) downwelling TIR radiation for nine canyon configurations (orientation 0°).The Root Mean Square Error (RMSE) is provided for the estimation of the total downwelling TIR radiation.

Figure 6 .
Figure 6.Comparison of estimated and reference total (black circles), atmospheric (dark grey dots) and scene (light grey dots) downwelling TIR radiation for nine canyon configurations (orientation 0 • ).The Root Mean Square Error (RMSE) is provided for the estimation of the total downwelling TIR radiation.

Figure 8 .
Figure 8. Differences in ↓ with LST varying from −5 to +5 °C around the reference LST value, expressed as percentage of reference value

Figure 8 .
Figure 8. Differences in R ↓ T with LST varying from −5 to +5 • C around the reference LST value, expressed as percentage of reference value.

20 Figure 9 .
Figure 9. Hourly estimates of ↓ for four clear days: one in March (blue), one in June (green), one in September (red) and one in December (yellow) and for the nine canyon configurations.The dots represent ↓ estimated using the proposed method and the error areas are computed using the minimum and maximum LST values observed over the scene.

Figure 9 .
Figure 9. Hourly estimates of R ↓ em for four clear days: one in March (blue), one in June (green), one in September (red) and one in December (yellow) and for the nine canyon configurations.The dots represent R ↓ em estimated using the proposed method and the error areas are computed using the minimum and maximum LST values observed over the scene.

Figure 10 .
Figure 10.Distribution of the differences between hourly scene emitted radiation estimated integrating the 3D geometry and assuming a flat surface, by month (a) and by canyon configuration (b).The red line being the median value, the upper and lower limits of the box the upper and lower quartile respectively and the whiskers representing the minimum and maximum values.

Figure 10 .
Figure 10.Distribution of the differences between hourly scene emitted radiation estimated integrating the 3D geometry and assuming a flat surface, by month (a) and by canyon configuration (b).The red line being the median value, the upper and lower limits of the box the upper and lower quartile respectively and the whiskers representing the minimum and maximum values.

Figure 11 .
Figure 11.Google maps image and 3D geometry of the scene extracted from the BD Topo to extend the validation tests.

Figure 11 .
Figure 11.Google maps image and 3D geometry of the scene extracted from the BD Topo to extend the validation tests.
presents the validation graph for the entire scene: for R ↓ T (black circles), R ↓ at (dark dots) and R ↓ s (light dots) along with their Remote Sens. 2018, 10, 746 14 of 21

Figure 12 .
Figure 12.Comparison of estimated and reference total (black circles), atmospheric (dark grey dots) and scene (light grey dots) downwelling TIR radiation for a whole district.

Figure 12 .
Figure 12.Comparison of estimated and reference total (black circles), atmospheric (dark grey dots) and scene (light grey dots) downwelling TIR radiation for a whole district.

Figure 13 .
Figure 13.Emissivity maps derived from ASTER band 12 for each district using total downwelling TIR radiation ( ↓ , 1st line), the top of canopy downwelling TIR radiation (TOC ↓ , 2nd line) and the difference between both (3rd line).

Figure 13 .
Figure 13.Emissivity maps derived from ASTER band 12 for each district using total downwelling TIR radiation (R ↓ T , 1st line), the top of canopy downwelling TIR radiation (TOC R ↓ , 2nd line) and the difference between both (3rd line).

Figure 14 .
Figure 14.LST maps derived from ASTER for each district using total downwelling TIR radiation (R

Table 1 .
Morphological characteristics of each district.

Table 3 .
RMSE of the estimation of↓ expressed in W•m −2 .

Table 3 .
RMSE of the estimation of R

Table 4 .
RMSE of the estimation of R

Table 4 .
RMSE of the estimation of↓ expressed in W•m −2 .

Table 5 .
Maximum, mean and standard deviation of the difference between emissivity estimated using total downwelling TIR radiation or top of canopy downwelling TIR radiation for each district (highest values in bold).