In this article, the contrast imaging problem in nuclear magnetic resonance is modeled as a Mayer problem in optimal control whose solutions can be parameterized using Pontryagin's maximum principle and analyzed using geometric optimal control. In particular, the optimal problem can be mainly reduced to the analysis of the Hamiltonian dynamics describing the singular trajectories, and encoding their optimality status.