Moment distance index

Last updated

The moment distance index (MDI) is a shape-based metric or shape index that can be used to analyze spectral reflectance curves and waveform LiDAR, proposed by Salas and Henebry in 2014. [1] In the case of spectral data, the shape of the reflectance curve should unmask fine points of the spectra usually not considered by existing band-specific indices. It has been used to identify spectral regions for chlorophyll and carotenoids, [2] detect greenhouses using WorldView-2, Landsat, and Sentinel-2 satellite data, [3] [4] identify greenhouse crops, [5] compute canopy heights, [6] estimate green vegetation fraction, [7] and optimize Fourier-transform infrared (FTIR) scans for soil spectroscopy. [8]

Various approaches have been devised to analyze medium and fine spectral resolution data and maximize their use to extract specific information for vegetation biophysical and biochemical properties. Combinations of spectral bands, called indices, have been used to diminish the effects of soil background and/or atmospheric conditions while highlighting specific spectral features associated with plant or canopy properties. Vegetation indices (VIs) use the concept of band ratio and differences or weighted linear combinations to take advantage of the visible and NIR bands, two important spectral bands for vegetation studies, in measuring the photosynthetic activity of the plant and explore vegetation dynamics. There is an extensive list of such indices, including the normalized difference vegetation index (NDVI), ratio-based indices such as the modified simplerRatio, soil-distance-based indices such as the modified soil adjusted vegetation index (MSAVI), and many others. Whereas most indices incorporate two-band or three-band relations – slope-based, distance-based on soil line or optimized (slope-based and distance-based concepts combined) – no approach deals with the raw shape of the spectral curve. MDI, however, investigates the shape of the reflectance curve using multiple spectral bands not considered by other indices, which could carry additional spectral information useful for vegetation monitoring.

A full-waveform Light Detection And Ranging (LiDAR) system has the ability to record many returns per emitted pulse, as a function of time, to reveal the vertical structure of the illuminated object, showing position of the individual targets, and finer details of the signature of intercepted surfaces or the proportion of the canopy complexity. Information associated with the illuminated object can be decoded from the generated backscattered waveform, as key features of the waveform such as the shape, area, and power are directly related to the geometry of the illuminated object. The richness of the LiDAR waveform holds a promise to address the challenge of characterizing in detail the geometric and reflection characteristics of vegetation structure, e.g., the vertical canopy volume distribution. MDI utilizes the raw waveform and place importance on its shape and its return power. MDI departs from the usual Gaussian modeling in detecting peaks (canopy and ground) for example in canopy height estimation and focus more on the full geometry (raw shape) and radiometry (raw power) of the LiDAR waveform to retain richness of the data.

The moment distance is a matrix of distances computed from two reference locations (pivots) to each spectral or waveform point within the specified range.

Assume that a curve (reflectance or absorption curve or backscattered waveform) is displayed in Cartesian coordinates with the abscissa displaying the wavelength λ or time lapse t and the ordinate displaying the reflectance ρ or the backscattered power p. Let the subscript LP denote the left pivot (located in a shorter wavelength for the spectral curve and earlier temporal reference point for the waveform) and subscript RP denote the right pivot (located in a longer wavelength for the spectral curve and later temporal reference point for the waveform). Let λLP and λRP be the wavelength locations observed at the left and right pivots for a reflectance data, respectively, where left (right) indicates a shorter (longer) wavelength. Let tLP and tRP be the time value observed at the left and right pivots for a waveform data, respectively, where left (right) indicates an earlier (later) time. The proposed MD approach can be described in a set of equations.

Moment distance index (MDI) schematic diagram for spectral reflectance curve (top) and LiDAR waveform (bottom). Note the number of points between pivots can vary. Moment Distance (MD) Framework.png
Moment distance index (MDI) schematic diagram for spectral reflectance curve (top) and LiDAR waveform (bottom). Note the number of points between pivots can vary.

For spectral data, the index is given as:

where is the moment distance from the right pivot, is the moment distance from the left pivot, is the wavelength location at left pivot, is the wavelength location at right pivot, is the spectral reflectance at a given wavelength, and is successive wavelength location.

For waveform LiDAR data, the index is given as:

where the moment distance from the left pivot (MDLP) is the sum of the hypotenuses constructed from the left pivot to the power at successively later times (index from tLP to tRP): one base of each triangle is difference from the left pivot (tLP) along the abscissa and the other base is simply the backscattered power at . Similarly, the moment distance from the right pivot (MDRP) is the sum of the hypotenuses constructed from the right pivot to the power at successively earlier times (index from tRP to tLP): one base of each triangle is the difference from the right pivot (tRP) along the abscissa and the other base is simply the backscattered power at .

MDI is an unbounded metric. It increases or decreases as a nontrivial function of the number of spectral bands or bins considered and the shape of the spectrum or waveform that spans those contiguous bands or bins. The number of bands or bins is a function of the spectral resolution of the imaging spectrometer or the temporal resolution of the LiDAR (digitization rate) and the length of the reference range (i.e., full extent or subsets of the curve) being analyzed.




Related Research Articles

<span class="mw-page-title-main">Refractive index</span> Ratio of the speed of light in vacuum to that in the medium

In optics, the refractive index of an optical medium is a dimensionless number that gives the indication of the light bending ability of that medium.

<span class="mw-page-title-main">Wavelength</span> Distance over which a waves shape repeats

In physics and mathematics, wavelength or spatial period of a wave or periodic function is the distance over which the wave's shape repeats. In other words, it is the distance between consecutive corresponding points of the same phase on the wave, such as two adjacent crests, troughs, or zero crossings. Wavelength is a characteristic of both traveling waves and standing waves, as well as other spatial wave patterns. The inverse of the wavelength is called the spatial frequency. Wavelength is commonly designated by the Greek letter lambda (λ). The term "wavelength" is also sometimes applied to modulated waves, and to the sinusoidal envelopes of modulated waves or waves formed by interference of several sinusoids.

<span class="mw-page-title-main">Optical depth</span> Physics concept

In physics, optical depth or optical thickness is the natural logarithm of the ratio of incident to transmitted radiant power through a material. Thus, the larger the optical depth, the smaller the amount of transmitted radiant power through the material. Spectral optical depth or spectral optical thickness is the natural logarithm of the ratio of incident to transmitted spectral radiant power through a material. Optical depth is dimensionless, and in particular is not a length, though it is a monotonically increasing function of optical path length, and approaches zero as the path length approaches zero. The use of the term "optical density" for optical depth is discouraged.

In physics, coherence length is the propagation distance over which a coherent wave maintains a specified degree of coherence. Wave interference is strong when the paths taken by all of the interfering waves differ by less than the coherence length. A wave with a longer coherence length is closer to a perfect sinusoidal wave. Coherence length is important in holography and telecommunications engineering.

<span class="mw-page-title-main">Reflectance</span> Capacity of an object to reflect light

The reflectance of the surface of a material is its effectiveness in reflecting radiant energy. It is the fraction of incident electromagnetic power that is reflected at the boundary. Reflectance is a component of the response of the electronic structure of the material to the electromagnetic field of light, and is in general a function of the frequency, or wavelength, of the light, its polarization, and the angle of incidence. The dependence of reflectance on the wavelength is called a reflectance spectrum or spectral reflectance curve.

In mathematics, the spectral radius of a square matrix is the maximum of the absolute values of its eigenvalues. More generally, the spectral radius of a bounded linear operator is the supremum of the absolute values of the elements of its spectrum. The spectral radius is often denoted by ρ(·).

The laser diode rate equations model the electrical and optical performance of a laser diode. This system of ordinary differential equations relates the number or density of photons and charge carriers (electrons) in the device to the injection current and to device and material parameters such as carrier lifetime, photon lifetime, and the optical gain.

<span class="mw-page-title-main">Covariant transformation</span>

In physics, a covariant transformation is a rule that specifies how certain entities, such as vectors or tensors, change under a change of basis. The transformation that describes the new basis vectors as a linear combination of the old basis vectors is defined as a covariant transformation. Conventionally, indices identifying the basis vectors are placed as lower indices and so are all entities that transform in the same way. The inverse of a covariant transformation is a contravariant transformation. Whenever a vector should be invariant under a change of basis, that is to say it should represent the same geometrical or physical object having the same magnitude and direction as before, its components must transform according to the contravariant rule. Conventionally, indices identifying the components of a vector are placed as upper indices and so are all indices of entities that transform in the same way. The sum over pairwise matching indices of a product with the same lower and upper indices are invariant under a transformation.

In radiometry, irradiance is the radiant flux received by a surface per unit area. The SI unit of irradiance is the watt per square metre (W⋅m−2). The CGS unit erg per square centimetre per second (erg⋅cm−2⋅s−1) is often used in astronomy. Irradiance is often called intensity, but this term is avoided in radiometry where such usage leads to confusion with radiant intensity. In astrophysics, irradiance is called radiant flux.

<span class="mw-page-title-main">Kirchhoff's law of thermal radiation</span> Law of wavelength-specific emission and absorption

In heat transfer, Kirchhoff's law of thermal radiation refers to wavelength-specific radiative emission and absorption by a material body in thermodynamic equilibrium, including radiative exchange equilibrium. It is a special case of Onsager reciprocal relations as a consequence of the time reversibility of microscopic dynamics, also known as microscopic reversibility.

<span class="mw-page-title-main">CIE 1931 color space</span> Color space defined by the CIE in 1931

The CIE 1931 color spaces are the first defined quantitative links between distributions of wavelengths in the electromagnetic visible spectrum, and physiologically perceived colors in human color vision. The mathematical relationships that define these color spaces are essential tools for color management, important when dealing with color inks, illuminated displays, and recording devices such as digital cameras. The system was designed in 1931 by the "Commission Internationale de l'éclairage", known in English as the International Commission on Illumination.

In astronomy, the spectral index of a source is a measure of the dependence of radiative flux density on frequency. Given frequency in Hz and radiative flux density in Jy, the spectral index is given implicitly by

<span class="mw-page-title-main">Young stellar object</span> Star in its early stage of evolution

Young stellar object (YSO) denotes a star in its early stage of evolution. This class consists of two groups of objects: protostars and pre-main-sequence stars.

The refractive index of water at 20 °C for visible light is 1.33. The refractive index of normal ice is 1.31. In general, an index of refraction is a complex number with real and imaginary parts, where the latter indicates the strength of absorption loss at a particular wavelength. In the visible part of the electromagnetic spectrum, the imaginary part of the refractive index is very small. However, water and ice absorb in infrared and close the infrared atmospheric window thereby contributing to the greenhouse effect ...

The linear attenuation coefficient, attenuation coefficient, or narrow-beam attenuation coefficient characterizes how easily a volume of material can be penetrated by a beam of light, sound, particles, or other energy or matter. A coefficient value that is large represents a beam becoming 'attenuated' as it passes through a given medium, while a small value represents that the medium had little effect on loss. The (derived) SI unit of attenuation coefficient is the reciprocal metre (m−1). Extinction coefficient is another term for this quantity, often used in meteorology and climatology. Most commonly, the quantity measures the exponential decay of intensity, that is, the value of downward e-folding distance of the original intensity as the energy of the intensity passes through a unit thickness of material, so that an attenuation coefficient of 1 m−1 means that after passing through 1 metre, the radiation will be reduced by a factor of e, and for material with a coefficient of 2 m−1, it will be reduced twice by e, or e2. Other measures may use a different factor than e, such as the decadic attenuation coefficient below. The broad-beam attenuation coefficient counts forward-scattered radiation as transmitted rather than attenuated, and is more applicable to radiation shielding. The mass attenuation coefficient is the attenuation coefficient normalized by the density of the material.

Free spectral range (FSR) is the spacing in optical frequency or wavelength between two successive reflected or transmitted optical intensity maxima or minima of an interferometer or diffractive optical element.

In mathematics, the spectral theory of ordinary differential equations is the part of spectral theory concerned with the determination of the spectrum and eigenfunction expansion associated with a linear ordinary differential equation. In his dissertation, Hermann Weyl generalized the classical Sturm–Liouville theory on a finite closed interval to second order differential operators with singularities at the endpoints of the interval, possibly semi-infinite or infinite. Unlike the classical case, the spectrum may no longer consist of just a countable set of eigenvalues, but may also contain a continuous part. In this case the eigenfunction expansion involves an integral over the continuous part with respect to a spectral measure, given by the Titchmarsh–Kodaira formula. The theory was put in its final simplified form for singular differential equations of even degree by Kodaira and others, using von Neumann's spectral theorem. It has had important applications in quantum mechanics, operator theory and harmonic analysis on semisimple Lie groups.

<span class="mw-page-title-main">Thin-film interference</span> Optical phenomenon

Thin-film interference is a natural phenomenon in which light waves reflected by the upper and lower boundaries of a thin film interfere with one another, increasing reflection at some wavelengths and decreasing it at others. When white light is incident on a thin film, this effect produces colorful reflections.

In optics, the Fraunhofer diffraction equation is used to model the diffraction of waves when the diffraction pattern is viewed at a long distance from the diffracting object, and also when it is viewed at the focal plane of an imaging lens.

<span class="mw-page-title-main">Dirac equation in curved spacetime</span> Generalization of the Dirac equation

In mathematical physics, the Dirac equation in curved spacetime is a generalization of the Dirac equation from flat spacetime to curved spacetime, a general Lorentzian manifold.

References

  1. Salas, Eric Ariel L.; Henebry, Geoffrey M. (2013-12-19). "A New Approach for the Analysis of Hyperspectral Data: Theory and Sensitivity Analysis of the Moment Distance Method". Remote Sensing. 6 (1): 20–41. Bibcode:2013RemS....6...20S. doi: 10.3390/rs6010020 .
  2. Salas, Eric Ariel L.; Henebry, Geoffrey M. (2012-01-01). "Separability of maize and soybean in the spectral regions of chlorophyll and carotenoids using the Moment Distance Index". Israel Journal of Plant Sciences. 60 (1–2): 65–76. doi:10.1560/IJPS.60.1-2.65 (inactive 2024-01-27). ISSN   0792-9978.{{cite journal}}: CS1 maint: DOI inactive as of January 2024 (link)
  3. Aguilar, Manuel A.; Nemmaoui, Abderrahim; Novelli, Antonio; Aguilar, Fernando J.; García Lorca, Andrés (2016-06-18). "Object-Based Greenhouse Mapping Using Very High Resolution Satellite Data and Landsat 8 Time Series". Remote Sensing. 8 (6): 513. Bibcode:2016RemS....8..513A. doi: 10.3390/rs8060513 . hdl: 10835/6421 .
  4. Zhang, Peng; Du, Peijun; Guo, Shanchuan; Zhang, Wei; Tang, Pengfei; Chen, Jike; Zheng, Hongrui (July 2022). "A novel index for robust and large-scale mapping of plastic greenhouse from Sentinel-2 images". Remote Sensing of Environment. 276: 113042. doi: 10.1016/j.rse.2022.113042 .
  5. Nemmaoui, Abderrahim; Aguilar, Manuel A.; Aguilar, Fernando J.; Novelli, Antonio; García Lorca, Andrés (2018-11-06). "Greenhouse Crop Identification from Multi-Temporal Multi-Sensor Satellite Imagery Using Object-Based Approach: A Case Study from Almería (Spain)". Remote Sensing. 10 (11): 1751. Bibcode:2018RemS...10.1751N. doi: 10.3390/rs10111751 . hdl: 10835/6417 .
  6. Salas, Eric Ariel; Henebry, Geoffrey (2016-11-24). "Canopy Height Estimation by Characterizing Waveform LiDAR Geometry Based on Shape-Distance Metric". AIMS Geosciences. 2 (4): 366–390. Bibcode:2016AIMSG...2..366S. doi: 10.3934/geosci.2016.4.366 .
  7. Dutta, D.; Das, P. K.; Alam, K. A.; Safwan, P.; Paul, S.; Nanda, M. K.; Dadhwal, V. K. (September 2016). "Delta Area at Near Infrared Region ($textDA_textNIR$) #x2014;A Novel Approach for Green Vegetation Fraction Estimation using Field Hyperspectral Data". IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing. 9 (9): 3970–3981. doi:10.1109/JSTARS.2016.2539359. ISSN   1939-1404. S2CID   31461877.
  8. Barra, I.; Haefele, S. M.; Sakrabani, R.; Kebede, F. (June 2021). "Optimizing setup of scan number in FTIR spectroscopy using the moment distance index and PLS regression: application to soil spectroscopy". Scientific Reports. 11 (1): 13358. doi:10.1038/s41598-021-92858-w. PMC   8233441 . PMID   34172802.