24小时热门版块排行榜

返回列表

zhuningzh

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 22.9
帖子: 36
在线: 31.5小时
虫号: 1680108
注册: 2012-03-10
专业: 有机分子功能材料化学

[求助] 已经有投影数据,怎么用Matlab重构Shepp-Logan图像

模型：Shepp-logan头模型，该模型限制在半径为1mm的圆内，该模型共包含了10个椭圆，每个椭圆内的数值大小（代表线性衰减系数），及椭圆分布见下图：
投影数据存储格式：float型数据，文件后缀为“.prj”，是没有文件头的纯投影数据，即：假设投影角度个数为M，探测器单元个数为N，则文件长度为M*N*4个字节。按照角度旋转采集的顺序依次排列，即每读入N个数据，表示一个投影角度下的投影，直到M个投影角度数据完全读入程序。

1、平行束投影数据，0-359度，每度一个投影，即角度间隔为1度，共360个投影；线阵列探测器的长度为2mm，共有256个探测单元，每个探测单元大小为2/256mm。

2、扇形束等距投影数据，0-359度，每度一个投影，即角度间隔为1度，共360个投影；X射线源到探测器的距离为8mm，X射线源到旋转中心的距离为4mm，线阵列探测器的长度为4mm，共有256个探测单元，每个探测单元大小为4/256mm。

3、扇形束等角投影数据，0-359度，每度一个投影，即角度间隔为1度，共360个投影；X射线源到探测器的距离为8mm，X射线源到旋转中心的距离为4mm，探测器阵列对应扇角30度，共有256个探测单元，每个探测单元对应的角度为30/256度。
哪位大侠帮忙写下matlab代码，比较急，谢啦！

回复此楼