We establish a self-dual version of Brenier's polar decomposition for non-degenerate vector fieldswhere S is a self-dual ("almost idempotent") measure preserving point transformation on Ω, and H : R N × R N → R is a globally Lipschitz anti-symmetric convex-concave Hamiltonian. Just like Brenier's polar decomposition, our representation unifies several seemingly unrelated classical results.