Defines functions temp2raw

Documented in temp2raw

#' @export
temp2raw <-
  # This function is complementary to the raw2temp command,
  # except it converts a temperature estimate (in deg C) into a raw value
  # The temperature estimate would already have been obtained during analysis 
  # from a thermal imaging program under known atmospheric conditions
  # Temp should be supplied in degrees Celsius
  # how to call this function: 
  # temp2raw(temp,E,OD,RTemp,ATemp,IRWTemp,IRT,RH,PR1,PB,PF,PO,PR2)
  # example with all settings at default/blackbody levels
  # temp2raw(25,1,0,20,20,20,1,50,PR1,PB,PF,PO,PR2)
  # example with emissivity=0.95, distance=1m, window transmission=0.96, all temperatures=20C, 50% relative hum
  # temp2raw(20,0.95,1,20,20,20,0.96,50)  # default calibration constants for my FLIR camera will be used
  # Note: PR1, PR2, PB, PF, and PO are specific to each camera and result from the calibration at factory
  # of the camera's Raw data signal recording from a blackbody radiation source
  # Calibration Constants                 (Glenn's FLIR, Ray's T300(25o) Ray's, T300(telephoto), Glenn's Mikron )
  # PR1: PlanckR1 calibration constant from FLIR file  21106.77       14364.633     14906.216       21106.77
  # PB: PlanckB calibration constant from FLIR file    1501           1385.4        1396.5          9758.743281
  # PF: PlanckF calibration constant from FLIR file    1              1             1               29.37648768
  # PO: PlanckO calibration constant from FLIR file    -7340          -5753         -7261           1278.907078
  # PR2: PlanckR2 calibration constant form FLIR file  0.012545258    0.010603162   0.010956882     0.0376637583528285
  # PR1=14906.216,PB=1396.5,PF=1,PO=-7261,PR2=0.010956882       
  #  Uncomment this line above to reveal calibration values for Ray's T300 camera
  # PR1=21106.77,PB=9758.743281,PF=29.37648768,PO=1278.907078,PR2=0.0376637583528285
  #  Uncomment this line above to reveal calibration values for Glenn's Mikron camera
  # Set constants below. Comment out those that are variables in this function
  # Keep those that should remain constants.  These are here to make troubleshooting calculations easier
  # temp<-25; E<-0.95; OD<-20; IRT<-0.96
  # RTemp<-20; IRWTemp<-RTemp; ATemp<-20; RH<-50
  # PR1<-21106.77; PB<-1501; PF<-1; PO<--7340; PR2<-0.012545258
  # Keep these humidity parameters in since they are set constants required for a calculation
  ATA1<-0.006569; ATA2<-0.01262; ATB1<--0.002276; ATB2<--0.00667; ATX<-1.9
  # Equations to convert to temperature
  # See,4898.60.html
  # Standard equation: temperature<-PB/log(PR1/(PR2*(raw+PO))+PF)-273.15
  # Other source of information: Minkina and Dudzik's Infrared Thermography: Errors and Uncertainties
  refl.wind<-0 # anti-reflective coating on window
  # converts relative humidity into water vapour pressure (I think in units mmHg)
  # transmission through atmosphere - equations from Minkina and Dudzik's Infrared Thermography Book
  raw.refl1<-PR1/(PR2*(exp(PB/(RTemp+273.15))-PF))-PO   # radiance reflecting off the object before the window
  raw.refl1.attn<-(1-E)/E*raw.refl1   # attn = the attenuated radiance (in raw units) 
  raw.atm1<-PR1/(PR2*(exp(PB/(ATemp+273.15))-PF))-PO # radiance from the atmosphere (before the window)
  raw.atm1.attn<-(1-tau1)/E/tau1*raw.atm1 # attn = the attenuated radiance (in raw units) 

Try the Thermimage package in your browser

Any scripts or data that you put into this service are public.

Thermimage documentation built on Nov. 30, 2017, 9:01 a.m.