A computerized facial approximation method for Homo sapiens based on facial soft tissue thickness depths and geometric morphometrics.

Facial approximation (FA) provides a promising means of generating the possible facial appearance of a deceased person. It facilitates exploration of the evolutionary forces driving anatomical changes in ancestral humans and can capture public attention. Despite the recent progress made toward improving the performance of FA methods, a limited understanding of detailed quantitative craniofacial relationships between facial bone and soft tissue morphology may hinder their accuracy, and hence subjective experience and artistic interpretation are required. In this study, we explored craniofacial relationships among human populations based upon average facial soft tissue thickness depths (FSTDs) and covariations between hard and soft tissues of the nose and mouth using geometric morphometrics. Furthermore, we proposed a computerized method to assign the learned craniofacial relationships to generate a probable facial appearance of Homo sapiens, reducing human intervention. A smaller resemblance comparison (an average Procrustes distance was 0.0258 and an average Euclidean distance was 1.79 mm) between approximated and actual faces and a greater recognition rate (91.67%) tested by a face pool indicated that average dense FSTDs contributed to raising the accuracy of approximated faces. Results of partial least squares (PLS) analysis showed that nasal and oral hard tissues have an effect on their soft tissues separately. However, relatively weaker RV correlations (<0.4) and greater approximation errors suggested that we need to be cautious about the accuracy of the approximated nose and mouth soft tissue shapes from bony structures. Overall, the proposed method can facilitate investigations of craniofacial relationships and potentially improve the reliability of the approximated faces for use in numerous applications in forensic science, archaeology, and anthropology.

Skeleton-based cerebrovascular quantitative analysis.

BACKGROUND: Cerebrovascular disease is the most common cause of death worldwide, with millions of deaths annually. Interest is increasing toward understanding the geometric factors that influence cerebrovascular diseases, such as stroke. Cerebrovascular shape analyses are essential for the diagnosis and pathological identification of these conditions. The current study aimed to provide a stable and consistent methodology for quantitative Circle of Willis (CoW) analysis and to identify geometric changes in this structure. METHOD: An entire pipeline was designed with emphasis on automating each step. The stochastic segmentation was improved and volumetric data were obtained. The L1 medial axis method was applied to vessel volumetric data, which yielded a discrete skeleton dataset. A B-spline curve was used to fit the skeleton, and geometric values were proposed for a one-dimensional skeleton and radius. The calculations used to derive these values were illustrated in detail. RESULT: In one example(No. 47 in the open dataset) all values for different branches of CoW were calculated. The anterior communicating artery(ACo) was the shortest vessel, with a length of 2.6mm. The range of the curvature of all vessels was (0.3, 0.9) ± (0.1, 1.4). The range of the torsion was (-12.4,0.8) ± (0, 48.7). The mean radius value range was (3.1, 1.5) ± (0.1, 0.7) mm, and the mean angle value range was (2.2, 2.9) ± (0, 0.2) mm. In addition to the torsion variance values in a few vessels, the variance values of all vessel characteristics remained near 1. The distribution of the radii of symmetrical posterior cerebral artery(PCA) and angle values of the symmetrical posterior communicating arteries(PCo) demonstrated a certain correlation between the corresponding values of symmetrical vessels on the CoW. CONCLUSION: The data verified the stability of our methodology. Our method was appropriate for the analysis of large medical image datasets derived from the automated pipeline for populations. This method was applicable to other tubular organs, such as the large intestine and bile duct.

A Comparison of Semilandmarking Approaches in the Visualisation of Shape Differences.

In landmark-based analyses of size and shape variation and covariation among biological structures, regions lacking clearly identifiable homologous landmarks are commonly described by semilandmarks. Different algorithms may be used to apply semilandmarks, but little is known about the consequences of analytical results. Here, we assess how different approaches and semilandmarking densities affect the estimates and visualisations of mean and allometrically scaled surfaces. The performance of three landmark-driven semilandmarking approaches is assessed using two different surface mesh datasets with different degrees of variation and complexity: adult human head and ape cranial surfaces. Surfaces fitted to estimates of the mean and allometrically scaled landmark and semilandmark configurations arising from geometric morphometric analyses of these datasets are compared between semilandmarking approaches and different densities, as well as with those from warping to landmarks alone. We find that estimates of surface mesh shape (i.e., after re-semilandmarking and then re-warping) made with varying numbers of semilandmarks are generally consistent, while the warping of surfaces using landmarks alone yields surfaces that can be quite different to those based on semilandmarks, depending on landmark coverage and choice of template surface for warping. The extent to which these differences are important depends on the particular study context and aims.

A Comparison of Semilandmarking Approaches in the Analysis of Size and Shape.

Often, few landmarks can be reliably identified in analyses of form variation and covariation. Thus, 'semilandmarking' algorithms have increasingly been applied to surfaces and curves. However, the locations of semilandmarks depend on the investigator's choice of algorithm and their density. In consequence, to the extent that different semilandmarking approaches and densities result in different locations of semilandmarks, they can be expected to yield different results concerning patterns of variation and co-variation. The extent of such differences due to methodology is, as yet, unclear and often ignored. In this study, the performance of three landmark-driven semilandmarking approaches is assessed, using two different surface mesh datasets (ape crania and human heads) with different degrees of variation and complexity, by comparing the results of morphometric analyses. These approaches produce different semilandmark locations, which, in turn, lead to differences in statistical results, although the non-rigid semilandmarking approaches are consistent. Morphometric analyses using semilandmarks must be interpreted with due caution, recognising that error is inevitable and that results are approximations. Further work is needed to investigate the effects of using different landmark and semilandmark templates and to understand the limitations and advantages of different semilandmarking approaches.

A PCA-Based method for determining craniofacial relationship and sexual dimorphism of facial shapes.

Previous studies have used principal component analysis (PCA) to investigate the craniofacial relationship, as well as sex determination using facial factors. However, few studies have investigated the extent to which the choice of principal components (PCs) affects the analysis of craniofacial relationship and sexual dimorphism. In this paper, we propose a PCA-based method for visual and quantitative analysis, using 140 samples of 3D heads (70 male and 70 female), produced from computed tomography (CT) images. There are two parts to the method. First, skull and facial landmarks are manually marked to guide the model's registration so that dense corresponding vertices occupy the same relative position in every sample. Statistical shape spaces of the skull and face in dense corresponding vertices are constructed using PCA. Variations in these vertices, captured in every principal component (PC), are visualized to observe shape variability. The correlations of skull- and face-based PC scores are analysed, and linear regression is used to fit the craniofacial relationship. We compute the PC coefficients of a face based on this craniofacial relationship and the PC scores of a skull, and apply the coefficients to estimate a 3D face for the skull. To evaluate the accuracy of the computed craniofacial relationship, the mean and standard deviation of every vertex between the two models are computed, where these models are reconstructed using real PC scores and coefficients. Second, each PC in facial space is analysed for sex determination, for which support vector machines (SVMs) are used. We examined the correlation between PCs and sex, and explored the extent to which the choice of PCs affects the expression of sexual dimorphism. Our results suggest that skull- and face-based PCs can be used to describe the craniofacial relationship and that the accuracy of the method can be improved by using an increased number of face-based PCs. The results show that the accuracy of the sex classification is related to the choice of PCs. The highest sex classification rate is 91.43% using our method.

The production of digital and printed resources from multiple modalities using visualization and three-dimensional printing techniques.

PURPOSE: Virtual digital resources and printed models have become indispensable tools for medical training and surgical planning. Nevertheless, printed models of soft tissue organs are still challenging to reproduce. This study adopts open source packages and a low-cost desktop 3D printer to convert multiple modalities of medical images to digital resources (volume rendering images and digital models) and lifelike printed models, which are useful to enhance our understanding of the geometric structure and complex spatial nature of anatomical organs. MATERIALS AND METHODS: Neuroimaging technologies such as CT, CTA, MRI, and TOF-MRA collect serial medical images. The procedures for producing digital resources can be divided into volume rendering and medical image reconstruction. To verify the accuracy of reconstruction, this study presents qualitative and quantitative assessments. Subsequently, digital models are archived as stereolithography format files and imported to the bundled software of the 3D printer. The printed models are produced using polylactide filament materials. RESULTS: We have successfully converted multiple modalities of medical images to digital resources and printed models for both hard organs (cranial base and tooth) and soft tissue organs (brain, blood vessels of the brain, the heart chambers and vessel lumen, and pituitary tumor). Multiple digital resources and printed models were provided to illustrate the anatomical relationship between organs and complicated surrounding structures. Three-dimensional printing (3DP) is a powerful tool to produce lifelike and tangible models. CONCLUSIONS: We present an available and cost-effective method for producing both digital resources and printed models. The choice of modality in medical images and the processing approach is important when reproducing soft tissue organs models. The accuracy of the printed model is determined by the quality of organ models and 3DP. With the ongoing improvement of printing techniques and the variety of materials available, 3DP will become an indispensable tool in medical training and surgical planning.

The role of three-dimensional printed models of skull in anatomy education: a randomized controlled trail.

Three-dimensional (3D) printed models represent educational tools of high quality compared with traditional teaching aids. Colored skull models were produced by 3D printing technology. A randomized controlled trial (RCT) was conducted to compare the learning efficiency of 3D printed skulls with that of cadaveric skulls and atlas. Seventy-nine medical students, who never studied anatomy, were randomized into three groups by drawing lots, using 3D printed skulls, cadaveric skulls, and atlas, respectively, to study the anatomical structures in skull through an introductory lecture and small group discussions. All students completed identical tests, which composed of a theory test and a lab test, before and after a lecture. Pre-test scores showed no differences between the three groups. In post-test, the 3D group was better than the other two groups in total score (cadaver: 29.5 [IQR: 25-33], 3D: 31.5 [IQR: 29-36], atlas: 27.75 [IQR: 24.125-32]; p = 0.044) and scores of lab test (cadaver: 14 [IQR: 10.5-18], 3D: 16.5 [IQR: 14.375-21.625], atlas: 14.5 [IQR: 10-18.125]; p = 0.049). Scores involving theory test, however, showed no difference between the three groups. In this RCT, an inexpensive, precise and rapidly-produced skull model had advantages in assisting anatomy study, especially in structure recognition, compared with traditional education materials.

Densely calculated facial soft tissue thickness for craniofacial reconstruction in Chinese adults.

Craniofacial reconstruction (CFR) is used to recreate a likeness of original facial appearance for an unidentified skull; this technique has been applied in both forensics and archeology. Many CFR techniques rely on the average facial soft tissue thickness (FSTT) of anatomical landmarks, related to ethnicity, age, sex, body mass index (BMI), etc. Previous studies typically employed FSTT at sparsely distributed anatomical landmarks, where different landmark definitions may affect the contrasting results. In the present study, a total of 90,198 one-to-one correspondence skull vertices are established on 171 head CT-scans and the FSTT of each corresponding vertex is calculated (hereafter referred to as densely calculated FSTT) for statistical analysis and CFR. Basic descriptive statistics (i.e., mean and standard deviation) for densely calculated FSTT are reported separately according to sex and age. Results show that 76.12% of overall vertices indicate that the FSTT is greater in males than females, with the exception of vertices around the zygoma, zygomatic arch and mid-lateral orbit. These sex-related significant differences are found at 55.12% of all vertices and the statistically age-related significant differences are depicted between the three age groups at a majority of all vertices (73.31% for males and 63.43% for females). Five non-overlapping categories are given and the descriptive statistics (i.e., mean, standard deviation, local standard deviation and percentage) are reported. Multiple appearances are produced using the densely calculated FSTT of various age and sex groups, and a quantitative assessment is provided to examine how relevant the choice of FSTT is to increasing the accuracy of CFR. In conclusion, this study provides a new perspective in understanding the distribution of FSTT and the construction of a new densely calculated FSTT database for craniofacial reconstruction.

A regional method for craniofacial reconstruction based on coordinate adjustments and a new fusion strategy.

Craniofacial reconstruction recreates a facial outlook from the cranium based on the relationship between the face and the skull to assist identification. But craniofacial structures are very complex, and this relationship is not the same in different craniofacial regions. Several regional methods have recently been proposed, these methods segmented the face and skull into regions, and the relationship of each region is then learned independently, after that, facial regions for a given skull are estimated and finally glued together to generate a face. Most of these regional methods use vertex coordinates to represent the regions, and they define a uniform coordinate system for all of the regions. Consequently, the inconsistence in the positions of regions between different individuals is not eliminated before learning the relationships between the face and skull regions, and this reduces the accuracy of the craniofacial reconstruction. In order to solve this problem, an improved regional method is proposed in this paper involving two types of coordinate adjustments. One is the global coordinate adjustment performed on the skulls and faces with the purpose to eliminate the inconsistence of position and pose of the heads; the other is the local coordinate adjustment performed on the skull and face regions with the purpose to eliminate the inconsistence of position of these regions. After these two coordinate adjustments, partial least squares regression (PLSR) is used to estimate the relationship between the face region and the skull region. In order to obtain a more accurate reconstruction, a new fusion strategy is also proposed in the paper to maintain the reconstructed feature regions when gluing the facial regions together. This is based on the observation that the feature regions usually have less reconstruction errors compared to rest of the face. The results demonstrate that the coordinate adjustments and the new fusion strategy can significantly improve the craniofacial reconstructions.

A novel skull registration based on global and local deformations for craniofacial reconstruction.

Craniofacial reconstruction is important in forensic identification. It aims to estimate a facial appearance for human skeletal remains using the relationship between the soft tissue and the underlying bone structure. Various computerized methods have been developed in recent decades. An effective way is to deform a reference skull to the discovered skull, and then apply the same deformation to the skin associated with the reference skull to provide an approximate face for the discovered skull. For this method, the better the two skulls match each other, the more face-like the reconstructed skin surface will be. In this paper, we present a novel skull registration method that can match the two skulls closely, so as to improve the accuracy of the reconstruction. It combines both global and local deformations. A generic thin-plate spline (TPS)-based deformation, which is global, is applied first to roughly align the two skulls based on two groups of manually defined landmarks. Afterwards, the two skulls are largely matched, except some regions, on which some new landmarks are automatically marked. A compact support radial basis functions (CSRBF)-based deformation, which is local, will then be performed on these regions to adjust the initial alignment of the two skulls. Such adjustment can be repeatedly implemented until the two skulls have optimal alignment. In addition, all the skulls and face involved in the registration are represented by their single outer surfaces to facilitate the reconstruction procedure. The experiments demonstrate that our method can create a plausible face even when the reference skull is very different from the discovered skull. As a result, we can make full use of our database to provide multiple estimates for a principle components analysis (PCA) for the final reconstruction.