Магнитно-резонансная томография с подавлением сигналов от нормальной ткани.


The Presentation inside:

Slide 0

Магнитно-резонансная томография с подавлением сигналов от нормальной ткани. Бабич П. В. Научные руководители: профессор, д.ф.-м.н. Ю.А. Пирогов, с.н.с., к.ф.-м.н. Н.В. Анисимов


Slide 1

Магнитно-резонансная томография с подавлением сигналов от нормальной ткани. Бабич П. В. Научные руководители: профессор, д.ф.-м.н. Ю.А. Пирогов, с.н.с., к.ф.-м.н. Н.В. Анисимов


Slide 2

МР-томограф Tomikon S50. (0.5 Тесла) Максимальный ток в катушке электромагнита - 260.64 Ампер. Охлаждение с использованием жидкого гелия (T=37.3 К). Резонансная частота для протонов – 21 МГц. Неоднородность поля < 5 м.д. в области до 40 см. Напряженность РЧ поля достигает 1.6 мкТл. Длительность 90° прямоугольного импульса – 325 мкс. Максимальное значение градиентного поля – 16 мТ/м. Время нарастания импульса – 0.5 мс. Программное обеспечение - пакет ParaVision ™ v.1.036 – надстройка к пакету XWINNMR 1.0.


Slide 3

Постановка задачи Задача исследования Разработка и реализация методик подавления сигналов от нормальных тканей в МРТ Направление исследования Задание режимов МРТ-сканирования и разработка методов преобразования данных от этих режимов для многокомпонентного подавления сигналов от нормальной ткани с целью лучшей визуализации патологически измененной ткани Основная цель исследования - хорошая визуализация патологически изменённой ткани. Основная стратегия оператора - подавление сигналов от нормальной ткани. Актуальность темы


Slide 4

Постановка задачи Задача исследования Разработка и реализация методик подавления сигналов от нормальных тканей в МРТ Направление исследования Задание режимов МРТ-сканирования и разработка методов преобразования данных от этих режимов для многокомпонентного подавления сигналов от нормальной ткани с целью лучшей визуализации патологически измененной ткани Актуальность темы Cтандартный набор средств включает : химсдвиговое подавление перенос намагниченности дифференциация по временам релаксации алгебраические преобразования комбинированные методы


Slide 5

Постановка задачи Задача исследования Разработка и реализация методик подавления сигналов от нормальных тканей в МРТ Направление исследования Задание режимов МРТ-сканирования и разработка методов преобразования данных от этих режимов для многокомпонентного подавления сигналов от нормальной ткани с целью лучшей визуализации патологически измененной ткани Актуальность темы Актуальность двух последних направлений следует из их основных преимуществ: Сокращение времени стандартного исследования. Разгрузка производительных мощностей томографа.


Slide 6

Постановка задачи Задача исследования Разработка и реализация методик подавления сигналов от нормальных тканей в МРТ Направление исследования Задание режимов МРТ-сканирования и разработка методов преобразования данных от этих режимов для многокомпонентного подавления сигналов от нормальной ткани с целью лучшей визуализации патологически измененной ткани Цели и задачи Произвести математический анализ А) МРТ-сигналов для стандартных режимов сканирования, включая режимы с подавлением одного или двух типов нормальных тканей. Б) МРТ-изображений, получаемых в результате алгебраических преобразований данных, получаемых от стандартных режимов сканирования.


Slide 7

Постановка задачи Задача исследования Разработка и реализация методик подавления сигналов от нормальных тканей в МРТ Направление исследования Задание режимов МРТ-сканирования и разработка методов преобразования данных от этих режимов для многокомпонентного подавления сигналов от нормальной ткани с целью лучшей визуализации патологически измененной ткани Цели и задачи Произвести математический анализ методики Провести экспериментальную проверку А) получить данные от сканирования с подавлением нормальных тканей. Б) разработать математический аппарат для алгебраических операций с данными от разных режимов МРТ-сканирования. В) провести сопоставление получаемых результатов для нескольких типов нормальной и патологически измененной тканей.


Slide 8

Постановка задачи Задача исследования Разработка и реализация методик подавления сигналов от нормальных тканей в МРТ Направление исследования Задание режимов МРТ-сканирования и разработка методов преобразования данных от этих режимов для многокомпонентного подавления сигналов от нормальной ткани с целью лучшей визуализации патологически измененной ткани Цели и задачи Произвести математический анализ методики Провести экспериментальную проверку Провести обработку данных МРТ-исследований А) разработать и реализовать алгоритм сегментации зоны повышенного МР-сигнала (зоны поражения) и методику оценки площади сегментируемой зоны. Б) произвести оценку изменения объёма злокачественного образования за продолжительный промежуток времени.


Slide 9

Постановка задачи Задача исследования Разработка и реализация методик подавления сигналов от нормальных тканей в МРТ Направление исследования Задание режимов МРТ-сканирования и разработка методов преобразования данных от этих режимов для многокомпонентного подавления сигналов от нормальной ткани с целью лучшей визуализации патологически измененной ткани Цели и задачи Произвести математический анализ методики Провести экспериментальную проверку Провести обработку данных МРТ-исследований Дополнительно Цели и задачи предэкспериментального этапа


Slide 10

Постановка задачи Дополнительно Цели и задачи предэкспериментального этапа Провести обработку данных МРТ-исследований Разработать пакет программ прикладного характера с открытым кодом и возможностью модификации и модернизации, для выполнения следующих функции Конвертация Конвертация данные исходных сигналов (К-пространство), матриц 3D изображений в матрицы для работы в среде Matlab, для визуализации и в пользовательские форматы изображений (JPEG, TIFF,PNG) как в прямом направлении так и в обратном, по необходимости вышепоставленных задач. К-пространство


Slide 11

Постановка задачи Цели и задачи предэкспериментального этапа Разработать пакет программ прикладного характера с открытым кодом и возможностью модификации и модернизации выполняющего следующие функции Конвертация Математический редактор Математический редактор и просмоторщик К-пространства (эксперименты по определению информативности К-пространства и математической корректировки К-пространства для создания алгоритмов устранения искажений изображений).


Slide 12

Постановка задачи Цели и задачи предэкспериментального этапа Разработать пакет программ прикладного характера с открытым кодом и возможностью модификации и модернизации выполняющего следующие функции Конвертация Математический редактор Модуль графического редактора Модуль для использования графического редактора Microsoft Paint в качестве редактора слоёв трёхмерного МРТ-изображения (для устранения искажений из изображения, связанных с аппаратурными несовершенствами).


Slide 13

Постановка задачи Цели и задачи предэкспериментального этапа Разработать пакет программ прикладного характера с открытым кодом и возможностью модификации и модернизации выполняющего следующие функции Сегментация Конвертация Математический редактор Модуль графического редактора Сегментация трёхмерных областей и определения их объёма.


Slide 14

Постановка задачи Цели и задачи предэкспериментального этапа Разработать пакет программ прикладного характера с открытым кодом и возможностью модификации и модернизации выполняющего следующие функции Сегментация Конвертация Математический редактор Модуль графического редактора Математические операции Математические операции над матрицами МРТ-изображений.


Slide 15

Постановка задачи Цели и задачи предэкспериментального этапа Разработать пакет программ прикладного характера с открытым кодом и возможностью модификации и модернизации выполняющего следующие функции Сегментация Конвертация Математический редактор Модуль графического редактора Математические операции Другие приложения Для достижения некоторых целей дипломной работы дополнительно разрабатывалось прикладное программное обеспечение.


Slide 16

Постановка задачи Цели и задачи предэкспериментального этапа Разработать пакет программ прикладного характера с открытым кодом и возможностью модификации и модернизации выполняющего следующие функции Сегментация Конвертация Математический редактор Модуль графического редактора Математические операции Другие приложения Параллельно с решением вышеперечисленных задач разрабатывалось прикладное программное обеспечение для коррекции аппаратурных артефактов на МРТ-изображениях


Slide 17

Постановка задачи Цели и задачи предэкспериментального этапа Разработать пакет программ прикладного характера с открытым кодом и возможностью модификации и модернизации выполняющего следующие функции Сегментация Конвертация Математический редактор Модуль графического редактора Математические операции Другие приложения Параллельно с решением вышеперечисленных задач разрабатывалось прикладное программное обеспечение для коррекции аппаратурных артефактов на МРТ-изображениях Данное прикладное программное обеспечения должно конвертировать МРТ-данные из специализированного фирменного формата, разработанного для ОС IRIX на платформе SiliconGraphics, в форматы, пригодные для графических редакторов в среде MS Windows на платформе PC. Благодаря этому предполагалось решать задачи: А) анализа и обработки исходных МРТ-сигналов Б) коррекции аппаратурных артефактов на МРТ-изображениях В) обратная конвертация данных для работы в среде IRIX


Slide 18

T2 (~150 сек) Практическая часть: математический анализ: МРТ-сигналов Т2 (T2 – взвешенного изображения) График зависимости максимально возможной интенсивности регистрации сигнала вещества от параметра T1: Mz(T1) (параметр спин-решёточной релаксации)


Slide 19

FLAIR(~210 сек) Практическая часть: математический анализ: МРТ-сигналов FLAIR (T2 – взвешенного изображения с подавлением жировой ткани) График зависимости максимально возможной интенсивности регистрации сигнала вещества от параметра T1: Mz(T1) (параметр спин-решёточной релаксации)


Slide 20

STIR (~210 сек) Практическая часть: математический анализ: МРТ-сигналов STIR (T2 – взвешенного изображения с подавлением свободной воды) График зависимости максимально возможной интенсивности регистрации сигнала вещества от параметра T1: Mz(T1) (параметр спин-решёточной релаксации)


Slide 21

DIR (~330 сек) Практическая часть: математический анализ: МРТ-сигналов DIR (одновременное подавление жировой ткани и воды) График зависимости максимально возможной интенсивности регистрации сигнала вещества от параметра T1: Mz(T1) (параметр спин-решёточной релаксации)


Slide 22

Графики рассчитаны по приближённым формулам: при значениях пар-в: TR/TIw /TIf = 5600 / 1300 / 80 ms . Практическая часть: математический анализ: МРТ-сигналов Основу методики анализа способов эмуляции составляют данные об операциях над режимами сканирования:


Slide 23

Практическая часть: математический анализ: МРТ-сигналов


Slide 24

Практическая часть: математический анализ: МРТ-сигналов


Slide 25

Практическая часть: математический анализ: МРТ-сигналов


Slide 26

Практическая часть: математический анализ: МРТ-сигналов


Slide 27

Практическая часть: математический анализ: МРТ-сигналов Более качественный анализ включает в себя: 1) Анализ информативности сигнала отрицательных значений пикселей после операции на примере экспериментальных данных


Slide 28

Практическая часть: математический анализ: МРТ-сигналов Более качественный анализ включает в себя: 1) Анализ информативности сигнала отрицательных значений пикселей после операции на примере экспериментальных данных. 2) анализ контрастности изображения.


Slide 29

Практическая часть: математический анализ: МРТ-сигналов Более качественный анализ включает в себя: 1) Анализ информативности сигнала отрицательных значений пикселей после операции на примере экспериментальных данных. 2) анализ контрастности изображения. 3) оценка возможности и целесообразности применения данной операции в медицинской диагностике. 4) сравнение результатов с результатами стандартного исследования.


Slide 30

Методика сегментации области


Slide 31

Увеличение точности определения координаты. производится путём повторного указания координаты точки на изображении окрестности точки, которая была зафиксирована при первом нажатии.


Slide 32

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


Slide 33

Фильтрация Для увеличения стабильности работы методики был дополнительно разработан алгоритм фильтрации сегментируемого изображения. Фильтр-матрица было подобрана с учётом следующих критерий: Стабильность работы основной программы. Быстродействие алгоритма. Сходимость алгоритма сегментации с данным фильтром.


Slide 34

Фильтр


Slide 35

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0


Slide 36


Slide 37

Критерий детектирования


Slide 38


Slide 39

Определение границы указанной области Для объяснения принципа работы алгоритма рассмотрим граф, имеющий следующие свойства: Все вершины строятся только на ячейках полученной матрицы. Все рёбра одинаковой длины, ориентированны, и могут принимать только одну из четырёх ориентаций (право, лево, верх, низ). Все вершины имеют нумерацию (1,2,3,4, и т.д.) Каждая вершина имеет одно входящее (из вершины с меньшим номером) и одно выходящее ребро (в вершину с большим номером).


Slide 40


Slide 41

Граница сегментации выделена и может быть скорректирована Есть возможность настройки яркости и контрастности Есть возможность настройки яркости и контрастности Исследование протоколируется Проверяется форма и размер сегментируемой области Проверяется форма и размер сегментируемой области


Slide 42


Slide 43

ОПИСАНИЕ ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ Эмуляция режима DIR


Slide 44

ОПИСАНИЕ ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ Эмуляция режима STIR


Slide 45

Эмуляция режима двухкомпонентного подавления сигналов (DIR) используя данные режимов STIR и Т2. Используя ранее введённые схемы: А) T2 - FLAIR = STIR Б) STIR x FLAIR = DIR Получаем: (T2 - FLAIR) x FLAIR = STIR x FLAIR = DIR Это означает, что существует методика эмуляции режима DIR из режимов FLAIR и Т2:


Slide 46

T2 Практическая часть: математический анализ МРТ-изображений FLAIR (под. воды) DIR (жира и воды) STIR (под. жира) Параметры сканирования: TSE: TR/TEeff=5600/100 ms, ETL=8, NS=22 (6mm) FOV=160x180 mm, Matrix=160(antialise) x180, (in plane res.=1x1 mm) Особенности сканирования: Необходимое условие – одинаковые параметры сканирования перечисленные выше.


Slide 47

T2 Практическая часть: математический анализ МРТ-изображений FLAIR (под. воды) DIR (жира и воды) STIR (под. жира) Два варианта: Операция вычитания Стандартное сканирование


Slide 48

Практическая часть: математический анализ МРТ-изображений


Slide 49

T2 Практическая часть: математический анализ МРТ-изображений FLAIR (под. воды) DIR (жира и воды) STIR (под. жира) Два варианта: Операция вычитания Стандартное сканирование


Slide 50

T2 Практическая часть: математический анализ МРТ-изображений FLAIR (под. воды) DIR (жира и воды) STIR (под. жира) Два варианта: Стандартное сканирование Операция умножения


Slide 51

Практическая часть: математический анализ МРТ-изображений


Slide 52

T2 Практическая часть: математический анализ МРТ-изображений FLAIR (под. воды) DIR (жира и воды) STIR (под. жира) Два варианта: Стандартное сканирование Операция умножения


Slide 53

T2 Практическая часть: математический анализ МРТ-изображений FLAIR (под. воды) DIR (жира и воды) STIR (под. жира) Эмуляция режимов: d=bx(a-b) c=a-b d=bxc


Slide 54

Практическая часть: сопоставление получаемых результатов Т2 FLAIR (TI=1300 ms)


Slide 55

Практическая часть: сопоставление получаемых результатов «Т2» - «FLAIR» (Эмуляция STIR) STIR (Истинный STIR)


Slide 56

Практическая часть: сегментация зоны повышенного МР-сигнала


Slide 57


Slide 58

Практическая часть: сегментация зоны повышенного МР-сигнала


Slide 59

Выводы Результаты проделанной работы сводятся к следующему: 1) Разработан новый метод обработки данных МРТ путём алгебраических операций с МРТ-изображениями, полученных при разных режимах сканирования. Метод сводится к получению новых изображений, представляющих собой эмуляцию особых режимов сканирования, включая режимы с подавлением сигналов от нормальной ткани. 2) Для теоретического обоснования метода проведен расчёт величин сигналов для исходных и эмулированных режимов сканирования.


Slide 60

Выводы Результаты проделанной работы сводятся к следующему: 4) Проведено сопоставление расчетных и экспериментально данных, благодаря чему выявлена корреляция между изменениями контраста на изображениях, полученных при разных режимах сканирования, проведенных для одной и той же зоны сканирования и одинаковой схемы регистрации МРТ-сигналов. Показана возможность использования эмулированных изображений в диагностических исследованиях. 3) Для экспериментального обоснования метода проведены алгебраические операции с МРТ-изображениями и дана оценка тканевого контраста на исходных и эмулированных МРТ-изображениях.


Slide 61

Выводы Результаты проделанной работы сводятся к следующему: 5) Разработан оригинальный алгоритм автоматической сегментации для обработки МРТ-изображений, полученных с подавлением сигналов от нормальных тканей, который обладает широкими возможностями для подбора параметров сегментации, что позволяет его использовать в условиях слабого контраста и малого соотношения сигнал/шум. 6) Проведён расчёт динамики роста зоны поражения при астроцитоме за пятилетний период. Для расчета использованы изображения, полученные с подавлением сигналов от нормальной ткани. Для объемной обработки данных применён оригинальный алгоритм сегментации и разработанная в рамках задач данной дипломной работы методика коррекции артефактов.


Slide 62

7) Разработан пакет программ для обработки МРТ-данных, который позволяет проводить следующие операции: а) сегментацию и измерение объёма исследуемой структуры; б) математические операции над МР-изображениями; в) проекционные преобразования с матрицами данных; г) коррекцию артефактов на МР-изображении; д) визуализацию и обработку (включая Фурье-преобразование) исходных МРТ-данных (К-пространства); е) пакетные операции по интерполяции МРТ данных и конвертации изображений из специализированного в форматы, приемлемые для широкого круга пользователей персональных компьютеров. С помощью данного пакета программ были решены задачи, поставленные в дипломной работе. Программы могут быть адаптированы для решения других задач, связанных с обработкой МРТ-данных. Выводы Результаты проделанной работы сводятся к следующему:


Slide 63


Slide 64


×

HTML:





Ссылка: