The enhanced horizontal derivative is usually used to interpret the potential field anomalies and is a high-resolution method to detection the horizontal boundaries of gravity sources. This method is highly flexible, and different conditions of signal-to-noise ratios can be treated by an appropriate selection of the terms of the summation. In this paper, the EHD method is used to recognition sub-structures and edge detection in synthetic and real data. In order to calculation enhanced horizontal derivative value of gravity anomalies, a code has been written in MATLAB software. The first we applied it on the synthetic data and then we used enhanced horizontal derivative method to interpret the gravity data from Mobrun sulfide body.