We demonstrate that the spatial profiles of both propagating and evanescent Bloch-modes in a periodic structure can be extracted from a single measurement of electric field at the specified optical wavelength. We develop a systematic extraction procedure by extending the concepts of high-resolution spectral methods previously developed for temporal data series to take into account the symmetry properties of Bloch-modes simultaneously at all spatial locations. We demonstrate the application of our method to a photonic crystal waveguide interface and confirm its robustness in the presence of noise. Periodically modulated optical waveguides offer new possibilities for controlling the propagation of light. Resonant scattering from periodic modulations can be used to tailor the dispersion, enabling in particular a dramatic modification of the group velocity and realization of slow-light propagation. Such fundamental effects can be directly visualized in experiment with near-field measurements, which can be used to recover the amplitude, phase, and polarization of the electric field at all spatial locations in the plane of the waveguide [1]. This information can then be used to extract the dispersion characteristics of the guided modes.A commonly used approach to the dispersion extraction is through the spatial Fourier-transform (SFT) of the field profiles, since peaks in the Fourier spectra correspond to the wavenumbers of guided modes [2,3,4]. However, there exists a fundamental limitation on results obtained with SFT: ∆k ≥ 2π/L, where ∆k is the resolution of the wavenumber, and L is the structure's length. Therefore, accurate dispersion results can only be obtained for long waveguides, extending over many periods of the underlying photonic structure. Another limitation of the SFT method is that it cannot provide information on the dispersion of evanescent waves, which may play an important role close to the structure boundaries or interfaces between different waveguides. For example, evanescent waves enable efficient excitation of slow-light waves without a transition region [5].Alternative methods for dispersion extraction have been developed to overcome the shortcoming of the SFT method. It was shown that an interference of two counter-propagating modes can be used to extract their wavenumbers [6], however this technique is not applicable under the presence of multiple propagating modes or evanescent waves. Recently, it was demonstrated that dispersion extraction in multi-mode waveguides with in principle unbounded resolution is possible even for short waveguide sections [7,8], using approaches based on an adaptation of high-resolution spectral methods previously developed for the analysis of temporal dynamics [9,10]. In this work, we introduce an important generalization of such methods taking into account the spatial symmetry properties of modes in periodic waveguides. We show that beyond the dispersion relations, it is possible to extract the spatial profiles of all guided modes. Our method is applicable to an...