Остеоартрит (ОА) представляет собой чрезвычайно распространённое заболевание [1], которое наносит существенный урон на здоровью и благосостоянию населения большинства развитых стран [2, 3]. Современная парадигма рассматривает целый ряд вопросов диагностики и лечения ОА. Одна из ключевых проблем диагностики ОА коленных суставов (КС) состоит в отсутствии надёжных визуальных маркеров, отражающих комплекс суставных изменений [3, 4].
Современные рекомендации предлагают использовать стандартную рентгенографию в качестве основного средства визуализации изменений, происходящих в суставе при ОА [5]. Вместе с тем рентгенография обладает ограниченной чувствительностью и специфичностью, которая зависит главным образом от квалификации рентгенолога и имеет описательный характер [4]. Это обстоятельство диктует потребность идентификации новых маркеров и способов оценки рентгенограмм при ОА.
Фрактальный анализ – перспективный метод обработки и интерпретации данных, активно использующийся в различных разделах биомедицины, таких как неврология, микробиология, нанотехнология, фармакология [6]. Фрактальная размерность применима в тех случаях, когда имеется математическое множество, обладающее свойством самоподобия, нетривиальной структурой во всех масштабах и обладает дробной метрической размерностью. Очевидно, что параметры суставной щели обладают признаками фракталов. В связи с этим измерение фрактального размера суставной щели представляет собой интерес с точки зрения надёжности [7] и независимости от целого ряда факторов, таких как вариации положения пациента [8], яркость и контрастность снимка [9], а также способа получения рентгеновского изображения [10].
Рентгеновские изображения обладают ещё одним нетривиальным свойством – неоднородностью. Существуют различные способы идентификации инвариантности, среди которых лакунарный анализ фрактальных текстур [11], градиентные текстурные показатели [12]. В доступной литературе нами не было найдено информации об использовании данных методов для оценки изображений суставной щели при ОА КС.
Цель исследования: дать оценку вариабельности фрактального размера и показателей текстурных параметров суставной щели при остеоартрите с данными стандартной рентгенографии коленного сустава.
Материалы и методы исследования
В исследование были включены 92 больных (78 женщин и 14 мужчин) в возрасте 66,1 ± 10,5 лет с ОА КС I-IV рентгенологических стадий по Kellgren-Lowrence. В группу сравнения были включены 24 здоровых добровольца (16 женщин и 8 мужчин) в возрасте 20–34 (29,6 ± 5,96) лет без признаков остеоартрита. Характеристика групп пациентов приведена в табл. 1. Всеми пациентами было подписано добровольное информированное согласие на участие в данном исследовании, протокол исследования был одобрен междисциплинарным комитетом по этике ФГБОУ ВО ТГМУ Минздрава России.
Критерии исключения: посттравматический гонартроз, длительная иммобилизация, тяжелая соматическая патология, отсутствие согласия на участие в исследовании.
Таблица 1
Клиническая характеристика групп пациентов
|
Показатель |
Группа ОА |
Группа сравнения |
|
Всего, n, абс. |
92 |
24 |
|
Пол ж/м, n, абс. |
78/14 |
16/8 |
|
Возраст, M ± SD |
66,1 ± 10,5 |
29,6 ± 5,96 |
|
Стадия ОА I/II/III–IV, абс. |
14/52/26 |
– |
В рамках исследования было проведено всестороннее обследование пациента опытным врачом-ревматологом, клиническое и лабораторное обследование, выполнялась рентгенография коленных суставов по общепринятой методике в прямой проекции на аппарате «КРТ ОКО Электрон». Рентгенограммы описывал опытный врач-рентгенолог, устанавливая стадию ОА. Рентгеновские изображения обрабатывали в режиме «pixel to pixel» в программах DICOM и ImageJ. Выбирали область интереса, включающую суставную щель с захватом
остеофитов, плато большеберцовой кости и мыщелков бедренной кости по межмыщелковое возвышение.
Сегментацию аналогового изображения осуществляли по порогу автоматической бинаризации. Бинарные изображения подвергали фрактальному анализу box-методом, получали фрактальный размер (Db). Для оценки инвариантности текстуры (Sλ) суставной щели использовали показатель графика наклона прямой коэффициента регрессии натуральных логарифмов лакунарности от масштаба измерения. Изучение текстурных вариаций суставной щели КС осуществляли с использованием аналитического протокола матриц совпадений градаций серого (GLCM – Grey Level Co-occurrence Matrix): получали цифровые выражения второго углового момента (ASM), обратного момента разностей (IDM), энтропии (E).
Статистическая обработка результатов проводилась в программе Stаtistica 6.0 (StatSoft, США). Для сравнении двух групп использовали Z-критерий Манна – Уитни, для сравнения более двух групп использовали H-критерии Краскела – Уоллиса. Результаты представлены в виде среднего и стандартного отклонения (M ± SD). Корреляционный анализ проводился по методу Спирмена. Различия считали статистически значимыми при p < 0,05.
Результаты исследования и их обсуждение
У больных ОА, включенных в исследование, ФР суставной щели коленного сустава обратно коррелировала со стадией заболевания (r = –0,4, p = 0,004) и большим размером краевых остеофитов (r = –0,3, p = 0,02), а также с возрастом пациентов (r = –0,3, p = 0,03). Показатель инвариантности суставной щели Sλ продемонстрировал прямую статистически значимую связь со стадией ОА (r = 0,5, p = 0,00007), наличием больших краевых остеофитов (r = 0,4, p = 0,004) и возрастом пациентов (r = 0,4, p = 0,004). Результаты текстурного анализа по протоколу GLCM не продемонстрировали корреляционных связей с рентгенологической стадией гонартроза, остеофитами и возрастом пациентов.
У больных ОА отмечено достоверное снижение Db по сравнению с группой лиц без ОА (Z = –2,0, p = 0,04). Наблюдали значимое повышение Sλ у пациентов с гонартрозом по сравнению с интактной группой (Z = 2,1, p = 0,04). С другой стороны, показатели текстурного анализа не продемонстрировали статистически значимых различий между группами ОА и сравнения (табл. 2). У лиц обоих полов в группе ОА текстурные параметры не имели статистически значимых различий.
Таблица 2
Фрактальные и текстурные характеристики изображений суставной щели коленного сустава в исследуемых группах
|
Параметр |
Группа ОА |
Группа сравнения |
||||
|
Общее n = 92 |
Женщины n = 78 |
Мужчины n = 14 |
Общее n = 24 |
Женщины n = 16 |
Мужчины n = 8 |
|
|
Db |
1,679 ± 0,060 |
1,676 ± 0,056 |
1,698 ± 0,081 |
1,703 ± 0,029* |
1,696 ± 0,033 |
1,716 ± 0,017 |
|
Sλ |
–0,245 ± 0,059 |
–0,245 ± 0,062 |
–0,240 ± 0,037 |
–0,279 ± 0,019* |
–0,284 ± 0,015 |
–0,268 ± 0,024 |
|
ASM |
0,479 ± 0,078 |
0,489 ± 0,038 |
0,423 ± 0,181 |
0,494 ± 0,021 |
0,498 ± 0,11 |
0,485 ± 0,034** |
|
IDM |
0,936 ± 0,099 |
0,951 ± 0,032 |
0,849 ± 0,237 |
0,959 ± 0,018 |
0,965 ± 0,006 |
0,947 ± 0,027** |
|
E |
0,984 ± 0,037 |
0,844 ± 0,084 |
1,767 ± 0,388 |
0,831 ± 0,053 |
0,810 ± 0,16 |
0,874 ± 0,079** |
Примечания: * – межгрупповые различия статистически значимы, p < 0,05; ** – внутригрупповые различия статистически значимы, p < 0,05.
У добровольцев без ОА Db и Sλ также не продемонстрировали значимых гендерных различий. В то время как показатели GLCM проявили дискриминационные признаки по половому признаку. Так, ASM и IDM были статистически значимо выше у женщин (соответственно: Z = 1,9, p = 0,04; Z = 2,5, p = 0,01). Энтропия была выше у лиц мужского пола, не имевших ОА (Z = –2,5, p = 0,01). Среди больных гонартрозом не удалось установить статистически значимых различий не по одному изучаемому параметру между группой сравнения и пациентами с I стадией ОА КС. Вместе с тем Db статистически значимо снижался по мере прогрессирования гонартроза (H = 8,2, p = 0,04, рис. 1), а инвариантность Sλ увеличивалась (H = 16,8, p = 0,0008, рис. 2). Больные ОА с большими краевыми остеофитами показатели Db и Sλ имели статистически значимые различия (соответственно Z = 2,1, p = 0,04; Z = –2,7, p = 0,007).
Рис. 1. Показатели фрактальной размерности суставной щели в исследуемых группах. Пояснения в тексте
Рис. 2. Показатели инвариантности суставной щели в исследуемых группах. Пояснения в тексте
Оценка рентгенограмм до настоящего времени носит преимущественно описательный характер, что нередко влечёт диагностические ошибки [13]. В настоящем исследовании показано, что суставная щель имеет нетривиальную и инвариантную форму, то есть обладает фрактальными свойствами, что может быть использовано для определения её размерностей.
Наши результаты показывают, что фрактальный размер суставной щели имеет обратную связь со стадией гонартроза. По мере прогрессирования рентгенологических симптомов гонартроза происходит уменьшение фрактального значения суставной щели, что позволяет отнести фрактальный размер к достоверным визуальным маркерам заболевания. Этим данным противоречат результаты исследования T. Miyamura и соавт. (2002), в котором фрактальные размерности несегментированных изображений не продемонстрировали значимых связей со стадией ОА коленных суставов [14]. В другом более раннем исследовании бинарных изображений фрактальный размер суставной щели был связан со степенью её сужения [15]. Это диктует необходимость разработки строгих условий протокола исследований с детализацией процедуры сегментации, например, бинаризации, которая была успешно применена в нашей работе.
Изображение суставной щели имеет ещё одно фрактальное свойство – неупорядоченность или лакунарность. Анализ инвариантных свойств бинарных изображений коленного сустава показал, что по мере усиления рентгенологических симптомов изменяется её неоднородность, что значимо отражает коэффициент наклона графика прямой лакунарности от масштаба измерения.
Показатели самоподобия и неоднородности позволяют дискриминировать пациентов с и без гонартроза. Более того, данные показатели склонны изменяться по мере прогрессирования симптомов и связаны с развитием краевых остеофитов. Полученные результаты согласуются с E.A. Messent и соавт. (2007), которые убедительно доказали возможность фрактальной размерности идентифицировать остеофиты у пациентов с ОА КС [15]. Однако в данном исследовании использовались не сегментированные изображения, а также не уточнялась взаимосвязь со стадией ОА. Данные обстоятельства позволяют наметить предпосылки для решения важных задач автоматической оценки рентгенологических симптомов ОА и принятия решения.
Текстурные характеристики суставной щели, вычисляемые по сопряжённости уровней яркости на изображении, не показали статистически значимых различий у пациентов исследуемых групп. По-видимому, это объясняется с тем, что сопряженность уровней градации серого определяется на основе пар смежных пикселей, разделенных заданным вектором и расстоянием. В нашем случае были использованы бинарные изображения, имеющие узкий спектр гистограммы оттенков серого, что в значительной степени ограничивало запись ячеек матрицы, большинство из которых были пустыми. Это может свидетельствовать и о том, что текстурные характеристики, основанные на оценке градаций оттенков серого в бинарном изображении, не отражают изменений пространственной организации суставной щели, что вполне предсказуемо. Параметры GLCM текстурного анализа показали возможность идентифицировать половые различия в характеристиках изображений суставной щели у лиц без ОА, что может быть использовано в решении узкоспециализированных исследовательских задач и, возможно, требует уточнения на больших выборках.
Выводы
Таким образом, нами установлены значимые показатели, отражающие изменения пространственной организации суставной щели при гонартрозе. Фрактальный размер и инвариантность суставной щели продемонстрировали возможность диагностировать гонартроз, наблюдать изменения, происходящие по мере прогрессирования сужения суставной щели и роста остеофитов. Данные обстоятельства ещё раз подтверждают самоподобие текстур суставной щели, которые возможно описать методами неевклидовой геометрии. С другой стороны, даются предпосылки для создания автоматизированных систем принятия решений и идентификации рентгенологических изменений при гонартрозе.