Обратный луч вентилятора преобразовывает
I = ifanbeam(F,D)I из данных о проекции луча вентилятора в F. Каждый столбец F содержит данные о проекции луча вентилятора при одном угле поворота. Угол между датчиками принят, чтобы быть универсальным и равным шагу между углами поворота луча вентилятора. D расстояние от вершины луча вентилятора до центра вращения.
I = ifanbeam(F,D,Name,Value)
Создайте демонстрационное изображение. phantom функция создает фантомное главное изображение.
ph = phantom(128);
Создайте преобразование луча вентилятора фантомного главного изображения.
d = 100; F = fanbeam(ph,d);
Воссоздайте фантомное главное изображение от представления луча вентилятора. Отобразите оригинальное изображение и воссозданное изображение.
I = ifanbeam(F,d); imshow(ph)
figure imshow(I);
Создайте демонстрационное изображение. Фантомная функция создает фантомное главное изображение.
ph = phantom(128);
Создайте преобразование радона изображения.
P = radon(ph);
Преобразуйте преобразование от параллельной проекции луча до проекции луча вентилятора.
[F,obeta,otheta] = para2fan(P,100,... 'FanSensorSpacing',0.5,... 'FanCoverage','minimal',... 'FanRotationIncrement',1);
Воссоздайте изображение из данных луча вентилятора.
phReconstructed = ifanbeam(F,100,... 'FanSensorSpacing',0.5,... 'Filter','Shepp-Logan',... 'OutputSize',128,... 'FanCoverage','minimal',... 'FanRotationIncrement',1);
Отобразите оригинал и преобразованное изображение.
imshow(ph)
figure imshow(phReconstructed)
F — Данные о проекции луча вентилятораДанные о проекции луча вентилятора в виде numsensors-by-numangles числовая матрица. numsensors является количеством датчиков луча вентилятора, и numangles является количеством углов поворота луча вентилятора. Каждый столбец F содержит выборки датчика луча вентилятора при одном угле поворота.
Типы данных: double | single
D — Расстояние от вершины луча вентилятора до центра вращенияРасстояние в пикселях от вершины луча вентилятора до центра вращения в виде положительного числа. ifanbeam принимает, что центр вращения является центральной точкой проекций, которая задана как ceil(size(F,1)/2). Фигура иллюстрирует D относительно вершины луча вентилятора для одной проекции луча вентилятора.
Типы данных: double | single
Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.
I = ifanbeam(F,D,'FanRotationIncrement',5)Угол поворота луча вентилятора постепенно увеличивается в градусах в виде разделенной запятой пары, состоящей из 'FanRotationIncrement' и положительная скалярная величина.
Типы данных: double
Расположение датчика луча вентилятора в виде разделенной запятой пары, состоящей из 'FanSensorGeometry' и одно из следующих значений.
Значение | Значение | Схема |
|---|---|---|
| Датчики расположены с интервалами под равными углами вдоль круговой дуги на расстоянии
|  |
| Датчики расположены с интервалами на равных расстояниях вдоль линии, которая параллельна x' ось. Самый близкий датчик является расстоянием
|  |
Filter фильтр'Ram-Lak' (значение по умолчанию) | 'Shepp-Logan' | 'Cosine' | 'Hamming' | 'Hann' | 'None'Отфильтруйте, чтобы использовать для фильтрации частотного диапазона в виде разделенной запятой пары, состоящей из 'Filter' и одно из значений в таблице. Для получения дополнительной информации смотрите iradon.
Значение | Описание |
|---|---|
| Обрезанный Поршень-Lak или фильтр пандуса. Частотная характеристика этого фильтра | |
| Умножает фильтр Поршня-Lak по |
| Умножает фильтр Поршня-Lak по |
| Умножает фильтр Поршня-Lak по Окну Хэмминга |
| Умножает фильтр Поршня-Lak по окну Hann |
'None' | Никакая фильтрация. ifanbeam возвращает неотфильтрованные данные. |
Типы данных: char | string
FrequencyScaling ScaleFactor Масштабный коэффициент для того, чтобы перемасштабировать ось частоты в виде разделенной запятой пары, состоящей из 'FrequencyScaling' и положительное число в области значений (0, 1]. Если ' меньше 1, затем фильтр сжат, чтобы поместиться в частотный диапазон FrequencyScaling'[0,FrequencyScaling], в нормированных частотах; все частоты выше FrequencyScaling установлены в 0. Для получения дополнительной информации смотрите iradon.
Типы данных: double
Interpolation — Тип интерполяции'Linear' (значение по умолчанию) | 'nearest' | 'spline' | 'pchip'Тип интерполяции используется между параллельным лучом и данными луча вентилятора в виде разделенной запятой пары, состоящей из 'Interpolation' и одно из следующих значений.
'nearest' nearestNeighbor
'linear' — Линейный (значение по умолчанию)
'spline' — Кусочный кубический сплайн
'pchip' — Кусочный кубический Эрмит (PCHIP)
Типы данных: char | string
OutputSize — Размер восстановленного изображенияРазмер восстановленного изображения в виде разделенной запятой пары, состоящей из 'OutputSize' и положительное целое число. Изображение имеет равное количество строк и столбцов.
Если вы задаете OutputSize, затем ifanbeam восстанавливает меньший или больший фрагмент изображения, но не изменяет масштабирование данных.
Примечание
Если проекции были вычислены с fanbeam функция, затем восстановленное изображение не может быть одного размера с оригинальным изображением.
Если вы не задаете OutputSize, затем размер вычисляется автоматически:
OutputSize = 2*floor(size(R,1)/(2*sqrt(2)))где R длина данных о проекции параллельного луча, используемых iradon. Для получения дополнительной информации см. Алгоритмы.
Типы данных: double
ifanbeam преобразует данные луча вентилятора, чтобы быть параллельными проекциям луча и затем использует отфильтрованный алгоритм задней проекции, чтобы выполнить обратный Радон, преобразовывают. Фильтр спроектирован непосредственно в частотном диапазоне и затем умножен на БПФ проекций. Проекции дополнены нулем к степени 2 прежде, чем отфильтровать, чтобы предотвратить пространственное доменное искажение и ускорить БПФ.
[1] Kak, A. C., и М. Слэни, принципы компьютеризированной томографической обработки изображений, Нью-Йорка, Нью-Йорка, нажатия IEEE, 1988.