EOS 840 Hyperspectral Imaging Applications

1
EOS 840 Hyperspectral Imaging Applications
October 27, 2004 Week 9
Ron Resmini v 703-735-3899 ronald.g.resmini_at_boein
g.com Office hours by appointment
2
Outline

• Review/context/a thread...
• A bit more RT (and one more item...)
• Finish algorithms (two more to go)
• Thermal infrared remote sensing
• Theory
• SEBASS
• Working with TIR in ENVI
• My semester project status
• Your semester project status

3
Review/Context/A Thread...
4
A Bit More RT...
5
A key point
Apply lots of simple algebra
6
One More Aside...

• I want to stress two points made by Dr. Barry
• 18 bands for application 1
• 18 bands for application 2
• etc... and
• Atmospheric compensation with HSI

7
Algorithms (continued)
8
Constrained Energy Minimization (CEM)

• The description of CEM is similar to that of
  OSP/DSR (previous slides)
• Like OSP and DSR, CEM is an Orthogonal Subspace
  Projection (OSP)family algorithm
• CEM differs from OSP/DSR in the following,
  important ways
• CEM does not simply project away the first n
  eigenvectors
• The CEM operator is built using a weighted
  combination of theeigenvectors (all or a subset)
• Though an OSP algorithm, the structure of CEM is
  equally readily observed bya formal derivation
  using a Lagrange multiplier

• CEM is a commonly used statistical spectral
  matched filter
• CEM for spectral remote sensing has been
  published on for over 10 years
• CEM has a much longer history in the
  multi-dimensional/array signalprocessing
  literature
• Just about all HSI tools today contain CEM or a
  variant of CEM
• If an algorithm is using M-1d as the heart of its
  filter kernel (where M is thedata covariance
  matrix and d is the spectrum of the target of
  interest), thenthat algorithm is simply a CEM
  variant

9

• The statistical spectral matched filter (SSMF)
• Derivation in detail
• Application of the filter
• Statistics
• Endmembers (FBA/MCEM)
• Interpretation of results
• Many algorithms are actually the basic SSMF
• Different ways to apply the filter/application
  strategies(i.e., in-scene spectra/library
  spectra)
• Matched filter in ENVI

10
Derivation taken from
Stocker, A.D., Reed, I.S., and Yu, X., (1990).
Multi-dimensional signal processing for
electro-Optical target detection. In Signal
and Data Processing of Small Targets 1990,
Proceedingsof the SPIE, v. 1305, pp. 218-231.
J of Bands
Form the log-likelihood ratio test of Hº and H1
11
Some algebra...
12
A trick...recast as a univariable problem
After lots of simple algebra applied to the r.h.s
Now, go back to matrix-vector notation
13
Take the natural log
14
Constrained Energy Minimization (CEM)
(Harsanyi et al., 1994)
15
(No Transcript)
16
Note...

• See also the Lagrange multiplier derivation
• Previous techniques exploit shape and albedo
• this can cause problems...
• Sub-classes of algorithms developed to mitigate
  this
• shape, only, operators
• MED, RSD of ASIT, Inc.
• MTMF of ENVI

17
Last Class of Algorithms

• Spectral feature fitting/derivativespectroscopy
• Spectral parameterizations
• Wavelets
• Band depth/band depth mapping
• Application strategies (i.e., in-scenespectra/lib
  rary spectra)
• Mixed pixels...

18
Another Note...

• Performance prediction/scoring/NP-Theory, etc...
• Hybrid techniques
• still some cream to be skimmed...
• Caveat emptor...
• lots of reproduction of work already accomplished
• who invented what? when?
• waste of resources
• please do your homework!read the lit.!

19
Thermal Infrared (TIR) Remote Sensing...
20
Reflective vs. Emissive
Reflective Band Emissive Band Physics
Reflected sunlight Direct thermal
emission Complex phenomena Simpler phenomena
Surface reflectivity Surface
temperature Illumination geometry Emissivity
Phenomena Familiar Familiar
21
Properties of MWIR/LWIR

• Day and Night Operations
• Material Identification/Quantification
• Gas Identification/Quantification
• Fingerprint Spectral Region
• Better T measurements (ideally)
• Physics-Based Algorithms
• True for all band regions

22
A Bit of Theory
The Planck or Blackbody Radiation Equation
23
An Energy Balance
Good absorbers are good emitters
Kirchhoffs Law
24
Stefan-Boltzmann Law
Two surfaces radiating at each other
What happens when T1 T2?
View Factor Algebra and Radiant Exchange...
25
BTW...M for an HSI band is
BTW...(again)...are we in irradiance or radiance?
26
Some More Interesting Stuff...
Wiens Displacement Law
A 2898 mm.K
27
The radiant exitance of the sun is
The total flux from the surface of the sun is
28
Some Values in the Previous Equations

• k Boltzmann Gas Constant 1.38 x 10-23 J/K
• s Stefan-Boltzmann Constant 5.67 x 10-8
  W/m2.K4
• c Speed of Light 2.9979 x 108 m/sec
• h Plancks Constant 6.6256 x 10-34 J.sec

29
The Basic TIR RT Expression
Just about all papers on TIR HSI will start with
this basic RT expression. Its not the only one to
use what about other terms from the Big
Equation?
30
Atmospheric Compensation
In-Scene Atmospheric Compensation (ISAC)
Young, S.J., Johnson, R.B., and Hackwell, J.A.,
(2002). An in-scene method for atmospheric
compensation of thermal hyperspectral data.
Journal of Geophysical Research, v. 107, no. D24,
4774, doi10.1029/2001JD001266, 20 p.
Start with
Let
Get
31
Atmospheric Compensation (continued)

• Calculate a brightness temperature for each
  spectrum
• Calculate Planck function radiance for
  that temperature
• Plot Ll vs. Planck radiance
• Interested in blackbodies thus e 1

Get

• Fit a line to tops of clusters where LS is large
  as is e

32
Atmospheric Compensation (continued)

• Note from linear equation that you get t and Lu
• Subtract Lu from original radiance spectra
• Divide by t
• Left with ground leaving radiance (GLR) spectra
• Must now apply temperature/emissivity
  separation(TES)
• Refer to Young et al., (2002) for more details

Other Theres also AAC and EELM
33
Temperature/Emissivity Separation (TES)
The Normalized Emissivity Method (NEM)

• On a pixel-by-pixel basis
• Find the maximum radiance value
• Assume e 0.97 (or some such value)
• Find T by inverting the Planck function
• Divide original GLR spectrum but thePlanck
  function just calculated

34
There are lots of TES routines.
See also
Kaiser, R. D. (1999), Quantitative comparison of
temperature / emissivity algorithm performance
using SEBASS data. SPIE Vol. 3717, pp. 47-57.
35
Exploitation of TIR HSI

• Atmospheric compensation/TES
• Otherwise, apply algos. already discussed!
• TIR HSI are also points in n-D hyperspace
• Radiance, GLR, emissivity, temp.
• Which one? When? Why?
• How you conceive of/think of your data...
• ENVI spectral library is in reflectance
• Apply Kirchhoffs Law (e 1 r)

36
(No Transcript)
37
The SEBASS TIR HSI Sensor...
SEBASS Radiance Units Are
38
TIR Exploitation with ENVI...
39
TIR HSI Applications

• Day/night remote sensing
• Gas phase remote sensing
• Solid material detection, identification, and
  quantification
• Precision thermometry T 0.01C
• Surveillance / security
• Geology / mineral mapping / volcanology
• Bowers and Resmini (2004)
• Others...many others!

40
MWIR HSI

• Reflected plus emittedduring the day
• MWIR at night is analogous to LWIR
• Reflected plus emittedvery complicated
• Not much more will be said now
• See the following excellent, recently
  publishedreview article

Boyd, D., and Petitcolin, F., (2004). Remote
sensing of the terrestrial environment using
middle infrared radiation (3.0-5.0 mm).
International Journal of Remote Sensing, v. 25,
no. 17, September, pp. 3343-3368.
41

• My semester project status
• slides from last week...
• Your semester project status?
• Got data?