Large-Scale Urban Impervous Surfaces Estimation Through Incorporating Temporal and Spatial Information into Spectral Mixture Analysis

With rapid urbanization, impervious surfaces, a major component of urbanized areas, have increased concurrently. As a key indicator of environmental quality and urbanization intensity, an accurate estimation of impervious surfaces becomes essential. Numerous automated estimation approaches have been developed during the past decades. Among them, spectral mixture analysis (SMA) has been recognized as a powerful and widely employed technique. While SMA has proven valuable in impervious surface estimation, effects of temporal and spectral variability have not been successfully addressed. In particular, impervious surface estimation is likely to be sensitive to seasonal changes, majorly due to the shadowing effects of vegetation canopy in summer and the confusion between impervious surfaces and soil in winter. Moreover, endmember variability and multi-collinearity have adversely impacted the accurate estimation of impervious surface distribution with coarse resolution remote sensing imagery. Therefore, the main goal of this research is to incorporate temporal and spatial information, as well as geostatistical approaches, into SMA for improving large-scale urban impervious surface estimation. Specifically, three new approaches have been developed in this dissertation to improve the accuracy of large-scale impervious surface estimation. First, a phenology based temporal mixture analysis was developed to address seasonal sensitivity and spectral confusion issues with the multi-temporal MODIS NDVI data. Second, land use land cover information assisted temporal mixture analysis was proposed to handle the issue of endmember class variability through analyzing the spatial relationship between endmembers and surrounding environmental and socio-economic factors in support of the selection of an appropriate number and types of endmember classes. Third, a geostatistical temporal mixture analysis was developed to address endmember spectral variability by generating per-pixel spatial varied endmember spectra. Analysis results suggest that, first, with the proposed phenology based temporal mixture analysis, a significant phenophase differences between impervious surfaces and soil can be extracted and employed in unmxing analysis, which can facilitate their discrimination and successfully address the issue of seasonal sensitivity and spectral confusion. Second, with the analyzed spatial distribution relationship between endmembers and environmental and socio-economic factors, endmember classes can be identified with clear physical meanings throughout the whole study area, which can effectively improve the unmixing analysis results. Third, the use of the spatially varying per-pixel endmember generated from the geostatistical approach can effectively consider the endmember spectra spatial variability, overcome the endmember within-class variability issue, and improve the accuracy of impervious surface estimates. Major contributions of this research can be summarized as follows. First, instead of Landsat Thematic Mapper (TM) images, MODIS imageries with large geographic coverage and high temporal resolution have been successfully employed in this research, thus making timely and regional estimation of impervious surfaces possible. Second, this research proves that the incorporation of geographic knowledge (e.g. phonological knowledge, spatial interaction, and geostatistics) can effectively improve the spectral mixture analysis model, and therefore improve the estimation accuracy of urban impervious surfaces

Uncertainty Assessment of Spectral Mixture Analysis in Remote Sensing Imagery

Spectral mixture analysis (SMA), a scheme of sub-pixel-based classifications, is one of the widely used models to map fractional land use and land cover information in remote sensing imagery. It assumes that: 1) a mixed pixel is composed by several pure land cover classes (endmembers) linearly or nonlinearly, and 2) the spectral signature of each endmember is a constant within the entire spatial extent of analysis. SMA has been commonly applied to impervious surface area extraction, vegetation fraction estimation, and land use and land cover change (LULC) mapping. Limitations of SMA, however, still exist. First, the existence of between- and within-class variability prevents the selection of accurate endmembers, which results in poor accuracy of fractional land cover estimates. Weighted spectral mixture analysis (WSMA) and transformed spectral mixture analysis (TSMA) are alternate means to address the within- and between- class variability. These methods, however, have not been analyzed systematically and comprehensively. The effectiveness of each WSMA and TSMA scheme is still unknown, in particular within different urban areas. Second, multiple endmember SMA (MESMA) is a better alternative to address spectral mixture model uncertainties. It, nonetheless, is time consuming and inefficient. Further, incorrect endmember selections may still limit model performance as the best-fit endmember model might not be the optimal model due to the existence of spectral variability. Therefore, this study aims 1) to explore endmember uncertainties by examining WSMA and TSMA modeling comprehensively, and 2) to develop an improved MESMA model in order to address the uncertainties of spectral mixture models. Results of the WSMA examination illustrated that some weighting schemes did reduce endmember uncertainties since they could improve the fractional estimates significantly. The results also indicated that spectral class variance played a key role in addressing the endmember uncertainties, as the better performing weighting schemes were constructed with spectral class variance. In addition, the results of TSMA examination demonstrated that some TSMAs, such as normalized spectral mixture analysis (NSMA), could effectively solve the endmember uncertainties because of their stable performance in different study areas. Results of Class-based MEMSA (C-MESMA) indicated that it could address spectral mixture model uncertainties by reducing a lot of the calculation burden and effectively improving accuracy. Assessment demonstrated that C-MEMSA significantly improving accuracy. Major contributions of this study can be summarized as follow. First, the effectiveness of addressing endmember uncertainties have been fully discussed by examining: 1) the effectiveness of ten weighted spectral mixture models in urban environments; and 2) the effectiveness of 26 transformed spectral mixture models in three locations. Constructive guidance regarding handling endmember uncertainties using WSMA and TSMA have been provided. Second, the uncertainties of spectral mixture model were reduced by developing an improved MESMA model, named C-MESMA. C-MESMA could restrict the distribution of endmembers and reduce the calculation burden of traditional MESMA, increasing SMA accuracy significantly

An evaluation of Multiple Endmember Spectral Mixture Analysis applied to Landsat 8 OLI images for mapping land cover in southern Africa\u27s Savanna.

Mapping land cover in southern Africa’s savannas using traditional pixel based remote sensing techniques can be very challenging due to the heterogeneity of its vegetation structure and the spectral difficulty in separating similar land covers across various land uses. In order to overcome these complications, a Multiple Endmember Spectral Mixture Analysis (MESMA) provides a potential remote sensing approach to quantify spectral variation in the physical environment at a subpixel level. The MESMA approach was applied in the study area of the Mayuni Conservancy, in Namibia. Results show that 32.3% of the study area is covered by photosynthetic vegetation (PV), 32.0% by non-photosynthetic vegetation (NPV), 25.2% by bare soil (B) and 10.6% by shade. Post-classification validation shows that MESMA presented a moderate performance in estimating the proportions of land cover types in the study area. However the validation process is limited to the available resources and carries great subjectivity. It is concluded that future research on the matter should include a more consistent investigation on the endmember selection methodology and expand the study area inside of the same ecosystem

Tree Cover Variability in the District of Columbia

Urban forests are increasingly a focus of interest as urbanized populations grow and urban areas expand. Urban forests change as trees are planted, grow, die, and are removed. These processes alter a city's tree cover over time, but this inherent dynamism is poorly understood. Better understanding of how tree cover is a variable land cover component will enhance knowledge of the urban environment and provide new perspectives for management of urban resources. In this study, tree cover variability within a major urban center was observed over a 20 year period. Changes in tree cover proportion were measured in the District of Columbia between 1984-2004 utilizing highly calibrated satellite remote sensing data. Testing of alternate methodologies demonstrated that an approach utilizing support vector regression provided most consistent accuracy across land use types. Tree cover maps were validated using aerial photography imagery and data from field surveys. Between 1984-2004, the city-wide tree cover remained between 22.1(+/-2.9)% and 28.8(+/-2.9)% of total land surface area. The District of Columbia did not experience an overall increase or decrease in total tree canopy area. Spatial patterns of tree cover variability were investigated to identify local scale changes in tree cover and connections with urban land use. Within the city, greatest variability was observed in low density residential zones. Tree cover proportion in these zones declined 7.4(+/-5.4)% in the years between 1990-1996 and recovered after 1996. Changes in tree cover were observed with high resolution aerial photography to determine relative contribution from fluctuation in the number of standing trees and changes in crown sizes. Land cover conversion removed dense tree cover from 50.2 hectares of the city's land surface between 1984-2004. The results demonstrate that tree cover variability in the District of Columbia occurred primarily within low population density residential areas. Neighborhoods within these zones were analyzed to identify factors correlated with tree cover. Implications of the results include enhanced understanding of the possible impact of urban forest management, and how a focus on low density residential zones is appropriate in setting goals for expansion of urban tree cover

Small-Area Population Estimation: an Integration of Demographic and Geographic Techniques

Knowledge of detailed and accurate population information is essential to analyze and address a wide variety of socio-economic, political, and environmental issues and to support necessary planning practices for both public agencies and the private sector. However, such important data are generally only available once every decade through the National Census. Moreover, populations in some rapidly-developing areas may increase quickly, such that this ten-year frequency does not meet the needs of these areas. Therefore, a cost-effective method for population estimation is necessary. To address this issue, this research integrated geographic, sociological, and demographic theories and exploited remotely sensed imagery and geographic information system (GIS) datasets to derive better population estimates at the census block level, the finest level of the national census. Specifically, three new approaches have been proposed in this dissertation to assist in the improvement of small-area population estimation accuracy. First, existing remotely sensed and GIS data have been adopted to estimate two major components of a demographic framework, including the redistribution of newly built dwelling units from the aggregated geographic level to the census block level and the estimation of persons per household (PPH) at such a fine scale. Second, in addition to the use of existing data, new urban environmental indicators were also extracted and employed to improve population estimation. In particular, to implement the automatic enumeration for individual housing units, a new spectral index, biophysical composition index (BCI), has been proposed to derive impervious surface information, a desirable urban environmental parameter. Third, using the extracted high-resolution urban environmental information and GIS data, a new bottom-up method was developed for small-area population estimation at the census block level by incorporating these high-resolution data into the demographic framework. Analyses of the results suggest three major conclusions. First, existing GIS spatial factors, together with demographic information, can assist in improving the accuracy of small-area population estimation. Second, the BCI has a closer relationship with impervious surface area than do other popular indices. Moreover, it was shown to be the most effective index of the four evaluated for separating impervious surfaces and bare soil, which consequently might assist in more accurately deriving fractional land cover values. Third, the use of the new environmental indicators extracted from remote sensing imagery and GIS data and the integration of demographic and geographic approaches has significantly improved the estimation accuracy of housing unit (HU) numbers, PPH, and population counts at the census block level. Therefore, this research contributes to both the remote sensing and applied demography fields. The contribution to the remote sensing field lies in the development of a novel spectral index to characterize urban land for monitoring and analyzing urban environments. This index provided more significant separability between impervious surfaces and bare soil than did other existing indices. Moreover, three major contributions have been made in the field of applied demography: 1) the generation of accurate HU estimates using high-resolution remote sensing and GIS datasets, 2) the development of a model to derive an accurate PPH estimate, and 3) the improvement of small-area population estimation accuracy through the integration of geographic and demographic approaches

Regression-based surface water fraction mapping using a synthetic spectral library for monitoring small water bodies

Small water bodies (SWBs), such as ponds and on-farm reservoirs, are a key part of the hydrological system and play important roles in diverse domains from agriculture to conservation. The monitoring of SWBs has been greatly facilitated by medium-spatial-resolution satellite images, but the monitoring accuracy is considerably affected by the mixed-pixel problem. Although various spectral unmixing methods have been applied to map sub-pixel surface water fractions for large water bodies, such as lakes and reservoirs, it is challenging to map SWBs that are small in size relative to the image pixel and have dissimilar spectral properties. In this study, a novel regression-based surface water fraction mapping method (RSWFM) using a random forest and a synthetic spectral library is proposed for mapping 10 m spatial resolution surface water fractions from Sentinel-2 imagery. The RSWFM inputs a few endmembers of water, vegetation, impervious surfaces, and soil to simulate a spectral library, and considers spectral variations in endmembers for different SWBs. Additionally, RSWFM applies noise-based data augmentation on pure endmembers to overcome the limitation often arising from the use of a small set of pure spectra in training the regression model. RSWFM was assessed in ten study sites and compared with the fully constrained least squares (FCLS) linear spectral mixture analysis, multiple endmember spectral mixture analysis (MESMA), and the nonlinear random forest (RF) regression without data-augmentation. The results showed that RSWFM decreases the water fraction mapping errors by ~ 30%, ~15%, and ~ 11% in root mean square error compared with the linear FCLS, MESMA unmixings, and the nonlinear RF regression without data-augmentation respectively. RSWFM has an accuracy of approximately 0.85 in R2 in estimating the area of SWBs smaller than 1 ha

Tensor-based Hyperspectral Image Processing Methodology and its Applications in Impervious Surface and Land Cover Mapping

The emergence of hyperspectral imaging provides a new perspective for Earth observation, in addition to previously available orthophoto and multispectral imagery. This thesis focused on both the new data and new methodology in the field of hyperspectral imaging. First, the application of the future hyperspectral satellite EnMAP in impervious surface area (ISA) mapping was studied. During the search for the appropriate ISA mapping procedure for the new data, the subpixel classification based on nonnegative matrix factorization (NMF) achieved the best success. The simulated EnMAP image shows great potential in urban ISA mapping with over 85% accuracy. Unfortunately, the NMF based on the linear algebra only considers the spectral information and neglects the spatial information in the original image. The recent wide interest of applying the multilinear algebra in computer vision sheds light on this problem and raised the idea of nonnegative tensor factorization (NTF). This thesis found that the NTF has more advantages over the NMF when work with medium- rather than the high-spatial-resolution hyperspectral image. Furthermore, this thesis proposed to equip the NTF-based subpixel classification methods with the variations adopted from the NMF. By adopting the variations from the NMF, the urban ISA mapping results from the NTF were improved by ~2%. Lastly, the problem known as the curse of dimensionality is an obstacle in hyperspectral image applications. The majority of current dimension reduction (DR) methods are restricted to using only the spectral information, when the spatial information is neglected. To overcome this defect, two spectral-spatial methods: patch-based and tensor-patch-based, were thoroughly studied and compared in this thesis. To date, the popularity of the two solutions remains in computer vision studies and their applications in hyperspectral DR are limited. The patch-based and tensor-patch-based variations greatly improved the quality of dimension-reduced hyperspectral images, which then improved the land cover mapping results from them. In addition, this thesis proposed to use an improved method to produce an important intermediate result in the patch-based and tensor-patch-based DR process, which further improved the land cover mapping results

Multiple Endmember Spectral Mixture Analysis (MESMA) Applied to the Study of Habitat Diversity in the Fine-Grained Landscapes of the Cantabrian Mountains

P. 1-19 ArtículoHeterogeneous and patchy landscapes where vegetation and abiotic factors vary at small spatial scale (fine-grained landscapes) represent a challenge for habitat diversity mapping using remote sensing imagery. In this context, techniques of spectral mixture analysis may have an advantage over traditional methods of land cover classification because they allow to decompose the spectral signature of a mixed pixel into several endmembers and their respective abundances. In this work, we present the application of Multiple Endmember Spectral Mixture Analysis (MESMA) to quantify habitat diversity and assess the compositional turnover at different spatial scales in the fine-grained landscapes of the Cantabrian Mountains (northwestern Iberian Peninsula). A Landsat-8 OLI scene and high-resolution orthophotographs (25 cm) were used to build a region-specific spectral library of the main types of habitats in this region (arboreal vegetation; shrubby vegetation; herbaceous vegetation; rocks–soil and water bodies). We optimized the spectral library with the Iterative Endmember Selection (IES) method and we applied MESMA to unmix the Landsat scene into five fraction images representing the five defined habitats (root mean square error, RMSE 0.025 in 99.45% of the pixels). The fraction images were validated by linear regressions using 250 reference plots from the orthophotographs and then used to calculate habitat diversity at the pixel ( -diversity: 30 30 m), landscape (-diversity: 1 1 km) and regional ("-diversity: 110 33 km) scales and thecompositional turnover ( - and -diversity) according to Simpson’s diversity index. Richness and evenness were also computed. Results showed that fraction images were highly related to reference data (R2 0.73 and RMSE 0.18). In general, our findings indicated that habitat diversity was highly dependent on the spatial scale, with values for the Simpson index ranging from 0.20 0.22 for -diversity to 0.60 0.09 for -diversity and 0.72 0.11 for "-diversity. Accordingly, we found -diversity to be higher than -diversity. This work contributes to advance in the estimation of ecological diversity in complex landscapes, showing the potential of MESMA to quantify habitat diversity in a comprehensive way using Landsat imageryS

Investigating impacts of natural and human-induced environmental changes on hydrological processes and flood hazards using a GIS-based hydrological/hydraulic model and remote sensing data

Natural and human-induced environmental changes have been altering the earth's surface and hydrological processes, and thus directly contribute to the severity of flood hazards. To understand these changes and their impacts, this research developed a GISbased hydrological and hydraulic modeling system, which incorporates state-of-the-art remote sensing data to simulate flood under various scenarios. The conceptual framework and technical issues of incorporating multi-scale remote sensing data have been addressed. This research develops an object-oriented hydrological modeling framework. Compared with traditional lumped or cell-based distributed hydrological modeling frameworks, the object-oriented framework allows basic spatial hydrologic units to have various size and irregular shape. This framework is capable of assimilating various GIS and remotely-sensed data with different spatial resolutions. It ensures the computational efficiency, while preserving sufficient spatial details of input data and model outputs. Sensitivity analysis and comparison of high resolution LIDAR DEM with traditional USGS 30m resolution DEM suggests that the use of LIDAR DEMs can greatly reduce uncertainty in calibration of flow parameters in the hydrologic model and hence increase the reliability of modeling results. In addition, subtle topographic features and hydrologic objects like surface depressions and detention basins can be extracted from the high resolution LiDAR DEMs. An innovative algorithm has been developed to efficiently delineate surface depressions and detention basins from LiDAR DEMs. Using a time series of Landsat images, a retrospective analysis of surface imperviousness has been conducted to assess the hydrologic impact of urbanization. The analysis reveals that with rapid urbanization the impervious surface has been increased from 10.1% to 38.4% for the case study area during 1974 - 2002. As a result, the peak flow for a 100-year flood event has increased by 20% and the floodplain extent has expanded by about 21.6%. The quantitative analysis suggests that the large regional detentions basins have effectively offset the adverse effect of increased impervious surface during the urbanization process. Based on the simulation and scenario analyses of land subsidence and potential climate changes, some planning measures and policy implications have been derived for guiding smart urban growth and sustainable resource development and management to minimize flood hazards

Ash Tree Identification Based on the Integration of Hyperspectral Imagery and High-density Lidar Data

Monitoring and management of ash trees has become particularly important in recent years due to the heightened risk of attack from the invasive pest, the emerald ash borer (EAB). However, distinguishing ash from other deciduous trees can be challenging. Both hyperspectral imagery and Light detection and ranging (LiDAR) data are two valuable data sources that are often used for tree species classification. Hyperspectral imagery measures detailed spectral reflectance related to the biochemical properties of vegetation, while LiDAR data measures the three-dimensional structure of tree crowns related to morphological characteristics. Thus, the accuracy of vegetation classification may be improved by combining both techniques. Therefore, the objective of this research is to integrate hyperspectral imagery and LiDAR data for improving ash tree identification. Specifically, the research aims include: 1) using LiDAR data for individual tree crowns segmentation; 2) using hyperspectral imagery for extraction of relative pure crown spectra; 3) fusing hyperspectral and LiDAR data for ash tree identification. It is expected that the classification accuracy of ash trees will be significantly improved with the integration of hyperspectral and LiDAR techniques. Analysis results suggest that, first, 3D crown structures of individual trees can be reconstructed using a set of generalized geometric models which optimally matched LiDAR-derived raster image, and crown widths can be further estimated using tree height and shape-related parameters as independent variables and ground measurement of crown widths as dependent variables. Second, with constrained linear spectral mixture analysis method, the fractions of all materials within a pixel can be extracted, and relative pure crown-scale spectra can be further calculated using illuminated-leaf fraction as weighting factors for tree species classification. Third, both crown shape index (SI) and coefficient of variation (CV) can be extracted from LiDAR data as invariant variables in tree’s life cycle, and improve ash tree identification by integrating with pixel-weighted crown spectra. Therefore, three major contributions of this research have been made in the field of tree species classification:1) the automatic estimation of individual tree crown width from LiDAR data by combining a generalized geometric model and a regression model, 2) the computation of relative pure crown-scale spectral reflectance using a pixel-weighting algorithm for tree species classification, 3) the fusion of shape-related structural features and pixel-weighted crown-scale spectral features for improving of ash tree identification