In a recent study, using the distribution of galaxies in the north galactic pole of the SDSS DR7 region enclosed in a 500 h(-1) Mpc box, we carried out our ELUCID simulation (ELUCID III). Here, we light the dark matter halos and subhalos in the reconstructed region in the simulation with galaxies in the SDSS observations using a novel neighborhood abundance matching method. Before we make use of the galaxy-subhalo connections established in the ELUCID simulation to evaluate galaxy formation models, we set out to explore the reliability of such a link. For this purpose, we focus on the following few aspects of galaxies: (1) the central-subhalo luminosity and mass relations, (2) the satellite fraction of galaxies, (3) the conditional luminosity function (CLF) and conditional stellar mass function (CSMF) of galaxies, and (4) the cross-correlation functions between galaxies and dark matter particles, most of which are measured separately for all, red, and blue galaxy populations. We find that our neighborhood abundance matching method accurately reproduces the central-subhalo relations, satellite fraction, and the CLFs, CSMFs, and biases of galaxies. These features ensure that galaxy-subhalo connections thus established will be very useful in constraining galaxy formation processes. We provide some suggestions for the three levels of using the galaxy-subhalo pairs for galaxy formation constraints. The galaxy-subhalo links and the subhalo merger trees in the SDSS DR7 region extracted from our ELUCID simulation are available upon request.