A new method is proposed for the study of crustal azimuthal anisotropy, and the lateral variation of other crustal aerogeophysical parameters, using a quantitative approach. The processing of a large, newly acquired, high-resolution set of aerogeophysical data is considered using the existing Radon transform and the geostatistical analysis approaches. The data set includes an area of the Guyana shield, in northern Brazil, which was not included in previous surveys. The area is covered by dense rain forest vegetation and thick soil layers. Parameterization was performed considering the possible anisotropic character of the geophysical 2D data. Application of the newly proposed geostatistical data processing yielded high-resolution images of the lateral variation of quantitative geophysical parameters, which indicate good correlation with previously determined seismic anisotropy in the area. Average anisotropy as measured by the Radon transform and variogram analysis is scale dependent. At scales greater than , the results seem to match those of the previous seismic studies. Images of the derived quantitative parameters from magnetic and radiometric data in the shield area indicate sharp changes that follow the known geologic changes observable in the survey area. Observed large-scale anisotropy in the area seems to be autocorrelated with two different geochronological provinces of the Amazonian craton. Application of the new technique to two other lower resolution aeromagnetic surveys located on the Amazon basin has produced two different results, including one successful data treatment of the aeromagnetic survey. The results obtained for this area indicate that the sedimentary cover seems to play almost no role in the observed magnetic anomaly anisotropy in the middle Amazon basin area. Most of the observed magnetic anomalies appear to be controlled by the continuing geologic features of the cratonic basement.