Метод извлечения и статистической оценки параметров походки из данных двумерного скелетного видеоанализа

Бесплатный доступ

Предложен метод извлечения биомеханических параметров походки из данных двумерного скелетного видеоанализа и их статистической оценки относительно стратифицированных популяционных норм, не требующий калибровки камеры и измерения антропометрических характеристик обследуемого. На выборке из 398 участников открытого набора Health&Gait установлено, что нормализация по полной высоте тела обеспечивает наименьший коэффициент вариации масштаба CV = 3,16 % и оптимальна для 88,4 % участников. Сформирован набор из 16 скорируемых параметров в шести клинических доменах, построены перцентильные нормативные интервалы для шести демографических групп по стандарту CLSI EP28-A3c и ГОСТ Р 53022.3-2008, разработана скоростно-условная оценка с трёхуровневой шкалой тяжести отклонений. Валидация параметров относительно инструментальной системы OptoGait показала индивидуальную корреляцию r = 0,365–0,388, соответствующую опубликованному авторами набора данных диапазону и отвечающую физическому пределу 2D-видеоанализа без калибровки; средняя абсолютная ошибка длины шага на 32 % ниже авторского результата. Нормативные медианы согласованы с медианами инструментальной системы OptoGait (средний знаковый сдвиг менее 1 %), а нормативные центры переносятся между двумя независимыми наборами данных (отношение медиан 0,94–0,98). Диагностическая валидация на клинической когорте KOA-PD-NM (187 записей пациентов с гонартрозом и болезнью Паркинсона и здоровых участников, зарегистрированных иной камерой) показала специфичность 98,9 % на внутренней leave-one-out выборке, 84,7 % — на внешней, чувствительность 40,6 % для значимых и 99,2 % для любых уровней отклонения. Предложенный метод применим для популяционного скрининга и мониторинга нарушений походки вне специализированной лаборатории.

анализ походки \ безмаркерный видеоанализ \ оценка позы \ пространственная нормализация \ безразмерные параметры \ статистические нормы \ перцентильные интервалы \ скоростно-условная оценка \ биомеханические параметры \ стратификация

Короткий адрес: https://sciup.org/148334346

IDS: 148334346   |   УДК: 004.932.2:519.25   |   DOI: 10.37313/1990-5378-2026-28-4-212-223

Method for Extraction And Statistical Assessment of Gait Parameters from 2D Skeletal Video Analysis Data

A method forinvariant extraction and statistical assessment of human gait parameters from twodimensional skeletal video analysis is proposed. The method requires no camera calibration or anthropometric measurements of the examined subject. On a sample of 398 participants from the open Health&Gait dataset, normalization by full body height is shown to provide the lowest coefficient of variation of the scale, CV = 3.16%, and is optimal for 88.4% of the participants. A set of 16 scored parameters within six clinical domains is constructed, percentile normative intervals are built for six demographic groups following the CLSI EP28-A3c and GOST R 53022.3-2008 standards, and a speed-conditional assessment with a three-level severity scale is developed. Validation against the instrumented OptoGait system yielded individual correlations r = 0.365–0.388, matching the range reported by the authors of the dataset and corresponding to the physical ceiling of 2D video analysis without camera calibration; the mean absolute error of step lengthis 32% lower than the original authors’ result. The normative medians agree with the OptoGait system medians (mean signed bias below 1%), and the normative centres transfer between two independent datasets (median ratio 0.94–0.98). Diagnostic validation on the KOA-PD-NM clinical cohort (187 records from patients with knee osteoarthritis and Parkinson’s disease and healthy controls acquired with a different camera) yielded a specificity of 98.9% on the internal leave-one-out cohort and 84.7% on the external cohort, a sensitivity of 40.6% for significant deviations and 99.2% for deviations of any level. The proposed method is applicable to population-level screening and monitoring of gait disorders outside specialized laboratory settings.

Текст научной статьи Метод извлечения и статистической оценки параметров походки из данных двумерного скелетного видеоанализа

Ходьба служит основным способом передвижения человека и одновременно одним из наиболее информативных индикаторов состояния здоровья. Нарушения походки сопровождают широкий спектр неврологических, ортопедических и нейродегенеративных заболеваний, включая болезнь Паркинсона, болезнь Альцгеймера, инсульт и остеоартроз крупных суставов [1, 2, 3], а вариабельность пространственно-временных характеристик шага, в частности коэффициент вариации длительности шага, признана цифровым биомаркером ранних стадий нейродегенеративных заболеваний [4, 5, 6]. Количественная оценка параметров походки обеспечивает раннюю диагностику нарушений опорно-двигательного аппарата и мониторинг реабилитационных мероприятий, однако массовое применение подобных методов сдерживается ограничениями существующего инструментария.

Золотым стандартом клинического анализа походки остаются маркерные оптоэлектронные системы (Vicon, Qualisys) и инструментальные дорожки (GAITRite, OptoGait), обеспечивающие точность до 1–2 мм для пространственных параметров и до 10 мс для временных [7, 8]. Стоимость оборудования, необходимость калибровки и стационарное размещение ограничивают использование подобных систем рамками клинических лабораторий, исключая возможность массового популяционного скрининга. Альтернативу представляет видеоанализ на основе нейросетевых моделей оценки позы (pose estimation); архитектуры AlphaPose, OpenPose и MMPose позволяют извлекать двумерные координаты суставов из записи одной камеры с погрешностью, сопоставимой с маркерными системами для крупных суставов в сагиттальной плоскости [9, 10, 11, 12]. В систематическом обзоре [13] показано, что пространственно-временные характеристики безмаркерного анализа демонстрируют хорошую и отличную надёжность (ICC > 0,80), а согласованность кинематических

переменных варьируется от умеренной до отличной для тазобедренного и коленного суставов. Безмаркерный видеоанализ уже применяется в клинической оценке походки, в том числе у постинсультных пациентов [14], и допускает трёхмерную реконструкцию движения по нескольким камерам [15]. В отечественной литературе развиваются подходы к анализу цикла шага на основе скелетных моделей, в том числе для задач идентификации личности [16, 17]. Вместе с тем координаты суставов в выходе нейросети представлены в абсолютных пикселях и зависят от расстояния до камеры, разрешения видеозаписи и антропометрических характеристик обследуемого, что исключает прямое сопоставление параметров между участниками без предварительной нормализации.

В работе Hof [18] предложена нормализация параметров походки к длине ноги с учётом ускорения свободного падения, получившая широкое распространение при обработке данных лабораторных систем. Muro-de-la-Herran и соавт. [19] предложили метод нормализации по диагонали «плечобедро» (Anatomical Shoulder-Hip, ASH) для оценки степени выраженности двигательных нарушений и показали его эффективность при фронтальной проекции. При латеральной проекции плечевой и тазобедренный суставы проецируются в близко расположенные точки, что приводит к вырождению ASH и делает его непригодным для бокового ракурса, преобладающего в клинической практике. Задача систематического сравнения методов нормализации координат скелетной модели для латеральной проекции до настоящего времени не решена.

Помимо инвариантного извлечения параметров, количественная интерпретация последних требует наличия референтных интервалов, стратифицированных по полу и возрасту. Ряд нормативных баз пространственно-временных параметров [7, 20, 21, 22, 23] построен на данных лабораторных систем и оперирует размерными величинами; нормативные интервалы для безразмерных параметров, извлечённых из двумерного видеоанализа без калибровки камеры, в литературе отсутствуют. Скоростно-условная оценка, устраняющая конфаундер темпа ходьбы, для безразмерных параметров двумерного видеоанализа ранее не реализована. Отсутствие согласованного набора «инвариантные параметры — стратифицированные нормы — скоростная коррекция» и сдерживает перенос видеоаналитических методов из лаборатории в повседневную клиническую практику.

Настоящая работа решает задачу извлечения параметров походки — их количественного измерения по видеозаписи, а не задачу распознавания образа походки или автоматической постановки диагноза. Под извлечением параметра понимается вычисление скалярной характеристики цикла шага (например, нормированной длины шага или амплитуды сустава) из последовательности координат скелетной модели; под пространственной нормализацией — приведение пиксельных координат к безразмерному виду делением на устойчивый базовый отрезок скелета; под скоростно-условной оценкой — сравнение параметра с нормой, зависящей от темпа ходьбы. В общем маршруте инструментальной диагностики нарушений походки предлагаемый метод занимает нишу количественного скринингового слоя: он формирует объективный профиль отклонений параметров от популяционной нормы, тогда как клиническая интерпретация этого профиля и постановка диагноза остаются в компетенции врача.

Цель работы. Разработка метода извлечения биомеханических параметров походки из данных двумерного скелетного видеоанализа и их статистической оценки относительно стратифицированных популяционных норм, инвариантного к условиям съёмки и не требующего калибровки камеры. Метод включает эмпирическое сравнение четырёх методов пространственной нормализации, формирование набора скорируемых параметров с доменной группировкой на основе стандартной биомеханической классификации [24, 25], построение перцентильных нормативных интервалов по стандарту CLSI EP28-A3c [26] и ГОСТ Р 53022.3-2008 [27] и разработку скоростно-условной оценки с трёхуровневой шкалой тяжести отклонений.

ОБЩАЯ СХЕМА МЕТОДА

Суть предлагаемого метода состоит в последовательном преобразовании исходных данных. Пиксельные координаты скелетной модели нормируются по полной высоте тела, что делает параметры безразмерными и устраняет зависимость от масштаба съёмки; из нормированных координат вычисляются 16 параметров, сгруппированных в шесть клинических доменов; для каждого параметра строятся перцентильные нормативные интервалы, стратифицированные по полу и возрасту; скоростно-зависимые параметры оцениваются относительно регрессионной нормы, учитывающей темп ходьбы, а результат выражается трёхуровневой шкалой тяжести отклонения. Отличие метода от известных подходов состоит в том, что он, во-первых, использует нормализацию по полной высоте тела вместо диагонали «плечо-бедро» (ASH [19]), вырождающейся при боковом ракурсе; во-вторых, строит нормативные интервалы для безразмерных параметров двумерного видеоанализа, отсутствующие в литературе; в-третьих, вводит скоростно-условную коррекцию нормы для таких параметров, ранее для видеоанализа не реализованную.

ПОДГОТОВКА СКЕЛЕТНЫХ ДАННЫХ И СЕГМЕНТАЦИЯ ШАГОВ

Скелетные данные извлечены моделью AlphaPose в формате COCO-17 и приведены к унифицированному представлению COMMON_13, содержащему двенадцать суставов конечностей и нос как верхнюю опорную точку (pose_dim = 26). Суставный состав включает шесть пар с обеих сторон: плечевые, локтевые, лучезапястные, тазобедренные, коленные и голеностопные. При значении confidence score отдельного сустава ниже 0,3 кадр помечается как неполный и исключается из расчёта масштаба нормализации; порог 0,3 согласуется с практикой работы [9]. Сегментация записи на отдельные шаги выполнена временной свёрточной сетью (Temporal Convolutional Network, TCN) из трёх свёрточных слоёв по 64 канала с размером ядра 15; совокупное рецептивное поле составляет 43 кадра, что соответствует примерно 1,4 с. Модель обучена на ручной разметке фаз шага, разделение выборки по участникам в соотношении 70/15/15. Макроусреднённая F1-мера на валидационной выборке составила 0,88.

Рис. 1. Общая схема конвейера обработки видеозаписи походки

ЭМПИРИЧЕСКОЕ СРАВНЕНИЕ МЕТОДОВ ПРОСТРАНСТВЕННОЙ НОРМАЛИЗАЦИИ

Для систематического сравнения реализованы четыре метода пространственной нормализации, различающиеся выбором базового отрезка скелетной модели. Метод full_height использует вертикальное расстояние от носа до голеностопного сустава опорной ноги; типичное значение ≈258 пикселей. Метод ASH основан на диагональном расстоянии от плечевого до тазобедренного сустава и ранее валидирован для фронтальной проекции [19]. Метод leg_length опирается на расстояние от тазобедренного до голеностопного сустава; теоретическое обоснование такой нормализации восходит к классической работе [18]. Метод torso_vertical использует вертикальное расстояние от носа до центра таза. Стабильность масштаба каждого метода оценивается коэффициентом вариации:

^(s) M^s} ’

cv

где σ (s) — стандартное отклонение покадровых значений масштаба s в пределах записи, M(s) – их медиана. Коэффициент CV вычислялся для каждой записи, далее для каждого участника определялась медиана CV по всем его записям, итоговая статистика агрегирована по 398 участникам.

Нормализация по полной высоте тела обеспечивает наименьший коэффициент вариации масштаба (3,16 %) и оказывается оптимальной для 88,4 % участников, превосходя ближайший метод

Таблица 1. Сравнительные показатели четырёх методов пространственной нормализации (398 участников, латеральная проекция, 30 кадров/с)

Метод

CV, %

Ранг

Доминирование, % участников

full_height

3,16

1

88,4

ASH

3,65

2

10,8

torso_vertical

3,92

3

0,8

leg_length

5,00

4

0,0

ASH на 13,4 % в относительном исчислении. Наблюдаемое превосходство обусловлено двумя факторами. Во-первых, расстояние от носа до голеностопного сустава приближается к реальной вертикальной протяжённости скелета и слабее зависит от фазы шага, поскольку нос и опорный голе-ностоп занимают противоположные вертикальные экстремумы. Во-вторых, максимальная базовая длина пропорционально уменьшает вклад случайного шума детекции координат в коэффициент вариации. Вырождение ASH при боковом ракурсе объясняется проекцией плечевого и тазобедренного сустава на близкие точки. Базовый отрезок становится коротким и шумным, что подтверждает нецелесообразность прямого переноса методики [19] на сагиттальную съёмку. Для построения нормативных таблиц и дальнейшей валидации принят метод full_height.

НАБОР БИОМЕХАНИЧЕСКИХ ПАРАМЕТРОВ И СТАТИСТИЧЕСКИЕ НОРМЫ

На нормализованных координатах рассчитывается 16 скорируемых параметров, сгруппированных в шесть клинических доменов (таблица 2). Доменная группировка опирается на стандартную биомеханическую классификацию, принятую в монографиях [24, 25], и на работы по контролю баланса, вариабельности и симметрии [5, 29]. Дополнительно в полной норм-таблице сохраняется нормированная скорость ходьбы для контекстного отображения без участия в скоринге, а также набор из 12 угловых перцентилей суставных дуг для описания референтных диапазонов.

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

Домен

Параметров

Список

Тип оценки

Пространственновременной

2

нормированная длина шага, каденция

двусторонняя

Амплитуда движений нижних конечностей

3

амплитуда коленного сустава в фазе переноса, амплитуда тазобедренного сустава в фазах переноса и опоры

двусторонняя

Туловище и голова

3

средний наклон туловища, размах наклона туловища, наклон головы вперёд

односторонняя

Стабильность

3

стабильность положения головы, стабильность положения туловища, вертикальная стабильность

односторонняя

Вариабельность

2

коэффициент вариации длины шага, коэффициент вариации каденции

односторонняя

Симметрия

3

индекс симметрии длины шага, индекс симметрии длительности шага, индекс симметрии амплитуды коленного сустава

односторонняя

Амплитуда коленного сустава в фазе опоры исключена из скоринга на этапе агрегации из-за высокого минимально детектируемого изменения (Minimal Detectable Change, MDC), достигающего 90 % диапазона параметра и превышающего клинически значимую разницу. Нормированная скорость ходьбы не скорируется, чтобы избежать циркулярности со скоростно-условной моделью, в которой скорость выступает предиктором.

Построение нормативных интервалов выполнено по стандарту CLSI EP28-A3c [26], совместимому с ГОСТ Р 53022.3-2008 [27]. Стандарт предписывает перцентильную модель с двусторонним интервалом [p2,5; p97,5], охватывающим 95 % центральной массы распределения и не требующим предположения о нормальности. Стратификация выполнена по шести демографическим группам, образованным сочетанием пола и трёх возрастных страт 18–29, 30–49 и 50–64 лет, что согласуется с возрастным диапазоном набора Health&Gait. Для исключения экстремальных выбросов, обусловленных ошибками детекции позы и артефактами сегментации, применён двухуровневый фильтр. На первом уровне отсекаются шаги с физиологически невозможными значениями: длиной менее 10 см, каденцией менее 40 шаг/мин или амплитудой сустава более 120°. На втором уровне применяется правило межквартильного размаха: значения вне интервала [Q1 - 1,5·IQR, Q3 + 1,5·IQR] исключаются из агрегации. Перцентильные интервалы рассчитываются по отфильтрованному распределению; винсоризация на уровне группы применяется отдельно при построении скоростно-условной модели. Агрегация шагов в параметр участника выполнена медианой, робастной к редким артефактам сегментации и корректно учитывающей статистическую зависимость шагов одного индивида.

Ряд параметров шага статистически связан с темпом ходьбы, при этом скорость выступает предиктором для длины шага, каденции и амплитуд тазобедренного сустава в фазах переноса и опоры. Для параметров с абсолютным значением коэффициента корреляции Пирсона со скоростью выше 0,30 применяется регрессионная модель предсказания. В каждой из шести групп обучается линейная регрессия параметра по нормированной скорости v , отклонение оценивается относительно регрессионной полосы p ( v ) ± 2 σ , где p ( v ) — регрессионное предсказание, σ — остаточное стандартное отклонение, вычисляемое отдельно для каждой записи. Нормализованное отклонение задаётся z-оценкой:

(x - p(v))

z = ------------- ,                                           (2)

ст где x – значение параметра у участника.

Итоговая шкала тяжести образована тремя уровнями. Норма соответствует |z| ≤ 2,0 для скоростно-зависимых параметров или значению внутри [p2,5; p97,5] для скоростно-независимых. Лёгкое отклонение фиксируется при 2,0 < |z| ≤ 2,33 или при попадании значения в интервалы [p1; p2,5] и [p97,5; p99]. Значимое отклонение соответствует |z| > 2,33 или значению вне [p1; p99]. Приведённые границы согласованы с рекомендациями CLSI для трёхуровневой клинической классификации. Одностороннее оценивание для параметров доменов «Стабильность», «Вариабельность», «Симметрия» и параметра наклона головы вперёд формализуется отсечением только верхнего хвоста распределения, поскольку их низкое значение не соответствует патологии. Корректность построенной модели проверялась в трёх направлениях — сопоставлением с инструментальным эталоном OptoGait на индивидуальном уровне, сравнением нормативных медиан с независимыми популяционными базами и оценкой диагностических характеристик на клинической когорте.

ВАЛИДАЦИЯ ПАРАМЕТРОВ ОТНОСИТЕЛЬНО OPTOGAIT

Для сопоставления безразмерных параметров с данными инструментальной системы OptoGait выполнена калибровка пиксельного пространства в физические единицы по формуле

H к = Г'                                              (3)

(h ■ с)

где k — калибровочный коэффициент (см/пиксель), H — рост участника (см), h — медианное значение нормализующего отрезка в пикселях, c = 0,92 — эмпирический коэффициент, учитывающий, что верхней опорной точкой в COMMON_13 служит нос, а не макушка. В рабочем режиме метода все параметры остаются безразмерными и не требуют калибровки камеры; формула (3) используется только для валидации.

Источники погрешности для разных параметров различны и оцениваются раздельно. Пространственные параметры (длина шага) наследуют погрешность оценки позы и калибровочного коэффициента c = 0,92, тогда как временные (каденция) определяются главным образом точностью сегментации моментов касания и от пространственной калибровки не зависят. Проекционные искажения бокового ракурса и покадровый шум детекции координат также вносят вклад, различный по каждому параметру, поэтому индивидуальная погрешность оценивалась отдельно для каждой величины: средняя абсолютная ошибка приведена в таблице 3, а надёжность параметров дополнительно контролировалась по минимально детектируемому изменению (MDC) на этапе отбора — так, амплитуда коленного сустава в фазе опоры, для которой MDC достигает 90 % диапазона, исключена из скоринга. Поскольку в рабочем режиме параметры безразмерны, инструментальная погрешность калибровки в результаты не входит и проявляется только при сопоставлении с OptoGait. Воспроизводимость параметров дополнительно оценена по повторным проходам (test-retest): основные параметры (нормированные длина шага и скорость, каденция, наклон головы вперёд, амплитуды тазобедренного сустава, вертикальная стабильность) показали хорошую воспроизводимость (внутриклассовый коэффициент корреляции 0,70–0,95; таблица 6); параметры вариабельности и симметрии показали более низкую воспроизводимость и используются как вспомогательные, интерпретируемые с осторожностью; полученные уровни надёжности согласуются с опубликованными оценками воспроизводимости безмаркерного захвата движения [30, 31].

Для повышения согласованности дисперсии параметров с эталоном авторами разработана процедура VP-коррекции (Variance Preserving): ортогональная регрессия метода и эталона с последующей коррекцией наклона так, чтобы стандартное отклонение скорректированного ряда совпадало с эталонным. Процедура исключает характерное для линейной регрессии сжатие дисперсии. Дополнительно применяется субкадровая интерполяция траектории голеностопного сустава в окрестности локальных минимумов на основе параболической аппроксимации по трём соседним точкам; интерполяция уточняет моменты касания на доли кадра и повышает корреляцию каденции с эталоном.

Таблица 3. Валидация извлекаемых параметров относительно OptoGait (n = 396, UGS). Покрытие p5–p95 рассчитывается как отношение ширины интервала p5–p95 метода к ширине аналогичного интервала эталона в процентах

Параметр

MAE

Относительная

MAE, %

r

Покрытие р5-p95, %

Длина шага

11,13 см

17,6

0,365

102,8

Каденция

13,36 шаг/мин

11,8

0,371

109,8

Скорость ходьбы

0,19 м/с

16,4

0,388

122,9

а) Длина шага, см                   б) Каденция, шаг/мин                в) Скорость ходьбы, м/с

Рис. 2 . Согласованность извлекаемых параметров с показаниями OptoGait.

Ось абсцисс — значение эталона OptoGait, ось ординат — значение метода.

Диагональ y = x — серая пунктирная линия, линия ортогональной регрессии — чёрная сплошная. Цветовая кодировка точек соответствует трём возрастным стратам 18–29, 30–49 и 50–64 лет

Для оценки соотношения полученной точности с результатами других безмаркерных методов видеоанализа походки выполнено сопоставление с опубликованными показателями ближайших по постановке работ (таблица 4).

Таблица 4. Сравнение MAE и коэффициента корреляции с опубликованными результатами безмаркерных методов

Работа

Условия

Pose

MAE (длина шага)

r (длина шага)

Настоящая работа

H&G, UGS, 2D лат.

AlphaPose

11,13 см

0,365

Zafra-Palma и соавт. [28]

H&G, UGS, 2D лат.

OpenPose/AlphaPose

16,47 см

0,39

Stenum и соавт.

[9]

Контрольная выборка, 2D сагит.

OpenPose

11,1 см

0,86

Ota и соавт. [33]

Treadmill

OpenPose/2D

9,0 см

0,82

Kanko и соавт.

[34]

3D markerless, overground

Theia3D

2,2 см

0,97

Вклад отдельных этапов обработки оценён аблационным экспериментом, в котором поочерёдно отключались VP-коррекция и субкадровая интерполяция траектории голеностопного сустава (таблица 5). Введение VP-коррекции повышает корреляцию на 15–23 % относительно базовой конфигурации, при этом наибольший прирост наблюдается для скорости. Дополнительное включение субкадровой интерполяции увеличивает корреляцию каденции с 0,352 до 0,371, что соответствует приросту 21 % относительно базовой конфигурации. Обе компоненты вносят независимый и взаимодополняющий вклад: VP-коррекция улучшает согласование дисперсий, субкадровая интерполяция повышает временную точность.

Таблица 5. Вклад этапов обработки в корреляцию с эталоном OptoGait (n = 396, UGS)

Конфигурация

r (длина шага)

r (каденция)

r (скорость)

Базовая (без коррекции и интерполяции)

0,298

0,307

0,301

С VP-коррекцией

0,344

0,352

0,371

С VP-коррекцией и субкадровой интерполяцией

0,365

0,371

0,388

Отдельно оценена воспроизводимость извлекаемых параметров при повторных проходах; показатели надёжности по каждому параметру приведены в таблице 6.

Таблица 6. Воспроизводимость извлекаемых параметров (test-retest, обычная ходьба): внутриклассовый коэффициент корреляции (ICC) и минимально детектируемое изменение (MDC95)

Параметр

ICC

MDC95

MDC95, % коридора

Нормированная скорость

0,95

0,19

19

Нормированная длина шага

0,92

0,07

24

Наклон головы вперёд, °

0,90

5,5

27

Каденция, шаг/мин

0,82

18,6

36

Амплитуда тазобедренного сустава (опора), °

0,79

3,7

39

Вертикальная стабильность

0,78

0,003

41

Амплитуда тазобедренного сустава (перенос), °

0,70

4,3

44

Амплитуда коленного сустава (перенос), °

0,63

6,1

50

Стабильность головы

0,53

0,002

60

Стабильность туловища

0,53

0,39

58

СОГЛАСИЕ И ПЕРЕНОСИМОСТЬ ПОСТРОЕННЫХ НОРМ

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

На популяционном уровне, в отличие от точности отдельного измерения (раздел 5), проверяются свойства самих построенных норм. Несмещённость центров установлена сопоставлением медиан нормативных параметров с медианами по данным инструментальной системы OptoGait на той же выборке. Для представительной группы мужчин 30–49 лет расхождение медиан не превышает 2,3 % по каждому параметру (таблица 7), а в среднем по всем шести демографическим группам систематический сдвиг не превышает 1 % (длина шага -0,9 %, каденция -0,2 %, скорость +0,6 %). Отсутствие систематического сдвига согласуется с разделом 5: умеренная индивидуальная корреляция (r = 0,365–0,388) и высокое согласие популяционных медиан суть два проявления одного свойства метода — воспроизведения распределения параметров в популяции при заметной случайной составляющей для отдельного участника. Дополнительно покрытие интервалов p5–p95 на уровне 102,8–122,9 % (раздел 5) показывает, что построенные нормы не оказываются ложно-узкими.

Таблица 7 . Согласие нормативных медиан группы мужчин 30–49 лет с инструментальной системой OptoGait (обычная ходьба)

Параметр

Наши

OptoGait

Расхождение, %

Длина шага, см

55,5

56,8

-2,2

Каденция, шаг/мин

114,3

114,0

+0,3

Скорость, см/с

157,5

154,0

+2,3

Переносимость норм проверена сравнением нормативных центров, независимо построенных по двум наборам данных — Health & Gait и клинической когорте KOA-PD-NM, зарегистрированной иной камерой и иной моделью оценки позы: отношение медиан каденции между двумя базами составило 0,98, безразмерной длины шага — 0,94. Поскольку сравнение выполнено на безразмерных величинах, оно не зависит от калибровки камеры и указывает, что построенные нормы отражают свойство походки, а не артефакт конкретной съёмки. Сравнение проведено на двух базах, поэтому подтверждает переносимость в пределах проверенных условий, а не универсальность норм; по порядку величины построенные нормы согласуются с опубликованными нормативами пространственно-временных параметров [20, 21, 22, 23].

ДИАГНОСТИЧЕСКАЯ ВАЛИДАЦИЯ НА НЕЗАВИСИМЫХ КОГОРТАХ

Таблица 8. Диагностическая валидация статистического критерия на независимых когортах

Когорта

n

Тип оценки

Уровень тяжести

Значение

Health&Gait LOO (здоровые, внутренняя)

379

FPR

значимое

1,1 % (4/379)

Health&Gait LOO (здоровые, внутренняя)

379

FPR

лёгкое и выше

51,4 % (195/379)

KOA-PD-NM NM (здоровые, внешняя)

59

FPR

значимое

15,3 % (9/59)

KOA-PD-NM, клинические патологии

128

Sens

значимое

40,6 % (52/128)

KOA-PD-NM, клинические патологии

128

Sens

лёгкое и выше

99,2 % (127/128)

Чувствительность на значимом уровне равна 40,6 %, что отражает консервативность порогов [p2,5; p97,5] по CLSI, интерпретирующих как значимое только отклонение за 95 % референтного диапазона. Существенно более информативна интегральная чувствительность: у 127 из 128 клинических пациентов (99,2 %) детектировано хотя бы одно лёгкое отклонение, то есть метод практически не пропускает патологию, но требует клинической интерпретации профиля для принятия решения о её значимости. Высокая интегральная чувствительность при умеренной специфичности на значимом уровне типична для скринингового инструмента и соответствует целевой области применения метода.

ПРИМЕР ИНТЕРПРЕТАЦИИ ПРОФИЛЯ ОТКЛОНЕНИЙ

Возможности структурированной интерпретации иллюстрирует клинический случай из когорты KOA-PD-NM. Мужчина 58 лет, диагноз «эссенциальный тремор» установлен независимой неврологической экспертизой; запись выполнена продолжительностью 12 с, после сегментации получено 27 валидных шагов. Итоговый профиль приведён в таблице 9; параметры доменов «Пространственновременные» и «Симметрия» остались в пределах нормативного интервала и в таблице не отображены.

Три параметра домена «Стабильность» классифицированы как значимые отклонения с z-оценками 2,5–3,1, что согласуется с клинической картиной усиления постуральных колебаний при эссенциальном треморе. Наклон головы вперёд и коэффициент вариации каденции получили статус лёгкого отклонения при сохранённой амплитуде движений в коленных суставах. В отличие от нейросетевых классификаторов типа «чёрный ящик» [35, 36], метод раскрывает, какие именно параметры отклонены и в каком направлении, что повышает интерпретируемость результатов для клинициста.

Таблица 9 . Профиль отклонений участника с эссенциальным тремором (группа мужчин 50–64 лет). Параметры стабильности — безразмерные величины, выражающие стандартное отклонение нормированной координаты за цикл шага; коэффициент вариации каденции приведён в процентах, углы — в угловых градусах

Домен

Параметр

Значение

Граница нормы

z-оценка

Тяжесть

Стабильность

стабильность положения головы

0,091

p97,5 = 0,065

+3,1

значимое

Стабильность

стабильность положения туловища

0,072

p97,5 = 0,054

+2,7

значимое

Стабильность

вертикальная стабильность

0,038

p97,5 = 0,029

+2,5

значимое

Амплитуда

амплитуда коленного сустава в фазе переноса

57,2°

p2,5 = 54,1°

-0,8

норма

Туловище и голова

наклон головы вперёд

12,4°

p97,5 = 9,8°

+2,1

лёгкое

Вариабельность

коэффициент вариации каденции

5,8 %

p97,5 = 4,9 %

+2,3

лёгкое

ОГРАНИЧЕНИЯ И ПЕРСПЕКТИВЫ

Индивидуальная корреляция r = 0,37–0,39 отвечает физическому пределу 2D-анализа при сравнении с контактными системами [28], поэтому метод применим для популяционного скрининга и мониторинга, но не для окончательной клинической диагностики. Валидация выполнена только для бокового ракурса; применение фронтальной проекции требует иной нормализации (например, ASH [19]) и отдельной валидации. Возрастной диапазон 18–64 лет определён границами выборки Health&Gait; построение нормативных интервалов для когорты старше 64 лет на независимой выборке составляет приоритетную задачу дальнейших исследований. Расширение набора параметров за счёт угловых скоростей и фазовых ошибок цикла шага затруднено шумом второй производной при 30 кадрах/с и потребует повышения частоты съёмки или адаптивных сглаживающих фильтров. Перспективным направлением развития метода выступает верификация нормативных интервалов на выборке с клинически подтверждённым статусом здоровья и расширение клинической когорты KOA-PD-NM за счёт других нозологий, характерных для поликлинической практики анализа походки.

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

ЗАКЛЮЧЕНИЕ

Предложенный метод извлечения биомеханических параметров походки из двумерного скелетного видеоанализа и их статистической оценки относительно стратифицированных популяционных норм объединяет пространственную нормализацию по полной высоте тела, набор из 16 скорируемых параметров в шести клинических доменах, перцентильные нормативные интервалы по CLSI EP28-A3c и ГОСТ Р 53022.3-2008 для шести демографических страт, а также скоростно-условную оценку с трёхуровневой шкалой тяжести отклонений.

Эмпирическое сравнение методов нормализации на представительной выборке подтвердило превосходство выбранного подхода по критерию стабильности масштаба. Валидация относительно инструментальной системы OptoGait показала, что индивидуальная точность находится в диапазоне, опубликованном составителями набора Health&Gait, при этом абсолютная погрешность длины шага существенно ниже авторской за счёт VP-коррекции и субкадровой интерполяции. Согласованность нормативных медиан с инструментальной системой OptoGait и переносимость норм между двумя независимыми наборами данных подтвердили корректность калибровки распределений на популяционном уровне. Диагностическая валидация на независимой клинической когорте продемонстрировала высокую интегральную чувствительность при сохранении приемлемой специфичности в условиях, отличающихся от условий сбора обучающей выборки. Отсутствие необходимости в калибровке камеры и антропометрических измерениях делает метод пригодным для популяционного скрининга и мониторинга нарушений походки вне специализированной лаборатории; структурированный выход метода, раскрывающий конкретные отклонённые параметры и направление отклонения, обеспечивает интерпретируемость результатов для клинициста.