The need for operational models describing the friction factor f in streams remains undisputed given its utility across a plethora of hydrological and hydraulic applications concerned with shallow inertial flows. For small-scale roughness elements uniformly covering the wetted parameter of a wide channel, the Darcy-Weisbach f = 8(u*/Ub)2 is widely used at very high Reynolds numbers, where u* is friction velocity related to the surface kinematic stress, Ub = Q/A is bulk velocity, Q is flow rate, and A is cross-sectional area orthogonal to the flow direction. In natural streams, the presence of vegetation introduces additional complications to quantifying f, the subject of the present work. Turbulent flow through vegetation are characterized by a number of coherent vortical structures: (i) von Karman vortex streets in the lower layers of vegetated canopies, (ii) Kelvin-Helmholtz as well as attached eddies near the vegetation top, and (iii) attached eddies well above the vegetated layer. These vortical structures govern the canonical mixing lengths for momentum transfer and their influence on f is to be derived. The main novelty is that the friction factor of vegetated flow can be expressed as fv = 4Cd(Uv/Ub)2 where Uv is the spatially averaged velocity within the canopy volume, and Cd is a local drag coefficient per unit frontal area derived to include the aforemontioned layer-wise effects of vortical structures within and above the canopy along with key vegetation properties. The proposed expression is compared with a number of empirical relations derived for vegetation under emergent and submerged conditions as well as numerous data sets covering a wide range of canopy morphology, densities, and rigidity. It is envisaged that the proposed formulation be imminently employed in eco-hydraulics where the interaction between flow and vegetation is being sought.