The complete iron vibrational spectrum of the five-coordinate high-spin complex [Fe(OEP)(2-MeHIm)], where OEP = octaethylporphyrinato and 2-MeHIm = 2-methylimidazole, has been obtained by oriented single-crystal nuclear resonance vibrational spectroscopy (NRVS) data. Measurements have been made in three orthogonal directions, which provides quantitative information for all iron motion. These experimental data, buttressed by density functional theory (DFT) calculations, have been used to define the effects of the axial ligand orientation. Although the axial imidazole removes the degeneracy in the in-plane vibrations, the imidazole orientation does not appear to control the direction of the in-plane iron motion. This is in contrast to the effect of the imidazolate ligand, as defined by DFT calculations, which does have substantial effects on the direction of the in-plane iron motion. The axial NO ligand has been found to have the strongest orientational effect (Angew. Chem., Int. Ed., 2010, 49, 4400). Thus the strength of the directional properties are in the order NO > imidazolate > imidazole, consistent with the varying strength of the Fe-ligand bond.