Virtual unrolling and deciphering of Herculaneum papyri by X-ray phase-contrast tomography

Abstract

A collection of more than 1800 carbonized papyri, discovered in the Roman 'Villa dei Papiri' at Herculaneum is the unique classical library survived from antiquity. These papyri were charred during 79 A.D. Vesuvius eruption, a circumstance which providentially preserved them until now. This magnificent collection contains an impressive amount of treatises by Greek philosophers and, especially, Philodemus of Gadara, an Epicurean thinker of 1st century BC. We read many portions of text hidden inside carbonized Herculaneum papyri using enhanced X-ray phase-contrast tomography non-destructive technique and a new set of numerical algorithms for 'virtual-unrolling'. Our success lies in revealing the largest portion of Greek text ever detected so far inside unopened scrolls, with unprecedented spatial resolution and contrast, all without damaging these precious historical manuscripts. Parts of text have been decoded and the 'voice' of the Epicurean philosopher Philodemus is brought back again after 2000 years from Herculaneum papyri.

Virtual unrolling and deciphering of Herculaneum papyri by X-ray phase-contrast tomography

R P

I. Bukreeva,, A. Mittone, A. Bravin, G. Festa,,, M. Alessandrelli, P. Coan,,V. Formoso,, R. G. Agostino,, M. Giocondo, F. Ciuchi, M. Fratini, L. Massimi,A. Lamarra, C. Andreani,,, R. Bartolino,,, G. Gigli, G. Ranocchia & A. Cedola

Herculaneum is the unique classical library survived from antiquity. These papyri were charred lies in revealing the largest portion of Greek text ever detected so far inside unopened scrolls, with unprecedented spatial resolution and contrast, all without damaging these precious historical

The Herculaneum papyri were discovered in the Roman Villa dei Papiri during modern excavations of the Herculaneum site in Campania (Italy) in the middle of 18th century. This unique library survived from antiquity contains valuable works by Greek philosophers, such as Epicurus, Chrysippus and Philodemus, in particular an impressive amount of extensive treatises by Philodemus of Gadara, an Epicurean philosopher of 1st century BC1,2.

Since their discovery, numerous eorts have been made to open Herculaneum papyri and read the precious works hidden inside them and several destructive and non-destructive techniques were employed. Synchrotron X-ray phase-contrast tomography (XPCT) non-destructive technique3,4, originally introduced mainly in bio-medical imaging, in the present study has been optimised and combined with versatile algorithms developed to deal with extracting distorted and corrugated papyri. These are summarized in Supplementary Information (SI). Our ndings provide new insights into the Herculaneum collection, going far signicantly beyond a previous feasibility test experiment which showed the benet of XPCT applied to rolled-up papyri in imaging letters5.

This paper shows our success in virtual-unrolling two unopened papyrus rolls, i.e. PHerc. 375 and PHerc. 495, both stored in Naples National Library. A phantom sample, composed of high quality sheets of bare papyrus paper, was constructed and carbonized in order to determine the capacity and reliability of XPCT technique and the viability of the numerical algorithms. The phantom sample was earlier examined as an archaeological object in both geometry and material using XPCT. Subsequently, both papyrus rolls, PHerc. 375 and PHerc. 495, were measured

Consiglio Nazionale delle Ricerche, Consiglio Nazionale delle Ricerche, Consiglio Nazionale delle Ricerche, Istituto

SCIENTIFIC REPORTS

1

www.nature.com/scientificreports/

Figure 1. PHerc. 375. (A) virtual-unrolling; (B) textual portion.

through the same technique. In order to reveal the hidden text inside the original unopened papyri and the phantom sample, a new set of numerical algorithms for the virtual-unrolling were developed in house (see SI). These were tailored to adapt to the complex and diversied amorphous 3D deformed materials composing the rolls. The test performed on the phantom guided us, allowing to optimize and integrate the set of algorithms.

The result was a digital virtual-unrolling of the unopened papyrus rolls with unprecedented resolution and contrast. Detailed tomographic images of the inner complex structure of the layers along the full longitudinal axis of the bookrolls were obtained (see SI). The optimized phase-contrast 3D tomography distinguished areas of dierent densities: those composed of amorphous carbon-based ink were dierentiated from that of the papyrus substrate, made of ordered carbon bres. The virtual-unrolling procedure, despite the complex organization in the inner portions of the scrolls, allowed to reveal Greek sequences of letters and words and, for the rst time, read and decipher expressions and textual portions as well as a coronis, a symbol used to mark the end of a textual section, a chapter or a book in classical papyri.

Dierent numerical procedures have recently been proposed to provide6,7 a digital version of the text hidden inside rolled-up papyri starting from their non-invasive acquisition. Nevertheless the structure of carbonized Herculaneum papyri is quite dierent from that of a regular roll where a standard virtual-unrolling algorithm is applicable. In recent years specic soware has been developed8 without achieving any reading of the text inside Herculaneum-papyri. Our success lies in combining XPCT with tailoring the basic principles of the virtual-unrolling to the Herculaneum papyri and, consequently, virtually open, read and decode the largest portion of text, hidden inside unopened papyri, ever detected so far. This success is based on three main facts: i) the experience gained with the experiment on the phantom; ii) the optimal experimental conditions which provided, to the best of our knowledge, 3D images of the inner structure in carbonized Herculaneum papyri of the best quality ever obtained; iii) the very intricate internal arrangement of chaotic bundles of layers still presented sets of similarly bent, splayed and twisted layers. Once the XPCT provided the 3D volume of the bookrolls, virtual-unrolling was performed. Upon partitioning the roll in subsets of similarly bent and twisted layers we were able to apply to each one attening procedures for each specic distortion of packed sheets (see SI).

Although the procedure is very elaborate and time-consuming, it is remarkable that the explorable partial area of the virtually unrolled portion of the text was eectively unlimited. This adaptive procedure, allowing the complete unrolling of the rolls, is the propaedeutic tool for the future reading of entire columns of text, possibily paving the way for the rst edition and interpretation of the works contained in unopened carbonized Herculaneum papyri.

Figure1A shows the sequence of the images obtained during the unrolling procedure of PHerc. 375. From le to right of Fig.1A, an increasingly large portion of at region appears, including an 11mm large text of more than three lines. Its position inside the bookroll is shown in Fig.1B. Figure2A shows the sequence of the images obtained from the unrolling procedure of PHerc. 495. Figure2B discloses the coronis with an exceptional image quality.

The revealed text was written on the papyrus side where the writing runs along the bres (recto). Both rolls were ascribed on palaeographical grounds to Philodemus of Gadara, a remarkable philosopher and poet who studied in Athens before moving to Rome and Campania. In terms of textual extension, this assignment goes far signicantly beyond that revealed in the previous XPCT study5. The handwriting was found to be a known typology in the Herculaneum library. As usual in classical antiquity, this was in capitals, with no word spacing,

SCIENTIFIC REPORTS

2

www.nature.com/scientificreports/

Figure 2. PHerc. 495. (A) virtual-unrolling; (B) image of coronis.

breathings, accents or lectional signs. The virtual-unrolling allowed i) to reconstruct a complete Greek alphabet for each bookroll (see Figure 5S in SI); ii) to detect several Greek sequences of letters and words, reported in Fig.3A-B, D-E le) to read for the rst time, expressions and portions of text in Fig.3C, E right-G and Fig.1; iv) to identify, also for the rst time, a textual sign (coronis), reported in Fig.3H le and 2. These all are clearly distinguishable from the papyrus substrate. Using the virtual-unrolling it was possible to perform the reconstruction of the largest portions of text ever detected in unopened papyrus rolls. Three distinct sequences of letters, i.e. ]

, ]c [ and ]c [ were detected in PHerc. 375. These are displayed from le to right in the top row of Fig.3A; three words, i.e. (tri), (perie) and ]c (peistheien), are shown in Fig.3B. In the case of PHerc. 495, the sequences ][ (ep), ]c[ (st) and ][ (eche) are shown in Fig.3D. Furthermore, it was possible to isolate the word ][ (teipomen) (le of Fig.3E) and the expressions] [ (en gar) (right of Fig.3E) and] [ (tel bma) (Fig.3F). This last expression is formed by the second component stem of a compound adjective in c, c or, less probably, the plural noun (from c, end or purpose) and the noun , tribune, in either the singular or the plural. The latter means the raised place from which orators spoke in a public assembly or in a lawcourt and is used twice in Philodemus On Rhetoric (col. 56, 8 Sudhaus II; col. 264, 21 Sudhaus I). The coronis is clearly shown in Fig.3H le beside a known example from an unrolled portion of PHerc. 1008 (Fig.3H right). This latter is imaged by means of multispectral technique.

An example of textual portion is reported in Fig.3C, for PHerc. 375, and Fig.3G, for PHerc. 495, together with a proposed textual reconstruction. In the former case, the verb ] might refer to a resolution taken by, e.g., a political assembly. In the latter, the focus seems to be on either civil or political matters. It must be stressed that the dierent size and the deformation observed in some letters as well as some of the images is due to the distortion produced by the irregular surface of the papyrus substrate. Finally a coronis, namely a typical sign used by scribes in antiquity to signal specic articulations of the text, was revealed in PHerc. 495 (Fig.3H le). Full details of the textual reconstructions can be found in SI. In terms of spatial resolution and contrast, the present virtual-unrolling permitted to fully exploit XPCT techniques.

The large number of letters and texts revealed by the virtual-unrolling made it possible to analyze the handwriting within each unopened roll. Despite the distortions of the letters caused by the irregular surface of the layers, detailed paleographical description was achieved. Full details of it can be found in SI. From a comparative analysis with dierent writing groups inside the Herculaneum collection, both unopened rolls could be attributed with a high degree of probability to the Epicurean philosopher Philodemus of Gadara. In particular, PHerc. 495 might possibly represent a book of his twenty-books large treatise On Rhetoric9, either one partially known to us or a completely unknown one. In addition, our study revealed unexpected events historically experienced by the rolls, equally reported in SI.

These results demonstrate the potential of XPCT technique in reading inside rolled-up carbonized Herculaneum papyri. The investigation will enable our current knowledge of Hellenistic philosophy and classical literature to be expanded, with expected impact in elds such as papyrology, palaeography, classical philology and history of ancient philosophy.

Methods.

X-ray Phase Contrast Tomography imaging. Standard X-ray tomography is based on absorption and it is a well known tool for imaging the internal structure of thick objects with hard X-rays. For low-absorption materials (like Carbon-based objects) small attenuation in the sample produces low contrast in the images. Thus standard absorption tomography is an unsuitable technique to discriminate details of similar densities in, for example, distinguishing carbon bre-based papyrus foil from carbon-based ink used in writing on the papyrus. In the latter case, as shown in this study, a better contrast was achieved by imaging the phase modulation induced by the object in a coherent or partially coherent beam3,4. Several experimental approaches exist for detecting X-ray phase contrast. A simple yet eective phase-contrast method for hard X-rays is based on in-line imaging aer free-space propagation.

Experiment. A series of experiments were carried out at ID17 of the European Synchrotron Radiation Facility (ESRF) in Grenoble (F), using a free-space propagation set-up. The imaging detector is a Fast Readout Low Noise (FReLoN) 2k charge-coupled device camera (CCD) connected with an X-ray optics determining an eective pixel size of 47 47 m2. The nal spatial resolution in the images is about 100 100 m2. The eld of

SCIENTIFIC REPORTS

3

www.nature.com/scientificreports/

Figure 3. Sequences of letters, words, textual portions and coronis, revealed in papyrus rolls through the virtual-unrolling. (AC) PHerc. 375 (top panel) and (DH) PHerc. 495 (bottom panel). Row H shows the coronis: on the le, the image obtained by our virtual-unrolling and, on the right, a known example from an unrolled fragment of PHerc. 1008. The latter is imaged by means of multispectral technique. By permissionof Ministero per i Beni e le Attivit Culturali e del Turismo (Steven W. Booras Biblioteca Nazionale Vittorio Emanuele III, Napoli Brigham Young University, Provo, USA). All rights reserved.

view (FOV) was limited to about 1 mm (vertical) 75 mm (horizontal). A double-silicon (111) crystal system was used to monochromatize the incident X-ray beam. Samples were located at about 10 meters from the CCD camera and they were rotated around the vertical axis parallel to the longitudinal axis of the rolls.a) multi-energy scans were performed on the phantom sample in order to distinguish between absorption and phase contributions: incident X-ray beam was monochromatized at three dierent energies: 80, 52 and 30keV. For each energy, 1400 radiograms were recorded for each scan, covering a total angle range of 180, with acquisition time of 0.05seconds per point. Results of these experiments allow to read and decipher letters and numbers previously written inside papyrus roll and to optimize the experimental conditions. Because of the ow of CO2 during the combustion severe and diversied deformations developed within the bulk of the phantom papyrus roll. This complex inner structure was ideal to develop and test the new numerical algorithms for the virtual unrolling and to tailor and optimize them to reveal the hidden writing accordingly. From the multi-energy experiment on the phantom sample it was also possible to record and to test the capability of distinguishing between the high signal due to the Pb in the ink from the low signal due to the other elements composing the ink.b) The XPCT experiments on the PHerc. 375 and PHerc. 495, were performed using a monochromatized incident X-ray beam with an energy of 73keV. The choice of the energy is one of the essential conditions to achieve a suitable contrast in the nal recorded images. The choice of the incident energy was selected following dierent

SCIENTIFIC REPORTS

4

www.nature.com/scientificreports/

guidelines: a) the experience gained with the phantom, which indicated the range of energies where contrast would have been higher; b) the measurement of the sample thickness which determined the minimum energy to be used; c) nally the optimizion was based on the maximization of the contrast transfer function3 taking into account the xed sample-detector distance and the pixel size of the CCD.

The two bookrolls, PHerc. 375 and PHerc. 495, were placed in cylindrical containers made of Plexiglas, specically fabricated for the experiment, and safely mounted on a rotating plate at the sample position. About 2000 radiograms were recorded for each scan, covering a total angle range of 180, with acquisition time of 0.05seconds per point. Due to the limited FOV, in order to scan the whole scroll, several tomographic measurements were acquired at dierent vertical positions. The total data-acquisition time for each bookroll was about 5hours.

In the rst step of data analysis, the phase retrieval algorithm proposed by Paganin et al.10 was applied to all projections of the tomographic scans using a modied version of the ANKAphase code11. Using a PyHST2 code12, the recorded images were tomographically reconstructed. As a result a stack of 2D images representing the density distribution at dierent depths inside the sample was obtained. By means of a commercial 3D rendering soware the 2D images of the stack were composed to build the whole 3D volume (see Figs 3S and 4S in SI). The inner structure of the whole reconstructed rolls shows the various papyrus layers to be diversely bent and twisted at the dierent depths inside the volume (see movie S1 in SI).

Image treatment. 3D rendering and segmentation of the impurities were carried out using the Visage Amira 6 soware (interactive system from FEI) and VolView, an interactive system for volume visualization from Kitware Inc. (see SI). The 3D structure of the papyri was obtained through a sort of segmentation procedure. The images were generated by associating colors to the dierent gray-scale ranges in the 3D model.

Virtual Unrolling. In the last years virtual unrolling problem using volumetric scanning in virtual restoration and preservation of ancient artefacts like parchments or papyri has been an extremely active area of research. The use of X-ray Computed Tomography (CT) and X-ray Micro Computed Tomography for the data digitalization has encouraged the development of restoration algorithms. Dierent numerical procedures have been proposed to provide a digital version of the text hidden inside rolled-up scrolls by starting from non-invasive acquisition.

The procedure of virtual unrolling can be divided into three main steps: volumetric scanning, segmentation, layered texture generation and restoration.

Volumetric scanning. In this stage the experimental data are acquired and their digitalization is performed.

Segmentation. The process of segmentation is used in order to segment a voxel set into two categories: material voxels versus empty space.

Layered texture generation. This is the stage where the surface modelling is performed.

Digital restoration. The surface attening and unrolling can be interpreted as an isometric mapping (i.e. it preserve distances, which minimises text distortion in the parchment) from 3D to 2D images.

To solve these problems dierent approaches were proposed. One of the most promising is described by Seales and colleagues6,8. The authors have developed soware that integrates functions of attening and unrolling based on mass-spring surface simulation. The algorithms proposed by O. Samko et al.7 allows to solve the problem of touching points between adjacent sheet layers.

The internal structure of the Herculaneum papyrus is quite dierent from the regular roll where a simple virtual unrolling was applicable. However the very intricate internal arrangement of chaotic bundles of layers still presented sets of similarly bent, splayed and twisted layers. The central part of the scroll is better preserved and the layers of papyrus are separated and loosely rolled up. The outer part of papyrus is tightly scrolled and stacked. It is evident that the virtual revealing of the text requires dierent approaches and has its own challenges. From the experiment performed on the prototype we learned that the letters can be deformed due to thermal and pressure impact and due to irregular surface of the bres. A priori it is not possible to know the letters aspect and even whether the text is present in the analysed portion of rolls. In principle the text could be destroyed or could not appear due to the wrong position in the papyrus foil.

We performed semi-manually procedure to reveal the hidden text based on computed algorithms exploiting Matlab codes ImageJ macros and dierent available commercial computer soware. Once the phase contrast

tomography provided the 3D volume of the bookrolls, to reveal the Greek text hidden inside, a multistep procedure has been settled:1. Accurate optimization of the range of grey-levels has been done to individuate the inner structure of papyri with preserving all information about text.

2. Selection of areas with packaging of sheets with similar orientation.3. Flattening of packed sheets, using dierent strategies including transformation of coordinates.4. Filtering and segmentation of individual sheets in 2D slices.5. Surface modelling.6. Segmentation based on threshold to reveal the text

References

1. Zarmakoupi, M. The Villa of the Papyri at Herculaneum: archaeology, reception, and digital reconstruction Vol. 1 (Walter de Gruyter, Berlin, Germany, 2010).

2. Sider, D. The library of the Villa dei Papiri at Herculaneum. (Getty Publications, Los Angeles, USA, 2005).

SCIENTIFIC REPORTS

5

www.nature.com/scientificreports/

3. Bravin, A., Coan, P. & Suortti, P. X-ray phase-contrast imaging: from pre-clinical applications towards clinics. Phys. Med. Biol 58, R1 (2012).

4. Lagomarsino, S. et al. Phase contrast hard x-ray microscopy with submicron resolution. APL 71, 25572559 (1997).5. Mocella, V., Brun, E., Ferrero, C. & Delattre, D. Revealing letters in rolled Herculaneum papyri by X-ray phase-contrast imaging. Nat. Commun. 6, 15, doi: 10.1038/ncomms6895 (2015).

6. Yun Lin and W. B. Seales, Opaque document imaging: building images of inaccessible texts. presented at Tenth IEEE International Conference on Computer Vision: Beijing, doi: 10.1109/ICCV.2005.182 (2005, October).

7. Samko, O., Lai, Y.-K., Marshall, D. & Rosin, P. L. Virtual unrolling and information recovery from scanned scrolled historical documents. Pattern Recogn. 47, 248259 (2014).

8. Seales, W. & Delattre, D. Virtual unrolling of carbonized Herculaneum scrolls: Research Status (20072012). Cronache Ercolanesi 43, 191208 (2013).

9. Ranocchia, G. Nuove acquisizioni sulla struttura del c di Filodemo. Un trattato in almeno venti libri. Zeitschri fr Papyrologie und Epigraphik, accepted for publication (2014).

10. Paganin, D., Gureyev, T. E., Pavlov, K. M., Lewis, R. A. & Kitchen, M. Phase retrieval using coherent imaging systems with linear transfer functions. Opt. Commun. 234, 87105 (2004).

11. Weitkamp, T., Haas, D., Wegrzynek, D. & Rack, A. ANKAphase: soware for single-distance phase retrieval from inline X-ray phase-contrast radiographs. J. Synchrotron Radiation 18, 617629 (2011).

12. Chilingaryan, S. et al. A GPU-based architecture for real-time data assessment at synchrotron experiments. IEEE Trans. Nucl. Sci. 58, 14471455 (2011).

Acknowledgements

The authors wish to thank Luigi Nicolais (CNR), Glenn Most (SNS, Pisa), Daniel Koger (Lindsey Wilson College, Kentucky, US), G. Neville Greaves (University of Cambridge, UK) for valuable discussions and revisions of the present study, the ESRF Directorate for having granted the beamtime, the Biblioteca Nazionale Vittorio Emanuele III of Naples, for lending us the samples (in particular, the officers Soa Maresca and Vincenzo Boni), Eugenio Amendola (CNR-IPCB), for valuable assistance in the preparation of the containers for the samples, Carlo Ionta for technical assistance in data-analysis, Luigi Verolino (University of Naples Federico II) and Gaetano Campi (CNR-IC), for preliminary discussions.

Author Contributions

A.C. proposed and coordinated the experiment, G.R. reconstructed the Greek texts and wrote the historical account, M.A. made the palaeographical analysis, R.B. coordinated the design and experimental measurements on the papyri puppet, R.G.A. and M.G. developed the concept of the papyri puppet and realized them. I.B. developed the algorithms for the virtual-unrolling procedure and A.M. performed the tomographic reconstruction. A.B, P.C., F.C. and V.F contributed to the experiments and the data collection at the ESRF beamline, G.F., M.F., L.M. designed the figures. C.A., R.B., A.L., G.G. contributed to the theoretical discussion and the writing of the manuscript. All authors discussed the results and contributed to the redaction of the manuscript.

Additional Information

Supplementary information accompanies this paper at http://www.nature.com/srep

Competing nancial interests: The authors declare no competing nancial interests.

How to cite this article: Bukreeva, I. et al. Virtual unrolling and deciphering of Herculaneum papyri by X-ray phase-contrast tomography. Sci. Rep. 6, 27227; doi: 10.1038/srep27227 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

SCIENTIFIC REPORTS

6