Математические преобразования в МР-томографии

Содержание

Слайд 2

План лекции Двумерное преобразование Фурье Преобразование Фурье Теорема о свертке Фильтрация

План лекции

Двумерное преобразование Фурье
Преобразование Фурье
Теорема о свертке
Фильтрация изображений
Преобразование Радона
Преобразование Радона
Теорема о

центральном слое
Filtered back projection
Слайд 3

Одномерное преобразование Фурье Прямое преобразование сигнал во времени в спектр по

Одномерное преобразование Фурье

Прямое преобразование сигнал во времени в спектр по частоте
Обратное

преобразование переводит спектр в сигнал по времени

 

 

Слайд 4

Двумерное преобразование Фурье Аналогично одномерному случаю – прямое преобразование И обратное преобразование

Двумерное преобразование Фурье

Аналогично одномерному случаю – прямое преобразование
И обратное преобразование

 

 

Слайд 5

Двумерное дискретное преобразование Фурье Для дискретного набора данных – прямое преобразование И обратное преобразование

Двумерное дискретное преобразование Фурье

Для дискретного набора данных – прямое преобразование
И обратное

преобразование

 

 

Слайд 6

Двумерное преобразование Фурье Как и в случае одномерного преобразования Фурье, двумерное

Двумерное преобразование Фурье

Как и в случае одномерного преобразования Фурье, двумерное преобразование

является сменой базиса разложения функций
Для одномерного преобразования – одномерные гармоники, для двумерного - двумерные
Слайд 7

Изображение и его пространственный спектр Изображение Спектр Фильтр

Изображение и его пространственный спектр

Изображение Спектр Фильтр

Слайд 8

Изображение и его пространственный спектр Изображение LP HP

Изображение и его пространственный спектр

Изображение LP HP

Слайд 9

Дискретная двумерная свертка Для двух дискретных функций свертка определяется, как Пределы

Дискретная двумерная свертка

Для двух дискретных функций свертка определяется, как
Пределы суммирования могут

варьироваться в зависимости от областей определения функций
Свертка функции представляет собой точки оригинальной функции, взвешенные ядром свертки h

 

Слайд 10

Связь свертки и преобразования Фурье Теорема о свертке

Связь свертки и преобразования Фурье

Теорема о свертке

 

 

 

 

Слайд 11

Свертка, как фильтр Так как результатом свертки является модификация каждого значения

Свертка, как фильтр

Так как результатом свертки является модификация каждого значения функции

f, то операцию свертки можно использовать для создания фильтров
Например, для размытия изображения, повышения резкости, поиска краёв изображения и других.
Слайд 12

Ядро усреднения Рассмотрим действие ядра Из определения свертки – сопоставит значению

Ядро усреднения

Рассмотрим действие ядра
Из определения свертки – сопоставит значению функции значение,

усредненное с 8 соседними
В общем случае

 

 

Слайд 13

Ядро размытия Примеры действия ядра 3х3 9х9

Ядро размытия

Примеры действия ядра 3х3 9х9

Слайд 14

Ядро усреднения по Гауссу Рассмотрим действие ядра При выборе большого σ

Ядро усреднения по Гауссу

Рассмотрим действие ядра
При выборе большого σ – фильтр

размытия (усреднения)
При малом σ – не влияет на изображение

 

Слайд 15

Усреднение по Гауссу

Усреднение по Гауссу

Слайд 16

Градиент (производная первого порядка) Ядро Превитта (в зависимости от направления взятия производной) Ядро Собеля

Градиент (производная первого порядка)

Ядро Превитта (в зависимости от направления взятия производной)
Ядро

Собеля

 

 

 

 

 

 

Слайд 17

Ядро Превитта Примеры действия ядра dx dy

Ядро Превитта

Примеры действия ядра dx dy

Слайд 18

Ядро Собеля Примеры действия ядра dx dy

Ядро Собеля

Примеры действия ядра dx dy

Слайд 19

Лапласиан Сочетание двух производных второго порядка по двум координатам Форма –

Лапласиан

Сочетание двух производных второго порядка по двум координатам
Форма – из численной

аппроксимации второй производной для дискретных функций

 

Слайд 20

Лапласиан

Лапласиан

Слайд 21

Преобразование Радона Прямое преобразование: переводит двумерную функцию в её интеграл вдоль произвольной оси r θ

Преобразование Радона

Прямое преобразование: переводит двумерную функцию в её интеграл вдоль произвольной

оси

 

r

θ

Слайд 22

Преобразование Радона Пример преобразования

Преобразование Радона

Пример преобразования

Слайд 23

Обратное преобразование Радона Обратное преобразование радона (алгоритм обратной проекции) И его дискретная модель

Обратное преобразование Радона

Обратное преобразование радона (алгоритм обратной проекции)
И его дискретная модель

 

 

Слайд 24

Обратное преобразование Радона Точность реконструкции зависит от числа проекций 5 12 180

Обратное преобразование Радона

Точность реконструкции зависит от числа проекций
5 12 180

Слайд 25

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

Обратное преобразование Радона

Однако, даже при большом числе проекций реконструкция получается неточной
Для

точечного источника реконструкция имеет вид 1/r
Для практической реконструкции используется алгоритм отфильтрованой обратной проекции, основанной на теореме центрального сечения
Слайд 26

Теорема о центральном сечении Фурье-преобразование проекции функции на ось является Фурье-образом

Теорема о центральном сечении

Фурье-преобразование проекции функции на ось является Фурье-образом функции

вдоль линии, проходящей через центр координат под углом проекции

 

 

Слайд 27

Алгоритм отфильтрованной обратной проекции Предполагает реконструкцию исходного изображения из проекций, прошедших фильтрацию в частотном пространстве

Алгоритм отфильтрованной обратной проекции

Предполагает реконструкцию исходного изображения из проекций, прошедших фильтрацию

в частотном пространстве

 

 

Слайд 28

Алгоритм отфильтрованной обратной проекции Пример преобразования

Алгоритм отфильтрованной обратной проекции

Пример преобразования