Volcanic activity is ultimately controlled by processes that take place many kilometers beneath the surface of a planet. The deeper processes are unlikely to reach the surface without some degree of modification at shallower levels. Nevertheless, traces of those deeper processes may still be found when examining the final products at the surface. In this work, it is shown that it is possible to gain insights concerning the integrated contribution of deep structures through the study of the spatial distribution of volcanic vents at the surface. The method here described relies on the systematic use of increasing smoothing factors in Gaussian kernel estimations. The sequences of probability density functions thus generated are equivalent to images obtained with an increasing wavelength, which therefore have the power to penetrate deeper below the surface. Although the resolution of this method is much smaller than the resolution provided by seismic or other geophysical surveys, it has the advantages of ease of implementation, extremely low cost, and remote application. Thus, the reported method has great value as a first-order exploration tool to investigate the deep structure of a planet, and it can make important contributions to our understanding of the volcano-tectonic relationship, not only on Earth, but also across the various bodies of the solar system where volcanic activity has been documented.