🌱 Brachytherapy Isodose Distribution Optimization

Регулярне шаблонне розміщення джерел — відправна точка планування

Кожен план брахітерапії починається з фізичного обмеження: джерела (голки, катетери або зерна) можуть бути введені лише через шаблонну решітку з фіксованим кроком, найчастіше 5 мм, поєднану з трансректальним УЗ-наведенням. На цьому етапі позиції ще не оптимізовані — це рівномірна початкова сітка, яка слугує стартовою точкою для подальшого ітеративного алгоритму, а не остаточним клінічним планом.

• 5 мм: Крок шаблону (стандартна решітка Mick/Vienna)
• 10–20: Типова кількість джерел (залежно від об’єму мішені)
• ~0.5–1 Ci: Активність джерела (Ir-192) (HDR after-loading)
• 15–30 хв: Час встановлення (під УЗ/КТ-контролем)

Чому починають саме з регулярної решітки

Регулярна решітка — не довільний вибір, а наслідок клінічної практики: голки або катетери вводяться через перфорований шаблон, жорстко закріплений відносно промежини пацієнта, тому доступні позиції дискретні й наперед визначені геометрією апарата. Побудова початкового плану починається з контурування об’єму-мішені (CTV) на серії УЗ- або КТ-зрізів, після чого система пропонує рівномірний набір позицій джерел, що покривають мішень із певним запасом.

Цей рівномірний розподіл свідомо ігнорує локальну форму мішені та розташування сусідніх органів ризику — саме тому він служить лише початковим наближенням. Без подальшої оптимізації рівномірна геометрія типово дає низький індекс конформності (CI ≈ 0.5–0.6): частина периферії мішені залишається недодозованою, тоді як у центральних ділянках, де щільність джерел найвища, формуються локальні гарячі точки з дозою, що суттєво перевищує приписану. Ці недоліки й мотивують перехід до розрахунку дозового поля, а потім — до формальної інверсної оптимізації.

Клінічно прийнятний CI (>0.85–0.90) майже ніколи не досягається на етапі початкової геометрії — оптимізаційний алгоритм є обов’язковим кроком, а не опційним покращенням плану.

Суперпозиція дозових ядер — формалізм AAPM TG-43

Дозове поле кожного окремого джерела описується протоколом AAPM TG-43: доза в точці спадає приблизно обернено пропорційно квадрату відстані від джерела, з поправками на геометричну функцію, радіальну дозову функцію g(r) та анізотропний фактор F(r,θ), що враховує самопоглинання випромінювання в капсулі джерела. Сумарна доза в будь-якій точці простору — це проста лінійна суперпозиція внесків усіх активних джерел плану.

• TG-43U1: Формалізм (AAPM, 2004 (оновлений))
• табличні дані: Радіальна функція g(r) (консенсусні дані виробника)
• 0.85–1.05: Анізотропний фактор F(r,θ) (залежить від кута й капсули)
• ~1–2 мм: Розмір розрахункової сітки (типовий крок вокселя)

Від точкового джерела до тривимірного дозового поля

Базове рівняння TG-43 записується як D(r,θ) = Sk · Λ · [G(r,θ)/G(r0,θ0)] · g(r) · F(r,θ), де Sk — потужність кермa джерела, Λ — константа потужності дози, G — геометрична функція (для точкового джерела зводиться до 1/r²), g(r) — радіальна дозова функція, що описує поглинання та розсіювання в середовищі, F(r,θ) — анізотропний фактор, що враховує несиметричність випромінювання відносно осі капсули джерела.

Практично розрахунок виконується на тривимірній вокселізованій сітці: для кожного джерела обчислюється його внесок у кожну точку сітки, після чого всі внески підсумовуються — принцип лінійної суперпозиції справедливий, оскільки дозові поля не взаємодіють одне з одним у тканині. Саме ця обчислювальна структура — швидкий перерахунок сумарного поля при зміні позиції чи ваги будь-якого джерела — і робить можливою подальшу ітеративну оптимізацію, де тисячі варіантів геометрії мають оцінюватися за прийнятний час.

Для сучасних HDR-планів повний перерахунок дозового поля для сотень точок сітки й десятків джерел займає долі секунди — це і дозволяє алгоритмам оптимізації виконувати сотні ітерацій у межах секунд-хвилин.

Симульований відпал — ітеративний пошук оптимальних ваг і позицій

Інверсне планування переформульовує задачу: замість того щоб вручну підбирати позиції та ваги джерел, алгоритм мінімізує функцію вартості, що штрафує недостатнє покриття мішені, надлишкові гарячі точки та перевищення дозових лімітів органів ризику. Симульований відпал (simulated annealing) — один з найпоширеніших методів: він випадково збурює параметри плану, приймаючи як покращення, так і — з поступово спадною ймовірністю — тимчасові погіршення, що дозволяє уникнути локальних мінімумів.

• Simulated Annealing: Алгоритм (або градієнтний спуск)
• покриття − λ·ОЗІ: Функція вартості (багатокритеріальна)
• 300–1000: Типова кількість ітерацій (до збіжності плану)
• експоненційне спадання: Температурний графік (T(k) = T0 · αᵏ)

Механіка алгоритму та фізичний сенс "температури"

На кожній ітерації алгоритм пропонує невелике випадкове збурення — зсув позиції джерела або зміну його вагового коефіцієнта (часу опромінення в HDR-плануванні) — і обчислює нове значення функції вартості. Якщо вартість зменшилась, зміна приймається завжди; якщо зросла, вона приймається з ймовірністю exp(−ΔE/T), де T — поточна "температура", що поступово знижується за заданим графіком. На початку (висока температура) алгоритм охоче досліджує невигідні напрямки, уникаючи передчасного застрягання в локальному мінімумі; до кінця процесу (низька температура) пошук звужується до точного локального доопрацювання найкращого знайденого рішення.

Функція вартості типово має вигляд зваженої суми: штраф за недопокриті вокселі мішені (доза нижче 100% приписаної), штраф за гарячі точки (доза вище 150–200%) та штраф за перевищення дозових обмежень на органах ризику. Баланс між цими складовими — саме те, що на практиці регулює параметр "вага оптимізації": зсув балансу в бік сильнішого захисту ОЗІ жертвує частиною периферичного покриття мішені заради зниження ризику ускладнень.

Візуально процес відпалу проявляється як тремтіння позицій джерел, що поступово згасає, тоді як ізодозні контури синхронно стягуються навколо мішені та відступають від органа ризику — це і є збіжність алгоритму до локально-оптимального плану.

Кількісна оцінка плану — CI, DHI та об’єм гарячих точок

Після завершення оптимізації план потребує об’єктивної кількісної оцінки, незалежної від візуального враження. Індекс конформності (CI) показує, наскільки щільно приписана ізодоза облягає мішень, не захоплюючи зайвого здорового об’єму. Індекс гомогенності дози (DHI) характеризує рівномірність розподілу всередині мішені, спираючись на співвідношення об’ємів, що отримують 100% і 150% приписаної дози.

• CI = TVRI / TV: Формула CI (RTOG) (об’єм мішені в приписаній ізодозі / об’єм мішені)
• 1 − V150/V100: Формула DHI (частка мішені без гарячих точок)
• ≥0.85–0.90: Клінічно прийнятний CI (для простати LDR/HDR)
• ≥0.75: Цільовий DHI (ASTRO/ABS консенсус)

Як читати CI, DHI та V150/V100 у клінічній практиці

Індекс конформності найчастіше обчислюється як CI = V100(мішень) / V(усього об’єму, охопленого приписаною ізодозою) — досконалий план мав би CI = 1.0, тобто ізодозна поверхня точно збігалася б з межами мішені. На практиці CI 0.85–0.95 вважається клінічно прийнятним для брахітерапії простати; нижчі значення вказують або на недопокриття периферії, або на надлишкове опромінення здорової тканини навколо мішені.

Індекс гомогенності дози (DHI) розраховується через об’єми, що отримують 150% (V150) та 100% (V100) приписаної дози: DHI = 1 − V150/V100. Значення, близькі до 1, означають рівномірний дозовий розподіл без виражених гарячих точок; типова мета в сучасних HDR-протоколах — DHI ≥ 0.75, тобто об’єм гарячих точок не повинен перевищувати чверті об’єму, що отримує приписану дозу. Ці два індекси разом — CI та DHI — формують стандартну пару метрик, за якою фізик-дозиметрист приймає або відхиляє план перед клінічним затвердженням.

Високий CI без контролю DHI може приховувати небезпечні внутрішньомішеневі гарячі точки — саме тому індекси завжди оцінюються разом, а не окремо один від одного.

Клінічне затвердження — дозові ліміти органів ризику

Останній крок оптимізації брахітерапії — перевірка фінального дозового розподілу проти жорстких дозиметричних обмежень на органи ризику, специфічних для анатомічної локалізації. Для брахітерапії простати ключовими ОЗІ є пряма кишка, сечовий міхур та уретра; план затверджується лише тоді, коли одночасно виконані вимоги покриття мішені та всі дозові ліміти на суміжні структури.

• <75%: Пряма кишка D2cc (приписаної дози (типовий ліміт))
• <80%: Сечовий міхур D2cc (приписаної дози)
• <118–125%: Уретра D10 (ризик стриктури/дизурії)
• ухвалено / доопрацювання: Статус плану (фінальне рішення фізика)

Клінічне ухвалення плану: баланс покриття та безпеки

Дозові обмеження на органи ризику формулюються як об’ємні дозові параметри — наприклад D2cc (доза, що отримують найбільш опромінені 2 см³ органа), або Dmax у безпосередній близькості до критичних структур, таких як уретра. Ці ліміти встановлюються на основі клінічних даних про токсичність: перевищення дозового ліміту на пряму кишку корелює з підвищеним ризиком проктиту та ректальної кровотечі, а перевищення уретрального ліміту — зі стриктурами й тривалою дизурією.

Якщо фінальний план не задовольняє хоча б одне з обмежень, цикл оптимізації повторюється зі зміщеним балансом вагового коефіцієнта — сильніший пріоритет захисту ОЗІ ціною часткової втрати периферичного покриття мішені, або додавання нових позицій джерел у ділянках недодозування, віддалених від органа ризику. Лише план, що одночасно задовольняє CI, DHI та всі обмеження ОЗІ, передається на клінічну реалізацію.

У 2023 консенсусні рекомендації American Brachytherapy Society підкреслюють, що жоден показник конформності не може виправдати перевищення дозового ліміту на орган ризику — безпека пацієнта завжди має пріоритет над формальною оптимальністю плану.