Studies of hyperbolic heat conduction have so far been limited mostly to one-dimensional frameworks. For two-dimensional problems, the reflection and interaction of oblique thermal waves and complicated geometries present a challenge. This paper describes a numerical solution of two-dimensional hyperbolic heat conduction by high-resolution schemes. First, the governing equations are transformed from Cartesian coordinates into generalized curvilinear coordinates. Then the dependent variables are cast in a characteristic form that decouples the original system equation into scalar equations. Two-dimensional high-resolution numerical schemes, suck as total variational diminishing ( TVD) are built up by forming symmetrical products of one-dimensional difference operators on each individual wave. Three examples are used to demonstrate the unique feature of complicated interaction of two-dimensional thermal waves.