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.

NOTE: All temperatures are assumed to be input in C and are converted to K for use with emissivity(), cal_rad(), irf_ir().

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] [,max_emiss=FLOAT] [,ratw=BOOL] [, wave1=FLOAT] [, wave2=FLOAT])

or (if using the LabView data acquisition software

Syntax: emcal2(data_qub, output_filename [, noise_free=BOOL] [, dehyd=BOOL] [, individual=BOOL] [, resample=BOOL][,warmbb=STRING] [,hotbb=STRING] [,max_emiss=FLOAT] [,ratw=BOOL] [, wave1=FLOAT] [, wave2=FLOAT])

data_qub - Full path to output data qub from LabView data acquisition software.

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 in C (from r2t_lo() or r2t_swri()) or Resistance values for each channel (Ch103//Ch104(asu) or Ch1lo//Ch2lo//Ch1hi//Ch2hi(swri))

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

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 ="bbw")

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

max_emiss - optional maximum emissivity (e.g. 0.98) (default = 1.)

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

wave1 - optional start wavenumber for brightness temperature calculation (Default is 500)

wave2 - optional end wavenumber for brightness temperature calculation (Default is 1700)

Examples

Either option shown below will produce the same results.

dv> out=emcal2("spectra/","output",70.1232,100.031)

or

dv> out=emcal2("spectra/","output",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: K_1
 25.6

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

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

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

 Enter the Chamber Temperature (Ch105 in degrees C) for
 Sample: P_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: K_1
       2: K_2
       3: K_3
       4: P_big_black
       5: P_medium_black
       6: P_red
       7: P_small_black
       8: R11066
       9: R11067
       10: R11069

Images

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