The magnetorotational instability (MRI) is a fundamental process of accretion disk physics, but its saturation mechanism remains poorly understood despite considerable theoretical and computational effort. We present a multiple-scales analysis of the non-ideal MRI in the weakly nonlinear regime—that is, when the most unstable MRI mode has a growth rate asymptotically approaching zero from above. Here, we develop our theory in a local, Cartesian channel. Our results confirm the finding by Umurhan et al. that the perturbation amplitude follows a Ginzburg-Landau equation. We further find that the Ginzburg-Landau equation will arise for the local MRI system with shear-periodic boundary conditions, when the effects of ambipolar diffusion are considered. A detailed force balance for the saturated azimuthal velocity and vertical magnetic field demonstrates that, even when diffusive effects are important, the bulk flow saturates via the combined processes of reducing the background shear and rearranging and strengthening the background vertical magnetic field. We directly simulate the Ginzburg-Landau amplitude evolution for our system, and demonstrate the pattern formation our model predicts on long scales of length- and timescales. We compare the weakly nonlinear theory results to a direct numerical simulation of the MRI in a thin-gap Taylor Couette flow.