This study applies a hybrid technique of the Laplace transform and finite-difference methods in conjunction with the least-squares method and experimental temperature data inside the test material to predict the unknown surface temperature, heat flux and absorptivity for various surface conditions in the laser surface heating process. In this study, the functional form of the surface temperature is unknown a priori and is assumed to be a function of time before performing the inverse calculation. In addition, the whole time domain is divided into several analysis sub-time intervals and then these unknown estimates on each analysis interval can be predicted. In order to show the accuracy of the present inverse method, comparisons are made among the present estimates, direct results and previous results, showing that the present estimates agree with the direct results for the simulated problem. However, the present estimates of the surface absorptivity deviate slightly from previous estimated results under the assumption of constant thermal properties. The effect of the surface conditions on the surface absorptivity and temperature is not negligible.