S U M M A R YThe Common-Reflection-Surface (CRS) stack method is a powerful tool to produce highquality stacked images of multicoverage seismic data. As a result of the CRS stack, not only a stacked section, but also a number of attributes defined at each point of that section, are produced. In this way, one can think of the CRS stack method as a transformation from data space to attribute space. Being a purely kinematic method, the CRS stack lacks amplitude information that can be useful for many purposes. Here we propose to fill this gap by means of a combined use of a zero-offset section (that could be a short-offset or amplitude-corrected stacked section) and common midpoint gather. We present an algorithm for an inverse CRS transformation, namely one that (approximately) transforms the CRS attributes back to data space. First synthetic tests provide satisfying results for the two simple cases of single dippingplane and single circular reflectors with a homogeneous overburden, and provide estimates of the range of applicability, in both midpoint and offset directions. We further present an application for interpolating missing traces in a near-surface, high-resolution seismic experiment, conducted in the alluvial plain of the river Gave de Pau, near Assat, southern France, showing its ability to build coherent signals, where recording was not available. A somewhat unexpected good feature of the algorithm, is that it seems capable to reconstruct signals even in muted parts of the section.The Common-Reflection-Surface (CRS) stack method is a recent data-driven time imaging process (see, e.g. Hubral 1999;Jäger et al. 2001 and also references therein) that has been originally proposed as an alternative to the classical normal moveout (NMO)-dip moveout (DMO) chain (Yilmaz 2000) to build seismic stacked, simulated zero-offset (ZO) time images of the subsurface. As already discussed elsewhere (Perroud & Tygel 2005), the CRS stack method has both advantages and disadvantages with respect to its alternative approaches. In fact, the adoption of the CRS stack method by the geophysical community has been until now only limited, because the classical NMO-DMO chain already provides good-quality robust results, so the need for a change is not obvious.However, the CRS stack method does not provide only ZO images, but also a set of wavefield attributes (also called CRS parameters: emergence angles and wave front curvatures) that have been exploited in several applications. These include, for example, velocity model building (Della-Moretta et al. 2001;Klüver 2006), multiple attenuation (Prüssmann et al. 2006) or residual statics correction (Koglin et al. 2006).The CRS stack method can be seen as a transformation from the data space (seismic amplitudes as a function of position and time) into attribute space (wavefield attributes as a function of position and time). Note that data space position variables include both midpoint and offset coordinates, while the attribute space position variables consist in the midpoint coordinat...