RWithRotAndTrans.m

function sinogram = RWithRotAndTrans( img, rotations, translations_pix, ...
  nDet, R )
  % rotations is an N element array specifying the rotation in radians
  %   corresponding to each sinogram projection
  % nDet is the number of detectors

  [M,~] = size(R);
  nThetas = M / nDet;
  sinogram = zeros(nThetas,nDet);

  [nTranslations,~] = size(translations_pix);
  if nTranslations ~= nThetas, error('Input error'); end;

  for i=1:nThetas
    thisRot = rotations(i);
    rotated = imrotate( img, thisRot, 'bilinear', 'crop' );

    thisTrans_pix = translations_pix(i,:);
    translated = translateImg( rotated, thisTrans_pix );

    %sinogram(i,:) = R(nDet*(i-1)+1:nDet*i,:) * translated(:);
    tmpSino = R * translated(:);
    tmpSino = reshape( tmpSino, [nThetas nDet] );
    sinogram(i,:) = tmpSino(i,:);
  end

end