A non-stationary model for simulating the dynamics of ocular aberrations

doi:10.1364/OE.18.021386

« Show journal navigation

A non-stationary model for simulating the dynamics of ocular aberrations

C. Leahy and C. Dainty »View Author Affiliations

Optics Express, Vol. 18, Issue 20, pp. 21386-21396 (2010)
http://dx.doi.org/10.1364/OE.18.021386

View Full Text Article

Acrobat PDF (1445 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Methods
3. Results
4. Discussion
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (32)
Cited By
Figures (7)
Metrics

Abstract

The time-evolution of ocular aberrations has been the subject of many studies, but so far there has been little discussion involving the modelling of the underlying temporal statistics. This paper presents a non-stationary modelling approach based on a coloured-noise generator, which can be applied to ocular aberration dynamics. The model parameters are computed from measured ocular aberration data. A custom-built aberrometer based on a Shack-Hartmann sensor was used for measurement. We present simulations based on our modelling approach, and validate them through comparison to real data. This work could be useful in areas such as the testing of ophthalmic devices and the development of improved algorithms for laser refractive surgery.

1. Introduction

Signal models of physiological processes are commonplace [1

E. N. Bruce, Biomedical Signal Processing and Signal Modeling (Wiley Series in Telecommunications and Signal Processing, 2001).

]; however there has been very little research regarding models of dynamic ocular aberrations other than defocus. Roberts [2

C. Roberts, “Future challenges to aberration-free ablative procedures,” J. Refract. Surg. 16, 623–629 (2000).

] asserted that measuring and understanding the dynamic performance of the eye is critical to the future development of advanced customised refractive surgeries, which is in itself important for generating an overall improvement of refractive surgery results. Therefore, any improvements in the measurement and modelling of dynamic aberrations could be beneficial in areas such as the development of customised ablation algorithms. Such models could also be of interest in the study of the relationship between dynamics and visual performance [3

A. Mira-Agudelo, L. Lundström, and P. Artal, “Temporal dynamics of ocular aberrations: Monocular vs binocular vision,” Ophthal. Physiol. Opt. 29, 256–263 (2009). [CrossRef]

], the development of customised contact lenses [4

R. Navarro, E. Moreno-Barriuso, S. Bará, and T. Mancebo, “Phase plates for wave-aberration compensation in the human eye,” Opt. Lett. 25, 236–238 (2000). [CrossRef]

], and in the design of ocular adaptive optics systems with real-time correction [5

S. O. Galetskiy, T. Yu. Cherezova, and A. V. Kudryashov, “Adaptive optics in ophthalmology: human eye wavefront generator,” Proc. SPIE 6849, 68490918 (2008).

Many dynamic properties of ocular aberrations have been established by previous authors [6

D. R. Iskander, M. Collins, M. Morelande, and M. Zhu, “Analyzing the dynamic wavefront aberrations in the human eye,” IEEE Trans. Biomed. Eng 51, 1969–1980 (2004). [CrossRef] [PubMed]

, 7

H. Hofer, P. Artal, and D. R. Williams, “Dynamics of the eyes wave aberration,” J. Opt. Soc. Am. A. 18, 497–506 (2001). [CrossRef]

]. In particular, the time evolutions of some Zernike aberration modes have been shown to be correlated with each other [8

K. Hampson, E. Mallen, and J. C. Dainty, “Coherence function analysis of the higher-order aberrations of the human eye,” Opt. Lett. 31, 184–186 (2006). [CrossRef] [PubMed]

], as well as with the cardiopulmonary system [9

K. Hampson, I. Munro, C. Paterson, and J. C. Dainty, “Weak correlation between the aberration dynamics of the human eye and the cardiopulmonary system,” J. Opt. Soc. Am. A. 22, 1241–1250 (2005). [CrossRef]

, 10

M. Zhu, M. J. Collins, and D. R. Iskander, “Microfluctuations of wavefront aberrations of the eye,” Opthal. Physiol. Opt. 24, 562–571 (2004). [CrossRef]

]. Previous work has shown that the power spectral density (PSD) of aberration signals tends to follow a power-law relationship for frequencies below approximately 10 Hz, when the subject is not actively accommodating [3

A. Mira-Agudelo, L. Lundström, and P. Artal, “Temporal dynamics of ocular aberrations: Monocular vs binocular vision,” Ophthal. Physiol. Opt. 29, 256–263 (2009). [CrossRef]

H. Hofer, P. Artal, and D. R. Williams, “Dynamics of the eyes wave aberration,” J. Opt. Soc. Am. A. 18, 497–506 (2001). [CrossRef]

]. The dynamics of Zernike modes have also been shown to exhibit features such as non-stationarity, harmonic distortion, and harmonic amplitude modulation [6

D. R. Iskander, M. Collins, M. Morelande, and M. Zhu, “Analyzing the dynamic wavefront aberrations in the human eye,” IEEE Trans. Biomed. Eng 51, 1969–1980 (2004). [CrossRef] [PubMed]

, 11

D. R. Iskander, M. R. Morelande, and M. J. Collins, “Estimating the dynamics of aberration components in the human eye,” In IEEE Signal Process. Workshop Statist. pages241–244 (2001).

Iskander et al. [11

D. R. Iskander, M. R. Morelande, and M. J. Collins, “Estimating the dynamics of aberration components in the human eye,” In IEEE Signal Process. Workshop Statist. pages241–244 (2001).

] proposed a methodology for modelling the dynamics of higher-order ocular aberrations using a parametric, amplitude-modulated, frequency-modulated (AM-FM) approach. This study focused in particular on the Zernike spherical aberration and coma coefficients, however measurements of just 5 s in duration (without taking into account the effects of subjects’ blinking) were considered. Other work has demonstrated a potential application of a model of dynamic aberrations, by employing a deformable mirror to generate aberrations in real time [5

S. O. Galetskiy, T. Yu. Cherezova, and A. V. Kudryashov, “Adaptive optics in ophthalmology: human eye wavefront generator,” Proc. SPIE 6849, 68490918 (2008).

In this work, we present a methodology for constructing power-law models of the temporal dynamics of ocular aberrations. The modelling approach is based on a numerical method for coloured-noise generation, developed by Billah and Shinozuka [12

K. Billah and M. Shinozuka, “Numerical method for colored noise generation and its application to a bistable system,” Phys. Rev. A 42, 7492–7495 (1990). [CrossRef] [PubMed]

]. We describe the connection between power-law relationships in the frequency domain and temporal stationarity, and how this may be exploited in the modelling of dynamic aberrations. We attempt to validate the performance of our model by comparing real and simulated data.

2. Methods

2.1. Data collection

All data was collected with a custom-built aberrometer based on a Shack-Hartmann wavefront sensor. The subject was a 26 year-old male. For each measurement, data was collected from the subject’s dominant eye, whereas the non-dominant eye was covered with an eye patch. Eye dominance was assessed through an eye examination by a qualified optometrist prior to our experiments. Data was collected over 40 s, with an active pupil diameter of 5.4 mm. The sampling frequency was 100 Hz, giving a total of 4,000 samples per trial. The defocus term had the greatest magnitude, but there were significant contributions from astigmatism, coma, spherical aberration, and other higher-order terms. The system noise was white, with an RMS value of approximately 3 nm, as evaluated with an artificial eye (consisting of an 18 mm singlet coupled with a screen in the focal plane).

Analysis of dynamic ocular wavefront sensor measurements often takes the form of analysis of the Zernike coefficients as time series [3

A. Mira-Agudelo, L. Lundström, and P. Artal, “Temporal dynamics of ocular aberrations: Monocular vs binocular vision,” Ophthal. Physiol. Opt. 29, 256–263 (2009). [CrossRef]

, 6

D. R. Iskander, M. Collins, M. Morelande, and M. Zhu, “Analyzing the dynamic wavefront aberrations in the human eye,” IEEE Trans. Biomed. Eng 51, 1969–1980 (2004). [CrossRef] [PubMed]

–9

, 13

N. Davies, L. Diaz-Santana, and D. Lara-Sucedo, “Repeatability of ocular wavefront measurement,” Optom. Vis. Sci. 80, 142–150 (2003). [CrossRef] [PubMed]

–17

K. Y. Li and G. Yoon, “Changes in aberrations and retinal image quality due to tear film dynamics,” Opt. Express 14, 12552–12559 (2006). [CrossRef] [PubMed]

]. The temporal reconstruction of ocular wavefronts using Zernike polynomials can be represented as given in Eq. (1) [6

D. R. Iskander, M. Collins, M. Morelande, and M. Zhu, “Analyzing the dynamic wavefront aberrations in the human eye,” IEEE Trans. Biomed. Eng 51, 1969–1980 (2004). [CrossRef] [PubMed]

W (ρ, θ; t) = \sum_{n = 1}^{N} \sum_{m = 0}^{n} c_{n}^{m} (t) Z_{n}^{m} (ρ, θ) + ɛ (ρ, θ; t)

(1)

where W(ρ, θ;t) is the dynamic wave aberration function, m and n are the radial and azimuthal orders respectively,

Z_{n}^{m} (ρ, θ)

is the Zernike polynomial expansion term,

c_{n}^{m} (t)

gives the corresponding Zernike coefficients, and ε(ρ, θ;t) represents the modelling error. For our studies of aberration signals, we used an expansion of Zernike terms up to and including the 8 ^th radial order.

Measurements of dynamic ocular aberrations of more than a few seconds in duration are inevitably degraded by the subject’s blinking. This leads to discontinuities and/or missing values in the measured signal [6

D. R. Iskander, M. Collins, M. Morelande, and M. Zhu, “Analyzing the dynamic wavefront aberrations in the human eye,” IEEE Trans. Biomed. Eng 51, 1969–1980 (2004). [CrossRef] [PubMed]

, 28

C. Leahy, C. Leroux, C. Dainty, and L. Diaz-Santana, “Temporal dynamics and statistical characteristics of the microfluctuations of accommodation: Dependence on the mean accommodative effort,” Opt. Express 18, 2668–2681 (2010). [CrossRef] [PubMed]

]. Though the blink interval is typically on the order of 250 ms [6

D. R. Iskander, M. Collins, M. Morelande, and M. Zhu, “Analyzing the dynamic wavefront aberrations in the human eye,” IEEE Trans. Biomed. Eng 51, 1969–1980 (2004). [CrossRef] [PubMed]

], the transient components of the measured aberration signal associated with the blink can last for up to a second or more [28

]. For all measurements, we systematically removed any part of the signal that was related to the blink interval: tear film break-up and build-up, eye movement, and increased noise in the reconstruction of the Zernike modes of the wavefront due to a vignetted pupil. To incorporate these phenomena would significantly increase the complexity of any modelling approach. Therefore, as an approximation, we did not include them as part of the modelling objectives in this study.

2.2. Modelling approach

The modelling method is based on the assumption of a linear or piecewise linear power-law structure of the power spectral density (PSD) of the aberration signals. Synthesis of signals is based on a modification of an algorithm for the generation of power-law processes (or “coloured noise”) [12

K. Billah and M. Shinozuka, “Numerical method for colored noise generation and its application to a bistable system,” Phys. Rev. A 42, 7492–7495 (1990). [CrossRef] [PubMed]

]. This algorithm has also been utilised as the basis for modelling approaches used in the field of EEG research [18

L. Rankine, N. Stevenson, M. Mesbah, and B. Boashash, “A nonstationary model of newborn EEG,” IEEE Trans. Biomed. Eng 54, 19–28 (2007). [CrossRef] [PubMed]

], which deals with the study of signals that have many similarities to ocular signals in terms of their non-stationarity, power-law form of the spectral density, and dynamic range [19

B. Boashash and M. Mesbah, “A time-frequency approach for newborn seizure detection,” IEEE Eng. Med. Biol. Mag. 20, 54–64 (2001). [CrossRef] [PubMed]

We initially assume that the power spectral density of the process to be modelled can be approximated by a power-law of the form given in Eq. (2):

P (f) \approx \frac{c}{{| f |}^{γ}}

(2)

where f is the cyclic frequency, γ is the power-law exponent, and c is a constant. In general, for a finite-length, continuous-time signal x(t), the power spectral density can be written as shown in Eq. (3) [1

E. N. Bruce, Biomedical Signal Processing and Signal Modeling (Wiley Series in Telecommunications and Signal Processing, 2001).

P (f) = 〈 \frac{1}{T} {| X (f) |}^{2} 〉

(3)

where the angular brackets denote ensemble average. X(f) is the Fourier transform of x(t), and can be written in the complex exponential form of Eq. (4):

X (f) = | X (f) | e^{i φ (f))}

(4)

where ϕ(f) denotes the phase spectrum. To synthesise a time-domain signal from X(f), we can specify the magnitude and phase spectra, and then take the inverse Fourier transform. We assume that the magnitude spectrum follows a power-law of slope γ/2.

In practice, the method described above can be implemented on a computer using an FFT algorithm. A detailed discussion of the discrete implementation of such an algorithm is given by Billah and Shinozuka [12

K. Billah and M. Shinozuka, “Numerical method for colored noise generation and its application to a bistable system,” Phys. Rev. A 42, 7492–7495 (1990). [CrossRef] [PubMed]

]. Rankine et al. [18

L. Rankine, N. Stevenson, M. Mesbah, and B. Boashash, “A nonstationary model of newborn EEG,” IEEE Trans. Biomed. Eng 54, 19–28 (2007). [CrossRef] [PubMed]

] suggested that the “roughness” of real-world power spectra (as opposed to the “smooth” power-law form of the model) can be approximated by performing the synthesis for several different realisations of ϕ(f), and then summing together the resultant signals in the time-domain.

2.3. Parameter estimation

For estimating the slope value to fit to the PSD, the well-known Hill estimator was used [22

B. M. Hill, “A simple general approach to inference about the tail of a distribution,” Ann. Statist. 3, 1163–1174 (1975). [CrossRef]

]. This method is simple to implement and gives an asymptotically normal and consistent estimate [23

A. Clauset, C. R. Shalizi, and M. E. J. Newman, “Power-law distributions in empirical data,” arXiv:0706.1062v1, URL http://arxiv.org/abs/0706.1062v1 (2007).

]. For the phase spectrum, ϕ(f), one may use the phase spectrum of a measured signal as a “surrogate” phase. Alternatively, Billah and Shinozuka [12

K. Billah and M. Shinozuka, “Numerical method for colored noise generation and its application to a bistable system,” Phys. Rev. A 42, 7492–7495 (1990). [CrossRef] [PubMed]

] as well as Kasdin [21

N. Kasdin, “Discrete simulation of colored noise and stochastic processes and 1/f power law noise,” Proc. IEEE 83, 802–827 (1995). [CrossRef]

] asserted that the phase spectrum of a power-law process can be assumed to be composed of random phase angles uniformly distributed over the interval [0,2π]. Artificial dynamic aberration signals can thus be generated by randomly choosing phase angles from a uniform distribution, to go with a PSD slope value arbitrarily chosen or empirically selected based on estimates from real data.

2.4. Validation of model

To justify the use of the model, we must assess its validity. Unfortunately, for non-stationary signal models, no well-established validation methodology exists [26

P. Celka and P. Colditz, “Nonlinear nonstationary Wiener model of infant seizures,” IEEE Trans. Biomed. Eng. 49, 556–564 (2002). [CrossRef] [PubMed]

]. Rankine et al. [18

L. Rankine, N. Stevenson, M. Mesbah, and B. Boashash, “A nonstationary model of newborn EEG,” IEEE Trans. Biomed. Eng 54, 19–28 (2007). [CrossRef] [PubMed]

] suggested using the Pearson product-moment correlation coefficient to validate simulated non-stationary data, in part due to its widespread familiarity. The application of this measure to processes that exhibit non-stationarity and non-linearity is however questionable. Muma et al. [20

M. Muma, D.R. Iskander, and M. J. Collins, “The role of cardiopulmonary signals in the dynamics of the eye’s wavefront aberrations,” IEEE Trans. Biomed. Eng. 57, 373–383 (2010). [CrossRef]

] described a method of using time-frequency coherence estimates to obtain a normalised estimate of the coupling between two non-stationary signals. We are grateful to the authors for providing us with scripts for the implementation of this method. In the absence of a more rigorous validation approach, we use the time-frequency coherence between real and simulated signals as a means to assess model performance. We do not claim this to be a conclusive or quantitative method of validation (our use of the computed phase spectrum of measured data as a model parameter means that a significant level of coupling is to be expected), but merely an indicator that the simulated data does indeed resemble measured data to a reasonable degree.

3. Results

Figure 1 (left) shows the measured time-evolution of several Zernike coefficients; Fig. 1 (right) shows the corresponding PSD estimates (via the Lomb-Scargle periodogram [24

N. R. Lomb, “Least-squares frequency analysis of unequally spaced data,” Astrophys. & Space Sci. 39, 447–462 (1975). [CrossRef]

,25

J. D. Scargle, “Studies in astronomical time series analysis ii. statistical aspects of spectral analysis of unevenly spaced data,” Astrophys. J. 263, 835–853 (1982). [CrossRef]

]). As mentioned earlier, the temporal behaviour of ocular aberrations is in general non-stationary. This is illustrated by Fig. 2, which gives a time-frequency representation of a measured Zernike astigmatism

(Z_{2}^{- 2})

signal. The time-frequency plot was obtained using the ZAM distribution, a member of Cohen’s class of distribution functions [32

L. Cohen, Time-Frequency Analysis , (Prentice-Hall, New York, 1995).

]. The plot shows low frequency content modulated in the < 0.5 Hz range over the full duration of the 40 s trial. There is significant power at approximately 1 Hz, which also varies in frequency and is likely to be related to cardiopulmonary activity [10

M. Zhu, M. J. Collins, and D. R. Iskander, “Microfluctuations of wavefront aberrations of the eye,” Opthal. Physiol. Opt. 24, 562–571 (2004). [CrossRef]

, 20

M. Muma, D.R. Iskander, and M. J. Collins, “The role of cardiopulmonary signals in the dynamics of the eye’s wavefront aberrations,” IEEE Trans. Biomed. Eng. 57, 373–383 (2010). [CrossRef]

Fig. 1 Dynamics of selected Zernike aberrations (left) and the associated PSD estimates (right).

Fig. 2 Time-frequency representation of Zernike astigmatism for subject CML (via the ZAM distribution).

To assess the validity of our model, we generated signals from the model using parameters obtained directly from real data. From a measured signal, we determined the measured phase spectrum and an estimate $\hat{γ}$ of the power-law slope, and used these as the model parameters ϕ(f) and γ respectively. We then synthesised a time-domain representation using these model parameters. For all measurements, our estimates of the power-law slope on a doubly logarithmic plot were greater than unity, and were similar to the PSD slope estimates of previous authors [3

A. Mira-Agudelo, L. Lundström, and P. Artal, “Temporal dynamics of ocular aberrations: Monocular vs binocular vision,” Ophthal. Physiol. Opt. 29, 256–263 (2009). [CrossRef]

, 7

H. Hofer, P. Artal, and D. R. Williams, “Dynamics of the eyes wave aberration,” J. Opt. Soc. Am. A. 18, 497–506 (2001). [CrossRef]

]. Figure 3 compares the dynamics of measured Zernike spherical aberration to a simulated signal generated using our modelling method, while Fig. 4 shows the corresponding PSD estimates. Figure 5 is a plot of the measured phase spectrum ϕ(f) that is used as the surrogate phase.

Fig. 3 Comparison of measured and simulated Zernike spherical aberration signals in the time-domain.

Fig. 4 Comparison of measured and simulated Zernike spherical aberration signals via their estimated power spectral densities.

Fig. 5 Discrete Fourier phase spectrum computed from measured data and used as the “surrogate” phase to produce the simulated aberration signal of Fig. 3. This Fourier phase spectrum representation was computed using an interpolation algorithm to replace missing data points.

Figure 6 shows the time-frequency coherence between the real and modelled signals of Fig. 1. The value of the time-frequency coherence function is considerably close to unity for most values of t and ω, which indicates that the simulated signal has spectral content that is close to that of the original measured signal from which the model parameters were estimated.

Fig. 6 Time-frequency coherence between real and simulated signals for Zernike spherical aberration.

4. Discussion

4.1. Performance of models

Rankine et al. [18

L. Rankine, N. Stevenson, M. Mesbah, and B. Boashash, “A nonstationary model of newborn EEG,” IEEE Trans. Biomed. Eng 54, 19–28 (2007). [CrossRef] [PubMed]

] suggested that greater flexibility in the modelling of non-stationary data can be achieved by allowing the estimated power-law exponent $\hat{γ}$ to be time-varying, i.e. $\hat{γ}$ = $\hat{γ}$ (t). We have observed that $\hat{γ}$ can vary significantly over a single epoch of data. However, by allowing $\hat{γ}$ to vary over time in our model, we observed only marginal improvement in the model performance at best. This could be partly due to the difficulties in validating model performance. However, we feel that the potential benefits in allowing $\hat{γ}$ to be time-varying must be weighed against some drawbacks. Firstly, the time-scales upon which γ varies significantly are difficult to determine. Secondly, by attempting to estimate γ over shorter time-scales, one must take into account the additional bias in the estimate of γ.

4.2. Two-slope model

The shape of the PSD varies with the level of accommodative effort [27

L.R. Stark and D.A. Atchison, “Pupil size, mean accommodation response and the fluctuations of accommodation,” Opthal. Physiol. Opt. 17, 316–323 (1997). [CrossRef]

, 28

]. When the eye actively accommodates, the shape of the PSD in this range often exhibits two distinct slopes with a “breaking frequency” between them [28

]. This effect is not observed when the eye is in its unaccommodated state, such as in Fig. 1. For the cases where the estimated PSD is better fitted by two slopes rather than one, the model can be adapted accordingly. A bilinear fit can be performed by adopting a regression approach and imposing a constraint on the regression function. In this case, $\hat{γ}$ = $\hat{γ}$ (f), with $\hat{γ}$ (f) given by Eq. (5):

\hat{γ} (f) = {\begin{matrix} γ_{1} & if f \leq f_{b r} \\ γ_{2} & if f > f_{b r} \end{matrix}

(5)

where γ ₁ and γ ₂ are constants representing estimates of the PSD slope for lower and higher frequencies respectively. The two slope regions are separated by the breaking frequency f_br .

4.3. Relationship Between Power-Law Slope and Stationarity

Kasdin [21

N. Kasdin, “Discrete simulation of colored noise and stochastic processes and 1/f power law noise,” Proc. IEEE 83, 802–827 (1995). [CrossRef]

] stipulated that for power-law processes, the stationarity of the process is related to the slope γ of the PSD. In particular, the slope of the PSD at low frequencies is indicative of whether the process is stationary or not. This can be summarised as follows:

A value of γ > 1 at low frequencies implies inherent non-stationarity of the underlying process.
Conversely, a value of 0 ≤ γ ≤ 1 at these lower frequencies implies that the process is stationary.

We illustrate this notion in Fig. 7, with some examples of signals synthesised using the two-slope model. The higher frequency slope was fixed to γ ₂ = 2.0. We then generated signals using five different values of γ ₁. In each case, an identical random phase spectrum (composed of uniform random variates distributed over the interval [0,2π]) was used. To employ a random phase spectrum (rather than a measured one), is potentially more useful as a general approach to modelling the temporal behaviour of the eye’s aberrations. Both methods can be used to generate either stationary or non-stationary signals with frequency content shaped by $\hat{γ}$ (f), as required. As the value of γ ₁ exceeds unity, the effect on the stationarity of the signals becomes clearly visible (for example, the mean of the process wanders significantly). We can consider the value of γ ₁ to relate to the “statistical stability” of the signal, rather than to stationarity in a more rigorous sense. It is useful in this context to think of stationarity in a relative sense, as truly stationary signals are mathematical constructs that do no not typically occur in real-world data [29

M. B. Priestley, Non-Linear and Non-Stationary Time Series Analysis (Academic Press, London, United Kingdom, 1988).

Fig. 7 Illustration of the relationship between statistical stability about a mean and the PSD slope at low frequencies. All signals were generated using the two-slope model, with an identical phase spectrum generated from a uniform distribution.

4.4. Application of models

The model for the dynamics of ocular aberrations was originally conceived as a way to extend static models of aberrations over a large population, such as the model developed by Thibos et al. [30

L. N. Thibos, A. Bradley, and X. Hong, “A statistical model of the aberration structure of normal, well-corrected eyes,” Ophthal. Physiol. Opt. 22, 427–433 (2002). [CrossRef]

]. If the dynamic model parameters could be similarly estimated over a large population, the authors’ idea of creating a database of “virtual eyes” could be extended to include the dynamics of aberrations.

A notable feature of the aberration signals is their dependence on the mean accommodative effort [28

]. The modelling technique described in this work could similarly be used to model dynamic accommodation, as measures of accommodation can be computed from Zernike terms [31

L. N. Thibos, X. Hong, A. Bradley, and X. Cheng, “Statistical variation of aberration structure and image quality in a normal population of healthy eyes,” J. Opt. Soc. Am. A 19, 2329–2348 (2002). [CrossRef]

]. A model of dynamic accommodation could be adapted to include the state-dependence by modelling the relationship between the power-law slope values and the accommodative effort.

It should be noted that in our modelling efforts we have not attempted to compensate for the movements of the eye during measurement. Given that the ultimate aim of this work is to simulate the temporal behaviour of the eye’s intrinsic aberrations (rather than those measured by the wavefront sensor), it follows that these movements should ideally be compensated for, e.g., by using some active tracking of the pupil.

Another property that we have not directly addressed in our modelling approach is the non-linearity of aberration and accommodation signals. We previously discussed the presence of harmonic distortion in our measured signals, which suggests the presence of non-linearity. This is not surprising for a physiological process. Non-linearity has been directly addressed in EEG signal models [26

P. Celka and P. Colditz, “Nonlinear nonstationary Wiener model of infant seizures,” IEEE Trans. Biomed. Eng. 49, 556–564 (2002). [CrossRef] [PubMed]

], and we feel that proper study of non-linearity in ocular signals is also warranted (e.g., using bispectrum or bicoherence methods) for the improvement of the modelling effort.

Acknowledgments

This research was funded by Science Foundation Ireland under grant number 07/IN.1/I906 and the Irish Research Council for Science, Engineering, and Technology. The authors wish to thank Dr. Charles Leroux at NUI Galway and Dr. Luis-Diaz Santana at City University, London for inspiring discussions and guidance during this project.

References and links

1.	E. N. Bruce, Biomedical Signal Processing and Signal Modeling (Wiley Series in Telecommunications and Signal Processing, 2001).
2.	C. Roberts, “Future challenges to aberration-free ablative procedures,” J. Refract. Surg. 16, 623–629 (2000).
3.	A. Mira-Agudelo, L. Lundström, and P. Artal, “Temporal dynamics of ocular aberrations: Monocular vs binocular vision,” Ophthal. Physiol. Opt. 29, 256–263 (2009). [CrossRef]
4.	R. Navarro, E. Moreno-Barriuso, S. Bará, and T. Mancebo, “Phase plates for wave-aberration compensation in the human eye,” Opt. Lett. 25, 236–238 (2000). [CrossRef]
5.	S. O. Galetskiy, T. Yu. Cherezova, and A. V. Kudryashov, “Adaptive optics in ophthalmology: human eye wavefront generator,” Proc. SPIE 6849, 68490918 (2008).
6.	D. R. Iskander, M. Collins, M. Morelande, and M. Zhu, “Analyzing the dynamic wavefront aberrations in the human eye,” IEEE Trans. Biomed. Eng 51, 1969–1980 (2004). [CrossRef] [PubMed]
7.	H. Hofer, P. Artal, and D. R. Williams, “Dynamics of the eyes wave aberration,” J. Opt. Soc. Am. A. 18, 497–506 (2001). [CrossRef]
8.	K. Hampson, E. Mallen, and J. C. Dainty, “Coherence function analysis of the higher-order aberrations of the human eye,” Opt. Lett. 31, 184–186 (2006). [CrossRef] [PubMed]
9.	K. Hampson, I. Munro, C. Paterson, and J. C. Dainty, “Weak correlation between the aberration dynamics of the human eye and the cardiopulmonary system,” J. Opt. Soc. Am. A. 22, 1241–1250 (2005). [CrossRef]
10.	M. Zhu, M. J. Collins, and D. R. Iskander, “Microfluctuations of wavefront aberrations of the eye,” Opthal. Physiol. Opt. 24, 562–571 (2004). [CrossRef]
11.	D. R. Iskander, M. R. Morelande, and M. J. Collins, “Estimating the dynamics of aberration components in the human eye,” In IEEE Signal Process. Workshop Statist. pages241–244 (2001).
12.	K. Billah and M. Shinozuka, “Numerical method for colored noise generation and its application to a bistable system,” Phys. Rev. A 42, 7492–7495 (1990). [CrossRef] [PubMed]
13.	N. Davies, L. Diaz-Santana, and D. Lara-Sucedo, “Repeatability of ocular wavefront measurement,” Optom. Vis. Sci. 80, 142–150 (2003). [CrossRef] [PubMed]
14.	L. Diaz-Santana, V. Guériaux, G. Arden, and S. Gruppetta, “New methodology to measure the dynamics of ocular wave front aberrations during small amplitude changes of accommodation,” Opt. Express 15, 5649–5663 (2007). [CrossRef] [PubMed]
15.	L. Diaz-Santana, C. Torti, I. Munro, P. Gasson, and C. Dainty, “Benefit of higher closed-loop bandwidths in ocular adaptive optics,” Opt. Express 11, 2597–2605 (2003). [CrossRef] [PubMed]
16.	S. Gruppetta, F. Lacombe, and P. Puget, “Study of the dynamic aberrations of the human tear film,” Opt. Express 13, 7631–7636 (2005). [CrossRef] [PubMed]
17.	K. Y. Li and G. Yoon, “Changes in aberrations and retinal image quality due to tear film dynamics,” Opt. Express 14, 12552–12559 (2006). [CrossRef] [PubMed]
18.	L. Rankine, N. Stevenson, M. Mesbah, and B. Boashash, “A nonstationary model of newborn EEG,” IEEE Trans. Biomed. Eng 54, 19–28 (2007). [CrossRef] [PubMed]
19.	B. Boashash and M. Mesbah, “A time-frequency approach for newborn seizure detection,” IEEE Eng. Med. Biol. Mag. 20, 54–64 (2001). [CrossRef] [PubMed]
20.	M. Muma, D.R. Iskander, and M. J. Collins, “The role of cardiopulmonary signals in the dynamics of the eye’s wavefront aberrations,” IEEE Trans. Biomed. Eng. 57, 373–383 (2010). [CrossRef]
21.	N. Kasdin, “Discrete simulation of colored noise and stochastic processes and 1/f power law noise,” Proc. IEEE 83, 802–827 (1995). [CrossRef]
22.	B. M. Hill, “A simple general approach to inference about the tail of a distribution,” Ann. Statist. 3, 1163–1174 (1975). [CrossRef]
23.	A. Clauset, C. R. Shalizi, and M. E. J. Newman, “Power-law distributions in empirical data,” arXiv:0706.1062v1, URL http://arxiv.org/abs/0706.1062v1 (2007).
24.	N. R. Lomb, “Least-squares frequency analysis of unequally spaced data,” Astrophys. & Space Sci. 39, 447–462 (1975). [CrossRef]
25.	J. D. Scargle, “Studies in astronomical time series analysis ii. statistical aspects of spectral analysis of unevenly spaced data,” Astrophys. J. 263, 835–853 (1982). [CrossRef]
26.	P. Celka and P. Colditz, “Nonlinear nonstationary Wiener model of infant seizures,” IEEE Trans. Biomed. Eng. 49, 556–564 (2002). [CrossRef] [PubMed]
27.	L.R. Stark and D.A. Atchison, “Pupil size, mean accommodation response and the fluctuations of accommodation,” Opthal. Physiol. Opt. 17, 316–323 (1997). [CrossRef]
28.	C. Leahy, C. Leroux, C. Dainty, and L. Diaz-Santana, “Temporal dynamics and statistical characteristics of the microfluctuations of accommodation: Dependence on the mean accommodative effort,” Opt. Express 18, 2668–2681 (2010). [CrossRef] [PubMed]
29.	M. B. Priestley, Non-Linear and Non-Stationary Time Series Analysis (Academic Press, London, United Kingdom, 1988).
30.	L. N. Thibos, A. Bradley, and X. Hong, “A statistical model of the aberration structure of normal, well-corrected eyes,” Ophthal. Physiol. Opt. 22, 427–433 (2002). [CrossRef]
31.	L. N. Thibos, X. Hong, A. Bradley, and X. Cheng, “Statistical variation of aberration structure and image quality in a normal population of healthy eyes,” J. Opt. Soc. Am. A 19, 2329–2348 (2002). [CrossRef]
32.	L. Cohen, Time-Frequency Analysis , (Prentice-Hall, New York, 1995).

OCIS Codes
(330.0330) Vision, color, and visual optics : Vision, color, and visual optics
(330.4060) Vision, color, and visual optics : Vision modeling

ToC Category:
Vision, Color, and Visual Optics

History
Original Manuscript: July 6, 2010
Revised Manuscript: September 9, 2010
Manuscript Accepted: September 22, 2010
Published: September 23, 2010

Virtual Issues
Vol. 5, Iss. 14 Virtual Journal for Biomedical Optics

Citation
C. Leahy and C. Dainty, "A non-stationary model for simulating the dynamics of ocular aberrations," Opt. Express 18, 21386-21396 (2010)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-18-20-21386

Sort: Author | Year | Journal | Reset

References

E. N. Bruce, Biomedical Signal Processing and SignalModeling (Wiley Series in Telecommunications and Signal Processing, 2001).
C. Roberts, "Future challenges to aberration-free ablative procedures," J. Refract. Surg. 16, 623-629 (2000).
A. Mira-Agudelo, L. Lundström, and P. Artal, "Temporal dynamics of ocular aberrations:Monocular vs binocular vision," Ophthal. Physiol. Opt. 29, 256-263 (2009). [CrossRef]
R. Navarro, E. Moreno-Barriuso, S. Bará, and T. Mancebo, "Phase plates for wave-aberration compensation in the human eye," Opt. Lett. 25, 236-238 (2000). [CrossRef]
S. O. Galetskiy, T. Yu. Cherezova, and A. V. Kudryashov, "Adaptive optics in ophthalmology: human eye wavefront generator," Proc. SPIE 6849, 68490918 (2008).
D. R. Iskander, M. Collins, M. Morelande, and M. Zhu, "Analyzing the dynamic wavefront aberrations in the human eye," IEEE Trans. Biomed. Eng. 51, 1969-1980 (2004). [CrossRef] [PubMed]
H. Hofer, P. Artal, and D. R. Williams, "Dynamics of the eyes wave aberration," J. Opt. Soc. Am. A. 18, 497-506 (2001). [CrossRef]
K. Hampson, E. Mallen, and J. C. Dainty, "Coherence function analysis of the higher-order aberrations of the human eye," Opt. Lett. 31, 184-186 (2006). [CrossRef] [PubMed]
K. Hampson, I. Munro, C. Paterson, and J. C. Dainty, "Weak correlation between the aberration dynamics of the human eye and the cardiopulmonary system," J. Opt. Soc. Am. A. 22, 1241-1250 (2005). [CrossRef]
M. Zhu, M. J. Collins, and D. R. Iskander, "Microfluctuations of wavefront aberrations of the eye," Opthal. Physiol. Opt. 24, 562-571 (2004). [CrossRef]
D. R. Iskander, M. R. Morelande, and M. J. Collins, "Estimating the dynamics of aberration components in the human eye," In IEEE Signal Process. Workshop Statist. pages 241-244 (2001).
K. Billah and M. Shinozuka, "Numerical method for colored noise generation and its application to a bistable system," Phys. Rev. A 42, 7492-7495 (1990). [CrossRef] [PubMed]
N. Davies, L. Diaz-Santana, and D. Lara-Sucedo, "Repeatability of ocular wavefront measurement," Optom. Vis. Sci. 80, 142-150 (2003). [CrossRef] [PubMed]
L. Diaz-Santana, V. Guériaux, G. Arden, and S. Gruppetta, "New methodology to measure the dynamics of ocular wave front aberrations during small amplitude changes of accommodation," Opt. Express 15, 5649-5663 (2007). [CrossRef] [PubMed]
L. Diaz-Santana, C. Torti, I. Munro, P. Gasson, and C. Dainty, "Benefit of higher closed-loop bandwidths in ocular adaptive optics," Opt. Express 11, 2597-2605 (2003). [CrossRef] [PubMed]
S. Gruppetta, F. Lacombe, and P. Puget, "Study of the dynamic aberrations of the human tear film," Opt. Express 13, 7631-7636 (2005). [CrossRef] [PubMed]
K. Y. Li and G. Yoon, "Changes in aberrations and retinal image quality due to tear film dynamics," Opt. Express 14, 12552-12559 (2006). [CrossRef] [PubMed]
L. Rankine, N. Stevenson,M. Mesbah, and B. Boashash, "A nonstationary model of newborn EEG," IEEE Trans. Biomed. Eng 54, 19-28 (2007). [CrossRef] [PubMed]
B. Boashash and M. Mesbah, "A time-frequency approach for newborn seizure detection," IEEE Eng. Med. Biol. Mag. 20, 54-64 (2001). [CrossRef] [PubMed]
M. Muma, D. R. Iskander, and M. J. Collins, "The role of cardiopulmonary signals in the dynamics of the eye’s wavefront aberrations," IEEE Trans. Biomed. Eng. 57, 373-383 (2010). [CrossRef]
N. Kasdin, "Discrete simulation of colored noise and stochastic processes and 1/f power law noise," Proc. IEEE 83, 802-827 (1995). [CrossRef]
B. M. Hill "A simple general approach to inference about the tail of a distribution," Ann. Statist. 3, 1163-1174 (1975). [CrossRef]
A. Clauset, C. R. Shalizi, and M. E. J. Newman, "Power-law distributions in empirical data," arXiv:0706.1062v1, URL http://arxiv.org/abs/0706.1062v1 (2007).
N. R. Lomb, "Least-squares frequency analysis of unequally spaced data," Astrophys. Space Sci. 39, 447-462 (1975). [CrossRef]
J. D. Scargle, "Studies in astronomical time series analysis ii. statistical aspects of spectral analysis of unevenly spaced data," Astrophys. J. 263, 835-853 (1982). [CrossRef]
P. Celka and P. Colditz, "Nonlinear nonstationary Wiener model of infant seizures," IEEE Trans. Biomed. Eng. 49, 556-564 (2002). [CrossRef] [PubMed]
L. R. Stark and D. A. Atchison, "Pupil size, mean accommodation response and the fluctuations of accommodation," Opthal. Physiol. Opt. 17, 316-323 (1997). [CrossRef]
C. Leahy, C. Leroux, C. Dainty, and L. Diaz-Santana, "Temporal dynamics and statistical characteristics of the microfluctuations of accommodation: Dependence on the mean accommodative effort," Opt. Express 18, 2668-2681 (2010). [CrossRef] [PubMed]
M. B. Priestley, Non-Linear and Non-Stationary Time Series Analysis (Academic Press, London, United Kingdom, 1988).
L. N. Thibos, A. Bradley, and X. Hong, "A statistical model of the aberration structure of normal, well-corrected eyes," Ophthal. Physiol. Opt. 22, 427-433 (2002). [CrossRef]
L. N. Thibos, X. Hong, A. Bradley, and X. Cheng, "Statistical variation of aberration structure and image quality in a normal population of healthy eyes," J. Opt. Soc. Am. A 19, 2329-2348 (2002). [CrossRef]
L. Cohen, Time-Frequency Analysis, (Prentice-Hall, New York, 1995).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper