This article describes a new method for 3D QEEG tomography in the frequency domain. A variant of Statistical Parametric Mapping is presented for source log spectra. Sources are estimated by means of a Discrete Spline EEG inverse solution known as Variable Resolution Electromagnetic Tomography (VARETA). Anatomical constraints are incorporated by the use of the Montreal Neurological Institute (MNI) probabilistic brain atlas. Efficient methods are developed for frequency domain VARETA in order to estimate the source spectra for the set of 10(3)-10(5) voxels that comprise an EEG/MEG inverse solution. High resolution source Z spectra are then defined with respect to the age dependent mean and standard deviations of each voxel, which are summarized as regression equations calculated from the Cuban EEG normative database. The statistical issues involved are addressed by the use of extreme value statistics. Examples are shown that illustrate the potential clinical utility of the methods herein developed.