Cortical Análisis Fuente de Alta Densidad EEG Grabaciones en niños

Summary

En los últimos años, ha habido un creciente interés en la estimación de las fuentes corticales de cuero cabelludo medidos actividad eléctrica para los experimentos de neurociencia cognitiva. En este artículo se describe cómo se adquiere EEG de alta densidad y cómo se procesan las grabaciones para la estimación de fuente cortical en los niños desde la edad de 2 años en el bebé Lab Londres.

Abstract

EEG se describe tradicionalmente como una técnica de neuroimagen de alta resolución espacial y temporal bajo. Los recientes avances en el modelado biofísico y procesamiento de señales hacen que sea posible la explotación de la información de otras técnicas de imagen como la resonancia magnética estructural que proporcionan una alta resolución espacial para superar esta restricción ^1. Esto es especialmente útil para las investigaciones que requieren alta resolución en lo temporal, así como el dominio espacial. Además, debido a la fácil aplicación y bajo costo de los registros de EEG, EEG suele ser el método de elección cuando se trabaja con poblaciones, como los niños pequeños, que no toleran IRM funcional analiza también. Sin embargo, con el fin de investigar lo que los sustratos neurales implicados, información anatómica de resonancia magnética estructural sigue siendo necesario. La mayoría de los paquetes de análisis de EEG trabajan con los modelos estándar de cabeza que se basan en la anatomía de adultos. La precisión de estos modelos cuando se utiliza para los niños está limitada ^2, debido a que el composition y la configuración espacial de la cabeza tissues cambios dramáticamente en desarrollo ^3.

En el presente trabajo, se proporciona un resumen de nuestro trabajo reciente en la utilización de modelos de cabeza sobre la base de imágenes por resonancia magnética estructural individuales o modelos de cabeza específicos de edad para reconstruir los generadores corticales de alta densidad de EEG. En este artículo se describe cómo se adquieren los registros de EEG, procesados ​​y analizados con la población pediátrica en el bebé Lab Londres, incluyendo la configuración de laboratorio, diseño de tareas, EEG preprocesamiento, procesamiento de resonancia magnética, y el nivel del canal EEG y análisis de origen.

Introduction

El presidente Barack Obama describió el cerebro humano como la próxima frontera de los descubrimientos científicos con una alta importancia para la salud y la economía de ³ (http://www.whitehouse.gov/share/brain-initiative). Sin embargo, como cualquier otro campo de las ciencias naturales, las ciencias neurológicas depende de los avances en las metodologías y técnicas de análisis para el progreso. Dos herramientas no invasivas de uso común en los estudios sobre la función cerebral en los seres humanos son la resonancia magnética (RM) y la electroencefalografía (EEG). Estos herramienta explotar diferentes propiedades físicas y proporcionar diferentes ideas sobre la función cerebral con ventajas y desventajas. MRI utiliza las propiedades magnéticas de las moléculas de agua dentro de los campos magnéticos para obtener imágenes de los tejidos vivos. El sujeto tiene que ser colocado en un imán con alta intensidad de campo. El movimiento de la participante está restringido durante estos procedimientos y el participante tiene que tolerar el ruido causado por los rápidos cambios en el campo magnéticocampo. Además de las imágenes estructurales, resonancia magnética también ofrece la posibilidad de medir los cambios en la oxigenación de la sangre para investigar la función cerebral (fMRI). En resumen, la RM ofrece relativamente alta resolución espacial de hasta 0,5 mm ³ con modernos altos campos escáneres y parámetros optimizados ^4. En contraste, la resolución temporal de resonancia magnética funcional se limita a la lenta cinética de la respuesta BOLD, que refleja sólo indirectamente la dinámica temporal altos de actividad neural ^5,6.

Por otro lado, el EEG mide los cambios en la actividad eléctrica causada por la actividad de las neuronas a través de electrodos colocados en el cuero cabelludo. Los recientes avances en la tecnología EEG permiten el uso rápido y fácil de los sensores de corto plazo o largo plazo y fijos, así como grabaciones ambulatorios. Debido a que el EEG es menos restrictiva, que también es el método de elección para ciertas poblaciones de los participantes que no toleran el medio ambiente así como la RM pediátrica y ciertageriatría y psiquiatría. Las propiedades de EEG muestran un patrón inverso a los de resonancia magnética: la resolución temporal es muy alta con una precisión de milisegundos, pero la resolución espacial es limitada. Las corrientes eléctricas pasan a través de diferentes tejidos entre su generador y los electrodos de EEG en la superficie del cuero cabelludo. Esto conduce a un mezclado y manchas espacial de actividad de la fuente conocida como el efecto de volumen de la conducción. Por lo tanto, la actividad medida por los electrodos en la superficie del cuero cabelludo refleja la actividad de múltiples fuentes que podrían ser distante a la posición del electrodo en la cabeza de ^1,7.

Gran parte del trabajo en los últimos años se ha dedicado a la fusión de la RM y EEG con el fin de aprovechar sus respectivos puntos fuertes. Una línea de trabajo se dedica a la adquisición simultánea de EEG y la RM en los estudios funcionales. Otro enfoque es usar la información espacial proporcionada por resonancia magnética estructural para tener en cuenta el volumen Conduction efecto a través del modelado biofísico. El uso de la información estructural para la reconstrucción fuente de registros de EEG es particularmente útil para estudios con una población pediátrica. La investigación sobre el desarrollo de la función cerebral es fundamental para la comprensión de cómo las habilidades cognitivas complejas se construyen en la parte superior de precursores simples ^8.

Estas investigaciones ayudan a poner de relieve los cambios en los sustratos neurales y propiedades de respuesta que se correlacionan con los cambios en el rendimiento de comportamiento. Sin embargo, la investigación de la función cerebral y la cognición durante el desarrollo también plantea retos específicos. En particular, la oportunidad para que los estudios de resonancia magnética funcional se limita como los niños pequeños y bebés, o bien tienen que estar dormido o sedado para obtener datos de MRI y sin artefactos de movimiento y el impacto negativo sobre el bienestar de los participantes. Además, el EEG se percibe como menos riesgoso e invasivo por los padres, lo que hace el reclutamiento de participantes en la investigación más fácil. Tor lo tanto, EEG es el método de elección para muchas investigaciones sobre la función del cerebro en los niños pequeños. Avances metodológicos en los sistemas de EEG permiten la aplicación de matrices de electrodos de alta densidad con 128 o más canales en cuestión de minutos. Facilidad de aplicación y la comodidad de uso son suficientes para permitir el registro EEG incluso en los niños más pequeños. Sin embargo, a menudo los investigadores no sólo están interesados ​​en la dinámica temporal de las respuestas a estímulos particulares, pero también le gustaría comparar los sustratos neurales que median las respuestas.

Una suposición que prevalece en el análisis de ERP de nivel de canal comparando diferentes grupos de edad es que los mismos sustratos neurales responden, pero que el momento o la amplitud de la respuesta varía entre las edades de ^9. Similar cuero cabelludo topografía se utiliza a menudo como un indicador de la actividad neuronal subyacente similar. Sin embargo, muchas diferentes configuraciones de origen pueden conducir a las topografías cuero cabelludo similares ^10. Mediante la aplicación de la estimación de la fuente, esta UNCErtainty se puede reducir y cuantificados. La independencia de las observaciones es crítico para las cuentas de la red de la función cerebral: si se mezclan las fuentes, las correlaciones estarán sesgados hacia una mayor conectividad local. Fuente de reconstrucción se puede aplicar para reducir este sesgo ^11. Alternativamente, las diferencias en tiempo y fase se pueden utilizar para el análisis de la conectividad, pero estos modelos matemáticos requieren supuestos que son difíciles de evaluar en datos no simulados ^12. En resumen, la estimación de fuente proporciona información adicional a nivel de los canales EEG y análisis ERP basado en el conocimiento sobre la anatomía y propiedades biofísicas de los tejidos.

Diferentes algoritmos se han ideado para encontrar soluciones al problema inverso. Estos algoritmos se dividen en dos categorías: paramétricas y no paramétricas ^13. Modelos paramétricos asumen una o múltiples dipolos que pueden variar en localización, orientación y fuerza. Por el contrario, no modelos paramétricos containa gran número de dipolos con la ubicación y orientación fija. En estos modelos, la actividad eléctrica del cuero cabelludo se explica como una combinación de activaciones en los dipolos fijos ^10,13,14. No paramétrica, modelos de fuentes distribuidas se pueden basar en el conocimiento sobre la anatomía y la conductividad en diferentes medios de comunicación. Modelos de Elementos de Contorno incorporan valores de conductividad de los principales tejidos de la cabeza con diferentes conchas para el cerebro, el líquido cefalorraquídeo y el cráneo. Esto se basa en la suposición de que la conductividad es principalmente constante dentro de cada compartimento, pero que los cambios marcados se produce en el límite de los diferentes compartimentos. Modelos de elementos finitos se basan en una mayor segmentación de los exámenes de IRM en la materia gris y el blanco de manera que los valores de conductividad se pueden asignar a cada voxel ^15.

En términos prácticos, los modelos no paramétricos son particularmente útiles para la reconstrucción de fuentes en las tareas cognitivas complejas, en las que el número de áreas involucradas puedenno ser conocido ^10. Modelos de elementos de contorno son más utilizados en la literatura actual, probablemente porque los más precisos modelos de elementos finitos representan comparativamente altas demandas computacionales. Además, existe una considerable variabilidad individual entre en la materia gris y el blanco de manera que los FEM debe basarse en imágenes por resonancia magnética individuales.

Modelos no paramétricos requieren un segundo paso para hacer coincidir la actividad medida el cuero cabelludo hasta las predicciones del modelo a seguir. Una vez más, los diferentes enfoques con diferentes ventajas y desventajas se han discutido en la literatura (véase Michel et al. 2004 para una visión general). Los algoritmos más utilizados se basan en la estimación de norma mínima (MNE), que coincide con la actividad medida cuero cabelludo para una distribución de corriente en el modelo hacia adelante con la intensidad total más bajo ^16. MNE está sesgada hacia las fuentes débiles y superficiales. Algoritmos EMN ponderados Profundidad tratan de reducir la desviación de superficies mediante la introducción de ponderaciónmatrices basadas en supuestos matemáticos ^10. El enfoque LORETA ampliamente utilizado se basa también en ponderada EMN, pero, además, reduce al mínimo el laplaciano de fuentes, que conduce a soluciones más suaves ^17,18. LORETA se ha encontrado que el mejor rendimiento para fuentes individuales en estudios de simulación ^19,20. Sin embargo, LORETA puede conducir a un exceso de suavizado de soluciones. Profundidad ponderada EMN es preferible cuando las fuentes son desconocidos o múltiples fuentes son propensos a estar presente ^{13, 16.} Se recomienda Comparando los resultados de diferentes algoritmos para evaluar la influencia de diferentes supuestos del modelo.

En resumen, la reconstrucción de fuentes a través de métodos de modelado se ha limitado a los niños hasta hace poco. Esto se debe a que la mayoría del software de análisis de EEG se basa en modelos de cabeza en base a la anatomía adulta que limita sustancialmente la precisión de las soluciones de código en los niños ^2,8. El acceso barato a la potencia de cálculo y la prestación desoftware fácil de usar para la reconstrucción de la fuente que sea posible superar estas limitaciones. La aplicación de la estimación de origen a la EEG proporciona dos ventajas importantes sobre el análisis sobre la base de observaciones del nivel del canal solas: una mejor resolución y la independencia de las observaciones espaciales.

Fuente de estimación puede no ser informativo en algunos casos: se requiere una buena cobertura de la cabeza para distinguir fuentes. Los sistemas de alta densidad con 128 o más electrodos se recomiendan ^10,15; una cobertura más escasa actuará como un filtro espacial que conduce a activación de la fuente propagación más amplia o resultados negativos falsos ^10. Además, sólo se ha informado de reconstrucción fuente basado en el método descrito en este artículo para generadores corticales. Por lo tanto, es menos adecuado para probar hipótesis acerca de sustratos subcorticales o interacciones subcorticales corticales. Por último, el análisis de origen debe basarse en hipótesis previas detalladas sobre los sustratos corticales,teniendo la literatura existente de otras modalidades de imágenes en cuenta. Técnicas de filtrado espacial también se pueden utilizar para mejorar la resolución espacial de la señal de EEG mediante la reducción de la mezcla espacial en el nivel del cuero cabelludo. Se utilizan métodos alternativos para reducir la influencia de los efectos del volumen de la conducción sin modelado de la cabeza, por ejemplo, el filtrado de Laplaciano ²¹ o análisis de la densidad Fuente de corriente ^22. Sin embargo, estos métodos no proporcionan más información acerca de los generadores neuronales como efectos volumen de la conducción son no sólo restringidos a los sensores en estrecha proximidad espacial ^1.

En las siguientes secciones, el artículo describe cómo los experimentos para la investigación de la función cerebral y cognitivo en los niños de 2 años de edad están diseñados en el bebé Lab Londres. A continuación, se analiza la adquisición de datos de EEG con los sistemas de baja impedancia de alta densidad con niños. A continuación, se presenta EEG preprocesamiento y análisis sobre el nivel del canal. Ultimoly, el artículo se centra en el tratamiento de los datos de resonancia magnética estructural para la reconstrucción fuente cortical y el análisis de las señales de nivel de fuente.

Protocol

Experimentos potenciales 1. Diseñar EEG y eventos relacionados para los niños

Nota: Un simple experimento fue diseñado para los fines de este artículo que se puede utilizar para investigar cara de procesamiento en los niños pequeños. La siguiente sección describe el experimento y explicar cómo implementarlo utilizando MATLAB R2012b y Psychtoolbox V3.0.11 ^23,24. Las imágenes tomadas desde el conjunto NimStim de la expresión facial emocional ²⁵ fueron utilizados para este ejemplo. Este conjunto de estímulo está disponible para fines de investigación a petición ( http://www.macbrain.org/resources.htm ).

1. Transfiera las fotografías RGB a escala de grises para reducir las diferencias entre los estímulos. Ver Tabla 1 Nota:. Estos comandos requieren la Imagen Processing Toolbox ( http://www.mathworks.co.uk/products/image/ ). Alternativas libres se pueden encontrar tediante el intercambio de archivos ( https://www.mathworks.co.uk/matlabcentral/fileexchange ).
2. Utilice el software de control experimental para poner en práctica el experimento con el tiempo exacto para disparadores de presentación de estímulos mediante una serie de comandos (consulte la Tabla 1 para un ejemplo).

2. Adquisición de Datos

1. Asegúrese de que el niño se sienta cómodo con el entorno de prueba. Permita que los niños más pequeños para sentarse en el regazo de su cuidador o en el asiento de un niño cómodo. Deje que el niño lo vea y sienta la red de sensores antes de aplicarlo a la cabeza del niño. Si hay una red adicional, haga que el padre también probar uno encima ni coloque uno en una muñeca o un peluche de peluche.
2. Mida la circunferencia máxima cabeza para seleccionar el tamaño neto correcto para el niño. Utilice una cinta de medición y manténgalo a la nasión. A continuación, mida alrededor de la cabeza alrededor de la circunferencia máxima (~ 1cm por encima del inión). Nota: mantener un registro de la circunferencia de la cabeza medida y la red de sensores utilizados para el análisis posterior ^26. Sirve de ayuda si la cabeza los padres también se mide para hacer que los niños más cómodo con la situación.
3. Identificar el vértice de la cabeza en la intersección de la distancia media entre nasión y inión y punto periauricular izquierda y derecha. Marque este punto con una pluma de China para asegurar que el canal de vértice se coloca correctamente cuando la aplicación de la red.
4. Aplique la red de sensores y asegúrese de que los canales principales están alineados con los puntos de referencia anatómicos (nasion, INION, vértices y mastoides izquierda / derecha). Nota: Para obtener resultados más precisos, la posición de los canales en la cabeza puede ser adquirido digitalmente con un equipo especial de la digitalización. Los investigadores que deseen adquirir la posición del sensor deben consultar en los manuales de hardware y software apropiados. Alternativamente, mapas de electrodos que asumen la colocación estándar de electrodos a lo largo anatómico LANdmarks se pueden utilizar. Estos mapas pueden ser ajustadas para envejecer modelos de cabeza adecuados como se describe en la sección de análisis.
5. Asegúrese de que los canales tienen un buen contacto con el cuero cabelludo por el posicionamiento de los sensores de forma individual; gire suavemente cada sensor de lado a lado para mover el pelo de la manera.
6. Medir ganancias de canal y las impedancias de canal. Haga clic en "Inicio" para iniciar la grabación en el software de grabación NetStation EEG y empezar la ganancia y la medición de la impedancia. Si la medición no se inicia automáticamente, utilice el "Amplificador Calibrar" y el botón "Medir Impedancias netos".
7. Compruebe el software de grabación para los canales con impedancias superiores a 50 kW que aparecerán en rojo. Aplicar la solución de electrolito adicional con una pipeta para bajar impedancias de canal. Compruebe la pantalla para los canales de EEG que muestran la actividad de alta frecuencia a pesar de baja impedancia o notablemente menos actividad que los canales que rodean (canales de línea plana). Estos canales pueden have contacto flojo con el cuero cabelludo y requerir un ajuste.
8. Con el fin de mantener a los niños cómodos durante la preparación EEG, permite al niño escuchar música, ver un dibujo apropiado para su edad o distraerlos utilizando otro experimentador, por ejemplo, soplando burbujas de jabón para los niños.

3. Análisis

1. Preprocesamiento
   1. Filtrar digitalmente los datos con un filtro de paso alto con un punto de corte a 0,1 Hz para eliminar el canal se desplaza ²⁷ (Tabla 1).
   2. Para el análisis de ERP, aplicar un filtro de paso bajo con un punto de corte a 30 Hz ²⁷ (Tabla 1).
   3. Época los datos continuos de acuerdo con los códigos de disparo establecidas durante la grabación. Para la mayoría de experimentos, utilizar una línea de base de 200 mseg antes de inicio del estímulo y un intervalo post estímulo de 600 mseg para cubrir el intervalo de tiempo de interés (Tabla 1).
   4. Retire épocas que contienen movimiento o parpadeo artefactos: marca channels con un pico a pico de amplitud superior a 150 mV - ajustar este umbral en función del grupo de participantes y la calidad de los datos. Para mantener la coherencia, utilice el mismo umbral para todos los participantes en un estudio. Si un canal está por encima de este umbral en más de 30% de las épocas, quitar el canal (la actividad del canal puede ser interpolada a partir de canales circundantes, si éstos contienen datos aceptables). Si más de 20% de los canales están marcados como no válidos en una época, eliminar la época. Si más de 20% de los canales son rechazados por el algoritmo o menos de 50% de épocas se conservan, considerar la eliminación del conjunto de datos del análisis adicional (Tabla 1).
      Nota: Los umbrales de porcentaje para la época y el rechazo de canal son cifras aproximadas que eliminan una cantidad suficiente de ruido en nuestra experiencia. La cantidad de artefacto en la grabación es probable que sea diferente el uso de otros grupos de participantes, paradigmas experimentales o sistemas de EEG. Los experimentadores pueden querer ajustar el porcentajeumbrales y comprobar si están satisfechos con el rechazo de artefactos. Alternativamente, los experimentadores pueden rechazar los ensayos que contienen artefacto mediante una inspección visual.
   5. Re referencia a la media de referencia restando la actividad Medio en todos los canales de cada canal (Tabla 1). Nota: El electrodo de vértice se utiliza típicamente como la referencia de grabación en NetStation.
2. Corrección de artefactos Usando Análisis de Componentes Independientes
   1. Importe los datos en la caja de herramientas más rápido ²⁸ y ejecutar el artefacto algoritmo de rechazo automático de los datos (Tabla 1).
   2. Utilice la interfaz gráfica de usuario (GUI) para más rápido; para abrir la interfaz gráfica de usuario, el tipo FASTER_GUI en la línea de comandos.
   3. Anule la selección de las opciones de filtrado en el menú de filtrado ya que los datos que ya se ha filtrado antes de la epoching.
   4. Especifique el número de canales: 126 canales de EEG con 2 electrooculagram (EOG) canales.
   5. Introduzca elmarcadores utilizados para epoching los datos como cadenas en una matriz de células. Para el caso que se presenta entre: {'cara', 'scra'} para la cara y las condiciones de cara revueltos.
   6. Seleccione los canales para el análisis de componentes independientes (ICA). Típicamente seleccionar todos los canales de registro, incl. canales de EEG no externos.
   7. Especificar la entrada y la carpeta de salida en el panel derecho de la interfaz gráfica de usuario.
   8. Seleccione el archivo de canal adecuado para las grabaciones. Nota: Los archivos de canal durante la mayor parte del sistema EEG se pueden descargar desde el fabricante, o se pueden descargar de la página web EEGLAB.
   9. Haga clic en Ejecutar para iniciar el procesamiento más rápido. Dependiendo de la duración de las grabaciones y el número de archivos, este tratamiento puede durar varias horas.
   10. Inspeccione visualmente las grabaciones, mapas independientes de componentes y sistemas ERP después del procesamiento.
3. Channel Level Análisis de Datos Potenciales relacionados con eventos
   1. Combine varios canales para formar un virtual cHannel con una mejor relación señal a ruido (Tabla 1). Nota: La selección de canales se debe basar en los informes anteriores en la literatura o hipótesis a priori. Selección de canales que muestran la mayor amplitud dentro de una ventana de tiempo determinado no se recomienda ^29.
   2. Obtener medidas como pico de amplitud, la media de latencia y amplitud pico para caracterizar la forma de onda y realizar pruebas estadísticas (Tabla 1).
4. Crear Modelos de Elementos de Contorno (BEM)
   1. Segmento de la resonancia magnética anatómica escanear con FreeSurfer. Nota: Para obtener los resultados más precisos, la base del modelo de los elementos de contorno a la resonancia magnética individuales para cada participante. Si esto no es posible, se deben utilizar las plantillas medias de resonancia magnética que coinciden con la edad de los participantes en la mayor medida posible. Por favor, tenga en cuenta que los BEM no se puede utilizar para los niños menores de 24 meses. Modelos de elementos de contorno suponen que cada shell (cerebro, el cráneo, la piel) se compone de una capa cerrada. Sin embargo, en yniños óvenes las fontanelas en el cráneo no están cerradas, lo que viola el supuesto de capa cerrada.
      1. Con el fin de instalar el software FreeSurfer, primero descargarlo de la página web FreeSurfer ( http://freesurfer.net/fswiki/DownloadAndInstall ). A continuación, configurar el entorno de shell que incluye un FreeSurfer; . para bashrc, incluir los siguientes comandos en el archivo bashrc.:
         1. Exportación FREESURFER_HOME = / Aplicaciones / freesurfer /
         2. Fuente $ FREESURFER_HOME / FreeSurferEnv.sh
            Nota: Estos comandos asumen que la carpeta FreeSurfer está en la carpeta Aplicaciones, en un sistema Unix. Hay más detalles sobre cómo configurar FreeSurfer con entornos alternativos de shell, por ejemplo, csh / tcsh, o sistemas operativos en el sitio web FreeSurfer ( http://freesurfer.net/fswiki/DownloadAndInstall ).
      2. A continuación, defina el Asunto direcciónhistoria, es decir, la carpeta que la salida se escribirá con el siguiente comando:
         1. exportación SUBJECTS_DIR = / ... / BEM /
            Nota: Los resultados se pueden escribir en cualquier carpeta del sistema.
      3. A continuación, cambie el directorio de trabajo a la carpeta que contiene el archivo de la RM para el Modelo de Elementos de Límites con el siguiente comando:
         1. cd / Users / joebathelt / Neurodevelopmental_MRI_database / Los niños / cerebro /
            Nota: Cualquier carpeta en el sistema puede ser especificado usando la sintaxis del comando cd. Este es un ejemplo de la estructura de archivos investigadores primaria.
      4. Por último, iniciar la reconstrucción utilizando los siguientes comandos:
         1. recon-all-i-<mri_file> subjid <subject_id>
         2. recon-todo-todo-subjid <subject_id>
            Nota: <mri_file> necesita ser reemplazado con el nombre del archivo de las imágenes por resonancia magnética deseados en el directorio actual. <subject_id> puede ser repmezclada con cualquier nombre. FreeSurfer creará una carpeta con este nombre en el directorio de temas. Dependiendo del sistema utilizado, los últimos comandos pueden requerir algún tiempo para funcionar.
   2. Compruebe la segmentación FreeSurfer para la segmentación incorrecta, por ejemplo. esferas superpuestas, anatómicamente improbables compartimientos etc. mediante la importación de los segmentos en BrainStorm y utilizar las herramientas de visualización en la interfaz gráfica de usuario:
      1. En BrainStorm, seleccione el panel de la anatomía. Importe el MRI segmentada haciendo clic derecho sobre el tema y seleccionando "Importar Anatomía carpeta". Asegúrese de que la carpeta con la salida FreeSurfer se selecciona. Inspeccione visualmente la segmentación haciendo clic derecho y seleccionando la opción "Display". Nota: Como alternativa, los comandos FreeSurfer pueden ser utilizados. Una descripción detallada se puede encontrar en el sitio web FreeSurfer: http://surfer.nmr.mgh.harvard.edu/fswiki/RecommendedRECONSTRUCCIÓN. Si se desea región de análisis de interés sobre la base de parcelación anatómica, la FreeSurfer funciones mris_ca_label , mri_annotation2label y mri_mergelabels se pueden utilizar. Consulte las publicaciones FreeSurfer y páginas de ayuda en línea para obtener información más detallada.
5. Estimar la actividad de origen en BrainStorm
   1. Comience BrainStorm escribiendo "lluvia de ideas" en la ventana de comandos.
   2. Crear un nuevo protocolo seleccionando Nuevo Protocolo en el menú Archivo.
   3. Añadir un nuevo subject al protocolo seleccionando Nuevo Tema en el menú Archivo.
   4. Importar datos de EEG para el participante con un clic derecho sobre el tema y seleccionando "Importar MEG / EEG".
   5. Importar un archivo de canal haciendo clic derecho y seleccionando "Archivo de canal Importar". Nota: El archivo de canal debe estar en consonancia con la RM para la reconstrucción de fuentes. BrainStorm utiliza un sistema de 4 puntos de referencia anatómicos que el usuario tiene que marcar en la resonancia magnética. Por favor, consulte los tutoriales BrainStorm para más información ( http://neuroimage.usc.edu/brainstorm/CoordinateSystems ). La posición estándar tal como se define por un archivo de canal para un sistema de EEG determinado o, idealmente, las posiciones de la cabeza que fueron digitalizados antes de la grabación de EEG se pueden utilizar.
   6. Compruebe que el BEM y los canales se alinean como se esperaba: Haga clic derecho sobre el archivo de canal para el tema y vaya a "registro MRI" y ", Hora ". Nota: Si las esferas dentro del modelo se superponen o si los canales están dentro del BEM, la reconstrucción de fuentes producirán resultados incorrectos. Ajuste la alineación con la opción "Editar" en el menú de "segmentación de IRM".
   7. Calcular la matriz de covarianza de ruido de la línea de base de cada época haciendo clic derecho sobre el participante y la elección de "Noise Covarianza Matrix" y "Cálculo de la grabación". Nota: Los autores de la caja de herramientas BrainStorm recomiendan el uso de una matriz diagonal ruido covarianza para grabaciones cortas (~ menos puntos de tiempo que los canales) y uno completo para grabaciones más largas. Por favor, consulte el tutorial BrainStorm Fuente Estimación para más información: http://neuroimage.usc.edu/brainstorm/Tutorials/TutSourceEstimation .
   8. Calcular el modelo de origen haciendo clic derecho sobre el objeto y seleccionando "Computmodelo de código e ".
   9. Calcular la solución inversa usando profundidad ponderada mínima Norma Estimación haciendo clic derecho sobre el objeto y seleccionando "Compute Fuente" y "estimación norma mínima". Nota: Otras opciones (dSpm, sLORETA) están disponibles. Cada opción tiene diferentes ventajas y desventajas. El algoritmo debe seleccionarse con base en consideraciones a priori y los informes previos en la literatura. Además, algunos algoritmos son mejores en la resolución de la activación focal en determinadas zonas, mientras que otros son más adecuados para la activación de amplia difusión. MNE se utilizó para este estudio basado en informes previos en la literatura ^16. Por razones de coherencia, el mismo algoritmo para la solución inversa se debe utilizar para todos los participantes en un estudio. Los investigadores también pueden desear comparar cómo conclusiones sólidas son para la aplicación de diferentes algoritmos de solución inversa.
   10. Repita la Sección 3 para todos los participantes en el estudio. Nota: Ninguno de utilizar la grainterfaz de procesamiento por lotes ficas o está para repetir los pasos de procesamiento para los participantes. Consulte la documentación para obtener instrucciones BrainStorm ( http://neuroimage.usc.edu/brainstorm/Tutorials/TutRawScript ).
   11. Actividad de la fuente promedio durante los ensayos por los participantes arrastrando las grabaciones al menú Proceso y seleccionando la opción "Normal" y "Por condición (media subject)".
   12. Contrasta la condición mediante la selección de los "Procesos de 2" y arrastrando cada condición en una sola ventana. A continuación, seleccione "Test" y "t de Student" o "la prueba t de Student (emparejado)", dependiendo del diseño. Para realizar comparaciones múltiples, amplitud establecido y los umbrales de la zona en la pantalla del mapa estadístico que resulta en el menú "Estadísticas". Nota: Como alternativa, los mapas de activación se pueden exportar a SPM ( http://www.fil.ion.ucl.ac.uk / spm /) para más análisis estadístico en profundidad ( http://neuroimage.usc.edu/brainstorm/ExportSpm12 HYPERLINK "http://neuroimage.usc.edu/brainstorm/ExportSpm12") .
   13. Calcular la respuesta relativa al evento para una región de interés. Para ROIs parcellation basado, cargue la parcelación FreeSurfer haciendo clic derecho sobre la superficie de la corteza en el menú de la anatomía y seleccione "etiquetas de importación". Navegar por el archivo correspondiente y cargarla. Ahora, seleccione un retorno de la inversión en el panel "Explorador" del menú de datos funcionales.
   14. Obtenga el retorno de la inversión relacionada con la actividad de eventos arrastrando los archivos a la ventana de proceso 1 y seleccione "Extraer Scouts Series de Tiempo" en el menú Fuentes. Nota: Varias regiones de interés se puede seleccionar de forma simultánea y las series de tiempo ROI puede ser exportada para su posterior trazado y análisis haciendo clic derecho en el momento del exploradordatos de series y seleccionando la opción "Exportar a Matlab".

Representative Results

El diseño de experimentos de potenciales evocados para lactantes y niños a menudo es un reto, debido a su limitada capacidad de tolerar los experimentos repetitivos largos ^30. Este problema se agrava aún más cuando el experimentador planea aplicar reconstrucción de fuentes, porque la reconstrucción fuente precisa requerirá una elevada relación señal-ruido ^1. Figura 1 muestra un protocolo experimental para la investigación de los mecanismos de procesamiento de la cara que se pueden utilizar con los niños muy pequeños. El paradigma está adaptado para a) reducir al mínimo parpadeo de los ojos y los movimientos oculares durante la presentación del estímulo, porque los niños van a ser menos capaces de controlar los movimientos de los ojos de los voluntarios adultos b) realizar el experimento más atractivo añadiendo acaparadores de la atención después de que el intervalo entre estímulos. Parpadeo de los ojos y los movimientos de los ojos son controlados mediante la presentación de una fijación cruz poco antes del estímulo con el fin de llamar la atención sobre el centro de la pantalla. Además, la duración del estímuloestá ajustado a 500 mseg, que permite la percepción consciente del estímulo y reducir al mínimo el tiempo para la digitalización de la imagen con los movimientos del ojo. El capturador de atención consiste en una imagen amigable para los niños se presenta con un sonido simultáneo. Una selección aleatoria de los diferentes estímulos capturador de atención puede ser usado para mantener la interesante tarea para el niño. La próxima prueba se puede iniciar por el experimentador, cuando es evidente desde el sistema de monitoreo que el niño está mirando el centro de la pantalla de nuevo. Además, historias pueden ser utilizados para ayudar a los niños mayores asisten a la pantalla. A menudo es útil para practicar las tareas con los niños antes de que el registro del EEG para asegurarse de que el niño entienda la tarea. Preguntas de selección o las puntuaciones obtenidas en la sesión de práctica se pueden utilizar como covariables en los análisis posteriores.

Al calcular el número de repeticiones necesarias para el experimento, es importante tener en cuenta que muchas pruebas se pueden perder debido a la falta de atención or movimiento artefacto cuando se trabaja con niños. Como regla general, el número necesario de repeticiones debería duplicarse en comparación con los estudios de adultos o contratar a un mayor número de participantes. La capacidad de atención y la cooperación es limitada en niños que en adultos. Por lo tanto, mantener a las necesidades específicas de los niños en mente al diseñar la tarea. Una tarea larga se puede dividir en varios bloques de tareas más cortos con descansos intermedios. Típicamente, el número de condiciones que pueden ser incluidos en el experimento es menor para los niños muy pequeños, ya que no serán capaces de cooperar para los períodos más largos necesarios para obtener ensayos suficientes para muchas condiciones de estímulo.

Las cifras que se presentan se basan en una grabación con un niño de 6 años (6 años 3 meses). El modelo de la cabeza se basa en una plantilla media de MRI de 6 años de edad ^31. Figura 5 muestra el nivel del canal de eventos relacionados con las respuestas potenciales (ERP) para enfrentar y revueltos cara estímulos. El waveform de la ERP a través de canales posteriores muestra el patrón esperado de una deflexión positiva seguida de una desviación negativa y una amplia desviación positiva posterior. Sobre la base de la topografía, curso de tiempo y la naturaleza del paradigma, estas desviaciones es probable que representan la P100, N1 y finales de componente potencial positivo. Además, la temprana desviación negativa es significativamente mayor para enfrentar los estímulos en comparación con caras revueltos. Por lo tanto, es probable que refleje la cara componente específico N170. Los mapas topográficos en la Figura 5 muestran la distribución de la tensión entre 250 y 300 mseg. Voltaje negativo, con un máximo a través de canales derecho occipito-temporal en la condición de las caras es evidente.

La figura 6 muestra la comparación estadística de actividad de la fuente proyectada en base a un modelo de la cabeza adulto estándar y un modelo de la cabeza adecuada para su edad. Fuente reconstrucción se basó en un modelo de elementos de contorno (BEM) con la profundidad ponderadaestimación mínima norma (wMNE) y el ruido total matriz de covarianza en Brainstorm v 3.1 ^32. El valor predeterminado MNI Colin27 BEM fue utilizado como el modelo adulto. Actividad Fuente fue promediado en el tiempo entre 250 y 300 ms en línea con respuestas específicas de cara en el nivel del canal.

El mapa muestra los resultados de una comparación t de Student entre las caras y revueltos condición caras corregido para comparaciones múltiples utilizando la tasa de falso descubrimiento (FDR). Los resultados muestran activación de la fuente significativamente más fuerte sobre el lóbulo temporal en las caras en comparación con la condición caras revueltos. La localización utilizando el modelo apropiado para su edad es más focal con fuertes diferencias en la superficie ventral de la corteza temporal. La localización basado en el modelo de la cabeza del adulto es más dispersa y muestra las diferencias de actividad de origen en el medial derecha y giro temporal superior, que es en su mayoría ausentes en el mapa basado en el modelo de la cabeza adecuada para su edad.

Figura 1. Ejemplo de una percepción de la cara experimentar adecuada en un amplio rango de edad El experimento consiste en la presentación visual de imágenes de rostros o caras revueltos. Los estímulos son físicamente idénticos, pero la disposición espacial se asignaron al azar en la condición revueltos. Cada ensayo se inicia con la presentación central de una cruz de fijación para minimizar los movimientos oculares durante la presentación del estímulo. La duración de la presentación cruzada de fijación se asignó al azar para evitar influencias de arrastre a través de múltiples repeticiones. El estímulo se presenta más de una duración de 500 mseg. La corta duración también a minimizar la posibilidad de movimientos de los ojos durante la ventana de presentación. Un estímulo atraiga la atención se presenta después de un intervalo de prueba, con una duración entre aleatorio entre 1 segundo y 2 segundos. El por locapturador de atención es especialmente útil para los participantes más jóvenes que no son propensos a asistir a muchas pruebas de material no atractivo en secuencia. El experimento puede comenzar la siguiente prueba, cuando el participante está mirando en la pantalla en respuesta a la llamada de atención.

Figura 2. Diagrama de flujo para el algoritmo de rechazo de umbral. El algoritmo compara el máximo de cada canal de EEG en cada época a un umbral establecido. Si un canal contiene una actividad máxima por encima del umbral, el canal se marca como mala. Si más de 20% de los canales son malos en una época dada, la época es rechazada. Después de rechazo época, la actividad máxima en cada canal se compara con el umbral de nuevo. Si un canal contiene la actividad por encima del umbral en más de 30% de todas las épocas, el canal se rechaza. Si más del 20% de los canales son rechazados por este procedimiento o menos del 50% de las épocas por condición se quedan después del rechazo época, el conjunto de datos debería ser excluido de mayor análisis. Haga clic aquí para ver una versión más grande de esta figura.

Figura 3. Corrección automática artefacto usando la caja de herramientas más rápido ^28. La figura muestra las opciones que necesitan ser cambiados para utilizar la caja de herramientas más rápido con la configuración y canalización de procesamiento descrito en este artículo:.. 1 Filtrado debe ser desactivado, debido a que el conjunto de datos que ya se ha filtrado 2 El número de canales debe ser ajustado. El sistema EEG se utiliza en este artículo tiene 126 canales con2 canales EOG. 3. Event-marcadores de tiempo de necesidad de ser especificado como una matriz de células de cadenas. 4 de bloqueo. La ventana de tiempo para la respuesta relacionados con el evento tiene que ser suministrado. Esto tiene que ser idéntica a la ventana utilizada en la etapa epoching anterior. 5. El usuario tiene que definir canales para el análisis de componentes independientes (ICA). En la mayoría de los casos, esto comprendería todos los canales de EEG y canales externos relevantes, como los canales oculares (EOG). 6. Los índices de los canales de los ojos también necesitan ser ajustados para el sistema de EEG utilizado. Para el sistema EEG descrito, estos serían los canales 125 y 128. Por favor, haga clic aquí para ver una versión más grande de esta figura.

Figura 4. análisis de código fuente en una lluvia de ideas ^32. 1. Después de importar el conjunto de datos de EEG y superficies FreeSurfer, el Modelo de Elementos de Contorno (MEC) se puede calcular mediante la selección de "modelo de la cabeza Compute" en el menú "Fuente". 2. La matriz de covarianza de ruido puede calcularse a partir de las grabaciones seleccionando "covarianza de ruido Compute". Si la grabación es lo suficientemente largo, es decir, más puntos de tiempo que los sensores, la matriz de covarianza total puede ser calculado, de otro modo se recomienda una matriz diagonal. 3. Después de calcular el modelo de la cabeza y la matriz de covarianza de ruido, es posible obtener la solución inversa. Diferentes algoritmos pueden ser utilizados. La profundidad mínima ponderada Norma Estimación (wMNE) algoritmo se utilizó para este artículo. 4. La evolución temporal de la actividad de la fuente en las regiones de interés (ROI) se puede extraer, seleccionando "Extraer Scouts Series de Tiempo" f esde el menú "Fuente". ROI de parcelación cortical automática en FreeSurfer fueron utilizados para este ejemplo. Haga clic aquí para ver una versión más grande de esta figura.

Relacionados con la Figura 5. Evento respuesta potenciales (ERP) para enfrentar y revueltos cara estímulos sobre derecho occipito canales temporales. El ERP muestra una deflexión más negativa entre 130 y 220 mseg después del inicio del estímulo en el lado derecho para las caras en comparación con caras revueltos. Estas propiedades están en consonancia con los informes anteriores sobre el componente N170 ^33. Por favor haga clic aquí para ver una versión más grande de esta figura.

ove_content "fo: keep-together.within-page =" always ">
Figura 6. Comparación de la localización de la fuente entre un modelo de la cabeza adulto por defecto y un modelo de la cabeza adecuada. La fila superior de la figura muestra el MNI límite adulta modelo de cabeza colin27 a la izquierda y un modelo de la cabeza apropiada para su edad sobre la base de una parcelación FreeSurfer de un promedio plantilla de la RM para los niños de 6 años de edad, a la derecha. También se presentan las ubicaciones de los lugares de electrodos coregistradas. La grabación fue obtenida de un niño de 6 años (6 años 3 meses). La segunda y la tercera fila muestran los resultados de una comparación estadística entre los mapas de activación de reconstrucción de fuentes sobre las EMN en las caras en comparación con la condición de rostros revueltos basada en una prueba t corregido para comparaciones múltiples utilizando Tasa de Falso Descubrimiento (FDR). El mapa de color ilustra la magnitud del efecto de color rojo que indica una mayor actividaden la condición de caras y azul que muestra una mayor actividad en la condición de rostros revueltos. Por favor, tenga en cuenta que no se aplicaron de amplitud o tamaño umbrales. La localización parece más focal en la superficie ventral del polo temporal utilizando el modelo de la cabeza adecuada para su edad en comparación con el adulto BEM. Haga clic aquí para ver una versión más grande de esta figura.

Figura 7. Fuente ERP de la circunvolución fusiforme derecha en respuesta a caras y trepó caras basado en la reconstrucción fuente de una grabación obtenida de un niño de 6 años (6 años 3 meses) con una edad BEM apropiado con MNE El ERP de código muestran una mayor deflexión negativa alrededor de 250 milisegundos después de la aparición de estímulo en la condición de caras en comparación con los scrambled se enfrenta a condiciones. Esta actividad es probable que refleje la contribución del giro fusiforme derecho en el componente N170 en la condición de las caras.

<td>

Entrada para el paso 1.1	Descripción
input_image = imread ('/ Users / some_user / images / example.jpeg');	% Lectura de la imagen
gray_image = rgb2gray (input_image);	% Transferencia de RGB a escala de grises
guardarcomo (gray_image, 'grey_image.tiff')	% Ahorra nueva imagen
Entrada para el paso 1.2	Descripción
Ejemplo de código:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Limpieza
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc	% variables de compensación de espacio de trabajo
borrar todos
cerca de todo	% Configuración de los drivers de audio al modo de baja latencia
InitializePsychSound ([0]);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Variables
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Input_folder = 'C: Documents and Settings usuario Lab ERP Mis documentos MATLAB Faces_Houses ';	% Definir el directorio de trabajo
netstation = 1;
	% Este conmuta la comunicación con el software de grabación NetStation EEG
	% Definir el número de ensayos
no_of_trials = 80;	% Unificar nombres de teclado para una fácil portabilidad entreen versiones de Unix y PC
KbName ('UnifyKeyNames')	% Definir una variable durante las pulsaciones de teclas de escape
ESCAPEKEY = KbName ('ESCAPE');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Cargando estímulos
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
face_directory = strcat ('C: Documents and Settings usuario Lab ERP Mis documentos MATLAB Faces_Houses Faces ');	% El suministro de la carpeta con las imágenes de los estímulos
files = dir (face_directory);	% Generando una lista con todos los estímulos y fotos
enfrenta = {files.name};
caras (01:02) = [];
capturadores = dir ('/ Usuarios / joebathelt / Dropbox / preterm_oddball / capturadores /');	</ Td>
Sujetadores = {grabbers.name};
capturadores (1:2) = [];
grabber_sounds = dir ('/ Usuarios / joebathelt / Dropbox / preterm_oddball / sonidos /');
grabber_sounds = {grabber_sounds.name};
grabber_sounds (01:02) = [];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% A partir del experimento
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
si netstation == 1; NetStation ('Connect', '194 .82.245.131 ', '55513')	% Sincronización con el software de grabación
NetStation ('Sincronizar', 10)	% De sincronización dentro de los 10 ms de exactitud
NetStation ('StartRecording')	% Inicia la grabación
final
% Configuración de la pantalla
screenNum = 0;	% Número de identificación de la pantalla de presentación
	% Definir una ventana para la presentación del estímulo
[Display.w, display.rect] = Pantalla ('OpenWindow', screenNum, 0);
[WPTR, rect] = Pantalla ('OpenWindow', screenNum);	% Ocultando el cursor del ratón
HideCursor;
negro = BlackIndex (WPTR);	% Estableciendo el fondo en negro
Pantalla ('fillRect', WPTR, negro);	% De fuente y tamaño de fuente para la definición de la pantalla de texto
Pantalla ('TextFont', display.w, 'Arial');
Pantalla ('TextSize', display.w, 32);
Pantalla ('TextStyle', display.w, 0);
para i = 1: no_of_trials
FlushEvents
randomizer = randi (2);
	% Seleccionando al azar una imagen capturador de atención
attention_grabber = imread (strcat (Input_folder, 'capturadores /', cell2mat (capturadores (randi (longitud (acaparadores))))));
[Sonidos, fs, nbits] = wavread (strcat (Input_folder, 'los sonidos /', cell2mat (grabber_sounds (randi (longitud (acaparadores))))));
audio_handle = PsychPortAudio ('Open', [], [], 0, fs, longitud (sonidos (1, :)));
sonidos = sonidos ';	% Con selección aleatoria de un sonido de llamada de atención
PsychPortAudio ('FillBuffer', audio_handle, sonidos);
	% La creación de la presentación del sonido
estímulo = strcat (face_directory, caras (randi (longitud (caras))));
estímulo = imread (cell2mat (estímulo));	% Este aleatoriza la presentación de caras y rostros revueltos
si randomizer == 2;	% De cargar el estímulo que fue seleccionada al azar de la lista de estímulos
estímulo = estímulo (randperm (longitud (estímulo (:, 1))), randperm (longitud (estímulo (1, :))));
final	% Revolver la matriz de estímulo, si el ensayo es un ensayo revueltos
stim = Pantalla ('MakeTexture', WPTR, doble (estímulo)); % Estímulo
	% Preparar el estímulo para la presentación
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Trial
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	% Este apartado sólo es relevante para el primer ensayo
si i == 1;
[Nx, ny, caja] = DrawFormattedText (WPTR, "Pulse cualquier tecla para iniciar", "centro", "centro", [255 255 255]);
Pantalla ('FrameRect', WPTR, 0, caja);
[VBLTimestamp StimulusOnsetTim] = Pantalla ('Flip', WPTR);	% En el primer juicio, el texto "Pulse cualquier tecla para iniciar" se presenta hasta que se pulse una tecla
KbWait ([], 2);
[Nx, ny, caja] = DrawFormattedText (display.w, '', 'centro', 'centro', 255); % Pantalla en blanco
Pantalla ('FrameRect', display.w, 0, caja);		[VBLTimestamp StimulusOnsetTim] = Pantalla ('Flip', WPTR);
final
% Fijación Cruz
[Nx, ny, BBOX] = DrawFormattedText (WPTR,, 'centro' '+', 'centro', 255);
Pantalla ('FrameRect', display.w, 0, BBOX);
[VBLTimestamp StimulusOnsetTim] = Pantalla ('Flip', WPTR);	% Preparar un signo más blanco en el centro de la pantalla como la cruz de fijación
si netstation == 1;	% Presentando la fijación cruz
NetStation ('Evento', 'fix +', StimulusOnsetTim) final
	% De enviar una señal con la marca de tiempo de la presentación cruzada de la fijación y el código "fix +" para el software de grabación de EEG
% Limitar el tiempo de presentación de la cruz de fijación a una duración aleatoria entre 0,15 y 0,2 seg
% de Estímulo	% Presentando la cara o revueltos cara estímulo
Pantalla ('DrawTexture', WPTR, stim);
[VBLTimestamp StimulusOnsetTim] = Pantalla (WPTR, 'Flip');
si netstation == 1;	% El envío de gatillo con indicación de la hora y el código "cara" o la cara 'scra' al software de grabación de EEG
Si aleatorizador == 1;
NetStation ('Event', 'cara', StimulusOnsetTim)
elseif randomizer == 2;
NetStation ('Evento', 'scra', StimulusOnsetTim)
final
final
WaitSecs (0.5)	% Límite la duración del estímulo a 500 ms
	% Presentando una pantalla en blanco
[Nx, ny, caja] = DrawFormattedText (WPTR, '', 'centro', 'centro', 255); % Pantalla en blanco
Pantalla ('FrameRect', WPTR, 0, caja);
[VBLTimestamp StimulusOnsetTim] = Pantalla ('Flip', WPTR);
WaitSecs (1 + rand (1))	% Intervalo de ensayo con una duración entre aleatorio entre 1 y 2 segundos
% Capturador de atención	% Presentará la llamada de atención con sonido
Pantalla ('DrawTexture', WPTR, capturador);
sonidos = PsychPortAudio ('Start', audio_handle, 15, 0, 1);
[VBLTimestamp StimulusOnsetTim] = Pantalla ('Flip';, WPTR);
si netstation == 1;
NetStation ('Evento', 'grbr', StimulusOnsetTim)	% Enviará la indicación de la hora y el código para el capturador de atención al software de grabación de EEG
final
	% Presentará la llamada de atención y sonido hasta que se pulse una tecla del teclado
KbWait ([], 2);
PsychPortAudio ('Stop', audio_handle);
Pantalla ('FrameRect', display.w, 0, caja);	% Abortar el experimento, si se ha pulsado la tecla de escape
[VBLTimestamp StimulusOnsetTim] = Pantalla ('Flip', WPTR);
[KeyIsDown, timeSecs, keyCode] = KbCheck; si keyIsDown
si keyCode (ESCAPEKEY)
sca
NetStation ('stoprecording')
retorno
final
final	% Cerrar el experimento y detener la grabación de EEG en el final del experimento
final
Pantalla ('CloseAll');
si netstation == 1;
NetStation ('stoprecording')
final
Entrada para el paso 3.1.1	Descripción
OUTEEG = pop_eegfilt (INEEG, 0.1, []);	% = OUTEEG EEG datos después de la filtración, es decir. salida de la función
	% = INEEG EEG de datos antes de la filtración, es decir, entrada de la función
	0.1%: de paso alto frecuencia de corte
	% []: De paso bajo de corte, sin definir porque se desea un filtro de paso alto
Entrada para el paso 3.1.2	Descripción
OUTEEG = pop_eegfilt (INEEG, [], 30);	% []: De paso alto corta, sin definir porque se desea un filtro de paso bajo
	30%: de paso bajo de frecuencia de corte
Entrada para el paso 3.1.3	Descripción
OUTEEG = pop_epoch (INEEG, "evento", {'cara', 'scra'}, [-0,2 0,6]);	% Epoching
	% 'Evento', {'cara', 'scra'}: la función se le dice que utilice los eventos de activación 'cara' y 'scra' como el tiempo de bloqueo de marcadores. Estos disparadores se definen en el guión experimento para marcar el inicio de la cara y revueltos sti carainicio presentación mulus.
	% [-0,2 0,6] = ventana de tiempo para el ERP de 0,2 segundos antes de que el evento de tiempo de bloqueo de 0.6 seg
OUTEEG = pop_rmbase (EEG, [-0,2 0]);	% Eliminación de la línea de base
	% [-0,2 0]: ventana de tiempo de línea de base, es decir, 0,2 segundos antes de que el evento de tiempo de bloqueo a la hora del evento de cierre
Entrada para el paso 3.1.4	Descripción
función [EEG] = threshold_rejection (EEG, umbral)	Definición de la función%, la función requiere una estructura de conjunto de datos de EEG y un umbral en mV
para j = 1:2
para i = 1: longitud (EEG.data (1,1, :))	% Este bucles pasa por todas las épocas en un conjunto de datos EEG dado
data = EEG.data (:,:, i);
datos = datos - media (datos, 2);	% Restando ªe actividad media para evitar la influencia de los cambios de amplitud
maxima = max (abs (datos '))';	% Identificar la actividad máxima absoluta en todos los canales
bad_channels = máximos> umbral;
channel_rejection (:, i) = bad_channels;
si suma (bad_channels)> 0.2 * 128	% Si más de 20% de los canales están por encima del umbral, la época está marcada para el rechazo
epoch_rejection (i) = 1;
más
epoch_rejection (i) = 0;
final
final
si j == 1;	% Seleccionando sólo las épocas que no están marcados como mala
EEG = pop_select (EEG, "juicio", de encontrar (epoch_rejection == 0))
más	Marcar canales% que son malas en más de 20% de épocas después del rechazo época para rechazo de canal
bad_channels = media (channel_rejection, 2)> 0,2;
EEG = pop_interp (EEG, encontrar (bad_channels == 1), 'esférica');	% Aplique interpolación esférica a malas canales
final
final
Entrada para el paso 3.1.5	Descripción
EEG = pop_reref (EEG, []);	% Calcula la media de referencia
Entrada para el paso 3.2.1	Descripción
Canales = {'E84', 'E89', 'E90', 'E91', 'E94', 'E95', 'E96'}; % N170 derecho	% Canales occipito-temporal del hemisferio derecho% para N170 respuestas utilizando las etiquetas de canal para el% 128 canales Geodésica Hydrocel Sensor Net
EEG = pop_select (EEG, "canal", los canales);	% Selección de canales
Virtual_channel = media (EEG.data, 1);	% La combinación de los canales individuales a un canal virtual
Entrada para el paso 3.3.2	Descripción
N170_peak = max (abs ((averaged_ERPs (0,13 * 0,2 * srate srate: 0.2 * 0.2 * srate srate), [], 2));	% De la amplitud máxima dentro de una ventana de latencia para cada participante en mV
N170_peaklatency = 1000 * (Encontrar (averaged_ERP (0,13 * 0,2 * srate srate: 0.2 * 0.2 * srate srate) == N170_peak) + 0,2 * EEG.srate + 0,13 EEG.srate). / EEG.srate	% Latencia pico en ms
N170_mean = medias (averaged_ERPs:, 0,13 * 0,2 * srate srate: 0.2 * 0.2 * srate SRATE), [], 2);	% Media de amplitud en mV

Tabla 1. Comandos de MATLAB para implementar ªe ejemplo experimento y analizar los registros de EEG de alta densidad en el canal y nivelado fuente. La tabla resume el código para implementar las caras vs revueltos caras experimento de ejemplo. Además, se presenta el código para el pre procesamiento de la EEG en bruto. Además, se muestran los métodos para la extracción de características de forma de onda para el análisis de nivel de canal de la respuesta relativa al evento.

Discussion

El presente artículo describe el registro y análisis de alta densidad de EEG para la reconstrucción de los generadores corticales utilizando modelos de elementos de frontera por razones de edad plantillas MRI medios apropiados y profundidad ponderada estimación mínima norma en un paradigma ERP estándar adecuado para los niños. En este paradigma, se presentan imágenes de rostros y rostros revueltos. Diferentes autores utilizaron este paradigma para investigar el desarrollo de mecanismos de procesamiento de la cara sobre el desarrollo ^35. En el nivel de canal, las deflexiones más negativos sobre derecho occipito canales temporales se describen para la condición de cara a la condición de cara revueltos. La topografía, la latencia y la característica de respuesta son consistentes con el componente N170 ^34. Fuente anterior y las investigaciones EEG fMRI simultáneas informan de que el giro fusiforme es un generador de probabilidad de la respuesta N170. Los resultados de los análisis actuales muestran que la inversión de origen con un límite ponderado profundidadelemento de modelo (BEM) se puede utilizar para localizar actividad de la fuente en el giro fusiforme en la cara revueltos cara paradigma en el nivel de los participantes individuales. El uso de modelos de cabeza en base a la anatomía del participante individual o el uso de apropiado para la edad promedio de las exploraciones anatómicas para estudios de desarrollo, en los que la anatomía individual no está disponible, permitirá a la localización de la fuente más exacta ^2. Además, las regiones de interés pueden ser identificados sobre la base de los conocimientos anatómicos o algoritmos parcellation automáticas para investigar la respuesta relacionada con eventos de todo, las regiones corticales.

Existen varias limitaciones a la reconstrucción de origen, sobre todo en muestras de desarrollo, en el momento. En primer lugar, la reconstrucción de origen sobre la base de las plantillas medias de diferentes grupos de edad se supone que el individuo muestra el desarrollo típico del cerebro para su edad cronológica, que puede no ser necesariamente el caso, sobre todo en los grupos de pacientes. Paraejemplo, diversos estudios describen trayectorias atípicas en el crecimiento del cerebro de los niños que nacen prematuros o ³⁶ niños con autismo ^37. Es difícil estimar cómo estas diferencias anatómicas influirá en la precisión de la solución inversa y sesgar los resultados de las comparaciones entre los grupos de control típicos y atípicos.

En segundo lugar, los modelos directos como el modelo de los elementos de contorno (BEM) no incorporan inhomogeneidades conductividad dentro de los compartimentos, por ejemplo, las diferencias entre materia gris y blanca. Por consiguiente, la precisión para las fuentes subcorticales está limitada. Soluciones Fuente estaban restringidos a las fuentes corticales por esa razón. Modelos de elementos finitos se pueden aplicar para la resolución más precisa de los generadores subcorticales. Con las soluciones restringidas a la corteza, es importante tener en cuenta que la activación en las regiones corticales puede reflejar mecanismos subyacentes subcorticales causantes, por ejemplo, la comunicación a través de lazos de retroalimentación del tálamo.Por lo tanto, las inferencias causales sobre la participación de las regiones corticales se limitan a menos que se utilizaron modelos más complejos que se encuentran actualmente disponibles sólo para adultos anatomía normal, por ejemplo, la dinámica causal Modelado ^38,39.

Además, los modelos de elementos de contorno asumen capas completas para cada compartimiento. Sin embargo, los niños pequeños tienen puntos blandos en el cráneo, donde las suturas entre los huesos del cráneo no se fusionaron plenamente ^15. Esta violación de supuestos BEM limita severamente la aplicabilidad de la reconstrucción con fuente de BEM en niños menores de 2 años de edad. Modelos de elementos finitos se pueden utilizar para la reconstrucción de origen en este rango de edad.

En tercer lugar, a pesar de que se utilizaron modelos de cabeza apropiadas para la edad para la reconstrucción de la fuente, se utilizaron los valores de conductividad basado en muestras de adultos para modelar la conductividad dentro de cada compartimento. Sin embargo, la conductividad del tejido es probable que cambiar con el desarrollo, por ejemplo., A través deaumentos en la densidad del hueso ^15. Los valores de conductividad para los tipos de tejidos utilizados en BEM para los bebés humanos y los niños no están actualmente disponibles en nuestro conocimiento.

Materials

Name	Company	Catalog Number	Comments
High-density EEG sensor net (128 or 256 channels)	HydroCel Geodesic Sensor Net 128	Electrical Geodesic Inc., Oregon, US
EEG high impendance amplifier	NetAmps 200	Electrical Geodesic Inc., Oregon, US
Data Acquisition Computer	PowerMac G4	Apple Inc, California, US
Stimulus Presentation Computer	Optiplex 745	Dell Computers Inc., Texas, US
Stimulus Presentation Software	MATLAB R2012b with PsychToolBox	Brainard et al. 1997
EEG recording software	NetStation 4.5.1	Electrical Geodesic Inc., Oregon, US
EEG analysis software	MATLAB R2012b	The Mathworks Inc.
	EEGLAB	Delorme et al. 2004
	BrainStorm	Sylvain et al. 2001
MRI processing software	FreeSurfer	Fischl et al. 2004
	OpenMEEG	Gramfort et. al. 2010
Delorme, A., & Makeig, S. EEGLAB: an open source toolbox for analysis of single trial EEG dynamics including independent component analysis. Journal of Neuroscience Methods, 134(1), 9–21, (2004).
Sylvain, B., John, C., Dimitrios, P., & Richard, M. Brainstorm: A User Friendly Application for MEG/EEG Analysis. Computational Intelligence and Neuroscience, 1–13, (2011).
Fischl, B. et al. Automatically parcellating the human cerebral cortex. Cerebral Cortex, 14(1), 11–22, (2004).
Gramfort, A., Papadopoulo, T., Olivi, E., & Clerc, M. OpenMEEG: opensource software for quasistatic bioelectromagnetics. BioMedical Engineering OnLine, 9(1), 45, (2010).
Brainard, D. H. The psychophysics toolbox. Spatial vision, (1997).