![](https://cer.blogs.auckland.ac.nz/files/2018/05/CeRTitleOtherWay-1jxwn5m.png)
A new ‘stratigraphy’: interpreting object relationships with 3D point densities
Dr Joshua Emmit, Rebecca Phillips, Simon Holdaway, Department of Anthropology; Sina Masoud-Ansari, Centre for eResearch
Introduction
Recording stratigraphic layers as a means for grouping archaeological phenomena into units dates back to the foundation of archaeology. Part of the craft of archaeology involves recognizing changes in sediment colour and composition, and the resultant need to change stratigraphic layer designation. Such changes are illustrated in stratigraphic diagrams that form a key part of the majority of archaeological site reports. Artefacts are then often assigned to these units which are used to understand artefact deposition, assemblage composition, and ultimately human behavior. Despite the history and centrality of stratigraphy to archaeological practice, stratigraphic recording has seen surprisingly little innovation and scrutiny. While some recent attempts have been made to rectify this, drawing stratigraphy continues to be as much an art as science. Learning how to draw stratigraphy is often seen as part of an ‘apprenticeship’ in archaeological methods. We collaborated with the Center for eResearch to propose an alternative to traditional stratigraphy by presenting a method for an object-oriented approach that utilizes 3D point densities to identify groups of objects.
Methods and materials
Some attempts have been made to apply three-dimensional (3D) point densities to archaeological interpretation, but have failed to be adopted as a mainstream analytical technique. In many cases, a lack of appropriate data may be the cause of this, but additional considerations are the relative difficulty of constructing 3D densities compared to 2D densities available in many GIS packages, and the inability to relate the resulting 3D densities to GIS data. Part of this issue is that 3D densities currently require volumetric rendering, which are not supported in GIS environments.
We developed an alternative approach which utilizes the computer-aided design program Blender with the Blender GIS add-on, and a custom 3D density add-on to conduct density analyses. The 3D density add-on works in conjunction with kernel based density estimation methods in the SciPy scientific computing package as well as GRASS GIS, an open source geospatial processing framework. The density add-on automates the import and density calculation process as well as the generation of 3D objects and components needed to visualise a semi-transparent data volume. The add-on allows filtering of data to view distributions of specific artefact types as well as settings to customise the grid resolution used for density statistics and coordinate-reference system used to register point positions.
An associated command line tool and Python library enables custom workflows for exporting the point densities into other formats such as image stacks and X3D compatible scenes for interactive web publication. The use of a library also provides researchers with the ability to reuse the developed components in novel ways that suit their particular analysis needs.
Application
We applied these methods to a case study from Waitapu, Ahuahu, New Zealand, an artefact rich site in a beach sand matrix. The number of artefacts from excavations at Waitapu provides the opportunity to consider how the 3D distribution of objects may be used to construct an alternative to grouping artefacts rather than conventional stratigraphic layers. These groupings have provided a new way of investigating artefact deposition and their relationship to other archaeological phenomena.
Through the use of this method we aim to remove the subjectivity around stratigraphic interpretation and construct a method that is reproducible in future research. The approach may be applied in any context with sufficient artefact recording, that is artefacts with x, y, and z data. The necessary scripts and add-ons will be published open source so that they may be further developed and can contribute to developments in 3D GIS.
![](https://www.eresearch.auckland.ac.nz/files/2020/01/image-for-Point-density-3840x2160-1.png)
Volume rendering of artefact distribution densities in Waitapu, Ahuahu, New Zealand. Image generated in Blender using GRASS GIS to compute spatial densities. Bright green regions indicate high artefact densities.
See more case study projects
![Our Voices: using innovative techniques to collect, analyse and amplify the lived experiences of young people in Aotearoa](https://www.eresearch.auckland.ac.nz/files/2024/11/Our-voice-banner-scaled-1-400x284.png)
Our Voices: using innovative techniques to collect, analyse and amplify the lived experiences of young people in Aotearoa
![Painting the brain: multiplexed tissue labelling of human brain tissue to facilitate discoveries in neuroanatomy](https://www.eresearch.auckland.ac.nz/files/2024/03/OFB-whole-section-images-32aca7cc70c4c921-400x284.png)
Painting the brain: multiplexed tissue labelling of human brain tissue to facilitate discoveries in neuroanatomy
![Detecting anomalous matches in professional sports: a novel approach using advanced anomaly detection techniques](https://www.eresearch.auckland.ac.nz/files/2024/03/Sport-fixing_NBA-646aaf7cc074ff84-400x284.jpg)
Detecting anomalous matches in professional sports: a novel approach using advanced anomaly detection techniques
![Benefits of linking routine medical records to the GUiNZ longitudinal birth cohort: Childhood injury predictors](https://www.eresearch.auckland.ac.nz/files/2024/03/Child-injury1-305c51cbbc7092bc-400x284.jpg)
Benefits of linking routine medical records to the GUiNZ longitudinal birth cohort: Childhood injury predictors
![Using a virtual machine-based machine learning algorithm to obtain comprehensive behavioural information in an in vivo Alzheimer’s disease model](https://www.eresearch.auckland.ac.nz/files/2024/03/DeepLab_one-side-692bb77ae6d4d336-400x284.jpg)
Using a virtual machine-based machine learning algorithm to obtain comprehensive behavioural information in an in vivo Alzheimer’s disease model
![Mapping livability: the “15-minute city” concept for car-dependent districts in Auckland, New Zealand](https://www.eresearch.auckland.ac.nz/files/2024/03/re_on_combine_micromobility-cdea073d4420e7c1-400x284.jpg)
Mapping livability: the “15-minute city” concept for car-dependent districts in Auckland, New Zealand
![Travelling Heads – Measuring Reproducibility and Repeatability of Magnetic Resonance Imaging in Dementia](https://www.eresearch.auckland.ac.nz/files/2023/03/Figure1_representative_image_final-400x284.jpg)
Travelling Heads – Measuring Reproducibility and Repeatability of Magnetic Resonance Imaging in Dementia
![Novel Subject-Specific Method of Visualising Group Differences from Multiple DTI Metrics without Averaging](https://www.eresearch.auckland.ac.nz/files/2023/03/Fig2_avg-value-visualsiation-400x284.png)
Novel Subject-Specific Method of Visualising Group Differences from Multiple DTI Metrics without Averaging
![Re-assess urban spaces under COVID-19 impact: sensing Auckland social ‘hotspots’ with mobile location data](https://www.eresearch.auckland.ac.nz/files/2021/11/K-road-400x284.jpg)
Re-assess urban spaces under COVID-19 impact: sensing Auckland social ‘hotspots’ with mobile location data
![Aotearoa New Zealand’s changing coastline – Resilience to Nature’s Challenges (National Science Challenge)](https://www.eresearch.auckland.ac.nz/files/2021/11/shoreline-change-rate_Nick-400x284.jpg)
Aotearoa New Zealand’s changing coastline – Resilience to Nature’s Challenges (National Science Challenge)
![Proteins under a computational microscope: designing in-silico strategies to understand and develop molecular functionalities in Life Sciences and Engineering](https://www.eresearch.auckland.ac.nz/files/2021/01/banner-for-shiga-toxin.jpg)
Proteins under a computational microscope: designing in-silico strategies to understand and develop molecular functionalities in Life Sciences and Engineering
![Coastal image classification and nalysis based on convolutional neural betworks and pattern recognition](https://www.eresearch.auckland.ac.nz/files/2020/10/2020_location-of-Tairua-beach-400x284.jpg)
Coastal image classification and nalysis based on convolutional neural betworks and pattern recognition
![Determinants of translation efficiency in the evolutionarily-divergent protist Trichomonas vaginalis](https://www.eresearch.auckland.ac.nz/files/2020/09/2020_featured-fig-4-400x284.jpg)
Determinants of translation efficiency in the evolutionarily-divergent protist Trichomonas vaginalis
![Measuring impact of entrepreneurship activities on students’ mindset, capabilities and entrepreneurial intentions](https://www.eresearch.auckland.ac.nz/files/2020/01/2019_cie-Dashboard-1-400x284.png)
Measuring impact of entrepreneurship activities on students’ mindset, capabilities and entrepreneurial intentions
![Using Zebra Finch data and deep learning classification to identify individual bird calls from audio recordings](https://www.eresearch.auckland.ac.nz/files/2020/01/Zebra-finch-group.jpg)
Using Zebra Finch data and deep learning classification to identify individual bird calls from audio recordings
![Automated measurement of intracranial cerebrospinal fluid volume and outcome after endovascular thrombectomy for ischemic stroke](https://www.eresearch.auckland.ac.nz/files/2020/01/segments-brain-400x284.jpg)
Automated measurement of intracranial cerebrospinal fluid volume and outcome after endovascular thrombectomy for ischemic stroke
![Using simple models to explore complex dynamics: A case study of macomona liliana (wedge-shell) and nutrient variations](https://www.eresearch.auckland.ac.nz/files/2020/01/parameterising-1.jpg)
Using simple models to explore complex dynamics: A case study of macomona liliana (wedge-shell) and nutrient variations
![Fully coupled thermo-hydro-mechanical modelling of permeability enhancement by the finite element method](https://www.eresearch.auckland.ac.nz/files/2018/04/2013-5-fem-velocity-15i9znp-400x284.png)
Fully coupled thermo-hydro-mechanical modelling of permeability enhancement by the finite element method
![Modelling dual reflux pressure swing adsorption (DR-PSA) units for gas separation in natural gas processing](https://www.eresearch.auckland.ac.nz/files/2018/04/2014-32-drpsa-exp-1qo75ef-400x284.jpg)
Modelling dual reflux pressure swing adsorption (DR-PSA) units for gas separation in natural gas processing
![Molecular phylogenetics uses genetic data to reconstruct the evolutionary history of individuals, populations or species](https://www.eresearch.auckland.ac.nz/files/2018/04/2014-28-phylogen-posterior-1j0qlkk-400x284.jpg)
Molecular phylogenetics uses genetic data to reconstruct the evolutionary history of individuals, populations or species
![Wandering around the molecular landscape: embracing virtual reality as a research showcasing outreach and teaching tool](https://www.eresearch.auckland.ac.nz/files/2018/04/2017-73-molecular-20voht3-400x284.png)
Wandering around the molecular landscape: embracing virtual reality as a research showcasing outreach and teaching tool
![Disposition of Microsoft HoloLenses for a Pop-Up Reality Shop to demonstrate the progress of a research project](https://www.eresearch.auckland.ac.nz/files/2018/04/2017-71-popup-banner-18wp6dq-400x284.jpg)