Fourier Tools for the Evaluation of Refractive Multifocal Designs

doi:10.21203/rs.3.rs-3332165/v1

Download PDF

Article

Fourier Tools for the Evaluation of Refractive Multifocal Designs

https://doi.org/10.21203/rs.3.rs-3332165/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 19 Dec, 2023

Read the published version in Scientific Reports →

You are reading this latest preprint version

Tools and methods for the theoretical evaluation of optical properties of refractive multifocal commercial designs are presented, which use segmentation techniques to divide the pupil and classical Fourier optics. They can be applied to refractive designs such as multifocal contact lenses and intraocular lenses as well as to the study of small aperture designs and corneal inlays. In this paper, the optical profiles of 12 multifocal, commercially available, contact lenses from four companies were measured with the help of the NIMO equipment (LAMBDA-X). The corresponding through-focus Visual Strehl ratios were computed in the frequency domain to quantify image quality as a single scalar value. The multifocal properties of each of the commercial contact lenses were evaluated for three different base prescriptions (-6.00 D, -3.00 D, and + 1.00 D) in 65 patients. The contact lenses of this study were also evaluated in conjunction with the higher-order aberrations obtained from 65 eyes collected by using a standard Shack-Hartmann sensor. The diameters used for the evaluation of the contact lenses varied according to each eye’s pupil size (from 2.45 to 6.27 mm). The proposed tools have allowed us to evaluate multifocal solutions in a theoretical manner without the necessity of creating prototypes. They helped us differentiate between lenses created for young and old presbyopes, instructed the proper utilization of hyperfocal distance in the field of contact lens design, and demonstrated the importance of the effective aperture effect. This is the first proposed method that converts sagittal powers of multifocal solutions into their equivalent optical quality metric. It can be applied to the evaluation of myopia control methods. It is the author’s hope that the reader appreciates the tools presented and uses them to further the field of refractive multifocality.

Health sciences/Health care

Health sciences/Medical research

Physical sciences/Optics and photonics

Several instruments capable of obtaining optical profiles of multifocal solutions are currently available to researchers ^1–5. The sagittal power of the optical profiles of multifocal contact lenses is increasingly made public, making it possible to infer the optical properties of commercial designs. It is anticipated that in the upcoming years a similar trend will occur with IOLs ⁶. The field of optical design has classically been interested in constructing optical quality metrics able to establish a connection between experimental measurements of Visual Acuity and optical characteristics of the human eye ^7–9. Zernike polynomials are often used in the assessment of the optical properties of the eye since they are a good match for the continuous wavefront of a healthy eye, and from them almost all the optical quality metrics of the eye can be constructed. This set of polynomial functions offers the advantage of accurately representing many of the most common eye refractive errors ^10,11. Second-order Zernike polynomials correctly represent myopia, hyperopia, and astigmatism. Third-order vertical and horizontal coma terms can be used to model misalignments occurring between cornea and crystalline lens, and fourth and sixth-order spherical aberration terms can account for asphericity values typically found in the optical structures of the eye.

Current multifocal refractive lenses use different areas of the pupil to provide clear vision at different distances, each of these areas potentially having different optical properties, such as radius of curvature ^1,12–16. Common terms used presently to refer to multifocal solutions include center-distance (CD) or center-near (CN), corresponding to designs in which the center and the surround areas of the pupil are matched with lens areas that have significantly different radius of curvature (aiming to provide clear vision either for distance or near). When trying to represent a multifocal wavefront for a CD or CN lens produced with Zernike coefficients, the fit is defective (i.e. in multifocal solutions two continuous zones can be discontinuous and present radically different sagittal powers that cannot be represented with a single set of Zernike polynomials). This is one of the major reasons for the lack of a procedure to derive an optical quality metric from optical profiles produced by optical metrology instruments such as NIMO ¹⁷.

During recent years, our group has explored techniques to segment the eye entrance pupil that allow introducing different sets of Zernike polynomials to each of the areas ^12–14. This pupil segmentation allows using multiple sets of Zernike polynomials (one per subarea) with several levels of defocus to represent one entrance pupil. With this technique, discontinuities representing radically different sagittal powers, can be properly represented. Furthermore, a full customization of Zernike coefficients for each Zernike polynomial on each subarea can be performed independently.

Figure 1 represents the basics of the procedure outlined in this paper. An optical profile is obtained from the NIMO, with a density of measurements of 28 points per millimeter across the diameter of the contact lens. This profile is a radial average of a sagittal (axial) power map, and it describes the intersection of light rays originating at any given point in the pupil with the optical axis. One wavefront is created per radial location measured in the optical profile and then a wavefront for a complete pupil is created with a single value of sagittal power, the one present at each radial location. For illustration purposes, in Fig. 1a only 10 radial locations (columns in grey and reddish backgrounds) are represented by their matching wavefront. The annuli that correspond to each of the radial values of sagittal power are cropped from the respective wavefront (lower line of insets in Fig. 1a) and then coupled with the rest of annuli (each representing the corresponding sagittal power for a particular pupil radius interval) to reconstruct a multifocal wavefront out of multiple monofocal wavefronts (Fig. 1b). It is important to note that spatial integration along the pupillary area of the annuli was performed, guaranteeing that spatial overlapping did not occur. Once a multifocal wavefront was created, a traditional Fourier optics procedure ¹⁸ was followed to obtain the Visual Strehl in the frequency domain (Fig. 1c) ¹⁹.

When evaluating multifocal corrections, it is of interest to investigate them for a wide range of working distances, as can be done through this method by adding the particular vergence power corresponding for a given working distance as a base power to the multifocal wavefront. Through-focus evaluation has been performed in this paper for object vergence values ranging from + 2.00 D (working distance located 0.50 m beyond infinity) to -6.00 D (working distance located − 16.6 cm in front of the eye). Distances beyond infinity have been included with dual purpose: 1- to have a full evaluation of the whole through-focus curve (distances beyond infinity could come into play if the eye has a certain degree of residual accommodation), and 2- to give dioptric room to correct for the eye’s spherical equivalent prescription.

The market offers a rich variety of multifocal contact lenses, and it was reasonably expected that separate manufacturers would have produced particular solutions for aiding the presbyopic patient to optimize their optical performance through-focus (known as defocus curve in clinical practice) ^15,16. Also, it was reasonable to expect that presbyopes of varying ages and therefore with distinct levels of residual accommodation would need unique addition prescriptions. Therefore, the study lenses cover four of the main manufacturers of contact lenses worldwide, and two-to-four models each of them manufactures for patients ranging from early presbyopes to patients who are > 65 years old within the same product line. This is the list of lenses in the study: two Bausch & Lomb lenses (Ultra high and Ultra Low), three Alcon lenses (Air Optix Aqua high, Air Optix Aqua medium, Air Optix Aqua low), three Johson and Johson lenses (Oasys high, Oasys Medium, Oasys Low), and four CooperVision lenses (Biofinity 2.50 N, Biofinity 2.00 N, Biofinity 1.50 N & Biofinity 1.00 N).

Potential variability across prescription powers has also been investigated in this study. It is expected that a given model of contact lens will preserve a fixed multifocal profile across all refractive prescription levels. However, a -3.00 D myope needs a lens with different radius of curvature and central thickness than a + 6.00 D hyperope. Therefore, changes in far contact lens prescription power could potentially produce changes in the multifocal profile. Three different far prescriptions (-6.00 D, -3.00 D, and + 1.00 D) have been evaluated with the NIMO instrument.

These 12 multifocal designs were evaluated in combination with the higher-order aberrations (HOAs) of a population of 65 eyes. By including 65 eyes we accounted for inter-subject variability and provided a reasonably solid average contact lens performance that allows for a fair comparison between lenses. Measurements of HOAs also have inherent intra-subject variability but since intra-subject variability is normally one order of magnitude lower than inter-subject, having 65 eyes on this study protects from artifacts arising from this unique source of error.

In summary, this paper presents the mathematical details for the implementation of the required pupil segmentation techniques and classical Fourier optics required to evaluate multifocal corrections. Also, a demonstration of the potential applications of this method is presented, through the evaluation of 12 commercial multifocal contact lenses, with three different far prescriptions (-6.00 D, -3.00 D, and + 1.00 D), and in combination with the HOAs of 65 eyes. All simulations presented include a total of 36 multifocal scenarios with 2340 individual subject cases.

2.1 Measurements of the Optical Profile of Contact Lenses

The NIMO equipment was used to measure the optical profiles of the 12 multifocal lenses in three different base corrections (-6.00 D, -3.00 D, and + 1.00 D as labelled by the manufacturer). Measurements with NIMO are based on the Phase-Shifting Schlieren technique. This technique combines the Schlieren principle with the phase-shifting technique that is used in interferometry by using a Schlieren filter and an adapted set-up. After the application of the phase shift technique, the Schlieren phase was calculated and converted into beam deviation values ¹⁷.

The sagittal powers as a function of the radial coordinate of the 12 multifocal lenses of the study are represented in Fig. 2: two Bausch & Lomb lenses (Ultra high and Ultra Low), three Alcon lenses (Air Optix Aqua high, Air Optix Aqua medium, Air Optix Aqua low), three Johson and Johson lenses (Oasys high, Oasys Medium, Oasys Low), and four CooperVision lenses (Biofinity 2.50 N, Biofinity 2.00 N, Biofinity 1.50 N & Biofinity 1.00 N).

2.2 Eye Population Description

A database with the optical aberrations of 65 healthy eyes provided by Bausch & Lomb was used in this study. This eye database was collected by using standard Shack–Hartmann wavefront sensor technology, this procedure is non-invasive and only requires the light of a super luminescent diode to enter and exit the eye ¹⁰. The procedure to collect the information on these patients adhered to the tenets of the Declaration of Helsinki and to all federal and state laws applicable in New York State. Aberrations included second, third, fourth, and fifth-order aberrations plus sixth-order spherical aberration (also known as secondary spherical aberration). The importance of HOAs in vision diminishes with their order. In this paper, secondary spherical aberration was included in the simulations for its known interactions with other radially-symmetrical aberrations such as defocus and fourth-order spherical aberration ^20,21, which could be of particular interest when elucidating the multifocal properties of ophthalmic devices. Figure 3 contains the description of the aberrations of the 65 eyes of this study ²².

For the simulations in this paper, the eye database had their equivalent sphere corrected. This was achieved in the image plane by shifting the maximum of the VSOTF obtained for each eye with the monofocal lens to 0.00 D vergence (distance correction) and performing the shift found necessary for the monofocal lens also for the 12 multifocal lenses. The amount of shift was particular to each eye but constant for all the lenses on the same eye.

2.3 Pupil Segmentation Techniques

The mathematical expressions that allow segmenting the eye’s pupil radially and angularly have been presented elsewhere ^12,13. In short, a multifocal wavefront with radial symmetry can be generated by creating monofocal wavefronts with different defocus values, cropping the areas of interest for each of the corresponding defocus values, and adding spatially all the cropped annuli. Equations 1 and 2 show the mathematical transcription of the procedure:

$${W}_{CL}=\sum _{j=1}^{n}{w}_{j}\left({D}_{j}\right)* {\sigma }_{j} \left(l\right) \left(1\right)$$

$${\sigma }_{j}\left(l\right)=\left\{\begin{array}{c}\begin{array}{c} 1, \forall \rho \in \left[\frac{Rl}{n},\frac{R\left(l+1\right)}{n}\right]\end{array}\\ 0, otherwise\end{array}\right. \left(2\right)$$

Where ${W}_{CL}$ represents a multifocal contact lens wavefront that is compounded by adding monofocal wavefronts ${w}_{j}$ with a defocus value (${D}_{j}$) obtained from the sagittal powers measured by the NIMO and $j$ denotes the radial coordinate. From each of these monofocal wavefronts, the appropriate radial areas are selected by multiplying ${w}_{j}$ by a masking function ${\sigma }_{j}$ that has a value of one for a particular annulus and takes the value of zero otherwise. Where $R$ represents the total radius measured of the lens and $\rho$ is the radial coordinate.

Where $l$ is linked to the number of radial zones and covers all integer values from (0 to$n-1$). $n$ is the number of radial areas used to represent the optical profile that is equal to the number of points measured by the NIMO in the multifocal lens and always assumed to be greater or equal than 2. Note that, within the same lens a particular $\rho$ and $l$ will have a uniquely defined ${D}_{j}$ obtained from the NIMO measurement. However, when changing from one lens to another the same combination of $\rho$ and $l$ will lead to a different value ${D}_{j}$. Also, to note, $n$ = 1 would represent a monofocal lens and a simple ${W}_{CL}={w}_{1}\left({D}_{1}\right)$.

The multifocal wavefront ${W}_{CL}$ was introduced in the expression of the generalized pupil from which different optical quality metrics can be calculated ¹⁸.

2.4 Metrics Used to Compare the Study Lenses

In this paper we have evaluated 12 multifocal and one monofocal designs by using the metric called Visual Strehl Ratio in the frequency domain as described in Eq. 3 ^7,19.

$$VSOTF=\frac{{\int }_{0}^{2\pi }{\int }_{0}^{60}CSF*\left|Re\left(OTF\right)\right| dw d\theta }{{\int }_{0}^{2\pi }{\int }_{0}^{60}CSF*{OTF}_{DL}dw d\theta } \left(3\right)$$

Where CSF stands for contrast sensitivity function, OTF represents optical transfer function, θ characterizes angles, and $w$ represents frequencies. Integration occurs for all angles (0–2$\pi$) and frequencies between 0 and 60 (VA = 20/10) cycles per degree. The contrast sensitivity function used for the calculations in this paper accounts only for the neural part of visual processing and was obtained from interferometric measurements available in the scientific literature ^24,25. The matrix size used in the simulations was 512 x 512 pixels.

Equation 4 shows the description of the procedure when eye aberrations are included in the simulations.

$${W}_{T}={W}_{CL} +{W}_{Eye} \left(4\right)$$

Where ${W}_{T}$ represents a wavefront superposition in the pupil plane of the multifocal profile ${(W}_{CL})$ of the contact lens and the eye’s aberrations${(W}_{Eye})$. This was achieved by using${ w}_{j(CL+Eye)}$ in Eq. 1 that contained both the sagittal power of the multifocal contact lens measured by the NIMO and the optical aberrations for each of the 65 eyes. From each of these wavefronts, the appropriate annuli are selected by multiplying ${w}_{j(CL+Eye)}$ by a masking function ${\sigma }_{j}$, equivalent to the one defined in Eq. 2. ${\sigma }_{j}$ has a value of 1 for a particular annulus and takes the value of zero otherwise. Where 𝑅 represents the maximum radius of the calculation (the smaller of the pair: pupil radius, contact lens radius measured by the NIMO). Where n, $\rho$,and l represent the same values as in Eq. 2.

The multifocal wavefront ${W}_{T}$loaded with the eye database HOAs was then introduced in the expression of the generalized pupil from which the VSOTF shown in Eq. 3 was calculated. The combination of eye plus the multifocal design has been considered as a single VSOTF because the image created by the multifocal correction keeps the phase information intact (i.e. no intermediate image is created on a screen) when interacting with the optics of the eye.

To establish comparisons between lenses, three metrics have been defined: area under the through-focus VSOTF, range above threshold, and peak performance. The definition of these three metrics can be seen in equations 5, 6, 7, and 8.

$$Area under {VSOTF}_{TF}=\sum _{-1}^{5}VSOTF \left(5\right)$$

$$Range above threshold= \sum _{-1}^{5}{x}_{d} \left(6\right)$$

$${x}_{d}\left(d\right)=\left\{\begin{array}{c}\begin{array}{c} \frac{6}{193} D, \forall VSOTF\left(d\right)>0.12\end{array}\\ 0, \forall VSOTF\left(d\right)<0.12\end{array}\right. \left(7\right)$$

Where 6/193 diopters for ${x}_{d}$ comes from the number of samples used to describe the through-focus curves (6.00 D of interval and 193 sample points). A VSOTF threshold of 0.12 predicts a level of logMAR visual acuity of 0.20 ²⁶, and it has been used abundantly in the literature ^26–28. Analogous calculations for other logMAR levels could be performed by varying the level of the threshold of the VSOTF. Values under 0.12 are not considered for the interval of acceptable vision but are still of some utility (corresponding logMAR values between 0.20 and 0.40, 6/9 to 6/15 or 20/30 to 20/50 Snellen)²⁶ and therefore counted on the area under the curve metric.

$$Peak Performance=\text{max} \left(Throughfocus VSOTF\right) \left(8\right)$$

In summary, through-focus (-1.00 D to + 5.00 D) monocular VSOTFs (Eq. 3) were calculated for the 12 multifocal commercial designs described in Fig. 2 in combination with the 65 eyes described in Fig. 3 and with one monofocal design that serves as a reference case.

As a result, a total of 2340 monocular through-focus VSOTFs curves were calculated. Figure 4 shows the average performance through-focus of each of the lenses, with the shaded areas representing ±½ of the standard deviation for each lens across the 65 eyes. We want to focus our interest on the multifocal nature of the profile more than in the refractive correction of the lens, therefore the nominal − 3.00 D base correction for distance has been removed previously to performing the calculations.

Upon direct visual inspection of Fig. 2 the reader can observe that no systematic remarkable differences are observed between distance prescription powers, with the potential exception of the − 6.00 D prescription offering more sagittal power than the − 3.00 D (0 D to 0.50 D), and the + 1.00 D prescription offering less sagittal power than the − 3.00 D (0.05 D to 0.40 D). This difference roughly translates into a shift of the overall VSOTF curves from one distance prescription to another. Therefore, in Fig. 4 only the through-focus VSOTFs results obtained for the − 3.00 D distance correction are presented. Data obtained for the − 6.00 D and + 1.00 D prescriptions will be made available to sensible direct requests to the author of this study.

Figure 5. Average performance of 13 lenses (one monofocal and 12 multifocal) for our population of 65 eyes. This graph shows area under the VSOTF through-focus, range above threshold, and peak performance as defined by equations 5, 6, and 8 respectively. All of them normalized by the corresponding value obtained for the monofocal lens. Dashed light blue line represents 100%; values larger than 100% indicate a larger value than the one obtained with the monofocal lens. Shaded blue indicates monofocal values, alternating gray and white shading is intended to help with visualization.

Nine out of the 12 multifocal designs offer a larger range above threshold than the monofocal design: Ultra High (129%), Ultra Low (112%), Air Optix Aqua High (126%), Air Optix Aqua Medium (122%), Air Optix Aqua Low (110%), Biofinity 2.00 N (104%), Biofinity 1.50 N (113%), Biofinity 1.00 N (106%), and Oasys High (109%). Three out of the 12 designs offer an area under the curve within a 2% of the one offered by a monofocal lens: Ultra High (102%), Ultra Low (101%), Biofinity 1.00 N (99.4%). All lenses offer a peak performance that is lower than the one provided by a monofocal lens, with the three best performers in this aspect being Air Optix Aqua Low (78%), Ultra Low (75%), and Biofinity 1.00 N (71%).

The optical performance as a function of the area under the VSOTF and the range above threshold of each of the 65 eyes of our study population averaged across the 12 multifocal lenses in reference to the eye’s pupil size (gray scale bar) is shown in Fig. 6.

4.1 Lens Characteristics

This paper clearly shows the existing tradeoff between peak performance and range above threshold through-focus provided by a particular lens. Early presbyopes can use residual accommodation to move the curves provided by multifocal corrections, therefore designs targeted to early presbyopes (Ultra Low, Air Optix Aqua low, Oasys Low, Biofinity 1.50 N & Biofinity 1.00 N), exhibit a higher peak performance that can be used in combination with residual accommodation to overcome the generally shorter interval through-focus above threshold that these lenses offer. On the contrary, lenses that are created for older presbyopes, whom have no residual accommodation left, offer larger intervals of optical quality above threshold, and sacrifice peak performance (Ultra high, Air Optix Aqua high, Oasys high, and Biofinity 2.50 N). Also, an intermediate version exists for presbyopes with some but not a large amount of residual accommodation (Air Optix Aqua medium, Oasys Medium, Biofinity 1.50 N & Biofinity 2.00 N).

The peak performance of all lenses in the study, with the exceptions of the Oasys Low and the Oasys Medium models, occurs at a distance from the patient between 80 cm and 2 meters (vergences of 1.25 D and 0.50 D respectively). Placing peak performance closer than infinity allows for better use of the interval of good quality of vision occurring on both sides of the work distance at which peak performance occurs (optimization of the concept of hyperfocal distance).

Of special interest are the Biofinity “N” presented in this study: their through-focus curves explain why they are to be used in combination with the “D” versions. The “N” lenses seem to disregard the optical quality obtained for distance correction, and focused on providing acceptable quality of vision at near; they are CN designs. Biofinity “D” lenses designed to be used in the dominant eye carry the distance correction on the center of the lens (CD); this is the opposite of the Biofinity lenses for the non-dominant eye “N” presented on this study. “D” lenses will provide the optical quality necessary for far vision. Together the “N” (CN) and “D” (CD) versions of the Biofinity lenses provide a modified monovision correction instead of a truly multifocal correction where both eyes will experience the same optical correction. Because of this monovision feature, and not having presented the characteristics of the “D” lenses in this manuscript, a direct comparison is difficult to be made with the rest of the study lenses and therefore the following discussion does not include them.

Oasys lenses show the importance of the effective aperture effect ²⁹. This effect basically predicts that the higher the level of change on the sagittal power across the radial coordinates in the entrance pupil, the lower the overall optical quality provided. Peak performance for Oasys Low lenses is 21% lower than the one offered by Ultra-Low and 23% lower than the one offered by Air Optix Aqua Low lenses. Interval through-focus above threshold for Oasys Low lenses is 18% lower than the one offered by Ultra-Low and 17% lower than the one offered by Air Optix Aqua Low lenses. Area under the curve for Oasys Low lenses is 19% lower than the one offered by Ultra-Low and 11% lower than the one offered by Air Optix Aqua Low lenses. Although the sagittal power shown in the third row of Fig. 2 shows two foci for the Oasys Lenses, there are three main reasons why their optical performance will not clearly display the expected two foci of clear vision. 1- The higher-order aberrations of the 65 eyes smooth the VSOTF of the lens alone, 2- As can be seen in Fig. 3, pupil sizes vary greatly between the 65 eyes, with a majority of our subject database having pupils under 4.25 mm, therefore the second increase of power in the profile is often not used, and 3- multiple transition areas between near and far powers contribute to smooth out the potential two foci design.

A clear trend appears between a smaller pupil size and a better performance through-focus on average across the 12 designs of the study. When dealing with a monofocal lens, smaller pupil sizes (up to 2 mm in diameter) reduce the amount of aberrations present in the eye and increase overall visual performance through-focus, for as long as diffraction does not become a dominant factor ³⁰. Figure 6 confirms that this effect is still present when multifocal lenses are used. This reinforces the predominant clinical practice whereby presbyopes are instructed by clinicians to maintain high levels of illumination when performing near tasks, and with natural evolution where senile myosis contributes to extend the depth of focus of elderly patients ^31,32. This phenomenon has allowed for the eruption of multiple technologies and pharmacological treatments in recent years that limit pupil size in order to increase a patient’s visual performance through-focus ^33–37.

4.2 Utility of the Proposed Method

There have been attempts in the past to evaluate depth of focus of healthy subjects from wavefront measurements ³⁸. However, these early attempts did not include multifocal corrections. As explained above, wavefront reconstruction with a multifocal correction would not be properly fitted with Zernike polynomials unless the process presented in this paper is used. The method we are describing arises from some of our work published during the last decade ^27,39,40. Also, the materials presented in this article have been defended in conferences by our group during the last six years ^41–43. However, this is the first time this method of converting sagittal power into the corresponding annuli of a multifocal wavefront (created with standard Zernike polynomials) is fully published in a paper format. The method allows the systematical evaluation of the optical performance of multifocal corrections through simulations, which do not require the implementation of optical benches. The central realization that makes this work valuable is the fact that sagittal power obtained with metrology devices can be represented locally with the corresponding wavefront obtained from Zernike polynomials. This was never presented in any of our previous papers ^27,39,40, and was not present in classical work ^7,8,44.

An approach to simulate multifocal corrections called temporal multiplexing has been recently introduced. Temporal multiplexing consists in presenting different refractive powers to the eye each of them ranging from a tenth of a microsecond to several microseconds ⁴⁵. Due to the short exposure time, each of the individual defocus values blends with others, the final integrated signal containing multiple defocus values. Temporal multiplexing does represent accurately the percentage of light carried by each refractive power. However, it uses the whole pupil for each refractive power represented and it is not capable of segmenting the pupil in different areas or creating any kind of asymmetries (neither radial nor angular). On the other hand, the method presented in this paper has the ability to represent any local distributions of power around the entrance pupil.

Other approaches to evaluate multifocal corrections were presented by Alarcon et al. ⁴⁶, Vega et al.⁴⁷, Armengol et al. ⁴⁸, Azor et al. ⁴⁹, and Torres-Sepulveda et al. ⁵⁰. In their works, the optical quality of intraocular lenses was evaluated by creating an optical system with an eye model. All these procedures do allow for testing multifocal designs but need a physical lens to be employed and experimental measurements to be taken in an optical bench. Furthermore, a recent study has shown an instrument capable of studying the optical capabilities of theoretical designs for spectacle lenses by using a spatial light modulator ⁵¹. Although this type of instrumentation can reproduce multifocal patterns, it requires an initial investment, and will have difficulties when evaluating smaller optical areas, such as the ones needed to represent contact lenses. On the other hand, our method represents the optics of contact lenses, and intraocular lenses already present in the market and those that do not yet exist. To reproduce lenses that do not exist yet, we could substitute the sagittal powers of the 12 lenses used as input on this article by any other combination of sagittal powers.

Bakaraju et al.⁵² have presented a method where by coupling a custom-made script in MATLAB with Zemax, the optical profiles obtained in a NIMO instrument can be evaluated. This is done through creating a super conic surface in the front surface of the correcting contact lens in their Zemax model. Their method requires the utilization of Zemax to be completed, which introduces extra steps and reduces control over the process. It is not clear from their publication if their method allows for inclusion of any given set of aberrations coming from a particular subject or whether they could have as output any desired optical quality metric ⁷.

Although radial asymmetries are not present on the set of multifocal contact lenses reported in this paper, in the simplest case of radial asymmetries, multifocal lenses with a correction for astigmatism could be studied. The most straightforward way to deal with the presence of astigmatism is to model the individual wavefronts of Fig. 1b, with a certain level of astigmatism already imbued in them that could easily be obtained from NIMO readings. Ultimately, the limits for representing asymmetries in a lens with this method are: 1- the precision with which the NIMO (or any other) equipment can measure sagittal powers, and 2 - the real-life utility of these constructions. Segmentation of the pupil could be done for any kind of sub-area considered of interest, by plugging a certain level of sagittal power, creating a wavefront representing it, and reconstructing the multifocal wavefront pixel by pixel if required. In short, the method proposed in this article is more cost effective, more flexible, and has more potential applications than the existing methods in the literature.

The algorithms presented provide a bridge that allows the evaluation of a multifocal design with Fourier optics and Zernike-based wavefronts. The work presented in this paper can be extended to problems related to myopia progression control where the central and peripheral optical quality provided by multifocal contact lenses specifically designed for myopia control (MYLO and MiSight) needs to be evaluated. It can also be extended to small aperture solutions such as the IC-8 Intraocular Lens, Kamra Inlay, or pharmaceutical approaches that aim to reduce pupil size.

Acknowledgements

This work was funded by Bausch and Lomb. I thank Ian Cox from Bausch and Lomb for providing the measurements for the multifocal contact lenses and sharing the database of aberrations for the 65 eyes on this study.

Disclosures

This work was funded by Bausch and Lomb. The company provided the database with 65 eyes and the measurements performed by with the NIMO instrument. The author did not receive any further directions on how to perform the study. The author has no current commercial interest in the results presented in this study.

Data availability

The data and algorithms used in this study are not publicly available. The author will consider on a case-by-case basis whether to share the information of this paper to reasonable requests. The author’s personal email is: [email protected].

Plainis, S., Atchison, D. A. & Charman, W. N. Power profiles of multifocal contact lenses and their interpretation. Optom Vis Sci 90, 1066-1077, doi:10.1097/opx.0000000000000030 (2013).
Kim, E., Bakaraju, R. C. & Ehrmann, K. Reliability of power profiles measured on NIMO TR1504 (Lambda-X) and effects of lens decentration for single vision, bifocal and multifocal contact lenses. Journal of Optometry 9, 126-136, doi:https://doi.org/10.1016/j.optom.2015.10.005 (2016).
Montés-Micó, R., Madrid-Costa, D., Domínguez-Vicent, A., Belda-Salmerón, L. & Ferrer-Blasco, T. In vitro power profiles of multifocal simultaneous vision contact lenses. Contact Lens and Anterior Eye 37, 162-167, doi:http://dx.doi.org/10.1016/j.clae.2013.09.014 (2014).
Ruiz-Alcocer, J. Analysis of the power profile of a new soft contact lens for myopia progression. Journal of optometry 10, 266-268, doi:10.1016/j.optom.2016.08.003 (2017).
Domínguez-Vicent, A., Marín-Franch, I., Esteve-Taboada, J. J., Madrid-Costa, D. & Montés-Micó, R. Repeatability of in vitro power profile measurements for multifocal contact lenses. Contact Lens and Anterior Eye 38, 168-172, doi:https://doi.org/10.1016/j.clae.2015.01.013 (2015).
Madrid-Costa, D. et al. Beam-Shaping Extended Depth of Focus Intraocular Lens: Optical Assessment With Corneas of Increasing Spherical Aberration. Journal of Refractive Surgery 39, 95-102, doi:doi:10.3928/1081597X-20221215-02 (2023).
Marsack, J. D., Thibos, L. N. & Applegate, R. A. Metrics of optical quality derived from wave aberrations predict visual performance. Journal of Vision 4, 322-328 (2004).
Applegate, R. A., Marsack, J. D. & Thibos, L. N. Metrics of retinal image quality predict visual performance in eyes with 20/17 or better visual acuity. Optometry and Vision Science 83, 635-640 (2006).
Applegate, R. A., Donnelly, W. J., Marsack, J. D., Koenig, D. E. & Pesudovs, K. Three-dimensional relationship between high-order root-mean-square wavefront error, pupil diameter, and aging. Journal of the Optical Society of America a-Optics Image Science and Vision 24, 578-587 (2007).
Liang, J. & Williams, D. R. Aberrations and retinal image quality of the normal human eye. Journal of the Optical Society of America A 14, 2873-2883 (1997).
Salmon, T., Thibos, L. & Bradley, A. Comparison of the eye´s wave-front aberration measured psychophysically and with the Shack-Hartmann wave-front sensor. J. Opt. Soc. Am. A 15, 2457-2465 (1998).
de Gracia, P., Dorronsoro, C. & Marcos, S. Multiple zone multifocal phase designs. Opt Lett 38, 3526-3529, doi:10.1364/ol.38.003526 (2013).
de Gracia, P. & Hartwig, A. Optimal orientation for angularly segmented multifocal corrections. Ophthalmic and Physiological Optics 37, 610-623, doi:https://doi.org/10.1111/opo.12402 (2017).
de Gracia, P. Optical properties of monovision corrections using multifocal designs for near vision. Journal of Cataract & Refractive Surgery 42, 1501-1510, doi:https://doi.org/10.1016/j.jcrs.2016.10.002 (2016).
Plainis, S., Ntzilepis, G., Atchison, D. A. & Charman, W. N. Through-focus performance with multifocal contact lenses: effect of binocularity, pupil diameter and inherent ocular aberrations. Ophthalmic and Physiological Optics 33, 42-50, doi:10.1111/opo.12004 (2013).
Kim, E., Bakaraju, R. C. & Ehrmann, K. Power Profiles of Commercial Multifocal Soft Contact Lenses. Optom Vis Sci 94, 183-196, doi:10.1097/opx.0000000000000998 (2017).
Joannes, L., Heraud, M., Ligot, R., Saoul, B. & Dupont, O. NIMO: a new tool for asphere and free-form optics measurement. Vol. 6341 SPE (SPIE, 2006).
Atchison, D. A. & Smith, G. Optics of the Human Eye. 1 edn, 288 (Butterworth-Heinemann, 2000).
Iskander, D. R. Computational Aspects of the Visual Strehl Ratio. Optometry & Vision Science 83, 57-59 (2006).
Benard, Y., Lopez-Gil, N. & Legras, R. Optimizing the subjective depth-of-focus with combinations of fourth- and sixth-order spherical aberration. Vision Research 51, 2471-2477, doi:10.1016/j.visres.2011.10.003 (2011).
Yi, F., Robert Iskander, D. & Collins, M. Depth of focus and visual acuity with primary and secondary spherical aberration. Vision Research 51, 1648-1658, doi:10.1016/j.visres.2011.05.006 (2011).
Spitzer, M., Wildenhain, J., Rappsilber, J. & Tyers, M. BoxPlotR: a web tool for generation of box plots. Nat Meth 11, 121-122, doi:10.1038/nmeth.2811 (2014).
Thibos, L. N., Applegate, R. A., Schwiegerling, J. T., Webb, R. H. & Members, V. S. T. Standards for reporting the optical aberrations of eyes. Vision Science and its Applications, OSA Trends in Optics & Photonics 35, 110-130 (2000).
Williams, D. R. Aliasing in human foveal vision. Vision Research 25, 195-205 (1985).
Williams, D. R., Brainard, D. H., McMahon, M. J. & Navarro, R. Double-pass and interferometric measures of the optical quality of the eye. J. Opt. Soc. Am. A 11, 3123-3135 (1994).
Cheng, X., Bradley, A. & Thibos, L. N. Predicting subjective judgment of best focus with objective image quality metrics. Journal of Vision 4, 310-321 (2004).
De Gracia, P., Dorronsoro, C. & Marcos, S. Multiple zone multifocal phase designs. Opt. Lett. 38, 3526-3529, doi:10.1364/OL.38.003526 (2013).
Radhakrishnan, A., Dorronsoro, C. & Marcos, S. Differences in visual quality with orientation of a rotationally asymmetric bifocal intraocular lens design. Journal of Cataract & Refractive Surgery 42, 1276-1287, doi:http://dx.doi.org/10.1016/j.jcrs.2016.06.034 (2016).
Xu, R., Bradley, A. & Thibos, L. N. Impact of primary spherical aberration, spatial frequency and Stiles Crawford apodization on wavefront determined refractive error: a computational study. Ophthalmic and Physiological Optics 33, 444-455, doi:10.1111/opo.12072 (2013).
Xu, R., Thibos, L. & Bradley, A. Effect of Target Luminance on Optimum Pupil Diameter for Presbyopic Eyes. Optometry and Vision Science 93, 1409-1419, doi:10.1097/opx.0000000000000963 (2016).
Sloane, M. E., Owsley, C. & Alvarez, S. L. Aging, senile miosis and spatial contrast sensitivity at low luminance. Vision Research 28, 1235-1246, doi:https://doi.org/10.1016/0042-6989(88)90039-9 (1988).
Suzuki, T.-a., Okajima, K. & Funai, T. Optical simulation of reduced retinal illuminance caused by senile miosis. Optical Review 19, 174-181, doi:10.1007/s10043-012-0029-3 (2012).
Grabner, G., Ang, R. E. & Vilupuru, S. The Small-Aperture IC-8 Intraocular Lens: A New Concept for Added Depth of Focus in Cataract Patients. American Journal of Ophthalmology 160, 1176-1184.e1171, doi:https://doi.org/10.1016/j.ajo.2015.08.017 (2015).
Dexl, A. K. et al. Reading performance after implantation of a small-aperture corneal inlay for the surgical correction of presbyopia: Two-year follow-up. Journal of cataract and refractive surgery, doi:10.1016/j.jcrs.2010.10.044 (2011).
Naroo, S. A. & Bilkhu, P. S. Clinical utility of the KAMRA corneal inlay. Clinical Ophthalmology 10, 913-919, doi:10.2147/OPTH.S89132 (2016).
Fenner, B. J., Moriyama, A. S. & Mehta, J. S. Inlays and the cornea. Experimental Eye Research 205, 108474, doi:https://doi.org/10.1016/j.exer.2021.108474 (2021).
Grzybowski, A., Markeviciute, A. & Zemaitiene, R. A Review of Pharmacological Presbyopia Treatment. Asia Pac J Ophthalmol (Phila) 9, 226-233, doi:10.1097/apo.0000000000000297 (2020).
Yi, F., Iskander, D. R. & Collins, M. J. Estimation of the depth of focus from wavefront measurements. Journal of Vision 10, 3-3, doi:10.1167/10.4.3 (2010).
De Gracia, P. & Hartwig, A. Optimal orientation for angularly segmented multifocal corrections. Ophthalmic & physiological optics : the journal of the British College of Ophthalmic Opticians (Optometrists) 37, 610-623, doi:10.1111/opo.12402 (2017).
De Gracia, P. Optical properties of monovision corrections using multifocal designs for near vision. Journal of cataract and refractive surgery 42, 1501-1510, doi:10.1016/j.jcrs.2016.10.002 (2016).
De Gracia, P. in American Academy of Optometry Meeting (American Academy of Optometry, Chicago, 2017).
De Gracia, P., Des Rosiers, C. & Whitescarver, T. D. Segmentation Tools for the Fourier Evaluation of Commercially Avalilable Multifocal Contact Lenses. Investigative ophthalmology & visual science 60, 6335-6335 (2019).
De Gracia, P., Purrier, N., Des Rosiers, C., Marcell, A. & Pashak, G. Study of the optical properties of two brands of commercial multifocal contact lenses (steep vs smooth addition increase) in combination with the higher order aberrations of 65 eyes. Investigative ophthalmology & visual science 59, 5809-5809 (2018).
Applegate, R. A., Marsack, J. D., Ramos, R. & Sarver, E. J. Interaction between aberrations to improve or reduce visual performance. Journal of cataract and refractive surgery 29, 1487-1495 (2003).
Vinas, M. et al. Optical and Visual Quality With Physical and Visually Simulated Presbyopic Multifocal Contact Lenses. Translational Vision Science & Technology 9, 20-20, doi:10.1167/tvst.9.10.20 (2020).
Alarcon, A. et al. Preclinical metrics to predict through-focus visual acuity for pseudophakic patients. Biomed. Opt. Express 7, 1877-1888, doi:10.1364/BOE.7.001877 (2016).
Vega, F. et al. Visual acuity of pseudophakic patients predicted from in-vitro measurements of intraocular lenses with different design. Biomed Opt Express 9, 4893-4906, doi:10.1364/boe.9.004893 (2018).
Armengol, J., Garzón, N., Vega, F., Altemir, I. & Millán, M. S. Equivalence of two optical quality metrics to predict the visual acuity of multifocal pseudophakic patients. Biomed. Opt. Express 11, 2818-2829, doi:10.1364/BOE.388531 (2020).
Azor, J. A., Vega, F., Armengol, J. & Millan, M. S. Optical Assessment and Expected Visual Quality of Four Extended Range of Vision Intraocular Lenses. J Refract Surg 38, 688-697, doi:10.3928/1081597x-20220926-01 (2022).
Torres-Sepúlveda, W., Mira-Agudelo, A., Barrera-Ramírez, J. F. & Kolodziejczyk, A. Objective method for visual performance prediction. J. Opt. Soc. Am. A 40, C138-C149, doi:10.1364/JOSAA.478022 (2023).
Arias, A., Ohlendorf, A., Artal, P. & Wahl, S. In-depth optical characterization of spectacle lenses for myopia progression management. Optica 10, 594-603, doi:10.1364/OPTICA.486389 (2023).
Bakaraju, R. C., Ehrmann, K. & Ho, A. Extended depth of focus contact lenses vs. two commercial multifocals: Part 1. Optical performance evaluation via computed through-focus retinal image quality metrics. Journal of Optometry 11, 10-20, doi:https://doi.org/10.1016/j.optom.2017.04.003 (2018).

Competing interest reported. This study was sponsored by Bausch & Lomb

Download PDF

Journal Publication

published 19 Dec, 2023

Read the published version in Scientific Reports →

Editorial decision: Revision requested
10 Nov, 2023
Reviews received at journal
06 Nov, 2023
Reviewers agreed at journal
21 Oct, 2023
Reviewers invited by journal
20 Sep, 2023
Editor assigned by journal
20 Sep, 2023
Editor invited by journal
08 Sep, 2023
Submission checks completed at journal
08 Sep, 2023
First submitted to journal
06 Sep, 2023

You are reading this latest preprint version

Fourier Tools for the Evaluation of Refractive Multifocal Designs

Status:

Journal Publication

Version 1

Abstract

Figures

1. Introduction

2. Methods

2.1 Measurements of the Optical Profile of Contact Lenses

2.2 Eye Population Description

2.3 Pupil Segmentation Techniques

2.4 Metrics Used to Compare the Study Lenses

3. Results

4. Discussion

4.1 Lens Characteristics

4.2 Utility of the Proposed Method

5. Conclusion and future work

Declarations

References

Additional Declarations

Status:

Journal Publication

Version 1