The Contour Integral Method (CIM) is a numerically efficient modeling technique for planar or infinitely extended two-dimensional (2-D) structures. In the optical regime, the CIM has already been adapted and applied for the modeling of TM0z-mode scattering in photonic crystals. In this work the dual case of TE0z-mode scattering is addressed. Making use of the duality principle, expressions for the behavior of the TE0z-mode can be derived from the system matrices associated with the TE0z-mode. This allows to reuse use formulas and program code written for the TE0z-mode with minimal adjustments. The results are validated by comparison to full-wave simulations.