ART迭代算法图像重建

%% ART based on FLA
clear all
close all
clc;
L = 128;%图像大小
nIter = 0;% 迭代次数
iterOK = 0;
iter0 = ones（4*L）*0.2;%元素为0.2的矩阵
theta0 = [0:14; 15:29; 30:44; 45:59; 60:74; 75:89; 90:104; 105:119; 120:134; 135:149; 150:164; 165:179];
theta = reshape （theta0 180 1）;%theta = [01530...16511631...16621732...167......]‘
nTheta = length（theta）;% 180

%% 生成投影，以及初始迭代矩阵
I = phantom （L）;% 原始图像

R1 = zeros（128128）;% 迭代初值为0

% P1 = radon （I 0:30:179）;
% R1 = iradon （P1 0:30:179 ‘linear‘ ‘hamming‘）;% 迭代初始值为R1（效果好）

for aa = 1:L %%
    for bb = 1:L
        iter0 （（4*aa-3）:4*aa （4*bb-3）:4*bb） = R1 （aa bb）/16;
        I1 （（4*aa-3）:4*aa （4*bb-3）:4*bb） = I （aa bb）/16;% 把原图扩大4倍
    end
end

%% 开始迭代  每考虑一个投影角