We demonstrate a new method of extracting parton distributions from lattice calculations. The starting idea is to treat the generic equal-time matrix element ${\cal M} (Pz_3, z_3^2)$ as a function of the Ioffe time $ u = Pz_3$ and the distance $z_3$. The next step is to divide ${\cal M} (Pz_3, z_3^2)$ by the rest-frame density ${\cal M} (0, z_3^2)$. Our lattice calculation shows a linear exponential $z_3$-dependence in the rest-frame function, expected from the $Z(z_3^2)$ factor generated by the gauge link. Still, we observe that the ratio ${\cal M} (Pz_3 , z_3^2)/{\cal M} (0, z_3^2)$ has a Gaussian-type behavior with respect to $z_3$ for 6 values of $P$ used in the calculation. This means that $Z(z_3^2)$ factor was canceled in the ratio. When plotted as a function of $ u$ and $z_3$, the data are very close to $z_3$-independent functions. This phenomenon corresponds to factorization of the $x$- and $k_\perp$-dependence for the TMD ${\cal F} (x, k_\perp^2)$. For small $z_3 \leq 4a$, the residual $z_3$-dependence is explained by perturbative evolution, with $lpha_s/\pi =0.1$.