function dQ = calculateDifferentialImages(Q);
        
    sizeQ = size(Q); N = sizeQ(3);
    dQ    = zeros(sizeQ(1), sizeQ(2), N-1);
    
    for tp = N:-1:3
        dQ(:,:,tp-1) = Q(:,:,tp) - Q(:,:,tp-1);
    end
        
    dQ(:,:,1) = Q(:,:,2);
    
  