Our investigation raises an important question that is of relevance to the wider turbomachinery community: how do we estimate the spatial average of a flow quantity given finite (and sparse) measurements? This paper seeks to advance efforts to answer this question rigorously. In this paper, we develop a regularized multivariate linear regression framework for studying engine temperature measurements. As part of this investigation, we study the temperature measurements obtained from the same axial plane across five different engines yielding a total of 82 datasets. The five different engines have similar architectures and therefore similar temperature spatial harmonics are expected. Our problem is to estimate the spatial field in engine temperature given a few measurements obtained from thermocouples positioned on a set of rakes. Our motivation for doing so is to understand key engine temperature modes that cannot be captured in a rig or in computational simulations, as the cause of these modes may not be replicated in these simpler environments. To this end, we develop a multivariate linear least-squares model with Tikhonov regularization to estimate the 2D temperature spatial field. Our model uses a Fourier expansion in the circumferential direction and a quadratic polynomial expansion in the radial direction. One important component of our modeling framework is the selection of model parameters, i.e., the harmonics in the circumferential direction. A training-testing paradigm is proposed and applied to quantify the harmonics.