emcal2

Description

Emissivity Calculation of Laboratory Spectra at ASU/SWRI

Computes an instrument response function from measured emission, spectra of hot and warm blackbody and known blackbody temperatures. Then calculates emissivity using previously calculated noise free response function and instrument radiance, hot bb spectrum, hot sample spectrum, and known environment temperature.

Arguments and Return Values

Arguments: A directory containing spectra files (*.CSV), *hot*.CSV (hot blackbody) and *warm*.CSV (warm blackbody)

Return Value: Writes either a single .hdf file or if individual is selected individual .txt files to be used with readvm()

Usage

Syntax: emcal2(dir, output_filename, warmBBT/warmBBRes, hotBBT/hotBBRes [, env_temp_array] [,lab="asu/swri"] [, noise_free=BOOL] [, dehyd=BOOL] [, individual=BOOL] [, resample=BOOL][,warmbb=STRING] [,hotbb=STRING] [,ratw=BOOL])

dir - A string defining the directory containing appropriate spectra

output_filename - A string containing the desired name of the output file

warmbbT/warmBBRes - Temperature of hot blackbody (from r2t_lo() or r2t_swri()) or Resistance values for each channel (Ch103//Ch104 or Ch1lo//Ch2lo//Ch1hi//Ch2hi)

hotbbT/hotBBRes - Temperature of warm blackbody (from r2t_hi() or r2t_swri()) or Resistance values for each channel (Ch103//Ch104 or Ch1lo//Ch2lo//Ch1hi//Ch2hi)

env_temp_array - temperature array containing the measured environmental temperature for each sample (sorted in alphabetical order)

lab - defines the lab where the spectrum was acquired (Default = "asu" or "swri"). This is only used if resistance values are input as it defines the use of r2t_lo()/rt2_hi() vs rt2_swri()

noise_free - computes a "noise free/reduced" instrument response function and instrument radiance. This helps increase the overall "Science" to Noise Ratio

dehyd - add dehydrated spectra to the output structure (Uses dehyd2())

individual - write spectra as individual text files (Uses write_vm())

resample - resample spectra to the response function xaxis (Uses resample())

warmbb - warm blackbody partial filename search (Default search ="warm")

hotbb - hot blackbody partial filename (Default search ="hot")

ratw - calibrate spectra using a "Rock Around the World" standard procedure

Examples

dv> out=emcal2("RATW_emcal2/","emcal2",1279.53//1281.20,1391.76//1390.63)

 Beginning manual environmental temperature entry
 Enter "q" to quit at any time

 Enter the Chamber Temperature (Ch105 in degrees C) for
 Sample: Kuiper_1
 25.6

 Enter the Chamber Temperature (Ch105 in degrees C) for
 Sample: Kuiper_2
 25.3

 Enter the Chamber Temperature (Ch105 in degrees C) for
 Sample: Kuiper_3
 25.5

 Enter the Chamber Temperature (Ch105 in degrees C) for
 Sample: Painter_big_black
 25.3

 Enter the Chamber Temperature (Ch105 in degrees C) for
 Sample: Painter_medium_black
 26.1
 ..................
 struct, 9 elements
   xaxis: 1x1x1971 array of float, bsq format [7,884 bytes]
   raw: struct, 15 elements
       ch1res: 1279.53
       ch2res: 1281.20
       ch3res: 1391.76
       ch4res: 1390.63
       warmbbt: 71.4075
       hotbbt: 100.268
       hdfname: "emcal2.hdf"
       inputdir: "RATW_emcal2/"
       speclist: Text Buffer with 20 lines of text
           1: RATW_emcal2/Kuiper_1.CSV
           2: RATW_emcal2/Kuiper_2.CSV
           3: RATW_emcal2/Kuiper_3.CSV
           4: RATW_emcal2/Painter_big_black.CSV
           5: RATW_emcal2/Painter_medium_black.CSV
           6: RATW_emcal2/Painter_red.CSV
           7: RATW_emcal2/Painter_small_black.CSV
           8: RATW_emcal2/RATW11066.CSV
           9: RATW_emcal2/RATW11067.CSV
           10: RATW_emcal2/RATW11069.CSV
       hotbbfile: "RATW_emcal2/bbhot.CSV"
       warmbbfile: "RATW_emcal2/bbwarm.CSV"
       hotbb: 1x1x1971 array of float, bsq format [7,884 bytes]
       warmbb: 1x1x1971 array of float, bsq format [7,884 bytes]
       data: 1x20x1971 array of float, bsq format [157,680 bytes]
       envtemp: 1x20x1 array of float, bsq format [80 bytes]
   irf_ri: struct, 12 elements
       xaxis: 1x1x1971 array of float, bsq format [7,884 bytes]
       inst: "spectrometer"
       t1: 373.418
       t2: 344.557
       irf: 1x1x1971 array of double, bsq format [15,768 bytes]
       ri: 1x1x1971 array of float, bsq format [7,884 bytes]
       irf1: 1x1x1971 array of double, bsq format [15,768 bytes]
       irf2: 1x1x1971 array of double, bsq format [15,768 bytes]
       noise_free_inst_temp: 298.954
       noise_free_method: "fit_bb"
       riw1: 300
       riw2: 2500
   cal_rad: struct, 3 elements
       xaxis: 1x1x1971 array of float, bsq format [7,884 bytes]
       inst: "spectrometer"
       radiance: 1x20x1971 array of double, bsq format [315,360 bytes]
   emissivity: struct, 27 elements
       xaxis: 1x1x1971 array of float, bsq format [7,884 bytes]
       dehyd: 0
       wave1: 500
       wave2: 1700
       w1_cold: 450
       w2_cold: 900
       threshold_t_cold: 160.000
       threshold_t_warm: 170.000
       filter: 3
       co2chan: 0
       co2lowchan: 0
       co2highchan: 0
       tb: 1x20x1971 array of double, bsq format [315,360 bytes]
       tfilt1: 1x20x1971 array of float, bsq format [157,680 bytes]
       tfilt2: 1x20x1971 array of float, bsq format [157,680 bytes]
       max_t1: 1x20x1 array of float, bsq format [80 bytes]
       max_t2: 1x20x1 array of float, bsq format [80 bytes]
       pos_t1: 1x20x1 array of int, bsq format [80 bytes]
       pos_t2: 1x20x1 array of int, bsq format [80 bytes]
       rad_bb: 1x20x1971 array of double, bsq format [315,360 bytes]
       max_emiss_low: 1.00000
       downwelling_e: 1.00000
       downwelling_t: 1x20x1 array of float, bsq format [80 bytes]
       downwelling_rad: 1x20x1971 array of double, bsq format [315,360 bytes]
       max_emiss: 1.00000
       emissivity: 1x20x1971 array of double, bsq format [315,360 bytes]
       target_temperature: 1x20x1 array of float, bsq format [80 bytes]
   radiance: 1x20x1971 array of double, bsq format [315,360 bytes]
   sample_temperature: 1x20x1 array of float, bsq format [80 bytes]
   data: 1x20x1971 array of double, bsq format [315,360 bytes]
   label: Text Buffer with 20 lines of text
       1: Kuiper_1
       2: Kuiper_2
       3: Kuiper_3
       4: Painter_big_black
       5: Painter_medium_black
       6: Painter_red
       7: Painter_small_black
       8: RATW11066
       9: RATW11067
       10: RATW11069

Images

The left spectrum is an example response function while the right spectrum is an example emissivity spectrum.