Bottomonia play a crucial role in our understanding of the quark gluon plasma. We present lattice non-relativistic QCD calculations of bottomonia at temperatures in the range $T\in\left[47,\,380\right]$ MeV using the FASTSUM Generation 2L anisotropic $N_f = 2 + 1$ ensembles. The use of a basis of smeared operators allows the extraction of excited-state masses at zero temperature and an investigation of their thermal properties at non-zero temperature. We dind that the ground state signal is substantially improved by this variational approach at finite temperature. We also apply the time-derivative moments approach to the projected or optimal correlation functions at finite temperature.