Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sparse sampling”

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.

At least 55 records · Page 3

Artificial Neural Networks for Determining Magnetospheric Conditions

This chapter presents a neural-network-based technique that allows for the reconstruction of the global, time-varying distribution of some physical quantity Q, that has been sparsely sampled at various locations within the magnetosphere, and at different times. We begin with a general introduction to the problem of prediction and specification, and why it is important and difficult to achieve with existing methods. We then provide a basic introduction to neural networks, and describe our technique using the specific example of reconstructing the electron plasma density in the Earth's inner magnetosphere on the equatorial plane. We then show more advanced uses of the technique, including 3D reconstruction of the plasma density, specification of chorus and hiss waves, and energetic particle fluxes. We summarize and conclude with a general discussion of how machine learning techniques might be used to advance the state-of-the-art in space weather prediction, and insight discovery.

Bortnik, Jacob↗

TPSAS-NF1676L-17943-DND

SAGE II used the solar occultation technique to measure atmospheric ozone for 21 years between 1984 and 2005 and has continued to play a key role in numerous international assessments. Traditionally, the determination of trends in ozone using SAGE II data has been done via time series analysis of monthly mean data within specific latitude bands. However, given the sparse sampling of SAGE II data, this will introduce temporal and/or spatial biases, particularly in the period after late 2000 when the instrument operated at half-duty cycle. Herein we present a new method for performing time series analysis of unevenly sampled data (in particular SAGE II); namely, we utilize the daily mean data and perform a simultaneous temporal and spatial analysis, making use of a Fourier series to constrain latitudinal variation. This methodology is applied to both the previous (v6.2) and current (v7.0) versions of the SAGE II data. Results are discussed that show the improvements to the data between versions 6.2 and 7.0 that make time series analysis produce reduced residuals as well as the preliminary results of the time series itself, which reveal some of the limitations in currently used proxies for the regression.

Robert P Damadeo↗

Remote Sensing of Lineage Functional Types for Modeling and Monitoring Biodiversity

Hyperspectral remote sensing has the potential to continuously scale plant function and plant diversity information from landscape to global extents. Numerous studies have indicated that VSWIR (400-2500 nm) reflectance properties of vegetation capture evolutionarily conserved biochemical, structural, and other functional attributes of plant species. Spectral properties conserved in plants provide the opportunity to both 1) aggregate species into lineages with improved classification accuracy and 2) link those lineages directly to plant traits. Full realization of this goal will enable parameterization of Land Surface Models (LSMs) with remotely sensed information, e.g., canopy nitrogen, and better representations of biodiversity and functional diversity in biogeographic studies. In this study, we use hyperspectral AVIRIS data from the 2013 HyspIRI campaign over the Southern Sierra Nevada, California flight box to investigate the potential for incorporating evolutionary thinking into landcover classification. We link the airborne hyperspectral data with vegetation plot data from roughly 1372 surveys and a phylogeny representing 1361 species. We aggregate species into lineages ranging from species level groups down to similar number of Plant Functional Types as often used in LSMs. We assessed the ability of Random Forest and Partial Least Squares Discriminant Analysis to discriminate across these different phylogenetic scales and determine the optimal number of lineages to classify. Although there are some temporal and spatial differences in our training data, our best approaches achieved moderate classification accuracy (Kappa > 0.65). Given an optimal number of lineages, we explored approaches to improve classifications including machine learning and unmixing approaches. This work suggests that lineage-based methods may be a promising way to leverage the huge amounts of data that will come from high resolution and high return interval hyperspectral data planned for the Surface Biology and Geology mission with sparsely sampled existing ground-based ecological data.

Hyperspectral↗

Assimilation of DAWN Doppler Wind Lidar Data During the 2017 Convective Processes Experiment (CPEX): Impact on Precipitation and Flow Structure

An improved representation of 3-D air motion and precipitation structure through forecast models and assimilation of observations is vital for improvements in weather forecasting capabilities. However, there are few independent data to properly validate a model forecast of precipitation structure when the underlying dynamics are evolving on short convective timescales. Using data from the JPL Ku/Ka-band Airborne Precipitation Radar (APR-2) and the 2 μmDoppler Aerosol Wind (DAWN) lidar collected during the2017 Convective Processes Experiment (CPEX), the NASA Unified Weather Research and Forecasting (WRF) Ensemble Data Assimilation System (EDAS) modeling system was used to quantify the impact of high-resolution sparsely sampled DAWN measurements on the analyzed variables and on the forecast when the DAWN winds were assimilated. Over-all, the assimilation of the DAWN wind profiles had a discernible impact on the wind field as well as the evolution and timing of the 3-D precipitation structure. Analysis of individual variables revealed that the assimilation of the DAWN winds resulted in important and coherent modifications of the environment. It led to an increase in the near-surface convergence, temperature, and water vapor, creating more favorable conditions for the development of convection exactly where it was observed (but not present in the control run). Comparison to APR-2 and observations by the Global Precipitation Measurement (GPM) satellite shows a much-improved forecast after the assimilation of the DAWN winds – development of precipitation where there was none, more organized precipitation where there was some, and a much more intense and organized cold pool, similar to the analysis of the dropsonde data. The onset of the vertical evolution of the precipitation showed similar radar-derived cloud-top heights, but delayed in time. While this investigation was limited to a single CPEX flight date, the investigation design is appropriate for further investigation of the impact of airborne Doppler wind lidar observations upon short-term convective precipitation forecasts

DAWN↗

Classifying Forest Type in the National Forest Inventory Context with Airborne Hyperspectral and Lidar Data

Forest structure and composition regulate a range of ecosystem services, including biodiversity, water and nutrient cycling, and wood volume for resource extraction. Forest type is an important metric measured in the US Forest Service Forest Inventory and Analysis (FIA) program, the national forest inventory of the USA. Forest type information can be used to quantify carbon and other forest resources within specific domains to support ecological analysis and forest management decisions, such as managing for disease and pests. In this study, we developed a methodology that uses a combination of airborne hyperspectral and lidar data to map FIA-defined forest type between sparsely sampled FIA plot data collected in interior Alaska. To determine the best classification algorithm and remote sensing data for this task, five classification algorithms were tested with six different combinations of raw hyperspectral data, hyperspectral vegetation indices, and lidar-derived canopy and topography metrics. Models were trained using forest type information from 632 FIA subplots collected in interior Alaska. Of the thirty model and input combinations tested, the random forest classification algorithm with hyperspectral vegetation indices and lidar-derived topography and canopy height metrics had the highest accuracy (78% overall accuracy). This study supports random forest as a powerful classifier for natural resource data. It also demonstrates the benefits from combining both structural (lidar) and spectral (imagery) data for forest type classification.

random forest↗

Trend Detection of Atmospheric Time Series: Incorporating Appropriate Uncertainty Estimates and Handling Extreme Events

This paper is aimed at atmospheric scientists without formal training in statistical theory. Its goal is to, 1) provide a critical review of the rationale for trend analysis of the time series typically encountered in the field of atmospheric chemistry; 2) describe a range of trend-detection methods; and 3) demonstrate effective means of conveying the results to a general audience. Trend detections in atmospheric chemical composition data are often challenged by a variety of sources of uncertainty, which often behave differently to other environmental phenomena such as temperature, precipitation rate, or stream flow, and may require specific methods depending on the science questions to be addressed. Some sources of uncertainty can be explicitly included in the model specification, such as autocorrelation and seasonality, but some inherent uncertainties are difficult to quantify, such as data heterogeneity and measurement uncertainty due to the combined effect of short- and long-term natural variability, instrumental stability, and aggregation of data from sparse sampling frequency. Failure to account for these uncertainties might result in an inappropriate inference of the trends and their estimation errors. On the other hand, the variation in extreme events might be interesting for different scientific questions, for example, the frequency of extremely high surface ozone events and their relevance to human health. In this study we aim to, 1) review trend detection methods for addressing different levels of data complexity in different chemical species; 2) demonstrate that the incorporation of scientifically interpretable covariates can outperform pure numerical curve fitting techniques in terms of uncertainty reduction and improved predictability; 3) illustrate the study of trends based on extreme quantiles that can provide insight beyond standard mean or median based trend estimates; and 4) present an advanced method of quantifying regional trends based on the inter-site correlations of multi-site data. All demonstrations are based on time series of observed trace gases relevant to atmospheric chemistry, but the methods can be applied to other environmental data sets.

Trace gas↗

Efficient Onboard Attitude Commanding for Fast Maneuvering of Lunar Reconnaissance Orbiter

To support Lunar Reconnaissance Orbiter’s extended science missions, an algorithm for autonomous optimization of fast occultation avoidance maneuvers was developed. The fast attitude maneuvers are inserted to the spacecraft attitude control system as trajectories to be tracked. Executing maneuvers in this way requires transmission and storage of a large number of time-tagged attitude waypoints on the spacecraft. For more efficient day-to-day operations, an interpolating filter was designed to perform onboard interpolation between sparse samples of commands. The interpolating filter – called FastXMan for efficient fast maneuvering – was patched into LRO’s flight software in late 2023 and is presently operational. This paper presents an overview of the issues related to the practical implementation of the filter and illustrates the flight performance of the new scheme.

Operations↗

Extracting Periodic Signals From Irregularly Sampled Data

Successive approximations formed in Fourier space. Algorithm extracts periodic signals from sparse, irregularly sampled sets of measurement data. Pertains to data processed via fast Fourier transforms (FFTs). Data represents signal components with initially unknown frequencies spanning large spectral range and includes frequencies not integer multiples of minimum FFT frequency.

Wilcox, Jaroslava Z.↗

Multiwavelength monitoring of the BL Lacertae object PKS 2155-304. 4: Multiwavelength analysis

Simultaneous X-ray, ultraviolet, optical, infrared, and radio monitoring data were used to test and constrain models of continuum emission from the BL Lacertae object PKS 2155-304. Intensively sampled ultraviolet and soft X-ray light curves showed a clear temporal correlation with the X-rays leading the ultraviolet by 2-3 hr. This lag was found to be significantly different from zero after an exhaustive comparison of four different techniques for measuring temporal correlations. Variations in the ultraviolet trough optical wave bands were also strongly correlated, with no measurable lag down to limiting timescales of approximately less than 1-2 hr. This strong correlation extends to the near-infrared, but the less intensive sampling precludes measurement of any lag beyomnd an upper limit of approximately less than 1 day. These lags and limits of the order of hours are much shorter than most rapid observed single-band variations. Because of the very sparse radio sampling, it was not possible to measure quantitatively the correlation and lag with shorter wavelengths, but the data do suggest that the radio may lag the optical/ultraviolet by approximately 1 week, with longer delays and weaker variations to longer radio wavelengths. The epoch-folding Q(exp 2) statistic was used to test for periodicity, and no evidence for strict or quasi-periodicity was found in any of the light curves. Because they lead the lower frequencies, the soft X-rays (approximately less than 1 keV) cannot arise from synchrotron self-Compton scattering. These results also rule out the accretion disk model, which predicts a measurable lag between ultraviolet/optical wavelength bands and a correlation between hardness and brightness, neither of which were seen. They are consistent with the entire radio through X-ray continuum arising from direct synchrotron emission from a relativistic jet. However, the tapered jet model, in which the X-ray emission is produced closer in, has problems explaining the magnitude of the ultraviolet/X-ray lag, because the X-ray-emitting electrons have very short lifetimes (t(sub 1/2) much less than 1 s). The result that the lag is much smaller than the variability timescale suggests instead that the radiation may be produced in a flattened region such as a shock front.

Edelson, R.↗

Using MODIS and GLAS Data to Develop Timber Volume Estimates in Central Siberia

The boreal forest is the Earth's largest terrestrial biome, covering some 12 million km2 and accounting for about one third of this planet's total forest area. Mapping of boreal forest's type, structure parameters and biomass are critical for understanding the boreal forest's significance in the carbon cycle, its response to and impact on global climate change. Ground based forest inventories, have much uncertainty in the inventory data, particularly in remote areas of Siberia where sampling is sparse and/or lacking. In addition, many of the forest inventories that do exist for Siberia are now a decade or more old. Thus, available forest inventories fail to capture the current conditions. Changes in forest structure in a particular forest-type and region can change significantly due to changing environment conditions, and natural and anthropogenic disturbance. Remote sensing methods can potentially overcome these problems. Multispectral sensors can be used to provide vegetation cover maps that show a timely and accurate geographic distribution of vegetation types rather than decade old ground based maps. Lidar sensors can be used to directly obtain measurements that can be used to derive critical forest structure information (e.g., height, density, and volume). These in turn can used to estimate biomass components using allometric equations without having to use out dated forest inventory. Finally, remote sensing data is ideally suited to provide a sampling basis for a rigorous statistical estimate of the variance and error bound on forest structure measures. In this study, new remote sensing methods were applied to develop estimates timber volume using NASA's MODerate resolution Imaging Spectroradiometer (MODIS) and unique waveform data of the geoscience laser altimeter system (GLAS) for a 10 deg x 10 deg area in central Siberia. Using MODIS and GLAS data, maps were produced for cover type and timber volume for 2003, and a realistic variance (error bound) for timber volume was calculated for the study area. In this 'study we used only GLAS footprints that had a slope value of less than 10 deg. This was done to avoid large errors due to the effect of slope on the GLAS models. The method requires the integration of new remote sensing methods with available ground studies of forest timber volume conducted in Russian forests. The results were compared to traditional ground forest inventory methods reported in the literature and to ground truth collected in the study area.

Ranson, K. Jon↗

Gravity and count probabilities in an expanding universe

The time evolution of nonlinear clustering on large scales in cold dark matter, hot dark matter, and white noise models of the universe is investigated using N-body simulations performed with a tree code. Count probabilities in cubic cells are determined as functions of the cell size and the clustering state (redshift), and comparisons are made with various theoretical models. We isolate the features that appear to be the result of gravitational instability, those that depend on the initial conditions, and those that are likely a consequence of numerical limitations. More specifically, we study the development of skewness, kurtosis, and the fifth moment in relation to variance, the dependence of the void probability on time as well as on sparseness of sampling, and the overall shape of the count probability distribution. Implications of our results for theoretical and observational studies are discussed.

Bouchet, Francois R.↗

The atmospheric SO2 budget for Pinatubo derived from NOAA-11 SBUV/2 spectral data

Spectral scan data from the NOAA-11 SBUV/2 instrument were used to derive SO2 for three days following the eruption of Mt. Pinatubo - June 19, July 1, and July 17, 1991. Band structure between 300 and 310 nm observed in the backscattered albedo uniquely identifies the presence of SO2. Band ratios are used to infer SO2 amounts to better accuracy (10-20 percent) and sensitivity (about 0.5 milli-atm-cm of SO2) than the TOMS retrieval, but with relatively poor spatial coverage because the measurement is nadir only. Only 7 scans showed detectable SO2 on June 19 when the cloud was still very localized. On July 1 there were 29 scans between 35N and 12S with SO2, with the highest concentration detected over the Atlantic, and on July 17 SO2 was detected in 30 scans around the world, but in decreased concentration. Estimates of the total SO2 budget made after the cloud had spread sufficiently for the sparse SBUV/2 sampling to be adequate indicated that there were 8.4 million metric tons (MMT) of SO2 in the stratosphere on July 1, 1991, and 4.1 MMT remaining on July 17. This corresponds to an e-folding time of about 24 days for the conversion of SO2 to aerosol, and is consistent with an initial injection into the stratosphere of 12-15 MMT of SO2.

Mcpeters, Richard D.↗

Characterization of the Crab Pulsar's Timing Noise

We present a power spectral analysis of the Crab pulsar's timing noise, mainly using radio measurements from Jodrell Bank taken over the period 1982-1989, an interval bounded by sparse data sampling and a large glitch. The power spectral analysis is complicated by nonuniform data sampling and the presence of a steep red power spectrum that can distort power spectra measurement by causing severe power 'leakage'. We develop a simple windowing method for computing red noise power spectra of uniformly sampled data sets and test it on Monte Carlo generated sample realizations of red power-law noise. We generalize time-domain methods of generating power-law red noise with even integer spectral indices to the case of noninteger spectral indices. The Jodrell Bank pulse phase residuals are dense and smooth enough that an interpolation onto a uniform time series is possible. A windowed power spectrum is computed revealing a periodic or nearly periodic component with a period of 568 +/- 10 days and a l/f(exp 3) power-law noise component in pulse phase with a noise strength S(sub infinity)=(1.24 +/- 0.067) x 10(exp 16) cycles(exp 2)/sec(exp 2) over the analysis frequency range f=0.003- 0.1 cycles/day. This result deviates from past analyses which characterized the pulse phase timing residuals as either l/f(sub 4) power-law noise or a quasiperiodic process. The analysis was checked using the Deeter polynomial method of power spectrum estimation that was developed for the case of nonuniform sampling, but has lower spectral resolution. The timing noise is consistent with a torque noise spectrum rising with analysis frequency as f implying blue torque noise, a result not predicted by current models of pulsar timing noise. If the periodic or nearly periodic component is due to a binary companion, we find a mass function f(M) = (6.8 +/- 2.4) x 10(exp -16) solar mass and a companion mass, M(sub c) is greater than or equal to 3.2 solar mass assuming a Crab pulsar mass of 1.4 solar mass.

Scott, D. M.↗

Timber Volume and Biomass Estimates in Central Siberia from Satellite Data

Mapping of boreal forest's type, structure parameters and biomass are critical for understanding the boreal forest's significance in the carbon cycle, its response to and impact on global climate change. The biggest deficiency of the existing ground based forest inventories is the uncertainty in the inventory data, particularly in remote areas of Siberia where sampling is sparse, lacking, and often decades old. Remote sensing methods can help overcome these problems. In this joint US and Russian study, we used the moderate resolution imaging spectroradiometer (MODIS) and unique waveform data of the geoscience laser altimeter system (GLAS) and produced a map of timber volume for a 10degx12deg area in Central Siberia. Using these methods, the mean timber volume for the forested area in the total study area was 203 m3/ ha. The new remote sensing methods used in this study provide a truly independent estimate of forest structure, which is not dependent on traditional ground forest inventory methods.

Ranson, K. Jon↗

The 2010 May Flaring Episode of Cygnus X-3 in Radio, X-Rays, and gamma-Rays

In 2009, Cygnus X-3 (Cyg X-3) became the first microquasar to be detected in the GeV gamma-ray regime, via the satellites Fermi and AGILE. The addition of this new band to the observational toolbox holds promise for building a more detailed understanding of the relativistic jets of this and other systems. We present a rich dataset of radio, hard and soft X-ray, and gamma-ray observations of Cyg X-3 made during a flaring episode in 2010 May. We detect a approx.3-d softening and recovery of the X-ray emission, followed almost immediately by a approx.1-Jy radio flare at 15 GHz, followed by a 4.3sigma gamma-ray flare (E > 100 MeV) approx.1.5 d later. The radio sampling is sparse, but we use archival data to argue that it is unlikely the gamma-ray flare was followed by any significant unobserved radio flares. In this case, the sequencing of the observed events is difficult to explain in a model in which the gamma-ray emission is due to inverse Compton scattering of the companion star's radiation field. Our observations suggest that other mechanisms may also be responsible for gamma-ray emission from Cyg X-3.

Williams, Peter K. G.↗

High-Resolution Infrared Spectroscopic Measurements of Comet 2PlEncke: Unusual Organic Composition and Low Rotational Temperatures

We present high-resolution infrared spectroscopic measurements of the ecliptic comet 2P/Encke, observed on 4-6 Nov. 2003 during its close approach to the Earth, using the Near Infrared Echelle Spectrograph on the Keck II telescope. We present flux-calibrated spectra, production rates, and mixing ratios for H2O, CH3OH, HCN, H2CO, C2H2, C2H6, CH4 and CO. Comet 2P/Encke is a dynamical end-member among comets because of its short period of 3.3 years. Relative to "organics-normal" comets, we determined that 2PlEncke is depleted in HCN, H2CO, C2H2, C2H6, CH4 and CO, but it is enriched in CH3OH. We compared mixing ratios of these organic species measured on separate dates, and we see no evidence of macroscopic chemical heterogeneity in the nucleus of 2P/Encke, however, this conclusion is limited by sparse temporal sampling. The depleted abundances of most measured species suggest that 2P/Encke may have formed closer to the young Sun, before its insertion to the Kuiper belt, compared with "organics-normal" comets - as was previously suggested for other depleted comets (e.g. C/1999 S4 (LINEAR)). We measured very low rotational temperatures of 20 - 30 K for H2O, CH3OH and HCN in the near nucleus region of 2P/Encke, which correlate with one of the lowest cometary gas production rates (approx. 2.6 x 10(exp 27) molecules/s) measured thus far in the infrared. This suggests that we are seeing the effects of more efficient radiative cooling, insufficient collisional excitation, and/or inefficient heating by fast H-atoms (and icy grains) in the observed region of the coma. Its extremely short orbital period, very low gas production rate, and classification as an ecliptic comet, make 2PlEncke an important addition to our growing database, and contribute significantly to the establishment of a chemical taxonomy of comets.

Radeva, Yana L.↗

Global Monitoring of Terrestrial Chlorophyll Fluorescence from Moderate-Spectral-Resolution Near-Infrared Satellite Measurements: Methodology, Simulations, and Application to GOME-2

Globally mapped terrestrial chlorophyll fluorescence retrievals are of high interest because they can provide information on the functional status of vegetation including light-use efficiency and global primary productivity that can be used for global carbon cycle modeling and agricultural applications. Previous satellite retrievals of fluorescence have relied solely upon the filling-in of solar Fraunhofer lines that are not significantly affected by atmospheric absorption. Although these measurements provide near-global coverage on a monthly basis, they suffer from relatively low precision and sparse spatial sampling. Here, we describe a new methodology to retrieve global far-red fluorescence information; we use hyperspectral data with a simplified radiative transfer model to disentangle the spectral signatures of three basic components: atmospheric absorption, surface reflectance, and fluorescence radiance. An empirically based principal component analysis approach is employed, primarily using cloudy data over ocean, to model and solve for the atmospheric absorption. Through detailed simulations, we demonstrate the feasibility of the approach and show that moderate-spectral-resolution measurements with a relatively high signal-to-noise ratio can be used to retrieve far-red fluorescence information with good precision and accuracy. The method is then applied to data from the Global Ozone Monitoring Instrument 2 (GOME-2). The GOME-2 fluorescence retrievals display similar spatial structure as compared with those from a simpler technique applied to the Greenhouse gases Observing SATellite (GOSAT). GOME-2 enables global mapping of far-red fluorescence with higher precision over smaller spatial and temporal scales than is possible with GOSAT. Near-global coverage is provided within a few days. We are able to show clearly for the first time physically plausible variations in fluorescence over the course of a single month at a spatial resolution of 0.5 0.5. We also show some significant differences between fluorescence and coincident normalized difference vegetation indices (NDVI) retrievals.

Fluorescence↗

Global Monitoring of Terrestrial Chlorophyll Fluorescence from Moderate-spectral-resolution Near-infrared Satellite Measurements: Methodology, Simulations, and Application to GOME-2

Globally mapped terrestrial chlorophyll fluorescence retrievals are of high interest because they can provide information on the functional status of vegetation including light-use efficiency and global primary productivity that can be used for global carbon cycle modeling and agricultural applications. Previous satellite retrievals of fluorescence have relied solely upon the filling-in of solar Fraunhofer lines that are not significantly affected by atmospheric absorption. Although these measurements provide near-global coverage on a monthly basis, they suffer from relatively low precision and sparse spatial sampling. Here, we describe a new methodology to retrieve global far-red fluorescence information; we use hyperspectral data with a simplified radiative transfer model to disentangle the spectral signatures of three basic components: atmospheric absorption, surface reflectance, and fluorescence radiance. An empirically based principal component analysis approach is employed, primarily using cloudy data over ocean, to model and solve for the atmospheric absorption. Through detailed simulations, we demonstrate the feasibility of the approach and show that moderate-spectral-resolution measurements with a relatively high signal-to-noise ratio can be used to retrieve far-red fluorescence information with good precision and accuracy. The method is then applied to data from the Global Ozone Monitoring Instrument 2 (GOME-2). The GOME-2 fluorescence retrievals display similar spatial structure as compared with those from a simpler technique applied to the Greenhouse gases Observing SATellite (GOSAT). GOME-2 enables global mapping of far-red fluorescence with higher precision over smaller spatial and temporal scales than is possible with GOSAT. Near-global coverage is provided within a few days. We are able to show clearly for the first time physically plausible variations in fluorescence over the course of a single month at a spatial resolution of 0.5 deg × 0.5 deg. We also show some significant differences between fluorescence and coincident normalized difference vegetation indices (NDVI) retrievals.

remote sensing↗