Strong field photoemission and electron recollision provide a viable route to
extract electronic and nuclear dynamics from molecular targets with attosecond
temporal resolution. However, since an {\em ab-initio} treatment of even the
simplest diatomic systems is beyond today's capabilities approximate
qualitative descriptions are warranted. In this paper, we develop such a
theoretical approach to model the photoelectrons resulting from intense
laser-molecule interaction. We present a general theory for symmetric diatomic
molecules in the single active electron approximation that, amongst other
capabilities, allows adjusting both the internuclear separation and molecular
potential in a direct and simple way. More importantly we derive an analytic
approximate solution of the time dependent Schr\"odinger equation (TDSE), based
on a generalized strong field approximation (SFA) version. Using that approach
we obtain expressions for electrons emitted transition amplitudes from two
different molecular centres, and accelerated then in the strong laser field.
One innovative aspect of our theory is the fact that the dipole matrix elements
are free from non-physical gauge and coordinate system dependent terms -- this
is achieved by adapting the coordinate system, in which SFA is performed, to
the centre from which the corresponding part of the time dependent wave
function originates. Our analytic results agree very well with the numerical
solution of the full three-dimensional TDSE for the H$_2^+$ molecule. Moreover,
the theoretical model was applied to describe laser-induced electron
diffraction (LIED) measurements of O$_2^+$ molecules, obtained at ICFO, and
reproduces the main features of the experiment very well. Our approach can be
extended in a natural way to more complex molecules and multi-electron systems.
The ability to directly follow and time resolve the rearrangement of the
nuclei within molecules is a frontier of science that requires atomic spatial
and few-femtosecond temporal resolutions. While laser induced electron
diffraction can meet these requirements, it was recently concluded that
molecules with particular orbital symmetries (such as {\pi}g) cannot be imaged
using purely backscattering electron wave packets without molecular alignment.
Here, we demonstrate, in direct contradiction to these findings, that the
orientation and shape of molecular orbitals presents no impediment for
retrieving molecular structure with adequate sampling of the momentum transfer
space. We overcome previous issues by showcasing retrieval of the structure of
randomly oriented O2 and C2H2 molecules, with {\pi}g and {\pi}u symmetries,
respectively, and where their ionisation probabilities do not maximise along
their molecular axes. While this removes a serious bottleneck for laser induced
diffraction imaging, we find unexpectedly strong back scattering contributions
from low-Z atoms.