In order to estimate the annual dose that the public receive from natural radioactivity, the identification of the potential risk areas is required which, in turn, necessitates understanding the relationship between the spatial distribution of natural radioactivity and the geogenic risk factors (e.g., rock types, presence of dikes, faults, physical conditions of soil, etc.). A detailed spatial analysis of outdoor ambient gamma dose equivalent rate was performed in the western side of Velence Hills, the largest outcropped granitic area in Hungary. In order to assess the role of local geology in the spatial distribution of gamma dose rates, field measurements were carried out at ground level at 300 sites along a 250 m x 250 m regular grid in a total surface of 19.8 km(2). Digital image processing methods were applied to identify anomalies, heterogeneities and spatial patterns in the measured gamma dose rates, including local maxima and minima determination, digital cross sections, gradient magnitude and gradient direction, second derivative profile curvature, local variability, lineament density, 2D auto correlation and directional variogram analyses. Statistical inference shows that different gamma dose rate levels are associated with the geological formations, with the highest level on the Carboniferous granite including outlying values. Moreover, digital image processing reveales that linear gamma dose rate spatial features are parallel to the SW-NE dike system and to the NW-SE main fractures. The results of this study underline the importance of understanding the role of geogenic risk factors influencing the ambient gamma dose equivalent rate received by public. The study also demonstrates the power of the image processing techniques for the identification of spatial pattern in field-measured geogenic radiation.