We present the results of a Lattice QCD computation of pion generalized parton distribution (GPD), employing perturbative matching up to next-to-next-to-leading order (NNLO). The computations are based on an ensemble of $N_f=2+1$ highly improved staggered quarks (HISQ) with a pion mass of 300 MeV and a lattice spacing of 0.04 fm. Centered on the zero-skewness limit, we utilize a recently proposed Lorentz-invariant definition of GPD, which is derived from Lorentz-invariant amplitudes. We analyze and compare these amplitudes in both Breit and non-Breit kinematic frames at comparable momentum transfers, validating their frame-independent nature. To obtain light-cone GPD, we integrate hybrid scheme renormalization with the large momentum effective theory (LaMET). Moreover, we determine the first three iso-vector generalized form factors (GFFs) of the pion using the ratio scheme renormalization and leading-twist factorization, achieving NNLO accuracy.