From a generic transformation, a $(3+1)$-dimensional nonlinear Schrödinger equation with spatially modulated nonlinearity is studied and exact spatiotemporal soliton cluster solutions are derived. When the azimuthal parameter $m = 0$, Gaussian solitons are constructed. For the modulation depth $q = 1$ and the azimuthal parameter $m \neq 0$, vortex solitons are obtained. In the limit $q = 0$ and $m \neq 0$, and for some specific values of the parameters $(m, l)$, multipole solitons are presented.
