We analyse the density field of galaxies observed by the Sloan Digital Sky Survey (SDSS)- III Baryon Oscillation Spectroscopic Survey (BOSS) included in the SDSS Data Release Nine (DR9). DR9 includes spectroscopic redshifts for over 400 000 galaxies spread over a footprint of 3275 deg2. We identify, characterize and mitigate the impact of sources of systematic uncertainty on large-scale clustering measurements, both for angular moments of the redshift-space correlation function, ξ0(s), and the spherically averaged power spectrum, P(k), in order to ensure that robust cosmological constraints will be obtained from these data. A correlation between the projected density of stars and the higher redshift (0.43 < z < 0.7) galaxy sample (the approximately constant stellar mass threshold ‘CMASS’ sample) due to imaging systematics imparts a systematic error that is larger than the statistical error of the clustering measurements at scales s > 120 h−1 Mpc or k < 0.01 h Mpc−1. We find that these errors can be ameliorated by weighting galaxies based on their surface brightness and the local stellar density. The clustering of CMASS galaxies found in the Northern and Southern Galactic footprints of the survey generally agrees to within 2σ.We use mock galaxy catalogues that simulate the CMASS selection function to determine that randomly selecting galaxy redshifts in order to simulate the radial selection function of a random sample imparts the least systematic error on ξ0(s) measurements and that this systematic error is negligible for the spherically averaged correlation function, ξ0.We find a peak in ξ0 at s ∼ 200 h−1 Mpc, with a corresponding feature with period ∼0.03 h Mpc−1 in P(k), and find features at least as strong in 4.8 per cent of the mock galaxy catalogues, concluding this feature is likely to be a consequence of cosmic variance. The methods we recommend for the calculation of clustering measurements using the CMASS sample are adopted in companion papers that locate the position of the baryon acoustic oscillation feature, constrain cosmological models using the full shape of ξ0 and measure the rate of structure growth.