Параллельные лучи - реконструкция фантома головы на основании проекционных данных.

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

Следующие три примера реконструкции (I1, I2 и I3) демонстрируют эффект изменения числа углов, по которым сделаны проекции. Для I1 и I2 некоторые особенности, которые четко заметны на исходном фантоме, видны не очень качественно. Что касается, элипсообразных предметов, то они достаточно хорошо заметны на каждом изображении. Результирующее изображение I3 наиболее точно соответствует исходному изображению.

Как уже отмечалось выше, на изображениях I1 и I2 есть некоторые артефакты, которые искажают действительную информацию. Для устранения этих артефактов используют большее число углов.

% Построение результирующей размерности

% при реконструкции с использованием

% размеров исходного изображения P.

output_size=max(size(P));

dtheta1=theta1(2)-theta1(1);

I1=iradon(R1, dtheta1, output_size);

figure, imshow(I1)

dtheta2=theta2(2)-theta2(1);

I2=iradon(R2, dtheta2, output_size);

figure, imshow(I2)

dtheta3=theta3(2)-theta3(1);

I3=iradon(R3, dtheta3, output_size);

figure, imshow(I3)

Веерные лучи - вычисление синтезированных проекций.

Рассмотрим вычисление синтезированных проекций с использованием веерно-лучевой геометрии и переменной 'FanSensorSpacing'.

D=250;

dsensor1=2;

F1=fanbeam(P, D, 'FanSensorSpacing', dsensor1);

dsensor2=1;

F2=fanbeam(P, D, 'FanSensorSpacing', dsensor2);

dsensor3=0.25;

[F3, sensor_pos3, fan_rot_angles3]=fanbeam(P, D, ...

'FanSensorSpacing', dsensor3);

Отобразим проекционные данные F3. Отметим, что веерный угловой поворот в диапазоне от 0 до 360 градусов может быть

преобразован в модель с использованием 180 градусов, поскольку некоторые свойства будут идентичны с обеих сторон.

Некоторые свойства на изображении при веерно-лучевых проекциях могут быть связаны с некоторыми свойствами на

изображениях параллельно-лучевых проекций и наоборот.

figure, imagesc(fan_rot_angles3, sensor_pos3, F3)

colormap(hot)

colorbar

Xlabel('Веерный поворот угла (в градусах) ')

ylabel('Веерное расположение датчиков (в градусах)')

Веерные лучи - реконструкция фантома головы на основании проекционных данных.

Каждому веерному расположению сенсоров в пространстве соответствует своя синтезированная проекция. При практической реализации реконструкции всегда известна геометрия источников и сенсоров изображения P. Изменение значения 'FanSensorSpacing' приводит к изменению числа сенсоров, которые используются при изменении угла. При каждой веерно-лучевой реконструкции используется поворот на все значения угла. Использование данных при различных поворотах угла влияет на качество результата обработки.

Отметим, что 'FanSensorSpacing' является только одним параметром из нескольких, которые можно контролировать при использовании функций fanbeam и ifanbeam. Также существует возможность конвертации между параллельными и веерно-параллельными данными. Для этого используют функций fan2para и para2fan.

Ifan1=ifanbeam(F1, D, 'FanSensorSpacing', dsensor1, 'OutputSize', output_size);

figure, imshow(Ifan1)

Ifan2=ifanbeam(F2, D, 'FanSensorSpacing', dsensor2, 'OutputSize', output_size);

figure, imshow(Ifan2)

Ifan3=ifanbeam(F3, D, 'FanSensorSpacing', dsensor3, 'OutputSize', output_size);

figure, imshow(Ifan3)