CT不完全投影数据重建算法与研究 ART与SART

CT不完全投影数据重建算法,ART与SART.CT不完全投影数据重建算法,ART与SART.CT不完全投影数据重建算法,ART与SART.CT不完全投影数据重建算法,ART与SART.CT不完全投影数据重建算法,ART与SART.

clear all
close all
tic
N = 256;
n = 0;%迭代次数
it0 = ones（4*N）;%数据矩阵扩增为原来四倍，保证对旋转过程的灵敏度
it1=it0;
sl=（0:1:100）;
slf=（0:1:99）;
angle=（0:180）‘;
nangle = length（angle）;
I = phantom（‘Modified Shepp-Logan‘256）;
%I=imread（‘C:\Users\adminraden\Desktop\医学成像系统大作业\CT图像\800px-CT_ScoutView.jpg‘）;
%I=rgb2gray（I）;
%I=double（I）;
%I = imresize（I[256 256]）;
%I=I./255;
P = radon （I angle）;
R0 = iradon （P angle）;
R0 = imresize（R0[256 256]）;
P1 = radon （I 0:10:180）;
%P1=P;
%P1（2001:181）=100;
%P1（2001:181）=（P1（1991:181）+P1（2011:181））/2;
R1 = iradon （P1 0:10:180）;
%R1 = iradon （P1 angle）;
for a = 1:N;
    for b = 1:N;
        it0 （（4*a-3）:4*a （4*b-3）:4*b） = R1 （a b）/16;
        I1 （（4*a-3）:4*a （4*b-3）:4*b） = I （a b）/16;
    end
end
while （n < 100）
    n = n + 1;
    for ii = 1:N
        pj = zeros （4*N）;
        pj （:（4*ii-3）:4*ii） = 1;
        pj1 = imrotate （pj angle （mod（n-1 nangle）+1） ‘crop‘）;%旋转采样
       
        P_pj = sum（sum（I1.*pj1））;
        R_pj = sum（sum（it0.*pj1））;

        it1 = it0 + pj1.*（P_pj-R_pj）/sum（sum（pj1））;
        it0 = it1;    
    end
%      for ii = N:2
%         pj = zeros （4*N）;
%         pj （:（4*ii-3）:4*ii） = 1;
%         pj1 = imrotate （pj angle （mod（n-1 nangle）+1） ‘crop‘）;%旋转采样
%        
%         P_pj = sum（sum（I1.*pj1））;
%         R_pj = sum（sum（it0.*pj1））;
% 
%         it1 = it0 + pj1.*（P_pj-R_pj）/sum（sum（pj1））;
%         it0 = it1;    
%     end

  for a = 1:N;
      for b = 1:N;
                R （a b） = sum （sum （it0 （ （4*a-3）:4*a （4*b-3）:4*b）））;
      end
  end
  dA=R-R0;
% slf（n）=max（max（abs（dA）））;
slf（n）=norm（dA1）;
end

% for a = 1:N;
%     for b = 1:N;
%                 R （a b） = sum （sum （it0 （ （4*a-3）:4*a （4*b-3）:4*b）））;
%     end
% end
figure
subplot（221）;imshow（I） title（‘原始图像‘）;stem（slf）;
subplot（222）;imshow（R0） title（‘完整的采集信号重构图像‘）;
subplot（223）;imshow（R1） title（‘不完整信号重构图像‘）;
subplot（224）;imshow（R） title（‘迭代重建图像‘）;

toc

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-06-15 19:39  大作业MATLAB程序\
     目录           0  2017-06-15 19:39  大作业MATLAB程序\art算法\
     文件        2071  2017-01-10 22:11  大作业MATLAB程序\art算法\main.m
     文件          45  2017-01-10 21:53  大作业MATLAB程序\art算法\Untitled111.m
     文件         684  2016-12-22 00:51  大作业MATLAB程序\art算法\Untitled2.m
     文件         202  2016-12-28 00:22  大作业MATLAB程序\art算法\Untitled3.m
     目录           0  2017-06-15 19:39  大作业MATLAB程序\sart算法\
     文件         719  2016-12-24 23:29  大作业MATLAB程序\sart算法\CWM.m
     文件         609  2017-01-10 22:29  大作业MATLAB程序\sart算法\main.m
     文件         445  2016-12-24 23:23  大作业MATLAB程序\sart算法\sart.m
     文件         737  2016-12-24 23:17  大作业MATLAB程序\sart算法\W.m