In large-momentum effective theory (LaMET), calculating parton physics starts from calculating coordinate-space-z correlation functions h?(z, a, Pz) in a hadron of momentum Pz 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 MS scheme, without introducing extra non-perturbative effects at large z. We analyze the effect of O(AQCD) 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 Pz, respectively. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP3.