Virtual anthropology: the digital evolution in anthropological sciences

Comparability of skeletal fibulae surfaces generated by different source scanning (dual-energy CT scan vs. high resolution laser scanning) and 3D geometric morphometric validation

This work aims to test accuracy and comparability of 3D models of human skeletal fibulae generated by clinical CT and laser scanner virtual acquisitions. Mesh topology, segmentation and smoothing protocols were tested to assess variation among meshes generated with different scanning methods and procedures, and to evaluate meshes-interchangeability in 3D geometric morphometric analysis. A sample of 13 left human fibulae were scanned separately with Revolution Discovery CT dual energy (0.625 mm resolution) and ARTEC Space Spider 3D structured light laser scanner (0.1 mm resolution). Different segmentation methods, including half-maximum height (HMH) and MIA-clustering protocols, were compared to their high-resolution standard generated with laser-scanner by calculating topological surface deviations. Different smoothing algorithms were also evaluated, such as Laplacian and Taubin smoothing. A total of 142 semilandmarks were used to capture the shape of both proximal and distal fibular epiphyses. After Generalized Procrustes superimposition, the Procrustes coordinates of the proximal and distal fibular epiphyses were used separately to assess variation due to scanning methods and the operator error. Smoothing algorithms at low iteration do not provide significant variation among reconstructions, but segmentation protocol may influence final mesh quality (0.09-0.24 mm). Mean deviation among CT-generated meshes that were segmented with MIA-clustering protocol, and laser scanner-generated ones, is optimal (0.42 mm, ranging 0.35-0.56 mm). Principal component analysis reveals that homologous samples scanned with the two methods cluster together for both the proximal and distal fibular epiphyses. Similarly, Procrustes ANOVA reveals no shape differences between scanning methods and replicates, and only 1.38-1.43% of shape variation is due to scanning device. Topological similarities support the comparability of CT- and laser scanner-generated meshes and validate its simultaneous use in shape analysis with potential clinical relevance. We precautionarily suggest that dedicated trials should be performed in each study when merging different data sources prior to analyses.

Retrodeformation of the Steinheim Cranium: Insights into the Evolution of Neanderthals

A number of different approaches are currently available to digitally restore the symmetry of a specimen deformed by taphonomic processes. These tools include mirroring and retrodeformation to approximate the original shape of an object by symmetrisation. Retrodeformation has the potential to return a rather faithful representation of the original shape, but its power is limited by the availability of bilateral landmarks. A recent protocol proposed by Schlager and colleagues (2018) overcomes this issue by using bilateral landmarks and curves as well as semilandmarks. Here we applied this protocol to the Middle Pleistocene human cranium from Steinheim (Germany), the holotype of an abandoned species named Homo steinheimensis. The peculiar morphology of this fossil, associated with the taphonomic deformation of the entire cranium and the lack of a large portion of the right side of the face, has given rise to different hypotheses over its phylogenetic position. The reconstruction presented here sheds new light on the taphonomic origin of some features observed on this crucial specimen and results in a morphology consistent with its attribution to the Neanderthal lineage.

Sex estimation in a Turkish population using Purkait’s triangle: a virtual approach by 3-dimensional computed tomography (3D-CT)

Sex estimation is considered one of the first steps in the forensic identification process. Morphological and morphometrical differences between males and females have been used as means for morphoscopic and metric methods on both cranial and postcranial skeletal elements. When dry skeletal elements are not available, virtual data can be used as a substitute. The present research explores 3-dimensional (3D) scans from a Turkish population to test a sex estimation method developed by Purkait (2005).

Overall, 296 individuals were used in this study (158 males and 138 females). Purkait’s triangle parameters were measured on computed tomography (CT) scans obtained from both right and left femora of each patient at the Bakirkoy Dr. Sadi Konuk Training Research Hospital (Istanbul, Turkey). Intra- and inter-observer errors were assessed for all variables through technical error of measurements analysis. Bilateral asymmetry and sex differences were evaluated using parametric and non-parametric statistical approaches. Univariate and multivariate discriminant function analyses were then conducted.

Observer errors demonstrated an overall agreement within and between experts, as indicated by technical error of measurement (TEM) results. No bilateral asymmetries were reported, and all parameters demonstrated a statistically significant difference between males and females. Fourteen discriminant models were generated by applying single and combined parameters, producing a total correct sex classification ranging from 78.4% to 92.6%. In addition, over 67% of the total sample was accurately classified, with 95% or greater posterior probabilities.

Our study demonstrates the feasibility of 3D sex estimation using Purkait’s triangle on a Turkish population, with accuracy rates comparable to those reported in other populations. This is the first attempt to apply this method on virtual data and although further validation and standardisation are recommended for its application on dry bone, this research constitutes a significant contribution to the development of population-specific standards when only virtual data are available.

Key points

CT analysis using Purkait’s triangle is a suitable tool for assessment of sex in unidentified individuals.

The best overall estimation rate was achieved with the F11 model, with around 92% of accuracy.

The results suggested 78.4% to 92.6% correct sex identification rates.

More research is needed to expand the sample set and verify the results.

Impact of 3D surface scanning protocols on the Os coxae digital data: Implications for sex and age-at-death assessment

The 3D imaging technologies have become of paramount importance for example in disciplines such as forensic anthropology and bioarchaeology, where they are being used more and more frequently. There are several new possibilities that they offer; for instance, the easier and faster sharing of data among institutions, the possibility of permanent documentation, or new opportunities of data analysis. An important requirement, however, is whether the data obtained from different scanning devices are comparable and whether the possible varying outputs could affect further analyses, such as the estimation of the biological profile. Therefore, we aimed to investigate two important questions: (1) whether 3D models acquired by two different scanning technologies (structured light and laser) are comparable and (2) whether the scanning equipment has an effect on the anthropological analyses, such as age-at-death estimation and sex assessment. 3D models of ossa coxa (n = 29) were acquired by laser (NextEngine) and structured light (HP 3D Structured Light Scanner PRO 2) scanners. The resulting 3D models from both scanners were subjected to age-at-death analyses (via the quantitative method of Stoyanova et al., 2017) and sex analyses (via Diagnose Sexuelle Probabiliste 2 of Brůžek et al., 2017). Furthermore, high quality scans of a small sample (n = 5) of pubic symphyseal surfaces with the RedLux Profiler device were acquired as reference surfaces to which the outputs from both scanners were compared. Small deviations between surfaces were more evident in more rugged surfaces (in areas of depression and protrusion). Even though small differences from the reference surfaces were found, they did not have a significant effect on the age and sex estimates. It never resulted in the opposite sex assignment, and no significant differences were observed between age estimates (with the exception of those with the TPS/BE model).

Forensic age estimation based on the pigmentation in the costal cartilage from human mortal remains

Age estimation is considered a crucial and challenging issue in forensic casework. Costal cartilage appears a potential mortal remain in age-at-death estimation attributable to its correlative alteration in color based on pigment accumulation with the advancing age. In this study, samples from the second costal cartilage were collected in a Chinese Han population, and the cross sections were subsequently scanned and digitalized in a standard way. Color change was quantified using mean gray value (MGV), which was measured by Photoshop CS5. After the exclusion of samples with factors which could impair the quality of images and the accuracy of values, a high correlation was demonstrated between age and MGV in samples. A linear regression model (AGE = 173.425-0.755*aveMGV) was established for age prediction, with its performance evaluated using both samples from the training set and the blind test set, in which a mean absolute deviation of 4.42 years and 3.57 years was obtained, respectively. Altogether, MGV could be reckoned as a precise quantification of pigmentation in costal cartilage and an excellent indicator of age prediction in the age interval from 20 to 60 years. Moreover, our strategy appears more user-friendly and accurate, thus exceedingly practical for age estimation in forensic anthropology.

Evaluating causes of error in landmark-based data collection using scanners

In this study, we assess the precision, accuracy, and repeatability of craniodental landmarks (Types I, II, and III, plus curves of semilandmarks) on a single macaque cranium digitally reconstructed with three different surface scanners and a microCT scanner. Nine researchers with varying degrees of osteological and geometric morphometric knowledge landmarked ten iterations of each scan (40 total) to test the effects of scan quality, researcher experience, and landmark type on levels of intra- and interobserver error. Two researchers additionally landmarked ten specimens from seven different macaque species using the same landmark protocol to test the effects of the previously listed variables relative to species-level morphological differences (i.e., observer variance versus real biological variance). Error rates within and among researchers by scan type were calculated to determine whether or not data collected by different individuals or on different digitally rendered crania are consistent enough to be used in a single dataset. Results indicate that scan type does not impact rate of intra- or interobserver error. Interobserver error is far greater than intraobserver error among all individuals, and is similar in variance to that found among different macaque species. Additionally, experience with osteology and morphometrics both positively contribute to precision in multiple landmarking sessions, even where less experienced researchers have been trained in point acquisition. Individual training increases precision (although not necessarily accuracy), and is highly recommended in any situation where multiple researchers will be collecting data for a single project.

Digital reconstruction of the Ceprano calvarium (Italy), and implications for its interpretation

The Ceprano calvarium was discovered in fragments on March 1994 near the town of Ceprano in southern Latium (Italy), embedded in Middle Pleistocene layers. After reconstruction, its morphological features suggests that the specimen belongs to an archaic variant of H. heidelbergensis, representing a proxy for the last common ancestor of the diverging clades that respectively led to H. neanderthalensis and H. sapiens. Unfortunately, the calvarium was taphonomically damaged. The postero-lateral vault, in particular, appears deformed and this postmortem damage may have influenced previous interpretations. Specifically, there is a depression on the fragmented left parietal, while the right cranial wall is warped and angulated. This deformation affected the shape of the occipital squama, producing an inclination of the transverse occipital torus. In this paper, after X-ray microtomography (μCT) of both the calvarium and several additional fragments, we analyze consistency and pattern of the taphonomic deformation that affected the specimen, before the computer-assisted retrodeformation has been performed; this has also provided the opportunity to reappraise early attempts at restoration. As a result, we offer a revised interpretation for the Ceprano calvarium’s original shape, now free from the previous uncertainties, along with insight for its complex depositional and taphonomic history.

Virtual anthropology - a brief review of the literature and history of computed tomography

Computed tomography (CT) has influenced numerous fields since its inception in the 1970s. The field of palaeoanthropology significantly benefited from this efficient and non-invasive medium in terms of the conservation, reconstruction and analysis of fossil human remains. Over the past decade, there has been a steady increase in the number of forensic anthropological studies incorporating virtual osteological analyses. Because of the increasing importance of these modern cross-sectional imaging techniques and the requirement for standardized parameters in forensic science, we deemed it important to outline the history and development of CT applications in these related academic areas. The present paper outlines the history of “virtual anthropology” and osteological multi-detector CT in the context of palaeoanthropology and forensic anthropology.

Sex Estimation from Human Cranium: Forensic and Anthropological Interest of Maxillary Sinus Volumes

Sex estimation is a key objective of forensic science. We aimed to establish whether maxillary sinus volumes (MSV) could assist in estimating an individual's sex. One hundred and three CT scans were included. MSV were determined using three-dimensional reconstructions. Two observers performed three-dimensional MSV reconstructions using the same methods. Intra- and interobserver reproducibility were statistically compared using the intraclass correlation coefficient (ICC) (α = 5%). Both intra- and interobserver reproducibility were perfect regarding MSV; both ICCs were 100%. There were no significant differences between right and left MSV (p = 0.083). No correlation was found between age and MSV (p > 0.05). We demonstrated the existence of sexual dimorphism in MSV (p < 0.001) and showed that MSV measurements gave a 68% rate of correct allocations to sex group. MSV measurements could be useful to support sex estimation in forensic medicine.

Determination of sex: interest of frontal sinus 3D reconstructions

Frontal sinuses (FSs) have been studied in radiology, anthropology, and forensic anthropology. This study aimed to determine whether it was possible to predict the age and sex of an individual using FS volume. Sixty-nine anonymized CT scans were imported to MIMICS 10.01(®) software (Materialise N.V.), and each FS volume was calculated in mm(3) . There was an absence of significant difference between right and left FS volume (p = 0.173) and an absence of correlation between age and FS volume (Pearson's test; p = 0.705). Sexual dimorphism was significantly different (p = 0.001). However, the most discriminant datum for determining sex was found to be the total FS volume (sum of an individual's right and left FS volumes) with linear discriminant Fisher's function coefficients of -2.759 for the male group and -1.275 for the female group. With this model, 72.5% of our sample was correctly classified according to sex.

Morphometric analysis of variation in the ribs with age and sex

Rib cage morphology changes with age and sex are expected to affect thoracic injury mechanisms and tolerance, particularly for vulnerable populations such as pediatrics and the elderly. The size and shape variation of the external geometry of the ribs was characterized for males and females aged 0-100 years. Computed tomography (CT) scans from 339 subjects were analyzed to collect between 2700 and 10 400 homologous landmarks from each rib. Rib landmarks were analyzed using the geometric morphometric technique known as Procrustes superimposition. Age- and sex-specific functions of 3D rib morphology were produced representing the combined size and shape variation and the isolated shape variation. Statistically significant changes in the size and shape variation (P < 0.0001) and shape variation (P < 0.0053) of all 24 ribs were found to occur with age in males and females. Rib geometry, location, and orientation varied according to the rib level. From birth through adolescence, the rib cage experienced an increase in size, a decrease in thoracic kyphosis, and inferior rotation of the ribs relative to the spine within the sagittal plane. From young adulthood into elderly age, the rib cage experienced increased thoracic kyphosis and superior rotation of the ribs relative to the spine within the sagittal plane. The increased roundedness of the rib cage and horizontal angling of the ribs relative to the spine with age influences the biomechanical response of the thorax. With the plane of the rib oriented more horizontally, loading applied in the anterior-posterior direction will result in increased deformation within the plane of the rib and an increased risk for rib fractures. Thus, morphological changes may be a contributing factor to the increased incidence of rib fractures in the elderly. The morphological functions derived in this study capture substantially more information on thoracic skeleton morphology variation with age and sex than is currently available in the literature. The developed models of rib cage anatomy can be used to study age and sex variations in thoracic injury patterns due to motor vehicle crashes or falls, and clinically relevant changes due to chronic obstructive pulmonary disease or other diseases evidenced by structural and anatomic changes to the chest.

To 3D or not to 3D, that is the question: do 3D surface analyses improve the ecomorphological power of the distal femur in placental mammals?

Improvements in three-dimensional imaging technologies have renewed interest in the study of functional and ecological morphology. Quantitative approaches to shape analysis are used increasingly to study form-function relationships. These methods are computationally intensive, technically demanding, and time-consuming, which may limit sampling potential. There have been few side-by-side comparisons of the effectiveness of such approaches relative to more traditional analyses using linear measurements and ratios. Morphological variation in the distal femur of mammals has been shown to reflect differences in locomotor modes across clades. Thus I tested whether a geometric morphometric analysis of surface shape was superior to a multivariate analysis of ratios for describing ecomorphological patterns in distal femoral variation. A sample of 164 mammalian specimens from 44 genera was assembled. Each genus was assigned to one of six locomotor categories. The same hypotheses were tested using two methods. Six linear measurements of the distal femur were taken with calipers, from which four ratios were calculated. A 3D model was generated with a laser scanner, and analyzed using three dimensional geometric morphometrics. Locomotor category significantly predicted variation in distal femoral morphology in both analyses. Effect size was larger in the geometric morphometric analysis than in the analysis of ratios. Ordination reveals a similar pattern with arboreal and cursorial taxa as extremes on a continuum of morphologies in both analyses. Discriminant functions calculated from the geometric morphometric analysis were more accurate than those calculated from ratios. Both analysis of ratios and geometric morphometric surface analysis reveal similar, biologically meaningful relationships between distal femoral shape and locomotor mode. The functional signal from the morphology is slightly higher in the geometric morphometric analysis. The practical costs of conducting these sorts of analyses should be weighed against potentially slight increases in power when designing protocols for ecomorphological studies.

Post-mortem computed tomography and 3D imaging: anthropological applications for juvenile remains

Anthropological examination of defleshed bones is routinely used in medico-legal investigations to establish an individual's biological profile. However, when dealing with the recently deceased, the removal of soft tissue from bone can be an extremely time consuming procedure that requires the presence of a trained anthropologist. In addition, due to its invasive nature, in some disaster victim identification scenarios the maceration of bones is discouraged by religious practices and beliefs, or even prohibited by national laws and regulations. Currently, three different radiological techniques may be used in the investigative process; plain X-ray, dental X-ray and fluoroscopy. However, recent advances in multi-detector computed tomography (MDCT) mean that it is now possible to acquire morphological skeletal information from high resolution images, reducing the necessity for invasive procedures. This review paper considers the possible applications of a virtual anthropological examination by reviewing the main juvenile age determination methods used by anthropologists at present and their possible adaption to MDCT.

GIS and paleoanthropology: incorporating new approaches from the geospatial sciences in the analysis of primate and human evolution

The incorporation of research tools and analytical approaches from the geospatial sciences is a welcome trend for the study of primate and human evolution. The use of remote sensing (RS) imagery and geographic information systems (GIS) allows vertebrate paleontologists, paleoanthropologists, and functional morphologists to study fossil localities, landscapes, and individual specimens in new and innovative ways that recognize and analyze the spatial nature of much paleoanthropological data. Whether one is interested in locating and mapping fossiliferous rock units in the field, creating a searchable and georeferenced database to catalog fossil localities and specimens, or studying the functional morphology of fossil teeth, bones, or artifacts, the new geospatial sciences provide an essential element in modern paleoanthropological inquiry. In this article we review recent successful applications of RS and GIS within paleoanthropology and related fields and argue for the importance of these methods for the study of human evolution in the twenty first century. We argue that the time has come for inclusion of geospatial specialists in all interdisciplinary field research in paleoanthropology, and suggest some promising areas of development and application of the methods of geospatial science to the science of human evolution.

Identification of group affinity from cross-sectional contours of the human midfacial skeleton using digital morphometrics and 3D laser scanning technology

Identifying group affinity from human crania is a long-standing problem in forensic and physical anthropology. Many craniofacial differences used in forensic skeletal identification are difficult to quantify, although certain measurements of the midfacial skeleton have shown high predictive value for group classifications. This study presents a new method for analyzing midfacial shape variation between different geographic groups. Three-dimensional laser scan models of 90 crania from three populations were used to obtain cross-sectional midfacial contours defined by three standard craniometric landmarks. Elliptic Fourier transforms of the contours were used to extract Fourier coefficients for statistical analysis. After cross-validation, discriminant functions based on the Fourier coefficients provided an average of 86% correct classifications for crania from the three groups. The high rate of accuracy of this method indicates its usefulness for identifying group affinities among human skeletal remains and demonstrates the advantages of digital 3D model-based analysis in forensic research.

Geometric morphometric methods for bone reconstruction: the mandibular condylar process of Pico della Mirandola

The issue of reconstructing lost or deformed bone presents an equal challenge in the fields of paleoanthropology, bioarchaeology, forensics, and medicine. Particularly, within the disciplines of orthodontics and surgery, the main goal of reconstruction is to restore or create ex novo the proper form and function. The reconstruction of the mandibular condyle requires restoration of articulation, occlusion, and mastication from the functional side as well as the correct shape of the mandible from the esthetic point of view. Meeting all these demands is still problematic for surgeons. It is unfortunate that the collaboration between anthropologists and medical professionals is still limited. Nowadays, geometric morphometric methods (GMM) are routinely applied in shape analysis and increasingly in the reconstruction of missing data in skeletal material in paleoanthropology. Together with methods for three-dimensional (3D) digital model construction and reverse engineering, these methods could prove to be useful in surgical fields for virtual planning of operations and the production of customized biocompatible scaffolds. In this contribution, we have reconstructed the missing left condylar process of the mandible belonging to a famous Italian humanist of the 15th century, Pico della Mirandola (1463-1494) by means of 3D digital models and GMM, having first compared two methods (a simple reflection of the opposite side and the mathematical-statistical GMM approach) in a complete human mandible on which loss of the left condyle was virtually simulated. Finally, stereolithographic models of Pico's skull were prototyped providing the physical assembly of the bony skull structures with a high fitting accuracy.

The feeding biomechanics and dietary ecology of Australopithecus africanus

The African Plio-Pleistocene hominins known as australopiths evolved a distinctive craniofacial morphology that traditionally has been viewed as a dietary adaptation for feeding on either small, hard objects or on large volumes of food. A historically influential interpretation of this morphology hypothesizes that loads applied to the premolars during feeding had a profound influence on the evolution of australopith craniofacial form. Here, we test this hypothesis using finite element analysis in conjunction with comparative, imaging, and experimental methods. We find that the facial skeleton of the Australopithecus type species, A. africanus, is well suited to withstand premolar loads. However, we suggest that the mastication of either small objects or large volumes of food is unlikely to fully explain the evolution of facial form in this species. Rather, key aspects of australopith craniofacial morphology are more likely to be related to the ingestion and initial preparation of large, mechanically protected food objects like large nuts and seeds. These foods may have broadened the diet of these hominins, possibly by being critical resources that australopiths relied on during periods when their preferred dietary items were in short supply. Our analysis reconciles apparent discrepancies between dietary reconstructions based on biomechanics, tooth morphology, and dental microwear.

Age and gender based biomechanical shape and size analysis of the pediatric brain

Injuries caused by motor vehicle crashes (MVCs) are the leading cause of head injury and death for children in the United States. This study aims to describe the shape and size (morphologic) changes of the cerebrum, cerebellum, brainstem, and ventricles of the pediatric occupant to better predict injury and assess how these changes affect finite element model (FEM) response. To quantify morphologic differences in the brain, a Generalized Procrustes Analysis (GPA) with a sliding landmark method was conducted to isolate morphologic changes using magnetic resonance images of 63 normal subjects. This type of geometric morphometric analysis was selected for its ability to identify homologous landmarks on structures with few true landmarks and isolate the shape and size of the individuals studied. From the resulting landmark coordinates, the shape and size changes were regressed against age to develop a model describing morphologic changes in the pediatric brain as a function of age. The most statistically significant shape change was in the cerebrum with p-values of 0.00346 for males and 0.00829 for females. The age-based model explains over 80% of the variation in size in the cerebrum. Using size and shape models, affine transformations were applied to the SIMon FEM to determine differences in response given differences in size and size plus shape. The geometric centroid of the elements exceeding 15% strain was calculated and compared to the geometric centroid of the entire structure. Given the same Haversine pulse, the centroid location, a metric for the spatial distribution of the elements exceeding an injury threshold, varied based on which transformation was applied to the model. To assess the overall response of the model, three injury metrics were examined to determine the magnitude of the metrics each element sustained and the overall volume of elements that experienced that value. These results suggested that the overall response of the model was driven by the variation in size, with little variation due to changes in shape. This study demonstrates a new methodology to quantify the shape and size variation of the brain from infancy to adulthood. The use of the changes in shape and size when applied to a FEM suggests that there are differences in the spatial distribution of the elements that exceed a specific threshold based on shape but the overall volume of elements experiencing the specified magnitude was more dependent on the changes in the size of the model with little change due to shape.

Virtual palpation of skeletal landmarks with multimodal display interfaces

The 3D location of skeletal landmarks on CT datasets is an important procedure, used in many research and clinical contexts. The standard procedure involves the segmentation of the CT images, the creation of a 3D surface bone model, and the location of the landmarks on this surface. However, the segmentation is time-consuming and requires skilled operators and sophisticated software. The aim of the present study is to evaluate the efficacy of a multimodal display interface to direct volumetric interactive visualization in performing a virtual palpation task. An expert operator used the CT dataset of a patient's thigh region to locate 14 femoral skeletal landmarks. This operation was repeatedly performed using different CT data representation; the accuracy and repeatability were compared to those achievable with the conventional procedure based on the segmented 3D surface. When a multimodal display interface (formed by an orthogonal slice, RXCT and interactive isosurface views) was used to perform the virtual palpation directly on the CT data, the average coordinates of the landmarks did not differ significantly from those located on the 3D surface, and the measurement repeatability was actually better with the multimodal display of the volumetric data than with the 3D surface. Thus, we can conclude that skeletal virtual palpation can be performed directly on the CT dataset, as far as the virtual palpation is performed with a multimodal display interface.

Facial surface analysis by 3D laser scanning and geometric morphometrics in relation to sexual dimorphism in cerebral--craniofacial morphogenesis and cognitive function

Over early fetal life the anterior brain, neuroepithelium, neural crest and facial ectoderm constitute a unitary, three-dimensional (3D) developmental process. This intimate embryological relationship between the face and brain means that facial dysmorphogenesis can serve as an accessible and informative index of brain dysmorphogenesis in neurological and psychiatric disorders of early developmental origin. There are three principal challenges in seeking to increase understanding of disorders of early brain dysmorphogenesis through craniofacial dysmorphogenesis: (i) the first, technical, challenge has been to digitize the facial surface in its inherent three-dimensionality; (ii) the second, analytical, challenge has been to develop methodologies for extracting biologically meaningful shape covariance from digitized samples, making statistical comparisons between groups and visualizing in 3D the resultant statistical models on a 'whole face' basis; (iii) the third, biological, challenge is to demonstrate a relationship between facial morphogenesis and brain morphogenesis not only in anatomical-embryological terms but also at the level of brain function. Here we consider each of these challenges in turn and then illustrate the issues by way of our own findings. These use human sexual dimorphism as an exemplar for 3D laser surface scanning of facial shape, analysis using geometric morphometrics and exploration of cognitive correlates of variation in shape of the 'whole face', in the context of studies relating to the early developmental origins of schizophrenia.