УДК 616.858-073.96-072.7-8
ПРИМЕНЕНИЕ ВИДЕОАНАЛИЗА ДВИЖЕНИЙ ДЛЯ ИССЛЕДОВАНИЯ ФУНКЦИИ РАВНОВЕСИЯ У ПАЦИЕНТОВ С БОЛЕЗНЬЮ ПАРКИНСОНА
С.А. ЛИХАЧЕВ, ВВ. ВАШИЛИН, И.С. ГУРСКИЙ
РНПЦ неврологии и нейрохирургии Ф. Скорины, 24 Минск, 220114, Беларусь
Поступила в редакцию 22 ноября 2016
Описан способ контроля постуральной функции у пациентов с болезнью Паркинсона, который не требует дорогостоящего оборудования и может быть использован в комплексе до- и послеоперационного обследования пациентов с болезнью Паркинсона.
Ключевые слова: болезнь Паркинсона, видеоанализ движений, глубокая стимуляция мозга, равновесие.
Введение
В развернутых стадиях болезни Паркинсона (БП), наряду с олигобрадикинезией и мышечной ригидностью, значимой и трудно поддающейся коррекции проблемой становится нарушение равновесия. Поддержание вертикальной позы и равновесия требует быстрой и точной коррекции тонуса определенных групп мышц в соответствии с информацией, поступающей от проприоцептивного, вестибулярного и зрительного анализаторов. С внедрением хирургических методов лечения БП, в частности - DBS, значимость проблемы возрастает. Усиление постуральных нарушений является одним из нежелательных эффектов DBS. Изучение эффективности и безопасности оперативного лечения требует объективного контроля постуральной функции. Известны способы оценки постуральной функции, основанные на использовании стабилоплатформы.
Целью настоящей работы явилась разработка способа исследования выраженности нарушения равновесия при БП, основанного на применения метода видеоанализа движений.
Методика эксперимента
Для изучения состояния функции поддержания вертикальной позы и равновесия пациент располагается перед цифровой камерой, которая подключена к персональному компьютеру. На передней поверхности туловища пациента в области передней брюшной стенки по средней линии и в области плечевых суставов закрепляются цветные маркеры (использовались окрашенные в красный цвет шарики, которые хорошо контрастировали с одеждой и кожей пациента, обстановкой помещения), которые облегчают решение задачи автоматического измерения движений. При использовании системы видеоанализа движений с одной камерой затруднительно провести триангуляцию для измерения положения маркера в пространстве. Возможно применение зеркал или призм, обеспечивающих обзор с двух точек, что, однако, значительно усложняет методику исследования. В данном исследовании для оценки расстояния от туловища пациента до камеры используется эффект перспективы (при приближении туловища пациента к камере его видимые размеры, и, в частности, расстояние между маркерами над плечевыми суставами увеличиваются, при удалении - наоборот, уменьшаются). Для оценки перпендикулярных к оптической оси камеры смещений используется видимое положение маркера на передней брюшной стенке.
Разработанная авторами программа обработки видеоизображений определяет
координаты центров маркеров на поступающих с камеры видеокадрах, и в реальном времени выводит на экран компьютера наглядную информацию для пациента о положении его тела. На экране отображаются две фигуры разного цвета, одна из которых отображает положение тела пациента (отметка пациента - ОП), а другая (мишень) перемещается по экрану в соответствии с заданной программой. Пациента инструктируют совмещать ОП с движущейся меткой-мишенью, совершая движения туловищем без наклонов в пояснице, не отрывая ноги от пола. ОП смещается в соответствии с перемещениями туловища: вправо и влево - при перемещении маркера на брюшной стенке вправо и влево соответственно; вверх - при смещении туловища вперед; вниз - при смещении туловища назад. Координаты мишени изменяются по синусоидальному закону, при этом частота изменения координаты y в 2 раза выше, чем координаты x:
x = хс + A■ sin(f ■ 2к-1), (1)
yt = yc + A ■ sin(2 ■ f ■ 2n ■ t), (2)
где xt и yt - координаты мишени, xc и yc - координаты центра окна на экране, A - настраиваемая амплитуда колебаний мишени, f - настраиваемая частота (в данном исследовании f = 0,05 Гц).
Для определения координат маркеров на кадре используется техника сравнения с шаблоном. Сравнение цвета пикселя цветом маркера проводится в цветовом пространстве HSV. Величину D, которая отражает степень соответствия цвета пикселя цвету маркера,
255 ■ S
рассчитывают по формуле D =-=-- , где S и H -
1 + 4 ■ (sin(H0) - sin( H) ) + 4 ■ ( cos( H0) - cos(H) )
значения насыщенности и оттенка цвета пикселя, H0 - значение оттенка цвета маркера. Чем больше величина D, тем ближе цвет пикселя по оттенку и насыщенности к цвету маркера. Значение (V) при сравнении не учитывается.
С целью уменьшения затрат времени на расчет величины D, ее значения для всех возможных значений R, G и B (256-256-256 возможных сочетаний) при заданном цвете маркера сохраняются в таблице. Степень соответствия изображению маркера с координатами центра
(x,y) и радиусом r оценивается функцией W(x, y, r), рассчитываемой по формуле:
( \
\ s(х, y, r) s(х, y, 2r) - s(x, y, r) W(x, y, r) = r ■ —-f—1-^-1-TL
Ц2r +1)2 (4r +1)2-(2r +1)2 , Здесь s(x, y, p) - сумма значений D в квадратной области изображения с координатами
х+р y+p
центра (xy) и стороной, равной 2p+1, определяемая как s(х, y, р) = ^ ^ Д. ., где i, j -
1=х-р j=y-p
целые числа, D- значение величины D для пикселя с координатами i, j; r - радиус маркера; x и y - координаты центра маркера. Для быстрого вычисления значения s(x,y,p) используется
i \ 1 j формула: s(x, y, р) = ^+р,у+р- ^х+р,у-р-1 ( 2х-р-1,у+р- 2х-р-1,у-р-1) , где Zi, j = ZZ Dk ,l , причем
k=0 l=0
таблица значений z,- быстро рассчитывается за время порядка O(N), где N - количество
пикселей в области интереса, с использованием формул:
Координаты (xm%ym) центра маркера рассчитывают по формулам: хт = argmax W(х, y, r) .
х
Ут = argmax W (х, У, r) , х ^ l^mm , хтах У Уш^ Уmax], r ^r^ rmax] , где Xmm и xmax -
c = c. , + D , c = D„
',] 1-1,j i,j ' 0,} 0,j zi, j = Zi, j-1 + ci, j , Zi,0 = ci,0
min? max^^ L^ min? У max J? L min? max J
У
минимальное и максимальное значение координаты х, задающие область поиска маркера на изображении, Ушт и утах минимальное и максимальное значение координаты у, задающие область поиска маркера на изображении, гшп и Гтах - минимальное и максимальное возможные значения радиуса маркера г. Координата х ОП в зависимости от видимого положения маркера на передней брюшной стенке рассчитывается по формуле: х = хс-ЕХ-(Х-Хо), где Ех -
настраиваемый масштаб, X - текущее положение центра маркера на кадре (после сглаживания), Xo - положение центра маркера на кадре в начале исследования. Координата y отметки на экране в зависимости от видимого расстояния между центрами маркеров над плечевыми суставами рассчитывается по формуле: y = yc-Fy-(d-do, где Fy - настраиваемый масштаб, d -текущее расстояние между центрами маркеров на кадре (после сглаживания), do - расстояние между центрами маркеров на кадре в начале исследования. Исследование выполняется у пациента в on-периоде и продолжается 10 минут, за это время 3 раза изменяется направление движения мишени. Критерием качества выполнения теста служат коэффициенты ранговой корреляции Spearman между положением мишени и отметки пациента на экране для движений в сагиттальной плоскости и во фронтальной плоскости.
Результаты и их обсуждение
Разработанная методика контроля постуральной функции была применена у пациента с БП, которому было проведено нейрохирургическое лечение - имплантации DBS. Оценка по 3 части шкалы UPDRS до операции составила 36 в off-периоде и 19 в оп-периоде, через 1 мес после операции - 7 в off-периоде, 0 в оп-периоде. Оценка по шкале Schwab-England on/off до операции составляла 90% / 80%, через 1 мес после операции - 100% / 90%. В дооперационном периоде при исследовании постуральной функции коэффициент корреляции для движений во фронтальной плоскости у данного пациента составил 0,964, для движений в сагиттальной плоскости коэффициент корреляции составил 0,906. При исследовании постуральной функции в послеоперационном периоде (через месяц после установки DBS), коэффициент корреляции для движений во фронтальной плоскости у данного пациента составил 0,905, для движений в сагиттальной плоскости коэффициент корреляции составил 0,755. Таким образом, данное исследование, несмотря на отсутствие субъективного ухудшения, демонстрирует возможное усиления нарушения функции поддержания равновесия у пациентов с БП в послеоперационном периоде DBS, поддержания равновесия пациентом. Планируется проспективное контролируемое исследование нарушения постуральной функции в до и послеоперационном периоде у пациентов с БП, которым будет имплантироваться DBS в РНПЦ неврологии и нейрохирургии.
Заключение
Разработан способ исследования выраженности нарушения равновесия при болезни Паркинсона на основе метода видеоанализа движений. Программа обработки видеоизображений определяет координаты центров маркеров на поступающих с камеры видеокадрах и в реальном времени выводит на экран компьютера наглядную информацию для пациента о положении его тела. Описанный способ исследования постуральной функции у пациентов с болезнью Паркинсона не требует дорогостоящего оборудования, и может быть использован в комплексе до- и послеоперационного обследования пациентов с болезнью Паркинсона.
UTILIZATION OF VIDEOMOTION ANALYSIS FOR ASSESSMENT OF BALANCE IN THE PATIENTS WITH PARKINSON'S DISEASE
S.A. LIKHACHEV, V.V. VASHCHYLIN, V.V. GOURSKY
Abstract
The objective of the research was to develop a method for the assessment of balance impairment in Parkinson's disease, based on videomotion analysis. The developed method for the assessment of balance function in the patients with Parkinson's disease does not require expensive equipment, and may be utilized while pre- and postoperative examination of patients with Parkinson's disease.
Keywords: Parkinson's disease, videomotion analysis, DBS, balance.