%%?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

%%?開始迭代??每考慮一個投影角