SUMMARYThe spectral domain approach for the solution of integral equations is especially attractive for planar circuits and antennas in multilayered environments characterized by analytically derived Green's functions. However, for circuits with large lateral dimensions the spectral integrands are strongly oscillatory and accurate evaluation of the integrals is extremely difficult and often time-consuming. In this paper, we propose higher-order Legendre quadrature rules optimized for the integration of the oscillating integrands using the Filon principle. The Legendre-Filon quadrature weights are computed for the specific oscillations by solving a small linear system of equations. The method is applied within a spectral domain solution approach evaluating the spectral integrals in cartesian wavenumbers and using asymptotic extraction techniques as well as complex integration path deformation for improved convergence. Since the asymptotic extraction causes a singularity in the origin of the wavenumber domain that is not annihilated by cartesian wavenumber integration, a special procedure for the treatment of this singularity is also introduced. The overall method needs very few sampling points for the numerical evaluation of the spectral integrals but is still applicable to very large structures as compared to the wavelength.