THEMIS deconvolution

Contents: Description, Procedure, Functions Used, Related Functions Description This procedure illustrates the method to produce a set of atmospherically corrected TES emissivity spectra to atmospherically correct THEMIS data using rematm2(). Procedure Load your THEMIS image: dv> thm_spec = load_pds("I10418002.atmcorr.cub") Load your THEMIS spectral library (downloaded from the ASU Spectral Library Tool output) with THEMIS 10-point spectra as the output. dv> thm_lib = read("your_spec_lib.hdf") Remove band 10 from library: dv> thm_lib.data = thm_lib.data[,,:9] dv> thm_lib.xaxis = thm_lib.xaxis[,,:9] Ideally, band 1 and 2 should be averaged in both library and THEMIS image. Now run deconvolution using SMA: dv> out = sma(thm_spec.surf[x1:x2,y1:y2],thm_lib,nn = 1,noerr = 1)

DavinciWiki Mini-Nav Bar Contents Functions Function Categories Common Procedures Glossary Download Davinci Change Log News Mailing List Contact Developers davinci-dev [AT] mars.asu.edu All other topics See navigation on the left Functions Used sma read Related Procedures Standard THEMIS IR Processing with ISIS3 Using TES to atmospherically correct THEMIS