Emulation as an Accurate Alternative to Interpolation in Sampling Radiative Transfer Codes
This work offers a more accurate and efficient method for reconstructing radiative transfer model spectral data, which is crucial for researchers and practitioners in remote sensing and atmospheric science who rely on these models.
This paper investigates using emulation as an alternative to interpolation for sampling radiative transfer codes, which are computationally expensive. The authors found that emulation methods consistently produced more accurate output spectra than classical interpolation methods, with Gaussian Process Regression (GPR) emulation performing up to ten times more accurately than the best interpolation method at a competitive speed.
Computationally expensive Radiative Transfer Models (RTMs) are widely used} to realistically reproduce the light interaction with the Earth surface and atmosphere. Because these models take long processing time, the common practice is to first generate a sparse look-up table (LUT) and then make use of interpolation methods to sample the multi-dimensional LUT input variable space. However, the question arise whether common interpolation methods perform most accurate. As an alternative to interpolation, this work proposes to use emulation, i.e., approximating the RTM output by means of statistical learning. Two experiments were conducted to assess the accuracy in delivering spectral outputs using interpolation and emulation: (1) at canopy level, using PROSAIL; and (2) at top-of-atmosphere level, using MODTRAN. Various interpolation (nearest-neighbour, inverse distance weighting, piece-wice linear) and emulation (Gaussian process regression (GPR), kernel ridge regression, neural networks) methods were evaluated against a dense reference LUT. In all experiments, the emulation methods clearly produced more accurate output spectra than classical interpolation methods. GPR emulation performed up to ten times more accurately than the best performing interpolation method, and this with a speed that is competitive with the faster interpolation methods. It is concluded that emulation can function as a fast and more accurate alternative to commonly used interpolation methods for reconstructing RTM spectral data.