Remote NASA Goddard’s LiDAR, Hyperspectral and Thermal (G-LiHT) Airborne Imager

: The combination of LiDAR and optical remotely sensed data provides unique information about ecosystem structure and function. Here, we describe the development, validation and application of a new airborne system that integrates commercial off the shelf LiDAR hyperspectral and thermal components in a compact, lightweight and portable system. Goddard’s LiDAR, 6,500 km 2 of data for NASA-sponsored studies across a broad range of ecoregions in the USA and Mexico. In this paper, we document G-LiHT design considerations, physical specifications, instrument performance and calibration and acquisition parameters. In addition, we describe the data processing system and higher-level data products that are freely distributed under NASA’s Data and Information policy.


Introduction
LiDAR, hyperspectral and thermal remote sensing are core areas of current and planned NASA remote sensing capability (e.g., Landsat Enhanced Thematic Mapper Plus, ETM+, and Operational Land Imager, OLI; Earth Observing-1, Hyperion; Earth Observing System's Advanced Spaceborne Thermal Emission and Reflection Radiometer, ASTER, Multi-angle Imaging SpectroRadiometer, MISR, and Moderate Resolution Imaging Spectroradiometer, MODIS; Suomi National Polar-orbiting Partnership Visible Infrared Imaging Radiometer Suite, NPP VIIRS; Ice, Cloud, and Land Elevation Satellite-2 Advanced Topographic Laser Altimeter System, ICESat-2; and Hyperspectral Infrared Imager, HyspIRI), and fusion of complementary data from different sensors offers the potential for improved global remote sensing of terrestrial ecosystems. Formation flying missions, such as NASA's A-Train satellite constellation, offer near-simultaneous coverage from a number of Earth-observing systems; however, the spatial resolution of the current generation of NASA sensors limits some terrestrial ecology applications. Furthermore, efforts to merge 3D structure information from LiDAR with imaging spectrometer and thermal data is currently hindered by the absence of a LiDAR designed for vegetation in space [1] and an EO-1 satellite that is approaching the end of its useful life.
Airborne platforms are more flexible than satellite missions for developing and testing data fusion at fine spatial (1 to 10 m) and spectral resolutions, which have prompted a new generation of LiDAR and imaging spectrometer instrument packages [2][3][4]. Airborne platforms offer specific advantages for the study of terrestrial ecosystem, including targeted acquisitions of seasonal and diurnal processes (e.g., wetland inundation, plant phenology, drought and fire impacts) and coordinated field and remote sensing data collection needed to scale process-level understanding to the scale of airborne and satellite remote sensing observations [5]. Furthermore, the addition of coincident thermal data to LiDAR and imaging spectrometer data broadens the range of potential terrestrial ecology applications to include research on evapotranspiration, hydrology, forest health and urban applications [6,7].
Here, we describe the development of a unique multi-sensor instrument, Goddard's LiDAR, Hyperspectral and Thermal (G-LiHT) airborne imager, which advances previous concepts for data fusion by integrating LiDAR, hyperspectral and thermal sensors in a lightweight and portable system for worldwide research applications. G-LiHT's single-solution GPS-INS (Global Positioning System and Inertial Navigation System) avoids multi-dimensional data effects that are introduced when data is collected at different time and observational scales [2]. Imaging spectroscopy provides quantitative information on vegetation cover, species composition and biophysical and chemical properties that can be derived from measurements of reflected sunlight in the visible through shortwave infrared wavelengths [8][9][10][11][12]. Light detection and ranging provides quantitative, 3D information on terrain and vegetation cover, height and distribution of canopy elements, which can be used to characterize biodiversity and habitat [13] and estimate light interception in plant photosynthesis and production models [14]. Land surface temperature provides data needed to estimate evapotranspiration and other surface energy fluxes [6], and can be an indicator of soil or vegetation moisture status [15]. By using commercial off-the-shelf instrumentation and general aviation aircraft, G-LiHT also reduces development and operational costs. Low-cost deployment of the G-LiHT system opens a wide range of applications for targeted, airborne remote sensing, including diurnal and seasonal processes in terrestrial ecosystems at ~1 m spatial resolution ( Figure 1).

Figure 1.
First coincident acquisition of passive optical, thermal and LiDAR data with G-LiHT (14 July 2011; 37.1839°N 76.5291°W) and key measurement characteristics of the instruments. Spectral and structural differences between a forest, river, golf course and buildings demonstrate the synergistic potential of data fusion for airborne remote sensing of ecosystem composition, structure, function and health. FOV, field of view; NETD, Noise Equivalent Temperature Difference.
The goal of this multi-sensor and data fusion effort is to characterize ecosystem form and function using remote sensing data, with a particular emphasis on the data products needed to develop a new generation of high resolution ecosystem and radiative transfer models. Analysis of G-LiHT data will be used to: • provide new insight into photosynthetic functionality and vegetation productivity, including new, spatially-explicit remote sensing indicators of key dynamic biological processes; • characterize fine-scale spatial and temporal heterogeneity in ecosystem structure and function under diverse environmental and climate conditions; and • create new methods for data fusion to monitor ecosystem health and the effects of climate and human-induced changes on these ecosystems.

Scientific Objectives and Design Considerations
G-LiHT was designed as a relatively inexpensive, robust and portable research tool for evaluating the potential benefits of data fusion for studies of terrestrial ecosystems. Table 1 lists the objectives and related measurement requirements for G-LiHT. One of the major obstacles to the development of science-based data fusion algorithms is the availability of accurately co-registered data of similar grain size for different information types. This is often the case when instruments are flown on different platforms and acquired on different dates. We believe that "instrument fusion" is a prerequisite to "data fusion" and conceived G-LiHT as a multi-sensor airborne imaging system that would simultaneously map the composition, structure and function of terrestrial ecosystems.
In addition, G-LiHT was designed to simplify worldwide deployment and minimize collection and data processing costs. As a result, G-LiHT features eye-safe lasers, a portable, low-power payload (37 kg; 30 × 30 × 60 cm; 210 W), a single-solution GPS-INS, compatibility with common, civilian-use aircraft (e.g., Cessna, Piper, Twin Otter; 12/28 VDC compatibility) using a standard camera port or custom wing-mounted pod and commercial off-the-shelf (COTS) instruments that are easy to replace and not regulated by ITAR (International Traffic in Arms Regulation). Both G-LiHT science and deployment objectives are traced to instrument and design requirements in Table 1.   Specifications of the individual G-LiHT instruments are provided in the sections below.

GPS and Inertial Navigation System (INS)
A single-solution GPS-INS (RT-4041, Oxford Technical Solutions, Oxfordshire, UK) is used to obtain high precision position and attitude measurements for all G-LiHT sensors. The unit is directly attached to the airborne laser scanning (ALS), which shares the same mounting plate as the optical breadboard for the imaging spectrometer and thermal camera ( Figure 2). The GPS-INS incorporates a six-axis inertial navigation system (three gyroscopes and three servo-grade accelerometers) and an L1/L2 GPS (Global Positioning System) and GLONASS (Global Navigation Satellite System) receiver with an OmniStar decoder to deliver 10 cm positioning, 0.1° heading and 0.03° roll and pitch accuracies. Measurements are acquired at 250 Hz, which is required to geolocate data from the optical imagers (25 to 50 Hz) and ALS (300 kHz laser). The real-time internal processing includes strap down algorithms, a WGS-84 and EGM96 (Earth Gravitational Model 1996) Earth model, Kalman filtering and in-flight alignment algorithms. The Kalman filter and in-flight algorithms monitor the performance of the system and update the measurements to correct for inertial sensor errors and maintain high positional accuracy.

Airborne Scanning LiDAR
The VQ-480 (Riegl USA, Orlando, FL, USA) airborne laser scanning (ALS) instrument was selected for use with G-LiHT, because it was affordable, compact, provided evenly distributed pulses on the ground and offered near-turnkey operation. The VQ-480 uses a high-performance laser rangefinder and a rotating polygon mirror with three facets to deflect a 1,550 nm Class 1 laser beam onto the ground. A user-selectable pulse repetition rate up to 300 kHz provides an effective measurement rate of up to 150 kHz along a 60° swath perpendicular to the flight direction ( Figure 3 ability to characterize fine scale disturbances, which are difficult to deconvolve from large footprint LiDAR waveforms ( Figure 4, [17]). The mirror speed is set to a maximum of 100 rotations s −1 during G-LiHT acquisitions, whose points are spaced 0.23 m apart within a line perpendicular to the flight direction and 0.57 m between lines with a nominal aircraft speed of 110 knots. For each laser shot, the return waveform is digitized, and online processing algorithms are used to provide ranging data for multiple targets. Up to eight discrete ranging returns may be identified and recorded for a given pulse. Each return target is time tagged and accurately synchronized with the GPS-INS using a 1 Hz Transistor-Transistor Logic (TTL) pulse. Riegl's software, RiACQUIRE, provides a graphical user interface for scanner control and near real-time monitoring of scanning LiDAR and GPS-INS data. Discrete returns from small footprint LiDAR are able to detect small gaps and characterize fine-scale disturbance (i.e., strip harvesting), which are challenging to deconvolve from large footprint LiDAR waveforms.

Profiling LiDAR
G-LiHT's profiling LiDAR is an LD321-A40 (Riegl USA, Orlando, FL, USA), multi-purpose laser distance meter that is similar to the LD90-3800-VHS used in the Portable Airborne Laser System (PALS) [16]. The LD321-A40 provides a continuous profile of canopy height measurements, which, from a regional sampling perspective, provides data similar to other space-based profilers that have been flown (Ice, Cloud, and Land Elevation Satellite Geoscience Laser Altimeter System, ICESat GLAS), proposed for flight (Deformation, Ecosystem Structure and Dynamics of Ice, DESDynI) or which may be flown in the near future (ICESat-2 ATLAS). In addition, profiling data is important for studying horizontal landscape patch structure [13] and ensuring continuity between ALS and PALS datasets that have been collected worldwide. Real-time digital echo signal processing with the LD321-A40 enables precise distance measurement for complex multi-target situations, resolving up to five target distances per pulse. Distance measurements parallel to the flight line are used to continuously measure vegetation height and structure along a sampling transect. A Class 1M laser diode emits a 905 nm beam with a divergence of 1.5 mrad to produce a 50 cm diameter footprint at the nominal operating altitude of 335 m.

Irradiance Spectrometer
Downwelling irradiance is measured with an Ocean Optics (Dunedin, FL, USA) USB4000-VIS-NIR spectrometer. Light energy is transmitted to the spectrometer through an upward looking opaline glass cosine diffuser with a 180° field of view (FOV). The cosine diffuser is mounted on top of the aircraft, where there is an unobstructed view of the sky. A custom mount was designed for the leading edge of the Cessna 206 ( Figure 5), which integrates both the GPS antenna and cosine diffuser. A 3 m long, 100 µm diameter optical fiber delivers the light energy through a 25 µm entrance slit and a multi-bandpass order-sorting filter, where it is dispersed with a fixed grating across a 3,648-element Toshiba linear Charge Coupled Device (CCD) array. The spectrometer covers the spectral range from 350 to 1,100 nm with a native optical resolution to ~1.5 nm (full width half maximum, FWHM). Power and communications are transferred through a USB 2.0 connection to the computer. The downwelling radiometer is operated continuously during flight with a 33 ms integration time, a 30 scan average and 3× spectral binning, which is equivalent to a 1 Hz data acquisition rate. The unit is cross-calibrated with the imaging spectrometer to enable atmospheric characterization of downwelling irradiance, which is used to compute the at-sensor reflectance product.

Imaging Spectrometer
The Hyperspec imaging spectrometer (Headwall Photonics, Fitchburg, MA, USA) enables high spectral and spatial resolution imaging by using f/2.0 telecentric optics and a high efficiency aberration-corrected convex holographic diffraction grating, providing an optical dispersion of 100 nm per mm over a 7.4 mm spatial by 6.0 mm spectral focal plane. The imaging spectrometer is based on the Offner form and is designed to operate in the 400-1,000 nm spectral region with a 50° full field of view. The Hyperspec imaging spectrometer accepts a C-mount objective lens (Cinegon f/1.4 8 mm, Schneider Optics, Hauppauge, NY, USA) with high optical performance using ultralow dispersion glass and broadband anti-reflection coating designed for the visible to near-IR spectrum. Coupled to a b c the spectrometer is the RA1000 m/D high speed rugged megapixel focal plane array that allows 50 progressive frames per second to be acquired through an EPIX (Buffalo Grove, IL, USA) PIXCI ECB1 PCI Express CameraLink interface. The camera uses a 1,004 × 1,004 pixel 2/3-inch format interline CCD with 7.4 µm square pixels and 12-bit radiometric resolution. Other camera features include a digital fine gain for adjustable camera sensitivity over a 60 dB dynamic range, electronic shuttering and low smear characteristics. The camera is controlled with the serial communication channel of the CameraLink interface. Each image frame is coded with a computer timestamp synchronized with a Time Tools (Dudley, UK) LC2750 GPS Timing Receiver for post-processing geolocation using data from the GPS-INS detailed in Section 2.2 ( Figure 6).

Thermal Imaging
G-LiHT's thermal remote sensing capability originates from the Gobi-384 thermal imaging camera (Xenics, Leuven, Belgium). The instrument measures long wave infrared (LWIR) radiation over a broad spectral range (8 to 14 µm) using an uncooled amorphous silica (ASi) microbolometer detector with 384 cross track pixels and a 30° FOV. This compact camera (72 × 60 × 50 mm) uses a Gigabit Ethernet TCP/IP (Transmission Control Protocol and the Internet Protocol) interface to deliver 16-bit radiometrically calibrated thermal imaging data at a 25 Hz frame rate to the computer running Xenics's Xeneth infrared camera software. Each image frame is coded with a computer timestamp synchronized to a GPS Timing Receiver (see Section 2.6).

Boresight Alignment
Georeferenced LiDAR returns and optical observations are computed as a function of lever arm offsets, boresight angles and optical measurements. Bias in any of these parameters and GPS-INS errors will result in offset point clouds and distorted images. Boresight biases are computed using data collected over buildings with peaked roofs, where overlapping data are collected with different headings. Boresight corrections are computed between campaigns and consistently result in an accuracy of ~10 cm (1 σ).

Radiometric Calibration
The absolute spectral response (ASR) of G-LiHT's imaging spectrometer was calculated using measurements made with a portable version of the US National Institute of Standards and Technology's (NIST) Spectral Irradiance and Radiance Responsivity Calibrations Using Uniform Sources (SIRCUS) [18], which is currently housed at NASA Goddard Space Flight Center (GSFC). SIRCUS uses continuously tunable lasers coupled to an integrating sphere as a radiance source for the calibration of detectors operating in the solar reflective spectrum. SIRCUS and appropriately characterized transfer radiometers allow absolute radiometric calibration with uncertainties of 0.5% and traceability to national standards [19]. The availability of lasers determines the spectral coverage on SIRCUS, and the uncertainties achievable are determined by the quality of the transfer radiometers and measurement technique. The high power and wavelength stability of the laser-based sources enable large aperture instruments to be characterized and calibrated. All irradiance and radiance responsivity calibrations are traceable to NIST's Primary Optical Watt Radiometer (POWR) through regularly characterized transfer radiometers. The short calibration chain from NIST's Primary Optical Watt Radiometer (POWR) to the transfer radiometer and imaging spectrometer minimizes the uncertainty of the Spectral Radiometer Facility (SRF). Spectral radiance calibrations at power levels typical of solar irradiance values allowed for quantitative measurements of the spectral response functions and under flight-like conditions over the focal plane from 415 to 1,020 nm. In addition to providing spectral and radiometric calibration, the combination of full-field, full-aperture and near-monochromatic method of radiometric calibration used here allows the unique capability of assessing and correcting for the total system stray light of the imaging spectrometer, not possible with traditional broadband calibration sources.

Wavelength and Radiometric Stability
Both the downwelling irradiance and upwelling imaging spectrometer are monitored for radiometric and wavelength stability. A portable, 10 cm Teflon integrating sphere with a Hg and Ar pen lamp illumination sources produce nine pronounced emission lines (Figure 7) that are used periodically to verify the spectral channel to wavelength relationship over the detector spectral range with sub-nanometer precision. A large aperture, 1 m integrating sphere and an NIST traceable uniform source is used to monitor the radiometric stability and performance of the spectrometers between campaigns.

Thermal Radiometric Calibration
G-LiHT's thermal imaging camera has factory radiometric calibration for surface temperatures between −20 °C and 120 °C, with a Noise Equivalent Temperature Difference (NETD) >50 mK. This wide dynamic range allows for airborne operations of the thermal image over a diverse range of surface temperatures. Thermal calibration stability as a function of a microbolometer operating temperature range of 25 °C to 50 °C was verified against a GSFC Calibration Facilities blackbody for stable operation over the target temperature range of −5 °C to 85 °C (Figure 8).  aircraft speed and frame acquisition rate. The instantaneous per pixel resolution of the thermal camera is approximately 3 times coarser than the imaging spectrometer, but oversampling with the 2D detector allows us to create a product of equal resolution (Section 5.2.4). Flight lines are uploaded onto a 696 GPS (Garmin International, Inc., Olathe, KS, USA) for pilot navigation and onto the instrument operator's computer, which is using Reigl's RiAcquire acquisition software. Each instrument is separately controlled and triggered to collect data using vendor-supplied software. Ejectable, solid state hard drives provide fast input/output and quick retrieval upon conclusion of the flight.

Data Products
NASA's Earth Science Data Systems program references data products based on their level of processing, ranging from Level 0 to Level 4 [20]. Level 0 (L0) data products include unprocessed instrument data; Level 1 (L1) products are time-referenced data that have been processed to at-sensor radiometric units; Level 2 (L2) products are geophysical variables derived from L1 products; and Level 3 (L3) products are geophysical variables mapped on a space-time grid scale. G-LiHT data products (L1 through L3) are listed in Table 2 for each of the instruments. Many of these products will help end users create additional higher-order products (e.g., solar and view angles are needed for computing bidirectional corrected reflectance; terrain and canopy heights are needed for orthorectification) and facilitate scientific data analysis (e.g., cloudiness and solar illumination conditions affect canopy photosynthesis) [21].

Data Processing System
An automated data processing system is key to releasing standardized products in a timely manner. Data processing begins with the GPS-INS and scanning LiDAR data, since some of the products are used in processing the imaging spectrometer and thermal data. Workflows for the scanning LiDAR and spectrometer data have been developed as illustrated in Figures 9 and 10. Manual pre-processing of L0 data with vendor supplied software is the first step in the workflows, followed by automated data processing algorithms that have been custom coded in the IDL-ENVI (Interactive Data Language and Environment for Visualizing Images) scientific programing language (Exelis Visual Information Solutions, Boulder, CO, USA). The following sections describe the specific data processing steps.

GPS and Inertial Data
GPS-INS data is stored on internal memory in a raw, unprocessed format. Lever arm and other offset coefficients are applied to the L0 data with Oxford's RT Post Process software, which converts the proprietary binary data to an ASCII format. During this step, the inertial data is processed forwards and backwards in time to minimize the effects of GPS data drift. During data acquisition, OmniStar HP differential correction service provides a real-time positioning error of <15 cm in extended periods of open sky. Increased precision can be achieved where base station GPS data is available.

Scanning LiDAR Data
Riegl's RiPROCESS software is used for managing, processing, analyzing and visualizing data acquired with the ALS system. RiPROCESS ingests raw laser scanner data and pre-processed GPS-INS data, applies calibration information, transforms the scan data into geographic coordinates and exports return data in LAS file format [22]. For mapping projects, individual swaths are co-aligned for small differences (typically <0.5 m) in elevation due to differences in position uncertainty, which is typically an effect of different satellite configurations and atmospheric conditions. A cross-track swath is typically collected at the conclusion of a mapping project to serve as a reference point for swath co-alignment.   Table; BIL, Band Interleaved by Line.
LiDAR returns are resampled and processed as smaller data volumes (~1 km tiles for mapping and ~7 km segments for sampling transects) for efficient processing. Tiles are buffered to prevent edge artifacts during ground classification and creation of the Digital Terrain Model (DTM). Classification of ground returns is performed with a progressive morphological filter [23]. Delaunay triangulation is used to create a Triangulated Irregular Network (TIN) of ground hits, and the TIN is used to linearly interpolate DTM elevations on a 1 m raster grid. Additionally, the TIN is used to interpolate the base elevation of every non-ground return, and vegetation heights are computed by difference. A Canopy Height Model (CHM) is created by selecting the greatest return height in every 1 m grid cell, using these points to create a TIN and interpolating canopy heights on a 1 m raster grid. An example of the DTM and CHM models are shown in Figure 11. For mapped areas, tiles are mosaicked, and terrain slope and aspect, canopy rugosity (standard deviation of height) and common LiDAR metrics [24,25] are derived from the resulting data at the scale of a US Forest Service Forest Inventory and Analysis (FIA) field subplot (7.32 m radius).

Imaging Spectrometer Data
Pre-processing of the image spectrometer data requires the use of XCAP software (EPIX, Inc., Buffalo Grove, IL, USA) to convert L0 data in a proprietary VIF format to ASCII time stamps and an image cube in standard Band Interleaved by Line (BIL) format. Calibration coefficients (see below) are applied to both the downwelling irradiance data and upwelling imaging spectrometer data to compute at-sensor radiance ( Figure 6). Prior to computing at-sensor reflectance, irradiance spectra are resampled spectrally to match the imaging spectrometer. Two reflectance products are computed as part of the workflow: at-sensor reflectance is computed as the fraction of irradiance detected by the imaging spectrometer, and surface reflectance is computed using the radiometric data and an atmospheric correction algorithm. Following the reflectance calculations, data are aggregated from the 1.5 nm native sampling interval using a Gaussian-based spectral resampling procedure to generate a distribution product with uniform 5 nm FWHM band spacing. This is a finer sampling and spectral resolution than is reported for NASA's Next Generation Airborne Visible/Infrared Imaging Spectrometer (AVIRISng) [26] and may be useful for measuring fine spectral features, such as solar induced fluorescence [27]. From the reflectance data, common vegetation indices for greenness (e.g., Normalized Difference Vegetation Index, (NDVI); Enhanced Vegetation Index, (EVI); red edge), light use efficiency (e.g., Photochemical Reflectance Index) and leaf pigments (e.g., carotenoid, anthocyanin) are computed. Each pixel in the image cube is georeferenced using Geo-Correction for Airborne Platforms (GCAP) [28], a software package developed at NASA Goddard Space Flight Center that provides the user the capability to georeference a raster image using the image time stamps and GPS-INS data. Pixel coordinates are used to build a geographic lookup table in ENVI that is distributed with the radiance data cube and is used to create georeferenced layers for the reflectance data, vegetation indices and ancillary data products.
Ancillary data products include valuable information for each pixel of the image cube, including: acquisition time, solar zenith angle, solar azimuth, incoming PAR, cloudiness index [21], aircraft elevation, aircraft altitude AGL, view angle and view azimuth. These time-space variables are computed from the LiDAR-derived DTM and a clear sky irradiance model [29]. The products will help end users create additional higher-order products (e.g., solar and view angles are needed for computing bidirectional corrected reflectance; terrain and canopy heights are needed for orthorectification) and facilitate scientific data analysis (e.g., cloudiness and solar illumination conditions affect canopy photosynthesis).
Future development will focus on improved co-registration through automated control points and implementing orthorectification and atmospheric correction algorithms specifically designed for low altitude, wide FOV airborne imagers [30].

Thermal Data
By selecting a single line of cross-track pixels from the thermal image array, the thermal data can be treated as a line imager in the same manner as the imaging spectrometer. However, the number of cross track pixels in a single line and frame rate limits the product to a coarse spatial resolution (~3 m) using this method. Since the field of view of the full array is much greater than the spacing between images, we can use spatial oversampling techniques to create a smoothed, fine resolution (1 m) temperature product. We might also use additional information from the LiDAR and imaging spectrometer (e.g., land cover) to resolve non-linear edge effects, such as abrupt land cover boundaries [31]. Both of these methods for deriving fine-resolution surface temperature products are currently under development.

Profiling LiDAR Data
Level 0 profiling LiDAR data is captured in ASCII format through a 10/100 Mbit TCP/IP port, and can be processed in the same manner as PALS data [16].

G-LiHT Data Distribution
NASA's Earth Science Program promotes the full and open sharing of data with all users in accordance with NASA's Data and Information Policy [32], and this includes G-LiHT data products. Every effort is taken by the G-LiHT instrument team to ensure that accurate, well-calibrated data is released in a timely manner, and data is distributed in common, readily usable file formats. Classified LiDAR returns and feature heights are made available in ASPRS LAS file format (American Society for Photogrammetry and Remote Sensing, LASer file format), a non-proprietary, binary file industry standard [22]. This open data file format allows for the raw data and information specific to data collection to be recorded (e.g., scan angle, return number, classification, AGL height). Gridded LiDAR products are made available as GeoTIFF files that conform to established Tagged Image File Format (TIFF) interchange format for georeferenced raster imagery [33]. GeoTIFF files include geographic metadata formally, using compliant TIFF tags and structures. Multiband products from the imaging spectrometer will be distributed in standard Band Interleaved by Line (BIL) format. In addition to the data files, quick look images of key spatial variables (e.g., aircraft trajectory, DTM, CHM, true color image) created in Keyhole Markup Language (KML) [34] for visualization in Google Earth.
The data archive for G-LiHT products can be accessed through the G-LiHT webpage, http://gliht.gsfc.nasa.gov or by directly connecting to the anonymous ftp site, fusionftp.gsfc.nasa.gov/ G-LiHT. The G-LiHT website contains additional instrument specifications, a description of various airborne campaigns, software tools and links to open source data analysis and visualization software. G-LiHT has been used to collect more than 6,500 km 2 of data for NASA-sponsored studies across a broad range of ecoregions in the USA during 2011-2012 ( Figure 12, [35][36][37]), with plans to collect data in Mexico and interior Alaska in 2013-2014. Our automated data processing and distribution system is designed to give scientists open access to both low-and high-level data products (http://gliht.gsfc.nasa.gov), which will stimulate community development of synergistic data fusion algorithms. Other collections in the USA include Forest Inventory and Analysis (FIA) plots and intensive study sites. The green color in (a) represents percent forest cover [35,36]; dark gray lines indicate country borders, and light gray lines delineate terrestrial ecoregions [37].

Conclusions
Goddard's LiDAR, Hyperspectral and Thermal (G-LiHT) Airborne Imager is a unique system that permits simultaneous measurements of vegetation structure, foliar spectra and surface temperatures at very high spatial resolution (~1 m). The complementary nature of LiDAR, optical and thermal data provide an analytical framework for the development of new algorithms to map plant species composition, plant functional types, biodiversity, biomass and carbon stocks and plant growth. G-LiHT data will enhance our ability to validate data from existing satellite missions, design new missions and produce data products related to biodiversity and climate change.
The scientific rationale and motivation for G-LiHT is similar to other multi-sensor systems (e.g., Carnegie Airborne Observatory, CAO, Alpha and Beta Systems [3]; CAO Airborne Taxonomic Mapping System, AToMS [2]; National Ecological Observatory's Airborne Observation Platform, NEON AOP [4]), and differences largely exist due to specific mission objectives. G-LiHT was specifically designed to simplify worldwide deployment and minimize collection and data processing costs by using commercial off-the-shelf instruments and local general aviation aircraft. In contrast, NEON AOP is an operational observatory that will acquire annual acquisitions over selected NEON sites in the US and additional in-country flights requested by the scientific community. CAO systems have largely focused on sustainable forest management and habitat conservation in global tropical forests, but unlike G-LiHT and NEON, the data is not openly distributed. Other differences between these systems include choice of instrumentation. NEON and CAO systems use a custom-built, high fidelity visible to shortwave infrared spectrometer (380 to 2,510 nm) that covers a wider range of the electromagnetic spectrum and has a greater signal-to-noise ratio than the G-LiHT spectrometer (400 to 1,000 nm), but the instrument cost is one-hundred-fold greater, and a more capable aircraft is needed to accommodate the greater size, mass and power requirements. Additionally, G-LiHT is the only system that currently acquires downwelling irradiance and surface temperature measurements. Comparisons are currently under way at locations where both NEON and G-LiHT data have been collected to evaluate the performance and cross-calibrate the LiDAR and spectrometer data products.
Additional G-LiHT instrument specifications, campaign information and access to more than 6,500 km 2 of data can be obtained from the through the G-LiHT webpage, http://gliht.gsfc.nasa.gov. not constitute endorsement, recommendation or favoring by the US government or National Aeronautics and Space Administration.