This paper proposes a recursive procedure, called the extended local rank factorization (elrf), that characterizes the order of the pole and the coefficients of the Laurent series representation of the inverse of a regular analytic matrix function around a given point. The elrf consists in performing a finite sequence of rank factorizations of matrices of nonincreasing dimension, at most equal to the dimension of the original matrix function. Each step of the sequence is associated with a reduced rank condition, while the termination of the elrf corresponds to a full rank condition; this last step reveals the order of the pole. The Laurent coefficients $B_{n}$ are calculated recursively as $B_{n} = C_{n} + \sum_{k=1}^{n} D_{k} B_{n-k}$, where $C_{n}$, $D_{k}$ have simple closed form expressions in terms of the quantities generated by the elrf. It is also shown that the elrf characterizes the structure of Jordan pairs, Jordan chains, and the local Smith form. The procedure is easily cast in an algorithmic ...