An automated profile-based depth estimation method wherein anomaly type and location are identified by cross-correlation of the observed magnetic profile with theoretical anomalies. The depth to a dike-like or plate-like source is then estimated from parameters relating source body half-width, depth, and data sampling interval.