Brain mapping data have been modelled as Gaussian random fields, and `activation', or local increases in mean, are detected by local maxima of a random field of test statistics derived from this data. Accurate P-values for local maxima are readily available for isotropic data based on the expected Euler characteristic of the excursion set of the test statistic random field. In this paper we give a simple method for dealing with non-isotropic data. Our approach has connections to the Sampson & Guttorp (1992) model for non-isotropy in which their exists an unknown mapping of the support of the data to a space in which the random fields are isotropic. Heuristic justification for our approach comes from the Nash Embedding Theorem. Formally we show that our method gives consistent unbiased estimators for the true P-values based on new results of Taylor & Adler (2003) for random fields on manifolds which replace the Euclidean metric by the variogram. The results are used to detect gender differences in the cortical thickness of the brain, and to detect regions of the brain involved in sentence comprehension measured by fMRI.