Passive infrared optical gas imaging (IOGI) is sensitive to toxic or greenhouse gases of interest, offers non-invasive remote sensing, and provides the capability for spatially resolved measurements. It has been broadly applied to emission detection, localization, and visualization.
|
Molecule |
Absorption Wavebands (µm) |
Molecule |
Absorption Wavebands (µm) |
||
|---|---|---|---|---|---|
|
CO |
2.29–2.48 |
4.36–5.05 |
NO |
2.63–2.78 |
5.03–5.78 |
|
CO2 |
2.66–2.81 |
4.07–4.43 |
NO2 |
3.38–3.53 |
5.92–6.57 |
|
CH4 |
3.07–3.71 |
6.67–9.09 |
SO2 |
3.94–4.07 |
6.94–9.44 |

where I is the radiative intensity (W/(m2⋅sr)), the subscripts v,o,e and t represent the frequency of the light wave, the output of the layer, emission of the layer, and transmission of the layer, respectively.
According to the Beer–Lambert law, transmissivity can be expressed as:
where t is the transmissivity, kv is the absorption coefficient (cm−1), and l is the light path length of the layer (cm). Meanwhile, the sum of transmissivity and absorptivity is 1 when scattering and reflection are neglectable, that is:
where αv is the absorptivity. Supposing all three layers are in thermal equilibrium, and following Kirchhoff’s law of thermal radiation that emissivity equals absorptivity, the researchers substitute Equations (2) and (3) into Equation (1):
where Iv,B is the black body radiation, Iv,i is the radiation intensity at the incident surface of the layer. Blackbody radiation can be modeled by Plank’s law as:
where h is the Plank constant, kB is the Boltzmann constant, c is the light speed, and T is the temperature. Then, the absorption coefficient, kv, of gas species j could be calculated from:
where s is the line intensity per molecule (cm−1/(molecule⋅cm−2)), which is a function of temperature, ϕv is the line shape function [36] (cm−1), which is a function of both pressure and temperature, P is the local pressure (Pa), and X is a mole fraction. The total absorption coefficient can be simplified as the sum of the absorption coefficients of each species:
In practice, a bandpass filter is used to select the target gas species, so Equation (7) can be simplified as follows:
where subscript t represents the target gas, therefore, combining Equations (4)–(8), Equation (4) can be rewritten as:
The mole fraction Xt is coupled with the light path length l. Using the ideal gas equation of the state, Xtl, can be transformed as the product of concentration Ct and light path length l, which is called column density [25,37][25][37] and can be expressed as follows:
where CLt, Mt and Ct are the column density, molecular mass, and concentration of target gas, respectively. The final intensity at each pixel is a function of the camera characteristics and double integration of the camera incident intensity at all wavelengths of the band and the covered surface of the pixel, i.e.,
where Ip is the pixel value, f symbolizes a functional relationship, Iv,c is the camera incident light strength at the wavelength inside the filtered waveband, dA is a finite element of the surface that a pixel covers, and D
represents the camera characteristics of transforming radiation into pixel charge and includes chip-sensing efficiency, transforming efficiency, and noise characteristics.
In most cases, pressure can be supposed to be constant along the optical path. After accumulating the transmission and absorption of three layers, the pixel value can be represented as follows:
where Rb is the background radiation, CLg and CLf are the column densities of the target gas substance in gas cloud and foreground, respectively, Tf and Tg are the temperatures of the gas cloud and foreground, respectively, D is the device characteristics, and ε is the noise that comes from the environment and devices, such as wind effect, scattering, etc., which can be neglected and regarded as measurement uncertainty. Parameters Rb, CLf, Tf, D and ε can be summarized as environmental factor εe as they are controlled by the environment and measurement devices and are in general considered constant for the experimental condition. Thus, for a given imager, Equation (12) can be simplified and rewritten as:
Consequently, column density is the function of Ip,Tg, εe, that is:
Equations (13) and (14) reveal two important insights. First, from the pixel value of an IR image, we cannot decouple the concentration and light path length can not be decoupled—they are represented by column density. Second, the fundamental quantitative parameter, column density, is a function of three parameters, i.e., pixel value, gas temperature, and environmental factors (Rb, CLf, Tf, D and ε ). Therefore, a single pixel value is not sufficient to retrieve column density.
The concept of the augmentation method is to add spectral information to each pixel, that is, spectral imaging, a combination of spectroscopy and imaging, so that the image conveys temperature information as well, in a latent way. The problem then transforms to column density quantification from the spectral information at each pixel, i.e.,
where v is the frequency, Bg is the characteristic band for the target gas, and SP is the spectrum obtained at a given pixel, which is a wavenumber function.