$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
В этом разделе представлен комплексный подход, разработанный для улучшения диагностики и классификации остеоартрита коленного сустава с использованием модифицированной модели XceptionNet. Представленная методология основана на тщательной предварительной обработке данных, тщательной настройке архитектуры модели и надежных методах оценки, которые предназначены для решения сложных проблем, связанных с визуализацией остеоартрита коленного сустава. На рисунке 2 проиллюстрирован ход модели.

Рисунок 2: Рабочий процесс модели. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Описание набора данных
Набор данных, использованный в этом исследовании, включает 9 786 рентгеновских изображений коленного сустава из набора данных OAI, которым была присвоена оценка OA на основе метода оценки KL. Этот набор данных, который предлагает широкий спектр примеров, охватывающих различные фазы KOA, имеет важное значение как для обучения, так и для оценки предлагаемой модели24. Для масштабирования всех изображений использовалась билинейная интерполяция. Этот метод был выбран потому, что он обеспечивает компромисс между вычислительной эффективностью и поддержанием качества изображения, которые необходимы для поддержания анатомических особенностей, важных для категоризации остеоартрита коленного сустава. В наборе данных оценки расположены от 0 до 4, где степень 0 означает здоровое колено, а степень 4 обозначает тяжелый остеоартрит. Эта степень предлагает сложный диапазон эволюции заболевания, что важно для модели для понимания мельчайших вариаций на разных стадиях.
Кроме того, набор данных был организован в тестовые, проверочные и обучающие наборы. Для оценки производительности модели и обеспечения того, чтобы убедиться, что она учится эффективно обобщать без переобучения обучающей выборки, такое разделение было необходимо. Тестовый набор позволяет объективно оценить эффективность окончательно доработанной модели, так как он также был использован для модификации весов модели. Кроме того, валидационный набор помогает в тонкой настройке гиперпараметров и оценке моделей на протяжении всего этапа обучения. Краткая сводка распределения данных представлена в таблице 2, а на рисунке 3 представлено ее визуальное изображение.
| Класс | Автотест | Тест | Поезд | Валь |
| 0 | 604 | 639 | 2286 | 328 |
| 1 | 275 | 296 | 1046 | 153 |
| 2 | 403 | 447 | 1516 | 212 |
| 3 | 200 | 223 | 757 | 106 |
| 4 | 44 | 51 | 173 | 27 |
Таблица 2: Краткое описание набора данных.

Рисунок 3: Распределение набора данных. Набор данных был организован в тестовый, валидационный и обучающий наборы. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этой цифры.
Предварительная обработка данных
Этот шаг имеет решающее значение и направлен на повышение производительности модели за счет повышения качества входных данных. Этапы предварительной обработки следующие.
Изменение размера изображения: К каждому рентгеновскому изображению применяется стандартный размер, чтобы обеспечить согласованность всего набора изображений. Этот этап необходим для того, чтобы модель могла последовательно и эффективно обрабатывать изображения. В этом процессе используется уравнение 1.
Новый размер изображения = Изменить размер (исходное изображение, стандартный размер) (1)
Нормализация: Значения пикселей на изображениях были стандартизированы таким образом, чтобы получить среднее значение, равное нулю, и стандартное отклонение, равное единице. Такой вид нормализации гарантирует, что входные значения попадают в один и тот же диапазон, что ускоряет сходимость модели во время обучения. Уравнение 2 используется для выполнения этого нормализующего шага.
(2)
Аугментация данных: Стратегии аугментации данных использовались для решения проблем, связанных с недостаточным объемом данных, и повышения способности модели к обобщению. Эти методы включают в себя перевороты, масштабирования, повороты и перемещения. Они вводят новые изображения, которые намеренно добавляются в обучающий набор данных. Уравнение 3 используется при реализации данной процедуры аугментации.
Дополненное изображение = Применить трансформацию (оригинальный iamge) (3)
Для обеспечения вариативности обучающей выборки с целью повышения надежности модели и минимизации переобучения были использованы различные методы аугментации. Диапазон сдвига 0,2, диапазон масштабирования от 0,8 до 1,2 (что соответствует коэффициенту масштабирования 0,2) и диапазон случайного вращения от -30° до +30° были включены в аугментацию обучающих данных. Для улучшения обобщения модели допускались как горизонтальные, так и вертикальные перевороты с вероятностью 0,5. Для вертикального и горизонтального сдвига также использовался коэффициент сдвига 0,1, или ± 10% искажения.
Балансировка классов: Алгоритмы взвешивания классов использовались из-за возможного дисбаланса классов в наборе данных (с большим количеством изображений одних классов открытого доступа, чем других). Чтобы модель не показывала смещение в пользу более распространенных классов, этот метод присваивает недопредставленным классам большие веса. В этом процессе используется уравнение 4. Для определения веса использовались следующие значения: Здоровый: 3857, Сомнительный: 1770, Минимальный: 2578, Умеренный: 1286 и Тяжелый: 295 - это распределения по классам 9786 всего выборок и 5 классов.
(4)
Разделение обучения и валидации: набор данных был организован в наборы для проверки и обучения с использованием традиционных соотношений 80:20. Разделение данных позволяет модели учиться на обучающих данных и регулярно оценивать ее производительность с помощью проверочного набора. Это позволяет выявить переобучение и оценить способность модели к обобщению.
Оптимизация конвейера данных: Был настроен улучшенный конвейер данных, гарантирующий эффективную обработку данных во время обучения модели. Чтобы уменьшить узкие места ввода-вывода и повысить производительность вычислений, в этом конвейере используются методы пакетной обработки, предварительной выборки и параллельного извлечения данных. Для повышения эффективности и производительности обучения был разработан конвейер входных данных. Чтобы гарантировать сбалансированное использование памяти и вычислительную нагрузку, был использован размер пакета 32. Для повышения эффективности ввода данных в функции map было использовано значение num_parallel_calls=4 для распараллеливания загрузки и предварительной обработки данных. Кроме того, чтобы свести к минимуму задержку ввода-вывода, используется prefetch(buffer_size=tf.data.AUTOTUNE) TensorFlow для автоматической регулировки размера буфера предварительной выборки. Это фактически перекрывало предварительную обработку данных и выполнение модели
На рисунке 4 показано несколько экземпляров различных классов изображений после базовой предварительной обработки.

Рисунок 4: Экземпляры из набора данных. На рисунке показано несколько экземпляров различных классов изображений после базовой предварительной обработки. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этой цифры.
Благодаря этим тщательным шагам предварительной обработки набор данных превращается в надежную основу для обучения модифицированной модели XceptionNet, создавая основу для точного и надежного обнаружения и классификации KOA.
Архитектура модели
Базовая модель: Благодаря своей уникальной конструкции, в которой используются разделяемые по глубине свертки для достижения наилучшего возможного баланса между вычислительной эффективностью и сложностью модели, модель XceptionNet, созданная Франсуа Шолле, выделяется25. Уравнение 5 было использовано для реализации сепарабельных сверток по глубине.
Глубинно разделимая свёртка(x) = Точечно (Глубина(x)) (5)
Это особенно актуально для данного исследования по нескольким причинам. Во-первых, его архитектура известна тем, что позволяет извлекать подробные и иерархические характеристики из изображений, что имеет решающее значение для анализа медицинских изображений, где тонкие особенности могут указывать на различные стадии состояния. Во-вторых, предложенная модель основана на архитектуре XceptionNet, которая продемонстрировала выдающуюся эффективность в нескольких задачах классификации изображений, направленных на обнаружение и классификацию остеоартроза коленного сустава26.
Представленный OsteoXceptionNet является специализированной версией архитектуры XceptionNet, созданной специально для автоматической классификации остеоартрита коленного сустава по рентгеновским снимкам. Важные изменения включают в себя добавление еще трех сверточных блоков (фильтры: 512, 256 и 128) для улучшения извлечения признаков, специфичных для структур коленного сустава, а также сверточных слоев с размерами фильтров 3 x 3 и переменным количеством фильтров в диапазоне от 32 до 1024 для входного, среднего и выходного потоков. Пакетная нормализация и активация ReLU происходят после каждого сверточного слоя, а для минимизации переобучения добавляются слои дропаута. Для ввода модели используется стандарт 224 x 224 пикселя.
Тонкая настройка и модификации
Модификации: Модель XceptionNet была стратегически модифицирована в нескольких направлениях, чтобы сделать ее пригодной для классификации и обнаружения KOA.
Тонкая настройка: Верхние слои предварительно обученной модели XceptionNet были разморожены, что позволило модели изучить высокоуровневые характеристики, специфичные для рентгеновских изображений коленного сустава. Тонкая настройка модели проводилась по уравнению 6.
(6)
где θfine-tuned - параметры после тонкой настройки; θпредварительно обученные параметры – это предварительно обученные параметры; λ – скорость обучения; ∇θпредварительно обученныйLoss — это градиент функции потерь, связанный с параметрами, которые были предварительно обучены.
Дополнительные сверточные слои: Дополнительные сверточные слои были введены на заключительном этапе модели. Эти слои имеют меньший размер ядра, чтобы сосредоточиться на извлечении более мелких деталей, относящихся к различным классам открытого доступа. Он был использован с использованием уравнения 7.
Сверточный выход = свертка (вход, ядро) (7)
Пакетная нормализация: После каждого добавленного сверточного слоя применяется пакетная нормализация для стабилизации обучения и увеличения скорости сходимости. Он используется с помощью уравнения 8.
(8)
где
— нормализованный выход; x — вход слоя пакетной нормализации; μ — среднее значение для пакета входных данных; σ2 — дисперсия пакета входных данных; ε — минимальная константа для численной стабильности.
Функции активации: Функции активации ReLU используются в дополнительных слоях для введения нелинейности, тем самым позволяя модели приобретать более сложные структуры в данных. ReLU рассчитывается с помощью уравнения 9.
ReLU(x) = max(0,x) (9)
где, x — вход функции активации ReLU.
Dropout: Для предотвращения переобучения включаются слои dropout, особенно после вновь добавленных слоев, чтобы обеспечить генерализацию модели даже на невидимых данных. Это было сделано с помощью уравнения 10.
Выход = Вход × Маска (10)
Конфигурация выходного слоя
Выходной слой: Исходная модель XceptionNet структурирована для многоклассовой классификации с функцией активации Softmax в выходном слое. В этой адаптированной модели выходной слой настроен таким образом, чтобы представить пять степеней остеоартрита коленного сустава, начиная от 0 до 4 степени, т.е. от здорового до тяжелого уровня. В частности, этот слой состоит из пяти нейронов, каждый из которых соответствует одной из степеней ОА. Функция активации Softmax была использована в этом слое для получения распределения вероятностей по пяти классам, что позволило предложенной модели прогнозировать соответствующую степень ОА для данного рентгеновского снимка коленного сустава. Этот подход согласовывает выходные данные модели непосредственно со шкалой клинической оценки, облегчая интуитивную и практическую интерпретацию прогнозов модели для медицинских работников.
Архитектура XceptionNet специально адаптирована, которая была предварительно обучена на наборе данных ImageNet, для целей представленного исследования по классификации тяжести остеоартрита коленного сустава. Модель использовала изученные признаки из большого и разнообразного набора данных с использованием предварительно обученных весов от ImageNet, что дало ей прочную основу. Чтобы настроить фундаментальную модель XceptionNet для предполагаемых исследований, было добавлено множество новых слоев.
Во-первых, двухмерные сверточные слои были объединены с различным количеством фильтров и размеров ядер, каждый из которых был активирован ReLU и пакетной нормализацией для добавления нелинейности. Эти дополнительные слои были предназначены для захвата и усиления особенностей, связанных с тяжестью остеоартрита коленного сустава. Затем модель была готова к последнему этапу классификации путем добавления слоя Global Average Pooling 2D (GAP), который сжимает карту признаков, сохраняя при этом важные данные. В частности, операция GAP обеспечивает уменьшение размерности с минимальными потерями данных за счет уменьшения пространственных размерностей каждой карты признаков до одного значения при сохранении глубины, которая может варьироваться от 7 x 7 x 1024 до 1 x 1 x 1024. Выходные данные слоя GAP подавались в функцию активации Softmax, которая обеспечивала многоклассовую классификацию путем преобразования логитов в вероятности. Уравнения 11, 12, 13, 14, 15, 16 и 17 используются для GAP, активации Softmax, вычисления карты признаков, снижения скорости обучения, компиляции модели, корректировки веса классов и потери многоклассовой классификации соответственно.
(11)
где, x,i,j – активация i-й строки карты признаков и j-го столбца; H — высота карты объектов; W – ширина карты объектов.
(12)
где, xi — вход в функцию Softmax для класса i; n — количество классов.
Карта признаков = σ(Свертка(Ввод, Ядро) + Смещение) (13)
где, Свёртка — операция свёртки; i - функция активации; Input — входной тензор к слою; Ядро — это сверточное ядро; Предвзятость — это термин предвзятости.
Новый коэффициент обучения = коэффициент × коэффициент обучения (14)
Model = Compile(Архитектура,Потери,Оптимизатор,Метрики) (15)
(16)
где,Класс веса — вес, присвоенный классу; Общее количество выборок — это total_number_of_samples набора данных; Number_of_Classes — количество отличительных классов в наборе данных; Выборка в классе — это количество выборок в определенном классе.
(17)
где, y,i — фактическое распределение вероятностей для класса I; pi — предсказание распределения вероятностей класса I ; N — количество классов.
Наконец, модель, теперь включающая модифицированную архитектуру XceptionNet с пользовательскими слоями, была скомпилирована для прогнозирования вероятностей для каждого класса. С помощью этих модификаций цель состояла в том, чтобы расширить возможности модели по выявлению нюансов, связанных с различными степенями тяжести остеоартрита коленного сустава, улучшив ее классификационные характеристики для конкретной задачи.
Для стабильной и эффективной сходимости оптимизатор Адама был использован для обучения модели с коэффициентом обучения 0,0001, при этом 0,5 использовался в качестве коэффициента отсева для минимизации переобучения. Для ограничения сложности модели и улучшения обобщения были использованы стратегии регуляризации L1 и L2. Для задач многоклассовой классификации с целочисленными метками подходящей функцией потерь была разреженная категориальная кроссэнтропия. Обучение проводилось на протяжении 50 эпох. Softmax, конечная функция активации, использовалась для генерации распределений вероятностей классов. Кроме того, для обучения использовался 64-ступенчатый размер.
Обучение: Процесс обучения был очень важным этапом. На этом этапе модифицированная модель XceptionNet обучается точному обнаружению и классификации KOA на рентгенографических изображениях. Для реализации модели XceptionNet были использованы Keras версии 2.6.0 и бэкенд TensorFlow версии 2.6.0. Ниже приведена подробная информация о ключевых компонентах этапа обучения, включая функцию потерь, оптимизатор, обратные вызовы, размер пакета и эпохи.
Loss_Function: Для задачи многоклассовой классификации использовалась категориальная cross_entropy. Эта функция потерь особенно подходит для задач, в которых ожидается, что каждому экземпляру будет присвоена одна и только одна метка из набора категорий27. Он оценивает производительность модели, получая оценку вероятности в пределах от нуля до единицы. Эта потеря была использована в качестве индикатора для эффективного обучения этой модели для точного прогнозирования серьезности KOA, поскольку она увеличивается, когда прогнозируемая вероятность отличается от фактической метки.
Оптимизатор: Был использован оптимизатор Адама, который хорошо известен своей эффективностью и динамическими характеристиками скорости обучения. Адам объединяет преимущества алгоритмов AdaGrad и RMSProp, предлагая алгоритм оптимизации, способный управлять разреженными градиентами в условиях зашумленной задачи28.
Основные параметры для оптимизатора Adam
Коэффициент обучения: Коэффициент обучения был использован в размере 0,0001, что позволило оптимизатору первоначально внести существенные корректировки в весовые коэффициенты, тем самым оптимизируя процесс обучения.
Beta1 и Beta2: Эти параметры регулируют скорость, с которой квадрат градиента и скользящие средние предыдущих градиентов уменьшаются соответственно. Используются значения по умолчанию 0,9 для beta1 и 0,999 для beta2.
Epsilon: Этот параметр предотвращает любое деление на ноль в реализации, установленное в маленькое число, близкое к нулю.
Обратные вызовы: Обратные вызовы используются во время обучения для контроля производительности модели и для корректировки. Были задействованы следующие обратные звонки:
Early_Stopping: Используется для мониторинга потери валидации модели, а также для остановки процесса обучения, если потери перестают уменьшаться в течение заданного количества эпох (это называется терпением). Когда данные валидации больше не показывают улучшения производительности модели, они останавливают процесс обучения, что помогает предотвратить переобучение.
Reduce LROn Plateau: Этот обратный вызов снижает скорость обучения, когда потери при проверке перестают улучшаться, что позволяет более тонко корректировать веса, что потенциально может привести к повышению общей производительности модели29.
Контрольная точка модели: Этот обратный вызов сохраняет модель через определенный интервал, чтобы после завершения процесса обучения можно было получить наилучшую версию модели. Как правило, он отслеживает точность или потери проверки и сохраняет вес модели при обнаружении улучшения.
Размер партии и эпохи:
Размер партии: 32 — это типичный размер партии, который использовался для балансировки требований к стабильности сходимости модели и эффективности вычислений. Размер партии 32 обеспечивает баланс, будучи достаточно большим для использования вычислительной оптимизации и в то же время достаточно маленьким, чтобы обеспечить стабильную оценку градиента.
Эпохи: модель настроена на обучение в течение до 50 эпох, хотя обучение может быть остановлено досрочно, если активируется обратный вызов EarlyStoping. 50 эпох обеспечивают достаточное количество итераций для корректировки весов и сходимости модели, в то время как EarlyStopping гарантирует, что обучение не будет продолжаться без необходимости.
Благодаря тщательной настройке этих параметров и использованию обратных вызовов процесс обучения был оптимизирован, чтобы модель эффективно обучалась и хорошо обобщала невидимые данные. Чтобы гарантировать справедливое и точное представление всех классов в прогнозах моделей, был использован ряд подходов для решения проблемы дисбаланса данных. Во-первых, чтобы компенсировать недостаточную представленность некоторых классов, взвешивание классов использовалось на протяжении всего этапа обучения. Классам с меньшим количеством выборок придавался больший вес. Чтобы еще больше уменьшить влияние дисбаланса классов во время обучения, также использовался ImageDataGenerator от Keras, чтобы убедиться, что модель подвергается воздействию различных классов в каждой партии30.
Метрики оценки
При оценке эффективности разработанной модели были использованы многочисленные параметры, которые описаны ниже.
Точность (ACC): Отношение точно предсказанных наблюдений ко всем наблюдениям количественно определяется этой статистикой. Когда целевые классы распределены поровну, это выгодно. Его расчет основан на уравнении 18.
(18)
Точность (PR): Точность модели измеряет, насколько хорошо она может различать все ожидаемые положительные и положительные случаи. В условиях, когда процент ложноположительных результатов существенен, это очень важно. Его расчет основан на уравнении 19.
(19)
Отзыв (R): Отзыв, также известный как чувствительность, количественно определяет процент правильно обнаруженных истинных положительных результатов. Это особенно важно в ситуациях, когда игнорирование хорошего примера может иметь серьезные последствия. Его вычисление включает в себя уравнение 20.
(20)
F1_Score: Это гармоническое среднее значение точности (PR) и полноты (R), оно обеспечивает справедливую оценку, особенно при наличии неравномерного распределения классов. Его расчет включает в себя уравнение 21.
(21)
ROC AUC: Area Under the Curve for Receiver Operating Characteristics (ROC AUC) количественно определяет способность модели различать классы. Повышенные значения AUC означают превосходную производительность модели. Его расчет основан на уравнении 22.
(22)
Каппа Коэна (CK): Эта метрика оценивает соответствие между двумя оценщиками, которые классифицируют N элементов по C взаимоисключающим классам. Он обеспечивает большую надежность по сравнению с точностью, особенно при работе с несбалансированными классами. Его расчет включает в себя уравнение 23.
(23)
Метрики Средняя абсолютная ошибка (MAE), Среднеквадратичная ошибка (RMSE) и Среднеквадратичная ошибка (MSE): Хотя эти параметры обычно используются для задач регрессии, они могут предоставить полезные сведения в сценариях классификации, особенно в порядковых классификациях. Они вычисляют разницу между реальным и прогнозируемым значениями. Эти метрики были рассчитаны с использованием уравнений 24, 25 и 26 соответственно.
(24)
(25)
(26)
Оценка F2: оценка F2 отдает приоритет запоминанию, а не точности, что ценно в сценариях, где упущение положительного прогноза влечет за собой большие затраты, чем создание ложного срабатывания. Его расчет включает в себя уравнение 27.
(27)
Кривая точности и полноты: Этот график показывает, как баланс полноты и точности на различных уровнях. Чем больше площадь под кривой, тем выше уровни как полноты, так и точности.
Вышеуказанные параметры были выбраны для того, чтобы обеспечить тщательную оценку производительности модели по всем параметрам, особенно в контексте дисбаланса набора данных и важности точного определения различных степеней остеоартрита коленного сустава.