An automatic depth estimation method in which the source parameters are estimated from the autocorrelation function of the magnetic anomaly. Like Werner deconvolution, the method uses a dike or contact model.