Преобразуйте проекции вентилятора-балки в параллельные
P = fan2para(F,D,Name,Value)
[ возвращает местоположение датчика параллельного луча в P,parallel_sensor_positions,parallel_rotation_angles] = fan2para(___)parallel_sensor_positions и углы поворота в parallel_rotation_angles.
Создайте синтетические данные параллельного луча.
ph = phantom(128);
Вычислите параллельное преобразование луча и отобразите его.
theta = 0:179; [Psynthetic,xp] = radon(ph,theta); imshow(Psynthetic,[],... 'XData',theta,'YData',xp,'InitialMagnification','fit') axis normal title('Synthetic Parallel-Beam Data') xlabel('\theta (degrees)') ylabel('x''') colormap(gca,hot), colorbar
Преобразуйте данные параллельного луча в диаграммы направленности вентилятора.
Fsynthetic = para2fan(Psynthetic,100,'FanSensorSpacing',1);Восстановите исходные данные параллельного луча.
[Precovered,Ploc,Pangles] = fan2para(Fsynthetic,100,... 'FanSensorSpacing',1,... 'ParallelSensorSpacing',1); figure imshow(Precovered,[],... 'XData',Pangles,'YData',Ploc,'InitialMagnification','fit') axis normal title('Recovered Parallel-Beam Data') xlabel('Rotation Angles (degrees)') ylabel('Parallel Sensor Locations (pixels)') colormap(gca,hot), colorbar
F - Данные проекции вентилятора-лучаДанные проекции вентилятора-луча, заданные как числовая матрица. Каждый столбец F содержит данные о пучке вентиляторов под одним углом поворота. Количество столбцов указывает количество углов поворота вентилятора-луча, а количество строк указывает количество датчиков вентилятора-луча.
Типы данных: double | single
D - Расстояние от вершины вентиляторной балки до центра вращенияРасстояние в пикселях от вершины балки вентилятора до центра вращения, заданное как положительное число. fan2para принимает, что центр поворота является центральной точкой проекций, которая определяется как ceil(size(F,1)/2). Рисунок иллюстрирует D относительно вершины вентилятора-балки для одной проекции вентилятора-балки.
Типы данных: double | single
Задайте необязательные разделенные разделенными запятой парами Name,Value аргументы. Name - имя аргумента и Value - соответствующее значение. Name должны находиться внутри кавычек. Можно задать несколько аргументов в виде пар имен и значений в любом порядке Name1,Value1,...,NameN,ValueN.
P = fan2para(F,D,'FanRotationIncrement',5)Шаг угла поворота вентилятора-луча в степенях, задается как разделенная разделенными запятой парами, состоящая из 'FanRotationIncrement' и положительная скалярная величина.
Типы данных: double
Позиционирование датчика-луча вентилятора, заданное как разделенная разделенными запятой парами, состоящая из 'FanSensorGeometry' и одно из следующих значений.
Значение | Значение | Диаграмма |
|---|---|---|
| Датчики расположены с равными углами вдоль дуги окружности на расстоянии
|  |
| Датчики расположены на одинаковых расстояниях вдоль линии, которая параллельна оси x '. Ближайшим датчиком является расстояние
|  |
'Interpolation' - Тип интерполяции'Linear' (по умолчанию) | 'nearest' | 'spline' | 'pchip'Тип интерполяции, используемой между данными параллельного луча и данных о вентиляторе, задается как разделенная разделенными запятой парами, состоящая из 'Interpolation' и одно из этих значений.
'nearest' - Ближайшая соседка
'linear' - Линейный (значение по умолчанию)
'spline' - Кусочно-кубический сплайн
'pchip' - Кусочно-кубический гермит (PCHIP)
'ParallelCoverage' - Область значений поворота параллельного луча'halfcycle' (по умолчанию) | 'cycleОбласть значений поворота параллельного луча, заданный как разделенная разделенными запятой парами, состоящая из 'ParallelCoverage' и одно из этих значений.
'cycle' - Параллельные данные охватывают полную область значений [0, 360) степеней.
'halfcycle' - Параллельные данные покрывают [0, 180) степеней.
'ParallelRotationIncrement' - Шаг угла поворота параллельного лучаШаг угла поворота параллельного луча в степенях, заданный как разделенная разделенными запятой парами, состоящая из 'ParallelRotationIncrement' и положительная скалярная величина k таким образом, что 180/k является целым числом. Если вы не задаете ParallelRotationIncrement, тогда значение по умолчанию равно FanRotationIncrement.
Типы данных: double
'ParallelSensorSpacing' - Расстояние между датчиками параллельного лучаДатчик параллельного луча интервала в пикселях, задается как разделенная запятой пара, состоящая из 'ParallelSensorSpacing' и положительная скалярная величина. Область значений расположения датчика параллельного луча вычисляется из области значений углов наклона вентилятора, fanangles, согласно: [.D* sin (min (fanangles)) D * sin (max (fanangles))]
Если вы не задаете ParallelSensorSpacingзатем интервал принимается равномерным и устанавливается на минимальное расстояние, подразумеваемое углами вентилятора и дискретизируемое в области значений, подразумеваемом углами вентилятора.
Типы данных: double