In this paper, we present a procedure to map subsidence at the regional scale by means of persistent scatterer interferometry (PSI). Subsidence analysis is usually restricted to plain areas and where the presence of this phenomenon is already known. The proposed procedure allows a fast identification of subsidences in large and hilly-mountainous areas. The test area is the Tuscany region, in Central Italy, where several areas are affected by natural and anthropogenic subsidence and where PSI data acquired by the Envisat satellite are available both in ascending and descending orbit. The procedure consists of the definition of the vertical and horizontal components of the deformation measured by satellite at first, then of the calculation of the “real” displacement direction, so that mainly vertical deformations can be individuated and mapped.