滤波反投影matlab实现matlab源程序

滤波反投影

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