The decomposition of an image into a linear combination of digitised basis functions is an everyday task in astronomy. A general method is presented for performing such a decomposition optimally into an arbitrary set of digitised basis functions, which may be linearly dependent, non-orthogonal and incomplete. It is shown that such circumstances may result even from the digitisation of continuous basis functions that are orthogonal and complete. In particular, digitised shapelet basis functions are investigated and are shown to suffer from such difficulties. As a result the standard method of performing shapelet analysis produces unnecessarily inaccurate decompositions. The optimal method presented here is shown to yield more accurate decompositions in all cases.

A self-consistent approach to nonequilibrium radiation temperature is introduced using the distribution of the energy over states. We begin rigorously with ensembles of Hilbert spaces and end with practical examples based mainly on the far from equilibrium radiation of lasers. We show that very high, but not infinite, laser radiation temperatures depend on intensity and frequency. Heuristic "temperatures" derived from a misapplication of equilibrium arguments are shown to be incorrect. More general conditions for the validity of nonequilibrium temperatures are also established.