We propose a stellar locus outlier (SLOT) method to determine the binary fraction of main-sequence stars statistically. The method is sensitive to neither the period nor mass ratio distributions of binaries and is able to provide model-free estimates of binary fraction for large numbers of stars of different populations in large survey volumes. We have applied the SLOT method to two samples of stars from the Sloan Digital Sky Survey (SDSS) Stripe 82, constructed by combining the recalibrated SDSS photometric data with the spectroscopic information from the SDSS and LAMOST surveys. For the SDSS spectroscopic sample, we find an average binary fraction for field FGK stars of 41% ± 2%. The fractions decrease toward late spectral types and are 44% ± 5%, 43% ± 3%, 35% ± 5%, and 28% ± 6% for stars with g − i colors in the range 0.3-0.6 mag, 0.6-0.9 mag, 0.9-1.2 mag, and 1.2-1.6 mag, respectively. A modest metallicity dependence is also found. The fraction decreases with increasing metallicity. For stars with [Fe/H] between −0.5 and 0.0 dex, −1.0 and −0.5 dex, −1.5 and −1.0 dex, and −2.0 and −1.5 dex, the inferred binary fractions are 37% ± 3%, 39% ± 3%, 50% ± 9%, and 53% ± 20%, respectively. We have further divided the sample into stars from the thin disk, the thick disk, the transition zone between them, and the halo. The results suggest that the Galactic thin and thick disks have comparable binary fractions, whereas the Galactic halo contains a significantly larger fraction of binaries. Applying the method to the LAMOST spectroscopic sample yields consistent results. Finally, other potential applications and future work with the method are discussed.