Standard THEMIS IR Processing

From DavinciWiki
Revision as of 23:23, 5 March 2006 by Cedwards (Talk | contribs)

Jump to: navigation, search

Description

This procedure is desinged to take an ISIS projected THEMIS Infrared image and run it through standard processing. The data output of this should be identical to what can be generated using thmproc.

NOTE: The example is shown for 10 bands, if less bands are present like in a nighttime image you can safely skip the autoradcorr step.

Procedure


Read in the ISIS map projected 10 band THEMIS infrared image.

dv>a=read("I10418002.irmc.cub")
..........577x1852x10 array of float, bsq format [42,744,160 bytes]

Rectify the image, which turns the projected image into a rectangle for easy manipulation.

dv>a=thm.rectfiy(a)
struct, 4 elements
   data: 301x1883x10 array of float, bsq format [22,671,320 bytes]
   leftedge: 1x1883x1 array of int, bsq format [7,532 bytes]
   width: 577
   angle: 3.052882433

Remove the line and row correlated noise.

dv> a.data=thm.deplaid(a.data)
301x1883x10 array of float, bsq format [22,671,320 bytes]

Perform an automatic radiance correction which allows for better separation of brightness temperature and emissivity.

dv> rc=autoradcorr(a.data)
..
struct, 4 elements
   avgdark: 1x1x10 array of float, bsq format [40 bytes]
   blackmask: 301x1883x1 array of byte, bsq format [566,783 bytes]
   dark: 7x38x10 array of float, bsq format [10,640 bytes]
   darkcol: 1x1883x10 array of float, bsq format [75,320 bytes]

To apply the correction the correction must be subtracted from the data.

dv>a.data=a.data-rc.darkcol

Return the data to its original map projection.

dv> a=thm.reconstitute(a)
577x1852x10 array of float, bsq format [42,744,160 bytes]


DavinciWiki Mini-Nav Bar

Contents


Contact Developers

  • davinci-dev [AT] mars.asu.edu

All other topics

  • See navigation on the left

Personal tools