OPPORTUNITIES COMPARISON OF FRACTAL EEG-ANALIZING METHODS APPLICABLE TO DETECT CHANGES IN HUMAN BRAIN ACTIVITY AT DIFFERENT EXTERNAL LIGHTING

Full Text

Введение Исходная задача применения математических фрактальных методов для анализа ЭЭГ заключалась в обнаружении изменения активности коры головного мозга при разной освещенности сетчатки глаз человека. Обнаружение разницы между состояниями активности коры больших полушарий головного мозга при воздействии на сетчатку светом разной интенсивности необходимо для исследования механизмов воздействия естественной освещенности на мозговые меха «Инфокоммуникационные технологии» Том 12, № 2, 2014 58 Антипов О.И., Захаров А.В., Пятин В.Ф. низмы контроля циркадианных ритмов человека. Циркадианные ритмы контролируются рецепторами освещенности сетчатки глаза посредством супрахиазматических ядер гипоталамуса, а функция зрачкового рефлекса при участии мла-нопсинпродуцирующих ганглиозных клеток сетчатки состоит в регуляции светового потока. В случае положительного эффекта от применения фрактальных методов к обнаружению различия в активности головного мозга будет косвенно подтверждена гипотеза о недостаточности бытовой освещенности для полноценной регуляции циркадных ритмов, что может являться причиной целого ряда распространенных неврологических расстройств, таких как инсомния, парасомния, нарколепсия и т.д. Постановка эксперимента С целью идентификации различий в функциональной активности коры больших полушарий головного мозга при разной естественной освещенности был проведен следующий эксперимент. С помощью прибора «Нейрон-спектр/4ВКМ» (Россия) регистрировалась ЭЭГ. Электроды для снятия ЭЭГ были наложены по системе «10-20». В общем виде схема расположения электродов представлена на рис. 1. Группы электродов можно разделить на пары, условные названия которых представлены в таблице 1. Как видно из таблицы и рис. 1, пары электродов симметрично расположены относительно сагиттальной линии Fpz-Oz, Использование данной системы отведений должно позволить определить квантующие механизмы мозга, которые активируются центральными циркадианными часами при участии фоточувст-вительных ганглиозных клеток сетчатки. В эксперименте снимались 10 серий ЭЭГ-сигналов при двух различных уровнях освещенности: - при естественной освещенности внутри помещения, соответствующей норме согласно СНиП 23-05-95, но, как свидетельствуют данные ряда авторов [32-33], недостаточной по величине для активации циркадианных часов. Условно такая величина естественной освещенности была обозначена как недостаточная освещенность; - при естественной освещенности, когда взгляд испытуемого был ориентирован на голубой безоблачный небосвод в околополуденное время, прямые солнечные лучи не освещали глаза добровольца. В первом случае имитировалась рабочая обстановка, соответствующая нормативам охраны труда ГОСТ ИСО 8995-2002 для выполнения «задания со средними требованиями к условиям зрительного восприятия», т.е. в пределах 300-500 лк. Во втором случае рассматривалось моделирование ситуации нахождения человека в условиях естественной освещенности, то есть не менее 3000 лк. Измерение освещенности проводилось с помощью прибора ТКА-ПКМ (02) (Россия), зарегистрированного в Реестре средств измерений 24248-09. Для объективности эксперимента серии со световой стимуляцией чередовались с сериями без таковой. Длительность каждой серии составляла не менее 5 мин., для того чтобы можно было в каждой серии выбрать эпоху длительностью 5 мин., которая была бы свободна от очевидных двигательных артефактов. Каждая пятиминутная серия, в свою очередь, разбивалась на 10 тридцатисекундных эпох, для каждой их которой были вычислены значения фрактальных величин в соответствии с применяемыми методами. Рис. 1. Биполярная схема расположения электродов для регистрации ЭЭГ (пунктирными линиями показана биполярные пары отведений электродов, а штрихпунктирной линией представлена сагиттальная линия головы) Обзор применяемых методов Применение методов моделирования зарекомендовало себя как научно обоснованный и высокоинформативный метод суждения о различных нелинейных процессах (в том числе заболеваний) [14-17; 21], на которые оказывают влияние множество факторов [18-19]. Являясь универсальными методами для описания различных биологических процессов, математическое «Инфокоммуникационные технологии» Том 12, № 2, 2014 Антипов О.И., Захаров А.В., Пятин В.Ф. 59 моделирование нашло применение не только при описание клинических состояний, но и при изучении физиологических процессов на животных моделях [23-24; 28; 30-31], в изучении иммунных ответов и протекании физиологических иммунных реакций [20; 22; 25; 29]. Всестороннее изучение патологического процесса в конечном счете позволяет подходить к решению задач терапии, реабилитации и профилактики [20; 26]. Наибольшими возможностями в решении нелинейных задач в описании биологических процессов обладают фрактальные методы. Поскольку человеческий мозг является классической синергетической неравновесной динамической системой [27], то соответственно применяемые методы должны быть фрактальными. Так как исследуемые временные ряды являются относительно короткими [9], то и применяемые методы должны быть адаптированы к решаемой задаче соответствующим образом. В данной работе применялись следующие методы: модифицированный метод ближайших ложных соседей (БЛС), метод нормированного размаха Херста (R/S анализ) и метод аппроксимационной энтропии. Таблица 1. Условные обозначения пар электродов ЭЭГ Модификация метода БЛС заключалась в том, что в отличие от стандартного метода рассчитывалось количество ближайший ложных соседей к длине исследуемого ряда PIN при переходе из двумерного в трехмерное псевдо-фазовое пространство. Данная методика хорошо зарекомендовала себя в авторских работах при решении различных задач не только в медицине [4; 8; 12], но и в экономике [2; 6; 9], и в технике [10-11; 13]. Поскольку существует связь между показателем Херста и показателем функции огибающей спектр мощности исследуемого сигнала, то метод нормированного размаха применялся по двум разным алгоритмам. Сначала он применял ся для всего ряда целиком (полученная величина обозначена как Н), а потом вычислялась величина < Н >, полученная путем усреднения значения показателей Херста отдельных участков равномерно разделенного исходного ряда. Это связано с результатами Анищенко [1], на основании которых он сделал вывод о лучшем приближении усредненных вышеприведенным способом спектрограмм к незашумленной исходной спектрограмме применительно к нелинейным неравновесным динамическим системам. Подобная методика хорошо зарекомендовала себя в авторских работах [7]. Последним из применяемых фрактальных методов являлся метод расчета аппроксимацион-ной энтропии, который ранее авторами успешно применялся в работах [3-5], где также производился анализ ЭЭГ-сигналов. Таким образом, для каждого пятиминутного опыта, характеризуемого 10 участками длиной в стандартную эпоху в 30 сек., получается 10 значений каждой фрактальной величины. Результаты представлены на рис. 2 в виде диаграмм, известных в статистике как ящики с усами (от англ. box-and-whiskers). Ящик с усами - компактное представление распределения дискретной случайной величины. Усы ящика - предельные значения, низ и верх ящика - нижняя и верхняя квартили (первая и третья соответственно), а линия внутри ящика - медиана. Этот график позволяет визуально выделить статистически значимые результаты, входящие внутрь квартильного размаха, а положение медианы относительно квартилей позволяет оценить смещение центра распределения. Для визуального удобства опыты без световой стимуляции окрашены в серый цвет и расположены в левой половине диаграммы, а результаты опыта со световой стимуляцией расположены в правой половине диаграммы и окрашены в светло-серый свет. Обсуждение результатов Поскольку статья посвящена фрактальным методам, то физиологическая и неврологическая стороны вопроса здесь освещены не будут. Применение метода БЛС не дало ожидаемых результатов. Возможно, это связано с тем, что не было правильно подобрано значение порога [9]. Также, возможно, методика усреднения для данного метода не является оптимальной. По применению метода R/S-размаха получены два вывода: результаты усредненных показателей отличаются от показателей целого ряда Области коры больших полушарий Локализация и обозначение пар электродов Левое полушарие Правое полушарие Лобная Fpl-F3 Fp2-F4 Височная F7-T3 F8-T4 Теменная СЗ-РЗ С4-Р4 Затылочная Т5-01 Т6-02 «Инфокоммуникационные технологии» Том 12, № 2, 2014 60 Антипов О.И., Захаров А.В., Пятин В.Ф. в) г) Рис. 2. Распределения фрактальных величин для канала F7-T3 (в виде ящиков с усами), полученных с помощью следующих фрактальных методов: а) метода ложных ближайших соседей, б) показателя Херста, в) усредненного за серию значения показателя Херста, г) аппроксимационной энтропии Рис. 3. Графики распределения аппроксимационной энтропии для каналов правого полушария: a) Fp2-F4, б) F8-T4, в) С4-Р4, г) Т6-02 «Инфокоммуникационные технологии» Том 12, № 2, 2014 Антипов О.И., Захаров А.В., Пятин В.Ф. 61 Таблица 2. Усредненные значения фрактальных величин по всем серям опытов Фрактальная величина <P/N> Н <Н> <АрЕл> Достаточная освещенность нет да нет да нет да нет да Значение для канала С4-Р4 0.48 0.48 0.31 0.32 0.61 0.60 0.97 0.98 Значение для канала Т6-02 0.50 0.52 0.31 0.35 0.62 0.61 1.02 1.11 Значение для канала F8-T4 0.42 0.45 0.35 0.36 0.68 0.66 0.74 0.88 Значение для канала Fp2-F4 0.34 0.42 0.34 0.32 0.62 0.66 0.79 0.88 Значение для канала СЗ-РЗ 0.47 0.50 0.32 0.30 0.62 0.60 0.95 1.02 Значение для канала Т5-01 0.50 0.51 0.31 0.35 0.60 0.60 1.04 1.10 Значение для канала F7-T3 0.41 0.43 0.33 0.36 0.67 0.66 0.69 0.84 Значение для канала Fpl-F3 0.46 0.46 0.32 0.33 0.67 0.66 0.61 0.72 не только качественно, но и количественно. То есть значение показателя Херста для всего ряда существенно отличается от среднего значения показателей его отдельных частей. Это говорит о локальности масштабной инвариантности ЭЭГ-сигналов, что, возможно, связано с ограничениями в частотном диапазоне работы головного мозга. Точность фиксирования сверхслабых сигналов, которыми являются ЭЭГ, ограничиваемая минимальным уровнем дискретизации, не позволяет различить фрактальные свойства на малых временных масштабах. Главный вывод заключается в том, что метод аппроксимационной энтропии позволил получить хорошо дифференцируемые результаты с хорошей статистической достоверностью, как показано на рис. 2 г. Также следует отметить, что подобный результат достигается только для отведений электродов, расположенных в проекции височных и затылочных областей коры больших полушарий головного мозга. Для иллюстрации данного вывода на рис. 3 представлены результаты применения метода аппроксимационной энтропии для различных пар отведений другого полушария. Общая картина применения методов фрактального анализа представлена в таблице 2 в виде упорядоченных усредненных значений для всех серий экспериментов отдельно для каждого условия эксперимента: недостаточная освещенность и естественная освещенность. Заключение Таким образом, в соответствии с задачей эксперимента - применение математических фрактальных методов для анализа ЭЭГ при интенсивной освещенности сетчатки глаз человека - было обнаружено изменение активности коры больших полушарий головного мозга в височных и затылочных областях симметрично относительно сагиттальной линии. Важной гипотезой, полученной в ходе проведенного анализа, является предположение, что на результаты применения фрактальных методик сильное влияние оказывают двигательные артефакты, лежащие в дельта-диапазонах ЭЭГ. Эта гипотеза, возможно, позволит повысить достоверность полученных результатов путем полосовой фильтрации в указанных диапазонах.

About the authors

O. I Antipov

Email: oleg1307@mail.ru

A. V Zakharov

Email: zakharov1977@mail.ru

V. F Pyatin

Email: pyatin_vf@list.ru

References

Supplementary files