One of the difficulties in interpreting magnetic anomaly data is the distortion caused by rock body magnetization. There have been several successful studies carried out on the estimation of body magnetization based on the correlation of the Bouguer and magnetic anomaly data. A number of other significant research projects have aimed at identifying the effect of remanent magnetization in magnetic data only. In this paper, we present a three-dimensional method for determining the body magnetization direction by means of comparing the amplitude of the analytic signal with the horizontal gradient of pseudogravity of magnetic anomaly using the correlation coefficient equation. The method has been successfully applied to test cases and is applied here to the aeromagnetic anomalies located in northern Central Turkey. Anomalies displaying individual polarities were separated out and correlated with themselves using the correlation coefficient equation. Counter-clockwise rotations in the range of 3 degrees and 174 degrees were determined. Clockwise rotations were determined in the range of 19 degrees and 153 degrees. The inclination of the body magnetization directions for one of the blocks was estimated to be 68.4 degrees, while all other blocks had a low inclination angle of body magnetization in the range of 32 degrees to 57 degrees. These results may imply that these blocks gained their magnetizations when Turkey was at low latitudes. Possible buried faults are also identified. The results reported here comply with the new mobilistic theory.