We propose a new method to determine the spatially or impact-parameter dependent nuclear parton distribution functions (nPDFs) using the double parton scattering (DPS) processes in highenergy heavy-ion (proton-nucleus and nucleus-nucleus) collisions. We derive a simple generic DPS formula in nuclear collisions by accommodating both the nuclear collision geometry and the spatially dependent nuclear modification effect, under the assumption that the impact-parameter dependence of nPDFs is only related to the nuclear thickness function. While the geometric effect is widely adopted, the impact of the spatially dependent nuclear modification on DPS cross sections has been overlooked so far, which can, however, be significant when the initial nuclear modification is large. In turn, the DPS cross sections in heavy-ion collisions can provide useful information on the spatial dependence of nPDFs. They can be, in general, obtained in minimum-bias nuclear collisions, featuring the virtue of independence of Glauber modeling.1 To the best of our knowledge, none of the existing DPS phenomenological applications in heavy-ion collisions [89][90][91][92][96][97][98][99][100][101][102][103][104] considers such an effect. 2 There are also a few attempts to obtain the spatial form of the nPDFs based on phenomenological models (see, e.g., ).