This Letter proposes a method to estimate phase derivatives of arbitrary order in digital holographic interferometry. Based on the desired order, the generalized complex-lag distribution is computed from the reconstructed interference field. Subsequently, the phase derivative is estimated by tracing the peak of the distribution. Simulation and experimental results are presented to validate the method’s potential.