Using these formulas, we can compute the orientation angle of both anterior and posterior corneal aberrations. Therefore, the orientation angle difference between anterior and posterior corneal aberrations can be calculated. In this article, Δβ, defined as the angular difference, is used to represent the orientation angle difference between the anterior aberration vector
(\mathbf|<||<|\FINLINE-FIGLINK LOCATOR="cjs1164"/INLINE-FIG/F|<|A|>||>||>|) and the posterior aberration vector
(\mathbf|<||<|\FINLINE-FIGLINK LOCATOR="cjs1164"/INLINE-FIG/F|<|P|>||>||>|) . The angular difference was calculated and rescaled in the range of 0° to 180° by the following formula:
\[{\Delta}{\beta}{=}\left\{\begin{array}{c}m\ {\cdot}\ {\vert}{\theta}_{\mathrm{A}}{-}{\theta}_{\mathrm{P}}{\vert}\ \mathrm{if}\ m\ {\cdot}\ {\vert}{\theta}_{\mathrm{A}}{-}{\theta}_{\mathrm{P}}{\vert}{\leq}180{^\circ}\\360{-}m\ {\cdot}\ {\vert}{\theta}_{\mathrm{A}}{-}{\theta}_{\mathrm{P}}{\vert}\ \mathrm{if}\ m\ {\cdot}\ {\vert}{\theta}_{\mathrm{A}}{-}{\theta}_{\mathrm{P}}{\vert}{>}180{^\circ}\ ,\end{array}\right.\]
where θ
A and θ
P are the orientation angles of anterior aberration vector
(\mathbf|<||<|\FINLINE-FIGLINK LOCATOR="cjs1164"/INLINE-FIG/F|<|A|>||>||>|) and posterior aberration vector
(\mathbf|<||<|\FINLINE-FIGLINK LOCATOR="cjs1164"/INLINE-FIG/F|<|P|>||>||>|) , respectively. The use of Δβ to represent the angular difference between the two vectors allowed uniformity of range from 0 to 180° for different sets of aberrations, thus facilitating easy comparison. As an example, for two aberrations with opposite directions, although |θ
A− θ
P| equals 90°, 60°, 180°, 45°, 90°, 36°, 60°, and 180° for astigmatism (
m = 2), trefoil (
m = 3), coma (
m = 1), quadrafoil (
m = 4), secondary astigmatism (
m = 2), pentafoil (
m = 5), secondary trefoil (
m = 3), and secondary coma (
m = 1) respectively, the value of Δβ is always 180°. To calculate the compensation effect between anterior and posterior corneal aberration, we decomposed the posterior corneal aberration
(\mathbf|<||<|\FINLINE-FIGLINK LOCATOR="cjs1164"/INLINE-FIG/F|<|P|>||>||>|) into a component (
F), defined as the compensation component, whose direction was parallel to the direction of the anterior corneal aberration vector
(\mathbf|<||<|\FINLINE-FIGLINK LOCATOR="cjs1164"/INLINE-FIG/F|<|A|>||>||>|) . The following formula was used to calculate the compensation component:
\[F{=}M_{\mathrm{P}}\ {\cdot}\ \mathrm{cos}\ ({\Delta}{\beta})\ \mathrm{if}\ m{\neq}0\]
and
\[F{=}(C_{n}^{0})_{\mathrm{p}}\ \mathrm{if}\ m{=}0,\]
where
F was the compensation component,
M P was the magnitude of the posterior corneal aberration vector, and (
C n 0)
p was the Zernike coefficients of the posterior corneal aberration. For spherical aberration (
m = 0), both anterior and posterior aberrations were not vectors and
F was equal to the Zernike coefficient of posterior corneal aberration. Using coma as an example,
Figure 2explains the method for calculating the compensation component (
F). For each patient and specific aberration, the percentage of the anterior corneal aberration compensated by posterior cornea was obtained by calculating the ratio of
F and the magnitude of the anterior corneal aberration,
M A (when
m ≠ 0), or the ratio of
F and the anterior corneal aberration Zernike coefficient, (
C n 0)
A (when
m = 0). To get the average compensation factor (
k), we computed the the linear regression (with 0
y-axis interception) between
F and the anterior corneal aberration
M A (for aberrations with
m ≠ 0) or between
F and the anterior corneal aberration Zernike coefficient (
C n 0)
A (for aberrations with
m = 0) for all subjects. The slope (
k) of the fitted line indicated that the average compensation factor and the determination factor (
R 2) was a measure of how strictly this percentage was followed in each subject in the investigated group.