Ультразвуковая симуляция для обучения нейросети в навигационной системе робот-ассистированного медицинского комплекса
Левин А.А., Климов Д.Д., Прохоренко Л.С., Лискевич Р.В., Мищенков Д.С., Коновалов М.Э., Тупикин К.А., Подураев Ю.В., Григорьева Е.В., Панченков Д.Н.
Введение: в статье описывается новый подход к организации навигационной системы робот-ассистированных медицинских комплексов, также находящий свое применение в методиках визуализации медицинских данных. Подход основывается на современных технологиях, таких как: применение нейронных сетей и глубокого обучения для сегментации медицинских данных, а также симуляции ультразвуковых изображений на основе диагностических данных других модальностей (КТ, МРТ). Применение данных технологий и, как следствие, комбинированного подхода на них основанного, стало возможным в настоящее время благодаря быстрому развитию видеокарт и специализированных вычислителей, что значительно снизило время, затрачиваемое на расчет данных ресурсоемких задач. Цель: апробация методики сегментации ультразвуковых изображений с помощью нейронных сетей, использующих предварительное персонализированное обучение. Материал и методы: для апробации данной методики было разработано программное обеспечение, симулирующее ультразвуковое изображение на основе данных компьютерной томографии. Полученные данные использовались для дообучения нейронной сети. Результаты: предложенный подход позволил улучшить точность сегментации в сравнении с актуальными на данный момент нейронными сетями, что положительно сказалось как на визуализации анатомических структур, так и на использовании полученных данных в навигационной системе.
Введение и методы
Ультразвуковая диагностика - популярный метод диагностики в современной медицине. Это обусловлено такими ее преимуществами, как безопасность в сравнении с методиками, которые используют рентгеновское или ионизирующее излучение; минимальная инвазивность, широкая доступность, как следствие невысокой стоимости и мобильности, а также визуализация в реальном времени. Кроме того, УЗИ можно проводить многократно и в любом возрасте, включая новорожденных и беременных пациентов. Однако данная технология обладает и рядом недостатков: ограниченная глубина проникновения, что может вызвать проблемы с визуализацией глубоко расположенных органов, особенно у пациентов с ожирением, и ограниченная визуализация вследствие того, что звуковые волны не проходят через воздух, газы в кишечнике и твердые ткани. Но самая важная проблема - это качество изображения. Ультразвуковые изображения отличаются большой зашумленностью и размытостью. В одной из работ авторы описывают показательный эксперимент: 10 специалистов проходили курс, который состоял из изучения 200 случаев ультразвуковой диагностики, при этом после каждых 50 случаев они сдавали экзаменационный тест по полученным компетенциям. При этом средние значения в определении анатомических структур выросли с 33% после 50 изученных случаев до 56,5% после 200 изученных случаев. Таким образом, интерпретация результатов во многом зависит от квалификации и опыта специалиста, проводящего исследование. В совокупности с не самым качественным изображением высока вероятность неверной интерпретации результатов и, как следствие, неверная постановка диагноза.
Вероятно, самой сложной задачей в процессе подготовки ультразвукового диагноста является умение визуально сегментировать анатомические структуры на ультразвуковом изображении. Очевидно, что данное умение приходит с опытом работы, однако современные технологии, в том числе глубокое обучение, готовы предложить различные алгоритмы сегментации изображений. Конкурс ImageNet Large Scale Visual Recognition Challenge - это ежегодное соревнование в области компьютерного зрения, которое проводится с 2010 года. Уровень точности в задаче классификации, достигнутый человеком и составляющий около 95,2%, был превзойден в 2015 году нейронной сетью ResNet от команды из Университета Торонто под руководством Джеффри Хинтона. С тех пор было представлено множество нейронных сетей, которые реализуют задачу сегментации изображений.
В области компьютерного зрения выделяют два вида сегментации: семантическая сегментация (Semantic Segmentation) и сегментация экземпляров (Instance Segmentation). Препятствием для использования нейронных сетей, обучаемых на изображениях окружающего мира, в решении задач сегментации медицинских диагностических данных является узкая специализация: такие сети изначально обучаются на изображениях, которые весьма далеки и от ультразвуковых изображений, и от изображений компьютерной и магнитно-резонансной томографии. Кроме того, изображения, получаемые с помощью средств медицинской диагностики, как правило, имеют один цветовой канал, что также усложняет работу для нейронных сетей.
Первые попытки сегментации диагностических медицинских данных были основаны на данных компьютерной и магнитно-резонансной томографии. Проект MONAI (Medical Open Network for Artificial Intelligence) достиг значимых успехов в этой области. Значимым прорывом в области сегментации именно медицинских изображений стала разработка исследователями из Фрайбургского Университета архитектуры U-Net в 2015 году. Модификации CMU-Net и CMUNeXt показывают на данный момент одни из самых лучших результатов в отношении семантической сегментации среди специализированных в данной области нейронных сетей. Метрика IoU (Intersection over Union) - это один из основных показателей, используемых для оценки качества сегментации изображений и обнаружения объектов; она измеряет степень совпадения между предсказанной и истинной областью. На специализированных ультразвуковых датасетах описанные выше сети достигают результата 60-70% точности.
Разработанная компанией Meta нейронная сеть Segment Anything Model 2 (SAM2) показывает очень хорошие результаты в отношении сегментации изображений окружающего мира, а возможность осуществления дообучения (Transfer Learning) открывает большие перспективы в плане улучшения качества ее работы в узкоспециализированных областях, в том числе и в отношении сегментации медицинских данных. На основе открытой сети SAM создано несколько открытых проектов для сегментации медицинских данных: MedSAM, MedSAM2 и UltraSAM. Оригинальная архитектура Segment Anything Model имеет три режима работы: автоматическое сегментирование и два режима работы с подсказками - в виде координат одного пикселя изображения, который принадлежит к искомому объекту, либо координат прямоугольной зоны поиска (ROI - Region of Interest). Данные подсказки решают сразу две проблемы: скорость работы и точность обнаружения различных анатомических структур сложной формы.
Принципиально все нейронные сети, описанные выше, в процессе своего дообучения используют открытые наборы медицинских данных, а потому их результаты крайне похожи по качеству работы - это еще одна проблема узкой специализации. Поскольку разметка датасета является затратной с точки зрения вычислительных ресурсов и труда специалиста, а для качественного обучения нейронной сети таких изображений должно быть в наличии как минимум несколько сотен, количество датасетов в свободном доступе конечно и весьма ограничено. Технология Data Augmentation позволяет обходить данное ограничение, однако является лишь вспомогательной и не в состоянии заменить полноценный репрезентативный набор данных.
Решением данной проблемы может являться применение синтетических данных - симулированных ультразвуковых изображений. Симуляция ультразвукового изображения - это процесс получения изображения, максимально приближенного по качеству к изображению, которое выдает аппарат ультразвуковой диагностики. Данная методика применяется как для разработки новых трансдьюсеров для ультразвуковых аппаратов, так и для обучения молодых специалистов. Современное программное обеспечение для симуляции УЗИ, такое как MUST (MatLAB UltraSound Toolbox) и Field II, симулирует ультразвуковое изображение на основе данных, полученных с помощью других средств медицинской диагностики, таких как, например, компьютерная томография, учитывая при этом такие физические характеристики, как частота работы трансдьюсера, количество элементов трансдьюсера и их физические размеры.
Результаты
Предлагаемый авторами подход заключается в создании персонализированной нейронной сети для сегментации ультразвукового изображения в реальном времени, при этом результат работы нейронной сети может использоваться как для удобства визуализации медицинских данных в режиме реального времени, так и для применения в навигационной системе робот-ассистированного комплекса. Задача персонализации решается путем дообучения нейронной сети на данных каждого конкретного пациента за счет формирования для него собственного набора обучающих данных с размеченными анатомическими структурами. Источниками таких персонализированных данных для дообучения нейронной сети могут быть: фактические данные, записанные во время предоперационной ультразвуковой диагностики и размеченные вручную квалифицированным специалистом; синтетические (симулированные) данные, полученные в результате предварительно проведенной компьютерной или магнитно-резонансной томографий.
Процесс разметки реальных данных - операция трудоемкая, однако на этапе предварительного планирования операции становится возможным определить углы обзора трансдьюсера и анатомические структуры, которые будут вызывать во время операции основной интерес, что в разы снижает объем данных, подлежащих разметке. Второй источник данных - симуляция ультразвуковых изображений на основе результатов компьютерной или магнитно-резонансных томографий. Такое предварительно проведенное исследование подлежит процессу ручной либо автоматизированной сегментации, а полученные результаты сегментации являются входными данными для симулятора ультразвуковых изображений. Данная технология позволяет получать практически неограниченные объемы размеченных данных, которые впоследствии станут обучающими для персонализированной нейронной сети. После обучения данная сеть становится готовой к использованию на реальном видеопотоке данных ультразвукового аппарата. При этом предлагаемый подход дополнительно может задействовать методику подсказок, описанную выше.
Описанная авторами ранее навигационная система комплекса робот-ассистированной хирургии имеет возможность смешивания изображений разных модальностей - наложение на реальное ультразвуковое изображение соответствующего среза данных компьютерной или магнитно-резонансной томографии. В виду того, что данные томографии были размечены ультразвуковым диагностом на этапе планирования операций, становится возможным определять границы зоны интереса, что, в свою очередь, позволяет задействовать систему подсказок для нейронной сети в виде координат пикселей определенной анатомической структуры, либо координат границ этой анатомической структуры, взятых с определенным допуском.
Результатом работы нейронной сети является набор координат для изображения-маски, которая позволяет выделить на изображении определенную структуру. Предлагаемый подход предусматривает два способа применения результатов: применение полученных изображений-масок для удобной визуализации анатомических структур на отдельных устройствах визуализации, в том числе и на устройствах дополненной реальности, и применение найденных координат контуров анатомических структур для системы навигации, в том числе для процесса отслеживания перемещений множества анатомических структур. При этом прототип системы, реализующий данный подход, имеет возможность отслеживания перемещений как с помощью конвенциональных алгоритмов технического зрения, так и с помощью алгоритмов отслеживания, использующих нейронные сети.
Эксперимент. Для проверки работоспособности способа был собран испытательный стенд. В качестве его элементов использовалось следующее оборудование: 1. УЗИ-аппарат Philips Affiniti 50; 2. Вычислительный сервер (Intel i7-13700K, ОЗУ 64 Гб, nVidia RTX 4090 24 Гб). Было вынесено предположение, что персонализированная нейронная сеть должна отличаться более высокой точностью распознавания на изображении контуров сложных анатомических структур.
Для подтверждения данного предположения был проведен эксперимент. На начальном этапе эксперимента было разработано соответствующее программное обеспечение: модуль симуляции ультразвуковых изображений, программное обеспечение для анализа массивов результатов, а также развернуты виртуальные окружения для разметки данных, дообучения и применения нейронных сетей. Для разработки программного обеспечения были использованы открытые библиотеки PyTorch, ArrayFire и OpenCV. На следующем этапе была обучена нейронная сеть. Для этого при помощи компьютерного томографа Siemens Definition AS 128 были получены DICOM-изображения многомодального фантома абдоминальной зоны человека. На основе сегментированных данных компьютерной томографии, с учетом такой информации, как глубина сканирования ультразвукового аппарата, физические характеристики используемого трансдьюсера и возможные его положения в пространстве, заранее определенные медицинским персоналом, был в автоматизированном режиме сохранен массив изображений соответствующих срезов компьютерной томографии. При этом каждое изображение из данного массива было дополнено изображением-маской, автоматически созданным на основе данных сегментации.
Полученный массив срезов компьютерной томографии был загружен в программный модуль симуляции ультразвукового изображения, что позволило получить массив симулированных ультразвуковых изображений с учетом физических параметров трансдьюсера. Массив изображений-масок компьютерной томографии был преобразован с учетом физических параметров трансдьюсера, таких как, например, угол обзора для конвексного трансдьюсера, что позволило иметь для каждого симулированного ультразвукового изображения соответствующие изображения-маски, содержащие информацию о расположении анатомических структур на симулированном изображении. Затем массив симулированных изображений вместе с соответствующими изображениями-масками был разбит на обучающий и проверочный датасеты, на которых происходило дообучение и проверка нейронной сети на базе описанного выше проекта SAM2.
На основном этапе эксперимента была осуществлена проверка работы обученной нейронной сети на записи потока фактических видеоданных ультразвукового аппарата. Кадры данной записи в количестве 50 штук были сегментированы квалифицированными диагностами, а затем проведено сравнение работы с нейронными сетями проектов UltraSAM, MedSAM2 и CMUNeXt.
Результаты эксперимента (метрика IoU, %): CMUNeXt - 64,52; UltraSAM - 71,93; MedSAM2 - 72,19; предложенный подход - 74,25. Данный эксперимент подтвердил возможность повышения качества результатов работы нейронных сетей за счет персонализации их обучения для конкретного пациента.
Обсуждение
Предложенный подход, комбинируя в себе несколько современных технологий обработки медицинских данных, позволяет решить несколько задач: использование результатов работы персонализированной нейронной сети для сегментации в системах навигации; повышение удобства визуализации анатомических структур; возможность применения алгоритмов отслеживания движения анатомических структур на видеопотоке ультразвуковых данных в режиме реального времени. Персонализация обучения нейронной сети на данных конкретного пациента, за счет комбинирования реальных и симулированных на основе КТ/МРТ данных, позволила предложенному подходу превзойти по точности сегментации (метрика IoU) как специализированную медицинскую сеть CMUNeXt, так и дообученные на медицинских данных версии общей архитектуры Segment Anything Model - UltraSAM и MedSAM2.
Заключение
Предложенный подход, комбинируя в себе несколько современных технологий обработки медицинских данных, позволяет решить несколько задач: использование результатов работы персонализированной нейронной сети для сегментации в системах навигации; повышение удобства визуализации анатомических структур; возможность применения алгоритмов отслеживания движения анатомических структур на видеопотоке ультразвуковых данных в режиме реального времени.
Иллюстрации
Список литературы
- Hertzberg B, Kliewer M, Bowie J, et al. Physician training requirements in sonography: how many cases are needed for competence? American Journal of Roentgenology. 2000; 174(5): 1221-1227. http://doi.org/10.2214/ajr.174.5.1741221
- Cardoso M, Li W, Brown R, et al. MONAI: An open-source framework for deep learning in healthcare. Springer Nature 2021 LATEX template. 2022; 1-25. http://doi.org/10.48550/arXiv.2211.02701
- Ronneberger O, Fischer P, Brox T. U-Net: Convolutional Networks for Biomedical Image Segmentation. Microscopy Today. 2024; 32(6): 13-19. http://doi.org/10.1093/mictod/qaae082
- Tang F, Wang L, Ning C, et al. CMU-Net: A Strong ConvMixer-based Medical Ultrasound Image Segmentation Network. Applied Sciences. 2022; 13(7): 4445. http://doi.org/10.3390/app13074445
- Tang F, Ding J, Wang L, et al. CMUNeXt: An Efficient Medical Image Segmentation Network based on Large Kernel and Skip Fusion. 2022. http://doi.org/10.48550/arXiv.2208.01456
- Facebook Research. Segment Anything Model 2. https://github.com/facebookresearch/sam2
- Jun M, Kim S, Li F, et al. Segment Anything in Medical Images and Videos: Benchmark and Deployment. 2024. https://doi.org/10.48550/arXiv.2408.03322
- Zhu J, Hamdi A, Qi Y, et al. Medical SAM 2: Segment medical images as video via Segment Anything Model 2. 2024. https://doi.org/10.48550/arXiv.2408.00874
- Meyer A, Murali A, Mutter D, et al. UltraSam: A Foundation Model for Ultrasound using Large Open-Access Segmentation Datasets. 2024. http://doi.org/10.48550/arXiv.2411.16222
- Teerath K, Muhammad T, Kislay R, et al. Advanced Data Augmentation Approaches: A Comprehensive Survey and Future directions. 2023. http://doi.org/10.48550/arXiv.2301.02830
- Rao M, Varghese T, Zagzebski J. Simulation of ultrasound two-dimensional array transducers using a frequency domain model. Medical physics. 2008; 35(7): 3162-3169. http://doi.org/10.1118/1.2940158
- Nicolau S, Vemuri A, Wu H, et al. A Cost Effective Simulator for Education of Ultrasound Image Interpretation and Probe Manipulation. Studies in Health Technology and Informatics. 2011; 163: 403-407. http://doi.org/10.3233/978-1-60750-706-2-403
- Damien GD. SIMUS: an open-source simulator for medical ultrasound imaging. Part I: theory & examples. Computer Methods and Programs in Biomedicine. 2022; 218(7): 106726. http://doi.org/10.1016/j.cmpb.2022.106726
- Cigier A, Varray F, Garcia D. SIMUS: an open-source simulator for ultrasound imaging. Part II: comparison with three popular simulators. 2021. http://doi.org/10.48550/arXiv.2103.05521
- Jensen J, Bk D. Transducer models in the ultrasound simulation program FIELD II and their accuracy. The Journal of the Acoustical Society of America. 2010; 127(3): 1828. http://doi.org/10.1121/1.3384240
- Field II. Fetus example. https://field-ii.dk/?examples/fetus_example/fetus_example.html
- Levin A, Klimov D, Nechunaev AA, et al. The comparison of the process of manual and robotic positioning of the electrode performing radiofrequency ablation under the control of a surgical navigation system. Scientific Reports. 2020; 10(1): 8612. https://doi.org/10.1038/s41598-020-64472-9
- Levin A, Klimov D, Nechunaev AA, et al. Assessment of experimental OpenCV tracking algorithms for ultrasound videos. Scientific Reports. 2022; 13(1). http://doi.org/10.1038/s41598-023-30930-3
- Heet T, Noopur T, Sanjana T, et al. Object Tracking by Detection using YOLO and SORT. International Journal of Scientific Research in Computer Science Engineering and Information Technology. 2020. http://doi.org/10.32628/CSEIT206256
- Konovalov M, Klimov D, Poduraev Y. Development of a Deformable Anthropomorphic Liver Phantom for Multimodal Imaging With Ultrasound and CT. IEEE Transactions on Radiation and Plasma Medical Sciences. 2023; (99): 1-1. http://doi.org/10.1109/TRPMS.2023.3303534