The STAR Collaboration introduces a novel method, collective flow-assisted nuclear shape imaging, to image the global shape of atomic nuclei by colliding them at ultrarelativistic speeds and analyzing the collective response of outgoing debris. This technique captures a collision-specific snapshot of the spatial matter distribution in the nuclei, which leaves imprints on the particle momentum distribution patterns observed in detectors. The method is benchmarked using ground-state Uranium-238 nuclei, known for its elongated, axial-symmetric shape. The findings confirm an overall deformation consistent with prior low-energy experiments but also indicate a small deviation from axial symmetry in the nuclear ground state. This approach provides a new way to image nuclei, especially those with uncertain shape characteristics, and refines initial conditions in high-energy nuclear collisions. It addresses the issue of nuclear structure evolution across various energy scales, highlighting the need to examine nuclear phenomena across different energy levels. The technique is applicable to a wide range of applications, including odd-mass nuclides, octupole and hexadecapole deformations, dynamic deformations in "soft" nuclei, and 0νββ decay.The STAR Collaboration introduces a novel method, collective flow-assisted nuclear shape imaging, to image the global shape of atomic nuclei by colliding them at ultrarelativistic speeds and analyzing the collective response of outgoing debris. This technique captures a collision-specific snapshot of the spatial matter distribution in the nuclei, which leaves imprints on the particle momentum distribution patterns observed in detectors. The method is benchmarked using ground-state Uranium-238 nuclei, known for its elongated, axial-symmetric shape. The findings confirm an overall deformation consistent with prior low-energy experiments but also indicate a small deviation from axial symmetry in the nuclear ground state. This approach provides a new way to image nuclei, especially those with uncertain shape characteristics, and refines initial conditions in high-energy nuclear collisions. It addresses the issue of nuclear structure evolution across various energy scales, highlighting the need to examine nuclear phenomena across different energy levels. The technique is applicable to a wide range of applications, including odd-mass nuclides, octupole and hexadecapole deformations, dynamic deformations in "soft" nuclei, and 0νββ decay.