function [rays,phis] = extractRays(I,start,increment,num)

szI = size(I);
phi = start;
sz = size(I);
rays = cell(num,1);
for i=1:num
  phis(i)=phi;
  rI = imrotate(I,phi);
  rays{i} = sum(rI,1);
  phi = phi+increment;
end

%ct(sz(1),sz(2))=0;
%for i=1:num
%  phis(i)=phi;
%  rI = imrotate(I,phi);
%  cropped = rI;
%  dimensions=size(rI);
  %Crop(rI,size(I));
%  rays = [rays;sum(cropped,1)];
%  ct = ct + Crop(backProject(rays,dimensions,phi),size(I));
%  phi = phi+increment;
%end

%for i=1:num
%  phis(i)=phi;
%  rI = imrotate(I,phi);
%  cropped = rI(1:szI(1),1:szI(2));
%  rays = sum(rI,1);
%  ct = ct + imrotate(ones(length(rays),1)*rays,-phi);
%  phi = phi+increment;
%end