Научная статья на тему 'Анализ данных и уменьшение размерности сигнала с помощью метода главных компонент для задачи удаления артефактов (шума)'

Анализ данных и уменьшение размерности сигнала с помощью метода главных компонент для задачи удаления артефактов (шума) Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
347
61
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЭЭГ / ИНТЕРФЕЙС МОЗГ-КОМПЬЮТЕР / ИМК / АНАЛИЗ ГЛАВНЫХ КОМПОНЕНТ / АГК / НЕЙРОУПРАВЛЕНИЕ / НЕЙРОКОНТРОЛЛЕРЫ / УДАЛЕНИЕ АРТЕФАКТОВ / СИГНАЛЫ / ОБРАБОТКА / EEG / BRAIN-COMPUTER INTERFACE / IMC / ANALYSIS OF MAIN COMPONENTS / AGC / NEURO MANAGEMENT / NEUROCONTROLLER / REMOVAL OF ARTIFACTS / SIGNALS / PROCESSING. _______________

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Калашников Евгений Александрович, Кондыбаева Алмагуль Бауржановна

Актуальность выбранной темы обусловлена необходимостью удаления артефактов получаемых с электроэнцефалографических данных для дальнейшей класстеризации по пространствам состояний в области задач построения интерфейса мозг компьютер (ИМК), класстеризации вызванных эмоциональных состояний.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Калашников Евгений Александрович, Кондыбаева Алмагуль Бауржановна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

DATA ANALYZYS AND THE PROCESS OF REDUCING THE DIMENSION OF THE DATA INPUT SIGNAL WITH THE PRINCIPAL COMPONENT METHOD FOR THE PROBLEM OF REMOVING ARTIFACTS (NOISE)

The relevance of the topic chosen due to the need of the detecting the artefact noise signals in the processing of the raw electroencephalogram (EEG) signals and their further use in the brain computer interface (BCI) and emotion classterization systems.

Текст научной работы на тему «Анализ данных и уменьшение размерности сигнала с помощью метода главных компонент для задачи удаления артефактов (шума)»

ИНФОРМАТИКА, ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА И УПРАВЛЕНИЕ

АНАЛИЗ ДАННЫХ И УМЕНЬШЕНИЕ РАЗМЕРНОСТИ СИГНАЛА С ПОМОЩЬЮ МЕТОДА ГЛАВНЫХ КОМПОНЕНТ ДЛЯ ЗАДАЧИ УДАЛЕНИЯ АРТЕФАКТОВ (ШУМА)

Калашников Евгений Александрович

канд. техн. наук, доцент, профессор, кафедра автоматизированных систем управления Национальный исследовательский технологический университет «МИСиС»,

119991, РФ, Москва, Ленинский пр-т, 4

Кондыбаева Алмагуль Бауржановна

магистр, кафедра автоматизированных систем управления Национальный исследовательский технологический университет «МИСиС»,

119991, РФ, Москва, Ленинский пр-т, 4 E-mail: ek. misis@gmail.com

DATA ANALYZYS AND THE PROCESS OF REDUCING THE DIMENSION OF THE DATA INPUT SIGNAL WITH THE PRINCIPAL COMPONENT METHOD FOR THE PROBLEM OF REMOVING ARTIFACTS (NOISE)

Evgeniy Kalashnikov

candidate of Technical Sciences, Associate Professor, Professor, Department of Automated Control Systems

National Research Technological University "MISiS", 119991, Russia, Moscow, Leninsky Prospekt, 4

Almagul Kondybayeva

master's Degree, Department of Automated Control Systems National Research Technological University "MISiS",

119991, Russia, Moscow, Leninsky Prospekt, 4

АННОТАЦИЯ

Актуальность выбранной темы обусловлена необходимостью удаления артефактов получаемых с электроэнцефалографических данных для дальнейшей класстеризации по пространствам состояний в области задач построения интерфейса мозг - компьютер (ИМК), класстеризации вызванных эмоциональных состояний.

ABSTRACT

The relevance of the topic chosen due to the need of the detecting the artefact noise signals in the processing of the raw electroencephalogram (EEG) signals and their further use in the brain computer interface (BCI) and emotion classter-ization systems.

Ключевые слова: ЭЭГ, интерфейс мозг-компьютер, ИМК, анализ главных компонент, АГК, нейроуправле-ние, нейроконтроллеры, удаление артефактов, сигналы, обработка.

Keywords: EEG, brain-computer interface, IMC, analysis of main components, AGC, neuro management, neuro-controller, removal of artifacts, signals, processing.

В данной статье рассматриваются вопросы обработки сигналов электроэнцефалографии (ЭЭГ) для интерфейса мозг-компьютер (ИМК) для удаления артефактов, в основе которых лежит анализ сенсомо-торных ритмов при помощи анализа главных компонент (АГК) сигналов получаемых с неинвазивных монополярных электродов. На момент написания статьи опубликовано множество значимых работ и достигнуты теоретические, практические и экспери-

ментальные успехи в области нейроуправления (технологии ИМК), так, например в работах [1, с. 167] подробно рассматриваются вопросы методологии снятия сигналов, их оцифровки и подавления шумов. В работах [3, с. 155] подробно описывается и перечисляются методы обработки сигналов и удаления артефактов.

Задача детектирования сигналов осложняется наличием так называемых артефактов. Артефакты по

Бибилографическое описание: Кондыбаева А.Б., Калашников Е.А. Анализ данных и уменьшение размерности сигнала с помощью метода главных компонент для задачи удаления артефактов (шума) // Universum: Технические науки: электрон. научн. журн. 2017. № 12(45). URL: http://7universum. com/ru/tech/archive/item/5417

их происхождению можно разделить на две группы: физические и физиологические [1, с. 78].

При анализе ЭЭГ в данной статье для удаления артефактов используются следующие методы:

а) Применение высочастотных и низкочастотных фильтров на уровне аналого-цифрового преобразователя (АЦП) [1, с. 34],

б) Перерасчет значений сигнала тока смещения относительного нового референтного потенциала (моноподного электрода),

в) Методы, основанные на линейном разложении сигнала ЭЭГ на составляющие компоненты, выделении артефакта и восстановлении ЭЭГ без компоненты артефакта. К методам данного типа относятся анализ главных компонентов (АГК) и нижеописанный анализ независимых компонентов (АНК) [1, с. 34].

ФИЛЬТРАЦИЯ СИГНАЛА

Пусть исходный сигнал X в каждом канале соответствует вектору размерности пространства соответственно. Тогда в векторной форме сигнал может быть представлен в виде:

X = (х(0),..., х(п),..., х^ - 1))

'Х1(0) ... Х1(М-1)

... хт(п) ... кхм(0) ... хм(Ы-1)/

где х(п) = (х1(п),...,хт(п),...,хм(п) )Т - значения сигнала в момент времени с номером п е {0,1, 1};

хт (п) значение т-й компоненты сигнала в указанный момент времени, т е {1,2,..., М}.

В общем случае дискретный фильтр описывается следующим образом[5, с. 45] и представлен на рисунке 1:

У(п) = %=о Ъ1х(п -I)- аку(п

к),пе {0,1,., N - 1},

где у(п) - отсчет выходного сигнала в текущий момент времени п;

х(п - С) - отсчеты входного сигнала в предыдущие моменты времени п - I, 1е {0,1, 1},

у(п- к) - отсчеты выходного сигнала в предыдущие моменты п - к, 1е {0,1, 1}, ЪI - коэффициенты входного сигнала, ак - коэффициенты обратной связи Р - порядок входного сигнала Q - порядок обратной связи

ш =

1 + 12=1

у(г),х(г) - Ъ - преобразования для сигналов у(г), х(г) соответсвенно; г = , / е [0,1] - относительная частота сигнала (отношение частоты сигнала к частоте дискретизации). С групповой задержкой:

а

Ф(/)=Агд(Н(г)).

Рисунок 1. Фильтрация ЭЭГ сигнала: а) Сигнал с высокочастотной помехой; б) Сигнал после проведенной

обработки (фильтрации) 0,5 - 49 Гц

МОНТАЖ ОТВЕДЕНИИ ЭЭГ

Представляет собой фиксированный набор отведений, благодаря которому без помощи селектора записываются специально подобранные или стандартные программы отведений. В данной работе используется референциальный — монтаж, состоящий из референциальных отведений;

Перерасчет и приведение к единой шкале значений потенциалов производятся по следующей формуле:

1 М

1=1

где Х1 - строки матрицы Х(формула (1)).

В идеальном случае, если бы электроды были равномерно расположены на сфере, общим средним числом для А^был бы нуль и переход, выполненный в соответсвии с выражением , что позволило бы аос-станавливать достаточно точные значения потенциалов ЭЭГ.

АНАЛИЗ ГЛАВНЫХ КОМПОНЕНТ

Вследствии движения человеческих глаз (глазных яблок) возникают артефакты, вызываемые изменением ориентации электрической оси глаза. Максимальная амплитуда артефактов наблюдается в лобных отведениях и уменьшаются в отведениях затылочной части. Для удаления из многоканального ЭЭГ сигнала артефактов, обусловленных непроизвольными движениями глаз испытуемого, применяется, так называемый метод главных компонент.

Данный метод позволяет находить ортогональные преобразования в новую систему координат, для которых выполняются следующие условия:

• Максимальна выборочная дисперсия первой координаты

• При условии ортогональности ^-п) первым координатам выборочная дисперсия данных вдоль k координаты максимальная.

Направления базисных векторов будут выбраны так, что коэффициент

ковариации между проекциями исходного набора данных на различные координатные оси будет равен нулю.

У = АХ,

где X - исходный многоканальный ЭЭГ сигнал размерности MxN,

A - матрица преобразования (нагрузок) размерности KxM,

Y - некоррелированный многомерный временной ряд (KxM), представляющий собой проекцию исходного сигнала на К главных компонентов.

Алгоритм метода анализа главных компонент:

Для исходного множества сигналов ЭЭГ X производится расчет ковариационной матрицы С =[с I ¿}мхм :

су = соу{х1,х]) ЕХд(Х3 -

ЕХ}) Л,] = [1,2, ...,М},

где Х^,Х] — строки матрицы X, E - оператор математического ожидания.

Ведется поиск собственных чисел Л и собственных векторов VI ковариационной матрицы С

Матрица преобразования А формируется из первых К собственных векторов VI и располагаются в порядке убывания соответсвующих собственных значений Л ковариационной матрицы С:

А=Ш=1 К> Л2> • > Лк > 0,0 < К < М

В задачах удаления артефактов из многоканального ЭЭГ сигнала в матрицу преобразования А включаются все собственные вектора/

При наличии записи ЭОГ процедура восстановления ЭЭГ сигнала без компоненты артефакта сводится к следующему:

Для каждой строки матрицы Y (формула (3)) рассчитывается коэффициент корреляции с каждым из каналов записи электроокулограммы.

Те строки, для которых коэффициент корреляции хотя бы с одним из каналов ЭОГ превышает по модулю некоторое пороговое значение, удаляются из матрицы Y. Соответствующие им строки матрицы А также удаляются.

В результате получают-ся матрицы У, Д размерностей (К^)ХМ соответственно, где d - общее число удаленных строк матрицы Y.

Восстановленный ЭЭГ сигнал без компоненты артефакта рассчитывается как:

{Х) = Ат?

РЕЗУЛЬТАТЫ ПРОВЕДЕННЫХ ОПЫТОВ

В результате исследования были замерены 450 единиц значений сигналов в дискретный момент времени в течении тридцати минут. Таким образом мы получили 15 измерений в минуту по 6 пятиминутных блоков. Блок состоит из семидесяти пяти элементов.

Всего исходных отведений: 40 единиц.

Рисунок 2. Результат исследования

Рисунок 3. Результат исследования

ВЫВОДЫ

В результате проведенной работы были произведены наглядные демонстрации по обработке сырого ЭЭГ сигнала для удаления артефактов посредством метода главных компонент. Каждый исходный сигнал был обработан и восстановлен в соответсвии с алгоритмом описанным в методолоогической части работы. На выходе после обработки получили очищенные от артефактов данные. Точность метода при

определении матрицы собственных чисел составила менее 1%, что является сравнительно невысокой степенью погрешности. При работе с методом было выяснено, что ранг полученной матрицы значительно уменьшается по сравнению с исходной и зависит от поставленной задачи. Таким образом в работе с выявлением характерных значений ЭЭГ был выбран ранг равный 75 при исходном количестве столбцов равному 450.

Список литературы:

1. Каплан А.Я., Кочетова А.Г., Шишкин С.Л., Басюл И.А., Ганин И.П., Васильев А.Н., Либуркина С.П. Экспериментально-теоретические основания и практические реализации технологии интерфейс мозг-компьютер. Бюллетень Сибирской медицины. 2013. Т.12, №2, С. 21-29.

2. Зенков Л.Р. Клиническая электроэнцефалография (с элементами эпилептологии). Руководство для врачей. 3-е изд. Москва: МЕДпрессинформ, 2004. 368 с.

3. Рабинер Л., Гоулд Б. Теория и применение цифровой обработки сигналов. Москва: Издательство «Мир», 1978.

4. Fukunaga K. Introduction to statistical pattern recognition. 2nd ed. San Diego: АР, 1990. 614 pp.

5. Kalcher J., Pfurtscheller G., "Discrimination between phase-locked and non-phase-locked event-related EEG activity," // Electroencephalography and Clinical Neurophysiology, Vol. 94(5), 1995. pp. 381-384.

6. Zhang D., Wang Y., and Gao X., "An Algorithm for Idle-State Detection in Motor-Imagery-Based Brain-Computer Interface," // Computational Intelligence and Neuroscience, Vol. 2007, 2007. 9 pp.

7. Brunner C., Billinger M., and Vidaurre C., "A comparison of univariate, vector, bilinear autoregressive, and band power features for brain-computer interfaces," // Med Biol Eng Comput, Vol. 49, 2011. pp. 1337-1346.

8. Ludwig K.A., Miriani R.M., and Langhals N.B., "Using a Common Average Reference to Improve Cortical Neuron Recordings From Microelectrode Arrays," // Journal of Neuro-physiology, Vol. 101, 2009. pp. 1679-1689.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

9. Кестер У. Проектирование систем цифровой и смешанной обработки сигналов. Москва: Техносфера, 2010. 328 с.

i Надоели баннеры? Вы всегда можете отключить рекламу.