ifanbeam

Обратный луч вентилятора преобразовывает

Синтаксис

I = ifanbeam(F,D)

I = ifanbeam(F,D,Name,Value)

[I,H] = ifanbeam(___)

Описание

I = ifanbeam(F,D) восстанавливает изображение I из данных о проекции луча вентилятора в F. Каждый столбец F содержит данные о проекции луча вентилятора при одном угле поворота. Угол между датчиками принят, чтобы быть универсальным и равным шагу между углами поворота луча вентилятора. D расстояние от вершины луча вентилятора до центра вращения.

пример

I = ifanbeam(F,D,Name,Value) пары "имя-значение" использования, чтобы управлять различными аспектами реконструкции.

[I,H] = ifanbeam(___) также возвращает частотную характеристику фильтра, H.

Примеры

свернуть все

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

Скрипт Open Live Script

Создайте демонстрационное изображение. phantom функция создает фантомное главное изображение.

ph = phantom(128);

Создайте преобразование луча вентилятора фантомного главного изображения.

d = 100;
F = fanbeam(ph,d);

Воссоздайте фантомное главное изображение от представления луча вентилятора. Отобразите оригинальное изображение и воссозданное изображение.

I = ifanbeam(F,d);
imshow(ph)

Figure contains an axes object. The axes object contains an object of type image.

figure
imshow(I);

Figure contains an axes object. The axes object contains an object of type image.

Сгенерируйте луч вентилятора с набором Fancoverage к минимальному

Скрипт Open Live Script

Создайте демонстрационное изображение. Фантомная функция создает фантомное главное изображение.

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 contains an axes object. The axes object contains an object of type image.

figure
imshow(phReconstructed)

Figure contains an axes object. The axes object contains an object of type image.

Входные параметры

свернуть все

`F` — Данные о проекции луча вентилятора
numsensors-by-numangles числовая матрица

Данные о проекции луча вентилятора в виде numsensors-by-numangles числовая матрица. numsensors является количеством датчиков луча вентилятора, и numangles является количеством углов поворота луча вентилятора. Каждый столбец F содержит выборки датчика луча вентилятора при одном угле поворота.

Типы данных: double | single

`D` — Расстояние от вершины луча вентилятора до центра вращения
положительное число

Расстояние в пикселях от вершины луча вентилятора до центра вращения в виде положительного числа. ifanbeam принимает, что центр вращения является центральной точкой проекций, которая задана как ceil(size(F,1)/2). Фигура иллюстрирует D относительно вершины луча вентилятора для одной проекции луча вентилятора.

Типы данных: double | single

Аргументы name-value

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Пример: I = ifanbeam(F,D,'FanRotationIncrement',5)

`FanCoverage` — Область значений вращения луча вентилятора
`'cycle'` (значение по умолчанию) | `'minimal'`

Область значений вращения луча вентилятора в виде разделенной запятой пары, состоящей из 'FanCoverage' и одно из этих значений.

'cycle' — Вращайтесь через полный спектр [0, 360), степени.
'minimal' — Вращайтесь через минимальную область значений, необходимую, чтобы представлять объект.

`FanRotationIncrement` — Шаг угла поворота луча вентилятора
1 (значение по умолчанию) | положительная скалярная величина

Угол поворота луча вентилятора постепенно увеличивается в градусах в виде разделенной запятой пары, состоящей из 'FanRotationIncrement' и положительная скалярная величина.

Типы данных: double

`FanSensorGeometry` — Расположение датчика луча вентилятора
`'arc'` (значение по умолчанию) | `'line'`

Расположение датчика луча вентилятора в виде разделенной запятой пары, состоящей из 'FanSensorGeometry' и одно из следующих значений.

Значение	Значение	Схема
`'arc'`	Датчики расположены с интервалами под равными углами вдоль круговой дуги на расстоянии `D` от центра вращения. `FanSensorSpacing` задает угловой интервал в градусах.
`'line'`	Датчики расположены с интервалами на равных расстояниях вдоль линии, которая параллельна x' ось. Самый близкий датчик является расстоянием `D` от центра вращения. `FanSensorSpacing` задает расстояние между лучами вентилятора на x' ось, в пикселях.

Значение

Схема

'arc'

Датчики расположены с интервалами под равными углами вдоль круговой дуги на расстоянии D от центра вращения.

FanSensorSpacing задает угловой интервал в градусах.

'line'

Датчики расположены с интервалами на равных расстояниях вдоль линии, которая параллельна x' ось. Самый близкий датчик является расстоянием D от центра вращения.

FanSensorSpacing задает расстояние между лучами вентилятора на x' ось, в пикселях.

`FanSensorSpacing` — Интервал датчика боба вентилятора
1 (значение по умолчанию) | положительная скалярная величина

Интервал датчика боба вентилятора в виде разделенной запятой пары, состоящей из 'FanSensorSpacing' и положительная скалярная величина.

Если FanSensorGeometry 'arc', затем FanSensorSpacing задает угловой интервал в градусах.
Если FanSensorGeometry 'line', затем FanSensorSpacing задает линейное расстояние между лучами вентилятора, в пикселях. Линейный интервал измеряется на x' ось.

Типы данных: double

`Filter` фильтр
`'Ram-Lak'` (значение по умолчанию) | `'Shepp-Logan'` | `'Cosine'` | `'Hamming'` | `'Hann'` | `'None'`

Отфильтруйте, чтобы использовать для фильтрации частотного диапазона в виде разделенной запятой пары, состоящей из 'Filter' и одно из значений в таблице. Для получения дополнительной информации смотрите iradon.

Значение	Описание
`'Ram-Lak'`	Обрезанный Поршень-Lak или фильтр пандуса. Частотная характеристика этого фильтра \| `f` \|. Поскольку этот фильтр чувствителен к шуму в проекциях, один из описанных ниже фильтров может быть предпочтительным. Эти фильтры умножают фильтр Поршня-Lak по окну, которое преуменьшает роль высоких частот.
`'Shepp-Logan'`	Умножает фильтр Поршня-Lak по `sinc` функция
`'Cosine'`	Умножает фильтр Поршня-Lak по `cosine` функция
`'Hamming'`	Умножает фильтр Поршня-Lak по Окну Хэмминга
`'Hann'`	Умножает фильтр Поршня-Lak по окну Hann
`'None'`	Никакая фильтрация. `ifanbeam` возвращает неотфильтрованные данные.

Типы данных: char | string

`FrequencyScaling` ScaleFactor
1 (значение по умолчанию) | положительное число в области значений (0, 1]

Масштабный коэффициент для того, чтобы перемасштабировать ось частоты в виде разделенной запятой пары, состоящей из 'FrequencyScaling' и положительное число в области значений (0, 1]. Если 'FrequencyScaling' меньше 1, затем фильтр сжат, чтобы поместиться в частотный диапазон [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

Выходные аргументы

свернуть все

`I` — Восстановленное изображение
2D числовая матрица

Восстановленное изображение в виде 2D числовой матрицы.

`H` — Частотная характеристика
числовой вектор

Частотная характеристика фильтра, возвращенного как числовой вектор.

Типы данных: double

Советы

Чтобы выполнить обратную реконструкцию луча вентилятора, необходимо дать ifanbeam те же параметры, которые использовались, чтобы вычислить данные о проекции, F. Если вы используете fanbeam чтобы вычислить проекцию, затем убедитесь, что параметры сопоставимы при вызове ifanbeam.

Алгоритмы

ifanbeam преобразует данные луча вентилятора, чтобы быть параллельными проекциям луча и затем использует отфильтрованный алгоритм задней проекции, чтобы выполнить обратный Радон, преобразовывают. Фильтр спроектирован непосредственно в частотном диапазоне и затем умножен на БПФ проекций. Проекции дополнены нулем к степени 2 прежде, чем отфильтровать, чтобы предотвратить пространственное доменное искажение и ускорить БПФ.

Ссылки

[1] Kak, A. C., и М. Слэни, принципы компьютеризированной томографической обработки изображений, Нью-Йорка, Нью-Йорка, нажатия IEEE, 1988.

Представлено до R2006a

Документация

ifanbeam

Синтаксис

Описание

Примеры

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

Сгенерируйте луч вентилятора с набором Fancoverage к минимальному

Входные параметры

`F` — Данные о проекции луча вентилятора
numsensors-by-numangles числовая матрица

`D` — Расстояние от вершины луча вентилятора до центра вращения
положительное число

Аргументы name-value

`FanCoverage` — Область значений вращения луча вентилятора
`'cycle'` (значение по умолчанию) | `'minimal'`

`FanRotationIncrement` — Шаг угла поворота луча вентилятора
1 (значение по умолчанию) | положительная скалярная величина

`FanSensorGeometry` — Расположение датчика луча вентилятора
`'arc'` (значение по умолчанию) | `'line'`

`FanSensorSpacing` — Интервал датчика боба вентилятора
1 (значение по умолчанию) | положительная скалярная величина

`Filter` фильтр
`'Ram-Lak'` (значение по умолчанию) | `'Shepp-Logan'` | `'Cosine'` | `'Hamming'` | `'Hann'` | `'None'`

`FrequencyScaling` ScaleFactor
1 (значение по умолчанию) | положительное число в области значений (0, 1]

`Interpolation` — Тип интерполяции
`'Linear'` (значение по умолчанию) | `'nearest'` | `'spline'` | `'pchip'`

`OutputSize` — Размер восстановленного изображения
положительное целое число

Выходные аргументы

`I` — Восстановленное изображение
2D числовая матрица

`H` — Частотная характеристика
числовой вектор

Советы

Алгоритмы

Ссылки

Смотрите также

Документация Image Processing Toolbox

Поддержка

Документация

ifanbeam

Синтаксис

Описание

Примеры

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

Сгенерируйте луч вентилятора с набором Fancoverage к минимальному

Входные параметры

F — Данные о проекции луча вентилятора numsensors-by-numangles числовая матрица

D — Расстояние от вершины луча вентилятора до центра вращения положительное число

Аргументы name-value

FanCoverage — Область значений вращения луча вентилятора 'cycle' (значение по умолчанию) | 'minimal'

FanRotationIncrement — Шаг угла поворота луча вентилятора1 (значение по умолчанию) | положительная скалярная величина

FanSensorGeometry — Расположение датчика луча вентилятора 'arc' (значение по умолчанию) | 'line'

FanSensorSpacing — Интервал датчика боба вентилятора1 (значение по умолчанию) | положительная скалярная величина

Filter фильтр 'Ram-Lak' (значение по умолчанию) | 'Shepp-Logan' | 'Cosine' | 'Hamming' | 'Hann' | 'None'

FrequencyScaling ScaleFactor 1 (значение по умолчанию) | положительное число в области значений (0, 1]

Interpolation — Тип интерполяции 'Linear' (значение по умолчанию) | 'nearest' | 'spline' | 'pchip'

OutputSize — Размер восстановленного изображения положительное целое число

Выходные аргументы

I — Восстановленное изображение 2D числовая матрица

H — Частотная характеристика числовой вектор

Советы

Алгоритмы

Ссылки

Смотрите также

Документация Image Processing Toolbox

Поддержка

`F` — Данные о проекции луча вентилятора
numsensors-by-numangles числовая матрица

`D` — Расстояние от вершины луча вентилятора до центра вращения
положительное число

`FanCoverage` — Область значений вращения луча вентилятора
`'cycle'` (значение по умолчанию) | `'minimal'`

`FanRotationIncrement` — Шаг угла поворота луча вентилятора
1 (значение по умолчанию) | положительная скалярная величина

`FanSensorGeometry` — Расположение датчика луча вентилятора
`'arc'` (значение по умолчанию) | `'line'`

`FanSensorSpacing` — Интервал датчика боба вентилятора
1 (значение по умолчанию) | положительная скалярная величина

`Filter` фильтр
`'Ram-Lak'` (значение по умолчанию) | `'Shepp-Logan'` | `'Cosine'` | `'Hamming'` | `'Hann'` | `'None'`

`FrequencyScaling` ScaleFactor
1 (значение по умолчанию) | положительное число в области значений (0, 1]

`Interpolation` — Тип интерполяции
`'Linear'` (значение по умолчанию) | `'nearest'` | `'spline'` | `'pchip'`

`OutputSize` — Размер восстановленного изображения
положительное целое число

`I` — Восстановленное изображение
2D числовая матрица

`H` — Частотная характеристика
числовой вектор