THEMIS deconvolution
(One intermediate revision by one user not shown) | |||
Line 1: | Line 1: | ||
− | {{DavinciWiki_Procedure|This procedure illustrates the method to | + | {{DavinciWiki_Procedure|This procedure illustrates the method to linearly unmix THEMIS spectrum using a spectral endmember library. It is highly advised that you follow the [[Standard_THEMIS_IR_Processing_with_ISIS3|THEMIS standard processing]] procedure, followed by [[rematm2|THEMIS atmospheric correction]] before attempting THEMIS deconvolution. You may also obtain a processed THEMIS image from [http://thmproc.mars.asu.edu ThmProc], but the atmospheric correction is (for the time being) up to you! |
| | | | ||
Load your THEMIS image: | Load your THEMIS image: | ||
Line 5: | Line 5: | ||
− | Load your THEMIS spectral library (downloaded from the ASU [[Spectral Library Tool]] output) with THEMIS 10-point spectra as the output. | + | Load your THEMIS spectral library (downloaded from the ASU [[Spectral Library Tool]] output) with THEMIS 10-point spectra as the output. You may also use the [[i2i]] function to resample an preexisting endmember library at a different instrument sampling (i.e., TES) to THEMIS spectral resolution. |
dv> thm_lib = read("your_spec_lib.hdf") | dv> thm_lib = read("your_spec_lib.hdf") | ||
Line 14: | Line 14: | ||
− | Ideally, band 1 and 2 should be averaged in both library and THEMIS image. | + | Ideally, band 1 and 2 should be averaged in both library and THEMIS image - or, just run SMA over bands 2-9, as below. |
− | + | ||
Now run deconvolution using SMA: | Now run deconvolution using SMA: | ||
− | dv> out = sma(thm_spec.surf[x1:x2,y1:y2],thm_lib,nn = 1,noerr = 1) | + | dv> out = sma(thm_spec.surf[x1:x2,y1:y2],thm_lib,nn = 1,noerr = 1, wave1 = 2, wave2 = 9, band = 1) |
| | | | ||
*[[sma]] | *[[sma]] |
Latest revision as of 12:06, 20 June 2014
Contents: Description, Procedure, Functions Used, Related Functions
Description This procedure illustrates the method to linearly unmix THEMIS spectrum using a spectral endmember library. It is highly advised that you follow the THEMIS standard processing procedure, followed by THEMIS atmospheric correction before attempting THEMIS deconvolution. You may also obtain a processed THEMIS image from ThmProc, but the atmospheric correction is (for the time being) up to you!
Procedure Load your THEMIS image: dv> thm_spec = load_pds("I10418002.atmcorr.cub")
dv> thm_lib = read("your_spec_lib.hdf")
dv> thm_lib.data = thm_lib.data[,,:9] dv> thm_lib.xaxis = thm_lib.xaxis[,,:9]
dv> out = sma(thm_spec.surf[x1:x2,y1:y2],thm_lib,nn = 1,noerr = 1, wave1 = 2, wave2 = 9, band = 1) |
DavinciWiki Mini-Nav Bar Contents
Contact Developers
All other topics
|