We present the thermal-magnetic relation in a simple, isolated sunspot deduced from the inversion of 1.56 μm spectropolarimetric data. Due to the high Zeeman sensitivity of the g = 3, Fe I 1.5648 μm line, we can study this relationship in the entire sunspot. An inversion technique based on response functions is used to derive various parameters, both as a function of location within the sunspot and of height in the atmosphere. In this paper we attempt to relate field strength, vertical and radial field components and the field inclination with temperature.