The pair correlation function (PCF) has proven an effective tool for analyzing many physical systems due to its simplicity and its applicability for simulated and experimental data. However, as an averaged quantity, the PCF can fail to capture subtle structural differences in particle arrangements, even when those differences can have a major impact on system properties. Here, we use Voronoi topology to introduce a discrete version of the PCF that highlights local interparticle topological configurations. The advantages of the Voronoi PCF are demonstrated in several examples including crystalline, hyperuniform, and active systems showing clustering and giant number fluctuations.