Преобразование веера-луча
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.
FposArc = D*tan(FposArcDeg*pi/180);
Вычислить проекции веерной балки для геометрии «линии».
[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., Принципы компьютеризированной томографической визуализации, IEEE Press, NY, 1988, стр. 92-93.