image=?phantom（‘Modified?Shepp-Logan‘256）;
subplot（221）;
imshow（image）;
title（‘原始圖像‘）;
[Rxp]=?radon（image0:179）;
theta?=?0:180;
I1?=?iradon（R0:179‘linear‘‘none‘）;
subplot（222）;
imagesc（I1）;colormap（gray）
title（‘無濾波重