Several binary systems that contain a massive star have been detected in both the radio band and at very high energies. In the dense stellar photon field of these sources, gamma-ray absorption and pair creation are expected to occur, and the radiation from these pairs may contribute significantly to the observed radio emission. We aim at going deeper in studying the properties, and in particular the morphology, of the pair radio emission in gamma-ray binaries. We applied a Monte Carlo code that computes the creation location, the spatial trajectory and the energy evolution of the pairs produced in the binary system and its surroundings. The radio emission produced by these pairs, with its spectral, variability and spatial characteristics, was calculated as it would be seen from a certain direction. A generic case was studied first, and then the specific case of LS 5039 was also considered. We found that, confirming previous results, the secondary radio emission should appear as an extended radio structure of a few milliarcseconds size. This radiation would be relatively hard, with fluxes of up to $\sim$ 10 mJy. Modulation is expected, depending on the gamma-ray production luminosity, system eccentricity, and wind ionization fraction, and to a lesser extent on the magnetic-field structure. In gamma-ray binaries in general, the pairs created due to photon-photon interactions can contribute significantly to the core, and generate an extended structure. In the case of LS 5039, the secondary radio emission is likely to be a significant fraction of the detected core flux, with a marginal extension.