Преобразование Вентилятор - Луч
F = fanbeam(I,D,Name,Value)
[ возвращает положение вентиляторно-лучевых датчиков в F,fan_sensor_positions,fan_rotation_angles]
= fanbeam(___)fan_sensor_positions и углы поворота, где проекции вентилятора-балки вычислены как fan_rotation_angles.
Установите выбор IPT, чтобы сделать оси видимыми.
iptsetpref('ImshowAxesVisible','on')
Создайте образец изображения и отобразите его.
ph = phantom(128); imshow(ph)
Вычислите проекции балки и отобразите их.
[F,Fpos,Fangles] = fanbeam(ph,250); figure imshow(F,[],'XData',Fangles,'YData',Fpos,... 'InitialMagnification','fit') axis normal xlabel('Rotation Angles (degrees)') ylabel('Sensor Positions (degrees)') colormap(gca,hot), colorbar
Вычислите проекции вентилятора-луча для геометрии 'дуги'.
I = ones(100); D = 200; dtheta = 45; [Farc,FposArcDeg,Fangles] = fanbeam(I,D,... 'FanSensorGeometry','arc',... 'FanRotationIncrement',dtheta);
Преобразуйте угловые положения в линейное расстояние вдоль оси X-prime.
FposArc = D*tan(FposArcDeg*pi/180);
Вычислите проекции вентилятора-луча для геометрии 'line'.
[Fline,FposLine] = fanbeam(I,D,... 'FanSensorGeometry','line',... 'FanRotationIncrement',dtheta);
Вычислите соответствующее преобразование Радона.
[R,Rpos]=radon(I,Fangles);
Отобразите три проекции под одним конкретным углом поворота. Обратите внимание, что три очень похожи. Различия связаны с геометрией выборки и численными приближениями, используемыми в вычислениях.
figure idx = find(Fangles==45); plot(Rpos,R(:,idx),... FposArc,Farc(:,idx),... FposLine,Fline(:,idx)) legend('Radon','Arc','Line')
I - Входное изображениеВходное изображение, заданное как 2-D числовая или логическая матрица.
D - Расстояние от вершины вентиляторной балки до центра вращенияРасстояние в пикселях от вершины балки вентилятора до центра вращения, заданное как положительное число. Центр вращения является центральным пикселем изображения, заданным как floor((size(I)+1)/2). D должен быть достаточно большим, чтобы гарантировать, что вершина веерной балки находится вне изображения при всех углах поворота. Инструкции по указанию Советов см. в разделе D. Рисунок иллюстрирует D по отношению к вершине вентилятора-балки для одной геометрии вентилятора-балки.
Задайте необязательные разделенные разделенными запятой парами Name,Value аргументы. Name - имя аргумента и Value - соответствующее значение. Name должны находиться внутри кавычек. Можно задать несколько аргументов в виде пар имен и значений в любом порядке Name1,Value1,...,NameN,ValueN.
F = fanbeam(I,D,'FanRotationIncrement',5)Шаг угла поворота вентилятора-луча в степенях, задается как разделенная разделенными запятой парами, состоящая из 'FanRotationIncrement' и положительная скалярная величина.
Типы данных: double
Позиционирование датчика-луча вентилятора, заданное как разделенная разделенными запятой парами, состоящая из 'FanSensorGeometry' и одно из следующих значений.
Значение | Значение | Диаграмма |
|---|---|---|
| Датчики расположены с равными углами вдоль дуги окружности на расстоянии
|  |
| Датчики расположены на одинаковых расстояниях вдоль линии, которая параллельна оси x '. Ближайшим датчиком является расстояние
|  |
Как руководство, попробуйте сделать D на несколько пикселей больше половины диагональной размерности изображения, вычисленного следующим образом.
sqrt(size(I,1)^2 + size(I,2)^2)
Значения, возвращенные в F являются численным приближением проекций вентилятор-балка. Алгоритм зависит от преобразования Радона, интерполированного в геометрию вентилятора-луча. Результаты варьируются в зависимости от используемых параметров. Вы можете ожидать более точных результатов, когда изображение больше, D больше, и для точек ближе к середине изображения, вдали от ребер.
[1] Kak, A.C., & Slaney, M., Principles of Computerized Tomographic Imaging, IEEE Press, NY, 1988, pp. 92-93.