A new scheme of sky pixelization is developed for CMB maps. The scheme is based on the Gauss--Legendre polynomials zeros and allows one to create strict orthogonal expansion of the map. A corresponding code has been implemented and comparison with other methods has been done.