CT-FBP 用MATLAB实现了CT平行束滤波反投影和扇束滤波反投影重建过程

用MATLAB实现了CT平行束滤波反投影和扇束滤波反投影重建过程，效果良好，其中radon变换的方法采用步长法求解，较新颖。

clc;
clear all;
I0=imread（‘lena.jpg‘）;
I=double（rgb2gray（I0））;
figureimshow（I0[]）;title（‘原图像‘）;
[mn]=size（I）;

% %% 使用内部函数
% theta=1:180;%计算投影角度向量
% [Rxp]=radon（Itheta）;%取投影
% figureimshow（R[]）;

%% 使用步长法
result=projection（I）;
figureimshow（result[]）;title（‘投影图像‘）;

%%
d=1;     %探测器的距离
n=301;   %探测器的个数
[M N]=size（result）;

%设计滤波器R_Lfilter（NL）
delta=1;   %对滤波函数进行离散化单位量（采样间隔）
for i=2:2:2*M;   %偶数项
    h（i）=0;
end
k=1/delta/delta;
h（M）=k/4;
for i=1:2:M-1    %奇数项
    down=-k/（i*i*pi*pi）;
    h（M+i）=down;
    h（M-i）=down;%%
end
g1=conv2（result‘h‘same‘）;    %得到固定角度的卷积180x301
g1=g1‘; %301x180
[PQ]=size（g1）;
figureimshow（g1[]）;title（‘投影滤波后图像‘）;

fbp = zeros（200200）;

delta_theta = pi / 180;
for i=1:200
    for j=1:200
        num=1;
      

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-07-02 10:19  考试用的程序\
     文件       17939  2015-06-19 10:34  考试用的程序\lena.jpg
     文件      105182  2015-06-24 13:13  考试用的程序\pingxingshu-chongjian.bmp
     文件        1523  2015-06-24 14:01  考试用的程序\pingxingshu_chongjian.m
     文件         844  2015-06-24 14:03  考试用的程序\projection.m
     文件      298280  2015-06-24 10:27  考试用的程序\result.txt
     文件        1475  2015-06-25 17:54  考试用的程序\shanshu_chongjian.m
     文件        1213  2015-06-25 17:33  考试用的程序\shanshu_projection.m
     文件      106326  2015-06-24 13:13  考试用的程序\yuantu.bmp
     文件         162  2015-06-25 21:18  考试用的程序\~$平行束与扇束重建算法实现.docx