CUMCM-2017A CT系统参数标定及成像matlab代码

CUMCM-2017A CT系统参数标定及成像前三题matlab代码。可直接运行。

R = xlsread（‘A题附件.xls‘2）; 
W = 512;   
theta = 1:180; 
fourier = fft（R W）; 
filter = 2*[0:（W/2-1） W/2:-1:1]‘/W;   
lvbo = zeros（W180）;  
for i = 1:180   
     lvbo（:i） = fourier（:i）.*filter;   
end

inversse = real（ifft（lvbo））;   
N = 360;
final = zeros（N）; 
for i = 1:180   
rad = theta（i）*pi/180; 
 for x = （-N/2+1）:N/2   
 for y = （-N/2+1）:N/2 
t = round（x*cos（rad+pi/2）+y*sin（rad+pi/2））;   
final（x+N/2y+N/2）=final（x+N/2y+N/2）+inversse（t+round（size（R1）/ 2）i）;
 end
 end
end
final=final/180;
save date final;
imshow（final）title（‘final‘）
load date final

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件     4231680  2017-08-09 10:09  A题附件.xls
     文件         594  2018-08-28 21:14  wen1.m
     文件        1197  2018-08-28 21:15  wen2.m
     文件        1225  2018-08-28 21:16  wen3.m