AbstractUnderstanding soil organic carbon (SOC) response to global change has been hindered by an inability to map SOC at horizon scales relevant to coupled hydrologic and biogeochemical processes. Standard SOC measurements rely on homogenized samples taken from distinct depth intervals. Such sampling prevents an examination of fine‐scale SOC distribution within a soil horizon. Visible near‐infrared hyperspectral imaging (HSI) has been applied to intact monoliths and split cores surfaces to overcome this limitation. However, the roughness of these surfaces can influence HSI spectra by scattering reflected light in different directions posing challenges to fine‐scale SOC mapping. Here, we examine the influence of prescribed surface orientation on reflected spectra, develop a method for correcting topographic effects, and calibrate a partial least squares regression (PLSR) model for SOC prediction. Two empirical models that account for surface slope, aspect, and wavelength and two theoretical models that account for the geometry of the spectrometer were compared using 681 homogenized soil samples from across the United States that were packed into sample wells and presented to the spectrometer at 91 orientations. The empirical approach outperformed the more complex geometric models in correcting spectra taken at non‐flat configurations. Topographically corrected spectra reduced bias and error in SOC predicted by PLSR, particularly at slope angles greater than 30°. Our approach clears the way for investigating the spatial distributions of multiple soil properties on rough intact soil samples.