clc;?clear?all;?close?all;
P=imread（‘xinzang.jpg‘）;
g2=rgb2gray（P）;
subplot（221）;?imshow（g2）;title（‘原始圖像‘）;
theta?=?0:0.5:179.5;%傳感器參數（角度數）
%執行雷登變換
[R2?xp2]?=?radon（g2?theta）;
%顯示投影圖像（正弦圖）
r2_show?=?flipud（R2‘）;
subplot（222）;imshow（r2_show?[]?‘XData‘?xp2（[1?end]）?‘YData‘?[179.5?0]）;title（‘原始圖像的雷登變換‘）;
axis?xy;
axis?on;
xlabel（‘\rho‘）;
ylabel（‘\theta‘）;
%從正弦圖得到反投影圖
w2=iradon（R2?theta?‘none‘）;?%原始圖像的反投影（未濾波）
subplot（223）;imshow（w2?[]）;title（‘原始圖像的反投影（未濾波）‘）;
f2?=?iradon（R2?theta?‘Hamming‘）;?%原始圖像的hamming濾波反投影
subplot（224）;imshow（f2?[]）;title（‘原始圖像的反投影（漢明濾波）‘）;
%使用扇形射線束
D?=?1.5*hypot（size（g21）size（g22））/2;
b2_line=fanbeam（g2?D?‘FanSensorGeometry‘?‘line‘?‘FanSensorSpacing‘?1?‘FanRotationIncrement‘?0.5）;
b2_line_s=flipud（b2_line‘）;
b2_src=fanbeam（g2?D?‘FanSensorGeometry‘?‘arc‘?‘FanSensorSpacing‘?.08?‘FanRotationIncrement‘?0.5）;
figure;
subplot（121）;imshow（g2）;%?原始圖像
subplot（122）;?imshow（b2_line_s?[]?‘XData‘?[0?850]?‘YData‘?[0?360]）;
axis?xy;
axis?on;
ylabel（‘扇形旋轉角度‘）;
xlabel（‘傳感器個數‘）;
%扇形反投影濾波重建
B1?=?fanbeam（g2?D）;
fB1=ifanbeam（B1?D）;
B2?=?fanbeam（g2?D?‘FanSensorSpacing‘?.05?‘FanRotationIncrement‘?.5）;
fB2=ifanbeam（B2?D?‘FanSensorSpacing‘?.05?‘FanRotationIncrement‘?.5?‘filter‘?‘Hamming‘）;
figure;
subplot（121）;imshow（fB1?[]）;
subplot（122）;imshow（fB2?[]）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1505??2017-04-10?16:35??fbprec.m
?????文件?????????621??2017-04-10?16:31??filter.m
?????文件???????15914??2017-04-09?15:16??xinzang.jpg