The computation of the many beam dynamical electron diffraction amplitudes or high resolution images can only be done numerically by using rather sophisticated computer programs so that the physical insight in the diffraction progress is often lost. Furthermore, it is not likely that in this way the inverse problem can be solved exactly, i.e. to reconstruct the structure of the object from the knowledge of the wavefunction at its exit face, as is needed for a direct method [1]. For this purpose, analytical expressions for the electron wavefunction in real or reciprocal space are much more useful. However, the analytical expressions available at present are relatively poor approximations of the dynamical scattering which are only valid either for thin objects ((weak) phase object approximation, thick phase object approximation, kinematical theory) or when the number of beams is very limited (2 or 3). Both requirements are usually invalid for HREM of crystals. There is a need for an analytical expression of the dynamical electron wavefunction which applies for many beam diffraction in thicker crystals. It is well known that, when a crystal is viewed along a zone axis, i.e. parallel to the atom columns, the high resolution images often show a one-to-one correspondence with the configuration of columns provided the distance between the columns is large enough and the resolution of the instrument is sufficient. This is for instance the case in ordered alloys with a column structure [2,3]. From this, it can be suggested that, for a crystal viewed along a zone axis with sufficient separation between the columns, the wave function at the exit face does mainly depend on the projected structure, i.e. on the type of atom columns. Hence, the classical picture of electrons traversing the crystal as plane-like waves in the directions of the Bragg beams which historically stems from the X-ray diffraction picture, is in fact misleading.