Медицинская диагностика — 08 — Конспект

$\global\def\at#1{\left. #1 \right\rvert}$ $\global\def\abs#1{\left\lvert #1 \right\rvert}$ $\global\def\norm#1{\left\lVert #1 \right\rVert}$ $\global\def\bvec#1{\mathbf{#1}}$ $\global\def\floor#1{\left\lfloor #1 \right\rfloor}$ $\global\def\limto#1{\underset{#1}{\longrightarrow}}$ $\global\def\prob#1{\mathbb{P} \left\{ #1 \right\}}$ $\global\def\mean#1{\mathbb{E} \left[ #1 \right]}$ $\global\def\disp#1{D \left[ #1 \right]}$ $\global\def\dp#1#2{#1 \cdot #2\,}$ $\global\def\vp#1#2{#1 \times #2\,}$ $\global\def\dv#1#2{\frac{d #1}{d #2}}$ $\global\def\rdv#1#2{\frac{d' #1}{d #2}}$ $\global\def\pd#1#2{\frac{\partial #1}{\partial #2}}$ $\global\def\pdv2#1#2{\frac{\partial^2 #1}{\partial #2^2}}$ $\global\def\pdvk#1#2#3{\frac{\partial^#1 #2}{\partial #3^#1}}$ $\global\def\ppdv#1#2#3{\frac{\partial^2 #1}{\partial #2 \partial #3}}$ $\global\def\pois#1{\left\{ #1 \right\}}$ $\global\def\paren#1{\left( #1 \right)}$ $\global\def\bydef#1{\overset{\mathrm{def}}{#1}}$ $\global\def\mbox#1{\text{#1}}$ $\global\def\div{\text{div}\,}$ $\global\def\dsum{\displaystyle\sum\,}$ $\global\def\grad{\text{grad}\,}$ $\global\def\rot{\text{rot}\,}$ $\global\def\vb#1{\textbf{#1}}$ $\global\def\op#1{\mathrm{#1}\,}$ $\global\def\Im{\text{Im}\,}$ $\global\def\Res{\text{Res}\,}$ $\global\def\Re{\text{Re}\,}$ $\global\def\argtg{\text{argtg}\,}$ $\global\def\ch{\text{ch}\,}$ $\global\def\const{\text{const}\,}$ $\global\def\degree{\text{degree}\,}$ $\global\def\proj{\mathrm{proj}}$ $\global\def\rank{\mathrm{rank}}$ $\global\def\res{\text{res}\,}$ $\global\def\sh{\text{sh}\,}$ $\global\def\sign{\text{sign}\,}$ $\global\def\tg{\mathrm{tg}\,}$

2024-02-13

Методы медицинской диагностики (визуализации)

Будем заниматься радионуклидными.

Главы курса:
  1. Физика визуализации
  2. Математическое моделирование
  3. Использование программного комплекса

Задачи радионуклидных методов

Реконструкция изображений по проекциям.
Восстановление функции $f(x, y)$ по проекциям $R f(\varphi, s)$: \[ R f(\varphi, s) = \int\limits_{-\infty}^{+\infty} f(s \cos \varphi - z \sin \varphi, s \sin \varphi + z \cos \varphi) dz. \]
Коррекция движения при томографических и планарных исследованиях.
Необходимо произвести коррекцию полученных в ходе радионуклидной диагностики данных на предмет движения пациента или его отдельных внутренних органов.
Моделирование транспорта индикатора
  1. Камерные линейные модели.
  2. Нелинейные камерные модели.
Трёхмерная визуализация.

Радионуклидная диагностика

Элементы ядерной физики в радионуклидной диагностики

Радиоактивный распад

Получение радионуклидов

Радиоактивный распад, сопровождающийся образованием дочернего радиоактивного продукта с иным атомным номером, позволяет простым химическим способом разделить дочерний и исходный продукт. Использование этого метода является характерной особенностью радионуклидного генератора.

Типы радиоактивного распада

  1. $\alpha$-распад;
  2. $\beta$-распад;
  3. электронный захват ($K$-захват);
  4. внутренняя конверсия (ВК);
  5. изомерные переходы (ИП);
  6. деление тяжёлых ядер;
  7. синтез лёгких ядер.
В медицине используются первые 5.

$\alpha$ и $\beta$ не подходят из-за их малого пробега в тканях, поэтому используют $\gamma$.

Выбор радиоизотопа для получения изображений

Всем этим требованиям удовлетворяет ${}^{99} Tc^m$. Он может быть получен из молибдена: \[ {}^{99} \mathrm{Mo} \limto{\beta^-} {}^{99} \mathrm{Tc}^m \limto{\text{ИП}} {}^{99} \mathrm{Tc} + \gamma. \]

Взаимодействие $\gamma$-излучения с веществом

2024-02-15

Доза излучения и её единицы

Зависимость: \[ D = f X, \qquad f = 0.873 \frac{\mu_\text{ткани}}{\mu_\text{воздуха}}, \qquad f = 33.8 \frac{\mu_\text{ткани}}{\mu_\text{воздуха}}. \] Первая $f$ — для рентгена, вторая — для Кл/кг. Мощность поглощённой дозы получается как производная: \[ D = \dv{D(t)}{t}, \] единица мощности — Гр/с. Эквивалентная доза $h$: \[ h = D k, \] где $k$ — коэффициент качества излучения.
Вид излучения Коэффициент качества ($k$)
1. рентгеновское и $\gamma$-кванты 1
2. $\beta$-частицы и электроны 1
3. $\alpha$-частицы и фотоны 10
4. тепловые нейтроны 2-5
5. быстрые нейтроны 5-10

Приборы для радионуклидной диагностики

in vivo — исследование на целом организме.
in vitro — исследование в биологических пробах.
схема 1
Установка для регистрации излучений всегда содержат 3 основных узла:
  1. один или несколько детекторов;
  2. устройства преобразования и отбора сигналов;
  3. одного или нескольких регистрирующих устройств.

Детекторы излучения. Сцинтилляционные детекторы

Типы детекторов:
ОФЭКТ и $\gamma$-томограф — одно и то же.
Современная $\gamma$-камера состоит из:
Схема 2 (функциональная блок-схема $\gamma$-камеры)
Схема 3 (блок детектирования)
Пространственное разрешение коллиматора: \[ R_c = \frac{d(L + z)}{L}, \] где Чем меньше $R_c$, тем лучше разрешение.

Геометрическая эффективность (чувствительность) коллиматора: \[ g \approx K \paren{\frac{d^2}{L(d + t)}}^2, \]

где

Если $z \gg L$, а $d \gg t$, то из формул можно получить приближённое соотношение между $g$ и $R_c$: \[ g \approx R_c^2. \] Видим, что пространственное разрешение коллиматора можно повысить только за счёт уменьшения геометрической эффективности (при постоянной толщине септы).

Алгоритмы вычисления координат

Алгоритмы можно разбить на 3 группы: Метод Энжера: используются все ФЭУ, вычисления проводятся по методу взвешенного среднего. \[ X = \frac{\sum\limits_{i=1}^{N} X_{ci} \cdot A_i}{\sum\limits_{i=1}^{N} A_i}, \qquad Y = \frac{\sum\limits_{i=1}^{N} Y_{ci} \cdot A_i}{\sum\limits_{i=1}^{N} A_i}, \] где $A_i$ (отклик сцинтилляционной вспышки) зависит от $X,Y$ (координат события). Что-то называют амплитудно-пространственными характеристиками (АПХ).

Метод максимального правдоподобия: считаем, что сигналы ФЭУ подчиняются стандартному нормальному распределению. \[ L(x,y,A_i,\sigma_i) = \prod_{i=1}^N \frac{1}{\sqrt{2\pi \sigma_i^2 (x,y)}} \exp \frac{(A_i - \overline{A_i(x,y)})^2}{2 \sigma_i^2(x,y)}. \]

Логарифмическая функция правдоподобия: \[ \begin{aligned} LL(x,y,A_i,\sigma_i) &= \ln L(x,y,A_i,\sigma_i) = \\ &= \sum\limits_{i=1}^{N} \paren{ \ln \left[ 2 \pi \sigma_i^2(x,y) \right]^{-1/2} - \frac{(A_i - \overline{A_i(x,y)})^2}{2 \sigma_i^2(x,y)} }. \end{aligned} \] Дифференцируя это выражение по $x,y$, получим систему уравнений, решение которой относительно $x,y$ и даёт оценку координат: \[ \pd{LL}{x} = 0, \qquad \pd{LL}{y} = 0. \]

Геометрические методы. Надо заранее знать АПХ ФЭУ. Полученное этими методами решение может использоваться как начальное приближение для других численных итерационных методов.

Основные метрологические характеристики сцинтилляционных $\gamma$-камер и $\gamma$-томографов

Подразделяются на: Основные характеристики.