We consider an inverse problem of time fractional heat conduction problem. It is shown that the problem is ill-posed. A method is investigated based on the finite difference to find heat distribution and boundary values. The discrete mollification regularisation is applied to obtain a stable numerical solution. Finally, some test problems are investigated to show the ability of the proposed scheme.