The main goal of this paper is the estimation of thermal resistive fields in multilayer samples using the classical front face flash method as excitation and InfRared Thermography (IRT) as a monitoring sensor. The complete inverse processing of a multilayer analytical model can be difficult due to a lack of sensitivity to certain parameters (layer thickness, depth of thermal resistance, etc.) or processing time. For these reasons, our present strategy proposes a Bayesian inference approach. Using the analytical quadrupole method, a reference model can be calculated for a set of parameters. Then, the Bayesian probabilistic method is used to determine the maximum likelihood probability between the measured data and the reference model. To keep the processing method robust and fast, an automatic selection of the calculation range is proposed. Finally, in the case of a bilayer sample, both the thickness and resistive 3D layers are estimated in less than 2 min for a space and time matrix of 50000 pixels by 4000 time steps with a reasonable relative error of less than 5%.