An enhanced three-dimensional Finite-Difference Time-Domain (FDTD) technique is applied to the characterization of the fundamental and higher order modes of single and coupled image guide structures. The basic FDTD algorithm is modified to include multidimensional changes of permittivity and an irregular mesh is used for computational efficiency. In this paper FDTD calculations are compared with experimental and theoretical results found in the literature and good agreement is foun