Many edge detection methods, based on horizontal and vertical derivatives, have been introduced to provide us with intuitive information about the horizontal distribution of a subsurface anomalous body. Understanding the characteristics of each edge detection method is important for selecting an optimized method. In order to compare the characteristics of the individual methods, this study applied each method to synthetic magnetic data created using homogeneous prisms with different sizes, the numbers of bodies, and spacings between them. Seven edge detection methods were comprehensively and quantitatively analyzed: the total horizontal derivative (HD), the vertical derivative (VD), the 3D analytic signal (AS), the title derivative (TD), the theta map (TM), the horizontal derivative of tilt angle (HTD), and the normalized total horizontal derivative (NHD). HD and VD showed average good performance for a single-body model, but failed to detect multiple bodies. AS traced the edge for a single-body model comparatively well, but it was unable to detect an angulated corner and multiple bodies at the same time. TD and TM performed well in delineating the edges of shallower and larger bodies, but they showed relatively poor performance for deeper and smaller bodies. In contrast, they had a significant advantage in detecting the edges of multiple bodies. HTD showed poor performance in tracing close bodies since it was sensitive to an interference effect. NHD showed great performance under an appropriate window.