  In large-momentum effective theory (LaMET), calculating parton physics starts
from calculating coordinate-space-$z$ correlation functions $\tilde h(z,
a,P^z)$ in a hadron of momentum $P^z$ in lattice QCD. Such correlation
functions involve both linear and logarithmic divergences in lattice spacing
$a$, and thus need to be properly renormalized. We introduce a hybrid
renormalization procedure to match these lattice correlations to those in the
continuum $\overline{\rm MS}$ scheme, without introducing extra
non-perturbative effects at large $z$. We analyze the effect of ${\cal
O}(\Lambda_{\rm QCD})$ ambiguity in the Wilson line self-energy subtraction
involved in this hybrid scheme. To obtain the momentum-space distributions, we
recommend to extrapolate the lattice data to the asymptotic $z$-region using
the generic properties of the coordinate space correlations at moderate and
large $P^z$, respectively.
