A method of constructing three-dimensional cubature formulae of high degree for three-dimensional regions is presented. The formulae possess full polyhedral symmetry. This means that every knot of the formula belongs to an orbit invariant with respect to the symmetries of the polyhedron and that all points of an orbit are knots of the formula with the same weight.