update projection function

Author: hanyoseob

demo_parallelbeam_ct.m

@@ -10,7 +10,7 @@
 % nView	: # of the views [unit]
 % dView	: Gap between view_(k) - view_(k-1) [degree]
 param.dAngle        = 180;  % degree
-param.nView         = 180; 	% # of unit
+param.nView         = 90; 	% # of unit
 param.dView         = param.dAngle/param.nView;  % degree
 param.pdView        = linspace(0, param.dAngle, param.nView + 1);
 param.pdView(end)   = [];

lib/projection.m

@@ -11,17 +11,17 @@
 dSizeImgX   = param.nImgX*param.dImgX;
 dSizeImgY   = param.nImgY*param.dImgY;
 
-dDiameter   = sqrt(dSizeImgX*dSizeImgX + dSizeImgY*dSizeImgY);
+dDiameter   = sqrt(dSizeImgY*dSizeImgY + dSizeImgX*dSizeImgX);
 dRadius     = 0.5*dDiameter;
 
-dSample     = min(param.dImgX, param.dImgY);
+dSample     = min(param.dImgY, param.dImgX);
 nSample     = ceil(dDiameter/dSample);
 
 pdY         = zeros(param.nDctX, param.nView, 'like', pdX);
 
 for iview = 0:param.nView-1
 
-    dBeta   = -iview*param.dView;
+    dBeta   = iview*param.dView;
 
     for idctx = 0:param.nDctX-1
 
@@ -30,15 +30,15 @@
 
         dGamma      = 0;
 
-        dOriNorDirY = dSample*cosd(dGamma);
         dOriNorDirX	= dSample*sind(dGamma);
+        dOriNorDirY = dSample*cosd(dGamma);
 
-        dNorDirX	= cosd(dBeta)*dOriNorDirX - sind(dBeta)*dOriNorDirY;
-        dNorDirY	= sind(dBeta)*dOriNorDirX + cosd(dBeta)*dOriNorDirY;
-        
-        dPosImgX    = cosd(dBeta)*dOriImgX - sind(dBeta)*dOriImgY;
-        dPosImgY    = sind(dBeta)*dOriImgX + cosd(dBeta)*dOriImgY;
+        dNorDirX	= cosd(dBeta)*dOriNorDirX + sind(dBeta)*dOriNorDirY;
+        dNorDirY	= -sind(dBeta)*dOriNorDirX + cosd(dBeta)*dOriNorDirY;
 
+        dPosImgX    = cosd(dBeta)*dOriImgX + sind(dBeta)*dOriImgY;
+        dPosImgY    = -sind(dBeta)*dOriImgX + cosd(dBeta)*dOriImgY;
+
         dCurY       = 0;
 
         for ismp = 0:nSample-1