We consider the interpretation of planetary anomalies of the gravitational fields of Mars and Moon in relation to their internal heterogeneous structure. As the source information we used the models of the gravitation fields of Mars and Moon up to the 20-th order and degree as well as parameters of three-layer model (crust, mantle, core) and planetary parameters. The models for three-dimensional density distribution for mars and Moon are constructed. The interpretation of the maps for horizontal heterogeneities of density at depths of 50, 100, and 1700 km for Mars and 60, 100, and 1400 km for Moon is presented.