A novel method is proposed for the direct and simultaneous estimation of multiple phase derivatives corresponding to strain and slope fields from a single moiré fringe pattern in digital holographic moiré. The interference field in a given row/column is a multicomponent complex exponential signal and is represented as a spatially-varying autoregressive (SVAR) process. The spatially-varying coefficients of the SVAR model are computed by approximating them as the linear combination of linearly independent basis functions. Further, the spatially varying poles of the transfer function corresponding to the SVAR model are computed which provide the accurate estimation of the multiple phase derivatives. The simulation and experimental results are provided to substantiate the effectiveness of the proposed method.