Quantifying spatiotemporal changes of the urban impervious surface of Dhaka District using Remote sensing Technology

Dhaka, the capital of Bangladesh, is one of the world's fastest-growing cities where imperviousness expanding in tandem. Therefore, accurate estimation of impervious surfaces is essential for urban planning and management. This paper attempts to quantify the changes of urban impervious surfaces in Dhaka district from 1990 to 2020 using remote sensing technology. Satellite images of 1990, 1995, 2000, 2005, 2010, 2015, and 2020 have been taken from the Landsat TM, ETM+, OLI sensor. Unsupervised classification with k-means clustering and three different RS indices NDVI, NDBI, and BUI was used to delineate the actual impervious area of Dhaka city. This study reveals that due to urbanization a net increase of 67.30 sq. miles impervious area is added to the existing amount over the study period. In 2020 total 300.749 sq. miles which contain 51.02% of the total land were occupied by impervious surfaces compared to the 233.446 sq. miles in 1990. Instantaneously taking appropriate strategies is crucial for sustainable urban growth.  

Remote sensing of impervious surface area and its interaction with land surface temperature variability in Pretoria, South Africa

Includes summary for chapter 1-5Pretoria, City of Tshwane (COT), Gauteng Province, South Africa is one of the cities that continues to experience rapid urban sprawl as a result of population growth and various land use, leading to the change of natural vegetation lands into impervious surface area (ISA). These are associated with transportation (paved roads, streets, highways, parking lots and sidewalks) and cemented buildings and rooftops, made of completely or partly impermeable artificial materials (e.g., asphalt, concrete, and brick). These landscapes influence the micro-climate (e.g., land surface temperature, LST) of Pretoria City as evidenced by the recent heat waves characterized by high temperature. Therefore, understanding ISA changes will provide information for city planning and environmental management. Conventionally, deriving ISA information has been dependent on field surveys and manual digitizing from hard copy maps, which is laborious and time-consuming. Remote sensing provides an avenue for deriving spatially explicit and timely ISA information. Numerous methods have been developed to estimate and retrieve ISA and LST from satellite imagery. There are limited studies focusing on the extraction of ISA and its relationship with LST variability across major cities in Africa. The objectives of the study were: (i) to explore suitable spectral indices to improve the delineation of built-up impervious surface areas from very high resolution multispectral data (e.g., WorldView-2), (ii) to examine exposed rooftop impervious surface area based on different colours, and their interplay with surface temperature variability, (iii) to determine if the spatio-temporal built-up ISA distribution pattern in relation to elevation influences urban heat island (UHI) extent using an optimal analytical scale and (iv) to assess the spatio-temporal change characteristics of ISA expansion using the corresponding surface temperature (LST) at selected administrative subplace units (i.e., local region scale). The study objectives were investigated using remote sensing data such as WorldView-2 (a very high-resolution multispectral sensor), medium resolution Landsat-5 Thematic Mapper (TM) and Landsat-8 OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) at multiple scales. The ISA mapping methods used in this study can be grouped into two major categories: (i) the classification-based approach consisting of an object-based multi-class classification with overall accuracy ~90.4% and a multitemporal pixel-based binary classification. The latter yielded an area under the receiver operating characteristic curve (AUROC) = 0.8572 for 1995, AUROC = 0.8709 for 2005, AUROC = 0.8949 for 2015. (ii) the spectral index-based approach such as a new built-up extraction index (NBEI) derived in this study which yielded a high AUROC = ~0.82 compared to Built-up Area Index (BAI) (AUROC = ~0.73), Built-up spectral index (BSI) (AUROC = ~0.78), Red edge / Green Index (RGI) (AUROC = ~0.71) and WorldView-Built-up Index (WV-BI) (AUROC = ~0.67). The multitemporal built-up Index (BUI) also estimated with AUROC = 0.8487 for 1993, AUROC = 0.8302 for 2003, AUROC = 0.8790 for 2013. This indicates that all these methods employed, mapped ISA with high predictive accuracy from remote sensing data. Furthermore, the single-channel algorithm (SCA) was employed to retrieve LST from the thermal infrared (TIR) band of the Landsat images. The LST overall retrieval error for the entire study generally was quite low (overall root mean square RMSE ≤ ~1.48OC), which signifies that the Landsat TIR used provided good results for further analysis. In conclusion, the study showed the potential of multispectral remote sensing data to quantify ISA and evaluate its interaction with surface temperature variability despite the complex urban landscape in Pretoria. Also, using impervious surface LST as a complementary metric in this research helped to reveal urban heat island distribution and improve understanding of the spatio-temporal developing trend of urban expansion at a local spatial scale.Rapid urbanization because of population growth has led to the conversion of natural lands into large man-made landscapes which affects the micro-climate. Rooftop reflectivity, material, colour, slope, height, aspect, elevation are factors that potentially contribute to temperature variability. Therefore, strategically designed rooftop impervious surfaces have the potential to translate into significant energy, long-term cost savings, and health benefits. In this experimental study, we used the semi-automated Environment for Visualizing Images (ENVI) Feature Extraction that uses an object-based image analysis approach to classify rooftop based on colours from WorldView-2 (WV-2) image with overall accuracy ~90.4% and kappa coefficient ~0.87 respectively. The daytime retrieved surface temperatures were derived from 15m pan-sharpened Landsat 8 TIRS with a range of ~14.6OC to ~65OC (retrieval error = 0.38OC) for the same month covering Lynwood Ridge a residential area in Pretoria. Thereafter, the relationship between the rooftops and surface temperature (LST) were examined using multivariate statistical analysis. The results of this research reveal that the interaction between the applicable rooftop explanatory features (i.e., reflectance, texture measures and topographical properties) can explain over 22.10% of the variation in daytime rooftop surface temperatures. Furthermore, analysis of spatial distribution between mean daytime surface temperature and the residential rooftop indicated that the red, brown and green roof surfaces show lower LST values due to high reflectivity, high emissivity and low heat capacity during the daytime. The study concludes that in any study related to the spatial distribution of rooftop impervious surface area surface temperature, effect of various explanatory variables must be considered. The results of this experimental study serve as a useful approach for further application in urban planning and sustainable development.Evaluating changes in built-up impervious surface area (ISA) to understand the urban heat island (UHI) extent is valuable for governments in major cities in developing countries experiencing rapid urbanization and industrialization. This work aims at assessing built-up ISA spatio-temporal and influence on land surface temperature (LST) variability in the context of urban sprawl. Landsat-5 Thematic Mapper (TM) and Landsat-8 OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) were used to quantify ISA using built-up Index (BUI) and spatio-temporal dynamics from 1993-2013. Thereafter using a suitable analytical sampling scale that represents the estimated ISA-LST, we examined its distribution in relation to elevation using the Shuttle Radar Topography Mission (SRTM) and also create Getis-Ord Gi* statistics hotspot maps to display the UHI extent. The BUI ISA extraction results show a high predictive accuracy with area under the receiver operating characteristic curve, AUROC = 0.8487 for 1993, AUROC = 0.8302 for 2003, AUROC = 0.8790 for 2013. The ISA spatio-temporal changes within ten years interval time frame results revealed a 14% total growth rate during the study year. Based on a suitable analytical scale (90x90) for the hexagon polygon grid, the majority of ISA distribution across the years was at an elevation range of between >1200m – 1600m. Also, Getis-Ord Gi* statistics hotspot maps revealed that hotspot regions expanded through time with a total growth rate of 19% and coldspot regions decreased by 3%. Our findings can represent useful information for policymakers by providing a scientific basis for sustainable urban planning and management.Over the years, rapid urban growth has led to the conversion of natural lands into large man-made landscapes due to enhanced political and economic growth. This study assessed the spatio-temporal change characteristics of impervious surface area (ISA) expansion using its surface temperature (LST) at selected administrative subplace units (i.e., local region scale). ISA was estimated for 1995, 2005 and 2015 from Landsat-5 Thematic Mapper (TM) and Landsat-8 OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) images using a Random Forest (RF) algorithm. The spatio-temporal trends of ISA were assessed using an optimal analytical scale to aggregate ISA LST coupled with weighted standard deviational ellipse (SDE) method. The ISA was quantified with high predictive accuracy (i.e., AUROC = 0.8572 for 1995, AUROC = 0.8709 for 2005, AUROC = 0.8949 for 2015) using RF classifier. More than 70% of the selected administrative subplaces in Pretoria experienced an increase in growth rate (415.59%) between 1995 and 2015. LST computations from the Landsat TIRS bands yielded good results (RMSE = ~1.44OC, 1.40OC, ~0.86OC) for 1995, 2005 and 2015 respectively. Based on the hexagon polygon grid (90x90), the aggregated ISA surface temperature weighted SDE analysis results indicated ISA expansion in different directions at the selected administrative subplace units. Our findings can represent useful information for policymakers in evaluating urban development trends in Pretoria, City of Tshwane (COT).Globally, the unprecedented increase in population in many cities has led to rapid changes in urban landscape, which requires timely assessments and monitoring. Accurate determination of built-up information is vital for urban planning and environmental management. Often, the determination of the built-up area information has been dependent on field surveys, which is laborious and time-consuming. Remote sensing data is the only option for deriving spatially explicit and timely built-up area information. There are few spectral indices for built-up areas and often not accurate as they are specific to impervious material, age, colour, and thickness, especially using higher resolution images. The objective of this study is to test the utility of a new built-up extraction index (NBEI) using WorldView-2 to improve built-up material mapping irrespective of material type, age and colour. The new index was derived from spectral bands such as Green, Red edge, NIR1 and NIR2 bands that profoundly explain the variation in built-up areas on WorldView-2 image (WV-2). The result showed that NBEI improves the extraction of built-up areas with high accuracy (area under the receiver operating characteristic curve, AUROC = ~0.82) compared to the existing indices such as Built-up Area Index (BAI) (AUROC = ~0.73), Built-up spectral index (BSI) (AUROC = ~0.78 ), Red edge / Green Index (RGI) (AUROC = ~0.71) and WorldView-Built-up Index (WV-BI) (AUROC = ~0.67). The study demonstrated that the new built-up index could extract built-up areas using high-resolution images. The performance of NBEI could be attributed to the fact that it is not material specific, and would be necessary for urban area mapping.Environmental SciencesD. Phil. (Environmental Sciences

Assessment of spatio-temporal direction of impervious surface area surface temperature in Pretoria, South Africa

Over the years, rapid urban growth has led to the conversion of natural lands into large man-made landscapes due to enhanced political and economic growth. This study assessed the spatio-temporal change characteristics of impervious surface area (ISA) expansion using its surface temperature (LST) at selected administrative subplace units (i.e., local region scale). ISA was estimated for 1995, 2005 and 2015 from Landsat-5 Thematic Mapper (TM) and Landsat-8 OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) images using a Random Forest (RF) algorithm. The spatio-temporal trends of ISA were assessed using an optimal analytical scale to aggregate ISA LST coupled with weighted standard deviational ellipse (SDE) method. The ISA was quantified with high predictive accuracy (i.e., AUROC = 0.8572 for 1995, AUROC = 0.8709 for 2005, AUROC = 0.8949 for 2015) using RF classifier. More than 70% of the selected administrative subplaces in Pretoria experienced an increase in growth rate (415.59%) between 1995 and 2015. LST computations from the Landsat TIRS bands yielded good results (RMSE = ∼1.44OC, 1.40OC, ∼0.86OC) for 1995, 2005 and 2015 respectively. Based on the hexagon polygon grid (90x90), the aggregated ISA surface temperature weighted SDE analysis results indicated ISA expansion in different directions at the selected administrative subplace units. Our findings can represent useful information for policymakers in evaluating urban development trends in Pretoria, City of Tshwane (COT).The University of South Africa Student Funding Directorate (UNISA, DSF) and GeoTerraImage (Pty) Ltd.https://www.tandfonline.com/loi/tgei20hj2023Geography, Geoinformatics and MeteorologyPlant Production and Soil Scienc

Review on Active and Passive Remote Sensing Techniques for Road Extraction

Digital maps of road networks are a vital part of digital cities and intelligent transportation. In this paper, we provide a comprehensive review on road extraction based on various remote sensing data sources, including high-resolution images, hyperspectral images, synthetic aperture radar images, and light detection and ranging. This review is divided into three parts. Part 1 provides an overview of the existing data acquisition techniques for road extraction, including data acquisition methods, typical sensors, application status, and prospects. Part 2 underlines the main road extraction methods based on four data sources. In this section, road extraction methods based on different data sources are described and analysed in detail. Part 3 presents the combined application of multisource data for road extraction. Evidently, different data acquisition techniques have unique advantages, and the combination of multiple sources can improve the accuracy of road extraction. The main aim of this review is to provide a comprehensive reference for research on existing road extraction technologies.Peer reviewe

Multisource Remote Sensing based Impervious Surface Mapping

Impervious surface (IS) not only serves as a key indicator of urbanization, but also affects the micro-ecosystem. Therefore, it is essential to monitor IS distribution timely and accurately. Remote sensing is an effective approach as it can provide straightforward and consistent information over large area with low cost. This thesis integrates multi-source remote sensing data to interpretate urban patterns and provide more reliable IS mapping results. Registration of optical daytime and nighttime lights (NTL) data is developed in the first contribution. An impervious surface based optical-to-NTL image registration algorithm with iterative blooming effect reduction (IS_iBER) algorithm is proposed. This coarse-to-fine procedure investigates the correlation between optical and NTL features. The iterative registration and blooming effect reduction method obtains precise matching results and reduce the spatial extension of NTL. Considering the spatial transitional nature of urban-rural fringes (URF) areas, the second study proposed approach for URF delineation, namely optical and nighttime lights (NTL) data based multi-scale URF (msON_URF).The landscape heterogeneity and development vitality derived from optical and NTL features are analyzed at a series of scales to illustrate the urban-URF-rural pattern. Results illustrate that msON_URF is effective and practical for not only concentric, but also polycentric urban patterns. The third study proposes a nighttime light adjusted impervious surface index (NAISI) to detect IS area. Parallel to baseline subtraction approaches, NAISI takes advantage of features, rather than spectral band information to map IS. NAISI makes the most of independence between NTL-ISS and pervious surface to address the high spectral similarity between IS and bare soil in optical image. An optical and NTL based spectral mixture analysis (ON_SMA) is proposed to achieve sub-pixel IS mapping result in the fourth study. It integrates characteristics of optical and NTL imagery to adaptively select local endmembers. Results illustrate the proposed method yields effective improvement and highlight the potential of NTL data in IS mapping. In the fifth study, GA-SVM IS mapping algorithm is investigated with introduction of the achieved urban-URF-rural spatial structure. The combination of optical, NTL and SAR imagery is discussed. GA is implemented for feature selection and parameter optimization in each urban scenario

FINE SCALE MAPPING OF LAURENTIAN MIXED FOREST NATURAL HABITAT COMMUNITIES USING MULTISPECTRAL NAIP AND UAV DATASETS COMBINED WITH MACHINE LEARNING METHODS

Natural habitat communities are an important element of any forest ecosystem. Mapping and monitoring Laurentian Mixed Forest natural communities using high spatial resolution imagery is vital for management and conservation purposes. This study developed integrated spatial, spectral and Machine Learning (ML) approaches for mapping complex vegetation communities. The study utilized ultra-high and high spatial resolution National Agriculture Imagery Program (NAIP) and Unmanned Aerial Vehicle (UAV) datasets, and Digital Elevation Model (DEM). Complex natural vegetation community habitats in the Laurentian Mixed Forest of the Upper Midwest. A detailed workflow is presented to effectively process UAV imageries in a dense forest environment where the acquisition of ground control points (GCPs) is extremely difficult. Statistical feature selection methods such as Joint Mutual Information Maximization (JMIM) which is not that widely used in the natural resource field and variable importance (varImp) were used to discriminate spectrally similar habitat communities. A comprehensive approach to training set delineation was implemented including the use of Principal Components Analysis (PCA), Independent Components Analysis (ICA), soils data, and expert image interpretation. The developed approach resulted in robust training sets to delineate and accurately map natural community habitats. Three ML algorithms were implemented Random Forest (RF), Support Vector Machine (SVM), and Averaged Neural Network (avNNet). RF outperformed SVM and avNNet. Overall RF accuracies across the three study sites ranged from 79.45-87.74% for NAIP and 87.31-93.74% for the UAV datasets. Different ancillary datasets including spectral enhancement and image transformation techniques (PCA and ICA), GLCM-Texture, spectral indices, and topography features (elevation, slope, and aspect) were evaluated using the JMIM and varImp feature selection methods, overall accuracy assessment, and kappa calculations. The robustness of the workflow was evaluated with three study sites which are geomorphologically unique and contain different natural habitat communities. This integrated approach is recommended for accurate natural habitat community classification in ecologically complex landscapes

A Fusion of Remotely Sensed Data to Map the Impervious Surfaces of Growing Cities of Punjab, Pakistan

Urban population is expected to exceed 70% of the world’s total by the middle of the 21st century. Thus, growth in number as well as the sizes of the cities are certain in the near future. The urbanization rates will be much higher in the developing countries than the developed. Such phenomena are accompanied by conversion of land cover from its natural use to built up environment to accommodation growing population. Built up surfaces include road networks, buildings, parking lots and pathways. They are permanently impervious and hydrologically active surfaces. Large volume and discharges of runoff characterize impervious surfaces with frequently occurring flash floods in cities. Besides, proliferating impervious surface is responsible for increasing surface temperature due to Urban Heat Island effect and are the major Nonpoint Source pollutants in the receiving water bodies. At the face of climate change, the consequence of urbanization and increasing impervious surface is exacerbating. Therefore, for sustainable development, spatial and temporal expansion information of impervious surfaces is essential to the planners. Thus, the overall objective of this thesis is to map the impervious surfaces and estimate the expansion rates in the growing cities of Punjab, Pakistan in the last four years. In this thesis, combined and individual datasets from Sentinel-1 and Sentinel-2 satellites were used to extract the amounts of impervious surfaces at city scale and to estimate the expansion rates of various cities of Punjab, Pakistan. The study period for the change analysis is from 2015-2021 based on the availability of satellite imagery. The satellite imageries were obtained from the Copernicus Services Data Hub. Information on different land covers in the form of reflectance, backscattering signal, and texture from a wide range of electromagnetic spectrum of light derived from Sentinels were used to map impervious surfaces. The following land covers were defined: barren soil, vegetation, water, and built-up surface. Four classification models were created from Random Forest algorithms and trained with land covers samples from Google Earth high resolution imagery. The 10 cities considered in this study were among the 50 cities extensively studied by the Urban Unit Pakistan covering the dynamics of Punjab in terms of urban extent, population distribution, area, and expansion. They make up the 21st largest cities in the province as well as represent spatial distribution from north to south. They include various climatic conditions ranging from arid in Multan to humid subtropical in Rawalpindi. They also represent different topographies of the cities such as plain and hilly. Validation samples for each land cover were also obtained from high resolution images to assess the classified land cover maps. Apart from validation of classified maps, quantitative comparison of resultant impervious surfaces was also conducted. For the purpose, the study used published datasets from Atlas of Urban Expansion and the Copernicus Land Service. If available, administrative boundaries of the cities were also used to define the urban extent. For other cities, coordinates were manually defined. The combined Sentinel datasets were able to improve the overall accuracy and kappa coefficient of the classified maps by up to 11% and 7% respectively. McNemar test revealed that the models trained with fused data performed better than the models trained with optical alone data for land cover classification. The cities were expanding at rates ranging from 0.5% to 2.5% annually. The highest rate was encountered in Rawalpindi-Islamabad which is also the capital city of Pakistan. At least for one of the study years (2015/6 or 2020/21) the area was being overestimated by the single optical data. For instance, the optical data overestimated the impervious area of Lahore by a factor of 1.12 times while that of Bahawalpur by a factor of 1.2 times. The incorrect original results attributed to misclassification of barren soil as built up. This conclusion emphasized that additional information on backscattering signal and texture derived from radar image aided to reduce the misclassified bare soil pixels into built up. Spectrum plots also showed that sigma db and variance bands from radar image added a distinct feature to the classifier to distinguish built-up surfaces from other non built-up surfaces. The built-up surface had the highest value in backscatter signals and variance texture bands. This study emphasized the usefulness of combining freely available remote sensing datasets for updating the city scale impervious surfaces cover information in developing countries. The contribution includes the assessment of suitability of combined Sentinel datasets to map the impervious surface at city scale. It also evaluates the rate of expansion of the cities. In conclusion, the combined radar and optical data can enhance the accuracy of classified maps for impervious cover mapping with benefits in complex topographies to update impervious surface information in developing countries. The results from this study could be used as inputs in hydrological and runoff models for urban studies. Other useful applications could be service allocation, drainage improvement, location determination for low impact development (LID) structures, stormwater utility fee determinations, flood control, and pollutants removal from runoff

Mapping urban surface materials using imaging spectroscopy data

Die Kartierung der städtische Oberflächenmaterialien ist aufgrund der komplexen räumlichen Muster eine Herausforderung. Daten von bildgebenden Spektrometern können hierbei durch die feine und kontinuierliche Abtastung des elektromagnetischen Spektrums detaillierte spektrale Merkmale von Oberflächenmaterialien erkennen, was mit multispektralen oder RGB-Bildern nicht mit der gleichen Genauigkeit erreicht werden kann. Bislang wurden in zahlreichen Studien zur Kartierung von städtischen Oberflächenmaterialien Daten von flugzeuggestützten abbildenden Spektrometern mit hoher räumlicher Auflösung verwendet, die ihr Potenzial unter Beweis stellen und gute Ergebnisse liefern. Im Vergleich zu diesen Sensoren haben weltraumgestützte abbildende Spektrometer eine regionale oder globale Abdeckung, eine hohe Wiederholbarkeit und vermeiden teure, zeit- und arbeitsaufwändige Flugkampagnen. Allerdings liegt die räumliche Auflösung der aktuellen weltraumgestützten abbildenden Spektroskopiedaten bei etwa 30 m, was zu einem Mischpixelproblem führt, welches mit herkömmlichen Kartierungsansätzen nur schwer zu bewältigen ist. Das Hauptziel dieser Studie ist die Kartierung städtischer Materialien mit bildgebenden Spektroskopiedaten in verschiedenen Maßstäben und die gleichzeitige Nutzung des Informationsgehalts dieser Daten, um die chemischen und physikalischen Eigenschaften von Oberflächenmaterialien zu erfassen sowie das Mischpixelproblem zu berücksichtigen. Konkret zielt diese Arbeit darauf ab, (1) photovoltaische Solarmodule mit Hilfe von luftgestützten bildgebenden Spektroskopiedaten auf der Grundlage ihrer spektralen Merkmale zu kartieren; (2) die Robustheit der Stichprobe von städtischen Materialgradienten zu untersuchen; (3) die Übertragbarkeit von städtischen Materialgradienten auf andere Gebiete zu analysieren.Mapping urban surface materials is challenging due to the complex spatial patterns. Data from imaging spectrometers can identify detailed spectral features of surface materials through the fine and continuous sampling of the electromagnetic spectrum, which cannot be achieved with the same accuracy using multispectral or RGB images. To date, numerous studies in urban surface material mapping have been using data from airborne imaging spectrometers with high spatial resolution, demonstrating the potential and providing good results. Compared to these sensors, spaceborne imaging spectrometers have regional or global coverage, high repeatability, and avoid expensive, time-consuming, and labor-intensive flight campaigns. However, the spatial resolution of current spaceborne imaging spectroscopy data (also known as hyperspectral data) is about 30 m, resulting in a mixed pixel problem that is challenging to handle with conventional mapping approaches. The main objective of this study is to perform urban surface material mapping with imaging spectroscopy data at different spatial scales, simultaneously explore the information content of these data to detect the chemical and physical properties of surface materials, and take the mixed-pixel problem into account. Specifically, this thesis aims to (1) map solar photovoltaic modules using airborne imaging spectroscopy data based on their spectral features; (2) investigate the sampling robustness of urban material gradients; (3) analyze the area transferability of urban material gradients