shepp-logan头模型锥束投影生成程序

基于matlab工具，对医学影像中最常用的shepp-logan模型进行投影仿真，射束为锥束，采用射线与椭球体求交的算法生成锥束射线下头模型的投影图像

clc;
clear all;
close all;
N = 256;
SOD = 44;%源到旋转中心的距离
SDD = 2 * SOD;%源到探测器中心的距离
vitual_detector_length = 2;%虚拟探测器长度
detector_length = vitual_detector_length * SDD / SOD;
theta_num = 360;%旋转角度数
detector_channel_size = detector_length / N;%探测器的单位长度
P = zeros（N N theta_num）;%存放投影数据
%%===================定义头模型参数=====================%%
% xe     ye     ze     ae     be     ce     phi    gamma    theta     rho
shep = [
  0      0      0  0.6900  0.9200 0.8100    0       0        0     1.0000
  0  -0.0184    0  0.6624  0.8740 0.7800    0       0        0    -0.8000
  0.2200 0      0  0.1100  0.3100 0.2200  -18.0000  0     10.0000 -0.2000
 -0.2200 0      0  0.1600  0.4100 0.2800   18.0000  0     10.0000 -0.2000
  0  0.3500 -0.1500 0.2100 0.2500 0.4100    0       0        0     0.1000
  0  0.1000 0.2500 0.0460  0.0460 0.0500    0       0        0     0.1000
  0 -0.1000 0.2500 0.0460  0.0460 0.0500    0       0        0     0.1000
-0.0800 -0.6050 0  0.0460  0.0230 0.0500  

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2019  2018-05-13 21:45  head_phantom.m

     文件       3123  2018-05-13 21:31  medfunc3DProjectHeadModel.m

----------- ---------  ---------- -----  ----

                 5142                    2