A numerical inversion approach to detect and localize inclusions in thick walls under quasi-periodic natural solicitations is presented. It is based on a preliminary analysis of surface temperature field evolution with time. This analysis is improved by taking advantage of a priori information provided by ground-penetrating radar reconstructions. In this way, it is possible to improve the accuracy of the images achievable with the stand-alone thermal reconstruction method.