 Methodology
 Open Access
 Published:
A 3D cephalometric protocol for the accurate quantification of the craniofacial symmetry and facial growth
Journal of Biological Engineering volume 13, Article number: 42 (2019)
Abstract
Background
Cephalometric analysis is used to evaluate facial growth, to study the anatomical relationships within the face. Cephalometric assessment is based on 2D radiographic images, either the sagittal or coronal planes and is an inherently inaccurate methodology. The wide availability of 3D imaging techniques, such as computed tomography and magnetic resonance imaging make routine 3D analysis of facial morphology feasible. 3D cephalometry may not only provide a more accurate quantification of the craniofacial morphology and longitudinal growth, but also the differentiation of subtle changes in occlusion. However, a reliable protocol for the computation of craniofacial symmetry and quantification of craniofacial morphology is still a topic of extensive research. Here, a protocol for 3D cephalometric analysis for both the identification of the natural head position (NHP) and the accurate quantification of facial growth and facial asymmetry is proposed and evaluated. A phantom study was conducted to assess the performance of the protocol and to quantify the ability to repeatedly and reliably align skulls with the NHP and quantify the degree of accuracy with which facial growth and facial asymmetry can be measured.
Results
The results obtained show that the protocol allows consistent alignment with the NHP, with an overall average error (and standard deviation) of just 0.17 (9.10e6) mm, with variations of 0.21 (2.77e17) mm in the frontonasal suture and 0.30 (5.55e17) mm in the most prominent point in the chin. The average errors associated with simulated facial growth ranged from 1.83 to 3.75% for 2 years’ growth and from − 9.57 to 14.69% for 4 years, while the error in the quantification of facial asymmetry ranged from − 11.38 to 9.31%.
Conclusions
The protocol for 3D skull alignment produces accurate and landmark free estimation of the true symmetry of the head. It allows a reliable alignment of the skull in the NHP independently of userdefined landmarks, as well as an accurate quantification of facial growth and asymmetry.
Introduction
Cephalometric analysis is used to evaluate facial growth, to study the anatomical relationships within the face and as a routinely used tool for treatment planning in orthodontics and craniomaxillofacial deformity surgery. A standard cephalometric assessment is based on 2D radiographic images taken in either the sagittal (lateral cephalogram) or coronal planes (posteroanterior cephalogram), where multiple landmarks, lines and angles are identified to quantify vertical and horizontal relationships in the face. The analysis is primarily dependent on the accurate and repeatable definition of a standardized head position, therefore patients are imaged in customized image acquisition systems [1]. Regardless of the screening protocol, cephalometric images may be affected by magnification artefacts, craniofacial asymmetry and the superimposition of anatomical structures, all contributing to the imprecise and inaccurate traditional evaluation methods of facial morphology.
2D cephalometrics is an inherently inaccurate methodology, regarding the quantification of small increments of facial growth and the analysis of the transverse dimensions, in particular the analysis of facial asymmetry. The accuracy of contemporary cephalometric analysis is further complicated by the variation, inconsistency and errors in image capture, image quality and landmark identification, as well as limited intra and inter observer reliability and reproducibility [2]. The inability to detect slight changes in craniofacial morphology limits the clinical application of this technique and hinders its research usefulness in areas such as the quantification of longitudinal growth or the differentiation of subtle changes in occlusion, such as intercuspal and retruded contact position. No technique to date offers a precise virtual reproduction of dental occlusion [3], with study models and facebow recordings being the current gold standard in dentistry for the reproduction of interocclusal relationships.
The wide availability of 3D imaging techniques, such as computed tomography (CT), conebeam CT (CBCT) and magnetic resonance imaging (MRI), makes routine 3D analysis of facial morphology feasible. 3D reconstructed cephalometry possesses several advantages when compared with traditional cephalometry. It avoids the need for standardized fixation of the head and allows an a posteriori definition of the anatomical head position. It also prevents structure magnification and distortion inherent to imaging process, allowing at the same time the full assessment of the 3D craniofacial morphology, which is important for preoperative evaluation, surgical planning and postoperative evaluation [1]. 3D imaging techniques have the potential to provide significant precision improvements, however the accurate translation of wellestablished 2D concepts to a 3D framework needs further development. Studies in 3D cephalometry show conflicting results, with some studies concluding that it provides a more accurate evaluation of the craniofacial anatomy [4,5,6], while others conclude that there is no statistical difference between the two techniques [7, 8], or that it is less accurate than its 2D counterpart [9].
The most critical aspect of 3D cephalometric assessment is the accurate positioning of the skull in space, requiring accurate definition of both the craniofacial symmetry plane and the appropriate anatomical frame for the head. In 2D analysis of the lateral cephalogram, the skull is positioned through the identification of the true horizontal or the true vertical direction, and the identification of the S, as the reference point for image superimposition and comparison. The Frankfurt horizontal (FH) plane is used to orient the lateral cephalometric image in some specific analyses, requiring the identification of two bilateral landmarks (porion, Po, and orbitale, Or), which do not lie in the same plane as assessed using 3D imaging due to the natural asymmetry of the human face. In addition, FH related landmarks are amongst the most subjective landmarks and therefore can be prone to erroneous identification [5, 10]. The alternative is to define the true vertical direction, such as the NasionPogonion line [11]. This approach has the advantage of being based in the identification of the nasion (N) and pogonion (Pog), which are also less subjective than Po and Or [5], and may be regarded as lying in the symmetry plane of the head. These landmarks provide direct information about the midline structures of the face and are not subject to anatomical variation due to facial asymmetry. In addition, during childhood and adolescence, craniofacial growth is characterized by a relatively stable cranial base and foramen magnum, which remain relatively unchanged, when compared with the significant expansion of the cranial vault and distal and forward growth of both maxilla and mandible [1, 12, 13]. In a 3D setting, the basicranium can provide a reliable reference structure for further craniofacial growth assessment [1].
The development of a consistent and reliable 3D cephalometric workflow will provide clinicians with a valuable tool for the 3D characterization and quantification of facial growth, for the evaluation of dentofacial relationships and subsequent anatomically accurate orthodontic and surgical planning, and for the evaluation of treatment outcomes. We present a protocol for the accurate estimation of craniofacial symmetry, evaluating this using random skull orientations, facial growth vectors and different types of facial asymmetry. The null hypothesis tested was that the 3D method for analysing facial growth and asymmetry was not reliably quantifiable. The alternative hypothesis was that both craniofacial symmetry plane, facial growth and facial asymmetry could be accurately quantified with the proposed 3D algorithm.
Materials and methods
A skull model was obtained from the GrabCAD database (https://grabcad.com/), and skull landmarking (Table 1) conducted in Avizo v9.4.0 (Visualization Science Group). The 3D landmark coordinates were exported for further analysis and morphing with MATLAB R2007b and RStudio, respectively.
To simulate different types of facial growth and asymmetry, different geometrical transformations were applied to the reference skull. Several definitions are used to distinguish the different models, namely:

raw model: the initial/raw skull geometry obtained from the medical image data (Fig. 1a);

gold standard model: the raw model aligned according to the natural head position (NHP) and landmarked by a clinical expert (Fig. 1b). The NHP is defined as the position that the head adopts when the patient is sitting or standing, looking to the horizon or at a distant object [14];

target model: the gold standard model after random rigid transformation (Fig. 1e);

warped model: the gold standard model after mesh warping and random rigid transformation (Fig. 1f).
Figure 1 provides a graphical representation of the workflow from the raw model, to the definition of the gold standard model, target models and the final alignment. The use of the NHP aims to provide a common ground between standard 2D cephalometric evaluation and the protocol proposed in this work. It is important to note that all target models and warped models are generated from the gold standard model, since this contains the correct anatomical orientation provided by the clinical expert.
Computation of the symmetry plane
To determine the correct symmetry plane of the skull Principal Component Analysis (PCA) [15] and Iterative Closest Point (ICP) methods [16, 17] are used interchangeably. The craniofacial symmetry procedure is as follows:

1.
Reposition the original model in space according to the mesh centre of mass;

2.
Perform PCA over the model and reorient it according to the calculated eigenvectors λ_{i} ∈ R^{3};

3.
Compute the reflection that minimizes the rootmeansquared error (RMSE) of the Euclidean distances between the nearest points in the original and reflected models, along the eigenplanes defined in step 2;

4.
Apply ICP to register the original and reflected mesh created in step 3;

5.
Perform the PCA over the combined models from step 4 and reorient the original model according to the newly obtained eigenvectors \( {\boldsymbol{\lambda}}_i^{new}\in {R}^3 \).
In steps 1 and 2 a local reference frame that is independent of the acquisition frame is computed. Step 3 finds the first approximation to the symmetry plane; however, it is affected by the natural asymmetry of the skull. In step 4 the original and mirrored models are matched to create a new perfectly symmetrical model, which is used in step 5 to compute a new estimate of the symmetry plane. The original model is then oriented according to the \( {\boldsymbol{\lambda}}_i^{new}\in {R}^3 \), which is independent of any existing asymmetry in the skull model. It is important to note that step 4 is key to the accuracy of the procedure and is dependent on the alignment strategy [18]. A preliminary sensitivity test of the ICP method was therefore conducted. These results showed that pointtoplane alignment provides very accurate alignment between equal meshes (average error over 15 trials was just 8.81e9 (9.09e15) mm). ICP registration with pointtoplane alignment was therefore set as the default in step 4 (for more details the interested reader is referred to [16, 17]).
Positioning of the skull in the NHP
Evaluation of craniofacial morphology and growth is quantified considering a standardized anatomical position and therefore recovery of the NHP of the raw model may provide a consistent framework for 3D cephalometric assessment in clinical practice. The NHP can be defined by two landmarks along the symmetry plane, for example N and Pog, with the subsequent definition of the NPog line as the true vertical direction. However, it is important to note that use of the NPog line is not mandatory and two other landmarks along the symmetry plane could also be used to define the true vertical direction (for instance, N and A if the mandible is not available).
To assess the method’s reliability in recovering the NHP, a sensitivity test was performed. Rigid transformations were randomly generated and applied to the gold standard model (Fig. 1c). The rigid transformations consisted of random displacements −50 ≤ 〈d_{x}, d_{y}, d_{z}〉 ≤ 50 mm and rotations along the three orthogonal directions −50 ° ≤ 〈θ_{x}, θ_{y}, θ_{z}〉 ≤ 50°. The target model was then realigned with the gold standard model (Fig. 1g and i). This process was repeated 15 times and tested considering full (with cranial vault) and partial (no cranial vault) skull data, as in a standard CBCT acquisition process.
Craniofacial growth and asymmetry
Craniofacial growth was simulated by warping the viscerocranium of the gold standard craniofacial model while keeping the neurocranium unchanged. Figure 2a shows the average craniofacial growth of male and female subjects between 12 to 15 years old (adapted from [13]). In normal craniofacial development, an anterior growth of approximately 5.0 mm and distal growth of typically 20.0 mm was observed in subjects aged between 10 to 18 years of age [19]. For the purpose of this study average anterior and distal growth rates of 0.625 mm/year and 2.50 mm/year were considered, respectively. The positions of Cg, N, So, Or, Ec and Ba were fixed, while Go was displaced downwards and the A, Pog and Me landmarks were displaced both forward and downward to match the desired age (Fig. 2b). Similarly, facial asymmetry was included by applying different displacements to the craniofacial landmarks, such as N, IN, A, Ec, Go and Or. Three types of facial asymmetry were considered: (1) protrusion of the midface where IN, A, left and right Or were displaced forward (Fig. 3a); (2) unilateral protrusion of the face where IN, A, left Ec were displaced forward and the right Or was displaced backward (Fig. 3b); and (3) lateral displacement of the mandible where A, Pog, left and right Go were displaced laterally (Fig. 3c).
Mesh warping for both facial growth and facial asymmetry transformations was applied to the gold standard model and computed using two geometric morphometric packages, namely Geomorph [20] and Morpho [21]. Since the basilar surface remains relatively unchanged from adolescence to adulthood [12] and is independent of the landmarks used in cephalometric assessment, it was used as the reference anatomical structure to guide craniofacial alignment of all skulls.
Quantification of growth patterns and facial asymmetries is clearly dependent on accurate superposition of the craniofacial scans corresponding to the different stages of development. Alignment of morphologically different scans can be accomplished through the following steps:

1.
Computation of the symmetry plane as described in section 2.1;

2.
Definition of the NHP as described in section 2.2;

3.
Alignment refinement using IPC with weights (w_{i} ∈ R) assigned to the basilar surface (Fig. 1h and i).
By using ICP with weights, the error along the basilar surface is minimised. In pointtoplane alignment with weights, the weights are set to w_{i} ∈ {0, 1}, thereby acting as a binary variable that only allows alignment of the points belonging to the basilar surface region. It is important to note that only translations along the symmetry plane are considered, in order to preserve any existing asymmetries computed with the procedure proposed in section 2.1.
Evaluation measures
Evaluation of the cephalometric alignment for NHP and craniofacial growth and asymmetry was performed considering several measures, namely the alignment and distances between cephalometric landmarks against their respective landmarks on the gold standard model and general alignment measures such as the Hausdorff distance (Hd), where:
and the mean symmetrical distance (MSD)
where d stands for the 3D Euclidean norm and m and n are the number of points in A and B, respectively [22].
In facial growth and asymmetry tests, there is no longer an exact correspondence between cephalometric landmarks. Since mesh warping was applied to the gold standard model, the ideal alignment between gold standard and warped mesh was considered to be the initial alignment after warping and before any rigid transformation (Fig. 1d). In both cases, the goodness of alignment between the gold standard and the warped models was quantified by comparing the realigned models with their corresponding nontransformed gold standard models (initial alignment after warping). In the evaluation of alignment of morphologically different skulls (section 2.3), the Hd, MSD and the deviation between the N and Pog were computed. To understand how much information about facial growth and asymmetry is preserved after transformation and realignment, the distances between gold standard landmarks and morphed landmarks were therefore also calculated.
Results
The reliability in the recovery of the NHP was tested by applying a random rigid transformation (Fig. 1c) to the target model (Fig. 1g) and comparing its calculated NHP with the gold standard model (Fig. 1i). The process was repeated 15 times and the results are summarized in Table 2. The MSD and Hd (with standard deviations) between the gold standard and target models are 0.17 (9.10e6) mm and 0.52 (4.65e5) mm for the full data, respectively. The true vertical was recovered with an average error of just 0.21 mm for the N and 0.30 mm for the Pog landmarks. For partial data, the MSD and Hd were 0.27 (2.23e5) mm and 0.52 (4.26e5), whereas the N and Pog landmark error were 0.32 (2.63e7) and 0.23 (1.76e7), respectively. Figure 4 shows the original transformation and subsequent realignment between the gold standard model and the target model after estimation of the NHP for the full skull (Fig. 4a) and partial skull data (Fig. 4b).
The ability to quantify facial growth and asymmetry were tested by considering theoretical growths of 2 and 4years (Fig. 2a) and three different types of facial asymmetry (Fig. 3ac). The algorithm was able to realign the warped models with their corresponding gold standard geometries (Fig. 5ab). The NPog line was also recovered accurately, with an error of only 0.09 mm in the N and 0.16 mm in the Pog landmarks for the 2year growth. For the 4 years data the error was 1.67 mm for the N and 1.79 mm for the Pog. The algorithm also performs well in the presence of partial skull data (Fig. 5b).
Table 3 shows the differences between the morphed landmarks from the gold standard and the 2year and 4year data before and after transformation. There is an overall overestimation of the distances between the gold standard and target model(s). For instance, before rigid transformation the distance AA* between point A in the gold standard and the point A in the 2year (warped gold standard) model is 4.97 mm. After rigid transformation and realignment, the distance AA* is 5.06 mm, which corresponds to an overestimation of 0.091 mm. In the simulated 4year facial growth, the alignment error increases and the distance AA* is again overestimated by 1.61 mm (Table 3). The error between the theoretical and the quantified facial growth ranges between < 5% error for 2years to < 15% for 4years facial growth (Table 3).
The alignment of the asymmetric models showed good agreement with the theoretical models (Fig. 6ac), and good agreement was obtained even when the two landmarks N and Pog did not lie in the symmetry plane (Fig. 3c and 6c). Table 4 summarizes the distances between cephalometric landmarks for each independent facial asymmetry model. The errors in the quantification of facial asymmetry range from − 11.34% for asymmetry 2, to < 10% for asymmetry 1 and < 5% for asymmetry 3 (Table 4).
Discussion
3D medical imaging techniques offer the opportunity for more accurate evaluation of the craniofacial anatomy, however most assessment protocols are still based on the concepts of 2D cephalometry. Here, a protocol for general 3D cephalometric analysis is proposed, which aims firstly to accurately estimate facial symmetry and the NHP and secondly to use the basicranium to superimpose other morphologically different skulls. A phantom study was conducted to quantify the level of accuracy and reliability with which small morphological differences within or between subjects can be evaluated. Since the measured differences between the gold standard and warped models is one to two orders of magnitude smaller than the initial perturbation, the facial growth and asymmetry assessments were assumed to be accurate and therefore the null hypothesis was rejected.
According to the results obtained, the protocol provides: (1) consistent definition of the craniofacial symmetry plane and the true vertical direction, and (2) a method for consistent alignment between different skulls, which is independent of the standard cephalometric landmarks. The true vertical direction is given by two userdefined landmarks (here the N and Pog), however these are internally corrected once the true symmetry plane is computed. To consistently align different skulls (for instance, of the same subject but at different time points), the basilar surface is used as the reference structure, instead of the geometrical centre of the sella turcica (Sella point (S)). In previously reported studies, the estimation of the symmetry plane was accomplished through the identification of multiple anatomical landmarks. Landmarks that were then mirrored about an arbitrary plane and matched using Procrustes superimposition [23,24,25]. Here, a global solution for the definition of the symmetry plane is proposed, simply by taking advantage of the 3D skull data and the properties of both PCA and ICP methods.
The PCA aligns the skull according to the main axes of variation and can be made relatively insensitive to any natural facial asymmetry by reflecting the skull along the initial estimate of the symmetry plane (section 2.1 step 2) and reapplying the PCA method to the matched original and mirrored skulls (section 2.1 step 5). It is important to stress that the definition of the symmetry plane is only as accurate as the ability to match the original and the mirrored models [18]. To obtain a better match between original and mirrored skulls the ICP method is applied. A similar approach was proposed in [26] to estimate the facial symmetry plane. Here, the procedure is landmark free by testing the eigenvector λ_{i} that minimizes the rootmeansquared error (RMSE) of the Euclidean distances between the nearest points in the original and reflected models. The algorithm performed consistently for both full and partial skull data (Table 2).
In the definition of the NHP and with the newly defined symmetry plane, new estimates of the N and Pog landmarks and NPog line/plane are computed by estimating the nearest node to the landmark projection along the symmetry plane. However, it is important to stress that in the presence of a mandibular asymmetry (Fig. 3c) the identification of the nearest node along the symmetry plane would not achieve accurate results in the identification of the new Pog position. In this extreme scenario, the algorithm still produces an accurate estimate of the facial symmetry. This shows not only that the independence of the algorithm from userdefined landmarks, but also its usefulness for the management of patients with chin deviation [27]. In [26] an average angular deviation of 0.028° (0.014) was obtained between the estimated and the true symmetry plane. With the proposed procedure the maximum coronal deviation from the true symmetry is 0.30° (NHP recovery 15 trials, Table 2), and an upper limit of 0.36° of error in the coronal plane obtained for facial asymmetry 2 (Table 2).
In facial growth and asymmetry quantification, the final skull alignment is defined by the facial symmetry procedure and the ICP with weights, which minimizes distances across the basilar surface. Since in adolescence and adulthood there is a stabilization of the craniofacial base and foramen magnum [12], this anatomical region was considered suitable for 3D skull alignment. Similarly to the NHP recovery, the final alignment was independent of the userdefined positions of both the N and Pog landmarks. The advantage of using the basilar surface becomes evident in the presence of both facial growth (Fig. 5ab) and facial asymmetry (Fig. 6ac), in particular when the Pog landmark is not in the symmetry plane (Fig. 6c).
The alignment protocol allowed an accurate quantification of facial growth, with an error that ranged from − 11.38 to 9.31% for 2years and from − 11.38 to 2.77% for 4years facial growth (Table 3). Good skull alignment was also obtained when considering facial asymmetry. The linear measurements between landmarks showed that the error ranges from − 1.50 to 9.31% for asymmetry 1, from − 11.38 to 2.77% for asymmetry 2 and from 0.89 to 3.65% for asymmetry 3 (Table 4). The results show that the protocol might be useful to quantify slight asymmetrical variations in craniofacial anatomy, especially in patients with cleftlip and palate. Facial asymmetry in patients with cleftlip and palate typically range between 0.79 (0.23) mm in the lower face up to 1.15 (0.44) mm in the midface [28]. All the errors in asymmetry estimation are below the aforementioned values (Table 4).
The new protocol to determine the true craniofacial symmetry plane has numerous clinical applications, namely (1) as a part of a true 3D cephalometric analysis protocol to assess craniofacial morphology, growth and asymmetry as demonstrated above; (2) determination of the orientation of dental occlusion plane relative to the craniofacial morphology; (3) quantification of mandibular asymmetry and biomechanical investigations in masticatory asymmetry [29], (4) it may also be useful for virtual surgical planning of mandibular reconstructions [30], and mandibular implant development [31]. The protocol is landmark free and allows the reliable determination of the true craniofacial symmetry against a gold standard model. It can be combined with other midplane cephalometric landmarks (such as N and Pog) to position the skull in a standard anatomical position. The alignment error was quantified for both NHP, facial growth and facial asymmetry and was evaluated against the ideal alignment defined by a clinical expert. The results obtained demonstrate that symmetry alignment is accurate and independent of errors associated with the userdefined landmarks. In addition, the correct alignment of morphologically disparate skulls can be accomplished using the basilar surface instead of the traditional sella turcina (S) landmark. This is an important advance alongside the contemporary trend towards 3D imaging and the need to be able to reliably quantify longitudinal changes in 3D, with landmark free superimposition producing greater precision and reliability than previously used methods. With the craniofacial images being derivable from MRI, it is also noninvasive, clinically useful in orthodontics and maxillofacial surgery, versatile and applicable to other anatomical areas and even other fields of interest.
Conclusions
The protocol proposed for 3D skull alignment is fully automated and produces an accurate and landmark free estimation of the true symmetry plane of the human skull. It also allows the accurate alignment of the skull in the NHP, taking into account the true symmetry plane of the skull and two useddefined landmarks along the facial profile that specify the true vertical direction. After the definition of the true symmetry plane, the userdefined landmarks are internally corrected and the final alignment is completely independent of any userdefined data. The algorithm produces accurate quantifications of both facial growth and asymmetry, which renders it useful for clinical applications.
References
 1.
Swennen GR, Schutyser FA, Hausamen JE. In: Gwen RJ, Swennen FAC, Schutyser JEH, editors. Threedimensional cephalometry: a color atlas and manual. Berlin Heidelberg: Springer Science & Business Media; 2005.
 2.
Mish EC. Contemporary implant dentistry. St. Louis: Mosby Elsevier; 2008.
 3.
Solaberrieta E, Mínguez R, Barrenetxea L, Otegi JR, Szentpétery A. Comparison of the accuracy of a 3dimensional virtual method and the conventional method for transferring the maxillary cast to a virtual articulator. J Prosthet Dent Mosby. 2015;113:191–7.
 4.
Gribel BF, Gribel MN, Frazão DC, McNamara JA, Manzi FR. Accuracy and reliability of craniometric measurements on lateral cephalometry and 3D measurements on CBCT scans. Angle Orthod. 2011;81:28–37.
 5.
Zamora N, Llamas JM, Cibrián R, Gandia JL, Paredes V. A study on the reproducibility of cephalometric landmarks when undertaking a threedimensional (3D) cephalometric analysis. Med Oral Patol Oral Cir Bucal. 2012;17:678–88.
 6.
Li N, Hu B, Mi F, Song J. Preliminary evaluation of cone beam computed tomography in threedimensional cephalometry for clinical application. Exp Ther Med Spandidos Publications. 2017;13:2451–5.
 7.
Zamora N, Llamas JM, Cibrián R, Gandia JL, Paredes V. Cephalometric measurements from 3D reconstructed images compared with conventional 2D images. Angle Orthod. 2011;81:856–64.
 8.
Pittayapat P, Bornstein MM, Imada TSN, Coucke W, Lambrichts I, Jacobs R. Accuracy of linear measurements using three imaging modalities: two lateral cephalograms and one 3D model from CBCT data. Eur J Orthod. Oxford University Press. 2015;37:202–8.
 9.
van Vlijmen OJC, Maal T, Bergé SJ, Bronkhorst EM, Katsaros C, KuijpersJagtman AM. A comparison between 2D and 3D cephalometry on CBCT scans of human skulls. Int J Oral Maxillofac Surg. 2010;39:156–60.
 10.
Olszewski R, Lambrichts I, Pittayapat P, Bornstein MM, Odri GA, Jacobs R, et al. Threedimensional Frankfort horizontal plane for 3D cephalometry: a comparative assessment of conventional versus novel landmarks and horizontal planes. Eur J Orthod. Oxford University Press. 2017;40:239–48.
 11.
Mitchell L. An introduction to orthodontics. Oxford: OUP Oxford; 2013.
 12.
Premkumar S. Textbook of craniofacial growth. New Delhi: JP Medical Ltd; 2011.
 13.
Broadbent BH, Broadbent BHJ, Golden WH. Bolton standards of dentofacial developmental growth. Saint Louis: Mosby; 1975.
 14.
Proffit WR, Fields HW Jr, Sarver DM. Contemporary Orthodontics. 4th ed. St. Louis: Elsevier Health Sciences; 2006.
 15.
Dryden IL, Mardia KV. Statistical shape analysis: with applications in R. second Edi. Chichester: Wiley; 2016.
 16.
Besl P, McKay N. A method for registration of 3D shapes. IEEE Trans Pattern Anal Mach Intell. 1992;14:239–56.
 17.
Segal A, Haehnel D, Thrun S. GeneralizedICP. In: Brock O, Trinkle FR JC, editors. Robot Sci Syst, vol. 2009. Seattle: MIT Press; 2009. p. 435.
 18.
Cicconet M, Hildebrand DGC, Elliott H. Finding Mirror Symmetry via Registration and Optimal Symmetric Pairwise Assignment of Curves. 2017 IEEE Int Conf Comput Vis Work. IEEE; 2017. p. 1749–58. https://doi.org/10.1109/ICCVW.2017.206.
 19.
Hayashi I, Sakuda M, Takimoto K, Miyazaki T. Craniofacial growth in complete unilateral cleft lip and palate: a roentgenocephalometric study. Cleft Palate J. 1976;13:215–37.
 20.
Adams DC, Collyer ML, Kaliontzopoulou A. Geomorph: Software for geometric morphometric analyses. R package version 3.0.6. 2018. Available from: https://cran.rproject.org/package=geomorph
 21.
Schlager S. Morpho and Rvcg  Shape Analysis in R. In: Zheng G, Li S, Szekely G, editors. Statistical Shape and Deformation Analysis. Orlando: Academic Press; 2017. p. 217–56.
 22.
Preim B, Bartz D. Visualization in medicine: theory, algorithms, and applications. San Francisco: Elsevier; 2007.
 23.
Hajeer MY, Ayoub AF, Millett DT. Threedimensional assessment of facial softtissue asymmetry before and after orthognathic surgery. Br J Oral Maxillofac Surg. 2004;42:396–404.
 24.
Damstra J, Fourie Z, De Wit M, Ren Y. A threedimensional comparison of a morphometric and conventional cephalometric midsagittal planes for craniofacial asymmetry. Clin Oral Investig. SpringerVerlag. 2012;16:285–94.
 25.
Shin SM, Kim YM, Kim NR, Choi YS, Park SB, Kim YI. Statistical shape analysisbased determination of optimal midsagittal reference plane for evaluation of facial asymmetry. Am J Orthod Dentofac Orthop. Mosby. 2016;150:252–60.
 26.
Hartmann J, MeyerMarcotty P, Benz M, Häusler G, StellzigEisenhauer A. Reliability of a Method for Computing Facial Symmetry Plane and Degree of Asymmetry Based on 3Ddata. J Orofac Orthop / Fortschritte der Kieferorthopädie. Urban Vogel. 2007;68:477–90.
 27.
Baek SH, Cho IS, Il CY, Kim MJ. Skeletodental factors affecting chin point deviation in female patients with class III malocclusion and facial asymmetry: a threedimensional analysis using computed tomography. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. Mosby. 2007;104:628–39.
 28.
MeyerMarcotty P, Alpers GW, Gerdes ABM, StellzigEisenhauer A. Impact of facial asymmetry in visual perception: A 3dimensional data analysis. Am J Orthod Dentofac Orthop. Mosby. 2010;137:168.e1–8.
 29.
Moroi A, Ishihara Y, Sotobori M, Iguchi R, Kosaka A, Ikawa H, et al. Changes in occlusal function after orthognathic surgery in mandibular prognathism with and without asymmetry. Int J Oral Maxillofac Surg. 2015;44:971–6.
 30.
Ren W, Gao L, Li S, Chen C, Li F, Wang Q, et al. Virtual planning and 3D printing modeling for mandibular reconstruction with fibula free flap. Med Oral Patol Oral y Cir Bucal. 2018;23:e359–66.
 31.
Pinheiro M, Alves JL. The feasibility of a custommade endoprosthesis in mandibular reconstruction: implant design and finite element analysis. J CranioMaxillofacial Surg. 2015;43:2116–28.
Acknowledgements
Not applicable. (Same as in Funding.)
Funding
This research was supported by a Wellcome Trust Seed Award in Science to Xinhui Ma for the project Three Dimensional Analysis of Maxillofacial Growth in Patients with Cleft Lip and Palate (grant 205923/Z/17/Z).
Availability of data and materials
The skull geometry used in this study is available from the open access GrabCAD database (https://grabcad.com/).
Author information
Affiliations
Contributions
XM, MJF, PAM and GTM conceived the research programme and secured the funding. MP conceived and developed the 3D cephalometric protocol, undertook the validation studies and drafted the manuscript. XM and MJF provided advice on development of the methods and interpretation of the results. PL provided advise on the mathematical methods. PAM, GTM and GS provided clinical input and advice on clinical application. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Authors’ information
Not applicable.
Ethics approval and consent to participate
Not applicable. The skull geometry used in this study was obtained from the open access GrabCAD database (https://grabcad.com/).
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Pinheiro, M., Ma, X., Fagan, M.J. et al. A 3D cephalometric protocol for the accurate quantification of the craniofacial symmetry and facial growth. J Biol Eng 13, 42 (2019). https://doi.org/10.1186/s1303601901716
Received:
Accepted:
Published:
Keywords
 Craniofacial morphology
 Facial symmetry
 Facial growth
 Cephalometry
 Phantom study