© 2024 authors. Published by the American Physical Society.We combine the Dyson-Schwinger/Bethe-Salpeter equations framework with modern numerical reconstruction methods to derive the three-dimensional and transverse two-dimensional charge distribution of an array of ground-state pseudoscalar and vector mesons from their elastic electromagnetic form factor in the low-momentum region. The charge radii obtained by averaging over the reconstructed charge distributions have been checked to be consistent with those calculated from the slope of the elastic electromagnetic form factor at zero transferred momentum. The capability of the reconstruction procedure for capturing a reliable low-distance charge distribution is discussed and argued to work down to distances of around 0.1 fm, such that it might be potentially applied to extract, e.g., mass densities from gravitational form factors.