SOFTWARE OF AUTOMATIC CLASSIFICATION OF A CHEST X-RAY BASED ON HYBRID CLASSIFIERS

Full Text

Стратегическая цель государственной политики в сфере здравоохранения - снижение уровня заболеваемости населения. Поэтому одним из возможных путей сохранения здоровья населения является раннее обнаружение и качественная диагностика заболеваний лёгких. В соответствии с Приказом № 124н Министерства здравоохранения РФ от 21 марта 2017 г. «Об утверждении порядка и сроков проведения профилактических медицинских осмотров граждан в целях выявления туберкулеза» [1] более половины взрослого населения нашей страны ежегодно проходят рентгенологическое обследование легких для диагностики различных заболеваний легких. Врачу-рентгенологу приходится просматривать большое количество рентгеновских снимков (РС), качество которых не всегда удовлетворительно, а размер сравнительно мал. Кроме того, яркость пикселей на изображении РС мало отличается у анатомических и патологических образований. Это затрудняет постановку диагноза врачом. Поэтому качественный анализ РС требует высокой квалификации врача-рентгенолога. При массовых обследованиях, например с помощью флюорографии, эта работа весьма утомительна и может привести к ошибкам. Широкое распространение компьютерной техники и информационных технологий создает условия для реализации автоматизированной обработки большого количества изображений. Поэтому повышение качества диагностики заболеваний легких на основе диагностических интеллектуальных систем анализа растровых изображений рентгенограмм грудной клетки и создание программного обеспечения автоматизированного рабочего места врача-рентгено-лога является важной практической инновационной задачей в области медицины. Обработка медицинских цифровых изображений на протяжении последних десятилетий является объектом пристального внимания исследователей. Множество работ посвящено математическим и программным методам в данной области, а также разработке аппаратных средств [2, 7, 8]. Интерес исследователей к современным средствам обработки рентгеновских медицинских изображений обусловлен возросшими требованиями к качеству и надежности разрабатываемых диагностических систем [3, 4, 9]. Для классификации изображений используются дескрипторы, получаемые посредством скользящих локальных окон, формируемых в координатах классифицируемых пикселей [3]. Одним из способов вычисления дескрипторов в локальным окне является нахождение двумерного спектра Фурье фрагмента изображения, попавшего в окно. Релевантность определенных таким образом дескрипторов (спектральных коэффициентов Фурье) определяют посредством нейросетевого моделирования [5, 6]. Для обработки цифровых биомедицинских изображений на кафедре биомедицинской инженерии ЮгоЗападного государственного университета (Курск) создано программное обеспечение, позволяющее получить набор дескрипторов, предназначенных для построения классификаторов изображений. Исходный интерфейс (главное окно) этого программного обеспечения, разработанного в среде Matlab 7.10, представлен на рис. 1. ІЬ^ІДифровая обработка изображений в среде MATLAB -ZWx 1 1 файл изображение методы обработки изображений помощь 1 Пневмония^вусторонняя-ТШ Рис. 1. Исходный интерфейс программного обеспечения для обработки биомедицинских изображений в среде Matlab 7.10 Цель исследования заключается в разработке инструментальной среды, предназначенной для автоматической классификации изображений рентгенограмм на основе дескрипторов Фурье. В соответствии с поставленной целью для реализации возможностей создаваемого программного обеспечения необходимо было решить следующие задачи: разработать алгоритмическое и математическое обеспечения для гибридных нейросетевых структур, выделить программные модули, реализующие спектральные методы обработки изображений. Методы На основе выдвинутой концепции построения классификатора рентгенограмм грудной клетки применялись аналитические методы исследования с целью выявления оптимальных значений качественных и количественных параметров. В качестве объекта исследования для обработки и классификации изображений были использованы рентгенограммы грудной клетки пятидесяти пациентов больницы скорой медицинской помощи г. Курска, больных пневмонией, и десять рентгенограмм с морфологическими образованиями, соответствующими онкологическим заболеваниям. В качестве референтного класса применялись рентгенограммы здоровых людей (не больных пневмонией), полученные из базы данных пациентов, имевших положительные заключения о состоянии здоровья при устройстве на работу. Эта база данных достаточно обширна, и из нее тоже было отобрано пятьдесят снимков. Контрольные и обучающие выборки формировались по рентгеновским снимкам здоровых и больных пневмонией. По рентгеновским снимкам онкологических больных формировались только контрольные выборки. Согласно этой концепции рентгенограмма грудной клетки сканируется окнами 60 Экология человека 2018.06 Медицинская экология переменного размера. При этом первоначальный размер выделенного окна полагали равным 1 % от размера исходного изображения флюорограммы грудной клетки. Из модулей отсчетов оконного двумерного преобразования Фурье определялись дескрипторы Фурье. Для каждого размера окна строится классификатор, который относят к «слабым» классификаторам. В основе построения «слабого» классификатора были использованы дескрипторы, которые определялись по гистограмме яркости амплитудного двумерного спектра Фурье [2, 10]. На рис. 2 в качестве примера представлены два фрагмента рентгенограммы грудной клетки, полученных в окнах одного и того же размера. Справа показаны соответствующие им амплитудные спектры Фурье. Рис. 2. Изображения оконных фрагментов рентгенограммы грудной клетки (слева) и их оконные двумерные преобразования Фурье (справа) Результаты При исследовании оконных спектров рентгенограмм грудной клетки были обнаружены эволюции амплитудных спектров Фурье в зависимости от наличия или отсутствия патологии. Эти эволюции могут быть установлены посредством «слабых» классификаторов, построенных на основе обучаемых нейронных сетей. В качестве информативных признаков были использованы гистограммы яркости изображения оконного амплитудного спектра Фурье. Примеры таких гистограмм показаны на рис. 3. Для сравнения, слева показаны гистограммы яркости соответствующих исходных изображений в окне. Анализ изображений, приведенных на рис. 2 и 3, показал, что форма гистограммы яркости оконных спектров имеет более высокую устойчивость по сравнению с соответствующими формами, получаемыми у изображений-оригиналов, что позволило на ее основе сформировать дескрипторы, воспользовавшись методикой, предложенной в [2, 5]. Второй «слабый» классификатор основан на дескрипторах, полученных в результате аппроксимации гистограмм в окне анализа. Для понимания этого метода построения «слабого» классификатора рассмотрим рис. 4. На рис. 4а слева показано окно анализа рентгенограммы в поле легкого без патологических изменений. Справа на этом рисунке приведена гистограмма яркости в этом окне. Анализ аналогичных гистограмм пятидесяти различных рентгеновских снимков показал, что при отсутствии патологических изменений в окне анализа гистограммы имеют многомодальную форму. На рис. 4б слева проиллюстрировано окно анализа, которое моделирует наличие заболевания пневмонией. Для моделирования в окне анализа флюорограммы здорового пациента использовалась низкочастотная фильтрация посредством двумерного преобразования Фурье. При наличии низкочастотной фильтрации исчезает модальность гистограммы яркости, и гистограмма приобретает форму, близкую к треугольной. На рис. 5 показаны окна анализа одного и того же изображения флюорограммы, у которого имеются сегменты, классифицируемые как пневмония. Рис. 3. Гистограммы яркости оконных фрагментов рентгенограммы грудной клетки (справа) и гистограммы их оконных двумерных амплитудных преобразований Фурье (слева) 61 Медицинская экология Экология человека 2018.06 Рис. 4. Окна анализа и гистограммы яркости рентгенограммы: окно без патологии - а; модель окна с патологией, полученная в результате низкочастотной фильтрации окна без патологии - б Рис. 5. Окна анализа и гистограммы яркости рентгенограммы: окно без патологии - а; окно с пневмонией - б Анализ представленных изображений позволяет сделать вывод, что и у реальных патологических образований с пневмонией гистограммы яркости в окне анализа приобретают треугольную форму. Поэтому дескрипторы второго классификатора формируются как примитивы, аппроксимирующие гистограммы яркости в окне анализа. «Слабые» классификаторы строятся для каждого размера окна. Для объединения «слабых» классификаторов в «сильные» используются обучаемые нейронные сети прямого распространения. Обучение классификаторов осуществляется для различной патологии и для каждого размера окна. Обсуждение результатов Таким образом, разработан трехуровневый классификатор и программное обеспечение для анализа и классификации рентгеновских снимков, позволяющий осуществлять анализ амплитудных спектров Фурье и гистограммы яркости пикселей изображения в скользящем окне. Рентгеновское изображение последовательно сканируется скользящими окнами различного размера. По мере продвижения окна анализа по изображению формируется амплитудный спектр Фурье и гистограмма яркости исходного изображения. На основе анализа гистограмм амплитудного спектра Фурье и гистограмм яркости в окне получают дескрипторы, которые используют в качестве пространства информативных признаков для «слабых» классификаторов, построенных на основе обучаемых нейронных сетей. В рамках данной работы классификатор построен на основе нейронных сетей прямого распространения [6, 11]. В качестве основного принципа обучения классификатора в работе реализован алгоритм обратного распространения ошибки [6, 12]. При этом следует подчеркнуть, что отличительной особенностью алгоритма настройки нейронной сети является то, что сначала нейронная сеть настраивается на тестовых образцах, моделирующих пневмонию. Затем на втором этапе настройки с реальными образцами фрагментов флюорограммы с пневмонией при необходимости осуществляется коррекция весовых коэффициентов нейронной сети. «Слабых» классификаторов, основанных на двух способах анализа, получают столько, сколько масштабов окон анализа было выбрано. На втором иерархическом уровне объединяются решения «слабых» классификаторов внутри каждого способа анализа первого иерархического уровня. Окончательное решение принимает финальный классификатор, агрегирующий решения двух классификаторов второго иерархического уровня. Алгоритмы настройки и функционирования гибридных моделей, а также 62 Экология человека 2018.06 Медицинская экология модели сегментации рентгенограмм апробированы на сегментации рентгенограмм больных пневмонией. В работе была произведена оценка качества классификации морфологических образований на изображениях рентгенограмм грудной клетки. С этой целью был сформирован критерий проверки качества классификации, основывающийся на количестве неправильно классифицированных пикселей данного œ{ класса к общему количеству пикселей этого класса на эталонном изображении: где n - количество классов, C - количество œi ’ œi правильно классифицированных пикселей œ класса, пщ - количество пикселей, в действительности /-] принадлежащих к œ€ классу. На рис. 6 представлены диаграммы, отражающие распределение этого критерия на моделях морфологических образований различных классов, полученные на 12 рентгеновских снимках. Рис. 6. Результаты критерия Mœ , полученные на 12 рентгенограммах грудной клетки: ЦЩ - модель активного контура; ЦЦ - модель сегментации на основе предварительной обработки и многооконного спектрального анализа; - модель сегментации без предварительной обработки посредством многооконного спектрального анализа На каждом из снимков формируются модели морфологических образований пневмонии (œj) и онкологии (œ2). Таким образом, сравнительный анализ показателей качества различных методов сегментации рентгеновских снимков показал, что предложенный метод многооконного спектрального преобразования позволяет осуществить дифференциальную диагно стику пневмонии и онкологических морфологических образований по критерию M ниже 15 %. Вместе с тем следует отметить ограничения автоматизированного анализа рентгеновских снимков, которые связаны с противоречием между качеством классификации морфологических образований на рентгеновских снимках и временем проведения исследований. При этом время исследования одного снимка может достигать двух-трех минут, что доставляет дискомфорт практикующему врачу. Поэтому необходимы дальнейшие исследования для преодоления этой проблемы. Результаты проведенных исследований могут быть использованы для построения интеллектуальных систем поддержки принятия решений по диагностике и прогнозированию социально значимых заболеваний. В перспективе дальнейшие исследования возможны для разработки аналогичных систем для ранней и донозологической диагностики заболеваний, а также связаны с развитием математического, алгоритмического и программного обеспечения автоматизированного рабочего места врача-рентгенолога. Международный опыт разработки подобных компьютерных систем диагностики цифровых рентгеновских снимков связан с автоматическим обнаружением областей легких на основе использования адаптивной модели легких, построенной на нечеткой логике и нейросетевом анализе [9-11].

About the authors

R A Tomakova

South-West State University

Email: rtomakova@mail.ru
Kursk

S A Filist

South-West State University

Kursk

I V Durakov

Municipal health care of emergency medical services

Kursk

References

Supplementary files