A formulation that combines the scalar Helmholtz equation with an integral equation radiation boundary condition is developed for electromagnetic scattering from heterogeneous geometries exhibiting one-dimensional periodicity. A finite element discretization is implemented using piecewise-linear interpolation functions. Numerical results are presented for the reflection and transmission coefficients associated with several diffraction gratings.