Estimating Ising Models in Total Variation Distance