clc;
clear?all;
close?all;
N?=?180;?%圖像大小
N2?=?N^2;
I?=?phantom（N）;%?產生頭模型圖像
theta?=?linspace（0180181）;
theta?=?theta（1:180）;
%%?=?=?=?=?=?=?產生投影數據?=?=?=?=?=?=?%%
P_num?=?260;?%探測器通道個數
P?=?ParallelBeam（theta?N?P_num）;?%產生投影數據
%P?=?radon（Itheta）;
%%?=?=?=?=?=?=?獲取投影矩陣?=?=?=?=?=?=%%
delta?=?1;%?網格大小，角度增量
[W_indW_dat]?=?SystemMatrix（thetaNP_numdelta）;
%%?=?=?=?=?=?=?進行ART迭代?=?=?=?=?=?=?%%
F?=?zeros（N21）;??%初始圖像向量
lambda?=?0.25;??%松弛因子
c?=?0;??%迭代計數器
irt_num?=?5;??
while（c????for?j?=?1:length（theta）
????????for?i?=?1:1:P_num
????????????%?取得一條射線所穿過的網格編號和長度
????????????u?=?W_ind（（j-1）*P_num?+?i:）;??%?編號
????????????v?=?W_dat（（j-1）*P_num?+?i:）;??%?長度
????????????
????????????if?any（u）?==?0
????????????????continue;
????????????end
????????????%恢復投影矩陣中與這一條射線對應的行向量?w
????????????w?=?zeros（1N2）;
????????????ind?=?u?>?0;
????????????w（u（ind））=v（ind）;
????????????%?圖像進行一次迭代
????????????PP?=?w?*?F;???%前向投影
????????????C?=?（P（ij）-PP）/sum（w.^2）?*?w‘;??%?修正項???????????
????????????F?=?F?+?lambda?*?C;
????????end
????end
????F（F<0）?=?0;???%?小于0的像素值置為0
????c?=?c+1;
end
F?=?reshape（FNN）‘;?%?轉換成N?x?N的圖像矩陣
%%?=?=?=?=?=?=?仿真結果顯示?=?=?=?=?=?=?%%
figure（1）;
imshow（I）;xlabel（‘（a）180x180頭模型圖像‘）;
figure（2）;
%A?=?imadjust（F）;
imshow（F[]）;xlabel（‘（b）ART算法重建的圖像‘）;
?????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-12-09?21:12??第一次實驗\
?????文件?????????889??2015-12-01?12:03??第一次實驗\Backprojection.m
?????文件????????1291??2015-12-03?10:21??第一次實驗\ParallelBeam.m
?????文件?????????345??2015-12-02?11:42??第一次實驗\RLfilter.m
?????文件?????????981??2015-12-01?12:04??第一次實驗\RLfilteredbackprojection.m
?????文件?????????268??2015-12-02?11:43??第一次實驗\SLfilter.m
?????文件?????????935??2015-12-01?12:04??第一次實驗\SLfilteredbackprojection.m
?????文件????????6830??2015-12-04?11:39??第一次實驗\SystemMatrix.m
?????文件????????1492??2015-12-05?10:39??第一次實驗\art.m
?????文件?????????575??2015-12-03?10:43??第一次實驗\lb_irad.m
?????文件????????1168??2015-12-04?16:42??第一次實驗\lb_irad_fh.m
?????文件?????????304??2015-12-03?09:24??第一次實驗\zj_irad.m
?????文件?????????119??2015-12-09?21:11??第一次實驗\說明.txt