A number of quantum Hall isospin ferromagnetic (QHIFM) states have been predicted in the "relativistic" zero Landau level (LL) of graphene monolayer. These states, especially the states at LL filling factor = 0 of charge-neutral graphene, have been extensively explored in experiment. To date, identification of these high-field broken-symmetry states has mostly relied on macroscopic transport techniques. Here, we study splitting of the zero LL of graphene at partial filling and demonstrate a direct approach by imaging the QHIFM states at atomic scale with a scanning tunneling microscope. At half filling of the zero LL ( = 0), the system is in a spin unpolarized state and we observe a linear magnetic-fieldscaling of valley splitting. Simultaneously, the spin degeneracy in the two valleys is also lifted by the magnetic fields. When the Fermi level lies inside the spinpolarized states (at = 1 or -1), the spin splitting is dramatically enhanced because of the strong many-body effects. At = 0, we direct image the wavefunctions of the QHIFM states at atomic scale and observe an interaction-driven density wave featuring a Kekulé distortion, which is responsible for the large gap at charge neutrality point in high magnetic fields.