Discos EEG 3D

Los modelos de los discos 3D son generados a partir de los datos procesados de los canales EEG por persona y por sesión. Se desarrollo una aplicación que  utiliza librería ModelbuilderMk2 sobre processing. La señales de EEG se encuentran por canales y cada canal es convertido en una de las capas verticales del modelo, cada capa tiene aproximadamente 250 mm de alto x 140 mm de ancho.

Screen Shot 2016-03-30 at 2.56.29 PM

 

Screen Shot 2016-03-30 at 2.59.46 PM 1 Screen Shot 2016-03-30 at 2.59.55 PM

Screen Shot 2016-03-30 at 2.58.52 PM Screen Shot 2016-03-30 at 2.59.06 PM

 

Diseño Tornamesa – Primera parte

Para el diseño de la tornamesa estamos utilizando un modelo basado en el sistema Lazy Susan de oomlout.co.uk. El sistema funciona con un motor stepper que va conectado a un engrane que mueve la parte superior. La rotación entre las capas se hace por medio de baleros.

Download (PDF, 63KB)

Parte superior con dos engranes, uno para el motor y otro para la base del motor

Download (PDF, 55KB)

Base con estrutura para montaje de motor, contenedores para valeros y base el lazy susan.

 

Sensor laser – primera parte

Para la lectura de los discos físicos estamos trabajando en el desarrollo de un sensor laser de una dimensión. Una array de 5 lasers con fotoresistencias que nos permitirán leer la topología de los discos 3D.  En el esquemático se muestran las conexiones generales de uno solo de los módulos. Aun estamos decidiendo si se trabajara sobre módulos individuales con una sola tarjeta con todos los conectores.

Screen Shot 2016-03-30 at 2.39.53 PM

El primer código nos permite leer los valores en average de todos los sensores funcionando como una sola señal. En este caso estamos trabajando con el microcontrolador Teensy 3.2.

 

 

Screen Shot 2016-03-30 at 2.40.38 PM

 

Screen Shot 2016-03-30 at 2.41.17 PM

El segundo código nos permite leer los laser independientemente así como controlar de forma global la intensidad del laser utilizando PWM.

 

Visualización de bandas EEG – Primer Grupo.

Para comenzar con el diseño de los discos EEG vimos la necesidad de primero conocer las características de las señales por cada una de las bandas que utilizamos para el registro. En este caso estamos trabajando con la zona frontal en los electrodos: AF3,F7,F3,F4,F8,AF4.

Las siguientes imagenes muestran las caracteristicas de señal de cada una de estas bandas por persona.  También se agregaron algunos filtros a las señales para mostrar las señales normalizadas.

Screen Shot 2016-03-30 at 12.34.13 PM

Screen Shot 2016-03-30 at 12.34.59 PM

Screen Shot 2016-03-30 at 12.35.10 PM

Screen Shot 2016-03-30 at 12.34.25 PM

Screen Shot 2016-03-30 at 12.34.35 PM

Screen Shot 2016-03-30 at 12.34.46 PM

El código completo puede encontrarse en el repositorio en el siguiente link : https://github.com/Lessnullvoid/ProyectoLiquido/tree/master/Processing/Galleta_Test

 

Reporte 01 – Análisis de datos del muestreo 01

image20

El objetivo del análisis es identificar patrones de comportamiento en la actividad cerebral asociados con manifestaciones del condicionamiento emocional/cultural al que México ha sido sometido históricamente.

La actividad cerebral es monitoreada mediante un dispositivo de EEG en sesiones donde se induce a cada sujeto a la meditación profunda enfocada en su práctica particular. Se espera que a partir de estas muestras se puedan identificar y caracterizar los patrones de actividad cerebral que están asociados a los estados emocionales específicos por los que atraviesa el sujeto en su meditación.

El muestreo es realizado siguiendo el procedimiento descrito en el documento “Guía de estímulos para experimento EEG”.

El método para el análisis está basado en los experimentos reportados en [1]. Aunque es necesario puntualizar las diferencias respecto a los objetivos y el método allí reportado junto con las condiciones que han sido tomadas en consideración en el análisis aquí descrito.

La conclusión de esta fase del análisis indica que existen y es posible caracterizar los patrones de actividad cerebral en las muestras tomadas, sin embargo son específicos para cada individuo y es difícil realizar una validación entre los distintos sujetos. Se propone la experimentación con métodos más sofisticados para obtener resultados más concluyentes al respecto.

1. Método de referencia

En [1] se aborda el problema de clasificar patrones de actividad cerebral asociados a emociones positivas, neutras y negativas, representando las muestras en un espacio de atributos basados en métricas de correlación estadística entre señales. Específicamente usan correlación, coherencia y un índice de sincronización de fase para cada banda de frecuencia entre cada uno los 171 pares de electrodos, y entre ellos seleccionan los atributos que ofrecen mejores resultados. Con ello muestran que la representación basada en patrones de conectividad ofrece una plataforma que se puede usar para clasificar (con un clasificador QDA) tres estados de valencia emocional: positivo, neutro y negativo.

Para el desarrollo de sus experimentos hacen las siguientes consideraciones:

  1. Muestreo
  • 40 participantes, no drogas
  • Exposición a 6 clips de video (2+, 2=, 2-)
  • Escáner de 64 canales
  • Medición de electrooculografía como señal de control
  • Bioamplificador 0-100 Hz, Conversión Analógica/Digital a 500Hz, Resolución de 12 bits
  • Indican con una tecla cuando experimentan un cambio emocional (eventos)
  1. Preprocesamiento de la señal
  • 16.386 segundos por segmento (8192 puntos a 500 Hz)
  • Filtrado con Pasabajas a 50Hz para eliminar la señal de AC (60Hz)
  • Remover componente de movimiento ocular
  • Remover artefactos, mismos parámetros independientemente de la emoción
  • Filtro Pasabanda Butterworth de orden 2 para acotar cada banda de frecuencia
  • Delta (0-4Hz), Theta (4-7 Hz), Alpha (8-12 Hz), Beta (13-30 Hz), Gamma (30-50Hz)

Después de caracterizar cada muestra (una muestra es el registro de EEG para un individuo y una emoción), la representan como un vector donde están contenidos, para cada de las cinco bandas de frecuencia: los índices de correlación para cada par de electrodos, los índices de coherencia para cada par de electrodos y los índices de sincronización de fase para cada par de electrodos. Esta representación resulta en un gran número de atributos para cada muestra por lo que seleccionan aquellos pares de electrodos que mejor caracterizan cada categoría. Estos son los patrones de conectividad que identifican y dan nombre al método. Luego prueban su representación en la tarea de entrenar un clasificador QDA con la mitad de los datos y usarlo para identificar la categoría de cada uno de los restantes, repitiendo el procedimiento con los conjuntos intercambiados (2Fold Cross Validation).

Su conclusión es que dado que la clasificación funciona, entonces tal representación basada en la conectividad entre canales es útil para caracterizar estados emocionales mediante dichos patrones de conectividad.

A diferencia del escenario reportado en [1], en este experimento se cuenta con un menor número de muestras. Aunque en éstas se puede identificar la entrada a la fase de meditación, no es posible determinar si ésta se refleja como un estado emocional positivo, neutro o negativo. Dado este desconocimiento tampoco es posible realizar la validación de la representación en una tarea como la clasificación ni entre distintos sujetos por lo que el método usando para analizar las muestras difiere del reportado. Otras diferencias se encuentran en el número de canales usados para registrar la EEG así como en la frecuencia de muestreo.

2. Método propuesto

A diferencia del método revisado, en nuestra tarea desconocemos si los patrones que representan el estado de meditación tienen características que se mantengan entre distintos sujetos por lo que no se puede recurrir a una tarea supervisada como la clasificación para validar los atributos propuestos.

Se propone entonces el uso de un método no supervisado para identificar las componentes en el estado de meditación de cada individuo, representándolo como un mapa de actividad conjunta para los distintos canales, que en otros términos es un patrón de conectividad. Esto restringe los resultados y no permite la comparación entre distintos sujetos pero proporciona evidencia del comportamiento de las señales y de cómo en ellas es posible caracterizar un estado emocional. Otros experimentos para distinguir entre un estado de meditación y cualquier otro deben ser realizados si se desea comprobar la hipótesis en ese sentido.

Los datos recolectados durante el primer muestreo tienen las siguientes características:

  • Un registro por sujeto con las características descritas en la “guía de estímulos para…”. De donde se observa que las muestras pueden dividirse en 3 tipos de regiones: artefactos, estado de no-meditación y estado de meditación. Este último corresponde a los últimos aproximadamente 5 minutos de la muestra, donde el comportamiento es en la mayoría de casos más estacionario.
  • Frecuencia de muestreo entre 50  y 51 Hz. Esto limita el ancho de banda que cada señal puede representar (por el teorema de Nyquist) a una frecuencia máxima de 25 Hz, por lo que la banda Gamma queda completamente descartada y la banda Beta se reporta incompleta

El método que se propone para identificar los patrones de actividad es como sigue:

  1. Identificar manualmente la región de cada muestra que representa el estado de meditación
  2. Filtrar la muestra para restringir el ancho de banda a frecuencias por debajo de 25 Hz
  3. Transformar la señal al dominio de la frecuencia para identificar los componentes de cada banda.
  4. Aplicar un algoritmo para identificar los componentes independientes principales que tienen un mayor aporte energético en cada banda de frecuencia.
  5. Aplicar una transformación a las señales para representarlas como una mezcla de las componentes independientes en lugar de los canales originales.
  6. Observar los patrones de actividad y determinar si los componentes hallados proporcionan información descriptiva.
  7. Una vez representados los estados de meditación y no-meditación, aplicar un algoritmo no supervisado para agrupamiento (clustering) y validar cada agrupamiento con una medida de calidad como el coeficiente silhouette o métricas de coherencia inter- e intra-clúster.

Durante el análisis reportado se realizaron los procedimientos descritos hasta el paso 6. Los resultados se muestran como una serie de gráficos e imágenes que es necesario interpretar para cada caso individual tomando en cuenta la naturaleza de la señal de EEG recuperada así como condiciones de recolección de los datos.

En el apartado siguiente se comentan los resultados obtenidos para un caso en específico, los datos para otros casos están incluidos en la carpeta de resultados correspondiente.

3. Resultados

La primera observación se realiza sobre las señales graficadas directamente. En ellas es posible observar tanto la calidad de la muestra como los distintos estados de actividad por los que atraviesa el sujeto durante el muestreo. Los resultados que aquí se comentan corresponden al sujeto “Marco”. Otros sujetos muestran características distintas. De aquí la importancia de verificar primero la señal obtenida.

Esta es la gráfica de la amplitud de la señal de cada canal en el tiempo. Se observa claramente como después de algunos artefactos a la mitad del rango de tiempo, todas las señales alcanzan un estado más o menos estacionario. De aquí proviene la decisión de analizar únicamente la región que corresponde al estado ‘meditación’, esto es la señal durante este estado estacionario más estable y que corresponde al punto de atención del experimento.

Cada señal es distinta y en el conjunto ocurre que casos como “Mauricio” o “Mireya” presentan una señal menos periódica y claramente más irregular, de donde se desprende la hipótesis de que no se trata de un solo estado de ‘meditación’ común a todos los sujetos, si no que más bien es necesario descubrir las características propias de cada uno de ellos.

Luego de la observación en el dominio del tiempo, se realiza un análisis en frecuencia, para lo que es necesario transformar las señales mediante la FFT, filtrar y aplicar el algoritmo para extracción de componentes que nos indicarán cómo es la actividad de acuerdo con la posición de los electrodos. Los resultados de cada sujetos están separados por carpetas e integrados como se describe a continuación.

             

La imagen channelLocations mostrada en primera posición indica la ubicación de los electrodos de cada canal de acuerdo con la distribución para un modelo esférico [2].

La imagen channelSpectraMaps mostrada en la segunda posición representa el espectro de frecuencia para los distintos canales (parte inferior) y la región de actividad en la frecuencia central para cada una de las bandas Delta (2Hz), Theta (6 Hz), Alpha (10 Hz), Beta (19 Hz) (parte superior) de acuerdo con el modelo y distribución de los electrodos/canales.

La imagen mostrada arriba y que lleva el nombre de componentMaps muestra las componentes calculadas, esto es, las áreas de activación que luego son usadas para representar cada una de las bandas de frecuencia. Estas componentes son identificadas por el algoritmo no supervisado referido en el paso 4 del método propuesto. Cada una está conformada como una combinación lineal e independiente de los distintos canales. Es posible eliminar componentes para los casos en que las señales de los canales son redundantes. Cada una de estas componentes tiene un espectro en frecuencia resultante que se muestra a continuación

 

Luego, cada una de estas componentes es integrada a una mezcla y usada para reconstruir las frecuencias centrales de cada banda. En cada imagen se muestran las distintas componentes señalando especialmente a aquella que hace la mayor aportación energética para una frecuencia determinada.

 

En este último grupo de imágenes se puede identificar la aportación de cada una de las componentes independientes a las distintas bandas de frecuencia. Esta es una caracterización análoga a la realizada mediante los patrones de actividad pero en lugar de poder distinguir entre estados emocionales entre distintas personas, pueden usarse para identificar el estado ‘meditación’ en el sujeto del que ha sido tomada la muestra.

Una observación a través de los resultados para distintos sujetos hace evidente que los patrones de actividad son distintos de sujeto a sujeto. La comparación se puede realizar a partir de los mapas de las componentes como se muestra a continuación:

 

Contesina

Jprico

Mafer

Marco

Mauricio

Mireya

Algunas propiedades como la similitud en el comportamiento de distintas señales en el caso de Jprico o Mafer se reflejan en que sus componentes independientes son muy parecidos entre sí (todos los de Jprico y todos los de Mafer respectivamente).

 

También a partir de esta comparativa podemos observar las distintitas zonas de activación para cada uno de los sujetos. La norma en general es que cada uno muestra un patrón de actividad distinto.

4 REFS

  • [1] Classifying Different Emotional States by Means of EEGBased Functional Connectivity Patterns. You-Yun Lee, Shulan Hsieh. PLOS One Vol 9, Issue 4. April 2014.
  • [2] Analyzing Neural Time Series Data. Mike X. Cohen. MIT PRess. 2014.

Guia de estimulos para experimento EEG.

Los primero 5 segundos de los datos son muestras de ruido previas al inicio de las indicaciones. Las indicaciones se presentaron como estimulo audio visual, utilizando un oculus rift como pantalla para simular un espacio aislado.

Preparación

0:00:00 – 0:00:10 (estimulo audio visual) Tome una posición relaja, respire profundo y trate de mantener su mente tranquila

0:00:12 – 0:00:38 (estimulo audio visual) Imagen para sincronización de la respiración

Registro de artefactos:

0:00:40 – 0:00:46 (estimulo audio visual) Por favor cierre y abra la boca como si estuviera masticando algo

0:00:47 – 0:00:52 (estimulo audio visual) Por favor, abra y cierre los ojos rápidamente.

0:00:54 – 0:00:58 (estimulo audio visual) Imagine que mueve sus extremidades.

0:01:00 – 0:01:03 (estimulo audio visual) Mueva su cabeza y cuello en círculos

Preparación

0:01:05 – 0:01:10 (estimulo audio visual) Una vez más respire profundo y trate de mantener su mente tranquila.

0:01:11 – 0:01:37 (estimulo audio visual) Imagen para sincronización de la respiración

Contexto

0:01:40 – 0:03:04 (estimulo audio visual)

  • Ahora te invitamos a hacer consciencia del lugar en el que te encuentras. El cerro del Tepeyac, en palabras de Ignacio Manuel Altamirano es el lugar en donde están todas las razas de la antigua colonia, todas las clases de la nueva República, todas las castas que viven en nuestra democracia, todos los trajes de nuestra civilización, todas las opiniones de nuestra política, todas las variedades del vicio y todas las máscaras de la virtud, en México.
  • La cuenca del valle de México es la cuna de una de las civilizaciones más complejas y poderosas de latinoamerica. En los últimos 700 años la cuenca ha sido objeto de un sin número de colonizaciones y conquistas, desde las disputas entre los pueblos indígenas hasta la llegadas de los españoles y los americanos con sus estrategias de enfoque capitalista. Cada una de estas conquista dejó su marca en nuestra psique a partir de los procesos socio culturales que llegaron a imponer.
  • A partir de estos proceso hemos sido receptores de una serie de condicionamientos, que han moldeado nuestra forma de convivir. Se trata de una programación histórica, que determina profundamente lo que deseamos.

Ejercicio

0:03:07 – 0:03:34 (estimulo audio visual)

A partir de este momento te invitamos a entrar en meditación profunda y, desde tu práctica, realizar un ejercicio que manifieste la intención a hacer conscientes nuestros deseos programados. Un ejercicio de intención a distancia que canalice energía hacia el lugar en el que te encuentras ahora.

inicie la meditación ahora

Meditación 0:03:36 – 0:08:45 (ojos cerrados y audio)

Mapa conceptual

topologiasMapa

 

I ¿Cuál es el papel de la historia en el deseo social?

  • Todas las conquistas/colonizaciones desde los mexicas hasta el salón de clases.
  • La influencia de los histórico
  • Deseo impuesto
  • Biopolíticas

II ¿Es posible pensar en una reprogramación de la memoria social?

  • Ejercicio de la conciencia
  • Metodología terapéutica
  • Hacerse cargo
  • Deseo como agente de cambio

III ¿Es posible hacer un registro fisiológico de esta intención?

  • Registros EEG no invasivo
  • ¿Arrojaría patrones entre personas en estados similares?
  • Intención a la transformación creativa/constructiva
  • Topologías
    • Modelo de mapeo de correlaciones, coherencia y fase.

IV ¿Es posible transmitir las propiedades de estos patrones a través de un acto colectivo?

  • Acto sonoro colectivo
  • Tecnoshamanismo
  • Materialización /Transducción
  • EEG Topografía 3D/
    • Self Organization
    • Graph Theory
    • Particle Spring
    • Software
  • Sonido
    • Bioacústica
    • Psicogeografía
    • Ritual sonoro
    • Psicoacústica

V ¿Es posible diseñar una interfase que expanda la conciencia de nuestra memoria?

  • Epigenetics