Актуальность проблемы. Создание автоматизированных систем медицинской и биологической диагностики в связи с современным внедрением в них цифровых средств обработки и анализа информации вошло в новую стадию своего развития. Диагностические системы базируются на вычислении разного рода показателей функционирования организма, которые получают из данных клинических обследований, в том числе получаемых из обработки электрофизиологических сигналов. Совокупность полученной информации характеризует так называемое функциональное состояние (ФС) биологического объекта.
Заключение
о заболевании осуществляются на базе оценки ФС биологического объекта (в том числе и искусственно вызванных ФС — функциональных проб) и сравнении ее с известными описаниями ФС в норме и патологии. Отсюда становится понятной актуальность поиска автоматизированных средств именно ФС биологических объектов.
Потребность в разработке современных средств и методов медицинской диагностики в условиях резкого ухудшения экологической обстановки становится все более выраженной. Быстро растущая заболеваемость и смертность указывают на необходимость поиска новых методов и систем мониторинга за изменением состояния здоровья людей и животных.
В области моделирования биологических сигналов существует множество моделей и алгоритмов с той или иной степенью точности описывающих реальный процесс, но чисто описательные характеристики процесса не позволяют его прогнозировать, что значительно снижает ценность таких моделей для экспериментальной медицины. Например, исследование деятельности мозга в реальном масштабе времени при воздействии внешних факторов, требует подхода, который позволяет прогнозировать развитие процесса по его начальным параметрам.
Такие модели необходимы как в экспериментальной так и в клинической медицине, (анестезиология, реаниматология, психопатология) где по электроэнцефалограммам определяется функциональное состояние головного мозга. Проблема восприятия времени представляет глубокий теоретический и практический интерес. В ряде исследований показано, что в отсчете времени могут принимать участие следовые процессы, уровень активации, ритмы ЭЭГ и т. д. Существует предположение о том, что ритмическая активность нервной системы играет роль шкалы отсчёта времени, но работы посвященные изучению связи между особенностями ЭЭГ и пространственно-временным восприятием недостаточно изучены. Очевидно, что особенности пространственно-временного восприятия определяются индивидуальными психофизиологическими особенностями, которые проявляются в объективных биоэлектрических показателях. ЭЭГ обладают свойствами особого рода множеств, которые принято называть фракталами. 1,3,4]. Информационно-динамические исследования ЭЭГ человека при различных умственных нагрузках получили широкие аналитические возможности благодаря использованию новых компьютерных технологий. Быстрый рост числа работ, применяющих эти технологии, разнообразие аналитических подходов, привели к получению как хорошо понятных, так и противоречивых результатов. В рамках развития новых методов системного анализа и диагностики исследования проводились с точки зрения определимости значимости ЭЭГ для количественной характеристики и объективной диагностики свойств головного мозга человека. Эти модели близки к реальному процессу, который не требует длительного времени обработки. Разработка таких алгоритмов даёт возможность их схемотехнической организации с минимальными аппаратными затратами. Такие алгоритмы могут использоваться для решения широкого круга задач цифровой обработки биологических сигналов, где требуется высокая точность и скорость. Эти алгоритмы и реализующий их процессор, могут быть применены для математического моделирования ЭЭГ, что даёт возможность выделять характерные признаки, проводить экспресс диагностику и классифицировать функциональное состояние организма человека. Разработка подобной математической программы решает важные аспекты экспериментальной медицины и является возможным средством применения в телемедицине.
Массовые формы мониторинга здоровья могут быть созданы на основе автоматических систем диспансеризации населения, которые в свою очередь могут быть созданы на базе разработок быстрых и эффективных алгоритмов поиска диагностических признаков для организма в норме и патологии. Одним из направлений такого поиска представляется анализ информации, получаемый из ритмической деятельности мозга. Интерес к такого рода системам постоянно возрастает в связи с появлением новых фактов, указывающих на связь ритмических показателей мозговой деятельности с функционированием организма в целом, с одной стороны, и специфической изменчивостью сигналов мозга, связанной с заболеванием конкретных функциональных узлов организма, с другой стороны. Установленные факты дали новый толчок к поиску системного подхода анализа информации, получаемой из сигналов электрической активности мозга, получившей название электроэнцефалограммы (ЭЭГ). Таким образом, необходимость разработки автоматизированных средств оценки ФС и экспресс контроля ФС биологических объектов по-прежнему остается актуальной темой научной и практической медицины и биологии.
Цель работы.
Целью диссертационной работы является разработка моделей, алгоритмов и программы оценки функционального состояния мозговой деятельности человека по данным электрофизиологических наблюдений и деятельности мозга, а также разработка методики для обеспечения сотрудников медицинских учреждений инструментом, позволяющим оказать помощь в определении состояния пациента и уточнения диагноза по результатам компьютерного анализа фрактальной динамики электроэнцефалограмм [1,2,3,4,5,7,8,9].
Задачи работы:
Разработка быстрого алгоритма вычисления показателей фрактальной динамики биологических сигналов.
Разработка математической модели и системы моделирования ЭЭГ для решения актуальных задач в медицине.
Разработка алгоритма экспресс-анализа по электрофизиологическим цифровым данным и решение задачи автоматической классификации ЭЭГ.
НАУЧНАЯ НОВИЗНА.
В процессе решения поставленных задач получены новые научные результаты выносимые на защиту:
Разработан быстрый алгоритм вычисления показателей (характеристик) фрактальной динамики биологических сигналов.
Выявлены параметры характеризующие функцию, которая описывает огибающую спектра мощности процесса-фрактала биологического сигнала.
Разработан алгоритм экспресс-анализа по электрофизиологическим цифровым данным.
Разработана программа для математической обработки файла с данными ЭЭГ, позволяющая получить информационные фрактальные характеристики, которые в свою очередь дают возможность провести классификацию конкретных пациентов по их состоянию.
Проведен сравнительный анализ алгоритмов дискриминантного анализа.
ПРАКТИЧЕСКАЯ ЗНАЧИМОСТЬ РАБОТЫ.
Разработан новый алгоритм обработки биомедицинской информации в реальном масштабе времени. Предложена математическая модель пространственно-временной динамики ЭЭГ позволяющая провести классификацию по состоянию адаптивности человека по группам. На основе этой классификации получена достаточная информация о состоянии испытуемого и проведено уточнение диагноза.
РЕАЛИЗАЦИЯ И ВНЕДРЕНИЕ.
Научные результаты работы с целью верификации предложенных положений и отладки разработанного алгоритма компьютерного анализа фрактальной динамики ЭЭГ были сравнены с результатами диагнозов полученных в институте экспериментальной медицины РАН под руководством доктора биологических наук Н. Б. Суворова. Исследования, проведённые им и его коллегами, показали, что каждый испытуемый может быть помещён в свою индивидуально-типологическую группу. Результаты полученные в данной диссертационной работе не расходятся с вышеупомянутыми.
АППРОБАЦИЯ РАБОТЫ.
Основные результаты диссертационной работы докладывались на: VI Санкт Петербургской Международной конференции «Региональная информатика — 98».
IV съезде Союза УИСП «Медицина, медтехника и фармация» 1998 г. Научно-практическом семинаре I Международной выставки «Интернет-Экспо», «Телемедицина — становление и развитие». Информатика и новые компьютерные технологии в образовании-1999;Санкт-Петербург, РГПУ им А. И. Герцена.
Информатика и новые компьютерные технологии в образовании-2000;Санкг-Петербург, РГПУ им А. И. Герцена.
ПУБЛИКАЦИИ.
По теме диссертации опубликовано 11 научных работ.
1. Ronald Orosco Albornos, La influencia de los campos magneticos debites y uitradebiles en los sistemas biologicos (resumen analitico). El diario, suplemento ciencia y computacion, La PazBolivia 1999, pag 6.
2. Polonnikov R. I., Orosco A. R., Modelacion de las senales debiles de importancia biologica. Bolivia. La Paz. y El diario, suplemento ciencia y computacion, La PazBolivia 1999 pag. 7,8.
3. Orosco A. R. Polonnikov R. I., Investigacion dinamica fractal de encefalogramas., El diario, suplemento ciencia y computacion, La PazBolivia 1999, pag. 8−12.
4. Polonnikov R. l, Orosco A. R, Research of the fractal dynamics of electroencephalograms-Tools for mathematical modeling. С-Петербург 1999, с 98−99.
5. Polonnikov R.I., Orosco A.R., Modeling and methods of treatment of electroencephalograms (eeg) // Информатика-исследования и инновации. Вып. З Санкт-Петербург 1999 с. 22−26.
6. Полонников. Р.И. ОроскоА. Р., Моделирование слабых биологически значимых сигналов // VI Санкт-петербургская международная конференция «Региональная информатика — 98». Тезисы докладов, часть 1, С-Петербург 1998 с. 147.
7. ОроскоА. Р., Влияние электромагнитного излучения слабых и сверхслабых полей на биологические объекты (обзор литературы) // Региональная информатика 98, СанкгПетербург 1998 с. 145 — 146.
8. ОроскоА. Р., Экспресс анализ ЭЭГ// IV съезд Союза УИСП «Медицина, медтехника и фармация"//Методы и средства анализа информации, СанктПетербург 1999 с.
3.6 Выводы по третьей главе.
Нами предложен принципиально новый алгоритм-метод анализа фрактальной динамики, свободный от анализа от отмеченных выше недостатков. В чем же состоят эти принципиальных отличия. Перечислим основные из них:
1. Нашей основной целью является получение информации, содержащейся в исследуемом процессе (развертывающейся ЭЭГ) о динамике изменения характеристик фрактала.
2. Нас интересуют, прежде всего информационные характеристики фрактала, к числу которых можно отнести:
• Фрактальную размерность (и функционального свяженные с нею величины). Форму огибающей спектра мощности ЭЭГ.
• Накопленное за определенный интервал времени количество информации по Хартли в каждом отведении.
• Форму пространственно распределения накопленного количества информации для левого и правого полушарий.
• Пространственно-временной код.
3. Все эти характеристики оцениваются в динамике от секунды к секунде всего «24 секунде и на базе полученных оценок синтезируются информативные признаки до 11 признаков, являющиеся переменными для классических процедур обучаемой автоматической классификации. Последние позволяет формулировать рекомендации для вынесения диагностического решения.
4. Все указанные процедуры требуют примерно 25 секунд и поэтому могут использоваться в медицине, телемедицине, ургентной медицине (медицине катастров), при мониторинге здоровья и т. д.
4.1.Примеры использования анализа фрактальной динамики электроэнцефалограмм в биологии и медицине.
В процессе анализа фрактальной динамики мы обнаружили характеристики дающие более точные оценки распознавания классов функционального состояния человека.
Эти характеристики приводят к необходимости выбора логической и физической структуры данных, учитывающей специфику решаемых системой задач и адекватной применяемым алгоритмам и физической организации (структуре) системы для правильного использования дискриминантного анализа.
Рассмотрим и оценим пригодность существующих подходов к представлению и организации данных:
4.1.1. Представление таблицы по шести интегральным характеристикам.
Винт" Bstdj ZBmaxj Кинт> Kstdj ZKmax).
Analysis Summary.
Classification variable: class Independent variables:
Bmean.
Bstd.
Kmean.
Kstd.
ZBmax.
Zkmax.
Number of complete cases: 35 Number of groups: 3.
Discriminant Eigenvalue Relative Canonical.
Function Percentage Correlation.
1 2,59 195 90,49 0,84 947.
2 0,272 456 9,51 0,46 273.
Functions Derived Wilks Lambda Chi-Square DF P-Value.
1 0,21 879 44,8295 12 0,0000.
2 0,785 882 7,1080 5 0,2127.
The StatAdvisor.
This procedure is designed to develop a set of discriminating functions which can help predict class based on the values of other quantitative variables. 35 cases were used to develop a model to discriminate among the 3 levels of class. 6 predictor variables were entered. The one discriminating function with P-value less than 0.05 is statistically significant at the 95% confidence level. To plot the discriminating functions, select Discriminant Functions from the list of Graphical Options. To predict new observations, select Classification Table from the list of Tabular Options.
Это процедура создана для разработки дискриминантных функций которые могут помочь определить класс, основанный на значениях других числовых переменных. 35 испытуемых было использовано для разработки дискриминантной модели на трех уровнях класса. Было введено 6 определяющих переменных. Одна дискриминантная функция со значением Р меньше 0.05 обычно является значимой в 95%.
В последующих пунктах будем увеличивать количество характеристик чтобы определить новые наблюдения.
Классификация Функциональных коэффициентов по классам.
Classification Function Coefficients for class.
1 3.
Bmean 1399.67 1461.39 1439.63.
Bstd 547.843 545.448 512.322.
Kmean 0.135 005 -0.654 901 -1.45 046.
Kstd -27.6029 -25.702 -24.4561.
ZBmax -1.19 251 -1.23 457 -1.21 083.
ZKmax 0.181 511 0.197 934 0.262 841.
CONSTANT -445.24 -480.427 -462.956.
The StatAdvisor.
This pane shows the functions used to classify observations. Ther is a function for each of the 3 levels of class. For example, the function used for the first level of class is.
— 445,24 + 1399,67*Bmean +• 547,843*Bstd + 0,135 005*Knean — 27,6029*Kst — l, 19 251*ZBinax + 0,181 511*ZKmax.
These functions are used to predict which level of class new observations belong to. For more detail, select Classification Table from the list of Tabular Options.
ЗАКЛЮЧЕНИЕ
.
Изложенные результаты являются первыми и представляют собой попытку наметить путь к созданию нового подхода в исследованиях электрической активности организма. Основной содержательной стороной такого подхода является наша гипотеза о том, что электрическая активность мозга — это некий язык, на котором организм обменивается сообщениями с внешней средой [1]. Если мы научимся понимать этот язык, то многие важные и трудные задачи диагностики и лечения будут решены. Элементами структуры этого языка являются фракталы, порождающие в динамике, которую мы исследовали, в том числе и коды.
В работе предпринята попытка применить накопленные научные знания и практический опыт работы со структурно-функциональными взаимоотношениями в головном мозге.
Нашей основной целью является получение информации, содержащейся в исследуемом процессе (развертывающейся ЗЭГ) о динамике изменения характеристик фрактала.
В рамках развития новых методов системного анализа и диагностики, исследования проводились с точки зрения определимости значимости ЗЭГ для количественной характеристики и объективной диагностики свойств головного мозга человека. Эти модели близки к реальному процессу, который не требует длительного времени обработки. Разработка таких алгоритмов даёт возможность их схемотехнической реализации с минимальными аппаратными затратами. Такие алгоритмы могут использоваться для решения широкого круга задач цифровой обработки биологических сигналов, где требуется высокая точность и скорость. Эти алгоритмы и реализующий их процессор, могут быть применены для математического моделирования ЭЭГ, что даёт возможность выделять характерные признаки, проводить экспресс диагностику и классифицировать фуйкЦибнальное состояние организма человека.
Одной из задач применения (АФДЭ) является оказания быстрой квалифицированной помощи человеку и уточнении диагноза для специалиста.
Разработка этих все более сложных аппаратно-программных функций привела к созданию новых технологий, называемых информационными и ориентированных на применение в различных прикладных областях, в том числе и в медицине. Информационными они названы потому, что главным действующим.
109 агентом в них является информация. В частности, в медицине с помощью подобных технологий осуществляется не только диагностика, но и лечение.
Создание обучающей выборки и накопление базы данных позволило определить принципы математического анализа фрактальной динамики ЭЭГ. Это попытка была реализована в программе «БАЙЕСОВСКИЙ КЛАССИФИКАТОР» для получения достоверного результата. Исследованием является накопление большой базы данных и ее последующая классификация.
Для детального описания методов, (АФДЭ) оказался полезным для экспери-ментальнь!х целей позволяющих объективно оценить степень выполнения заданной классификации.
Подобные математические программы решают важные аспекты экспериментальной медицины и являются возможным средством применения в телемедицине.