zz 图像数据投影&投影数据重建图像 radon

Radon变换

图形处理工具箱中的radon函数计算指定方向上图像矩阵的投影。二维函数f(x,y)的投影是一组线积分。radon函数计算一定方向上平行光束的线积分。光线间隔1个像素单位。为了表示图像，radon函数通过围绕图像的中心旋转光源来从不同的角度获得图像的平行光投影。
%下面的命令行计算和显示一幅包含一个方形对象的图像在0度和45度上的radon变换
I=zeros(100,100);
I(25:75,25:75)=1;
imshow(I);
%用下面的命令行进行变换
[R,xp]=radon(I,[0 45]);
figure;plot(xp,R(:,1));title('R_{0^o}(x/prime)');
figure;plot(xp,R(:,2));title('R_{45^o}(x/prime)');
%当角度很多时，radon变换的结果常常用图像显示。
theta=0:180;
[R,xp]=radon(I,theta);
imagesc(theta,xp,R);
title('R_{/theta}(X/prime)');
xlabel('/theta(degrees)');
ylabel('/X/prime');
set(gca,'XTick',0:20:180);
colormap(hot);
colormap;
利用投影数据重建图像
下面代码演示了如何利用平行投影数据重建图像。测试图像是Shepp-Logan头部剖视图，可以用工具箱中的phantom函数生成。剖视图中外部的白色椭圆形壳可看着头骨，内部的许多椭圆形可看着脑部组织。
%生成Shepp-Logan头部剖视图图像
P=phantom(256);
imshow(P);
%计算三套theta值对应剖视图的radon变换
theta1=0:10:170;[R1,xp]=radon(P,theta1);
theta2=0:5:175;[R2,xp]=radon(P,theta2);
theta3=0:2:178;[R3,xp]=radon(P,theta3);
%显示一幅Shepp-Logan头部剖视图的radon变换图。显示第三套theta值对应的图像
figure;imagesc(theta3,xp,R3);colormap(hot);colorbar;
xlabel('/theta');ylabel('x/prime');
%根据第二部生成的投影数据重建头部试图图像
I1=iradon(R1,10);
I2=iradon(R2,5);
I3=iradon(R3,2);
figure;imshow(I1);
figure;imshow(I2);
figure;imshow(I3);

fanbeam变换

扇形光束投影数据以及用扇形光束投影数据重建图像
P=phantom(256);
imshow(p);
%用FanSensorSpacing参数修改传感器间隔，本例采用弧形布置传感器，所以通过指定光束的角度间隔来指定传感器之间的间隔（单位度）
D=256;
dsensor1=2;
F1=fanbem(P,D,'FanSensorSpacing',dsensor1);
dsensor2=1;
F2=fanbem(P,D,'FanSensorSpacing',dsensor2);
dsensor3=0.25;
F3=fanbem(P,D,'FanSensorSpacing',dsensor3);
%用ifanbeam函数根据投影数据重建图像
out_size=max(size(P));
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);