Search NASA⌕ Search

SEARCH · Search NASA

Results for “bin mapping”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

63 records · Page 4

Detecting Change in Urban Areas at Continental Scales with MODIS Data

Urbanization is one of the most important components of global environmental change, yet most of what we know about urban areas is at the local scale. Remote sensing of urban expansion across large areas provides information on the spatial and temporal patterns of growth that are essential for understanding differences in socioeconomic and political factors that spur different forms of development, as well the social, environmental, and climatic impacts that result. However, mapping urban expansion globally is challenging: urban areas have a small footprint compared to other land cover types, their features are small, they are heterogeneous in both material composition and configuration, and the form and rates of new development are often highly variable across locations. Here we demonstrate a methodology for monitoring urban land expansion at continental to global scales using Moderate Resolution Imaging Spectroradiometer (MODIS) data. The new method focuses on resolving the spectral and temporal ambiguities between urban/non-urban land and stable/changed areas by: (1) spatially constraining the study extent to known locations of urban land; (2) integrating multitemporal data from multiple satellite data sources to classify c. 2010 urban extent; and (3) mapping newly built areas (2000-2010) within the 2010 urban land extent using a multi-temporal composite change detection approach based on MODIS 250 m annual maximum enhanced vegetation index (EVI). We test the method in 15 countries in East-Southeast Asia experiencing different rates and manifestations of urban expansion. A two-tiered accuracy assessment shows that the approach characterizes urban change across a variety of socioeconomic/political and ecological/climatic conditions with good accuracy (70-91% overall accuracy by country, 69-89% by biome). The 250 m EVI data not only improve the classification results, but are capable of distinguishing between change and no-change areas in urban areas. Over 80% of the error in the change detection can be related to definitional issues or error propagation, rather than algorithm error. As such, these methods hold great potential for routine monitoring of urban change, as well as for providing a consistent and up-to-date dataset on urban extent and expansion for a rapidly evolving region.

Urban areas↗

Multiple GISS AGCM Hindcasts and MSU Versions of 1979-1998

Multiple realizations of the 1979-1998 time period have been simulated by the Goddard Institute for Space Studies Atmospheric General Circulation Model (GISS AGCM) to explore its responsiveness to accumulated forcings, particularly over sensitive agricultural regions. A microwave radiative transfer postprocessor has produced the AGCM's lower tropospheric, tropospheric and lower stratospheric brightness temperature (Tb) time series for correlations with the various Microwave Sounding Unit (MSU) time series available. MSU maps of monthly means and anomalies were also used to assess the AGCM's mean annual cycle and regional variability. Seven realizations by the AGCM were forced by observed sea surface temperatures (sst) through 1992 to gather rough standard deviations associated with internal model variability. Subsequent runs hindcast January 1979 through April 1998 with an accumulation of forcings: observed ssts, greenhouse gases, stratospheric volcanic aerosols. stratospheric and tropospheric ozone and tropospheric sulfate and black carbon aerosols. The goal of narrowing gaps between AGCM and MSU time series was complicated by MSU time series, by Tb simulation concerns and by unforced climatic variability in the AGCM and in the real world. Lower stratospheric Tb correlations between the AGCM and MSU for 1979-1998 reached as high as 0.91 +/-0.16 globally with sst, greenhouse gases, volcanic aerosol, stratospheric ozone forcings and tropospheric aerosols. Mid-tropospheric Tb correlations reached as high as 0.66 +/-.04 globally and 0.84 +/-.02 in the tropics. Oceanic lower tropospheric Tb correlations similarly reached 0.61 +/-.06 globally and 0.79 +/-.02 in the tropics. Of the sensitive agricultural areas considered, Nordeste in northeastern Brazil was simulated best with mid-tropospheric Tb correlations up to 0.75 +/- .03. The two other agricultural regions, in Africa and in the northern mid-latitudes, suffered from higher levels of non-sst variability. Zimbabwe had a maximum mid-tropospheric correlation of 0.54 +/- 0.11 while the U.S. Cornbelt had only 0.25 +/- .10. Precipitation and surface temperature performance are also examined over these regions. Correlations of MSU and AGCM time series mostly improved with addition of explicit atmospheric forcings in zonal bands but not in agricultural regional bins each encompassing only six AGCM gridcells.

Shah, Kathryn Pierce↗

Mapping Impervious Surfaces Globally at 30m Resolution Using Global Land Survey Data

Impervious surfaces, mainly artificial structures and roads, cover less than 1% of the world's land surface (1.3% over USA). Regardless of the relatively small coverage, impervious surfaces have a significant impact on the environment. They are the main source of the urban heat island effect, and affect not only the energy balance, but also hydrology and carbon cycling, and both land and aquatic ecosystem services. In the last several decades, the pace of converting natural land surface to impervious surfaces has increased. Quantitatively monitoring the growth of impervious surface expansion and associated urbanization has become a priority topic across both the physical and social sciences. The recent availability of consistent, global scale data sets at 30m resolution such as the Global Land Survey from the Landsat satellites provides an unprecedented opportunity to map global impervious cover and urbanization at this resolution for the first time, with unprecedented detail and accuracy. Moreover, the spatial resolution of Landsat is absolutely essential to accurately resolve urban targets such a buildings, roads and parking lots. With long term GLS data now available for the 1975, 1990, 2000, 2005 and 2010 time periods, the land cover/use changes due to urbanization can now be quantified at this spatial scale as well. In the Global Land Survey - Imperviousness Mapping Project (GLS-IMP), we are producing the first global 30 m spatial resolution impervious cover data set. We have processed the GLS 2010 data set to surface reflectance (8500+ TM and ETM+ scenes) and are using a supervised classification method using a regression tree to produce continental scale impervious cover data sets. A very large set of accurate training samples is the key to the supervised classifications and is being derived through the interpretation of high spatial resolution (approx. 2 m or less) commercial satellite data (Quickbird and Worldview2) available to us through the unclassified archive of the National Geospatial Intelligence Agency (NGA). For each continental area several million training pixels are derived by analysts using image segmentation algorithms and tools and then aggregated to the 30m resolution of Landsat. Here we will discuss the production/testing of this massive data set for Europe, North and South America and Africa, including assessments of the 2010 surface reflectance data. This type of analysis is only possible because of the availability of long term 30m data sets from GLS and shows much promise for integration of Landsat 8 data in the future.

global land survey↗

An Automated Algorithm to Screen Massive Training Samples for a Global Impervious Surface Classification

An algorithm is developed to automatically screen the outliers from massive training samples for Global Land Survey - Imperviousness Mapping Project (GLS-IMP). GLS-IMP is to produce a global 30 m spatial resolution impervious cover data set for years 2000 and 2010 based on the Landsat Global Land Survey (GLS) data set. This unprecedented high resolution impervious cover data set is not only significant to the urbanization studies but also desired by the global carbon, hydrology, and energy balance researches. A supervised classification method, regression tree, is applied in this project. A set of accurate training samples is the key to the supervised classifications. Here we developed the global scale training samples from 1 m or so resolution fine resolution satellite data (Quickbird and Worldview2), and then aggregate the fine resolution impervious cover map to 30 m resolution. In order to improve the classification accuracy, the training samples should be screened before used to train the regression tree. It is impossible to manually screen 30 m resolution training samples collected globally. For example, in Europe only, there are 174 training sites. The size of the sites ranges from 4.5 km by 4.5 km to 8.1 km by 3.6 km. The amount training samples are over six millions. Therefore, we develop this automated statistic based algorithm to screen the training samples in two levels: site and scene level. At the site level, all the training samples are divided to 10 groups according to the percentage of the impervious surface within a sample pixel. The samples following in each 10% forms one group. For each group, both univariate and multivariate outliers are detected and removed. Then the screen process escalates to the scene level. A similar screen process but with a looser threshold is applied on the scene level considering the possible variance due to the site difference. We do not perform the screen process across the scenes because the scenes might vary due to the phenology, solar-view geometry, and atmospheric condition etc. factors but not actual landcover difference. Finally, we will compare the classification results from screened and unscreened training samples to assess the improvement achieved by cleaning up the training samples. Keywords:

Tan, Bin↗

Put X-Ray Polarimetry on the MAP!

With Prof. R. Novick and others at the Columbia Astrophysics Laboratory I help found the field of X-ray polarimetry in the early 1970s. Currently I have more experience with the design, construction, calibration, and space flight of such instruments than anyone on the planet. The early probing beginnings saw only one definitive measurement (that of the integrated low-energy X-ray emission from the Crab Nebula sans pulsar) and a number of upper limits. The limited success did nevertheless inspire a number of detailed theoretical calculations, concentrating at first on neutron stars and black holes showing how precise measurements (e.g. degree of polarization and position angle as a function of pulse phase) would provide definitive limitations on otherwise equally plausible theoretical models. Over time the theoretical foundation has grown (e.g. the proceedings of the X-Ray Polarimetry Workshop held at SLAC in 2004). I will outline these foundations. It is important to understand the history of X-ray polarimetry beyond the early excitement. A polarimeter was at the focus of the original Einstein mission but was dropped during the restructuring. A polarimeter was successfully proposed (R. Novick PI, I was a Co-I) and built for the original Spectrum-X mission. During the years before the cancellation of Spectrum-X, the potential flight of this device stood in the way of other space flights for polarimeters --- "let us wait and see what it finds". This was unfortunate as there were a number of reasons why that polarimeter should not have been flown on the mission. Perhaps the most significant (but not only) reason was that a shared focal plane provided very little observing time. This is an extremely important point in considering the Roadmap. It is doubtful that there many 100%-polarized sources and so the "signal" is more typically a small fraction of the source flux. Thus, the source itself provides a substantial background, making continuum polarimetry even more difficult than narrow-line spectroscopy. X-ray polarimetry thus requires a dedicated mission that can, without programmatic pressures from other instruments/users, devote the integration time to perform meaningful measurements. The recently cancelled GEMS might have been such a mission. At least it was dedicated to polarimetry. Performing meaningful measurements is not going to be easy. In part because of the long hiatus and lack of experience, there appears to be too much pressure to "sell" polarimetry missions by the number of sources for which one might answer the simple question is, or is not, the integrated and time averaged emission from the source polarized at some confidence level? This was a fine question for the 1970 s but, I maintain, it is not today. One simply doesn't want to measure the time averaged polarization of the Crab s pulsar, but one wants to know the polarization as a function of energy and pulse phase to compare, e.g. to optical and radio measurements which divide even the primary pulse into dozens of phase bins. Such observations can distinguish amongst competing theories for the pulsed emission. The Roadmap should define what meaningful experiments are. I will pose some suggestions. Note that , because X-rays are usually believed originate in either non-thermal or highly aspherical situations we expect X-ray polarimetry to be much more important and rich in astrophysical information as opposed to the visible, where starlight often dominates the emission. One has often dreamt about an instrument that does polarimetry whilst it does other things, and I will discuss this. Even in this case, one needs to realize that the observing time will be driven by the polarimetry, otherwise no useful polarization measurements will be made. Finally, I will discuss some misconceptions that have appeared in the literature and at conferences which indicate to me that certain fundamental principles of polarimeter design and performance are not clearly understoodt is important for the Roadmap that we understand these principles and not advocate instruments that will not (or even worse cannot) accomplish the promised scientific return. I will discuss these principles.

Weisskopf, Martin C.↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗

Transmission Electron Microscopy (TEM) Sample Preparation of Si(1-x)Gex in c-Plane Sapphire Substrate

The National Aeronautics and Space Administration-invented X-ray diffraction (XRD) methods, including the total defect density measurement method and the spatial wafer mapping method, have confirmed super hetero epitaxy growth for rhombohedral single crystalline silicon germanium (Si1-xGex) on a c-plane sapphire substrate. However, the XRD method cannot observe the surface morphology or roughness because of the method s limited resolution. Therefore the authors used transmission electron microscopy (TEM) with samples prepared in two ways, the focused ion beam (FIB) method and the tripod method to study the structure between Si1-xGex and sapphire substrate and Si1?xGex itself. The sample preparation for TEM should be as fast as possible so that the sample should contain few or no artifacts induced by the preparation. The standard sample preparation method of mechanical polishing often requires a relatively long ion milling time (several hours), which increases the probability of inducing defects into the sample. The TEM sampling of the Si1-xGex on sapphire is also difficult because of the sapphire s high hardness and mechanical instability. The FIB method and the tripod method eliminate both problems when performing a cross-section TEM sampling of Si1-xGex on c-plane sapphire, which shows the surface morphology, the interface between film and substrate, and the crystal structure of the film. This paper explains the FIB sampling method and the tripod sampling method, and why sampling Si1-xGex, on a sapphire substrate with TEM, is necessary.

Kim, Hyun Jung↗

Large Eddy Simulation of Cirrus Clouds

The Regional Atmospheric Modeling System (RAMS) with mesoscale interactive nested-grids and a Large-Eddy Simulation (LES) version of RAMS, coupled to two-moment microphysics and a new two-stream radiative code were used to investigate the dynamic, microphysical, and radiative aspects of the November 26, 1991 cirrus event. Wu (1998) describes the results of that research in full detail and is enclosed as Appendix 1. The mesoscale nested grid simulation successfully reproduced the large scale circulation as compared to the Mesoscale Analysis and Prediction System's (MAPS) analyses and other observations. Three cloud bands which match nicely to the three cloud lines identified in an observational study (Mace et al., 1995) are predicted on Grid #2 of the nested grids, even though the mesoscale simulation predicts a larger west-east cloud width than what was observed. Large-eddy simulations (LES) were performed to study the dynamical, microphysical, and radiative processes in the 26 November 1991 FIRE 11 cirrus event. The LES model is based on the RAMS version 3b developed at Colorado State University. It includes a new radiation scheme developed by Harrington (1997) and a new subgrid scale model developed by Kosovic (1996). The LES model simulated a single cloud layer for Case 1 and a two-layer cloud structure for Case 2. The simulations demonstrated that latent heat release can play a significant role in the formation and development of cirrus clouds. For the thin cirrus in Case 1, the latent heat release was insufficient for the cirrus clouds to become positively buoyant. However, in some special cases such as Case 2, positively buoyant cells can be embedded within the cirrus layers. These cells were so active that the rising updraft induced its own pressure perturbations that affected the cloud evolution. Vertical profiles of the total radiative and latent heating rates indicated that for well developed, deep, and active cirrus clouds, radiative cooling and latent heating could be comparable in magnitude in the cloudy layer. This implies that latent heating cannot be neglected in the construction of a cirrus cloud model. The probability density function (PDF) of w was analyzed to assist in the parameterization of cloud-scale velocities in large-scale models. For the more radiatively-driven, thin cirrus case, the PDFs are approximately Gaussian. However, in the interior of the deep, convectively unstable case, the PDFs of w are multi-modal and very broad, indicating that parameterizing cloud-scale motions for such clouds can be very challenging. The results of this research are described in detail in a paper submitted to the Journal of Atmospheric Science (Wu and Cotton, 1999), which is enclosed as Appendix 2. Using soundings extracted from a mesoscale simulation of the November 26, 1991 cirrus event, the radiative effects on vapor deposition/sublimation of ice crystals was studied using a two-dimensional cloud-resolving model (CRM) version of RAMS, coupled to an explicit bin-resolving microphysics. The CRM simulations of the November 26, 1991 cirrus event demonstrate that the radiative impact on the diffusional growth (or sublimation) of ice crystals is significant. In this case, the ice particles experienced radiative warming. Model results show that radiative feedbacks in the diffusional growth of ice particles can be very complex. Radiative warming of an ice particle will restrict the particle's diffusional growth. In the case of radiative warming, ice particles larger than a certain size will experience so much radiative warming that surface ice saturation vapor pressures become large enough to cause sublimation of the larger crystals, while smaller crystals are growing by vapor deposition. However, ice mass production can be enhanced in the case of radiative cooling of an ice particle. For the November 26, 1991 cirrus event, radiative feedback results in significant reduction in the total ice mass, especially in the production of large ice crystals, and consequently, both radiative and dynamic properties of the cirrus cloud are significantly affected. A complete description of this research has been submitted as a paper to the Journal of Atmospheric Science (Wu et al., 1999), and included as Appendix 3.

Wu, Ting↗

Integrated Topographic Corrections Improve Forest Mapping Using Landsat Imagery

In mountainous environments, topography strongly affects the reflectance due to illumination effects and cast shadows, which introduce errors in land cover classifications. However, topographic correction is not routinely implemented in standard data pre-processing chains (e.g., Landsat Analysis Ready Data), and there is a lack of consensus whether topographic correction is necessary, and if so, how to conduct it. Furthermore, methods that correct simultaneously for atmospheric and topographic effects are becoming available, but they have not been compared directly. Our objects were to investigate (1) the effectiveness of two topographic correction approaches that integrate atmospheric and topographic correction, (2) improvements in classification accuracy when analyzing topographically corrected single-date imagery (14 July 2016 and 2 October 2016), versus a full Landsat time series from 2014 to 2016, and 3) improvements in classification accuracy when including additional terrain information (i.e., topographic slope, elevation, and aspect). We developed a physical based model and compared it with an enhanced C-correction, both of which integrate atmospheric and topographic correction. We compared classification accuracies with and without topographic correction using combinations of single-date imagery, image composites and spectral-temporal metrics generated from the full Landsat time series, and additional terrain information in the Caucasus Mountains. We found that both the enhanced C-correction and the physical model performed very well and largely eliminated the correlation (Pearson’s correlation coefficient r ranges from 0.06 to 0.24) between surface reflectance and illumination condition, but the physical model performed best (r ranges from 0.05 to 0.11). Both image composites, and spectral-temporal metrics generated from corrected imagery, resulted in significantly (p ≤ 0.05) higher classification accuracies and better forest classifications, especially for the mixed forests. Adding terrain information reduced classification error significantly, but not as much as topographic correction. In summary, topographic correction remains necessary, even when analyzing a full Landsat time series and including a digital elevation model in the classification. We recommend that topographic correction should be applied when analyzing Landsat satellite imagery in mountainous region for forest cover classification.

Atmospheric correction↗