matlab中自带的CT投影与反投影函数(扇束、平行束及其之间的转化)

最新推荐文章于 2026-06-23 14:06:10 发布

原创

最新推荐文章于 2026-06-23 14:06:10 发布 · 1.5w 阅读

119

本文深入探讨MATLAB中CT成像的关键函数，包括扇束和平行束投影与反投影，详细介绍函数特点及参数设置，如fanbeam、ifanbeam、radon和iradon等，展示不同滤波器对图像质量的影响。

2020-1-27
作者：老李

matlab中自带的CT投影与反投影函数(扇束、平行束及其之间的转化)

我参考了这篇文章，以及mathwork文档，这些给我帮助很大
链接: linkhttps://blog.csdn.net/nanhuaibeian/article/details/95234813

matlab中的函数的特点

matlab里面自带的反投影函数中，可以去选择滤好几种波器，这是非常棒的
在这里插入图片描述
而且有着投影模式转换的函数（完成了数据重排的工作）

扇束投影

函数 fanbeam 可以生成扇形射线束投影：

%% function test fanbeam
%1-26
% B = fanbeam(g,D,param1,val1,parm2,val2,...)
% g：是包含被投影的物体的图像
% 
% D：是从扇形射线束的顶点到旋转中心的距离（单位为像素）
%    规定D大于g的直径的一半，D = K*sqrt(size(g,1)^2+size(g,2)^2)/2
%    其中，K 是大于1的常数(例如，K=1.5到2是合理的值)
% B 的每一列包含扇形射线束传感器每旋转一个角度得到的样本。
% B 中的列数由扇形旋转增量决定。默认情况下B 有 360 列，B中的行数由传感器的数目决定。
% 函数 fanbeam 通过计算对于任意旋转角度覆盖全部图像所需要的射线条数，来求出传感器的数量
% 这个数字强烈的依赖于指定的几何形状(直线或圆弧)
close all;clear;clc;

对于扇束投影方式，我们可以分别生成等弧度的探测器和等长度的探测器。
在这里插入图片描述
对于等长度的探测器，他的光子接收器的间距长度是由像素长度进行度量的

我们生成图像并用两种模式来进行投影，而投影的方式是由’FanSensorGeometry’这一个参数来决定的（‘line’,‘arc’），剩下的参数可以在上图中对应找到

%% fanbeam test
% 生成 Shepp-Logan 头部幻影
g2 = phantom('Modified Shepp-Logan',600);
%

最低0.47元/天解锁文章