FDKFDK 三维图像重建算法

clc;
clear all;
clear all figure;

load image2.mat
Img = （image2）;
Img2 = permute（Img[3 2 1]）;
Img3 = permute（Img[3 1 2]）;

% Vplan=[256256256];
% Img=single（phantom3d（Vplan））;



figure（1）; imagesc（Img（::end/2））; axis off; axis equal; colormap gray; colorbar;
% figure（2）; imagesc（fliplr（Img（::end/2）））; axis off; axis equal; colormap gray; colorbar;
% figure（3）; imagesc（rot90（Img（::end/2）2））; axis off; axis equal; colormap gray; colorbar;
% figure（4）; imagesc（Img2（::end/2））; axis off; axis equal; colormap gray; colorbar;
% figure（5）; imagesc（Img3（::end/2））; axis off; axis equal; colormap gray; colorbar;

Voxelsize = size（Img）;
fprintf（‘the size of Voxel is %d %d %d\n‘Voxelsize）;

nx = Voxelsize（1）;
ny = Voxelsize（2）;
nz = Voxelsize（3）;

sx = 256; % in mm
sy =256; % in mm
sz =32*1.25; % in mm

% Detector setting according to Varian Trilogy OBI （real size）
su = 256*2;	% in mm
sv = 256*2;	% in mm

%The real detector panel pixel density （number of pixels）
nu = 256;		
nv = 256;

%