update

Author: hanyoseob

demo_parallelbeam_ct.m

@@ -14,7 +14,6 @@
 param.dAngle        = 360;  % degree
 param.nView         = 360; 	% # of unit
 param.dView         = param.dAngle/param.nView;  % degree
-param.pdView        = linspace(0, param.dAngle - param.dAngle/param.nView, param.nView);
 param.DSO           = 400 ; % mm
 param.DSD           = 800;  % mm
 
@@ -45,7 +44,7 @@
 param.nImgX         = 256;  % # of elements
 
 param.dOffsetImgY	= 0;    % # of elements
-param.dOffsetImgX   = 0;    % # of elements
+param.dOffsetImgX   = 0;    % # of elements4
 
 %% Load image
 load('XCAT512.mat');
@@ -59,7 +58,7 @@
 prj         = projection(input, param);
 toc;
 
-disp ('filtering - implementation of Ch.3 Equation (3.29) & (3.30)');
+disp ('filtering - implementation of Ch.3 Equation (3.28) & (3.29) & (3.30)');
 tic;
 prj_flt  	= filtering(prj, param);
 toc;

lib/filtering.m

@@ -1,6 +1,6 @@
 % Please read the Ch. 3 Image Reconstruction
 
-% Implementation for filtering operator based on Ch. 3 Equation (3.29) & (3.30)
+% Implementation for filtering operator based on Ch. 3 Equation (3.28) & (3.29) & (3.30)
 % Filtering operator is implemented by both convolution and FFT versions.
 function [pdFltY, pdFlt] = filtering(pdY, param)