In this paper we develop the theory of clusterization of peaks in a Gaussian random field. We have obtained new mathematical results from this theory and the theory of percolation and have proposed a topological method of analysis of sky maps based on these results. We have simulated $10^otimes10^o$ sky maps of the cosmic microwave background anisotropy expected from different cosmological models with $0.5^o-1^o$ resolution in order to demonstrate how this method can be used for detection of non-Gaussian noise in the maps and detection of the Doppler-peak in the spectrum of perturbation of $gD T/T$.