In this paper the application of an edge detection technique to gravity data is described. The technique is based on the tilt angle map (TAM) obtained from the first vertical gradient of a gravity anomaly. The zero contours of the tilt angle correspond to the boundaries of geologic discontinuities and are used to detect the linear features in gravity data. I also present that the distance between zero and +/-pi/4 pairs obtained from the TAM corresponds to the depth to the top of the vertical contact model. Alternatively, the half distance between -pi/4 and +pi/4 radians is equal to the depth to the same model. I illustrate the applicability of the present method by gravity data due to buried vertical prisms, imaging the positions of the edges of the prisms. The results obtained from the theoretical data, with and without random noise, have been discussed. The analysis of the TAM has been demonstrated on a field example from the KozaklA +/--Central Anatolian region, Turkey, and the location and depth of the edges of the structural uplifts of the KozaklA +/- graben are imaged. The results indicated that depth values from these sources have ranged between 0.2 and 0.6 km. I have also compared the Euler deconvolution technique with the TAM images obtained from the first vertical gradient of residual gravity anomaly. Both techniques have agreed closely in detecting the horizontal location and depth of the uplift edges in the subsurface with good precision.