In large-momentum effective theory (LaMET), calculating parton physics starts from calculating coordinate-space-z correlation functions in a hadron of momentum 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 scheme, without introducing extra non-perturbative effects at large z. We analyze the effect of 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 propertie...