|
|
||||||||
Scientific Section |
2nd Department of Oral and Maxillofacial Surgery, Faculty of Dentistry, Osaka University, 18, Yamadaoka, Suita, Osaka 565-0871, Japan
Abstract
In recent years, methods have been developed that calculate three-dimensional (3D) co-ordinates of orthodontic landmarks from lateral and frontal cephalograms. However, precise measurement has been impossible with these methods because, although they corrected the magnification of the image, they did not correct 3D cephalic malpositioning that occurs during the measurement of human subjects. In this study, we developed a 3D cephalometric system that corrected not only for magnification of the image, but also 3D cephalic malpositioning during cephalogram exposure.
Magnification of the image was corrected for first. Cephalic revolution was then sequentially corrected and divided into elements of x-, y-, and z-axes. The origin was parallelly translated to the mid-point of bilateral porion. In order to examine the accuracy of this system, seven human dry skulls were measured.
The accuracy unaffected by the cephalic revolution in any direction and standard errors was within 0·8 mm in any orthodontic landmarks. It was suggested that this measurement system would have sufficient accuracy for clinical application.
The results indicated that precise cephalometric measurement was possible with this system and it was suggested that its clinical application would be possible.
Key words: Three-dimensional analysis, Cephalometry, Stereograms, Mathematical Correction.
Introduction
A cephalogram assists orthodontic diagnosis and treatment planning. A variety of methods have been established for the analysis of lateral cephalograms (Tweed, 1946
; Downs, 1948
; Graber, 1952
; Ricketts, 1966
). However, frontal cephalograms have been limited in their use to the examination of facial asymmetry. This is largely due to distortion, mainly around the ear-rod, caused by cephalic revolution. In recent years, methods have been developed by which one can calculate three-dimensional (3D) co-ordinates of orthodontic landmarks from lateral and frontal cephalograms, for the quantitative evaluation of dento-facial deformities including facial asymmetry (Savara, 1965
; Yamazaki et al., 1981; Grayson et al., 1988
; Tasman et al., 1989; Bookstein et al., 1991
; Spolyar et al., 1993
). In comparison to computed tomography (CT), the cephalogram has advantages of low levels of exposure to radioactivity and no requirements for special equipment. However, if cephalic malpositioning and magnification of the images are not corrected for, precise measurements cannot be expected. In previous reports, accuracy has been examined by measuring human dry skulls. The correction of 3D cephalic malpositioning, which occurs during exposing cephalograms of human subjects, has not been described.
In this study, a 3D cephalometric system has been developed, which corrects not only magnification of the image, but also cephalic malpositioning. The accuracy of this system was evaluated by measuring dry skulls.
Methods
Development of the system
Composition of the system (Figure 1
).
Lateral and frontal cephalograms were exposed and traced respectively. Following input into an image processor (NEXUS6800, Kashiwagi Laboratory Co., Tokyo, Japan) mediated by a Charge Coupled Device (CCD) camera (DXC-930, Sony, Tokyo, Japan), landmarks were observed on a monitor, and located manually using a digitizer. The programme, which ran on a personal computer (PC-9801 As21, Nippon Electric Co., Tokyo, Japan) was developed in a BASIC language (N88-BASIC, Nippon Electric Co., Tokyo, Japan). Three-dimensional co-ordinates of landmarks were calculated from the data by this program.
|
In this study, we have devised a new methodology and constructed a measurement system that corrects for cephalic revolution in frontal and horizontal planes, in addition to correcting for cephalic revolution in the sagital plane as described by Yamazaki (1981).
Establishment of the co-ordinate system
The co-ordinate system used in this study was established as follows (Figure 2
). The plane determined by the midpoint (PoC) of bilateral porion (right side: PoR; left side: PoL) and bilateral orbita (right side: OrR; left side: OrL), is defined as the standard plane. The right hand co-ordinate system is decided by PoC fixed as the origin, by the straight line, which passes through PoC and OrC, as the x-axis (posteroanterior direction as positive), by the straight line, which passes through the origin and vertical to the standard plane as y-axis (downward as positive), and by the straight line, which is vertical to x- and y-axis as the z-axis (left direction as positive). This co-ordinate system is defined as the standard co-ordinate system.
|
|
The co-ordinate (x0, y0), which is corrected for the magnification of the image, is calculated by following numerical formulae:
![]() | (1) |
![]() | (2) |
yPA is calculated by the following numerical formula:
![]() | (3) |
A guide line y = yPA is indicated on the frontal cephalogram on the monitor and zPA is detected with reference to this line. Z0 is calculated by the following numerical formula:
![]() | (4) |
In order to calculate K for each landmark, an initial value K0 (Table 1
) is substituted for K, and z0 is calculated from numerical formulae above. Here, K is nothing else but z0. In addition, K0 is referred to as the average value of Japanese human dry skulls.
|
Second processing: method A (correction for cephalic revolution around the z-axis). Cephalic revolution around the z-axis or around the ear-rod often occurs. Yamazaki's correction method uses nasion as a reference point to calculate the revolution angle around the z-axis. However, the detection of nasion on the frontal cephalogram of a human subject is difficult. For this reason, our method adopts the use of upper incisal point (U1), which can be clearly detected from both lateral and frontal cephalograms as the reference. While Yamazaki's method adopted the cephalogram co-ordinate system, our method uses the above-mentioned standard co-ordinate system.
1. Calculation of cephalic revolution angle around the z-axis,
iS (Figure 4
).
x And y co-ordinate of U1 on an exposed subject are defined as x0(U1) and y0(U1), and the y co-ordinate of U1 on the frontal cephalogram as yPA(U1). The revolution angle around z-axis, iS (clockwise direction as positive), is summarized by the following numerical formula:
|
![]() | (5) |
However,
![]() |
2. Calculation of the z co-ordinate on the frontal cephalogram. yPA of point A is calculated using numerical formulae 1, 2, 3, and 5 and summarized by the following numerical formula:
![]() | (6) |
A guide line y = yPA is indicated on the frontal cephalogram on a monitor and zPA is detected with reference to this line.
3. Correction for magnification of the z co-ordinate. z0 is calculated using numerical formula 4 and summarized by the following numerical formula:
![]() | (7) |
4. Parallel translation of the origin. While the origin in the standard co-ordinate system is PoC, in order to minimize the influence of cephalic revolution around the y-axis on accuracy, OrC is used as a reference point for parallel translation for the z co-ordinate.
OrC on the frontal cephalogram is decided by detecting the intersection point of the straight line passing bilateral Or and the straight line passing the midpoint (PmC) of the most inside point of the bilateral orbita (right side: PmR, left side: PmL) and vertical to the straight line passing bilateral Or (Figure 5
). z0(OrC) is calculated by substituting zPA(OrC), which is the z co-ordinate on the frontal cephalogram, into numerical formula 7. The co-ordinate (X0, Y0, Z0) of point A on the standard co-ordinate system is calculated by the following numerical formula:
|
![]() |
Second processing: method B sequential correction for cephalic revolution, which is divided into elements of x-, y-, and z-axes. The method corrects for cephalic revolution around x- and y-axes sequentially, in addition to around the z-axis. After correcting for cephalic revolution around the z-axis as described in method A, the following calculation is carried out.
1. Correction for cephalic revolution around the x-axis.
The revolution angle around the x-axis, iF (Figure 6
anticlockwise direction as positive), is the angle between the straight line passing bilateral Or and the z-axis on the frontal cephalogram. It is calculated by the following numerical formula:
|
![]() |
The co-ordinate (x', y', z'), after correcting for revolution, is calculated by the following numerical formula:
![]() |
2. Correcting cephalic revolution around the y-axis.
The revolution angle around the y-axis, iH (Figure 7
anticlockwise direction as positive), is the angle between the straight line passing PoC and OrC, which is projected in the xz plane and the x-axis. It is described by the following numerical formula:
|
![]() |
In addition, as Po cannot be detected on the frontal cephalograms when measuring human subjects, the midpoint of the bilateral processus mastoideus is detected as a substitute point for zPA (PoC). This point is substituted into numerical formula 7 and z0 (PoC) is calculated. Moreover, z0 (OrC) is calculated in the same way as method A.
The co-ordinate (x'', y'', z'') after correcting revolution is described by the following numerical formula:
![]() |
3. Parallel translation of the origin. The co-ordinate (X0, Y0, Z0) of point A on the standard co-ordinate system is calculated by the following numerical formula:
![]() |
Examination of accuracy
Examination of the influence of cephalic revolution on measurement accuracy.
Lead leaves were affixed to 12 landmarks on a human dry skull (Table 1
). Measurement of the cephalogram of the human dry skull, which was exposed at positions revolved around x-, y-, and z-axes, was performed. The 3D co-ordinate of each landmark obtained with our system was compared with the true value obtained with the 3D digitizer (Tristation 400CNC, Nikon Co, Tokyo, Japan). Spatial accuracy within 5 µm was recorded. The accuracy of method A and B was compared.
Further examination of measurement accuracy
1. Calculation of the standard error.
Error is defined as the difference between a true value and a measured value. It includes both the systematic error and the random error. The systematic error is defined as a difference between the true value and the mean value. It is caused by the measurement system itself, such as characteristics of the cephalogram image, inputting of traced images into the image processor as mediated by a camera, the pixel size of the monitor, calculations in co-ordinates in the computer, and so on. A random error is unevenness in a measurement and it caused by human errors, such as tracing of the cephalogram, detection of landmarks, and so on. When accuracy is examined, both types of errors need to be considered.
Lead leaves were affixed to twelve landmarks on seven human dry skulls (Table 1
). The 3D co-ordinate of each landmark obtained with our system was compared with the true value. The standard error was calculated, as the linear sum of the systematic error and the random error, the confidence coefficient being 0·84.
![]() |
where u = square root of unbiased variance, n = number of measurements, u/
= random error, and
= systematic error.
Results
Construction of the 3D cephalometric system
A 3D cephalometric system was constructed. Magnification of the image and 3D cephalic malpositioning during cephalogram exposure were corrected and 3D co-ordinates of landmarks were calculated.
Accuracy
The influence of cephalic revolution on measurement accuracy.
Measurement of cephalograms of human dry skulls, which were exposed at positions revolved around x-, y-, and z-axes, were taken and the differences from true values calculated. The accuracy of methods A and B was compared (Figure 8
).
|
When the skull was revolved around the y-axis (5, 3, 0, 3, and 5 degrees), the difference from the true value in the x co-ordinate using method A increased with the degree of revolution, and became about 3 mm at 5 degrees revolution. Moreover, when tendencies among all landmarks were investigated, the accuracy of measurement of x and z co-ordinates of bilateral Or, Po, and Go, apart from the mid-sagittal plane, decreased with method A. However, when method B was applied, differences were within 0·7 mm for all co-ordinates.
When the skull was revolved around the z-axis 10_10 degrees (5-degree increments), differences from true values were within 0·8 mm for all co-ordinates with both methods.
Although the tendency differed with the various landmarks, the accuracy of method A was markedly influenced by cephalic revolution. However, the accuracy of method B was not influenced by cephalic revolution.
Further examination of measurement accuracy.
Further examination of the accuracy of method B was subsequently performed. The standard error was within 0·8 mm for the x co-ordinate, within 0·5 mm for the y co-ordinate, and within 0·8 mm for the z co-ordinate (Table 2
).
|
Since the introduction of cephalograms by Broadbent in 1931, various analytical methods have been established, most usually for lateral cephalograms. In recent years, methods that calculate 3D co-ordinates of orthodontic landmarks from lateral and frontal cephalograms have been developed. Since cephalic positioning is not usually accurate during cephalogram exposure, precise measurement requires correction for magnification of the image caused by the extension of X-ray bundles and 3D cephalic revolution.
A new methodology for the correction of measurement errors caused by 3D cephalic revolution (method B) has been developed and its measurement accuracy examined. Cephalograms of a human dry skull positioned to revolve around x-, y-, and z-axes were measured. Measurement accuracy was unaffected by cephalic revolution around any of the axes. Further examination of the accuracy was performed by evaluating standard errors in the measurement of human dry skulls. Seven human dry skulls were used and each subject was exposed seven times. The data obtained can be used for statistical procedures in accordance with a normal distribution based on the central limit theorem.
The error in this system is attributed to following factors: the error associated with tracing and manual inputting of landmarks, and distortion of the image when traced images is input into the image processor using a camera and displayed on a monitor. Although our system is considered to possess sufficient accuracy for clinical applications, it necessitates certain changes in the hardware, an image scanner substituted for the camera, an image processor and a monitor with smaller pixels, and the development of software capable of detecting landmarks automatically in order to improve measurement accuracy.
Conclusions
A 3D cephalometric system, which corrects magnification of the image and cephalic malpositioning in cephalograms three-dimensionally, was developed to improve the accuracy of cephalometric measurements for the purpose of clinical application. By this method, magnification of the image was corrected for first and cephalic revolution was subsequently sequentially corrected. In order to examine the accuracy of this system, measurements were made of seven human dry skulls. Results indicated that precise cephalometric measurement became possible with this new system and it was suggested its clinical application would be possible.
References
Bookstein F. L., Grayson B. H., Cutting C. B., Kim H. and McCarthy J. G. (1991) Landmarks in three dimensions: reconstruction from cephalograms versus direct observation, American Journal of Orthodontics and Dentofacial Orthopedics, 100(2), 133132.[Medline]
Broadbent B. H. (1931) A new X-ray technique and its application to orthodontia, Angle Orthodontist, 1, 4566.
Downs W. B. (1948) Variations in facial relationship: their significance in treatment and prognosis, American Journal of Orthodontics, 34, 812840.
Graber T. M. (1952) New horizons in case analysisclinical cephalometrics, American Journal of Orthodontics, 38, 603624.
Grayson B. H., Cutting C. B., Bookstein F. L., Kim H. and McCarthy J. G. (1988) The three-dimensional cephalogram: theory, technique, and clinical application, American Journal of Orthodontics and Dentofacial Orthopedics, 94(4), 327337.[Medline]
Ricketts R. M. (1966) Cephalometric synthesis, American Journal of Orthodontics, 46, 647673.
Savara B. S. (1965) A method for measuring facial bone growth in three dimensions, Human Biology, 37, 245255.[Medline]
Spolyar J. L., Vasileff W. and Macintosh R. B. (1993) Image corrected cephalometric analysis (ICCA): design and evaluation, Cleft Palate-Craniofacial Journal, 30(6), 528539.
Tasman B. and Amanda H. A. (1989) Computer assisted location of reference points in three dimensions for radiographic cephalometry, American Journal of Orthodontics and Dentofacial Orthopedics, 95(6), 490498.[Medline]
Tweed C. H. (1946) The Frankfort-mandibular plane angle in orthodontics, American Journal of Orthodontics, 32, 175230.
Yamazaki T. and Iwasawa T. (1981)A study on the three-dimensional position and morphology of the mandible in cephalic regionwith reference to comparison between the normal occlusion and true mandiblar protrusion groups,Journal of Japan Orthodontic Society, 40(2), 205227. [In Japanese]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |