In this paper, we present a quite simple recursive method for the construction of classical tensor product Hermite spline interpolant of a function defined on a rectangular domain. We show that this function can be written under a recursive form and a sum of particular splines that have interesting properties. As application of this method, we give an algorithm which allows to compress Hermite data. In order to illustrate our results, some numerical examples are presented.