CT成像的代码

压缩包里有着关于CT成像的代码，是利用matlab进行实现。

clear  
load data.mat ;                                             %下载数据
H = [0:floor（size（data1）/2）floor（size（data1）/2）-1:-1:0]; %补偿高频的滤波器
H = H（:） ;
l = pow2（nextpow2（size（data1））-1） ;                        %重构图像的大小
F = zeros（ll） ;                                            %用于存放重构后的图像
f = zeros（size（data））;
for i = 1 : size（data2） 
    x = data（:i） ;                                         
    y = fft（x） ;                                         
    f（:i） = real（ifft（y）） ;                                      
    tmp = imrotate（ repmat（f（:i）1size（data1））i-1‘bilinear‘ ） ; %将每一行的傅立叶反变换由一维扩展到二维
    tmp = tmp（floor（size（tmp1）/2-l/2）+1:floor（size（tmp1）/2+l/2）floor（size（tmp2）/2-l/2）+1:floor（size（tmp2）/2+l/2）） ; 
    F = F + tmp ;                           

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件      844368  2011-05-14 08:31  data.mat
     文件       48129  2011-05-22 23:13  xct.jpg
     文件         924  2011-05-22 23:27  xct.m