Biomedical Signal Processing and Modeling Complexity of Living Systems 2014
View this Special IssueResearch Article  Open Access
Craniofacial Reconstruction Evaluation by Geodesic Network
Abstract
Craniofacial reconstruction is to estimate an individual’s face model from its skull. It has a widespread application in forensic medicine, archeology, medical cosmetic surgery, and so forth. However, little attention is paid to the evaluation of craniofacial reconstruction. This paper proposes an objective method to evaluate globally and locally the reconstructed craniofacial faces based on the geodesic network. Firstly, the geodesic networks of the reconstructed craniofacial face and the original face are built, respectively, by geodesics and isogeodesics, whose intersections are network vertices. Then, the absolute value of the correlation coefficient of the features of all corresponding geodesic network vertices between two models is taken as the holistic similarity, where the weighted average of the shape index values in a neighborhood is defined as the feature of each network vertex. Moreover, the geodesic network vertices of each model are divided into six subareas, that is, forehead, eyes, nose, mouth, cheeks, and chin, and the local similarity is measured for each subarea. Experiments using 100 pairs of reconstructed craniofacial faces and their corresponding original faces show that the evaluation by our method is roughly consistent with the subjective evaluation derived from thirtyfive persons in five groups.
1. Introduction
Craniofacial reconstruction [1] aims to estimate an individual’s face appearance from its skull using the relationship between soft tissues and the underlying bone structure. It has a widespread application in many areas, such as forensic medicine, archeology, medical cosmetic surgery, and public safety. With the development of 3D digitalization technology, the research on computer aided craniofacial reconstruction has widely received attention. The evaluation of the craniofacial reconstruction has an important significance in improving the craniofacial reconstruction methods. However, nearly all researches on craniofacial reconstruction focus on the reconstruction method itself, and little attention is paid to the evaluation method of the reconstruction results.
The craniofacial face is one of the most complex geometric objects in the natural world. How to evaluate the results of the craniofacial reconstruction is still a challenging problem. Existing evaluation methods of the craniofacial reconstruction can be divided into three types: subjective qualitative evaluation method, objective quantitative evaluation method, and the combination method of subjective and objective evaluation. Subjective evaluation methods evaluate the craniofacial reconstruction results subjectively by designing different evaluation strategies. Quatrehomme et al. [2] invited 25 subjects to evaluate and got excellent or good to middle resemblance in 9 out of 25 cases. Snow et al. [3] invited more than 200 respondents to compare a manually recovered craniofacial model with 7 photos; 68% of men and 26% of women gave the correct results. Stephan and Henneberg [4] invited 37 respondents to verify 16 craniofacial reconstruction results; their experimental results also indicated that the recognition rate of the men was higher than that of the women. Helmer et al. [5] invited 24 respondents to compare 24 reconstructed craniofacial faces with their real photos; their results showed that 38% of the reconstruction results were very similar, 17% were similar, 42% were mildly similar, and only one result was considered to be unrelated to the original photo. Although subjective evaluation method is consistent with human cognitive theory, the evaluation process requires a lot of manpower and time, and the accuracy of the evaluation results is influenced by human subjective factors.
Some scholars made a preliminary exploration on evaluating craniofacial reconstruction results by objective method. Feng et al. [6] used a relative anglecontext distribution (RACD) to compare the two craniofacial faces. They defined the probability density function of the relative anglecontext distribution by counting the number of the relative angles in different intervals. Considering the calculation instability and high time complexity of RACD, Zhu et al. [7] extended the RACD to bendingrelative anglecontext distribution (BRACD) algorithm by measuring the bending of a reference hemisphere to a craniofacial model. Duan et al. [8, 9] analyzed the correlation between the skull and face shape and measured the craniofacial reconstruction error by the distance of corresponding points between the reconstructed craniofacial face and the original face.
Several researchers combined subjective and objective evaluations. For example, VaneZis [10] invited 20 assessors to choose the three best matches from 10 reconstructed craniofacial faces of one skull with the original face. They also computed the correlation between the subjective results with the objective evaluation results by mathematical shape analysis assessment using Procrustes Analysis. The results are not statistically significant but indicate that the objective method does seem to capture some perceptual similarity in human observers. Moorthy et al. [11] also proposed a combination evaluation method. In subjective aspect, they conducted a subjective study, where a set of human subjects (12 subjects on 180 3D faces) rated the similarity of pairs of faces (a total of 5490 pairs of similarity scores). In objective aspect, they extracted Gabor features from automatically detected feature points on the range and texture images from 3D faces. Finally, they demonstrated that these features correlated well with human judgment of similarity.
In this paper, we propose a new global and local evaluation method of craniofacial reconstruction based on the geodesic network. We define the weighted average of the shape index value in a neighborhood as the feature of one vertex. The absolute value of the correlation coefficient of the feature of all corresponding geodesic network vertices between two models is taken as the similarity. It lays a foundation for qualitative and quantitative analysis of craniofacial reconstruction results and provides guidance for improvement of the craniofacial reconstruction methods.
The rest of this paper is organized as follows. The materials and methods are presented in Section 2. Section 3 presents experimental results. Some conclusions are provided in Section 4.
2. Materials and Methods
2.1. Materials
This study was approved by the Institutional Review Board (IRB) of Image Center for Brain Research, National Key Laboratory of Cognitive Neuroscience and Learning, Beijing Normal University. The study was carried out on a database of 208 whole head CT scans on voluntary persons that mostly come from Han ethnic group in North of China, aged from 19 to 75 years. The CT scans were obtained by a clinical multislice CT scanner system (Siemens Sensation16) in Xianyang Hospital located in western China. As Figure 1 shows, we firstly extract skull and face borders from the original CT slice images and then reconstruct the triangle mesh models of the 3D skull and skin surfaces by a marching cubes algorithm [12]. More details on the procedure of the data processing can be found in [13]. The obtained threedimensional craniofacial mesh models often contain defects such as holes, gaps, degeneracies, or nonmanifold configurations. We need to fill holes and gaps and remove the scattered points to make the 3D face model become a full and wellstructured manifold. To eliminate the effects of data acquisition, posture, and scale, all 3D craniofacial data are transformed into a unified Frankfurt coordinate system [14]. Finally, as Figure 2 shows, we select a craniofacial data as a reference template and cut away the back part of the reference craniofacial model considering that there are too many vertices in the whole head and the face features mainly concentrate on the front part of the head. All of the craniofacial models are registered with the reference model automatically by the nonrigid data registration method in [14].
(a)
(b)
(c)
(d)
(e)
2.2. Overview of the Method
The craniofacial reconstruction evaluation is to assess the similarity between the reconstructed craniofacial face and the original face. We propose an objective method of the craniofacial reconstruction evaluation based on the shape index of the geodesic network vertices. The method evaluates the similarity both globally and locally. The geodesic network of the reconstructed craniofacial face and the original face are built, respectively, by geodesics and isogeodesics, whose intersections are geodesic network vertices. The weighted average of the shape index values within a neighborhood of each geodesic network vertex is computed for each craniofacial face model. The weighted average is considered as the feature of the vertex. The absolute value of the correlation coefficient of the features of all corresponding geodesic network vertices is taken as the holistic similarity between two models. To evaluate the reconstruction locally, these geodesic network vertices of each model are divided into six subareas, forehead, eyes, nose, mouth, cheeks, and chin. The absolute value of the correlation coefficient of the features of corresponding network vertices in each subarea between two craniofacial models is taken as the local similarity.
2.3. Geodesic Network Construction
Geodesic is the curve with geodesic curvatures equal to zero. It is the shortest path between two points on a surface. Due to its intrinsic invariance, geodesic can be applied to face recognition that is insensitive to expression. According to geodesic, a geodesic polar coordinate system can be constructed by exponential mapping [15, 16]. Inspired by this idea, we construct the geodesic network according to geodesics and isogeodesics, respectively, on two craniofacial facial models and take the intersection points of geodesics and isogeodesics, that is, geodesic network vertices, as the corresponding points between the two craniofacial face models. Thus, we can compare the similarity of the two models through the features of the corresponding geodesic network vertices.
The detailed steps for constructing the geodesic network are described as follows. Firstly, we find the nose tip of the craniofacial face model and make it as the source point of geodesics. Because the craniofacial models are frontal, the nose tip point is the highest point of the whole model, that is, the point which has the biggest value (or value) in Frankfurt coordinate. Secondly, isogeodesics are extracted according to the geodesic distance to the source point. The extracted isogeodesics are welldistributed with the same interval between two of them. The interval is one out of of the shortest geodesic distance from the nose tip point to the boundary point set . So the isogeodesics set can be obtained by the following formula: where represents the isogeodesic that the geodesic distance from the source point to the point on the curve is equal to . Geodesic distance can be solved by the existing geodesic algorithms, such as MMP [17], ICH [18], PCH [19], and SVG [20]. Here, we adopt the classical MMP algorithm, which has been realized and described in detail in [21]. Thirdly, geodesics are computed from the tip nose according to an equal angular interval. The equal angle can be found as follows: the direction of the nose tip and the middle point of eyebrows is selected as the initial direction; the outermost isogeodesic is projected into the tangent plane of the nose tip, and the cross points of the equal angular line from the nose tip with the projection of the outermost isogeodesic on the tangent plane are computed; project the cross point into original outermost isogeodesic, and these equal division points are taken as the target points of geodesics; geodesic from the nose tip to the target points can be computed as follows: Finally, the intersection points of the geodesics and isogeodesics are calculated, and all intersection points constitute the intersection point set , that is, geodesic network point set:
Figure 3 shows the geodesic networks of the reconstructed and the original craniofacial face models. Let and be the geodesic network point sets of the craniofacial models and , respectively. The geodesic network point of and of are the corresponding points of the two craniofacial face models since they have the same initial directions and the same source (the nose tip).
(a) The reconstructed craniofacial face model
(b) The original craniofacial face model
2.4. Extracting Shape Index Features
Due to the intrinsic invariance of geodesic, the features of the geodesic network vertices can be used to evaluate the results of craniofacial reconstruction. Shape index is the feature generated by the principal curvatures and . It can capture the “local” shape of a surface. Thus, we evaluate the craniofacial reconstruction based on the correlation coefficient of the shape index value of the corresponding geodesic network vertices between the original and reconstructed craniofacial faces. Shape index is defined by the curvednessorientationshape map on sphere (COSMOS) representation [22] and can quantitatively measure the shape of a surface at a point . Every distinct surface shape corresponds to a unique value of . At a point , the shape index is defined as
For the corresponding network vertexes and at the original and the reconstructed craniofacial models, the shape index values are and , respectively. To avoid the interference of the noise and obtain the robust shape index value, we define the weighted average of the shape index value in a neighborhood as the feature of one vertex. Owing to the important role of the geodesic network vertex in evaluating the two craniofacial models, the weight of the geodesic network vertex itself should be the highest, and the weights of the 1ring neighborhood of the vertex should be higher than those of the 2ring neighborhood. In this paper, we choose the weights of 100, 10, and 1 for the network vertex, 1ring neighborhood, and 2ring neighborhood, respectively. The features of the corresponding network vertex and of the two craniofacial models are and , respectively. For the two whole models, the features of all geodesic vertices constitute, respectively, the vectors and where is equal to , that is, the product of the number of isogeodesics and the number of geodesics.
2.5. Evaluating the Craniofacial Reconstruction Globally and Locally
Generally, the correlation coefficient is used to measure how two sets of variables are linearly related. Here, we use the correlation coefficient to measure the similarity of two craniofacial models. For evaluating the craniofacial reconstruction result globally, we compute the correlation coefficient of the features of all corresponding geodesic network vertices between the original and the reconstructed craniofacial models by the following formula: The value of correlation coefficient ranges in the interval . We take the absolute value as the similarity measure. The smaller the value is, the weaker the similarity is.
The global evaluation reflects only the whole similarity between the reconstructed and the original craniofacial models and cannot discern whether the region is well reconstructed or not. In order to evaluate the reconstructed craniofacial model locally, we divide the geodesic network vertices into six parts, that is, forehead, eyes, nose, mouth, cheeks, and chin (as shown in Figure 4). The features of the geodesic network vertices in each subarea are used to compute the similarity measure of the subarea by (5). According to the local evaluation, we can make a statistical analysis of the local similarity measure and give some feedback to promote the craniofacial reconstruction methods.
(a) The reconstructed craniofacial face model
(b) The original craniofacial face model
3. Experiments and Discussion
To acquire the experimental data, we use the partial least squares regression (PLSR) method [8] to reconstruct the craniofacial models, where 108 pairs of skulls and face skins among the 208 CT scans are used as the training data, and the other 100 skulls are used as the test data for the craniofacial reconstruction. Thus, we obtain 100 pairs of the reconstructed face models and the corresponding original craniofacial models. To compare with the subjective evaluation, we first introduce the subjective evaluation procedure we designed. Then, the reconstruction results are evaluated globally and locally by the proposed objective method.
3.1. Subjective Evaluation
In order to evaluate the proposed objective method, we invited 35 subjects to evaluate the 100 reconstructed craniofacial face models. 100 pairs of craniofacial face models were divided into five groups and each group had 20 pairs. 35 subjects were divided into five groups, and each group had seven subjects to evaluate twenty pairs of reconstructed craniofacial faces and corresponding original craniofacial faces. Each subject was asked to observe every pair of faces on the screen and choose the overall similarity degree from the following five degrees as Figure 5 shows: sufficiently (above 90%), highly (70%~90%), somewhat (50%~70%), slightly (30%~50%), and lowly (0%~30%). They were also asked to select the most similar area and the least similar area from the following six subareas: nose, eyes, mouth, forehead, cheeks, and chin. To avoid visual fatigue, each subject was only in charge of evaluating twenty pairs of craniofacial faces. According to the subjective assessing from the five groups, the similarities of 100 pairs of craniofacial faces were gained, and each pair had seven similarity degrees by seven different subjects. We computed the mean maximum similarity score according to the upper limits of the seven similarity degrees and the mean minimum similarity score by the lower limits. Thus, we obtained a similarity interval for each pair of craniofacial faces by the subjective evaluation. In the following, if the similarity obtained by the objective method is in the interval for one pair of craniofacial faces, we think the evaluation is consistent with the subjective evaluation.
3.2. Global Evaluation
The global evaluation is to compare the features of all geodesic network vertices to get the objective similarity score between the reconstructed craniofacial face and the corresponding original face. Table 1 shows the similarity scores by the objective and subjective assessments. From Table 1, we can see that the similarity scores by the objective assessment are within the similarity interval of the subjective evaluation. The assessment results show that the objective assessment is consistent with the human subjective evaluation.

3.3. Local Evaluation
Local evaluation is the evaluation of the similarity in six subareas: forehead, nose, eyes, mouth, cheeks, and chin between the reconstructed craniofacial face and the corresponding original craniofacial face. We find the most and the least similar areas and compare them with the subjective assessments. We take three pairs of craniofacial faces in Table 1; the local similarity scores are computed by the features of the geodesic network vertices in each subarea. Table 2 shows the evaluation results by the objective method. We can see that the nose area is the most similar area and the eyes area is the least similar area for the case 0401; the mouth area is the most similar area and the forehead area is the least similar area for the case 0012354; the cheeks area is the most similar area and the mouth area is the least similar area for the case 3718. These results are consistent with the subjective evaluation.

We evaluate all of the 100 reconstructed craniofacial faces locally and compare the local evaluation of six subareas with the global evaluation. Figure 6 shows the comparisons. From Figure 6, we can see that the nose similarity is roughly consistent with the global similarity and the similarity of eyes is almost unrelated with the global similarity. We also compute the absolute value of correlation coefficients between the local similarities of the six subareas with the global similarities using the l00 cases. Table 3 shows the absolute value of correlation coefficients. From Table 3, we can also see that the similarity of the nose area is highly interrelated with the global similarity and the eyes area is the region of the lowest correlation with the global similarity. Compared with other subareas, the nose area has more of the geodesic network vertices and the curvature of the nose area changes more significantly; so the similarity of nose is closer to the global similarity. The eyes area has fewer geodesic network vertices; thus the similarity of eyes has a little correlation with the global similarity.

(a) Nose similarity and global similarity
(b) Eyes similarity and global similarity
(c) Chin similarity and global similarity
(d) Forehead similarity and global similarity
(e) Cheeks similarity and global similarity
(f) Mouth similarity and global similarity
We compute the mean similarity score of each subarea of the 100 pairs of faces. The mean similarity scores are shown in Table 4. We can see that the maximum is in the chin area and the minimum is in the eyes area. This indicates that the eyes area is not well reconstructed. These objective evaluation results are consistent with the subjective evaluation results.

4. Conclusions
Craniofacial reconstruction has a widespread application in forensic medicine, archeology, medical cosmetic surgery, and so forth. However, nearly all researches on craniofacial reconstruction focus on the reconstruction method itself, and little attention is paid to the evaluation methods of the reconstruction. This paper proposes an objective method to evaluate globally and locally the reconstructed craniofacial faces based on the shape index of geodesic network vertices. The geodesic networks of the craniofacial faces are built by geodesics and isogeodesics. For each geodesic network vertex, we define the weighted average of the shape index value in a neighborhood as the feature of the network vertex. The absolute value of the correlation coefficient of the feature of all corresponding geodesic network vertices between two models is taken as the similarity. We used 100 pairs of reconstructed craniofacial faces and their corresponding original faces to evaluate our method. To compare with the subjective evaluation, we also invited 35 subjects to evaluate visually the reconstructed craniofacial faces. Experimental results show that the evaluation by our method is roughly consistent with the subjective evaluation. By evaluating the craniofacial reconstruction effects both globally and locally, we can provide guidance for improvement of the craniofacial reconstruction methods. In addition, since small face expression can be regarded as the isometric deformation, under which the geodesic distance is invariant, the proposed method is also fit for the craniofacial faces of small expression variation.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors gratefully appreciated the anonymous reviewers for all of their helpful comments. They also thank Surazhsky et al. and Rusinkiewicz for their public part code of geodesic or curvature. This work is partially supported by the National Natural Science Foundation of China (nos. 61170170 and 61272363) and the Program for New Century Excellent Talents in University (NCET130051).
References
 P. Claes, D. Vandermeulen, S. De Greef, G. Willems, J. G. Clement, and P. Suetens, “Computerized craniofacial reconstruction: Conceptual framework and review,” Forensic Science International, vol. 201, no. 1–3, pp. 138–145, 2010. View at: Publisher Site  Google Scholar
 G. Quatrehomme, T. Balaguer, P. Staccini, and V. AlunniPerret, “Assessment of the accuracy of threedimensional manual craniofacial reconstruction: a series of 25 controlled cases,” International Journal of Legal Medicine, vol. 121, no. 6, pp. 469–475, 2007. View at: Publisher Site  Google Scholar
 C. C. Snow, B. P. Gatliff, and K. R. McWilliams, “Reconstruction of facial features from the skull: an evaluation of its usefulness in forensic anthropology,” American Journal of Physical Anthropology, vol. 33, no. 2, pp. 221–228, 1970. View at: Google Scholar
 C. N. Stephan and M. Henneberg, “Building faces from dry skulls: are they recognized above chance rates?” Journal of Forensic Sciences, vol. 46, no. 3, pp. 432–440, 2001. View at: Google Scholar
 R. P. Helmer, S. Rohricht, D. Petersen et al., “Assessment of the reliability of facial reconstruction,” in Forensic Analysis of the Skull, WileyLiss, New York, NY, USA, 1993. View at: Google Scholar
 J. Feng, H. H. S. Ip, L. Y. Lai, and A. Linney, “Robust point correspondence matching and similarity measuring for 3D models by relative anglecontext distributions,” Image and Vision Computing, vol. 26, no. 6, pp. 761–775, 2008. View at: Publisher Site  Google Scholar
 X.Y. Zhu, G.H. Geng, and C. Wen, “Craniofacial similarity measuring based on BRACD,” in Proceedings of the 4th International Conference on Biomedical Engineering and Informatics (BMEI ’11), vol. 2, pp. 942–945, Shanghai, China, October 2011. View at: Publisher Site  Google Scholar
 F. Duan, S. Yang, D. Huang et al., “Craniofacial reconstruction based on multilinear subspace analysis,” Multimedia Tools and Applications, 2013. View at: Publisher Site  Google Scholar
 F. Duan, Y. Yang, Y. Li et al., “Skull Identification via Correlation Measure Between Skull and Face Shape,” IEEE Transactions on Information Forensics and Security, vol. 9, no. 8, pp. 1322–1332, 2014. View at: Publisher Site  Google Scholar
 M. VaneZis, Forensic Facial Reconstruction Using 3D Computer Graphics: Evaluation and Improvement of Its Reliability in Identification, University of Glasgow, 2008.
 A. K. Moorthy, A. Mittal, S. Jahanbin, K. Grauman, and A. C. Bovik, “3D facial similarity: Automatic assessment versus perceptual judgments,” in Proceedings of the 4th IEEE International Conference on Biometrics: Theory, Applications and Systems (BTAS '10), pp. 1–7, IEEE, September 2010. View at: Publisher Site  Google Scholar
 W. E. Lorensen and H. E. Cline, “Marching cubes: a high resolution 3D surface construction algorithm,” ACM Siggraph Computer Graphics, vol. 21, no. 4, pp. 163–169, 1987. View at: Publisher Site  Google Scholar
 Q. Deng, M. Zhou, W. Shui, Z. Wu, Y. Ji, and R. Bai, “A novel skull registration based on global and local deformations for craniofacial reconstruction,” Forensic Science International, vol. 208, no. 1–3, pp. 95–102, 2011. View at: Publisher Site  Google Scholar
 Y. Hu, F. Duan, B. Yin et al., “A hierarchical dense deformable model for 3D face reconstruction from skull,” Multimedia Tools and Applications, vol. 64, no. 2, pp. 345–364, 2013. View at: Publisher Site  Google Scholar
 X. Ying, S. Q. Xin, Q. Sun, and Y. He, “An intrinsic algorithm for parallel poisson disk sampling on arbitrary surfaces,” IEEE Transactions on Visualization and Computer Graphics, vol. 19, no. 9, pp. 1425–1437, 2013. View at: Publisher Site  Google Scholar
 Q. Sun, L. Zhang, M. Zhang et al., “Texture brush: an interactive surface texturing interface,” in Proceedings of the 17th ACM SIGGRAPH Symposium on Interactive 3D Graphics and Games (I3D '13), pp. 153–160, March 2013. View at: Publisher Site  Google Scholar
 J. S. B. Mitchell, D. M. Mount, and C. H. Papadimitriou, “The discrete geodesic problem,” SIAM Journal on Computing, vol. 16, no. 4, pp. 647–668, 1987. View at: Publisher Site  Google Scholar  MathSciNet
 S. Xin and G. Wang, “Improving chen and han’s algorithm on the discrete geodesic problem,” ACM Transactions on Graphics, vol. 28, no. 4, 2009. View at: Publisher Site  Google Scholar
 X. Ying, S. Xin, and Y. He, “Parallel ChenHan (PCH) algorithm for discrete geodesics,” ACM Transactions on Graphics, vol. 33, no. 1, article 9, 2014. View at: Publisher Site  Google Scholar
 X. Ying, X. Wang, and Y. He, “Saddle vertex graph (SVG): a novel solution to the discrete geodesic problem,” ACM Transactions on Graphics, vol. 32, no. 6, article 170, 2013. View at: Publisher Site  Google Scholar
 V. Surazhsky, T. Surazhsky, D. Kirsanov et al., “Fast exact and approximate geodesics on meshes,” ACM Transactions on Graphics, vol. 24, no. 3, pp. 553–560, 2005. View at: Google Scholar
 C. Dorai and A. K. Jain, “COSMOSa representation scheme for 3D freeform objects,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 19, no. 10, pp. 1115–1130, 1997. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Junli Zhao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.