  (Abridged) We have performed a comprehensive multiwavelength analysis of a
sample of 20 starburst galaxies that show the presence of a substantial
population of massive stars. The main aims are the study of the massive star
formation and stellar populations in these galaxies, and the role that
interactions with or between dwarf galaxies and/or low surface companion
objects have in triggering the bursts. We completed new deep optical and \NIR\
broad-band images, as well as the new continuum-subtracted H$\alpha$ maps, of
our sample of Wolf-Rayet galaxies. We analyze the morphology of each system and
its surroundings and quantify the photometric properties of all important
objects. All data were corrected for both extinction and nebular emission using
our spectroscopic data. The age of the most recent star-formation burst is
estimated and compared with the age of the underlying older low-luminosity
population. The \Ha-based star-formation rate, number of O7V equivalent stars,
mass of ionized gas, and mass of the ionizing star cluster are also derived. We
found interaction features in many (15 up to 20) of the analyzed objects, which
were extremely evident in the majority. We checked that the correction for
nebular emission to the broad-band filter fluxes is important in compact
objects and/or with intense nebular emission to obtain realistic colors and
compare with the predictions of evolutionary synthesis models. The estimate of
the age of the most recent star-formation burst is derived consistently. With
respect to the results found in individual objects, we remark the strong \Ha\
emission found in IRAS 08208+2816, UM 420, and SBS 0948+532, the detection of a
double-nucleus in SBS 0926+606A, a possible galactic wind in Tol 9, and one
(two?) nearby dwarf star-forming galaxies surrounding Tol 1457-437.
