This thesis presents a robust method for estimating the relaxations of a metallic object from its electromagnetic induction (EMI) response. The EMI response of a metallic object can be accurately modeled by a sum of real decaying exponentials. However, it is diffcult to obtain the model parameters from measurements when the number of exponentials in the sum is unknown or the terms are strongly correlated. Traditionally, the time constants and residues are estimated by nonlinear iterative search that often leads to unsatisfactory results. In this thesis, a constrained linear method of estimating the parameters is formulated by enumerating the relaxation parameter space and imposing a nonnegative constraint on the parameters. The resulting algorithm does not depend on a good initial guess to converge to a solution. Using tests on synthetic data and laboratory measurement of known targets the proposed method is shown to provide accurate and stable estimates of the model parameters.