Xylella fastidiosa is one of the most destructive plant pathogenic bacteria worldwide, affecting more than 500 plant species. In Apulia region (southeastern Italy), X. fastidiosa subsp. pauca (Xfp) is responsible for a severe disease, the olive quick decline syndrome (OQDS), spreading epidemically and with dramatic impact on the agriculture, the landscape, the tourism, and the cultural heritage of this region. An early detection of the infected plants would hinder the rapid spread of the disease. The main objective of this paper was to define a geostatistical approach of data fusion, which combines remote (radiometric), and proximal (geophysical) sensor data and visual inspections with plant diagnostic tests, to provide probabilistic maps of Xfp infection risk. The study site was an olive grove located at Oria (province of Brindisi, Italy), where at the time of monitoring (September 2017) only few plants showed initial symptoms of the disease. The measurements included: 1) acquisitions of reflected