The reconstruction of volcanic landforms is an effective way to assess the construction history and growth rates of volcanic structures and is a useful tool to estimate erosion rates. We propose a simple geometrical approach to reveal the location of the long-term eruption centers and the general shape of eroded volcanoes within the altitudinal range of the input data. This DEM-based approach consists of finding the radial topographic profile of an eroded volcanic landform from its remaining surfaces, by coupling two independent and simple mathematical methods (linear regression and Spearman rank correlation). Here we demonstrate that this simple approach is suitable for circular to elliptical (up to an ellipticity of 0.5) and to some extent for asymmetrical volcanoes, and we apply it to Réunion Island. Our results confirm previous conclusions from geological observations and data requiring more complex instrumentation. We confirm that the center of long-term volcanic activity of Piton des Neiges remained stable for the past 1.4 Myr while Piton de la Fournaise has been migrating eastward for~300 kyr. Our approach also reveals that Piton des Neiges is the dominant volcano in Réunion and that Piton de la Fournaise is a smaller edifice built on its edge. Moreover, despite the contribution of two other volcanoes in the construction of Réunion, we suggest that the Piton des Neiges has rapidly become the main structure of the whole island and has subsumed the previous volcanic centers. Finally, we use our approach to model the magmatic activity of Piton des Neiges, and suggest a deceleration of its magma emission rate since 2Myrs.