An automated profile-based depth estimation method derived from S. Werner's analysis of magnetic anomalies from sheet-like bodies. Polynomials representing a total field anomaly or its derivative (horizontal gradient) can be simultaneously solved to estimate the depth, dip, horizontal location, and susceptibility of the source body (thin sheet or interface).