CT反投影重建算法

function [ image ] = Backprojection（ sinogram thetas ）
%% 直接反投影法功能描述
% 输入参数
% sinogram 正弦图坐标系统数据
% thetas 扫描角度向量
% 输出参数
% image 重建图像
%该函数先建立初始化重建图像，即全零矩阵image，然后取正弦图每个角度的投影，原路反投影到每条射线通过的像素点上。

%% 输入正弦图数据，测试直接反投影法函数
% close all;clear all;clc;
% folderpath=‘E:\MATLAB\CT成像_NLT\2014-09-05 23-24-24高斯‘;
% name=‘\z3‘;
% sinogram=load（strcat（folderpathname‘.txt‘））;
% thetas=linspace（0360size（sinogram1））;

%% 初始化重建图像
[numOfParallelnumOfAngular] = size（sinogram）;
image = zeros（numOfParallel numOfParallel）;                                             %初始化重建图像为全零矩阵，行数和列数等于正弦图平移数
midindex = floor（numOfParallel / 2） + 1;
[xCoords yCoords] = meshgrid（ceil（-n