This work presents the reconstruction of the pressure head map of a water distribution system (WDS). This approach relies on historical data collected from a reduced number of sensors placed at some nodes of the WDS. Thus, a Gaussian regression process is then applied to estimate the pressure head at those nodes without a sensor, which allows the reconstruction of the pressure map for the whole network. Then, for leak diagnosis purposes, a dataset of pressure head maps of the WDN is created considering leaky scenarios, and a correlation method is applied to estimate the leak location. Then, for clarity, the Hanoi network is used to evaluate the performance of this leak diagnosis strategy in a simulation environment, assuming the availability of only three sensors. The results showed the potential for pressure head map reconstruction and leak localization.