Reconstruction of current density in electron beams is complicated by distortion of phase space which can require very fine discretization of the beam into trajectories. An efficient discretization of phase space is exploited, using conservation of charge and current in hypertriangle patches, to reconstruct the current density by fitting Gaussians through the distorted hypertriangles. Advantages and limitations are discussed.