Преобразование проекций параллельной балки в вентиляторную балку
F = para2fan(P,D)P к данным вентиляторной балки F. Каждый столбец P содержит образцы датчика параллельного луча под одним углом поворота. D - расстояние от вершины веерной балки до центра вращения.
Предполагается, что датчики с параллельным лучом имеют интервал в один пиксел. Углы поворота параллельного луча расположены на равном расстоянии друг от друга для покрытия [0,180] градусов. Вычисленные углы поворота веерной балки имеют такое же расстояние, как углы поворота параллельной балки, и перекрывают [0,360) градусов. Вычисленные углы пучка вентиляторов разнесены на равные расстояния с интервалом, установленным на наименьший угол, подразумеваемый интервалом между датчиками.
F = fan2para(P,D,Name,Value)
[ возвращает расположение датчика луча вентилятора в F,fan_sensor_positions,fan_rotation_angles] = fan2para(___)fan_sensor_positions и углы поворота в fan_rotation_angles.
Создание проекций с параллельными балками
ph = phantom(128); theta = 0:180; [P,xp] = radon(ph,theta); imshow(P,[],'XData',theta,'YData',xp,'InitialMagnification','fit') axis normal title('Parallel-Beam Projections') xlabel('\theta (degrees)') ylabel('x''') colormap(gca,hot), colorbar
Преобразовать в вентиляторно-балочные проекции
[F,Fpos,Fangles] = para2fan(P,100); figure imshow(F,[],'XData',Fangles,'YData',Fpos,'InitialMagnification','fit') axis normal title('Fan-Beam Projections') xlabel('\theta (degrees)') ylabel('Sensor Locations (degrees)') colormap(gca,hot), colorbar
P - Данные проекции параллельного лучаДанные проекции параллельного луча, заданные как числовая матрица. Каждый столбец P содержит данные параллельного луча под одним углом поворота. Число столбцов указывает количество углов поворота параллельного луча, а число строк - количество датчиков параллельного луча.
Типы данных: double | single
D - Расстояние от вершины балки вентилятора до центра вращенияРасстояние в пикселях от вершины луча вентилятора до центра вращения, заданное как положительное число. para2fan предполагает, что центр вращения является центральной точкой проекций, которая определяется как ceil(size(F,1)/2). Значение D должно быть больше или равно ParallelSensorSpacing*(SIZE(P,1)-1)/2
Рисунок иллюстрирует D относительно вершины веер-балки для одной проекции веер-балки.
Типы данных: double | single
Укажите дополнительные пары, разделенные запятыми Name,Value аргументы. Name является именем аргумента и Value - соответствующее значение. Name должен отображаться внутри кавычек. Можно указать несколько аргументов пары имен и значений в любом порядке как Name1,Value1,...,NameN,ValueN.
F = para2fan(P,D,'FanRotationIncrement',5)'FanRotationIncrement' - Приращение угла поворота балки вентилятораПриращение угла поворота балки вентилятора в градусах, указанное как пара, разделенная запятыми, состоящая из 'FanRotationIncrement' и положительный скаляр.
Если FanCoverage является 'cycle', то 360/FanRotationIncrement должно быть целым числом.
Если не указать FanRotationIncrement, то значение по умолчанию равно углу поворота параллельного луча.
Типы данных: double
Расположение датчика луча вентилятора, указанное как пара, разделенная запятыми, состоящая из 'FanSensorGeometry' и одно из следующих значений.
Стоимость | Значение | Диаграмма |
|---|---|---|
| Датчики разнесены под равными углами по дуге окружности на расстояние
|  |
| Датчики разнесены на равные расстояния вдоль линии, параллельной оси x '. Ближайший датчик - расстояние
|  |
'FanSensorSpacing' - Расстояние между датчиками вентиляторовРасстояние между датчиками вентиляторов, указанное как пара, разделенная запятыми, состоящая из 'FanSensorSpacing' и положительный скаляр.
Если FanSensorGeometry является 'arc', то FanSensorSpacing определяет угловой интервал в градусах.
Если FanSensorGeometry является 'line', то FanSensorSpacing определяет линейное расстояние между веерными лучами в пикселях. Линейный интервал измеряется на оси x '.
Если не указать FanSensorGeometry, то значение по умолчанию FanSensorSpacing наименьшее значение, подразумеваемое ParallelSensorSpacing чтобы:
Если FanSensorGeometry является 'arc', то FanSensorSpacing является 180/pi*ASIN(ParallelSensorSpacing/D)
Если FanSensorGeometry является 'line', то FanSensorSpacing является D*ASIN(ParallelSensorSpacing/D)
Типы данных: double
'Interpolation' - Тип интерполяции'Linear' (по умолчанию) | 'nearest' | 'spline' | 'pchip'Тип интерполяции, используемой между данными параллельного луча и веерного луча, определяемый как разделенная запятыми пара, состоящая из 'Interpolation' и одно из этих значений.
'nearest' - Ближайший сосед
'linear' - Линейный (по умолчанию)
'spline' - Кусочно-кубический сплайн
'pchip' - Кусочно-кубический эрмит (PCHIP)
'ParallelCoverage' - Диапазон вращения параллельной балки'halfcycle' (по умолчанию) | 'cycleДиапазон вращения параллельного луча, определяемый как разделенная запятыми пара, состоящая из 'ParallelCoverage' и одно из этих значений.
'cycle' - Параллельные данные охватывают весь диапазон [0, 360) градусов.
'halfcycle' - Параллельные данные охватывают [0, 180) градусов.
'ParallelSensorSpacing' - Расстояние между параллельно-балочными датчиками1 | положительный скалярРасстояние между датчиками параллельного луча в пикселях, указанное как пара, разделенная запятыми, состоящая из 'ParallelSensorSpacing' и положительный скаляр.
Типы данных: double