Pulsed-field-gradient nuclear magnetic resonance (NMR) combined with magic echo decoupling is applied to study anisotropic diffusion in samples with strong static dipolar spin interactions. The approach, due to its moderate demands on the NMR hardware, can be implemented on standard commercial equipment for routine diffusion studies of liquid crystals. Using a microimaging probe, measurement of diffusion in arbitrary spatial direction is possible. Hence, the principal components of the diffusion tensor are directly obtained. Anisotropic diffusion is investigated in the thermotropic mesophases of a homologous series of nOCB liquid crystals and an analogous compound with hydroxyl groups. The geometric average diffusion coefficient changes continuously at the isotropic-nematic phase transition. Experimental data are described in terms of the molecular translation models in the nematic phase and for the second-order nematic-smectic A phase transition. The diffusion anisotropy is higher for the sample with terminal hydroxyl groups suggesting significant molecular association via hydrogen bonding.