New Approaches, Findings, and Diagnostics in Medical and Surgical Retina
View this Special IssueResearch Article  Open Access
Advance in ERG Analysis: From Peak Time and Amplitude to Frequency, Power, and Energy
Abstract
Purpose. To compare time domain (TD: peak time and amplitude) analysis of the human photopic electroretinogram (ERG) with measures obtained in the frequency domain (Fourier analysis: FA) and in the timefrequency domain (continuous (CWT) and discrete (DWT) wavelet transforms). Methods. Normal ERGs were analyzed using traditional peak time and amplitude measurements of the a and bwaves in the TD and descriptors extracted from FA, CWT, and DWT. Selected descriptors were also compared in their ability to monitor the longterm consequences of disease process. Results. Each method extracted relevant information but had distinct limitations (i.e., temporal and frequency resolutions). The DWT offered the best compromise by allowing us to extract more relevant descriptors of the ERG signal at the cost of lesser temporal and frequency resolutions. Followups of disease progression were more prolonged with the DWT (max 29 years compared to 13 with TD). Conclusions. Standardized time domain analysis of retinal function should be complemented with advanced DWT descriptors of the ERG. This method should allow more sensitive/specific quantifications of ERG responses, facilitate followup of disease progression, and identify diagnostically significant changes of ERG waveforms that are not resolved when the analysis is only limited to time domain measurements.
1. Introduction
The electroretinogram (ERG) identifies the electrical signal that is generated by the retina in response to a light stimulus. It is the first biopotential ever recorded from a human subject, namely, by Dewar in 1877 [1]. However, despite significant (r)evolution in the recording technologies (essentially from the string galvanometer to the digital amplifier and supporting computer software) and, consequently, the significantly enhanced quality of the ERG signal thus obtained, analysis of the ERG remains for the most part limited to amplitude and peak time measurements of its major components, namely, the a and bwaves. This is at least what is recommended in the ERG standard of the International Society for Clinical Electrophysiology of Vision (ISCEV) [2]. The a and bwaves of the ERG are said to reflect the activity generated by the photoreceptors and the bipolarMüller cell complex, respectively [3–5]. These components are usually referred to as the slow waves of the ERG. Also identified in the ERG signal are the small, highfrequency, oscillations that are often seen riding on the ascending limb of the bwave [6, 7]. These components, referred to as oscillatory potentials (OPs), are most probably generated by the retinal cells of the inner retina (i.e., bipolar, amacrine, or horizontal cells) although their exact origin remains debated [8, 9]. The OPs appear to be major contributors to the shaping of the ERG waveform [10] and there is an abundant literature attesting to the clinical value of including the OPs when analyzing pathological ERGs [6, 11, 12]. Unfortunately, in order to optimize the visualization of the OPs one must modify the recording bandwidth of the ERG from a broadband (e.g., 1–1000 Hz) to a narrower band (e.g., 100–1000 Hz) that removes the lowfrequency components of the ERG (i.e., a and bwaves) and consequently selectively enhances the highfrequency components (i.e., OPs) [2]. However, when doing so one must always keep in mind the possibility of introducing artifactual components (such as ringing artifacts and phase lags) to the ERG thus obtained.
It is clear from the above that the ERG waveform results from the amalgamation of several frequency components. Is it possible to monitor the frequency composition of the ERG signal without altering the signal as it is done with the bandwidth restriction approach? Would the use of such an approach significantly improve analysis of the ERG beyond what is accomplished when using time and amplitude measures of the ERG only? Although advanced analytical approaches are now frequently used when studying biopotentials, such as the electroencephalogram [13], the electrocardiogram [14], and the electromyogram [15], to date they have only been sporadically applied to the ERG [8, 16, 17]. The purpose of this study was therefore to compare peak time and amplitude measurements of human photopic ERGs with measures obtained in the frequency domain using the Fourier analysis as well as in the timefrequency domain using the continuous (CWT) and discrete (DWT) wavelet transforms. For the sake of brevity, our study was limited to the photopic ERG only.
2. Materials and Methods
Normal photopic ERGs were obtained from 40 healthy subjects (26 females and 14 males, average age years) using a protocol that was approved by the Institutional Review Board of the Montreal Children’s Hospital and in accordance with the Declaration of Helsinki.
According to a previously published method of ours, the ERGs were recorded with both eyes dilated (tropicamide 1%) using an active electrode (DTL fiber electrode) placed in the inferior conjunctival bag, with reference and ground electrodes pasted at the external canthi and forehead, respectively [18–21]. The ERGs were evoked to flashes of white light (flash duration: 20 μs; interstimulus interval: 1.5 s; and average of at least 10 flashes per recording) of 0.64 log cd·s·m^{−2} in intensity that were delivered against a rod desensitizing background light of 30 cd·m^{−2} (measured using a research radiometer IL1700; International Light, Newburyport, MA, USA). ERG waves from both eyes were averaged to yield a single waveform of 150 ms in length (sampling rate: 3413.33 Hz) that included a prestimulus baseline of 20 ms.
2.1. ERG Analysis
The amplitude of the awave was measured from the prestimulus baseline to the most negative trough of the ERG, while the amplitude of the bwave was measured from the trough of the awave to the most positive peak of the ERG that followed the awave [2]. Peak times were measured from flash onset to the peak of the a and bwaves [2]. Given that these measures of the ERG are taken in the time domain, they will be referred to as time domain (TD) measurements.
Frequency domain analysis (or Fourier analysis (FA)) of the ERG was carried out using the fast Fourier transform (FFT) algorithm implemented in MATLAB R2013b (Mathworks, Natick, MA, USA) as follows: where represents the FFT coefficients, denotes the raw ERG timeseries, and denotes the number of data points in . Each FFT coefficient weighs the energetic contribution of a single frequency component to the signal so that a frequency spectrum can be illustrated by tracing . Considering the size (512 data points) and sampling frequency (3413.33 Hz) of our ERG waveforms, we were able to compute the FFT coefficients for frequencies ranging between 0 and 1706.66 Hz in increments of 6.66 Hz (i.e., frequency resolution). However, given the limitation imposed by our recording bandwidth (1–1000 Hz) we limited our analysis to frequencies ranging between 0 and 300 Hz to safely avoid artifactual contamination (such as that predicted by the NyquistShannon sampling theorem [22]).
In order to localize the energy content of the ERG in both time and frequency we computed, using MATLAB, the continuous wavelet transform (CWT) of selected ERGs as follows: where represents the wavelet coefficients localized at scales (frequency) and moments (time), denotes the unprocessed ERG timeseries, and denotes the complex conjugate of the Morse wavelet [23], which was chosen for its good frequency resolution [24]. To illustrate the timefrequency scalogram of the CWT, we took the absolute value of and normalized it to its maximal value so that the timefrequency scalograms of the ERG are shown as colored twodimensional plots of in which minimum energy values are displayed in blue and maximal values in red.
Use of the CWT approach allowed us to analyse the ERG, continuously, at every possible scale and translation . This approach, however, requires extensive computation time and also yields a lot of redundant information (i.e., since each coefficient has similar neighboring values) that will remain unused in the set of coefficients [25]. Interestingly, if the scale and translation parameters of the wavelet are taken at discrete values, we then obtain a discrete wavelet transform (DWT), where the scales and translations are based on powers of two (i.e., , ) so that (2) can be discretized as follows: where represents the wavelet coefficients localized at discrete scales (frequency) and discrete moments (time), designates the raw ERG timeseries, and designates the Haar wavelet [25]. was computed using the fast wavelet transform algorithm of Mallat [25, 26] implemented in MATLAB. We chose the Haar wavelet, for its simplicity (i.e., simplest wavelet available) and its orthonormal basis, which allows the wavelet coefficients to be reconstructed accurately and efficiently without any loss of information (using the inverse DWT) even if all the redundant information contained in the CWT is discarded [25, 27]. Similar to the CWT scalograms, the most prominent energy component of the DWT scalogram will appear as a deep dark red (high energy) rectangle in the region of the DWT scalogram where it is located (i.e., located in time and frequency), and, conversely, the absence of any component at given locations will appear as deep dark blue (no energy) rectangles.
2.2. Statistical Analyses
Mean value, standard deviation (SD), and coefficient of variation (CV) were computed for all ERG parameters that were identified using the different analytical approach. scores were used to evaluate the significance of selected descriptor changes. All tests were set to a level of significance of 5%.
3. Results
3.1. Time and Amplitude Measurements in the Time Domain
As reported in Table 1, time domain analysis allowed the identification of two major ERG components, one peaking at ms (mean ± SD obtained in our 40 subjects) with an amplitude of μV (identified as the awave at Figure 1(a)) and another one which peaks at ms with an amplitude of μV (identified as the bwave at Figure 1(a)).

(a)
(b)
(c)
(d)
3.2. Fourier Analyses (FA)
The frequency components contributing to the genesis of the ERG can be identified using the FA, such as that obtained with the FFT. This is best exemplified in Figure 1 and Table 1, where three major frequency components are identified in the normal ERGs. As shown with the black arrows, the lowfrequency components of the ERG (presumably a and bwaves) usually formed a smooth peak of large magnitude, culminating at Hz (i.e., mean ± SD obtained from our 40 subjects) on the frequency spectrum. However, in some instances, two distinct peaks could be identified (see double black arrows in panels (a) and (d)). In the later cases, only the peak of highest magnitude was considered for further analysis. In contrast, the higher frequency components of the ERG (probably the OPs) usually formed two distinct peaks (see grey arrows) of low magnitude located at Hz and Hz, respectively. However, given that FA only looks at the frequency content of the ERG without taking into consideration if those frequencies are timelocked or not to the stimulus, its use can lead to erroneous interpretations of the ERG. This is best illustrated with the ERGs shown in Figures 1(c) and 1(d) where the noise contaminants (such as 60 Hz in Figure 1(c)) appear to contribute more to the making of the ERG than the retinal evoked components themselves. These limitations of FA can be overcome by adding a temporal resolution to the frequency domain.
3.3. Continuous Wavelet Transform (CWT) Analyses
As shown in Figure 2, use of the CWT approach allowed us to more precisely localize (in both time and frequency, as reported in Table 1) the abovementioned frequency components of the ERG signal and even in the presence of significant noise contaminants (as seen in panels (c) and (d)). In each scalogram, the a and bwaves formed a cluster of hot (dark red) coefficients (see white arrows) centered at Hz (i.e., mean ± SD obtained in our 40 subjects) and peaking at ms (i.e., timelocked to the peak time of the ERG bwave). Similarly, the OPs formed two distinct clusters (see grey arrows) centered at Hz and Hz and peaking at ms and ms, respectively. As shown in the scalogram of panels C and D, the high (i.e., OPs) and lowfrequency (i.e., corresponding to the a and bwaves) components continued to remain the major lightevoked (i.e., timelocked to the stimulus) components of the ERG response in spite of significant noise contamination. The latter contrasts with results obtained using the FA where one cannot dissociate evoked from nonevoked frequency components (compare results shown in Figures 2(c) and 2(d) with corresponding FA results shown in Figures 1(c) and 1(d)).
(a)
(b)
(c)
(d)
In each scalogram of Figure 2, the value of the coefficients that identified the hot clusters pertaining to the a and bwaves or OPs was centered on a section of the clusters where all coefficients had the same value (i.e., equipotential regions). This attribute of the CWT will generate redundant information (i.e., similar or equal coefficient values) that will complicate the accurate identification of the timefrequency coordinates (i.e., illposed coordinates) of the ERG components, by introducing an uncertainty factor. For example, given that the red clusters of Figure 2 extend over a large area of the scalograms, this prevents an accurate quantification of the awave energy, which is most probably hidden by the higher energy bwave. These limitations can be overcome if we impose a discretization over the possible frequencies and times at which the information is computed.
3.4. Discrete Wavelet Transform (DWT) Analysis
As illustrated in Figure 3, the DWT scalograms decomposed the signals into seven contiguous frequency bands (20, 40, 80, 160, 320, 640, and 1280 Hz), each including a range of frequencies (reported in Table 1) around their respective central frequency (CF). For example, the 20 Hz band quantified the energy oscillating between 13.33 and 26.66 Hz (i.e., /3, that is, CF ± CF/3). Use of this expansion simplifies the choice of relevant coefficients (i.e., seen as rectangles of different sizes in the scalograms of Figure 3) that may be used as energy descriptors of the ERG.
3.4.1. Identification of DWT Descriptors
As seen with the DWT scalograms of Figure 3, the major frequency components are confined to the timefrequency region that is surrounded by white borders and where six major components can be identified (i.e., rectangles of various colors shown in panel (a) and magnified in panel (b)) and thus quantified (as reported in Table 1). The DWT also removed the redundancy so that the bwave energy, quantified with the 20b and 40b descriptors (panel (b)), is now represented by two rectangles (i.e., identified as 20b and 40b) confined to an area of the scalogram that is limited to the vicinity of the bwave peak rather than spread across most of the CWT scalogram. This allowed the quantification of two descriptors, timelocked to the awave, which we identified as 20a and 40a. The 20a and 40a were of lower energy compared to the 20b or 40b (see Table 1), indicating that, as expected, the bwave energy is greater than that of the awave. Finally, the highfrequency components, indicated as 80ops and 160ops, were also easily identified (i.e., maximal value in the 80 and 160 Hz bands, resp.).
3.4.2. Improvement of ERG Segregation Using the DWT
As indicated above, the a and bwave components were seen on both the 20 Hz (20a and 20b descriptors) and 40 Hz bands (40a and 40b descriptors). These descriptors can be used to segregate ERGs of different morphologies. This is better illustrated in Figure 3(c), where two ERGs of distinct morphologies were similar () in terms of a and bwave peak time and amplitude but were significantly different () on the basis of their DWT bwave descriptors (20b and 40b). In one example, the ascending limb of the bwave has a sharp morphology (blue tracing) and showed a lower 20b descriptor (), compared to the broader ascending limb of the bwave (compare first half of ascension of tracing in red) which disclosed an attenuated 40b parameter (). Similarly, we were also able to segregate ERGs that differed in OPs prominence (such as what is shown in Figure 3(d)). Although these ERGs were indistinguishable () on the basis of peak time and amplitude measurements of the a and bwaves, they were significantly different () on the basis of their 80ops and 160ops energy content which were higher (blue tracing) or lower (red tracing) compared to average.
3.5. Applications of Refined Analytical Approach to Clinical ERGs
In Figure 4(a) the ERG waveforms obtained from a patient diagnosed with retinitis pigmentosa (RP) that was followed up for more than 30 years are illustrated. As shown, the low amplitude (and low signaltonoise ratio or low SNR) of these pathological ERGs, especially those obtained later in the disease process (tracings 23 and 29), seriously compromises an accurate measurement of these waveforms. This is best exemplified in Figure 4(c), where an accurate measurement of the bwave amplitude could only be achieved for ERGs obtained within the first 13 years, due to the highly contaminated ERGs recorded subsequently. However, use of the DWT still permitted the extraction of bwave descriptors (i.e., 20b and 40b) as shown in the scalograms of Figure 4(b) and therefore allowed us to monitor progression of the disease process for an additional period of 16 years, as shown in Figure 4(d). Furthermore, as revealed in Figure 4(d), both eyes followed the same degeneration pattern, which appeared to follow an exponential decay function that correlated well with that which characterized (using bwave amplitude measures) the first 13 years. Interestingly, the use of the inverse DWT of the lowfrequency bands (i.e., 20 and 40 Hz bands) allowed us to reconstruct noisefree ERGs (Figure 4(e)), that were nearly identical in both eyes, as it was also the case for the ERGs (measurable, high SNR) recorded in earlier exams (tracings 0, 1, 3, and 9 of Figure 4(a)). The validation of this denoising approach (i.e., inverse DWT) is further demonstrated in Figure 5, where we reconstructed the 20 consecutive singleflash recordings obtained from a RP patient that were used to generate the average waveform. As shown, the 20 denoised singleflash waveforms (red tracings) are nearly identical (mean (±SD) Pearson coefficients = ) to the averaged response (i.e., blue tracing obtained by averaging the 20 consecutive noisy responses (gray traces)). In contrast, the mean Pearson coefficients obtained between the single sweeps and the averaged response were of .
4. Discussion
To date, analysis of the ERG relies mostly on time domain (TD) measurements (peak time and amplitude) of its two major components, namely, the a and bwaves. However, as shown with the examples illustrated in Figures 4 and 5 and as previously suggested [28–31], TD measurements are subjected to noise contamination. These contaminants can arise from numerous factors such as the subject (e.g., eye blinks, head/eye movements, etc.), external sources (e.g., mechanical vibrations, electromagnetic coupling with the 50/60 Hz power lines, computer monitors, electrical lighting, etc.), and, if the data is digitized, the digitization process itself (e.g., digitization artifacts, aliasing, etc.). Therefore, limiting the ERG analysis to TD measures only could jeopardize the detection of subtle functional changes.
4.1. From Fourier to Wavelets
FA methods, such as the one presented in Figure 1 (i.e., FFT), performed well in identifying the three major frequency components of the normal photopic ERG response. However, when a noise contaminant was distributed over the entire ERG response, the resulting frequencypower distribution was misleading (as shown in Figures 1(c) and 1(d)). This is due to the fact that the FA assumes that all the frequency components that compose a signal are periodic and, consequently, ignores the possibility that some frequencies could be found at precise poststimulus time locations only. This explains why the amplitude of a component can be over or underestimated when relying solely on FA analysis. In other words, while Fourier analyses are wellsuited to identify the frequency components that compose the ERG signal, they are of no use to determine their respective magnitude and temporal location within the signal. Such information is of crucial importance if one wishes to define the signal with all its subtleties. This can be accomplished with a timefrequency domain analysis of the ERG.
Use of the CWT (Figure 2) approach allowed us to clearly identify the temporal as well as the frequency coordinates of the three major frequency components of the ERG, thus remedying the FA limitations alluded to above. We have shown that, with the CWT scalograms, each of these components was timelocked to the largest wave of the ERG (i.e., the bwave). Furthermore, compared to corresponding FA estimates, their respective weights (i.e., relative energy levels) were also more accurately determined, even in noisy ERG recordings.
4.2. The DWT: An Optimal Compromise
Interestingly, it seems, from FA and CWT analyses, that the components of the ERG cannot be associated with single frequency values but rather to a range of values. For example, in the FA power spectrums of Figure 1, each component had a broad Gaussianlike distribution (i.e., suggestive of a band of frequencies) rather than a sharp peak (i.e., suggestive of a single frequency) as it would be the case for a pure sinusoid. Similar broad distributions of individual frequency components were also observed in the CWT of the same ERG signals (see Figure 2), the major difference being that the magnitude of these frequency components can now be timecorrelated to specific events of the ERG signal. Consequently, since these components contain a band of frequencies rather than a single frequency, analysis of the ERG at different frequency bands using the DWT scalogram (rather than at each possible frequency with the CWT) offers a simplified scheme (exempt of redundancy) to identify relevant ERG descriptors (see Figures 3(a) and 3(b)).
In the DWT scalograms the a and bwaves were characterized by distinct components located in the 20 Hz (20a and 20b descriptors) and 40 Hz (40a and 40b descriptors) bands. It was difficult to accurately determine these distinct frequency components using the CWT, although in the FA some ERGs did show both the 20 and 40 Hz components (i.e., identified as doublepeaks in Figures 1(b) and 1(d)). Furthermore, quantifying the ERG waveforms using the DWT augmented the specificity (i.e., ability to discriminate between distinct ERG morphologies, as shown in Figures 3(c) and 3(d)) and sensitivity (i.e., reduced the variability, as shown by the CV reported in Table 1) of the measures obtained.
Finally, at the endstage of severe degenerative retinopathies (such as RP), nearly extinguished ERGs (e.g., lowSNR, such as the one shown in Figures 4(a) and 5) are often the last measurable signs of functional vision [32]. As shown in Figure 4, extracting relevant ERG descriptors from these residual ERGs using a TD approach (e.g., peak time and amplitude measures) becomes nearly impossible as the SNR decreases. Use of the DWT permitted the quantification of such responses, thus extending the length of the followup period of disease progression by an additional 16 years. This allowed us to demonstrate the exponential decay known to characterize the longterm course of the cone ERG amplitude in patients affected with RP [33]. Furthermore, use of the inverse DWT of selected frequency bands allowed us to reconstruct noisefree ERG waveform thus confirming the presence of a residual biological response in signals that were reported as nonmeasurable using the TD approach. The validation of this DWTdenoising approach was demonstrated in Figure 5, where each of the 20 denoised ERGs was highly similar to the averaged response (obtained by averaging the 20 consecutive responses).
4.3. Limitations of the Study
In this study we limited the TD approach to its most widespread descriptors (i.e., amplitude and peak time of the a and bwaves), but other unusual descriptors (e.g., areaunderthecurve of the a or bwave, time to reach a certain percentage of the a or bwave amplitudes, steepness of the rising or descending flank of the bwave, filtered OPs measurements, etc. [6, 31, 34]) could also be of use to identify subtle morphological changes, albeit similarly sensitive to noise contaminant errors.
5. Conclusions
In this paper, we have presented a brief overview of the different analytical approaches that can be used to quantify the ERG waveform. As long as the response remains measurable, the traditional measurements of the a and bwaves can be used to monitor the peak time and the amplitude of the ERG signal. However, these measurements only look at the ERG signal as a whole, instead of looking at the different frequency components (possibly of distinct cellular origin) separately. The discrete wavelet transform offers the possibility to extract more components of the ERG signal, even in very poor SNR responses. Standardized time domain analysis of retinal function should thus be complemented with advanced DWT descriptors of the ERG. The latter should allow more sensitive/specific quantifications of ERG responses, facilitate followup of disease progression, and identify diagnostically significant changes of ERG waveforms that are not resolved when the analysis is only limited to time domain measurements, thus bringing the analysis and interpretation of the ERG signal in the 21st century, as it is already the case with other biopotentials such as the electroencephalogram and electrocardiogram.
Conflict of Interests
The authors of the paper (Mathieu Gauvin, JeanMarc Lina, and Pierre Lachapelle) do not have any financial disclosure or conflict of interests to report.
Acknowledgment
This study was funded by Grantsinaid from the Canadian Institutes for Health Research (MOP126082) and the Vision Health Research Network of the Fonds de Recherche du QuébecSanté.
References
 A. F. de Rouck, “History of the electroretinography,” in Principles and Practice of Clinical Electrophysiology of Vision, J. R. Heckenlively and G. B. Arden, Eds., pp. 139–185, The Mit Press, London, UK, 2006. View at: Google Scholar
 M. F. Marmor, A. B. Fulton, G. E. Holder, Y. Miyake, M. Brigell, and M. Bach, “ISCEV Standard for fullfield clinical electroretinography (2008 update),” Documenta Ophthalmologica, vol. 118, no. 1, pp. 69–77, 2009. View at: Publisher Site  Google Scholar
 P. A. Sieving, K. Murayama, and F. Naarendorp, “Pushpull model of the primate photopic electroretinogram: a role for hyperpolarizing neurons in shaping the bwave,” Visual Neuroscience, vol. 11, no. 3, pp. 519–532, 1994. View at: Google Scholar
 D. C. Hood and D. G. Birch, “Human cone receptor activity: the leading edge of the awave and models of receptor activity,” Visual neuroscience, vol. 10, no. 5, pp. 857–871, 1993. View at: Publisher Site  Google Scholar
 R. F. Miller and J. E. Dowling, “Intracellular responses of the Müller (glial) cells of mudpuppy retina: their relation to bwave of the electroretinogram,” Journal of Neurophysiology, vol. 33, no. 3, pp. 323–341, 1970. View at: Google Scholar
 P. Lachapelle, “The human suprathreshold photopic oscillatory potentials: method of analysis and clinical application,” Documenta Ophthalmologica, vol. 88, no. 1, pp. 1–25, 1994. View at: Publisher Site  Google Scholar
 L. Wachtmeister, “Oscillatory potentials in the retina: what do they reveal,” Progress in Retinal and Eye Research, vol. 17, no. 4, pp. 485–521, 1998. View at: Publisher Site  Google Scholar
 W. Zhou, N. Rangaswamy, P. Ktonas, and L. J. Frishman, “Oscillatory potentials of the slowsequence multifocal ERG in primates extracted using the Matching Pursuit method,” Vision Research, vol. 47, no. 15, pp. 2021–2036, 2007. View at: Publisher Site  Google Scholar
 H. Heynen, L. Wachtmeister, and D. van Norren, “Origin of the oscillatory potentials in the primate retina,” Vision Research, vol. 25, no. 10, pp. 1365–1373, 1985. View at: Publisher Site  Google Scholar
 P. Lachapelle, “Oscillatory potentials as predictors to amplitude and peak time of the photopic bwave of the human electroretinogram,” Documenta Ophthalmologica, vol. 75, no. 1, pp. 73–82, 1990. View at: Publisher Site  Google Scholar
 P. Lachapelle, S. Rousseau, M. Mckerral et al., “Evidence supportive of a functional discrimination between photopic oscillatory potentials as revealed with cone and rod mediated retinopathies,” Documenta Ophthalmologica, vol. 95, no. 1, pp. 35–54, 1998. View at: Publisher Site  Google Scholar
 L. Wachtmeister, “Basic research and clinical aspects of the oscillatory potentials of the electroretinogram,” Documenta Ophthalmologica, vol. 66, no. 3, pp. 187–194, 1987. View at: Publisher Site  Google Scholar
 K. Gadhoumi, J. Lina, and J. Gotman, “Discriminating preictal and interictal states in patients with temporal lobe epilepsy using wavelet analysis of intracerebral EEG,” Clinical Neurophysiology, vol. 123, no. 10, pp. 1906–1916, 2012. View at: Publisher Site  Google Scholar
 P. S. Addison, “Wavelet transforms and the ECG: a review,” Physiological Measurement, vol. 26, no. 5, pp. R155–R199, 2005. View at: Publisher Site  Google Scholar
 A. Phinyomark, C. Limsakul, and P. Phukpattaranont, “Application of wavelet analysis in EMG feature extraction for pattern classification,” Measurement Science Review, vol. 11, no. 2, pp. 45–52, 2011. View at: Publisher Site  Google Scholar
 J. D. Forte, B. V. Bui, and A. J. Vingrys, “Wavelet analysis reveals dynamics of rat oscillatory potentials,” Journal of Neuroscience Methods, vol. 169, no. 1, pp. 191–200, 2008. View at: Publisher Site  Google Scholar
 J. M. MiguelJiménez, L. Boquete, S. Ortega, J. M. RodríguezAscariz, and R. Blanco, “Glaucoma detection by waveletbased analysis of the global flash multifocal electroretinogram,” Medical Engineering and Physics, vol. 32, no. 6, pp. 617–622, 2010. View at: Publisher Site  Google Scholar
 M. Rufiange, S. Rousseau, O. Dembinska, and P. Lachapelle, “Conedominated ERG luminanceresponse function: the Photopic Hill revisited,” Documenta Ophthalmologica, vol. 104, no. 3, pp. 231–248, 2002. View at: Publisher Site  Google Scholar
 M. Rufiange, J. Dassa, O. Dembinska et al., “The photopic ERG luminanceresponse function (photopic hill): method of analysis and clinical application,” Vision Research, vol. 43, no. 12, pp. 1405–1412, 2003. View at: Publisher Site  Google Scholar
 M. Rufiange, M. Dumont, and P. Lachapelle, “Modulation of the human photopic ERG luminanceresponse function with the use of chromatic stimuli,” Vision Research, vol. 45, no. 17, pp. 2321–2330, 2005. View at: Publisher Site  Google Scholar
 M. Garon, M. Rufiange, R. Hamilton, D. L. McCulloch, and P. Lachapelle, “Asymmetrical growth of the photopic hill during the light adaptation effect,” Documenta Ophthalmologica, vol. 121, no. 3, pp. 177–187, 2010. View at: Publisher Site  Google Scholar
 H. Nyquist, “Certain topics in telegraph transmission theory,” Proceedings of the IEEE, vol. 90, no. 2, pp. 280–305, 2002. View at: Publisher Site  Google Scholar
 J. M. Lilly and S. C. Olhede, “On the analytic wavelet transform,” IEEE Transactions on Information Theory, vol. 56, no. 8, pp. 4135–4156, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 G. A. Worrell, K. Jerbi, K. Kobayashi, J. M. Lina, R. Zelmann, and M. Le Van Quyen, “Recording and analysis techniques for highfrequency oscillations,” Progress in Neurobiology, vol. 98, no. 3, pp. 265–278, 2012. View at: Publisher Site  Google Scholar
 S. G. Mallat, A Wavelet Tour of Signal Processing: The Sparse Way, Academic Press, Houston, Tex, USA, 3rd edition, 2009. View at: MathSciNet
 J. Buckheit, S. Chen, D. Donoho, and I. Johnstone, About WaveLab, Stanford University Press, Stanford, Calif, USA, 1995.
 S. G. Mallat, “Theory for multiresolution signal decomposition: the wavelet representation,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 11, no. 7, pp. 674–693, 1989. View at: Publisher Site  Google Scholar
 M. Bock, M. Andrassi, L. Belitsky, and B. Lorenz, “A comparison of two multifocal ERG systems,” Documenta Ophthalmologica, vol. 97, no. 2, pp. 157–178, 1998. View at: Publisher Site  Google Scholar
 M. Gur and Y. Zeevi, “Frequencydomain analysis of the human electroretinogram,” Journal of the Optical Society of America, vol. 70, no. 1, pp. 53–59, 1980. View at: Publisher Site  Google Scholar
 A. C. Fisher, R. P. Hagan, and M. C. Brown, “Automatic positioning of cursors in the transient pattern electroretinogram (PERG) with very poor SNR using an Expert System,” Documenta Ophthalmologica, vol. 115, no. 2, pp. 61–68, 2007. View at: Publisher Site  Google Scholar
 H. Jägle, J. Heine, and A. Kurtenbach, “L: Mcone ratio estimates of the outer and inner retina and its impact on sex differences in ERG amplitudes,” Documenta Ophthalmologica, vol. 113, no. 3, pp. 105–113, 2006. View at: Publisher Site  Google Scholar
 E. Rispoli, A. Iannaccone, and E. M. Vingolo, “Lownoise electroretinogram recording techniques in retinitis pigmentosa,” Documenta Ophthalmologica, vol. 88, no. 1, pp. 27–37, 1994. View at: Publisher Site  Google Scholar
 E. L. Berson, “Longterm visual prognoses in patients with retinitis pigmentosa: the Ludwig von Sallmann lecture,” Experimental Eye Research, vol. 85, no. 1, pp. 7–14, 2007. View at: Publisher Site  Google Scholar
 J. J. vos Hzn., M. A. J. M. Coemans, and J. F. W. Nuboer, “No evidence for polarization sensitivity in the pigeon,” The Journal of Experimental Biology, vol. 198, part 2, pp. 325–335, 1995. View at: Google Scholar
Copyright
Copyright © 2014 Mathieu Gauvin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.