Measuring the spatial distribution of heavy metal contaminants is the basis of pollution evaluation and risk control. Considering the cost of soil sampling and analysis, spatial interpolation methods have been widely applied to estimate the heavy metal concentrations at unsampled locations. However, traditional spatial interpolation methods assume the sample sites can be located stochastically on a plane and the spatial association between sample locations is analyzed using Euclidean distances, which may lead to biased conclusions in some circumstances. This study aims to analyze the spatial distribution characteristics of copper and lead contamination in river sediments of Daye using network spatial analysis methods. The results demonstrate that network inverse distance weighted interpolation methods are more accurate than planar interpolation methods. Furthermore, the method named local indicators of network-constrained clusters based on local Moran’ I statistic (ILINCS) is applied to explore the local spatial patterns of copper and lead pollution in river sediments, which is helpful for identifying the contaminated areas and assessing heavy metal pollution of Daye.