Cleaning the files

Author: hanyoseob

lib/backprojection.m

@@ -1,6 +1,6 @@
 % Please read the Ch.3 Image Reconstruction
 
-% Implementation for backprojection operator based on Ch. 3 Equation (3.22)
+% Implementation for backprojection operator based on Ch.3 Equation (3.22)
 % Backprojection operator is implemented using pixel-driven method
 function pdX = backprojection(pdY, param, bfig)
 
@@ -15,7 +15,7 @@
 pdX         = zeros(param.nImgY, param.nImgX, 'like', pdY);
 dCurX       = zeros(param.nImgY, param.nImgX, 'like', pdY);
 
-% Ch. 3 Equation (3.22)
+% Ch.3 Equation (3.22)
 % Backprojection operator
 for iview = 0:param.nView-1
 

lib/filterProjections.m

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

lib/filtering.m

@@ -23,7 +23,7 @@
 
 pdY         = zeros(param.nDctX, param.nView, 'like', pdX);
 
-% Ch. 3 Equation (3.6)
+% Ch.3 Equation (3.6)
 % Projection operator
 for iview = 0:param.nView-1
 
@@ -55,7 +55,7 @@
         dCurY       = 0;
 
 
-        % Ch. 3 Equation (3.6)
+        % Ch.3 Equation (3.6)
         % Line Integration along the s-axis
         for ismp = 0:nSample-1
 

lib/id2pos.m

@@ -7,4 +7,5 @@
         y(i+1) = y(i+1) + ker(((n-1) + i+1) - (j+1) + 1)*x(j+1);
     end
 end
+
 end

lib/pos2id.m

@@ -0,0 +1,6 @@
+% Transform from Index in [0, nDct] to Position in [-(dDct*nDct)/2, +(dDct*nDct)/2]
+function pos = id2pos(id, d, n)
+
+pos = (id - (n - 1)/2)*d;
+
+end

lib/projection.m

@@ -1,3 +1,6 @@
+% Reference
+% Bilinear Interpolation
+% https://en.wikipedia.org/wiki/Bilinear_interpolation
 function dst  = interpolation1d(src, curid)
 
 preid	= floor(curid);

lib/convolution1d.m renamed to util/convolution1d.m

@@ -1,4 +1,6 @@
-
+% Reference
+% Linear Interpolation
+% https://en.wikipedia.org/wiki/Linear_interpolation
 function dst  = interpolation2d(src, curidy, curidx)
 
 preidy	= floor(curidy);