function ROT = Rotationnel(II);
  [N1,N2,N3] = size(II);
  T1  = [2:N1 1]; T2 = [2:N2 1]; 
  ROT = (II(:,T2,2)-II(:,:,2))- (II(T1,:,1)-II(:,:,1)); 
endfunction
