In the paper, we develop three new methods for estimating unknown initial temperature in a backward time-fractional diffusion problem, which is transformed to a space-dependent inverse source problem for a new variable in the first method. Then, the initial temperature can be recovered by solving a second-order boundary value problem. The boundary functions and a unique zero element constitute a group symmetry. We derive energetic boundary functions in the symmetry group as the bases to retrieve the source term as an unknown function of space and time. In the second method, the solution bases are energetic boundary functions, and then by collocating the governing equation we obtain the expansion coefficients for retrieving the entire solution and initial temperature. For the first two methods, boundary fluxes are over-specified to retrieve the initial condition. In the third method, we give two boundary conditions and a final time temperature to construct the bases in another symmetry group; the governing equation is collocated to a linear system to obtain the whole solution (initial temperature involved). These three methods are assessed and compared by numerical experiments.