Using the Minnaert-k parameter derived from CHRIS/PROBA data for forest heterogeneity mapping

CHRIS/PROBA is capable of sampling reflected radiation at five viewing angles over the visible and near-infrared regions of the solar spectrum with a relatively high spatial resolution (~17m). We exploited both the spectral and angular domain of CHRIS data in order to map the surface heterogeneity of an Alpine coniferous forest during winter. In the spectral domain, linear spectral unmixing of the nadir image resulted in a canopy cover map. In the angular domain, pixelwise inversion of the Rahman–Pinty–Verstraete (RPV) model at a single wavelength at the red edge (722 nm) yielded a map of the Minnaert-k parameter that provided information on surface heterogeneity at subpixel scale. Merging both maps resulted in a forest cover heterogeneity map, which contains more detailed information on canopy heterogeneity at the CHRIS subpixel scale than can be obtained from a single-source data set

Merging the Minnaert-k parameter with spectral unmixing to map forest heterogeneity with CHRIS/PROBA data

The Compact High Resolution Imaging Spectrometer (CHRIS) mounted onboard the Project for Onboard Autonomy (PROBA) spacecraft is capable of sampling reflected radiation at five viewing angles over the visible and near-infrared regions of the solar spectrum with high spatial resolution. We combined the spectral domain with the angular domain of CHRIS data in order to map the surface heterogeneity of an Alpine coniferous forest during winter. In the spectral domain, linear spectral unmixing of the nadir image resulted in a canopy cover map. In the angular domain, pixelwise inversion of the Rahman-Pinty-Verstraete (RPV) model at a single wavelength at the red edge (722 nm) yielded a map of the Minnaert-k parameter that provided information on surface heterogeneity at a subpixel scale. However, the interpretation of the Minnaert-k parameter is not always straightforward because fully vegetated targets typically produce the same type of reflectance anisotropy as non-vegetated targets. Merging both maps resulted in a forest cover heterogeneity map, which contains more detailed information on canopy heterogeneity at the CHRIS subpixel scale than is possible to realize from a single-source optical data set

Multiscale Imaging of Evapotranspiration

Evapotranspiration (ET; evaporation + transpiration) is central to a wide range of biological, chemical, and physical processes in the Earth system. Accurate remote sensing of ET is challenging due to the interrelated and generally scale dependent nature of the physical factors which contribute to the process. The evaporation of water from porous media like sands and soils is an important subset of the complete ET problem. Chapter 1 presents a laboratory investigation into this question, examining the effects of grain size and composition on the evolution of drying sands. The effects of composition are found to be 2-5x greater than the effects of grain size, indicating that differences in heating caused by differences in reflectance may dominate hydrologic differences caused by grain size variation. In order to relate the results of Chapter 1 to the satellite image archive, however, the question of information loss between hyperspectral (measurements at 100s of wavelength intervals) laboratory measurements and multispectral (≤ 12 wavelength intervals) satellite images must be addressed. Chapter 2 focuses on this question as applied to substrate materials such as sediment, soil, rock, and non-photosynthetic vegetation. The results indicate that the continuum that is resolved by multispectral sensors is sufficient to resolve the gradient between sand-rich and clay-rich soils, and that this gradient is also a dominant feature in hyperspectral mixing spaces where the actual absorptions can be resolved. Multispectral measurements can be converted to biogeophysically relevant quantities using spectral mixture analysis (SMA). However, retrospective multitemporal analysis first requires cross-sensor calibration of the mixture model. Chapter 3 presents this calibration, allowing multispectral image data to be used interchangeably throughout the Landsat 4-8 archive. In addition, a theoretical explanation is advanced for the observed superior scaling properties of SMA-derived fraction images over spectral indices. The physical quantities estimated by the spectral mixture model are then compared to simultaneously imaged surface temperature, as well as to the derived parameters of ET Fraction and Moisture Availability. SMA-derived vegetation abundance is found to produce substantially more informative ET maps, and SMA-derived substrate fraction is found to yield a surprisingly strong linear relationship with surface temperature. These results provide context for agricultural applications. Chapter 5 investigates the question of mapping and monitoring rice agricultural using optical and thermal satellite image time series. Thermal image time series are found to produce more accurate maps of rice presence/absence, but optical image time series are found to produce more accurate maps of rice crop timing. Chapter 6 takes a more global approach, investigating the spatial structure of agricultural networks for a diverse set of landscapes. Surprisingly consistent scaling relations are found. These relations are assessed in the context of a network-based approach to land cover analysis, with potential implications for the scale dependence of ET estimates. In sum, this thesis present a novel approach to improving ET estimation based on a synthesis of complementary laboratory measurements, satellite image analysis, and field observations. Alone, each of these independent sources of information provides novel insights. Viewed together, these insights form the basis of a more accurate and complete geophysical understanding of the ET phenomenon

Hyperspectral Imaging System Model Implementation and Analysis

In support of hyperspectral imaging system design and parameter trade-off research, an analytical end-to-end model to simulate the remote sensing system pipeline and to forecast remote sensing system performance has been implemented. It is also being made available to the remote sensing community through a website. Users are able to forecast hyperspectral imaging system performance by defining an observational scenario along with imaging system parameters. For system modeling, the implemented analytical model includes scene, sensor and target characteristics as well as atmospheric features, background spectral reflectance statistics, sensor specifications and target class reflectance statistics. The sensor model has been extended to include the airborne ProspecTIR instrument. To validate the analytical model, experiments were designed and conducted. The predictive system model has been verified by comparing the forecast results to ones obtained using real world data collected during the RIT SHARE 2012 collection. Results include the use of large calibration panels to show the predicted radiance consistent with the collected data. Grass radiance predicted from ground truth reflectance data also compare well with the real world collected data, and an eigenvector analysis also supports the validity of the predictions. Two examples of subpixel target detection scenario are presented. One is to detect subpixel wood yellow painted planks in an asphalt playground, and the other is to detect subpixel green painted wood planks in grass. To validate our system performance, the detection performance are analyzed using receiver operating characteristic (ROC) curves in a comprehensive scenario setting. The predicted ROC result of the yellow planks matches well the ROC derived from collected data. However, the predicted ROC curve of green planks differs from collected data ROC curve. Additional experiments were conducted and analyzed to discuss the possible reasons of the mismatch including scene characterization inaccuracy. Several subpixel target detection parameter trade-off analyses are given, including relative calibration error vs SNR, the relationship among probability of detection, meteorological range, pixel fill factor, relative calibration error and false alarm rate. These trade-off analyses explain the utility of this model for hyperspectral imaging system design and research

Master of Science

thesisHyperspectral remote sensing of wildfires combines principles of emitted radiation with advanced spectrometry to model wildfire area and temperature, as well as background land cover classification, at the subpixel level. Yet airborne hyperspectral sensors face problems of inconsistent spatial resolutions and have limited spatial and temporal coverage. A proposed hyperspectral/thermal infrared satellite, the Hyperspectral InfraRed Imager (HyspIRI), will provide hyperspectral data over a spectral range of 350- 2500 nm at a spatial resolution of 60.0 m. Hyperspectral radiance data have previously been shown to allow fire detection and retrieval of fire temperature, although these abilities have not been demonstrated at spatial resolutions coarser than 16.1 m. For this study, four hyperspectral images containing active fires were acquired by the Airborne Visible-Infrared Imaging Spectrometer (AVIRIS), with spatial resolutions ranging from 3.8 to 16.1 m. By resampling these AVIRIS images to coarser spatial resolutions and by modeling fire area, fire temperatures and background land cover, the impacts of spatial resolution on fire detection and temperature retrieval were simulated. Multiple endmember spectral mixture analysis (MESMA) methods were used to model fire temperature and background land cover types. Modeling at coarser spatial resolutions produced larger areas of low fire temperatures with lower modeling error than modeling at finer spatial resolutions. Modeling results comparing 60.0 m data with and without a Gaussian point spread function validated pixel aggregation resampling as a suitable approximation of coarser spatial resolution imagery. Coarser spatial resolution hyperspectral data, such as that collected by the future HyspIRI sensor, are likely to model more fire area and lower temperatures when compared against simultaneously acquired higher spatial resolution data. Increasing the saturation thresholds of SWIR channels could greatly improve the fire detection and temperature modeling capabilities of a HyspIRI-like sensor

Using Lidar to geometrically-constrain signature spaces for physics-based target detection

A fundamental task when performing target detection on spectral imagery is ensuring that a target signature is in the same metric domain as the measured spectral data set. Remotely sensed data are typically collected in digital counts and calibrated to radiance. That is, calibrated data have units of spectral radiance, while target signatures in the visible regime are commonly characterized in units of re°ectance. A necessary precursor to running a target detection algorithm is converting the measured scene data and target signature to the same domain. Atmospheric inversion or compensation is a well-known method for transforming mea- sured scene radiance values into the re°ectance domain. While this method may be math- ematically trivial, it is computationally attractive and is most e®ective when illumination conditions are constant across a scene. However, when illumination conditions are not con- stant for a given scene, signi¯cant error may be introduced when applying the same linear inversion globally. In contrast to the inversion methodology, physics-based forward modeling approaches aim to predict the possible ways that a target might appear in a scene using atmospheric and radiometric models. To fully encompass possible target variability due to changing illumination levels, a target vector space is created. In addition to accounting for varying illumination, physics-based model approaches have a distinct advantage in that they can also incorporate target variability due to a variety of other sources, to include adjacency target orientation, and mixed pixels. Increasing the variability of the target vector space may be beneficial in a global sense in that it may allow for the detection of difficult targets, such as shadowed or partially concealed targets. However, it should also be noted that expansion of the target space may introduce unnecessary confusion for a given pixel. Furthermore, traditional physics-based approaches make certain assumptions which may be prudent only when passive, spectral data for a scene are available. Common examples include the assumption of a °at ground plane and pure target pixels. Many of these assumptions may be attributed to the lack of three-dimensional (3D) spatial information for the scene. In the event that 3D spatial information were available, certain assumptions could be levied, allowing accurate geometric information to be fed to the physics-based model on a pixel- by-pixel basis. Doing so may e®ectively constrain the physics-based model, resulting in a pixel-specific target space with optimized variability and minimized confusion. This body of work explores using spatial information from a topographic Light Detection and Ranging (Lidar) system as a means to enhance the delity of physics-based models for spectral target detection. The incorporation of subpixel spatial information, relative to a hyperspectral image (HSI) pixel, provides valuable insight about plausible geometric con¯gurations of a target, background, and illumination sources within a scene. Methods for estimating local geometry on a per-pixel basis are introduced; this spatial information is then fed into a physics-based model to the forward prediction of a target in radiance space. The target detection performance based on this spatially-enhanced, spectral target space is assessed relative to current state-of-the-art spectral algorithms