We present lattice QCD calculation of pion valence quark distribution using next-to-next-to-leading order perturbative matching. We use three lattices with fine lattice spacings of $a$ = 0.04, 0.06 and 0.076 fm and two valence pion mass $m_\pi$=300 MeV and 140 MeV to have both discretization and pion mass effect under control. The bare quasi-PDF matrix elements of pion extracted from three-point functions are renormalized in the hybrid scheme and matched to the $\overline{\rm MS}$ scheme. With boost momentum as large as 2.42 GeV, we present reliable lattice QCD prediction of Bjorken-$x$ dependent valence quark distribution of pion in middle $x$ region with controlled systematics.