The reversal of magnetic bubble helicity through topologically trivial transient states provides an additional degree of freedom that promises the development of multidimensional magnetic memories. A key requirement for this concept is the stabilization of bubble states at ambient conditions on application-compatible substrates. In the present work, we demonstrate a stabilization routine for remanent bubble states in high perpendicular magnetic anisotropy ${\text{[Co(0.44 nm)/Pt(0.7 nm)]}}_{X},$ $X=48$, 100, and 150 multilayers on $\mathrm{Si}/{\mathrm{SiO}}_{2}$ substrates by exploring the effect of external magnetic fields $({H}_{\mathrm{m}})$ of different strength and angles $(\ensuremath{\theta})$ with respect to the film surface normal. By systematic variation of these two parameters, we demonstrate that remanent bubble density and mean bubble diameter can be carefully tuned and optimized for each sample. Our protocol based on magnetometry only reveals the densest remanent bubble states at ${H}_{\mathrm{m}}=0.87{H}_{\mathrm{s}}$ (${H}_{\mathrm{s}}$ is the magnetic saturation field) and $\ensuremath{\theta}={60}^{\ensuremath{\circ}}$--${75}^{\ensuremath{\circ}}$ for all $X$ with a maximum of 3700 domains/100 $\ensuremath{\mu}{\mathrm{m}}^{2}$ for the $X=48$ sample. The experimental observations are supported by micromagnetic simulations, taking into account the nanoscale lateral grain structure of multilayers synthesized by magnetron sputter deposition, and thus helping to understand the different densities of the bubble states found in these systems.