We present an automatic method for extracting the geologic features of gravity and magnetic anomaly bodies. Combining the curvature, analytic signal and upward continuation processing together, we present a depth-estimation method for interpretation of sources. The method only uses the first derivatives of the gravity or magnetic anomaly reducing to the pole, which is insensitive to noise and has higher precision comparing with the process of second derivatives. Using depth-estimation, the geometric types of sources can also be inferred. Tests on synthetic and real data excellent results in extracting geologic features.