You have to specify what you measure when measuring a temperature.
A nice example is the temperature of the cosmic microwave background radiation. We always quote a temperature ##T \simeq 2.725 \; \text{K}##. So what is this temperature and how is it measured?
What's really measured by, e.g., the Planck satellite is the spectrum of the radiation. To understand what's measured it's sufficient to know that the cosmic background radiation (locally) has the spectrum of a black body and the satellite is moving, together with our entire solar system, relative to the restframe of this thermal radiation with a certain velocity ##\vec{v}##. So which spectrum does the satellite measure?
In the local rest frame of the CMBR the spectrum, i.e., the number of photons per unit volume and per unit volume in momentum space (i.e., the phase-space distribution function) is given by
$$f^{*}(\vec{k})=\frac{2}{\exp(\beta |\vec{k}^*|)-1},$$
where the ##2## comes from the two polarization degrees of freedom of photons. Here and in the following I use natural units with ##\hbar=c=k_{\text{B}}=1## and ##\beta=1/T## with ##T## the temperature as defined in modern terms as a scalar quantity. A comoving CMBR satellite, i.e., a satellite at rest relative to the restframe of the CMBR would just measure the isotropic temperature ##T##, which by definition is what we call temperature in the modern relativistic formulation of thermodynamics and statistical physics in terms of manifestly covariant quantities. In this formulation the phase-space distribution function of photons is a scalar field.
It's thus easy to find the phase-space distribution function for the real Planck satellite. We just have to write ##f^*## (the distribution function as seen from the point of view of an observer at rest relative to the restframe of the black-body radiation) in a manifestly covariant way. The restframe observer's four-velocity is simply ##u^*=(1,0,0,0)##. Now the momentum dependence in ##f^*## is to be understood to be valid for the on-shell four-momentum of a photon obeying ##k_{\mu} k^{\mu}=0##, i.e., ##|\vec{k}|=k^0## in any frame of reference. For the rest frame we thus have ##|\vec{k}^*|=u^{*\mu} k_{\mu}^*##, and that's a manifestly covariant expression, i.e., we can immediately write the phase-space distribution function in any frame as
$$f(k)=\frac{2}{\exp(\beta u^{\mu} k_{\mu})-1}.$$
So what's measuring the Planck satellite? Since ##u=\gamma (1,\vec{v})## with ##\gamma=1/\sqrt{1-\vec{v}^2}## we have
$$u^{\mu} k_{\mu}=\frac{|\vec{k}|}{\sqrt{1-\vec{v}^2}}(1-|\vec{v}|\cos \vartheta),$$
where ##\cos \vartheta=\vec{k} \cdot \vec{v}/(|\vec{k}||vec{v}|)##.
So for any fixed direction you observe a black-body spectrum with an effective temperature given by
$$\tilde{\beta}=\frac{1}{\tilde{T}}=\frac{\beta}{\sqrt{1-\vec{v}^2}}(1-|\vec{v}| \cos \vartheta)$$
or
$$\tilde{T}=T \frac{\sqrt{1-\vec{v}^2}}{1-|\vec{v}| \cos \vartheta}.$$
That's also easily interpreted. For ##\vartheta=0## what the satellite observes are photons coming from a source of black-body radiation coming towards the detector, i.e., you have a blue shift and thus you measure a higher temperature
$$\tilde{T}(\vartheta=0)=T \sqrt{\frac{1+|\vec{v}|}{1-|\vec{v}|}}.$$
For ##\vartheta=\pi## it's measuring photons from a source moving away from the detector with the corresponding red shift
$$\tilde{T}(\vartheta=\pi)=T \sqrt{\frac{1-|\vec{v}|}{1+|\vec{v}|}}.$$
For values of ##\vartheta## somewhere in between you get the corresponding blue and red shifted photons leading to a temperature given by the general formula. For ##\vartheta=\pi/2## the redshift is due to the transverse Doppler effect, where the observed frequency change from the frequency as measured in the restframe of the source is solely due to relativistic time dilation.
One should note that indeed the CMBR satellites measure this "dipole contribution" to the variation of the temperature with direction and what's shown as pictures of the CMBR temperature variations is corrected by this dipole contribution, such that you get a nearly isotropic temperature of ##T=2.725 \; \text{K}## quoted above. The slight variations in temperature in dependence of direction of ##\Delta T/T \simeq 10^{-5}##, i.e., the deviations to a isotropic black-body spectrum are then the really interesting signal for cosmologists. From the dipole contribution to the derivation we can deduce that our solar system is moving with about ##370 \; \text{km}/\text{s}## towards the constellation Leo relative to the rest frame of the CMBR. For details see
https://en.wikipedia.org/wiki/Cosmic_microwave_background#Features
and the quoted reference from the PLANCK collaboration:
https://arxiv.org/abs/1807.06205
Another example, where the meausurement of temperatures of a moving medium becomes relevant is in my own field of research, ultrarelativistic heavy-ion collisions: If you collide heavy nuclei together (as done at the Relativistic Heavy Ion Collider at the Brookhaven National Lab and at the Large Hardon Colider at CERN) a very hot and dense medium of strongly interacting matter is formed which can be described with astonishing precision as a blob of collectively streaming almost ideal fluid, which implies that it is, most of the time, in a state close to local thermal equilibrium and thus a (local) temperature is, in principle, well defined. So how to measure "the temperature"? That's a not so simple story, and indeed you have to account for the flow of the medium, and you need different measurements to get a (rough) picture of the temperature evolution of this rapidly expanding and cooling fireball.
One meausurement is simply the abundance of all kinds of hadrons produced in the collision, and since the particle numbers are scalar quantities you can indeed infer a temperature from it by just fitting the measured abundancies of different hadrons to a model, where all hadrons are in thermal equilibrium at the time the "chemistry" of the fireball is fixed. That's called the chemical freezeout. The so deduced temperature at the highest collision energies turns out to be around ##155 \; \text{MeV}##, close to the pseudocritical transition temperature for a zero net-baryon-number (baryo-chemical potential ##mu_{\text{B}}=0##) medium of strongly interacting particles between a state where quarks and gluons are the relevant degrees of freedom ("quark gluon plasma") and a "hadron resonance gas", which is known from lattice-QCD calculations at finite temperature. That means that shortly after the fireball undergoes the (cross-over) transition from a QGP to a hadron-resonance-gas state the inelastic collisions cease and almost only elastic collisions stay relevant.
The other signal you can use to measure a temperature from hadronic observables are the transverse-momentum spectra of the hadrons. These hadrons have a momentum distribution as given after "thermal freeze-out", i.e., it's the spectrum at the time when also the elastic collisions don't play a role anymore in the fireball distribution and the particles are freely streaming towards the detector (the decay of unstable hadrons to stable ones has to be taken into account as well). This spectrum can be fitted well with a socalled "blast wave model", assuming that one has a collectively streaming medium at local thermal equilibrium. To get the true temperature you have to take into account the flow of the particles to get a true "scalar" temperature of about ##100 \; \text{MeV}##.
Another way to measure an space-time weighted average temperature over the entire fireball evolution is to measure the spectrum of emitted photons. Since they don't strongly interact they come out of the fireball nearly undisturbed, i.e., they are emitted from a flowing source in thermal equilibrium, and again you have to take into account the blue shift from the radial motion of this source to get a true "scalar" average temperature.
Last but not least another related signal are "dileptons", i.e., the measurement of electron-positron or muon-antimuon pairs which come mostly from quark-antiquark annihilation in the QGP phase and from decays of the light vector mesons in the hot and dense hadron-resonance-gas phase of the medium. There you can, in addition to the transverse-momentum spectrum, also measure the invariant-mass spectrum of the dileptons. Since invariant mass is a scalar in principle you get directly the "scalar temperature" by measuring the slope of the invariant-mass spectrum, but there you also have the shape of the vector-meson resonance spectrum, and this shape is majorly changed in the hot and dense medium as compared to the vacuum. There's however a "mass window" between the ##\phi##-meson mass and the ##J/\psi## mass where you have a pretty flat spectral function of the here relevant autocorrelation function of the electromagnetic current, such that you get a pretty much undistorted thermal spectrum from the part of the dileptons which is due to quark-antiquark annihilations in the QGP phase of the fireball expansion. In addition the mass range of about ##1 \; \text{GeV}## to ##2 \; \text{GeV}## which is quite a bit larger than the usual temperatures of the fireball (of around a few 100 MeV) the emission from the early hot phases of the fireball evolution is the dominant contribution. The only obstacle here is that in this mass region also another dilepton contribution is relevant, namely, the semileptonic decay of correlated charm-anticharm pairs (or D and anti-D mesons in the hadronic phase), which have to be somehow subtracted to get a clean spectrum from the thermal dilepton radiation from the QGP phase.