I think you are on the right track in thinking about something *like* Zernike polynomials. As noted in the answer by jwimberly, these are an example of a  system of [orthogonal basis functions][1] on a disk. I am not familiar with Zernike polynomials, but *many* other families of orthogonal functions (including Bessel functions) arise naturally in classical mathematical physics as [eigenfunctions][2] for certain partial differential equations (at the time of this writing, the animation at the top of that link even shows an example of a vibrating drum head). 

Two questions come to my mind. First, if all you are after is the radial profile ($\theta$ averaged), then how much constraint on the spatial pattern do you need? Second, what types of variability occur in the spatio-temporal data?

In terms of the first question, there are two concerns that come to mind. Due to the polar coordinates, the support-area for each sensor has a trend with $r$. The second concern would be the possibility of [aliasing][3], essentially a mis-alignment of your sensors relative to the phase of the pattern (to use a Fourier analogy). Note that aliasing will likely be the primary uncertainty in constraining the peak temperatures (i.e. $T_{95}$).

In terms of this second question, data variability could actually help with any aliasing issues, essentially allowing any mis-alignment to average out over the different measurements. (Assuming no systematic bias ... but that would be a problem for any method, without e.g. a physical model to give more information). 

So one possibility would be to define your spatial orthogonal functions purely at the sensor locations. These ["Empirical Orthogonal Functions"][4] could be computed via [PCA][5] on your spatiotemporal data matrix. (Possibly you could use some weighting to account for the variable sensor support areas, but given the uniform polar grid and target of radial averages, this may not be required.) 

Note that if there **is** any physical modeling data available for "expected" variations in the temperature, available on a dense spatiotemporal computational grid, then the same PCA procedure could be applied to *that* data to derive orthogonal functions. (This would typically called "[Proper Orthogonal Decomposition][6]" in engineering, where it is used for model reduction, e.g. an expensive computational fluid dynamics model can be distilled for use in further design activities.)

A final comment, if you were to weight the sensor data by support area (i.e. polar cell size), this would be a type of diagonal covariance, in framework of [GLS][7]. (That would apply to your prediction problem more, although weighted PCA would be closely related.)

I hope this helps!

  [1]: https://en.wikipedia.org/wiki/Orthogonal_functions
  [2]: https://en.wikipedia.org/wiki/Eigenfunction
  [3]: https://en.wikipedia.org/wiki/Aliasing
  [4]: https://en.wikipedia.org/wiki/Empirical_orthogonal_functions
  [5]: http://stats.stackexchange.com/questions/tagged/pca?sort=votes
  [6]: https://scholar.google.com/scholar?q=proper%20orthogonal%20decomposition
  [7]: https://en.wikipedia.org/wiki/Generalized_least_squares