$\langle x\rangle$ and $\langle x^2\rangle$ of the pion PDF from lattice QCD with $N_f=2+1+1$ dynamical quark flavors