Using TES to atmospherically correct THEMIS

Description

This procedure illustrates the method to produce a set of atmospherically corrected TES emissivity spectra to atmospherically correct THEMIS data using rematm2(). The TES data must overlap the THEMIS training area and should be downloaded from JMARS as a .txt file. When possible, select pixels from a single ock and from the same season and similar atmospheric opacity as your THEMIS image.

Procedure

Load your spectral library (downloaded from the ASU Spectral Library Tool output) with tes 73-point spectra as the output.

dv> lib = read("your_spec_lib.hdf")

Load the 73 channel atmospheric library, ether from the Spectral Library Tool or the example file below.

dv> atm = read($DV_EX+"/TES_atmlib_73channel.hdf")

Load your TES emissivity data and split the data by orbit using split_ock() so that an assessment can be made of the atmospheric conditions from each orbit.

dv> tes = split_ock("tesfile.txt")

Use sma() to model TES atmospherically corrected surface emissivity spectra.

dv> spec = sma(tes.####avgcat, lib, atm, wave1 = 200, wave2 = 1305, surface = 1, nn = 1)

Resample the TES atmospherically corrected data to the same spectral sampling as THEMIS.

dv> tes73x = make_band(tes73 = 1)
dv> surf_em = i2i(spec, from = tes73x, to = 'themis')

Add a discussion of rematm2