Performance Evaluation of Urban Canopy Parameters Derived from VHR Optical Stereo Data †

: Urban canopy parameters (UCPs) are parameters which are utilized to define the thermal, radiative, and roughness properties of urban areas, which have a significant impact on the urban microclimate. The rapidly growing urbanization, especially in developing regions, leads to the modification of urban geometry, which calls for the characterization of UCPs in the countries of such regions to account for high population pressure, heterogeneous urban environments


Introduction
Urban canopy parameters (UCPs) define the thermal, radiative, and roughness properties of urban areas and impact the urban microclimate significantly.UCPs, like street height or type of buildings, etc., strongly affect the urban microclimate [1].The availability of UCPs is very important for urban climate studies, as recent urban canopy models (UCMs) ingest detailed UCPs characterizing urban heterogeneity to obtain quantified information on the urban climate.Numerical weather simulations with UCMs have also indicated that the availability of accurate information on UCPs provide improved simulation results for understanding urban heat island (UHI), air pollution, and climate change scenarios [2][3][4][5][6][7][8].The remote sensing datasets of airborne LiDAR and aerial photography have been widely used in developed countries to obtain information on UCPs [9].However, due to the Environ.Sci.Proc.2024, 29, 62 2 of 7 limited availability of airborne LiDAR and aerial photography in developing regions, ref. [10] has employed very-high-resolution satellite (VHRS) optical stereo datasets for the retrieval of building heights and UCPs in the complex urban environment of Delhi due to their cost effectiveness and easy availability [10][11][12].Although the study evaluated key input parameters, such as building height, building surface area, and land use and land cover (LULC), and found acceptable accuracies for the generated UCPs, it was felt that UCPs generated from VHRS optical stereo data needed to be evaluated against field-based benchmark data.Hence, this study chooses a representative urban area to evaluate the accuracies of generated UCPs derived from VHRS optical stereo data against reference data collected from the ground.

Study Area and Datasets
Bhubaneswar has been selected for this study because of the geographic location of the city (Figure 1).It is located near the Bay of Bengal, making it vulnerable to intense cyclones, precipitation, and flood events.The city has a warm and humid climate.The city is a planned urban climate; evaluation of the study area involves fewer constraints than other cities with very complex urban structures.Ward number 27 (also known as Nayapalli) of the Bhubaneswar Municipal Corporation has been selected for ground data collection to validate UCPs derived from satellite images.The criteria for selecting the area were based on heterogeneity in building heights and variation in LULC categories.
have been widely used in developed countries to obtain information on UCPs [9].However, due to the limited availability of airborne LiDAR and aerial photography in developing regions, ref. [10] has employed very-high-resolution satellite (VHRS) optical stereo datasets for the retrieval of building heights and UCPs in the complex urban environment of Delhi due to their cost effectiveness and easy availability [10][11][12].Although the study evaluated key input parameters, such as building height, building surface area, and land use and land cover (LULC), and found acceptable accuracies for the generated UCPs, it was felt that UCPs generated from VHRS optical stereo data needed to be evaluated against field-based benchmark data.Hence, this study chooses a representative urban area to evaluate the accuracies of generated UCPs derived from VHRS optical stereo data against reference data collected from the ground.

Study Area and Datasets
Bhubaneswar has been selected for this study because of the geographic location of the city (Figure 1).It is located near the Bay of Bengal, making it vulnerable to intense cyclones, precipitation, and flood events.The city has a warm and humid climate.The city is a planned urban climate; evaluation of the study area involves fewer constraints than other cities with very complex urban structures.Ward number 27 (also known as Nayapalli) of the Bhubaneswar Municipal Corporation has been selected for ground data collection to validate UCPs derived from satellite images.The criteria for selecting the area were based on heterogeneity in building heights and variation in LULC categories.Material used for the study were obtained from field surveys undertaken to collect building heights, hemispherical photographs, building footprints, and reference points for land use and land cover information.Around 954 building height samples and footprints were collected using the Leica Distometer Instrument for the creation of a threedimensional (3D) reference geodatabases for validation.The acquired remote sensing data for the study include stereo pairs from Pleiades.Pleiades is a French satellite and has a resolution of 0.5 m.Pleiades stereo pairs with both multi-spectral and panchromatic sensors are used for this study.

Methodology
This study was performed in a number of stages to achieve the research aim (Figure 2).The first step involved the generation of a digital surface model (DSM) from VHRS optical stereo datasets using ground control points (GCPs), which were obtained from the differential ground positioning system (DGPS survey).A digital terrain model (DTM) was generated by applying morphological filters on the generated DSM.By utilizing a DSM and a DTM of the optical stereo pair, a normalized digital surface model (nDSM) was derived by subtracting the DTM from the DSM.Accuracy assessments of all generated DSMs, DTMs, and nDSMs were performed using GCPs and field data of building heights.Material used for the study were obtained from field surveys undertaken to collect building heights, hemispherical photographs, building footprints, and reference points for land use and land cover information.Around 954 building height samples and footprints were collected using the Leica Distometer Instrument for the creation of a three-dimensional (3D) reference geodatabases for validation.The acquired remote sensing data for the study include stereo pairs from Pleiades.Pleiades is a French satellite and has a resolution of 0.5 m.Pleiades stereo pairs with both multi-spectral and panchromatic sensors are used for this study.

Methodology
This study was performed in a number of stages to achieve the research aim (Figure 2).The first step involved the generation of a digital surface model (DSM) from VHRS optical stereo datasets using ground control points (GCPs), which were obtained from the differential ground positioning system (DGPS survey).A digital terrain model (DTM) was generated by applying morphological filters on the generated DSM.By utilizing a DSM and a DTM of the optical stereo pair, a normalized digital surface model (nDSM) was derived by subtracting the DTM from the DSM.Accuracy assessments of all generated DSMs, DTMs, and nDSMs were performed using GCPs and field data of building heights.The second step involved the generation of orthorectified images using multi-spectral Pleiades images.
The second step involved the generation of orthorectified images using mu Pleiades images.The generated orthorectified multi-spectral Pleiades images were further obtain the distribution of LULC, which was followed by an accuracy assess generated classes.Supervised classifications of MX orthorectified images hav ried out for six classes: built-up, vegetation, agricultural land, water body, ba river bed.The third step of the study involved the computation of UCPs.The required for the generation of UCPs, such as building footprints and buildi were derived from the nDSM.The most widely used UCPs for climate studie calculated for the selected study area using VHRS optical stereo data.The fou volved a 3D reference geodatabase creation for the study area (ward number was obtained through an extensive field survey.The final step consisted of th tive evaluation of computed UCPs from a 3D reference geodatabase, with UCP from VHRS optical stereo data at a grid resolution of 30 m.The statistical ev UCPs was carried out by computing the bias, mean absolute error (MAE), square error (RMSE), and correlation coefficient (r).

Results
The VHRS DSM, generated at 1 m spatial resolution through photogram cessing, was found to have an RMSE of 3.08 m, while the DTM, after the filt DSM, had an RMSE of 2.8 m.The nDSM obtained after the subtraction of the the DSM displayed an RMSE of 0.38 m with respect to height observations col the ground (Figure 3a-c).Land use and land cover characteristics are very im urban climate studies.Accuracy assessment has been performed on classified selecting sample points using a stratified random sampling method.All the c classified with an overall accuracy of 91.79%, which is more than the recomme racy of 85% [13], and with an overall kappa index of 0.9042 (Figure 3d).The generated orthorectified multi-spectral Pleiades images were further classified to obtain the distribution of LULC, which was followed by an accuracy assessment of the generated classes.Supervised classifications of MX orthorectified images have been carried out for six classes: built-up, vegetation, agricultural land, water body, bare soil, and river bed.The third step of the study involved the computation of UCPs.The input data required for the generation of UCPs, such as building footprints and building heights, were derived from the nDSM.The most widely used UCPs for climate studies [13] were calculated for the selected study area using VHRS optical stereo data.The fourth step involved a 3D reference geodatabase creation for the study area (ward number 27), which was obtained through an extensive field survey.The final step consisted of the comparative evaluation of computed UCPs from a 3D reference geodatabase, with UCPs generated from VHRS optical stereo data at a grid resolution of 30 m.The statistical evaluation of UCPs was carried out by computing the bias, mean absolute error (MAE), root mean square error (RMSE), and correlation coefficient (r).

Results
The VHRS DSM, generated at 1 m spatial resolution through photogrammetric processing, was found to have an RMSE of 3.08 m, while the DTM, after the filtering of the DSM, had an RMSE of 2.8 m.The nDSM obtained after the subtraction of the DTM from the DSM displayed an RMSE of 0.38 m with respect to height observations collected from the ground (Figure 3a-c).Land use and land cover characteristics are very important for urban climate studies.Accuracy assessment has been performed on classified images by selecting sample points using a stratified random sampling method.All the classes were classified with an overall accuracy of 91.79%, which is more than the recommended accuracy of 85% [13], and with an overall kappa index of 0.9042 (Figure 3d).

Building Footprints and Building Height
Since no authentic 2D and 3D GIS databases were available for study the area, a 3D geodatabase (consisting of 954 buildings) for the study area was prepared after the collection of building footprints digitized from VHRS orthoimages from Pleiades (0.5 m spatial resolution) and the collection of building height data using the Leica Distometer Instrument, as explained in earlier sections (Figure 4a,b).Further, this dataset was utilized to assess the accuracy of the building footprints and building heights obtained from VHRS optical stereo data.

Building Footprints and Building Height
Since no authentic 2D and 3D GIS databases were available for study the area, a 3D geodatabase (consisting of 954 buildings) for the study area was prepared after the collection of building footprints digitized from VHRS orthoimages from Pleiades (0.5 m spatial resolution) and the collection of building height data using the Leica Distometer Instrument, as explained in earlier sections (Figure 4a,b).Further, this dataset was utilized to assess the accuracy of the building footprints and building heights obtained from VHRS optical stereo data.

Building Footprints and Building Height
Since no authentic 2D and 3D GIS databases were available for study the area, a 3D geodatabase (consisting of 954 buildings) for the study area was prepared after the collection of building footprints digitized from VHRS orthoimages from Pleiades (0.5 m spatial resolution) and the collection of building height data using the Leica Distometer Instrument, as explained in earlier sections (Figure 4a,b).Further, this dataset was utilized to assess the accuracy of the building footprints and building heights obtained from VHRS optical stereo data.Building footprints and building height were also retrieved from the VHRS nDSM (generated in the previous section) by following the methodology adopted by [10].The completeness analysis of the retrieved building footprints shows that 82.72% of buildings were correctly extracted from VHRS optical stereo data.A false positive rate of 17.27% shows that a much smaller percentage of buildings had been incorrectly classified as other features.Similarly, building heights retrieved from satellite data show a positive bias of 0.09 m.The MAE is 0.57 m, while the RMSE is 0.70 m.A correlation coefficient of 0.99 has been observed between the reference and retrieved building heights.

Urban Canopy Parameters
By utilizing the extracted input parameters-building footprints, building heights, and land use land cover-a set of UCPs has been derived at a spatial resolution of 30 m by following the methodology developed by [10].All twelve (12) UCPs, as shown in Figure 5a-l-mean building height (MBH) and standard deviation of building height (SDBH), building surface fraction (BSF), pervious surface fraction (PSF), impervious surface fraction (ISF), frontal area index (FAI), building-surface-to-plan-area ratio (BSPAR), heightto-width ratio (H/W ratio), complete aspect ratio (CAR), zero-plane displacement height (Z d ), roughness length (RL), and plan area density (PAD-were derived pixel by pixel.The computed UCPs from satellite data were compared against reference-data-derived UCPs using statistical measures such as bias, MAE, RMSE, and correlation coefficients, as shown in Tables 1 and 2. It can be seen that most of the satellite-data-derived UCPs have shown high correlation coefficients (above 0.7) with reference-data-derived UCPs.

Conclusions
The study has presented the evaluation of retrieved UCPs from VHRS optical stereo data for a part of Bhubaneshwar city, India, against ground-based field observations and a ground-based reference 3D geodatabase.The assessment of building footprints and building heights retrieved from VHRS optical stereo data has revealed the high performance of the Pleiades dataset in the extraction of building footprints, with a true positive rate of 82.7%, and building heights with an RMSE of 0.7 m and a correlation coefficient of 0.99.Further, performance evaluation of the most widely used UCPs for climate studies shows that at 30 m resolution, all UCPs derived from VHRS optical stereo data have demonstrated very high accuracy, with RMSE values beneath 2 and a high correlation coefficient of more than 0.7.The generated UCPs can be utilized for ventilation studies, generating urban climatic maps and local climate zonation, climatic modeling, numerical weather prediction over urban areas, shading analysis, urban heat island studies, and much more.

Figure 2 .
Figure 2. Methodology of the study.

Figure 2 .
Figure 2. Methodology of the study.

Table 1 .
Accuracy evaluation of satellite-derived UCPs with reference UCPs.

Table 2 .
Accuracy evaluation of satellite-derived UCPs with reference UCPs.