function [ sino2 ] = manual_COR_correction( sino, dcentre )
%MANUAL_COR_CORRECTION manual centre of rotation correction
% Jelle van der Horst - januari 2016
[a, b]=size(sino);
p=1:a;
[XX, YY]=meshgrid(1:b,1:a);

background=mean(sino(1,:)).*ones(size(sino));

sino2=interp2(XX,YY,sino,XX,YY+dcentre);
ind=isnan(sino2);
sino2(ind)=mean(mean(background));

end

