The unprecedented spatiotemporal resolution of the sixty-four dual electron and ion spectrometers comprising the Fast Plasma Investigation (FPI) onboard the four Magnetospheric Multiscale (MMS) spacecraft enables us to compute terms of the Vlasov equation and thus compare MMS measurements directly to kinetic theory. Here we discuss our techniques for determining spatial and velocity-space gradients of the ion and electron distribution function from the skymaps provided by FPI. We present initial results validating and comparing the contributions from ∂f/∂t, v·∇f, and a·∇vf for a variety of kinetic plasma contexts including both reconnection and non-reconnection events, such as electron diffusion regions (EDRs), magnetic holes, and thin electron-scale current sheets. The ability to resolve gradients of the distribution function motivates comparison of MMS observations to predictions from gyro-kinetic theory and particle-in-cell (PIC) simulations as an approach for determining which physical mechanism is responsible for generating ion and electron crescent distributions observed in association with EDRs and thin boundary layers.