Distribution of the sum of symmetric lattice vectors with supports on coordinate axes is approximated by multivariate compound Poisson distribution. The characteristic function and Kerstan’s methods are used to obtain local estimates and estimates in total variation.