Standard THEMIS IR Processing with ISIS3
Line 13: | Line 13: | ||
− | Remove across track temperature gradients using [[rtilt|Rtilt]]. This step could be skipped if the cube is being processed for spectral analysis, but is highly recommended for | + | Remove across track temperature gradients using [[rtilt|Rtilt]]. This step could be skipped if the cube is being processed for highly rigorous spectral analysis, but is highly recommended for typical spectral analysis, DCS, and aesthetics. See Edwards et al. (2011) mosaicking paper for details. |
dv> a.cube = rtilt(a.cube) | dv> a.cube = rtilt(a.cube) | ||
Revision as of 12:39, 16 April 2012
Contents: Description, Procedure, Functions Used, Related Functions
Description This procedure is designed to take an ISIS projected THEMIS Infrared image and run it through "standard processing", correcting for noises and producing an emissivity array. NOTE: The example is shown for data containing 10 bands. If less bands are present, such as in a nighttime image, you can safely skip the autoradcorr step.
Procedure Run isis3setup(). dv> isis3setup()
Undrift and Dewobble the unprojected RDR data. There are several options to uddw that require access to the ASU Mars Space Flight Facility Database. Http = 1 is recommended for general use. If the http option is selected the file will be downloaded for you and all database calls are now made with the themis3db() function. dv>a = uddw("I10418002",isis3struct = 1, http = 1)
Remove across track temperature gradients using Rtilt. This step could be skipped if the cube is being processed for highly rigorous spectral analysis, but is highly recommended for typical spectral analysis, DCS, and aesthetics. See Edwards et al. (2011) mosaicking paper for details. dv> a.cube = rtilt(a.cube)
Write out the corrected ISIS cube dv> write(a,"I10418002.isis3.uddw.rtilt.cub",type = isis3, force = 1)
Outside of davinci (in the terminal), attach SPICE data to corrected cube:
Back in davinci, run isis3setup() again. dv> isis3setup()
Read in the ISIS map projected 10 band THEMIS infrared image. dv> a = load_pds("I10418002.isis3.uddw.rtilt.sinu.cub")
dv> a.cube[where min(a.cube,axis = z)< -2] = -32768
dv> b = thm.rectify(a.cube,ignore = -32768)
dv> b.data = thm.deplaid(b.data, ignore = -32768)
dv> rc = autoradcorr(b.data, ignore = -32768)
dv> b.data[where b.data != -32768] = b.data-rc.darkcol
dv> b.data = destreak(b.data, ignore = -32768)
dv> c = rrmnoise_pca(b.data, ignore = -32768)
dv> a.cube = thm.reconstitute(b)
dv> write(a,"I10418002.cub",isis3,force = 1)
dv> a.cube = c dv> write(a,"I10418002.rmnoise.cub",isis3,force = 1)
|
DavinciWiki Mini-Nav Bar Contents
Contact Developers
All other topics
Functions Used
|