Accuracy and reliability of orthogonal polynomials in representing corneal topography

Wang, Junjie ORCID: 0000-0003-4708-9447, Li, Xuefei, Wang, Zheng, Davey, Pinakin G, Li, Yiyu, Yang, Lanting, Lin, Mao, Zheng, Xiaobo, Bao, Fangjun and Elsheikh, Ahmed
(2022) Accuracy and reliability of orthogonal polynomials in representing corneal topography. Medicine in Novel Technology and Devices, 15. p. 100133.

Access the full-text of this item by clicking on the Open Access link.
[thumbnail of Junjie - Orthogonal polynomials in corneal topography 2022.docx] Text
Junjie - Orthogonal polynomials in corneal topography 2022.docx - Author Accepted Manuscript

Download (750kB)


Fitting of corneal topography data to analytical surfaces has been necessary in many clinical and experimental applications, yet absolute superiority of fitting methods was still unclear, and their overfitting risks were not well studied. This study aimed to evaluate the accuracy and reliability of orthogonal polynomials as fitting routines to represent corneal topography. Four orthogonal polynomials, namely, Zernike polynomials (ZPs), pseudo-Zernike polynomials (PZPs), Gaussian-Hermite polynomials (GHPs) and Orthogonal Fourier-Mellin polynomials (OFMPs), were employed to fit anterior and posterior corneal topographies collected from 200 healthy and 174 keratoconic eyes using Pentacam topographer. The fitting performance of these polynomials were compared, and the potential overfitting risks were assessed through a prediction exercise. The results showed that, except for low orders, the fitting performance differed little among polynomials with orders ≥10 regarding surface reconstruction (RMSEs ∼0.3 ​μm). Anterior surfaces of normal corneas were fitted more efficiently, followed by those of keratoconic corneas, then posterior corneal surfaces. The results, however, revealed an alarming fact that all polynomials tended to overfit the data beyond certain orders. GHPs, closely followed by ZPs, were the most robust in predicting unmeasured surface locations; while PZPs and especially OFMPs overfitted the surfaces drastically. Order 10 appeared to be optimum for corneal surfaces with 10-mm diameter, ensuring accurate reconstruction and avoiding overfitting. The optimum order however varied with topography diameters and data resolutions. The study concluded that continuing to use ZPs as fitting routine for most topography maps, or using GHPs instead, remains a good choice. Choosing polynomial orders close to the topography diameters (millimeters) is generally suggested to ensure both reconstruction accuracy and prediction reliability and avoid overfitting for both normal and complex (e.g., keratoconic) corneal surfaces.

Item Type: Article
Uncontrolled Keywords: 40 Engineering, 4003 Biomedical Engineering, Eye Disease and Disorders of Vision
Divisions: Faculty of Science and Engineering > School of Engineering
Depositing User: Symplectic Admin
Date Deposited: 15 Mar 2022 16:13
Last Modified: 20 Jun 2024 20:59
DOI: 10.1016/j.medntd.2022.100133
Open Access URL:
Related URLs: