Standard THEMIS IR Processing
|
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. The data output of this should be identical to what can be generated using thmproc. 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
dv>a = read("I10418002.irmc.cub")
..........577x1852x10 array of float, bsq format [42,744,160 bytes]
Rectify the image, which rotates and slightly warps the projected image into a rectangle for easy manipulation and to minimize storage requirements. dv>a = thm.rectify(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 plaid. dv> a.data = thm.deplaid(a.data) 301x1883x10 array of float, bsq format [22,671,320 bytes] Call the automatic radiance correction function to remove the atmospheric radiance from surface data providing 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] Subtract the column of correction values from the data where the array is not the ignore value. dv>a.data[where a.data != -32768] = 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
All other topics
Functions Used
Related Procedures |
