Грубая карта груди: доза в 22 раза быстрее ✛

Оригинал: "Fast Patient-Specific Breast CT Dosimetry: 22-Fold Acceleration of Monte Carlo MGD Estimation"
arXiv:2609.03263v1 · 2026-09-03 · CC BY 4.0 · 1 мин
Монте-Карло оценка дозы для конкретной пациентки ускорилась с 27 часов до 72 минут без потери точности.

Картограф, который рисует маршрут через город, не вычерчивает каждый кирпич в стене. Ему хватает улиц, перекрёстков и расстояний. Так же и с расчётом дозы облучения при томографии груди: не нужна каждая клетка, достаточно верной общей картины.

Раньше для каждой пациентки строили цифровую модель груди из миллионов крошечных кубиков (вокселей), и компьютер перебирал их часами. Расчёт занимал около 27 часов. Это как если бы для поездки на другой конец города вы сначала изучали план каждой квартиры в каждом доме по пути.

Авторы работы проверили: насколько можно огрубить карту, чтобы маршрут остался верным?

Оказалось, что размер кубика в плоскости снимка влияет на дозу сильнее, чем его толщина. Если сделать кубики крупнее в 2,5–3 раза по ширине и высоте, доза меняется на 5–9%. А если увеличить только толщину среза в те же 2,5–3 раза, отклонение составит лишь 1–3%. Поэтому разрешение снижали неравномерно: сильнее по толщине, слабее в плоскости. Итоговый размер вокселя 0,30 × 0,30 × 3,00 мм³ удержал точность дозы в пределах 3,6–4,4% для всех четырёх образцов. Число кубиков при этом упало в 20–31 раз.

Следом сократили число моделируемых фотонов с 3 миллиардов до 200 миллионов. Статистическая погрешность осталась ниже 0,3%.

Полное время расчёта сжалось с 1600 минут до 72 минут. Это как если бы картограф вместо детального атласа взял туристическую схему и всё равно не заблудился.

Раньше было неясно, можно ли так огрублять модель груди для синхротронной фазово-контрастной КТ. Теперь ясно: можно. Метод пригодится в клиниках, где дозу нужно оценить быстро, до основного сканирования. Грубая карта оказалась не хуже подробной, зато в 22 раза быстрее.

Исходный фантом одной груди содержал до 95 миллионов вокселей, а после оптимизации их осталось около 3 миллионов. При этом разница в рассчитанной дозе составила меньше половины процента от целевого уровня облучения.

MGD = DgN \times K_{air}
MGD — средняя железистая доза, DgN — нормализованный коэффициент дозы, K_air — воздушная керма
MGD = \frac{\sum E_i}{\sum m_i}
E_i — энергия, поглощённая в железистом вокселе i, m_i — масса этого вокселя
\sigma \propto \frac{1}{\sqrt{N}}
σ — статистическая неопределённость, N — число первичных фотонных историй
Понятия
сегментация опухолей фантом компьютерная томография сегментация медицинских изображений