The computation of electromagnetic fields in planar, multilayered, uniaxial media is discussed. In particular, a complete set of electric and magnetic dyadic Green functions (DGFs) is derived for such media in a computationally efficient and convenient to use form, based on a spectral domain transmission line equivalent network along the axis normal to the stratification. Efficient methods for the computation of the Sommerfeld integrals that arise in the DGFs are also presented.