Gobierno de la ciudad de Buenos Aires

Hospital Neuropsiquiátrico "Dr. José Tiburcio Borda"

Laboratorio de Investigaciones Electroneurobiológicas

y Revista

Electroneurobiología

ISSN: 0328-0446

 

 

Do’s and don’ts in Fourier analysis of steady-state potentials

Michael Bach1 and Thomas Meigen2

1Univ.-Augenklinik, Freiburg, Germany

2Univ.-Augenklinik, Würzburg, Germany


Electroneurobiología vol. 14 (1), pp. 67-70, Dec. 2005; URL http://electroneubio.secyt.gov.ar/Fourier_analysis.htm

Copyright © 2005 del autor / by the author. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's full citation and original URL (above).

 

Purpose 

Fourier analysis is a powerful tool in signal analysis that can be fruitfully applied to steady-state evoked potentials (flicker ERG, pattern ERG, VEP etc.). However, there are some inherent theoretical assumptions in the discrete Fourier transform (DFT) that are not necessarily fulfilled in real-world applications. Furthermore, software packages developed for general engineering may be ill-suited and may not fully exploit the information of steady-state recordings.



Table of contents

  1. Purpose (above)
  2. What is special?
  3. Overspill
  4. 256 – No!
  5. Trend (“sawtooth artifacts”)
  6. Windowing? – No!
  7. Number of stimuli per sweep
  8. Response & Noise
  9. Conclusions
  10. Literature


What is special in clinical electrophysiology? 

Stimuli are periodic, and there is control over their frequency and absolute phase. This a-priori information can increase the power of the response measures as compared to textbook examples on Fourier analysis.

What we will ignore here:

 

The discrete Fourier analysis/synthesis (DFT) formulas, which contain all the effects demonstrated below. The term “FFT” (“Fast Fourier Transform” describes a special rapid implementation to do the DFT. The “fast” aspect is a trivial methodical detail, eternalized by an “FFT” button on some commercialized equipment, and should be replaced by the generic term “DFT”, whether or not performed by an FFT. Furthermore, the FFT has a minor technical disadvantage (see below).

 


“Overspill” 

If the analysis interval comprises an exact integer number of stimulus periods (here 4), a sharp frequency will be found (top). However, if there is a non-integer number of stimulus periods (here 3.5), “overspill” into neighbouring frequencies will occur (center graph). Fourier analysis assumes periodic progression outside the analysis interval; the resulting vertex generates artifactal high frequencies (center and bottom graphs). As we know the stimulus frequency exactly, we can avoid this artifact. (When employing a raster-based CRT, the stimulus frequency has to be an integer divider of the frame frequency; Bach et al 1997.)

Thomas Heinrich recently spotted an error in this figure: For the lowest curve, I concatenated the wrong traces. Luckily, it doesn’t affect the conclusion. In the forthcoming paper [Bach & Meigen (in press) Do’s and don’ts in Fourier analysis of steady-state potentials. Doc Ophthalmol] this is corrected.

 


256 – No! 

The experimental parameters principally determine the number of data points per sweep. This will often lead to an arbitrary number (e.g. 40, 200), not necessarily a power of 2 (e.g., 64, 128, 256, ...). However, the fast algorithms for the DFT need such a number (FFT by Coley & Tokey 1965; FHT (Fast Hartley transform), O’Neill 1988). Textbooks suggest to fill up with zeros. This plays havoc with the spectrum of periodic data (center graph). An easy solution is to ‘resample’: Interpolation to the next higher 2n number of data points (bottom graph). This can be done with linear or non-linear interpolation. A very small amount of spurious frequencies will be generated that can be neglected against other noise sources.

As an alternative to FFT/FHT, one can use straight numerical integration procedures to determine the Fourier coefficients, these methods do not require power-of-two data points. Typical flicker-ERG, PERG or VEP sweeps will have less than 1000 data points, so a numerical integration to calculate the Fourier coeffiecients will only take milliseconds on current computers. To calculate the entire spectrum, however, the calculation time will be noticable.

 


Trend (“sawtooth artifacts”) 

Often records will show a “trend”, a slow upwards or downwards shift of the sweep (see “PERG, 16 rps” in the adjacent figure).

This can be caused by electrode drift or through the following mechanism: After eye blinks, there will be a huge offset on the signal, and the averager rejects the sweeps until the baseline enters its acceptance windowwith a speed determined by the time constant (or highpass setting) of the amplifier. Then averaging is resumed, but the baseline still drifts toward zero, causing the “trend”.

When a Fourier analysis is performed on such records, the trend will cause a typical “sawtooth spectrum”: The Fourier analysis assumes that the analysis interval repeats periodically and thus sees the trend as a sawthoot waveform with its typical spectrum (every harmonic present with diminishing amplitude).

 


Windowing? – No! 

Textbooks suggest to “window” the data. This makes sense because it reduces somewhat the influence of overspill due to a non-integer number of responses per analysis interval (see “overspill” above); it also reduces the “sawtooth” or trend artifacts (see previous paragraph).

The neighbouring figure shows that indeed many lines from the sawtooth spectrum are suppressed through windowing. But note also that the line at 1 Hz remains. Furthermore, new lines are created adjacent to the response (“overspill”).

 

Windowing generally leads to overspill near the response frequency. This is here more clearly shown: At the top there is a “clean” sinusoid with a corresponding single line in the spectrum. After windowing (bottom), adjacent lines (so called “overspill”, “leaking” or “sidebands”) occur. This is an unavoidable consequence of the modulation of the signal by windowing.

When (as in evoked responses) the stimulus frequency is known exactly, windowing is not necessary and trend artifacts can be otherwise dealt with (see above).

 


Number of stimuli per sweep 

A high number of stimulus periods per analysis interval allows to recognize “trend artifacts” that may arise from blink rejection and electrode drift. A linear trend will generate a sawtooth-spectrum (see above). If the response frequency is high as compared to the lowest possible frequency (representing the fundamental of the sawtooth), the response will not be markedly affected by sawtooth harmonics. If there is only one response period per sample interval, the entire low frequency noise will “fold” into the response frequency.

Increasing the sweep length thus has two advantages:

  • The response is increasingly disambiguated from drift artifacts
  • The frequency resolution increases.

A long sweep duration has one major disadvantage, though: If the seep needs to be rejected due to a blink artifact, too much valuable recording time is lost.

As a compromise between number of averages and sweep length (keeping total recording time constant) we suggest to choose about 8 responses per sweep.

 


Response & Noise 

Noise and true response (magnitude in the frequency domain) do not add linearly (Strasburger 1987, Norcia et al 1989).

Bad news: small amplitudes will likely be noise only. Good news: if the amplitude is ~3x the noise estimate, there is virtually no noise in that line.

 

The signal-to-noise ratio snr can be estimated either by the tCIRC statistic (Victor & Mast 1991), or by comparing the response frequency with the average of neighbouring frequencies (Meigen & Bach, ISCEV 1996).

The latter authors calculated analytically the effects and provide the following rule of thumb:

 snr    p  
2.00.12
2.80.05
3.00.04

 


Conclusions 

Discrete Fourier analysis (DFT) can greatly enhance analysis of our results. As we have much a-priori knowledge about the signal in evoked potentials, we can:

  • Choose an integer number of stimuli per sweep, avoiding overspill
  • Choose an unconstrained number of points per sweep and interpolate for fast HFT or FFT
  • Choose an integer relationship between stimulus period and frame rate
  • Avoid windowing
  • Disambiguate the response from trend artifacts
  • Make use of noise estimates
  • Estimate the significance of our responses
  • Use specific hypotheses relating the discrete frequency spectrum to the current physiological problem.

 


Literature 

  1. Bach M, Meigen T, Strasburger H (1997) Raster-scan cathode ray tubes for vision research – limits of resolution in space, time and intensity, and some solutions. Spatial Vision 10:403-414.
  2. Cooley JW, Tukey J (1965) An algorithm for the machine calculation of complex fourier series. Math Comput 19:297-301
  3. Meigen T & Bach M (1996) ISCEV poster [now in press: Meigen T & Bach M () On the statistical significance of electrophyisiological steady-state responses. Doc Ophthalmol ]
  4. Norcia AM, Tyler CW, Hamer RD, Wesemann W (1989) Measurement of spatial contrast sensitivity with the swept contrast VEP. Vision Res 29:627-637
  5. O’Neill MA (1988) Faster than fast Fourier. Byte 4:293-300
  6. Papoulis A (1984) Probability, random variables, and stochastic processes. McGraw-Hill, New York
  7. Strasburger H (1987) The analysis of steady state evoked potentials revisited. Clin Vision Sci 1:245-256
  8. Victor JD, Mast J (1991) A new statistic for steady-state evoked potentials. Electroenceph Clin Neurophysiol 78:378-388
  9. Winograd S (1978) On Computing the Discrete Fourier Transform. Mathematics of Computation Vol 48 Number 141:175-199

The present work, based on a poster at the ISCEV meeting 1997 in Asilomar, is expounded in more detail in

 

• Statistical response/noise estimate Meigen T & Bach M (1999) On the statistical significance of electrophysiological steady-state responses. Doc Ophthalmol 98:207-232

 

•Fourier stuff: Bach M, Meigen T (1999) Do’s and don’ts in Fourier analysis of steady-state potentials. Doc Ophthalmol 99:69-82

•Addition 2005-05-14: there is a nice demo of Fourier analysis/synthesis at Dan Russel's website.

_______________________________

 

 

 revista

Electroneurobiología

ISSN: 0328-0446

Croquis de este sitio - Outline of this site

Índice - Table des matières - Inhaltsverzeichnis - Table of Contents

Some downloadable articles and documents in this site are in PDF format (.pdf). To view or print these articles, you must have Adobe Acrobat Reader installed. Acrobat Reader can be downloaded from Adobe's web site free of charge by clicking here  /  Se recomienda mucho leer o imprimir algunos de los trabajos mas extensos con el programa Acrobat Reader,
que puede obtenerse gratuitamente pulsando aqui


Haga doble "click" en el título de cualquier artículo, para leerlo ahora - Double-click on any article to read it now:

 

 

 

SOCIOLOGÍA DE LAS NEUROCIENCIAS

 Pour comprendre l'enjeu :  L’anthropologie ganglionnaire, un psychovirus démasqué (français)

Puede leer, imprimir o guardar en su disco duro esta investigación en versión .PDF (190 kB: recomendada) o .DOC (76 kB).

 

 L'aliénisme en Argentine :  Diego Alcorta (1827) : Dissertation sur la manie... aiguë? (français)

Puede leer, imprimir o guardar en su disco duro esta investigación en versión .PDF (600 kB: recomendada) o .DOC (320 kB).

 

 

ELECTRONEUROBIOLOGÍA

 

Efectos relativísticos en biofísica cerebral:

  An ESSENTIAL preprint:  Effects of relativistic motions in the brain and their physiological relevance (To be published in Helmut Wautischer, ed., Ontology of Consciousness: A Modern Synthesis) (English) (Tema: funcionamiento del cerebro y psiquismo)

Puede leer, imprimir o guardar en su disco duro esta investigación en versión .PDF (496 kB: recomendada) o .DOC (227 kB)

MBYKYHÁPE GUARANÍME 

  SUMARIO Y PÁRRAFOS INICIALES EN CASTELLANO   

 SUMÁRIO EM PORTUGUÊS  

 ABSTRAKTI SUOMEKSI

 

 

Diversificación de recursos electroneurobiológicos en la evolución del sistema nervioso:

 

  Anticipo de impreso:  ¡Alma 'e reptil! Los contenidos mentales de los reptiles y su procedencia filética (Introducción a una próxima edición de Biomedicina de reptiles por Juan Carlos Troiano y colaboradores) (Castellano with a VERY LONG English abstract) (Tema: origen del sistema nervioso)

Puede leer, imprimir o guardar en su disco duro esta investigación en versión .PDF (733 kB: recomendada) o .DOC (406 kB)

SUMÁRIO EM LÍNGUA PORTUGUESA

  SUMARIO CASTELLANO

 

  Anticipo de impreso:   Electroencefalograma y cerebro en reptiles (Fragmentos del Capítulo 2 del libro mencionado, Biomedicina de reptiles por Juan Carlos Troiano y colaboradores) (Tema: el electroencefalograma de los reptiles invierte una relación básica en el electroencefalograma de los mamíferos: con frecuencia el de los reptiles pasa a presentar sincronías cuando están despiertos y a desincronizarse cuando están dormidos, mientras que en mamíferos es al revés. Así, esta sincronización del electroencefalograma reptil no concurre con la desconexión sensitiva del entorno, que encontramos en el sueño mamífero.) (Castellano)

Puede leer, imprimir o guardar en su disco duro esta investigación en versión .PDF (354 kB: recomendada) o .DOC (236 kB).

 

 

Cálculo de potenciales dentro de las células

 Calcule intensidades eléctricas y magnéticas en cada compartimiento neuronal:  The nervous principle: active versus passive electric processes in neurons (Explains how to calculate electric and magnetic field strengths inside different neuronal compartments) (LONG FILE IN ENGLISH with Bulgarian, Russian and Spanish abstracts/TOCs)

Podrá  leer, imprimir o guardar en su disco duro esta investigación en versión .PDF (2 Mb): recomendada) o .DOC (1,5 Mb). También como .html comprimido (compressed .HTML folder: 0,5 Mb) .ZIP.

ENGLISH ABSTRACT   

 AБСТРАКТ НА БЪЛГАРСКИ

SUMARIO CASTELLANO

РЕЗЮМЕ НА РОССИЙСКОМ ЯЗЫКЕ

 

 

NOCIONES GENERALES

Conceptos:

  Noticia general -- ¿Qué es electroneurobiología? -- La atmósfera intelectual (all in Spanish) -- Main Technical Ideas / Conceptos técnicos principales (English and Spanish) -- El descubrimiento de la Doppelrinde (German and Spanish)

 

 

Historia de las experimentaciones:

 1883: el primer mapeo con electricidad en todo el mundo ­¡durante ocho meses! de un cerebro humano consciente, realizado en San Nicolás, Provincia de Buenos Aires (Spanish)

Table of Contents (partial) of "Sensing: a new fundamental action of nature" (English) -- Índices

 

 

Recepción de los aportes de Chr. Jakob en la neurobiología germana

 

  De "La citoarquitectonía de la corteza cerebral":  Constantin von ECONOMO y Georg KOSKINAS (1925) - Kommentar über die Forschungsergebnissen von Prof. Christfried Jakob (Deutsch)

You can download a .PDF (700 kB: recommended) or .DOC (1.3 MB) file for printing this work

 

   The comments on Professor Christfried Jakob's contributions made in 'The Cytoarchitectonics of the Adult Human Cortex' by Professors Constantin, Baron von Economo, and Georg N. Koskinas, rendered in English by Prof. H. Lee Seldon (English)

You can download a .PDF (673 kB: recommended) or .DOC (1.29 MB) file for printing this work.

 

Los comentarios de Economo y Koskinas sobre las contribuciones de Christofredo Jakob, traducidos por M. Szirko (Castellano: todavía no disponible)

 

  Morphogenetic versus morphofunctional theory -- Franz J. Irsigler's intervention in the Behavioral and Brain Sciences' discussion on the implications of the "initial brain" concept for brain evolution in Cetacea (1988) (English)

Puede obtener un archivo .PDF (256 kB: recommended) o .DOC (106 kB) para imprimir este trabajo.

 

 

Recepción de los aportes de Chr. Jakob en la neurobiología y la sociedad argentina

 

  Piva y Virasoro - Christofredo Jakob, neurobiólogo: científico en diálogo filosófico (Spanish)

Puede obtener un archivo .PDF (391 kB: recommended) o .DOC (466 kB) para imprimir este trabajo.

 

 

  Kurowski - Animal ótico, o homem: aportes de Christofredo Jakob incoporados à "biopsicosociología" de Juan Cuatrecasas na ciência argentina (Portuguese)

Puede obtener un archivo .PDF (385 kB: recommended) o .DOC (813 kB) para imprimir este trabajo.

 

  Kurowski - Animal óptico, el hombre: aportes de Christofredo Jakob incorporados a la "biopsicosociología" de Juan Cuatrecasas en la ciencia argentina (Spanish)

Puede obtener un archivo .PDF (384 kB: recommended) o .DOC (814 kB) para imprimir este trabajo.

 

 

Panorama evolutivo:

 Table of Contents of "Sensing: a new fundamental action of nature" (English) -- Humour... -- Comment l’ hylozoïsme scientifique contemporain aborde-t-il la sélection naturelle du parenchyme neurocognitif? (français)

 

 

 Vom Tierhirn zum Menschenhirn - Die Organisation der grauen Substanz des menschlichen Zentralnervensystems (German and Spanish) -- Christofredo Jakob: su vida y obra (Spanish)

 

 

 

FILOSOFÍA DE LA CIENCIA - CONCEPTO DE TIEMPO EN NEUROBIOFÍSICA

  Why is Time Frame-dependent in Relativity? Minkowski's spacetime as a Kantian 'condition of possibility' for relativistic calculations  (English) -- Abstract: Minkowski spacetime is a condition of possibility for relativistic calculations. It keeps special relativity kinematic (i.e. avoids introducing it it causal dynamics) but in exchange forces relativistic observers to describe at once past, present, and future states along the length of the observed moving entities. This has been considered a proof that past and future components of real entities enjoy an unremitting mode of existence, because, if the world described by relativity were three-dimensional, the kinematic consequences of special relativity and the experiments confirming them would be impossible. The latter is acknowledged as exact but the assumed proof is viewed as incorrect because, such imposition being a Kantian condition of possibility, it cannot reveal anything about what Kant called noumenon, namely extramentality. ( Traducción castellana en preparación ).

You can download a .PDF (1 MB: recommended) or .DOC (571 kB) file for printing this work.

 

RESUMEN DIVULGATORIO Y PARA ESTUDIANTES EN CASTELLANO

Puede leer, imprimir o guardar en su disco duro el resumen en castellano en versión .PDF (204 kB: recomendada) o .DOC (96 kB)

 

 

 

“ANTAGONISMO ENTRE CIENCIAS DURAS Y HUMANIDADES BLANDAS”

 Hubert Markl: Dementia dichotoma: Der Zweikulturenwahn - La ilusión de las dos culturas  (Spanish and German) -- Vista del Sur: una nota a la American Association for the Advancement of Science, Washington (Presidenta Jane Lubchenco) sobre prioridades científicas (English)

 

 

MALFORMACIONES Y PAPEL DEL ÓRGANO CEREBRAL

 Christofredo Jakob:  Los Monstruos Anencéfalos”  (Spanish)

Puede leer, imprimir o guardar en su disco duro esta investigación en versión .PDF (346 kB: recomendada) o .DOC (280 kB).

 

 

 

BIOÉTICA

  Éthique de la Bio-Éthique  (français)  ¡Nuevo! 

Puede leer, imprimir o guardar en su disco duro esta investigación en versión .PDF (323 kB: recomendada) o .DOC (161 kB).

 

 

 

EL PRESUNTO DUALISMO CUERPO - ALMA

  Soigner et Guérir ? Le rapport soma-psyché et l’archéologie du divorce culturel entre Médecine et Philosophie (français) ¡Nuevo!

Puede leer, imprimir o guardar en su disco duro esta investigación en versión .PDF (373 kB: recomendada) o .DOC (201 kB)

SUMARIO EN CASTELLANO

 

 

 

PSICOANÁLISIS Y FACILITACIóN PSICOSOMÁTICA DE LA ENF. DE ALZHEIMER

 CELOS: Olvido, homosexualidad, desmemoria  [Temas: Celos y olvido. Los celos de los muertos. Celos y homosexualidad. Los significados de ‘Untergang’. La ‘destrucción’ del complejo de Edipo. Una hipótesis acerca de la desmemoria (destrucción de la memoria). La enfermedad de Alzheimer.] (Spanish)

Puede leer, imprimir o guardar en su disco duro esta investigación en versión .PDF (502 kB: recomendada) o .DOC (110 kB)

 

 

 

PERICIAS JUDICIALES Y CASUÍSTICA

 Pericia médica en un caso de lisencefalia [Temas: Derecho de incapaces. Fundamento científico de la tutela jurídica de la adjudicación de recursos para enriquecimiento experiencial (calidad de vida) de oligofrénicos con retardo profundo del desarrollo intelectual. Acuerdo de juspositivismo y jusnaturalismo al respecto.] (Spanish)

Puede leer, imprimir o guardar en su disco duro esta investigación en versión .PDF (300 kB: recomendada) o .DOC (150 kB)

 

 

 

NUESTRA GENTE

 Reseña biográfica: Ramón Carrillo, el Gran Sanitarista Argentino (Spanish)