Multi layer finite difference method and other methods were used to calculate the surface temperature changes of different parts of city buildings and their surrounding parts with time in different seasons in a year considering the effect of outer meteorological factors. Using methods of Computer Graphics, modelings and realistic IR images of buildings in different seasons and at different time in a day were realized and drawn.