Introduction

From the viewpoint of classical field praxis, before the impact of computer applications, the documentation of a site with prehistoric rock art or with Prehistoric Graphic Expression (hereafter PGE)Footnote 1 had as its primary objective to obtain ‘reliable tracings’ and some good photos of the icons and the surrounding rock support. Within this analogical context of the discipline, successful completion of the task involved producing a cumbersome cellophane filled with coloured lines, which was then scaled with a pantograph or photocopied in parts to create a visually pleasing final image.

The inertia of that procedure justifies the reflection of Peter Ucko in the 1980s, in which he warned us that, with these routines, we unconsciously assume that we all ‘know what we are doing and that we have good theoretical support to do it so’ (Ucko 1989: 285). Indeed, this apparently consensual documentation pattern of the PGE, being more or less technical and with an equally common vocabulary, seemed congruent until faced with results generated by other researchers on the same documentation space. If the lines on the cellophane did not coincide (a common occurrence), in what terms was the discussion established? Was one documentation more objective than another? Who was right in drawing the lines?

‘Original image’ and ‘tracing’. These were the two defining elements of the method. They were essential and, in a certain way, axiomatic, which in many cases obviated the need for definition.

However, from the discipline’s early moments, there has been a framework for reflection on these methodological principles. For example, on the role of the image (Lorblanchet 1993: 330), the concept of fidelity to an ‘original’ (Lemozi 1929: 42), or on objectivity in the tracing process (Aujoulat 1987: 9; Raphael 1986: 133). The issues addressed by these and other authors are still relevant to current research. However, the development of new digital strategies for field work and cabinet processing has helped us to understand, from an unprecedented perspective, the theoretical problems dealt with in that analogical context of the discipline. We can say, then, that reflection born of the digital discussion behaves as a metalanguage which orders the theoretical issues dealt with by classical historiography in another dimension. In other words, the current ability to quantify documentation variables that were previously unapproachable has let us transfer to the field of ‘data’, documentation gestures that belonged, in the classical method, to interpretations. Although, in this context, we have not been aware of this change of role due to a lack of perspective.

Let us take an example: If we analyse a publication on PGE, especially if it predates digital archaeology, although it is easy to identify the internal order of its text (we can isolate an introduction, a descriptive chapter of the observable evidence and a section of interpretative conclusions), it is not unusual for doubts to arise when referring to the role of the graphic set that goes along with this text, i.e. the illustrations, the tracings. Do these belong to the exposition of data, or are they the product of the interpretative conclusions? Compared with what we could show for other archaeological processes, where was the data found in our study? In the cellophane tracings?

What does the current technological scenario contribute to the dialectic of classical documentation (original image / tracing)? During the first 20 years of the digital evolution of the discipline, as a general rule, the specialists, trained in a pre-digital era, saw in the ‘new technologies’ the possibility of being more efficient in the method, ‘the Method’ that is as ever defined by this ‘original image / tracing’ duality. Thus, vector graphic apps effectively replaced cumbersome cellophane and markers.

Three-dimensional models made it possible to create grazing shadows synthetically, using light and shadow maps. With this, our brain understood what was engraved, and the outline of its edges was safely delineated using a vector graphic application with its elegant B-splines. However, the method remained locked in its duality: ‘original image / tracing’, and the tracings still did not correspond between different researchers.

The digital world has gone up exponentially and overwhelmingly our capacity to calculate, store and analyse data. As we say, this has made it possible to quantify previously unattainable variables and multiply the points of view or the documentation strategies of the prehistoric graphic original. Suppose we order these circumstances in a methodological reflection; in that case, the elements that articulated the classical duality of the original image / tracing scenario have been changed in their defining parts and have confirmed the role of a novel third part, i.e. the data.

Let us dwell for a moment on each defining element of the documentation process. On the one hand, we must abandon the idea that at the outset of the process, there is an ‘original’ understood as an image. Any image is a construct of our brain, whether we refer to direct visualisation of the object of study or if we use intermediaries such as specialised radiometric or geometric sensors. From a different perspective, computing has taught us that when we document PGE, we face a complex physical reality from which we extract some particular variable that let us reconstruct the graphic message, i.e. the image. This elective act is subjective, but it can be quantified.

A second intrinsically digital circumstance is that we can now present an abundant and ordered array of data in our method. This factor begins with subjective choices; however, these are quantifiable and reproducible by other researchers once their features have been defined. At the same time, the same observer/researcher can reproduce the variables that define the method in different documentation scenarios. This new method, which we have noted is intrinsically digital, makes the documentation modes in the PGE concurrent with any other general archaeology or science research process.

Finally, we must conclude that the interpretive graphic is inescapable. The prehistorian must interpret his data, and that interpretation, in the PGE, will be a graphic document. The difference is that now, we can support this interpretive document in a structured data section (Fig. 1).

Fig. 1
figure 1

The ‘method’ and the impact of digital environments

However, the structure of this new section (the data) is far from being normalised. The main problem is that in our new habits of digital research, we are reproducing the tics that Peter Ucko talked about in the 1980s of the twentieth century. We use digital terms and concepts we do not define because we carry the false idea that we all know what we are talking about. Electronic tracing, data capture, 3D, 4D, virtual reality, orthoimage and giga-image are frequent terms in current documentation processes.

Do we all have the same understanding when we use them? How does this lack of definition interfere with the distinction between data and interpretation?

Undoubtedly, one of the most common terms and one of the most conflicting is 3D. Let us pause momentarily and reflect on the intrinsically related idea of ‘virtual reality’. Beyond its Euclidean definition (length, width, depth), most meanings associated with this term concern ‘visual experience’ (Soska et al. 2010; Storrs and Fleming 2021), if we understand as such our ability to discern between sizes and distances. This ability is a complex process based on the maturation of our visual system, which in this case is reproduced in two-dimensional settings through images of successive reproduction in time (4D). Therefore, we are talking about the individual apprehension of a subjective experience. Can we use these virtual reality scenarios to augment and structure our quantifiable data section, making it interchangeable between other observers? On many occasions, researchers/observers claim to understand three-dimensional space within these virtual scenarios as they have never experienced it before. However, when an observation differs from the others, even using this visual tool, how does that affect the discussion? Where is the quantifiable data in virtual reality?

Let us take this theme further, imagining that virtual reality worked ideally, as some hope, and that we could not distinguish the features of physical reality from its virtual version.

Would the problems be over for the researcher? Would the documentation objectives have been achieved? We believe not. Virtual reality would have taken us to the starting point, i.e. the visualisation of a panel. In this virtual or real scenario, we will probably hear interjections such as ‘fantastic’ and ‘impressive’… Everything will be fine until one observer asks: ‘Okay, but what do you see?’ In response, someone will take the initiative and begin to outline a line in the air with their finger. That is the core of the problem in the PGE: to transmit to another observer what one sees or interprets.

Insofar as it immerses us in an individual viewing experience, virtual reality moves us away from the notion of objective and quantifiable knowledge of the data. It invites us to interpret (Guertin-Lahoud et al. 2023: 264). Rather, what we need to create in our documentation process are metric documents (Domingo Sanz et al. 2013), where can these be found in virtual reality?

Let us look at the cases in a general sense. To build a data section in our documentation process, we must move away from graphic solutions that appeal to our ‘visual experience’, which is true for 3D scenarios, such as virtual reality. It also applies to two-dimensional documents, such as synthetic light and shadow maps that, in the same way, appeal to our visual experience to understand what is concave, what is flat, or what is convex, and to be able, based on that, to draw our interpretive B-splines in a vector graphic environment. Here, we would return to the rock-cellophane loop, i.e. physical reality / tracing.

However, this section which is key in bringing our documentation methods closer to other scientific processes that exchange data should be structured from graphic solutions far removed from the ‘visual experience’. These will generate actual interchangeable metric documents. We refer, for example, to the principal component analysis (PCA) (Cerrillo-Cuenca et al. 2014) of raster files generated by specialised radiometric sensors (Bayarri et al. 2021b). We are also referring to the analysis of geometric peculiarities of a 3D model (facet-edges, discrete curvatures, Hausdorff distances, morphological residual models etc.) (Cerrillo-Cuenca et al. 2019; Pires et al. 2014; Rivera Jiménez et al. 2021; Rogerio-Candelera and Linares 2015; Teira Mayolini and Ontañón Peredo 2016). They are not definitive, incontestable solutions. As we said, the variables have been chosen subjectively but are quantifiable.

Consistent with the above, in our methodological proposal, we aim to present an individualised and coherent data structure. From this data, we will create a ‘third order’ of interpretive images (Montero Ruiz et al. 1998).

We have chosen a classic case study, an obligatory reference in Late Peninsular Prehistory synthesis literature which has been visited by multiple researchers throughout the twentieth century and which brings together, on the same panels, engraved and painted icons. The most significant ones contain both techniques of expression. This is the main challenge of the study. The combined use of radiometric and geometric techniques to obtain metric documents lets us analyse this graphic narrative from quantifiable variables.

Regarding radiometric exploration, we will demonstrate the suitability of using a hyperspectral sensor with a working band between 350 and 1050 nm. With it, we have isolated several phases of pigment additions from subtly differentiated spectral signatures.

The geometric exploration of the rock surface is aimed at distinguishing between natural roughness and engraving grooves traced by the human hand. This implies an adequate data capture, i.e. accurate and high resolution, and processing of these data that differentiates between a natural or artificial origin of the micro-reliefs. To achieve this goal, we will rely on graphic solutions outside the visual experience (hidden geometries, search for edges, discrete curvatures etc.) and others within the visual experience (light and shadow maps), which are more intuitive but difficult to quantify.

Case study: Peña Tu’s rock (Puertas de Vidiago, Asturias, Spain)

Peña Tu is in the north of Spain, in the central longitudes of the Cantabrian Region, at the eastern end of the region of Asturias. It is a prominent, isolated rock at the end of one of the Sierras Planas de Llanes, erosive witnesses of ancient levels of marine terraces, carved on quartzite and sandstone from the Cambrian and Ordovician periods, also affected by phylogenetic tectonic movements, and which define two altitudinal levels, around 150 m and 220 m amsl (Domínguez-Cuesta et al. 2019; López-Fernández et al. 2020).

The Sierras Planas form a plateau between the Bay of Biscay, 1.6 km to the north, and the mountain barrier of the Sierra de Cuera, 4 km to the south, which, in this longitude, reaches altitudes above 1100 m asl. All these mountainous features, running parallel to the coast, are the prelude to the huge massif of the Picos de Europa, with several peaks above 2600 m. This rugged terrain forms a narrow coastal communication corridor intensely inhabited by humans during prehistoric times. It is one area with the highest density of archaeological sites in the Cantabrian region (Fano Martínez 1998: 13) (Fig. 2).

Fig. 2
figure 2

Location of Peña Tu and other nearby prehistoric sites

The rock or ‘Peña’ is a residual outcrop of the quartz arenites and white sandstones of the ‘Barrios’ formation, which form these flat mountain ranges (Gutiérrez-Alonso et al. 2007). It has a height on its southern side (the most prominent, in favour of the slope of the hillside) of about 10 m. The opposite side, to the north, is only 5.8 m above the ground. Between the two opposing sides, there is an average thickness of about 4.8 m. It is on the eastern side of this thickness where the decorated panel is located (Figs. 3 and 4). The rock has an east-west length of 7.8 m.

Fig. 3
figure 3

View of the Peña Tu site from the east in the evening of one of the documentation days. The pre-coastal mountain ranges around Llanes can be seen in the background

Fig. 4
figure 4

Plan and elevations of the decorated rock of Peña Tu

Peña Tu was identified in the summer of 1913 By Eduardo Hernández Pacheco, the Count of Vega del Sella y Juan Cabré, and published by its discoverers the following year, once a first tracing of the set was composed in the autumn of that year 1913 (Hernández-Pacheco et al. 1914) (Fig. 5a). This was done freehand. It describes all the major elements that frame the formal and chronological discussions on the decorated panel throughout the twentieth century and the beginning of the 21st. That is, the authors describe (A) a larger anthropomorphic figure, engraved and painted with ochre strokes, slightly more than one metre high and about 60 cm wide, arranged at the extreme right of the panel (to the north); (B) the representation of a weapon arranged vertically next to it on the left, also engraved and with remains of paint; and (C) various depictions painted in red, always to the left of the larger anthropomorphic representation, among which we can distinguish groups of anthropomorphs of a classical schematic convention, series of points distributed in a nuclear and linear way, and other signs difficult to interpret due to their state of conservation. The entire decorated panel spans an area of just over 4 m wide and 1.3 m high.

Fig. 5
figure 5

Comparative image of the Peña Tu tracings published in 1914 (a) and 2007 (b). In the case of the 1014 tracing, as it only has a numerical scale (1:8), the print size of the original figure has been checked to be able to reproduce the two tracings at the same scale

Only certain icons carved on the wall using a picketing technique are omitted in the tracing, probably because, considering the authors, these motifs were recent; they should not be part of the archaeological analysis.

In that publication, the tracing is printed on a fold-out sheet in two inks (Red and black), which refer to the two techniques used in the graphic narrative (painting and engraving). There is no graphic reference to ridges of the rock support. A graphic scale was also not drawn, although a numerical fraction scale of 1:8 is indicated. There is no reference to the horizontal level of the composition. As we say, etched grooves are rendered in black with different thicknesses, using a shading convention. The painted areas are shown in red ink with different intensities, thus suggesting different conservation and sharpness of the icons. This early 20th-century tracing was the canonical image of the site for over six decades until the review was published in the 1980s.

With a different scientific paradigm and renewed reproduction techniques, a team led by P. Bueno and M. Fernández Miranda carried out a new tracing of the panel in 1981 (Bueno and Fernández-Miranda 1981). In this case, it was done using cellophane in direct contact with the rocky surface. Unfortunately, the graphic work of these authors was largely ruined by the faulty edition of the book in which it was published. The lack of tinting of the sheets on which the tracing was shown made most icons illegible (Fig. 6). For this reason, the narrative of the decoration of the panels can only be followed from the descriptions in the text, which are themselves extensive and detailed.

Fig. 6
figure 6

Partial image of the tracing of El Peñatu after (Bueno and Fernández-Miranda 1981)

One author of this review, Primitiva Bueno, researched the site again twenty-six years later. It is the last major revision of the panel and was published as a diachronic context for the discovery of human remains in the nearby Fuentenegroso cave (Barroso Bermejo et al. 2007) (Fig. 5b). In this case, the authors tell us that for the elaboration of their tracing, they have not resorted to methods that imply contact with the surface, but rather that they have used more sophisticated photographic techniques. This change in technical strategy brought with it some notable variations to the interpretation of the icons, as we will see in the ‘Discussion: critical and comparative analysis of data and interpretations’ section.

On this occasion, the tracing incorporates graphic ‘shading’ conventions that try to represent the natural roughness of the rock support and the pecked areas that represent recent cruciforms. Neither of these related aspects of the feature was discussed in the 1914 tracing. Nevertheless, the same graphic convention is used for the grooves delineating the greater anthropomorphic figure and the weapon. The icons and areas with pigment are represented in red ink with slight variations in intensity, which appear to represent different degrees of preservation of the original lines.

From a different viewpoint, three years later, M.A. de Blas published a new interpretation of the site (Blas Cortina 2010). This author invites us to consider, with reference to the symbolic meaning of Peña Tu, not only the surfaces decorated with paintings and engravings but also the rock support itself: It is a stone volume with a supposed zoomorphic appearance. The work does not discuss in detail the appearance of the panel icons, nor does it reproduce a complete view of them. However, it does expand on the description of the weapon arranged next to the main anthropomorphic figure using graphics and text, which is interesting considering our methodological purpose. In addition to this icon, he also reproduces a tracing of one of the small-sized anthropomorphic figures painted in red, the interpretation of which has interested researchers the most throughout the last century (Fig. 7). In none of the cases is any detail of the methodology followed for the elaboration of the tracings specified. However, his peculiar interpretation of both figures will help us show; in the ‘Discussion: critical and comparative analysis of data and interpretations’ section, the abundance and variety of explanations Peña Tu has offered throughout the history of the investigation.

Fig. 7
figure 7

Representation of the anthropomorphic figure with crozier (A) and the weapon (B) of Peña Tu after Blas Cortina (2010)

Technical scenario and methodology: approach of a dual strategy for the documentation of the archaeological site

The study of the painting and engraving techniques seen at the site involves two documentation strategies. The first, concerning the identification of pigment groups, is based on specialised electromagnetic spectrum sensors. The second, which defines the geometric reality of the engravings and the natural rock support, uses photogrammetry and laser scanning techniques.

Both strategies are carried out concurrently in the field, but their joint and complementary analysis is only possible through accurate georeferencing.

The prehistoric graphic expression needs comprehensive georeferencing to guarantee spatial coherence of the data in a time series, which can then be acquired to review its condition.

Every measurement has an error: No measurement is ever exact and must be rigorously adjusted, with the results undergoing statistical analysis. In geomatics, integrating global navigation satellite systems (GNSS), topographic total station (TTS), 3D terrestrial laser scanner (3DTLS), digital metric cameras and hyperspectral imaging systems all allow an accurate, reliable and rapid recording of information (Bayarri et al. 2021a, 2023). This information can be integrated into a geographic information system (GIS), which can be used extensively for research, management or decision-making (Bayarri et al. 2019; Ontañón et al. 2014).

Integrating 3DTLS and GNSS allows a rapid, accurate and reliable recording of complex features such are the caves or cavities (Bayarri and Castillo 2009). It plays a pivotal role in the successful hybridization of results obtained through the integration of methodologies such as 3DTLS, GNSS and photogrammetry (Bayarri 2020). The accurate alignment of spatial data is essential for combining information derived from these diverse sources seamlessly. In data collection, the 3DTLS captures high-resolution point cloud data, providing detailed and accurate (up to 1–2 mm in phased-based systems for high reflectances) geometric information of the surveyed area. Simultaneously, GNSS ensures accurate positioning by leveraging satellite signals, offering global coverage for geolocation with an accuracy of around a few centimetres. Photogrammetry complements these techniques by using images to reconstruct three-dimensional scenes. The 3DTLS is adjusted within a free network linked to GNSS coordinates, allowing for the propagation of errors (Bayarri et al. 2021a). This interconnected system makes sure the laser scanner’s measurements are aligned with accurate geographical references, enhancing the reliability and spatial coherence of the hybridised results. The ground control points (GCP) of photogrammetry are obtained from the adjusted 3DTLS point cloud, reaching a few millimetres of accuracy in their determination.

Hyperspectral imaging is typically acquired through scanning by using a rotary scanner. To effectively use these datasets alongside with spatial data, there is often a need to make a registration or georeferencing process to set a specific map coordinate system, in our case European Terrestrial Reference System 1989 (ETRS89). An essential step involves determining the precise position from which hyperspectral images were captured. This is achieved through the registration of GCP, which serve as reference markers in both the 3D model and the hyperspectral imagery. By aligning these control points, it is possible to establish the spatial relationship between the two datasets, helping with an accurate fusion of hyperspectral information into the three-dimensional model. This process not only enhances the spatial context of the hyperspectral data but also enables the reprojection of 2D hyperspectral information onto the 3D model’s surface. Integrating hyperspectral and 3D spatial data offers a holistic understanding of the rock art site, allowing for in-depth analysis of material composition, pigment distribution and hidden details.

The fusion of hyperspectral imaging with photogrammetry yields several useful outcomes. First, the union provides an enriched spatial context by integrating hyperspectral data with a 3D model. This synergy empowers researchers to better interpret the complex relationship between materials and the surface of rock art. Second, this combination enhances material discrimination capabilities, enabling the precise identification of pigments, binders and other materials used to make rock art. Third, the integration offers a holistic approach to analysis, leveraging the strengths of both techniques. Finally, the comprehensive dataset obtained could aid in developing effective preservation and conservation plans. It helps identify areas of deterioration, understands material vulnerabilities and guides targeted conservation efforts, thus helping with the long-term preservation of these culturally significant artworks.

Image registration serves as the initial phase, aligning hyperspectral and photogrammetric datasets for spatial coherence by correlating corresponding points between the two sets of images. Following registration, the integration stage combines the datasets, merging hyperspectral and photogrammetric data into a cohesive whole. This integration involves overlaying hyperspectral information onto the 3D model, associating spectral data with precise spatial coordinates. Then, a comprehensive analysis of the integrated dataset ensues, exploring the hyperspectral information within the 3D model. This in-depth examination allows for the identification of materials, pigments and subtle features present on the rock art surface, establishing correlations between spectral signatures and specific areas. Finally, the interpreted findings are contextualised within an archaeological framework. The integrated dataset enables a more nuanced understanding of rock art, aiding in pigment identification, unveiling concealed details and unveiling patterns that might remain obscured when using individual techniques in isolation.

Hyperspectral documentation of the rock art panel

Hyperspectral imaging is a powerful technology that has found applications in various fields, including rock art studies. It refers to the acquisition and analysis of images captured at multiple wavelengths in a region of the electromagnetic spectrum. The significance of hyperspectral imaging in rock art studies lies in its ability to reveal details and information that may be invisible or difficult to discern with the naked eye or traditional imaging methods. Unlike conventional photography or even multispectral imaging, which captures images at a few wavelengths, hyperspectral imaging records data across a continuous and wide range of wavelengths, often spanning the visible and near-infrared regions of the spectrum.

This spectral range lets researchers analyse the unique spectral signatures of different materials present in rock art, such as pigments, binders and underlying rock surfaces. By processing the reflectance or absorption patterns at various wavelengths, hyperspectral imaging can provide valuable insights into the composition, layering and even faded or obscured elements and contribute to analyse, interpret and preserve such valuable artwork.

Hyperspectral sensor

The hyperspectral system used for the study was formed by a Specim V10E spectrograph (Spectral Imaging Ltd., Oulu, Finland) covering the spectral range of 400–1000 nm, an sCMOS camera (CL-30, effective resolution 1312 × 728 pixels by 12 bits), an objective lens (Cinegon 2.4/30 mm focal length) and a rotating scanner light and data acquisition (DAQ) spectral software (Fig. 8).

Fig. 8
figure 8

View of the Specim V10E spectrograph scanning the lower bench of the site

The bandwidth was reduced to 5.6 nm by applying binning to improve the signal-to-noise ratio to produce 214 spectral bands from each pixel to form a relatively continuous spectral curve.

The data were recorded as a hypercube with two-dimensional spatial images and wavelength bands as the third dimension. Two hundred fourteen spectral bands were extracted from each pixel to form a relatively continuous spectral curve.

Illumination sources

Data was recorded in two ways. The first used natural direct sunshine incident light, and the second used a Philips Photolita Lamp Bulb 240 V 500 W. The spectral signature of the lights was obtained from the 1 nm wide narrowband FieldSpec Pro FR spectroradiometer made by analytical spectral devices, measuring spectra over a range of 350–1050 nm (Fig. 9).

Fig. 9
figure 9

Spectral signatures of lights at working distance. The vertical axis represents the relative spectral power measured, and the horizontal axis is the wavelength in nanometres where it was measured

Pre-processing: reflectance calibration (Fig. 10)

Fig. 10
figure 10

General WorkFlow, adapted from Bayarri (2020) and Bayarri et al. (2019, 2021a, b)

During the field campaign, apart from taking the hyperspectral data, a white target using the light set-up and a black image by closing the camera’s shutter were measured. These measurements were taken to eliminate the effects of black (B) and white (W) backgrounds on the hyperspectral images. The calibrated reflectance values (R) of raw sample images (I) were calculated by the following formula:

$${R}_{i}=\frac{{I}_{i}-{B}_{i}}{{W}_{i}-{B}_{i}}$$

where i is the pixel index.

Processing: geometric correction

This process aims to eliminate errors or distortions of the image to adapt the image to an established 3D reference system that enables comparison.

The steps to follow in the geometric correction of an image are

  • control points selection,

  • transformation calculation,

  • image resampling, and

  • process verification.

Here, the correction was made by projecting the hyperspectral image onto the 3D model obtained by photogrammetry, which was subsequently re-projected following an orthogonal plane. The hyperspectral imaging camera had certain limitations regarding the imaging range and resolution. The pixel size was set to 1 mm, and each recording had about 1 m2.

Pigment analysis

The results can be divided into two categories.

  • Fully interpreted information: the end is to create a cartographic representation of pigment by classifying data.

  • Fully uninterpreted information: the end is to create false colour compositions to enhance paintings hardly appreciated by the naked eye.

Minimum noise fraction transformation

Original data had 214 spectral bands; the minimum noise fraction (MNF) transformation determines the dimensionality of the image by separating the noise from the image (Green et al. 1988). It produces orthogonal bands ordered by their information content.

Pixel purity index

The pixel purity index (PPI) aims to locate the purest spectral points of the hyperspectral image. The method is based on the assumption that the most extreme points in the scatter plots are the best candidates to be used as end members. This model offers satisfactory results when the components residing at a sub-pixel level seem spatially separated (Boardman and Kruse 1994).

n-D display

This locates, identifies and groups the purest pixels and the most extreme spectral responses in a data set. The n-D viewer helps to visualise the shape of a data cloud resulting from the plotting of image data in spectral space (with image bands as raster axes). The maximum number of bands to be displayed is 54.

Unsupervised ISODATA (iterative self-organising data analysis technique) classification

The unsupervised classification methods consist of the search for natural groupings of pixels (clusters) in the spectral space defined by the p variables (spectral bands), with no auxiliary information from the analyst, i.e. no learning samples are used.

A new adaptation of the ISODATA (iterative self-organising data analysis technique) method has been applied, according to Bilius and Pentiuc (2020). Parallel factor analysis (Parafac) is a method to decompose multi-way data and analyse the data using non-supervised classification tools. It is an iterative method repeated until some convergence is reached; it is self-organising in locating clusters from minimal initial data.

Visualisation improvement

It consists of applying specific techniques to enhance and group pigments that are difficult to see with the naked eye.

Decorrelation adjustment

Decorrelation adjustment is implemented for generic binary band data stored in 16 bits, as described by Sabins (1986), and is used to enhance the image. It consists of converting the data to the space defined by the main components of the original bands, followed by equalising the data according to the new axes, and finally converting the data to the initial space and combining the bands with the primary colours, red/green/blue (RGB). In this way, the dots are more evenly distributed in the RGB space so that the image will show much higher contrast.

Principal component analysis

Two principal component algorithms were used, one adapted from Richards (1986) and the Karhunen–Loeve transformation, as described by Loeve (1978), which were programmed in the IDL language. Both methods were carried out using correlation matrices, and all bands were included in the transformation.

Minimum noise fraction transformation

The minimum noise fraction (MNF) transformation was implemented as described by Green et al. (1988). It was created to enhance the ability of principal component analysis (PCA) in separating signal and noise components in multispectral images. Unlike PCA, MNF maximises each component’s noise content instead of the data’s variance, resulting in a higher signal-to-noise ratio when reversed. The MNF transformation effectively removes noise from the image, determining its actual size and reducing computational requirements. This linear transformation involves two separate rotations using the main components of the noise covariance matrix to decorrelate and rescale noise in the data and then using the main components from the original image data after the noise has been whitened and rescaled. The final eigenvalues and associated images determine the inherent dimensionality of the data for subsequent spectral processing.

Independent component analysis

Independent component analysis (ICA) can be used in multispectral datasets to transform a set of randomly mixed signals into mutually independent components, according to Hyvarinen (1998, 1999).

These transformation analyses are used for deriving a set of low-dimensional features from a larger group of variables. However, another application is to visualise higher dimensional data and clustering, one of the key data mining methods for discovering knowledge in multivariate data sets. The goal is to identify groups (i.e. similar coverages) of similar objects within a data set of interest.

Geometric and textural characterisation of the rock support and the engravings

Photogrammetry

To create the 3D model, a Sony ILCE 7 S Mark ii was used, with a 35 mm Sony Distagon T* FE 35 mm f/1.4 ZA previously calibrated. The mean flying altitude was 45 cm, and the ground resolution was around 100 microns (Fig. 11). The coverage area was 5 m2, and ten GCP taken from the 3DTLS data scanned targets were used with a mean error of 0.43 cm.

Fig. 11
figure 11

Density of triangles in the photogrammetric model

Stereoscopic coverage has been ensured over the entire area. In the flight to obtain the 3D model, the overlap in the longitudinal direction reached 80%, and the transversal overlap varied according to the coincidence of the axes of the flight (Fig. 12).

Fig. 12
figure 12

View of the main panel with the photogrammetry equipment

Illumination source

An LED Ring CNR-332 with a colour rendering index (CRI) of 95% was used to ensure the faithful reproduction of the colour of the rock art panel. The illuminance is 1900 lx/1 m, dimmable from 3 to 100% and was spectrally calibrated (Fig. 13).

Fig. 13
figure 13

Spectral signature of light at working distance. The vertical axis represents the relative spectral power measured, and the horizontal axis is the wavelength in nanometres where it was measured

Processing

In this stage, the photogrammetric shots are aligned as described by Bayarri et al. (2015). Traditional digital methods were executed in VisualSFM (http://ccwu.me/vsfm/).

During the data capture stage, a control process was applied to the incident light to ensure colour management, keeping the colour and its perception constant across all the relevant processes. A Datacolor SpyderCHECKR colour chart was used to record colour (Fig. 14).

Fig. 14
figure 14

Three measurements of light with Datacolor SpyderCHECKR colour chart: (a) natural incident light, (b) frontal LED ring light, and (c) Philips Photolita lamp

The colour management was based on the consensus achieved through the International Color Consortium (ICC), which focuses on the following:

  • Using colour profiles to describe colour transformations. The colour profile is characterised by gamut, dynamic range and tonal reproduction.

  • The profile connection space (PCS) is a strategy for the transformation of colour between an input and output space.

  • The colour management module is a colour transformation tool that feeds on colour profile information.

The radiometry of the processed images has made effective use of all the bits according to each case, avoiding the appearance of empty digital levels in the case of the 12-bit image.

Analysis strategy and visualisation

A 3D model is several points in space connected by complex geometric entities (i.e. vertex, edges, triangles, curved surfaces etc.) with a texture assigned. Textures consist of overlaying a two-dimensional colour-calibrated image on the surface of the 3D model.

The result of the photogrammetric model has been a 3D model with 12 million polygons and four textures with 8192 × 8192 size (Fig. 15). Three projection planes have been applied to the generated model to have as orthogonal a projection of the elements as possible (Figs. 19 and 20).

Fig. 15
figure 15

Frontal (a), isometric (b) and lateral (c) projection of the main panel based on the photogrammetric model

In addition to the 3D model, a georeferenced orthoimage with 100-micron pixel size and a digital elevation model (DEM) has been generated as a basis for further studies.

As we have indicated in the ‘Introduction’ section, we can group the different ways of visualising these 3D geometric models into two groups:

  1. a)

    graphical solutions within the visual experience

  2. b)

    graphical solutions outside the visual experience

Graphic solutions within the visual experience

These are solutions that deal with our understanding of three-dimensional space from an ‘appearance of volume’. Most of the time, this space materialises in maps of light and shadow. Originally, the hillshade represented shadows and solar radiance levels on the ground. Using Digital 3D, models can be simulated by locating a spotlight or sun in a specific place, which allows us to create a sensation of depth with shadows projected on the model (Fig. 16). They generally improve the visual quality of the 3D model by assigning clarity values to each triangle of the model according to the angle of incidence of the simulated spotlight or sun. So, these maps allow the adjustment of light, opacity, reflectivity, shading and other parameters. This process also allows the derivation of products such as polynomial texture mapping, also known as reflectance transformation imaging (Mudge et al. 2006), for interactive visualisation of objects under different lighting conditions to reveal surface phenomena.

Fig. 16
figure 16

3D hillshaded model of the panel

Graphic solutions outside the visual experience

These graphical solutions let us isolate or analyse certain geometric features of the 3D model and express them in maps of arbitrary scalar values. Such maps are independent of the observer’s viewpoint or their understanding of space. They are accurate metric documents, interchangeable between observers. In fact, both types of solutions complement each other and are often used in the same document.

In our analysis, we have used these two filters.

  • Curvature map of measured data: curvature is the amount by which a geometric object deviates from being flat (Nielsen 2020). The maximum curvature value at each polygon vertex has been calculated and represented.

  • Angle differences between each normal of a polygon vertex and a reference direction normal to the panel (Figs. 17 and 18).

Fig. 17
figure 17

Top: curvature map of measured data. Bottom: angle differences between each normal of a polygon vertex and the normal to the panel

Fig. 18
figure 18

Example of a combination of a hillshaded map and a curvature analysis (left) versus the curvature map alone (right)

Projection planes selection

A set of 3 points for each working area has been selected to define precisely each projection plane.

The zones have been divided (Figs. 13 and 14):

  • anthropomorphic and weapon

    The main anthropomorphic figure and the weapon next to it are engraved on the far right-hand side of the decorated wall. Its ‘median plane’ is slightly oblique to the median major plane of the rest of the decorated panel on its upper bench. It is a subplane oriented towards the outside of the rock. It is appropriate to apply its own projection to it.

  • upper bench

    This is an overhanging rocky bench bordered, at its base, by a natural sub-horizontal crack. It is a surface of 3.4 × 1.9 m, more or less homogeneous on its face side, so it is appropriate to establish a single median projection plane. However, the right end of the surface moves progressively away from this plane of projection, so the outer limit of the larger anthropomorphic figure is slightly distorted (Fig. 19).

  • lower bench

    Unlike the upper bench, with which it borders, the face side of this bench has a homogeneous ‘positive’ or convex inclination and a general pseudo-vertical dip, which allows for a projection onto a common plane without notable distortions (Fig. 20).

Fig. 19
figure 19

View of the projection planes with the location of the points that defined these planes

Fig. 20
figure 20

Front and side views of the three projections used in all images. The position of the projection points is indicated

Parallelism of radiometric spectra in geometric models

The surface represented in the hyperspectral images results from a conical projection that suffers anomalies or deformations regarding the position it should actually occupy in comparison with the true orthoimage obtained from integrating the 3D laser scanner and terrestrial photogrammetry. The geometric corrections are intended to reduce these errors, so the resulting image after these processes keeps, as far as possible, the radiometric values of the initial image and adapts to a chosen reference.

In this case, the process of taking an image as a reference is known as registration.

Basically, the two main reasons for geometrically correcting an image are to eliminate errors or distortions and to adapt the image to a reference system, allowing it to be compared with other elements in the same projection.

Geometric corrections refer to any change in the position of pixels in an image. The changes can be explained by a numerical transformation, the expression of which may be

$${\text{F}}1\left(x,y\right)={\text{G}}1\left(f,c\right)$$
$${\text{F}}2\left(x,y\right)={\text{G}}2\left(f,c\right)$$

where x and y refer to the coordinates of the corrected image, and f and c are the input image coordinates.

The sequence of phases in the geometric correction of an image is as follows.

  1. 1.

    Hyperspectral image acquisition: to obtain RAW data

  2. 2.

    Pre-processing: the result is the hyperspectral reflectance image.

  3. 3.

    Resampling of the image: to theoretically match the cell size in x and Y.

  4. 4.

    Calibration: obtaining reflectances (I) from raw data (I0) by using dark (D) and white (W) references by using the expression I = (I_0 − D)/(W − D)

  5. 5.

    Georeferencing: to convert hyperspectral data into spatial data. The result is a georeferenced 3D-calibrated hyperspectral model. This model will be used for analysis and visualisation to create cartographies and false colour compositions. The stages are as follows:

  6. 6.

    selection of control points,

  7. 7.

    calculation of the transformation,

  8. 8.

    resampling of the image, and

  9. 9.

    verification of the process.

Here, the correction has been carried out by projecting the hyperspectral image onto the 3D model obtained by photogrammetry, which has been re-projected following an orthogonal plane (Fig. 21).

Fig. 21
figure 21

Parallax process between radiometric data and geometric models

Visualising and analysing data derived from hyperspectral imaging and photogrammetry in rock art studies requires the use of specialised software tools and algorithms. The typical steps and tools used for this purpose are as follows.

Photogrammetric 3D models can be generated using various photogrammetric software options such as Agisoft Metashape (Agisoft LLC., St. Petersburg, Russia), RealityCapture (Epic Games, Inc., Cary, NC, USA) or Visual SFM (Developer: Changchang Wu). These models are saved in the Wavefront.obj file format, accompanied by several 8192 × 8192 images in JPEG, PNG or TIFF formats, in the function of the rock art’s extension and the predetermined resolution. Visualisation of the 3D models can be done using tools like Blender (Blender Foundation, Amsterdam, the Netherlands), MeshLab (ISTI-CNR, Pisa-Rome, Italy) or Unity (Unity Software Inc., San Francisco, USA).

Hyperspectral remote-sensing data is stored in the BSQ generic binary format (.BSQ), and the algorithms produce data stored in the same format. Remote-sensing software such as ERDAS Imagine, Exelis ENVI or Fiji can be employed for data calibration and processing. Tools like Matplotlib in Python or specialised hyperspectral visualisation software contribute to the creation of informative visual representations. Only the outcomes of pigment analysis, namely Cartography and False Colour compositions, are stored in standard JPG or TIF formats indexed to the 3D model using GCP and the transformation method. A texture is then generated and indexed to the geometric model, and the results can be visualised with the same tools.

Orthocorrection of 3D data can be performed to generate 2D representations, which can be managed using GIS software such as ArcGIS (Esri, California, USA), QGIS (QGIS Development Team) or GvSIG (Asociación gvSIG, Valencia, Spain).

Integrating hyperspectral and photogrammetric data often requires a multidisciplinary approach, involving collaboration among researchers in fields like remote sensing, archaeology and computer vision (Barrera et al. 2009). The selection of specific tools and algorithms is contingent on the research objectives, the features of the rock art site and the available resources.

Discussion: critical and comparative analysis of data and interpretations

Hyperspectral imaging is now a widely adopted tool in heritage studies and art history, used for the examination of paintings, objects, manuscripts, wall paintings etc. The diverse pigments used, often derived from minerals in ancient artworks, are easily discernible through spectroscopy. Hyperspectral imaging proves valuable in mapping the distribution of pigments and inks, characterising them, identifying binders, recognising techniques, visualising overlaid sketches or figure outlines and studying palimpsests (Fischer and Kakoulli 2006; Le Mouélic et al. 2013; Daniel and Mounier 2015; Daniel et al. 2017; Mulholland et al. 2017; Daveri et al. 2018; Cucci et al. 2019; de Viguerie et al. 2020; Li et al. 2020).

In archaeology, hyperspectral imaging has been applied to investigate the polychromy of Greek sculpted friezes, Egyptian tombs and Roman baths (Alfeld et al. 2018, 2019; Cortea et al. 2021). Despite the advancements of hyperspectral imaging in heritage studies, its application in rock art archaeology is relatively new and remains largely uncommon (Bayarri et al. 2016; Bayarri et al. 2019, 2021a, b). This is generally associated with the challenges with the outdoor implementation of hyperspectral instruments. The conditions pose difficulties in obtaining images that adhere to the prerequisites of uniform illumination and the minimization of optical and geometrical distortions (Alexopoulou et al. 2021).

Integrating hyperspectral imaging and photogrammetry has emerged as a powerful and multifaceted approach, delivering significant contributions to archaeological research and enriching our understanding of diverse sites (Schmitt et al. 2023). This combined method proves instrumental in various parts of analysis, including the characterisation of pigment composition and artistic techniques. It helps with the discernment of faded or weathered elements, unveiling hidden motifs and layers that may be obscured. The integration enables the mapping of mineralogy and subsurface features, shedding light on the geological context of archaeological sites. Beyond the technical aspects, it contributes to a deeper understanding of cultural context and symbolism, unravelling layers of meaning embedded in the artefacts (Defrasne et al. 2023). Additionally, the documentation capabilities of this integrated approach play an important role in conservation efforts, providing detailed records for preservation initiatives. Furthermore, it helps with chronological and stylistic analysis, offering a comprehensive toolkit for archaeologists to explore and interpret the intricate tapestry of human history encapsulated within archaeological remains.

In this paper, we will not detail the chrono-cultural implications of our reading of the panel from an archaeological viewpoint. As can be deduced from the previous sections, the contribution of this study is eminently theoretical and methodological. In other texts, we plan to develop these other interpretative questions. Still, from an eminently descriptive viewpoint, we will address the evidence that our method contributes to the archaeological discussions raised in Peña Tu because of the successive documentation of its icons throughout the twentieth century and the first decade of the twenty-first century. Above all, we will focus on certain ‘hot spots’ that have been described differently according to authors (or according to specific publications by these authors) or whose very existence has been denied. Although it may seem that a large part of our scientific effort merely discusses previous proposals, beyond these validation-refutation exercises, our documentation can detail original descriptions of the panel’s most important and ‘classic’ motifs and show the existence of decorated areas never previously considered. No less important are the possibilities that our method reveals in terms of detecting later additions of pigment and evaluating the state of conservation of the decorated ensemble.

General view

From the viewpoint of the panel as a whole, in terms of general proportions, we start from the disadvantage of lacking data on how the three major revisions of the site were projected onto the plane (Fig. 5). Whether a single projection plane was used or several. As the tracing techniques have also been threefold, freehand (Hernández-Pacheco et al. 1914), copying on a cellophane adhered to the rock (Bueno and Fernández-Miranda 1981) and photographic techniques (Barroso Bermejo et al. 2007), it is possible that, in some proposals, this question has not been considered. Nor do we have references to the horizontal plane to help us explain the inclination of any of the motifs or the angular relationship between motifs. Thus, in the 1914 proposal, the angular relationship between the larger anthropomorphic-sword duo and the rest of the representations to their left does not make sense in a common projection. In this proposal, there is also an error in the position of the motifs on the far left in the possible representation of a quadruped or quadrangular geometric figure. Clearly, this group on the far left has been placed in a lower position than in reality. These anomalies probably arise because the tracer, using a freehand technique, has worked from various positions, partially reproducing the motifs and with little concern for angular rectifications of the whole.

In this sense, the 2007 reproduction is more accurate if we imagine a single projection plane. However, the long-width proportions of the larger anthropomorphic figure are striking, as well as the marked vertical asymmetry of the borders and arches that define its upper part. The authors of the publication consciously reproduce this last detail. They indicate that the motif has been adapted to the edge of the right end of the rocky plane, and for this reason, this area has been delineated in a more ‘ovoid’ shape than the arches on the left side (Barroso Bermejo et al. 2007: 134). We agree in part with this novel assessment, although not to the degree of deformation in their tracing. Although we have no data on the photographic techniques used, we can only imagine that optical aberrations in the lenses used to capture the data explain this marked deformation.

The major anthropomorphic figure

In addition to this disparity in the general proportions of each tracing, we have observed other differences in the internal decoration of this motif. It is difficult to extract the peculiarities of the engraved traces from the 2007 tracing since it uses a shading convention that is sometimes indistinguishable from other conventions dedicated to the natural roughness of the rock. But, the texts of all these publications hardly describe the details and scope of this technique. From our data, we can see that the engraved grooves delineate a shape enveloped by two larger borders down to the foot of the image. Here, only the inner border has horizontal continuity. In addition to these bands, the upper area (which we identify as ‘the face’) is bounded by two other inner bands that stop at the first groove, defining the horizontal levels. In this sense, we do not observe a supposed extension downwards and to the right of one of the two borders defined in ‘the face’, as the 2007 interpretation implies. Referring to the same tracing, we observe neither the lower seventh line of this area organised in horizontal stripes nor the transversal grooves of the outermost border on the left side.

As for the painted areas, our spectral classification leads us to consider several clusters in the ochre tones. To put it colloquially, we observe several instances of repainting. In an effort to clarify, we have grouped all the spectral ‘signatures’ of these tones into three classes. Two of them are associated with thick-stroked, dark-toned varieties of thinned applied paint. We suppose they were applied in prehistoric times. One of them, the clearest, is the one on which the descriptions of the classic tracings have been based. The other tone, which is less obvious, enriches the description of some of the most significant elements of the design. These are features that have remained unpublished until the present day (Fig. 22).

Fig. 22
figure 22

Different results for the area of the largest anthropomorph and weapon: (a) pigment classification, (b) digital elevation model with hillshade overlayed, (c) pigment classification over hillshade model, (d) hillshade model and (e) pigment absorption band

There is a third, more orangey type, which was probably applied with chalk. Although clearly, in some of the strokes, it is difficult to isolate it from other ‘stains’, which we believe to be natural exudations from the rock.

This last type takes us back to an unfortunate recent intervention, specifically from the 1980s. R. de Balbín described it shortly after it took place. A teacher from the surrounding area, with the intention of ‘clarifying’ for his pupils the faded images documented there, traced with chalk several coloured lines in some grooves (Balbín-Behrmann, 1989:29) (Fig. 23). Until now, the extent of this unfortunate intervention was not known. With this work, we have been able to map the areas of deterioration.

Fig. 23
figure 23

RGB photograph of the detail of the engraved and painted grooves of the larger anthropomorph (Sony αSII—FE 35mm F1.4 ZA). The chalk repainting affecting the protuberances inside the grooves can be seen

The most important discrepancies in the classical tracings are in the decorative conventions of what we identify as ‘the dress’, specifically the number of zigzags and the number of parallel, serial strokes that define its design. Thus, while in the 1914 tracing, the inner horizontal stripes are decorated with parallel serial lines; in the 2007 interpretation, the decoration of this area is based on zigzags. The same is true of one of the arched borders on the head.

Although these areas are difficult to interpret, especially those involving lines not highlighted with an engraving groove, our data lead us to believe these areas were decorated with parallel lines.

The significance of these discrepancies concerns a question of ‘archaeological context’, i.e. formal parallels seen in other examples of anthropomorphic representations from the north of the Iberian Peninsula. Thus, researchers are asking to what extent the conventions of representation of these icons are uniform from one region to another. How the graphic models of representation of ‘power’ were exchanged at the dawn of social portraiture. In this sense, the Peña Tu example is integrated into a larger group of anthropomorphic images: the well-known Peña Tu/Sejos/Tabuyo group (Bueno Ramirez et al. 2005) (Briard 1993).

As for the elements defining ‘the face’, the most significant discrepancy between earlier tracings and our data must do, once again, with the painted strokes that are not underlined by an engraving groove. All the reproductions identify two engraved and painted circles, which were interpreted as the ‘eyes’ of the figure. Until now, between these two circular shapes, a vertical line painted in red was depicted, which was interpreted as the ‘nose’ (Hernández-Pacheco et al. 1914) (Bueno and Fernández-Miranda 1981) (Barroso Bermejo et al. 2007). We have been able to document two other painted lines parallel to this one, which reach above and below the circles of the eyes mentioned above. They are interrupted by the circles. They are barely perceptible traces and belong to the second of the clusters of prehistoric red paint established in our tonal classification. It is an area with multiple flaking, preventing clear identification of the remaining pigment. We believe it is also possible to identify the extension to the upper arch of the central line, which classical tracings interpret as the ‘nose’. So, in this space of ‘the face’, we can see a total of three parallel vertical lines of similar thickness that reach, at the top, the inner arch of the first border and, at the bottom, the first of the horizontal lines of the ‘dress’.

The weapon

The image of the weapon is one of the key elements in giving chrono-cultural context to the decoration of the panel or part of it. Although, in some publications, its shape has been interpreted as the silhouette of an anthropomorphic tomb (Fernández-Menéndez 1925), most authors agree that its design concerns a weapon. It is a metal weapon with a more or less developed hilt, probably made of perishable organic materials. Sometimes, it is called a dagger (Bueno and Fernández-Miranda 1981, p. 457), and other times, a sword (Almagro Gorbea 1972, p. 67). Sometimes, physical parallels are found in the early stages of the Copper Age metallurgy in the north of the Iberian Peninsula. Other times, it is shifted to the early stages of the Bronze Age. That means we are talking about a period between 2500 and 2000 BC. In technical terms, we would talk about a copper dagger, the hilt of which is embedded in a metal tang, or a weapon (copper or bronze) with a slightly more complex hilt, in which a system of fastening the hilt with metal rivets cannot be ruled out. This motif has always been observed and described in great detail and has given rise to various discussions about the elements that characterise its design.

The engraving groove has not been the subject of debate in any discussion. At least, successive observers have paid little attention to its description. The type of groove, homogeneous in thickness and U-shaped in section, is similar, although somewhat narrower, to the group of engravings that make up the carved design of the anthropomorphic figure. However, from a purely archaeological perspective, the size ratio (in length and width) between the hilt and the metal blade indicates a type of weapon of relatively large dimensions compared to what is known in the northern third of the Iberian Peninsula (Teira-Mayolini and Ontañón-Peredo 2016). By this, we mean it is more appropriate to speak about a ‘sword’ rather than a ‘dagger’, which would distance us from the early Copper Age weapons mentioned above.

While the technique and features of the line of the groove are not the subjects of debate, the situation is different if we refer to possible traces of paint in the groove itself or inside the surface silhouetted by the groove. Thus, while the groove is unpainted in the 1914 and 1981 publications, the 2007 tracing indicates that it was also painted. However, this circumstance is not mentioned in the text of this publication. In this tracing, the stroke painted in red is represented with a colour convention of a fainter tone than that applied to other areas of the anthropomorphic figure. As a third option, in Blas Cortina’s interpretation, the entire inner surface of the motif in addition to the groove was painted red (Blas Cortina 2010).

We have been able to observe that this engraving groove was also painted several times using different colouring media. We believe that almost all belong to very recent interventions, probably applied with dry-coloured pencils or with coloured chalk. As in the case of the larger anthropomorph, several orange strokes are superimposed on all the others. Undoubtedly, this is once again the mark of the teacher in his regrettable intervention in the 1980s (Balbín-Behrmann, 1989: 29). The spectral traces of recent interventions are so confused that we cannot be sure whether the groove was painted in prehistoric times using a diluted painting technique. It is possible to distinguish at least two types of colour made with chalk (Fig. 24).

Fig. 24
figure 24

Spectral signatures of the traces of pigment outlined with coloured chalks and their distribution in the engraving of the weapon in false colour

However, as previously described, this orange tone is close to various natural exudations of the rock, so as in the mapping of its tonal family, other stains seem to have nothing to do with the action of the human hand. These stains are associated with small natural concavities, which has led to less exposure of these surfaces to weathering. In the engraving grooves, these recent repaintings can be seen as globular remains adhering mainly to the protuberances of the inner micro-relief without impregnating the bottom of the groove. They have been traced with a dry-coloured pencil.

We also have different descriptions regarding the existence of red dots, which could be interpreted as the rivets fixing the hilt to the metal blade of the weapon. The 1914 and 1981 publications describe five points in an arc (Hernández-Pacheco et al. 1914) (Bueno and Fernández-Miranda 1981). The 2007 revision, however, states these points never existed (Barroso Bermejo et al. 2007:137). For Blas Cortina, as mentioned above, the supposed dots are only relics of paint on the inner surface of the engraved design, which was covered with pigment. Other authors, such as Almagro Gorbea, accept the interpretation of the five dots drawn in an arc (Almagro Gorbea 1972) (Fig. 25).

Fig. 25
figure 25

Different tracings of the weapon of the larger anthropomorph compared with the radiometric data obtained in this work. All cases are shown at the same size, ignoring the inconsistencies of the published scales. a Orthoimage with RGB texture. b Absorption band for ochre (600–700 nm). c (Hernández-Pacheco et al. 1914). d (Almagro Gorbea 1972). e (Bueno and Fernández-Miranda 1981). f (Barroso Bermejo et al. 2007). g (Blas Cortina 2010)

We have been able to confirm that there are at least three painted dots in the area corresponding to the grip of the handle, in a position coinciding with that described in the 1914 and 1981 versions, but reduced to three elements. The pigment is similar to that found in the grooves of the anthropomorphic figure; their placement is symmetrical regarding the axis of the engraved silhouette of the sword, and they are in the centre and at the ends of the metal blade. However, this area is dotted with multiple recent, intentional chipping or punctures. We do not rule out the existence of the other two ‘dots’, although we find it curious that only these three have survived, apparently arranged in axial symmetry (Fig. 26).

Fig. 26
figure 26

a, b Detail of the weapon next to the larger anthropomorph. The absorption band for ochre in (b) shows three dots of pigment that we interpret as rivets on the handle

In addition to the remains described above, other pigment stains remain on the inside of the silhouette and on the outside, next to the anthropomorphic figure. We could recognise no structured form in them, but their spectral signature is similar to that of other motifs painted to the left of these two larger elements. In this sense, Blas Cortina’s interpretation is correct, as he sees the carving and painting of these two figures on the far right of the panel as belonging to a second phase of decoration after most of the working surface had been used (Blas Cortina 2010).

Anthropomorphic figures in the central area of the panel

Between the groups of punctuations on the left of the panel and the large anthropomorphic figure with a sword on the far right, we can see a group of anthropomorphic motifs of standard convention, painted in red: a vertical fusiform trunk, limbs drawn in straight, oblique, downward lines and the head as a mere prolongation of the line of the trunk. It is a ‘scene’. This is how the authors of the 1914 tracing interpret it (a dance) and those of the 2007 version (an aggregation scene). However, the two major differences lie in the identification (or not) of lower limbs in almost all the figures and in the different interpretations of the ‘anthropomorphic figure with staff’ or ‘crook’ in the lower part of the group.

The 1914 tracing reproduces these as complete figures, i.e. with all their limbs, upper and lower, almost like a repeated pattern. The 2007 interpretation is very different, showing no lower limbs in any of the icons. Nor is it easy to agree on the number of anthropomorphs. The ‘Dance’, in 1914, consists of six figures, plus the figure with the crook. These are visible in his graphic proposal, and their number is also detailed in the text of the article (Hernández-Pacheco et al. 1914). As mentioned above, this tracing has no graphic convention regarding natural relief or what the authors consider non-prehistoric. Thus, there are no references to the punctuated motifs of the Christianisation crosses nor to other signs executed with the same technique of more dubious interpretation. Their position is not insignificant regarding the problems of interpretation in this area. The two crosses with orbs (globus cruciger) and another simpler one on the right seriously affect the painted remains of four of the six anthropomorphs. Their location in the tracing could have served as a reference for locating these painted motifs. Thus, in the same way as with the distribution of figures on the far left of the panel, some of these anthropomorphs seem to have been displaced from their actual position. This is noticeable in those on the vertical to the weapon, which have been drawn immediately next to it when they are over 20 cm away. Their location is correct in the 2007 tracing (Fig. 27).

Fig. 27
figure 27

Anthropomorphs in the central area of the panel. a Ochre absorption band (600–700nm). b Tracing from 1914. c Tracing from 2007

We can interpret the five figures to the right of the crosses with orbs as anthropomorphs because they all have lower limbs. There are undoubtedly more among the amalgam of traces in that area and those between the abovementioned crosses. The clearest figure is to the group’s left, before the figure with the staff. Both the 1914 and 2007 interpretations describe it as an anthropomorphic figure without lower limbs. From our spectral selections, this figure can be interpreted as a repainting of an earlier one that had ‘legs’. The panel was decorated at various times. Not only new figures were added, but also older motifs seem to have been repainted. We can even say there is an overlapping between figures. This can be seen in the far right anthropomorphic figure, whose strokes interrupt a rectangular sign composed of five vertical strokes and another horizontal one that joins them at the top (Fig. 28).

Fig. 28
figure 28

a, b Central area on the upper bench of the panel. In the ochre absorption band, b two families of red pigments can be seen, which may indicate two different moments of execution of the decoration

Interpreting the character with the staff or crook has also been conflicting. It is perceived as such in the tracings of 1914 and 1981, unlike in 2007, where the staff becomes an anthropomorphic figure in ‘phi’, who comes with another similar one to his right, with arms appearing as a handle (Barroso Bermejo et al. 2007). Both figures are affected, on the left and on the right, by two punctured motifs, a cross and a possible zoomorphic figure, which prevent the respective ends of these figures from being analysed.

We believe it is a human figure in a dynamic pose, with its legs open and arms almost horizontal. In this sense, although it reproduces a standard convention, it is a variant for this graphic group. The right arm seems to cling to a vertical linear element, with a circle or ring at its end. This would identify the crook referred to in the 1914 and 1981 interpretations. However, in reality, the design does not end with this linear line. Next to it, towards the left, up to the alteration of the picket of a cross, we can see other lines clearly linked to the ‘staff’. They seem to have three points arranged vertically and other lines with a horizontal tendency. In reality, this reasoning related to the recognition of staff is based on describing what we think we understand while ignoring the elements for which we find no context in our visual experience (Fig. 29).

Fig. 29
figure 29

The ‘figure with staff’, in the central area of the panel. a Ochre absorption band. b 1914 tracing. c 2007 tracing. d 2010 tracing

Groups of dots and other elements in the centre and far right of the panel

The central and left-hand area of the panel is where we can identify the fewest narrative ideograms. It is an area full of simple geometric elements, mere accumulations of dots, in which researchers have sought to recognise a specific formal order. Two motifs of a somewhat more elaborate design stand out, which we can describe as quadrangular geometric figures and anchor-shaped figures. Perhaps it is the most neglected area in terms of elaborate tracings or textual interpretations. In the 1914 version, we see sets of roughly misplaced strokes. From the quadrangular motif (or quadruped) to the far left of the panel, all the figures depicted are displaced downwards, leaving the horizontal crack decorated with dots above them. In reality, these figures are above this line.

This area at the end is correctly placed in the 2007 interpretation. However, the groupings of dots here become dotted ‘stelliform motifs’, small in size but similar in silhouette to the large anthropomorphic figure at the other end (Barroso Bermejo et al. 2007). Likewise, the quadrangular motif, identified as a ‘schematic animal’ in 1914 and 1981, disappears from his tracing. Our data reveal that this area on the far left has only a nuclear set of dots with no recognisable geometric form. Below this, the end of the crack is decorated, confirmed in this case, with a line of dots. We can also see the ‘animal’ figure, which we prefer to describe as a quadrangular geometric figure (Fig. 30).

Fig. 30
figure 30

Group of dots and other ideomorphs on the far left of the panel. a 1914 tracing. b 2007 tracing. c Ochre absorption band

The most numerous grouping of dots is found between this group of figures, at the end of the panel, and the area of the figure with the staff. For this area, both the tracings published in 1914 and 2007 and the descriptions in the texts of their publications, in addition to that of 1981, differ notably. The 1914 interpretation sees internal orders based on curved dotted lines. It even mentions a concave shape (Hernández-Pacheco et al. 1914) which is recognisable in the tracing but which, in reality, acquires its appearance from chipping off the wall. It seems that the writers of the text have described these areas based on Cabré’s tracing (the real draughtsman) and not based on direct observation of the panel. As we have indicated, this tracing shows no modern alterations, whether natural or man-made.

In the 2007 tracing, the density of dots is lower and irregular, with abundant ‘bald spots’ between groups. This is an altered area due to micro-organisms, which have generated a dark biofilm on the rocky surface. This makes it difficult to identify the remains of paint with the naked eye or with commercial RGB sensors. Here, the equipment used in our research has made it possible to document a circular area of about 45 cm in diameter, in which more than a hundred red dots of homogeneous size and dispersion are distributed. They do not seem to be arranged in recognisable shapes. This circular ‘patch’ of dots is delineated by the horizontal cracks in the sandstone strata, which are also decorated with dots. To the right of this area, there is an anchor-shaped figure, damaged on the lower part by the puncturing of a modern cross. This motif was already recognised in 2007 with this label (anchor-shaped), assigning it an anthropomorphic character. This figure comes with another, on the right, which the authors of the 2007 tracing also interpret as anthropomorphic. The lower part of this motif is also changed by more recent flaking. If we accept this anthropomorphic character of the two figures, their graphic convention will differ from the rest of the figures given as such on the panel. They are two designs with a larger side and vertical symmetry, like the other anthropomorphic figures, but with the upper area resolved as an inverted anchor. At least, this is clear in the figure on the left. The immediate one, on the right, is very deteriorated (Fig. 31).

Fig. 31
figure 31

Dots and ancoriform on the left side of the panel. a Ochre absorption band. b Tracing from 1914. c Tracing from 2007

Traces of paint in the lower back

Finally, the inspection of the lower rock bed of the panel has brought us interesting insights. This has never been reproduced in the various tracings. There is only one brief reference to it in the 1914 publication in which the existence of ten dots below the figure with the crook is noted and is also reflected in his tracing. This publication also refers to ‘vague spots’ (Hernández-Pacheco et al. 1914).

On this lower bed, which is slightly more convex than the upper stratum, almost 1 m high, we have observed various coloured patches. Some are sharper than others. Some have a specific delineated geometric appearance (Fig. 32). The dots referred to in the 1914 publication, below the figure with the staff, are now seen as two small parallel lines between two horizontal cracks decorated with dots. Below these two cracks, a closed oval shape can be seen, 8cm in length, and perhaps subdivided by other traces inside it, although we may read these details by using some flaking to our advantage (Fig. 33).

Fig. 32
figure 32

Lower bench of the panel. a Projection to the plane of the RGB texture. b Ochre absorption band (600–700 nm)

Fig. 33
figure 33

Detail of the lower bench of the panel. Ochre absorption band

In a somewhat riskier reading, we could glimpse a closed rectangular shape, about 70-cm long (vertically), with a horizontal internal division located slightly to the left of the abovementioned shape. Other patches of less well-defined shapes can be seen spreading to the foot of the bed.

Another circumstance which supports the idea that this lower bed had also been decorated is the presence of picket crosses from the Christianisation of the site, as seen among the motifs on the upper panel. Here, there are at least two simple crosses, one of which is of the ‘Jerusalem’ type. This leads us to believe that the figures painted in this area were visible to the naked eye in relatively recent historical times.

Conclusions

The documentation strategy undertaken in Peña Tu shows two main points of interest. One is methodological, and the other is technical.

On the one hand, we have organised a coherent set of quantifiable data, which enables other observers to verify our results. This also lets us reproduce the model in successive examples for comparative purposes. These two circumstances (verification and controlled reproduction) separate this eminently digital model from classical or ‘analogue’ strategies, in which we can distinguish only two single milestones that define its documentation process: the initial physical reality and the final, intrinsically interpretative, tracing.

In this sense, the documentation of sites with Prehistoric Graphic Expression is on par with other processes in archaeology. In recent years, the generalisation of computer applications has been aimed to a large extent, at measuring those variables of analysis which, in an earlier ‘analogue’ environment, were impossible to systematise. This is especially true in studies where geometry or space are relevant, such as landscape analysis (Figueiredo Daves and Barrocá Faccio 2022; Nakoinz and Knitter 2016) and in micro spatial analysis, in research on human behaviour ‘intrasite’ (Clark and Gingerich 2022; Gopher et al. 2016).

However, research on Prehistoric Graphic Expression has a differentiating circumstance because the form itself (its geometry) has a semantic content. They are designs constructed to be seen and to be read.

This means that the documentation process is confronted by a ‘metalanguage’ belonging to the data (May 2016: 172; Morphy 2005:51), sometimes forcing us to begin our process with subjective choices related to our ability to recognise icons or signs (Beaune 2015). This is more evident with paintings (pictorial stains) than engravings, which we can characterise, with solid criteria, the grooves made by the human hand.

However, we have correctly established two documentation strategies for painting and engraving. The first is eminently radiometric, and the other is geometric. This has also involved an exhaustive parallax control method between the two data sets.

Using a hyperspectral sensor has made it possible to differentiate several families of pigments that can be interpreted as diachronic moments of repainting, at least two of which belong to prehistoric times. With this technique, we have also provided data and obtained complex graphic descriptions of specific icons, which have been conflicted throughout the history of research at the site. This concerns the number and shape of the anthropomorphs, the distribution of the sets of dots and the ideographic interpretation of some of the motifs (Fig. 34). Although it is not the aim of this paper, some of these new descriptions have chronological implications, for example, the combination of painting and engraving found in the depiction of the weapon next to the larger anthropomorph. The dots of red paint on the proximal end of the weapon’s blade can be interpreted as handle attachment rivets, indicating a weapon engineering technique that can be chronologically referenced. Some historical ‘readings’ of this icon refer to five painted dots, but in others, this particular detail does not appear. The aim is to determine whether the weapon belongs to the time horizon of the earliest known copper daggers in the northern third of the Iberian Peninsula or to a later development, in the early Bronze Age, of larger weapons made of alloyed metals. The latter used rivets at the proximal end of the blade to fix the metal blade to the handle. The chronological discussion is very rich and suggestive in terms of the typologies concerned, but poor in terms of the actual weapons found in archaeological sites in these central Cantabrian valleys.

Fig. 34
figure 34

General view of the upper bench of the Peña Tu panel. The projection plane of this view is detailed in the ‘Case study: Peña Tu’s rock (Puertas de Vidiago, Asturias, Spain)’ section. The false colour shows the ochre spectral band (600–700 nm), including the recent alterations made with chalk. The rocky support and the engravings are shown as a light/shadow map

In addition, we have documented new areas with pigment in the lower bench of the quartzite outcrop, whose existence, although mentioned merely as a possibility in the first descriptions (Hernández-Pacheco et al. 1914, p. 7), has not been noted in more recent publications (Fig. 35).

Fig. 35
figure 35

General view of the lower bench of the Peña Tu panel. The false colour shows the ochre spectral band (600–700 nm). The rocky support and the engravings are shown as a light/shadow map

No less important are the parts that refer to the state of the panel’s conservation. With the technique used, we have mapped the extent of the chalk repainting from the 1980s. Although this issue is mentioned in some works (Balbín Behrmann, 1989, p. 29), it has never been possible to obtain a concrete image of such deterioration (Fig. 36).

Fig. 36
figure 36

Image with the dispersion of the dyes attributed to recent repainting with chalk. The rocky support is expressed as a light/shadow map

We have established the relationship between painted strokes and engraving grooves on the larger anthropomorph, which is poorly defined in earlier versions.

To logically conclude the chapter and in our position as prehistorians, we cannot end this study without providing an interpretative or ‘third order’ image (Montero Ruiz et al., 1998, p. 157) of, at least, the major anthropomorphic figure. An image that brings together measured but subjectively selected variables.

Here, it is a matter of tonally matching the two families of ochre paint and their parallax with the rock relief, expressed as a light/shadow map. This convention of the support is the most intuitive for understanding the relief (concave and convex areas). But, as we have explained in the ‘Introduction’ section, as it is a graphic solution within the visual experience, it is difficult to measure (Fig. 37).

Fig. 37
figure 37

Interpretative or ‘third order’ image (Montero Ruiz et al. 1998:157) of the larger anthropomorph and the weapon. Projection of the ochre absorption band shown in false colour, with arbitrary cleaning of the backgrounds. This pigment band has been referenced on the photogrammetric model of the rocky support, expressed as a light/shadow map

The documentation of the Peña Tu site has consolidated, in our experience as archaeologists, a methodological model and a technical praxis. Given the evolution of data capture equipment and analysis software, the latter will obviously have an expiration date. We hope that the work will interest other observers/analysts.