Abstract-The application of the finite-difference time-domain method with rectangular periodic boundary conditions to the analysis of a hexagonal photonic crystal results in a folded bandgap diagram. The aim of this paper is to introduce a new unfolding method, which allows unambiguously determining the position of the modes in a wavevector space by taking the advantage of the fast Fourier transform of modal field distributions. Unlike alternative solutions, it does not require any modifications of the FDTD method and is based solely on the post-processing of the simulation results. The proposed method can be applied to any non-rectangular lattice types, such as hexagonal, face-centered cubic or body-centered cubic.