thm.radcorr

From DavinciWiki
Jump to: navigation, search


Description

Calculate the radiance correction and emissivity offset values for a THEMIS radiance array.

Produces a 2x1x10 float array containing the radiance correction and em offset values for a given array. Values are minimized RMS differences of radiance from blackbody of same temperature. Uses inverse parabolic interpolation for minimization algorithm

This is the minimization function for autoradcorr and radcorr2 used in standard THEMIS processing.

Arguments and Return Values

Arguments: A small (100x100x1) THEMIS radiance array

Return Value: A 2x1xN float array, where N is the number of bands

Usage

Syntax: thm.radcorr(radiance = VAL [, bandlist = STRING] [, ignore = VAL] [, temp_rad_path = STRING] [, b1 = INT] [, b2 = INT] [, space = 1])

'radiance' - a small (50x50x10 pixels) array selected from the entire radiance array

'bandlist' - a list of THEMIS bands in radiance array. Default is ordered bands 1-10

'ignore' - the non-data value. Default is -32768

'temp_rad_path' - the path and filename of the temperature vs radiance lookup table

'b1' - the first band in the ordered list to search for the maximum brightness temperature. Default = 3

'b2' - the last band in the ordered list to search for the maximum brightness temperature. Default = 9

'space' - a binary flag to return a complete, bounded solution space

Examples

dv> a
345x500x10 array of float, bsq format [6,900,000 bytes]

dv> rc = thm.radcorr(a[20:120,450:500])
.2x1x10 array of double, bsq format [160 bytes]
-2.93331e-06    0.874220

-3.38604e-07    0.884298

2.01230e-06     1.00053

4.89527e-06     0.906034

3.45005e-06     0.852073

2.51699e-06     0.864984

1.87052e-06     0.899031

8.77423e-07     0.940794

4.50610e-06     0.926166

0.00000         0.00000

or, a more verbose example:

dv> rc = thm.radcorr(a[20:120,450:500],ignore=-32768,b1=4,b2=6)
.2x1x10 array of double, bsq format [160 bytes]
-4.34115e-06    0.983762

-1.71475e-06    0.994783

1.74100e-07     1.10584

1.88402e-06     0.998468

5.69619e-07     0.931107

-5.51093e-07    0.938657

-1.66467e-06    0.971046

-2.65994e-06    1.01136

1.08694e-06     0.991761

0.00000         0.00000

DavinciWiki Mini-Nav Bar

Contents


Contact Developers

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

All other topics

  • See navigation on the left

Related Functions

Recent Core Changes

Modified On: 04-15-2006

Personal tools