A Novel Three-Dimensional Vector Analysis of Axial Globe Position in Thyroid Eye Disease
Purpose. To define a three-dimensional (3D) vector method to describe the axial globe position in thyroid eye disease (TED). Methods. CT data from 59 patients with TED were collected and 3D images were reconstructed. A reference coordinate system was established, and the coordinates of the corneal apex and the eyeball center were calculated to obtain the globe vector . The measurement reliability was evaluated. The parameters of were analyzed and compared with the results of two-dimensional (2D) CT measurement, Hertel exophthalmometry, and strabismus tests. Results. The reliability of measurement was excellent. The difference between and 2D CT measurement was significant (), and was more consistent with Hertel exophthalmometry than with 2D CT measurement (). There was no significant difference between and Hirschberg test, and a strong correlation was found between and synoptophore test. When one eye had a larger deviation angle than its fellow, its corneal apex shifted in the corresponding direction, but the shift of the eyeball center was not significant. The parameters of were almost perfectly consistent with the geometrical equation. Conclusions. The establishment of a 3D globe vector is feasible and reliable, and it could provide more information in the axial globe position.
Globe position is a key element in the diagnosis and treatment of orbital diseases, especially thyroid eye disease (TED). Exophthalmos is the most widely used indicator in TED. Although different measurement methods have been proposed, all of them assess the axial position of the corneal apex and represent the change in the anterior-posterior direction only. Additionally, the computed tomography (CT) exophthalmos measurement, which is considered the reference standard, needs further evaluation for its reliability and accuracy [1–5]. Furthermore, the indicators to assess the axial globe position in the superior-inferior and medial-lateral directions are also important, but in TED, they are not uniform. Alsuhaibani et al., Takahashi and Kakizaki, and Fichter et al. have all reported that the globe position changes horizontally in TED patients, based on different measurement methods [6–8].
Strabismus is not uncommon in TED, meaning that the eyeballs may have a tilt change in these patients, in addition to protrusion or translation. On the one hand, there might be some deviation if we define the position of the whole eyeball by the corneal point only; on the other hand, it could not describe the orientation of the globe. As mentioned above, it might be better to evaluate the globe position in three-dimensional (3D) space using a vector, which could show the location as well as the orientation.
Describing the globe position with a 3D vector method could be also useful for orbital decompression surgery in TED. Having a more accurate measurement method could help evaluate postoperative changes of the axial globe position in different directions. What is more, new onset strabismus is an important complication after orbital decompression surgery, and the 3D vector method might help in studying the eyeball motion and tilt change.
The 3D reconstruction technique provides a more comprehensive view of the orbit and improves the ability to analyze in orbital diseases [9, 10]. In the following study, we established a coordinate system based on 3D CT images, defined the globe vector, and discussed the reliability of this method. We also compared the globe vector with the results of traditional two-dimensional (2D) CT measurement, Hertel exophthalmometry, and clinical strabismus tests.
2. Patients and Methods
Eighty-one sets of CT data from 59 patients diagnosed with TED were collected between 2012 and 2016 at the Eye and ENT Hospital of Fudan University. There were 30 females and 29 males, and the average age was 48.6 ± 11.3 years (ranging from 21 to 68 years). The exclusion criteria included the patients with orbital rim or zygomatic change from fracture or surgery, orbital tumor, and orbital inflammation. The research adhered to the tenets of the Declaration of Helsinki and was approved by the Ethics Committee of Fudan University. Informed consent was obtained from all the patients.
2.2. CT Scan and 2D Exophthalmos Measurement
CT scans were obtained at contiguous 0.75 mm thickness (SOMATOM Sensation, Siemens AG, Germany). The patients were instructed to look straight with fixed eyes during scanning. CT data were recorded in DICOM format (Digital Imaging and Communications in Medicine) and archived on compact discs, then imported to MIMICS 16.0 software (Materialise Dental, Leuven, Belgium) for image analysis.
The axial image including the thickest lens was chosen, and a line between the lateral orbital rims was drawn; then, the perpendicular distance between the corneal apex and this line was defined as the 2D CT exophthalmos.
2.3. Establishment of Globe Vector
The volume-rendering image was reconstructed in MIMICS, and then a 3D coordinate system was established based on the landmarks and reference planes described in Table 1. The landmarks were pointed to the volume-rendering image and adjusted on the axial, coronal, and sagittal images. The intersection of the three reference planes was set as the origin (0, 0, 0), and positive x, y, and z coordinate values indicated the front, left, and superior orientations, respectively (Figure 1).
The points of the corneal apex (C) and the eyeball center (E) were marked using the following procedure. The cursor was adjusted and positioned on the middle point of the cornea section in the axial and sagittal images simultaneously; this was then named as point C. A sphere was created with the CAD (computer-aided design) function and was made to overlap with the eyeball wall on the axial, sagittal, and coronal images; then, the center of this sphere was named point E. A 3D vector was defined to present the position of the globe in the reference coordinate system.
The coordinate of in the reference coordinate system could be calculated with mathematical methods and expressed as an actual size in millimeters: the STL coordinates of the reference landmarks were recorded, and the equations of the three reference planes could be calculated based on these coordinates as in Table 2. The distance (d) from the points C and E to the coronal, midsagittal (MS), and Frankfort horizontal (FH) planes was equal to their coordinates in the reference coordinate system and recorded as (xC, yC, and zC) and (xE, yE, and zE), respectively: here, , , , and represent the constant coefficients of the plane equation, and , , and represent the STL coordinates of points E and C.
The angle of and the FH plane () or the MS plane () was calculated as follows: here, represents the normal vector of the reference plane.
The angle of and axes , , and (, , and ) was calculated as follows:
All calculations were preset using Microsoft Excel software, and the results could be automatically acquired after inputting the original STL coordinate data. All the measurements were made by the same oculoplastic specialist.
2.4. Clinical Examination
Hertel exophthalmometry was performed by an experienced ophthalmologist using the standard method. Hirschberg corneal reflex test was performed before CT scanning by the same ophthalmologist; additionally, the objective strabismus angles were measured by synoptophore and recorded.
2.5. Statistical Analysis
Statistical analysis was performed using SPSS. Mean values with standard deviations (SD) were given. Twenty patients were randomly selected for interobserver and intraobserver reliability analysis. 2D CT measurement, landmark positioning, and calculation were repeated by the same observer as well as by another observer without referring to the previous results. The absolute difference between the two measurements was analyzed. Intraclass correlation coefficient (ICC) and Bland-Altman plot were used to evaluate method error. A paired t-test was used to assess the paired data, and a nonparametric test was used to compare the other data. Pearson correlation was used to analyze the correlation between data. Linear regression analysis was used to analyze the relativity. was considered statistically significant.
3.1. Measurement Reliability Analysis
For repeated measurements on the same scan, the intraobserver and interobserver agreements are listed in Table 3. Among the patients, 13 underwent unilateral orbital decompression surgery and were given a CT examination pre- and postoperatively within 2 weeks. All 13 of these patients were in the inactive phase, and we assumed that their nonoperated eyes would remain in the same positions; the average absolute difference in coordinates between the two scans ranged from 0.24 to 0.60 mm, and the average absolute difference was 1.32 mm for 2D exophthalmos measurements. The absolute difference of the 2D exophthalmos between the two scans was significantly larger than (p ranged from 0.001 to 0.008).
3.2. Parameter Analysis
The vector of the two eyes could be visualized with software and observed and compared from different orientations (Figure 2). For the 61 cases without a postsurgery change, the coordinates and deviation angles of are listed in Table 4.
3.3. and 2D CT Measurement
The coordinate xC of represented the exophthalmos theoretically. The difference between xC and the 2D exophthalmos measurement ranged from −2.5 to 3.2 mm, and there was a significant difference between them (). The absolute difference was 1.0 ± 0.73 mm, and there were 29 of 81 (35.8%) cases which had an absolute difference of 1 mm or more.
3.4. and Hertel Exophthalmometry
The exophthalmos measured by Hertel exophthalmometer was 19.1 ± 3.6 mm (ranged from 9 to 28 mm). The absolute difference between the Hertel measurement and 2D CT exophthalmos was 1.39 ± 0.97 mm, significantly larger than the difference between the Hertel measurement and xC (0.95 ± 0.67 mm, ), and the correlation between the Hertel measurement and xC was slightly stronger () than the correlation between the Hertel measurement and 2D CT exophthalmos (). The Bland-Altman plots showed that the difference between the Hertel measurement and 2D CT exophthalmos was more dispersed than that between the Hertel measurement and xC and had much fewer plots within −1~1 mm (Figure 3).
3.5. and Clinical Strabismus Angle
Hirschberg test was performed by an experienced ophthalmologist, and the results ranged from 0 to 45° (15.3° ± 12.4°). The direction of strabismus was compared with the vector orientation, and no significant difference was found between them (). The values of were also compared with the angles of Hirschberg test, and there was no significant difference between them (); the correlation between them was strong ().
The synoptophore test was performed in 37 cases to record the objective deviation angles: the horizontal deviation degrees ranged from −10 to 40° (10.1° ± 13.2°), with a positive value representing esotropia and a negative value representing exotropia; the vertical deviation ranged from 0 to 35° (9.2° ± 10.2°), with a positive value representing hypertropia and a negative value representing hypotropia. The values of were compared with the vertical deviation angles tested by synoptophore, and there was no significant difference between them (); the correlation between them was strong (). The correlation between the values of and the horizontal deviation angles was strong (), but the angles tested by synoptophore were significantly larger than those tested by (including both positive and negative values, ).
3.6. and Eyeball Movement
The values between the two eyes were compared, and it was found that in the eyes with a larger (including positive and negative values), the absolute values of yC were significantly smaller than the values of their fellows (); in the eyes with a smaller , the absolute values of yC were significantly larger than those of their fellows (), and there was a strong negative correlation between the difference of and the difference of absolute values of yC (, ). The values of were also compared: in the eyes with a larger , the values of zC were significantly larger than those of their fellows (); in the eyes with a smaller , the values of zC were significantly smaller than those of their fellows (), and there was a strong correlation between the difference of and the difference of zC (, ). However, no significant difference was found in the absolute values of yE or the values of zE between the two eyes (p ranged from 0.08 to 0.96).
In theory, the relation between points E and C should be as follows: ; ; and (Figure 4). We calculated the theoretical values of xC, yC, and zC based on these equations, and there was no significant difference between the theoretical xC and xC (, ), the theoretical yC and yC (, ), or the theoretical zC and zC (, ). The linear regression equations of the theoretical xC, yC, and zC and their actual values were all y = x ( with no constant, ) (Figure 5).
In previous research on TED, there have been few references to the change in the axial globe position in vertical and lateral orientations, and measurement methods have not been uniform: horizontal eye positions are sometimes evaluated using the interpupillary distance  and sometimes using the distance between the most medial part of the eyeballs measured on an axial CT image . Fichter et al. examined the horizontal position of the eyeball by measuring the distance between the nasion and the pupil using a digital pupillometer . Alsuhaibani et al. measured the vertical difference between the two eyeballs on coronal CT images to evaluate the vertical position change . However, nearly all of these methods only evaluated the relative position between the two eyes.
Based on 3D reconstructed images, we tried to describe the axial globe position in more directions. The establishment of a craniofacial 3D coordinate system has been described in many reports [11–14]; in this study, we chose the lateral orbital rim to set the coronal plane because it could make the coordinates consistent with the clinical exophthalmos measurement. The reproducibility of the reference coordinate system was proven stable, which is consistent with previous studies [11, 15]. The coordinates of were obtained by preset mathematical calculation rather than measuring images, which made the process more accurate and convenient. Besides the same CT scan, we assessed the measurement reliability of different CT scans, which has not been discussed before. Our results indicated that the reproducibility of was significantly better than the 2D CT measurement. Additionally, the coordinate xE, which represented the exophthalmos, was significantly different from exophthalmos values obtained from the 2D CT images, and in more than one third of cases differed by more than 1 mm, which was considered an accepted error between the repeated exophthalmometry by different observers . The difference between Hertel exophthalmos and 2D CT measurement was found to be variable in previous reports [3, 17, 18]. In our study, the consistency between and Hertel exophthalmos was better than consistency between 2D CT measurement and Hertel exophthalmos. This might be because the stability of the 3D coordinate system reduces head position change and measurement error, especially when following up patients using different CT scans.
Vector analysis has been used in some research [11, 19, 20], but this study represents the first time that a vector has been used to define the position of the eyeball in TED. A 3D vector can give a more comprehensive and intuitive view of the globe position. In previous reports, the mean values of 2D CT exophthalmos ranged from 21.2 mm to 22.5 mm in TED [1, 3, 21]; in our study, the mean value was 19.5 mm. The mean distance between the nasion and the pupillary center was 33.8 mm in TED , and the mean distance between the corneal apex and the MS plane was 32.2 in our study. The mean vertical height of the corneal apex from the FH plane was 20.8 mm in our study. Additionally, we have also listed the location of the eyeball center in Table 4. These results might give a new reference for the axial globe position research in TED.
The vector can describe the deviation of the eyeball. We assumed that approximated the ocular axis and that its normal position was approximately parallel to the x-axis; the angles of then indicated the severity of eyeball deviation in different directions. The orientation and angle of were largely in accordance with Hirschberg test; however, depending on the experience of the ophthalmologist, the Hirschberg test gave a rough value while could give a more exact deviation angle. The vertical deviation of was also in accordance with the synoptophore results. However, the horizontal deviation tested by synoptophore was significantly larger than that of the in the esotropia patients and smaller in the exotropia patients, although there was a very strong correlation between them. This difference might be because the synoptophore was affected by the artificial viewing conditions, resulting in overconvergence . Although there were some other factors that might have affected the clinical strabismus examination and measurement, such as the angle kappa, convergence, test distance, and measurement error, our study still indicated that, to some extent, the deviation of the eyeball as measured by morphological method was consistent with that found in clinical examination and might better reflect the anatomical position of the eyeball.
When an eye had a larger than its fellow, meaning that the eyeball had a more severe medial tilt; its corneal apex was closer to the nasal side than its fellow. When an eye had a larger , meaning that the eyeball had a more severe superior tilt, its corneal apex was higher than its fellow. However, there was no significant change in the positions of the eyeball center. The results indicated that the tilt of eyeball mainly depended on a rotation motion around its center. The fibrotic extraocular muscles pulled the eyeball [23, 24] and made it rotate around the eyeball center, making the corneal apex change its position in addition to protruding, but the lateral or vertical shift of the eyeball center was not obvious. Of course, the influence came from extrusion of the expanded orbital tissue, and the shift during ocular ductions could not be excluded completely . Additionally, due to changes in orbital bone, fat, and extraocular muscles, eyeball motion can be complicated during TED surgery, possibly leading to changes in strabismus or exophthalmos [26, 27]. The vector method might help us to better study the effects of surgery better in TED patients, including eyeball shift and rotation.
To further discuss the problem mentioned above, the relativity between the corneal apex and the eyeball center was analyzed. The relative position of the corneal apex to the eyeball center was influenced by the deviation angles of , suggesting that the bilateral corneal apexes might be asymmetrical even if the eyeball centers were symmetrical. The deviation of the corneal apex could disturb the exophthalmos measurement both using Hertel exophthalmometer and 2D CT images; the vector would improve this problem.
Although there were some drawbacks to the vector method, for example, it could not reflect the globe rotation around the ocular axis, it still supplied a novel and useful way of measuring the axial globe position.
In conclusion, the establishment of a 3D coordinate system and globe vector is a feasible and reliable method of depicting the axial globe position in TED. The vector method can depict the globe position in a more intuitive and comprehensive way and can provide more information about the eyeball location and its deviation than traditional methods. Additionally, vector analysis might help evaluate eyeball motion more accurately after decompression surgery in TED.
Conflicts of Interest
No conflicting relationship exists for any author.
The project was supported by the Shanghai Committee of Science and Technology, China (Grant no. 14411961800), and Shanghai Health and Family Planning Commission (Grant no. 20164Y0144).
C. M. Bingham, J. A. Sivak-Callcott, M. J. Gurka, J. Nguyen, J. P. Hogg, and S. E. Feldon, “Axial globe position measurement: a prospective multicenter study by the International Thyroid Eye Disease Society,” Ophthalmic Plastic and Reconstructive Surgery, vol. 32, no. 2, pp. 106–112, 2015.View at: Publisher Site | Google Scholar
E. Nkenke, T. Maier, M. Benz et al., “Hertel exophthalmometry versus computed tomography and optical 3D imaging for the determination of the globe position in zygomatic fractures,” International Journal of Oral and Maxillofacial Surgery, vol. 33, no. 2, pp. 125–133, 2004.View at: Publisher Site | Google Scholar
M. Segni, G. B. Bartley, J. A. Garrity, E. J. Bergstralh, and C. A. Gorman, “Comparability of proptosis measurements by different techniques,” American Journal of Ophthalmology, vol. 133, no. 6, pp. 813–818, 2002.View at: Google Scholar
N. Fichter, H. Krentz, and R. F. Guthoff, “Functional and esthetic outcome after bony lateral wall decompression with orbital rim removal and additional fat resection in Graves’ orbitopathy with regard to the configuration of the lateral canthal region,” Orbit, vol. 32, no. 4, pp. 239–246, 2013.View at: Publisher Site | Google Scholar
F. Borumandi, B. Hammer, H. Noser, and L. Kamer, “Classification of orbital morphology for decompression surgery in Graves’ orbitopathy: two-dimensional versus three-dimensional orbital parameter,” The British Journal of Ophthalmology, vol. 97, no. 5, pp. 659–662, 2013.View at: Publisher Site | Google Scholar
A. C. D. O. Ruellas, C. Tonello, L. R. Gomes, M. S. Yatabe, L. Macron, and J. Lopinto, “Common 3-dimensional coordinate system for assessment of directional changes,” American Journal of Orthodontics and Dentofacial Orthopedics, vol. 149, no. 5, pp. 645–656, 2016.View at: Publisher Site | Google Scholar
K. Wang, Y. Yang, F. Shen, J. Tao, H. Xu, and J. Portnof, “Utilization of three-dimensional computed tomography for craniofacial phenotypic analysis in children with velocardiofacial syndrome,” The Journal of Craniofacial Surgery, vol. 20, no. 6, pp. 2013–2019, 2009.View at: Publisher Site | Google Scholar
A. Muramatsua, H. Nawa, M. Kimura, K. Yoshida, M. Maeda, and A. Katsumata, “Reproducibility of maxillofacial anatomic landmarks on 3-dimensional computed tomographic images determined with the 95% confidence ellipse method,” The Angle Orthodontist, vol. 78, no. 3, pp. 396–402, 2007.View at: Publisher Site | Google Scholar
S. Thapa, A. K. Gupta, A. Gupta, V. Gupta, P. Dutta, and R. S. Virk, “Proptosis reduction by clinical vs radiological modalities and medial vs inferomedial approaches comparison following endoscopic transnasal orbital decompression in patients with dysthyroid orbitopathy,” JAMA Otolaryngology. Head & Neck Surgery, vol. 141, no. 4, pp. 329–334, 2015.View at: Publisher Site | Google Scholar
R. A. Clark, J. M. Miller, and J. L. Deme, “Location and stability of rectus muscle pulleys muscle paths as a function of gaze,” Investigative Ophthalmology & Visual Science, vol. 38, no. 1, pp. 227–240, 1997.View at: Google Scholar