  We define in the framework of light-cone collinear factorization method, the
chiral-odd generalized parton distributions (GPDs) of a pseudoscalar hadron
(such as the $\pi^0$) up to twist 6. For that, we introduce the relevant matrix
elements for 2-parton non-local operators, as well as matrix elements for
3-parton non-local correlators. Their detailed parametrization is fixed based
on parity, charge conjugation and time reversal invariance. This leads to the
introduction of 28 real GPDs, which are subject to constraints coming from the
QCD equations of motion.
