clc;
clear?all;
I0=imread（‘lena.jpg‘）;
I=double（rgb2gray（I0））;
figureimshow（I0[]）;title（‘原圖像‘）;
[mn]=size（I）;

%?%%?使用內部函數
%?theta=1:180;%計算投影角度向量
%?[Rxp]=radon（Itheta）;%取投影
%?figureimshow（R[]）;

%%?使用步長法
result=projection（I）;
figureimshow（result[]）;title（‘投影圖像‘）;

%%
d=1;?????%探測器的距離
n=301;???%探測器的個數
[M?N]=size（result）;

%設計濾波器R_Lfilter（NL）
delta=1;???%對濾波函數進行離散化單位量（采樣間隔）
for?i=2:2:2*M;???%偶數項
????h（i）=0;
end
k=1/delta/delta;
h（M）=k/4;
for?i=1:2:M-1????%奇數項
????down=-k/（i*i*pi*pi）;
????h（M+i）=down;
????h（M-i）=down;%%
end
g1=conv2（result‘h‘same‘）;????%得到固定角度的卷積180x301
g1=g1‘;?%301x180
[PQ]=size（g1）;
figureimshow（g1[]）;title（‘投影濾波后圖像‘）;

fbp?=?zeros（200200）;

delta_theta?=?pi?/?180;
for?i=1:200
????for?j=1:200
????????num=1;
??????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-07-02?10:19??考試用的程序\
?????文件???????17939??2015-06-19?10:34??考試用的程序\lena.jpg
?????文件??????105182??2015-06-24?13:13??考試用的程序\pingxingshu-chongjian.bmp
?????文件????????1523??2015-06-24?14:01??考試用的程序\pingxingshu_chongjian.m
?????文件?????????844??2015-06-24?14:03??考試用的程序\projection.m
?????文件??????298280??2015-06-24?10:27??考試用的程序\result.txt
?????文件????????1475??2015-06-25?17:54??考試用的程序\shanshu_chongjian.m
?????文件????????1213??2015-06-25?17:33??考試用的程序\shanshu_projection.m
?????文件??????106326??2015-06-24?13:13??考試用的程序\yuantu.bmp
?????文件?????????162??2015-06-25?21:18??考試用的程序\~$平行束與扇束重建算法實現.docx