We compute the pion quark Generalized Parton Distribution H and quark Double Distributions in a coupled Bethe-Salpeter and Dyson-Schwinger approach in terms of quarks flavors or isospin states. We use simple analytic expressions inspired by the numerical resolution of Dyson-Schwinger and Bethe-Salpeter equations. We obtain an analytic expression for the pion Generalized Parton Distribution at a low scale. Our model fulfils most of the required symmetry properties and compares very well to experimental pion form factor or parton distribution function (PDF) experimental data. In addition, we have highlighted gluons contributions which play a significant role when computing the PDF and which were neglected before.