Plant viral diseases, and especially the ones transmitted by aerial vectors, can cause considerable yield losses. A good knowledge of the distances of spread is key to the understanding of disease dynamics. Exploratory approaches aiming at characterizing the spatiotemporal distribution of diseased plants are often used to get an insight into the distances of spread. A more powerful approach is based on stochastic spatiotemporal modelling in order to estimate the dispersal function of the disease (probability density function describing the probability for an infectious plant to infect a healthy plant at distance d). In this study, we implemented a method for estimating the dispersal function of the sharka disease. Sharka is one of the most serious diseases of stone fruit trees (Prunus sp.). It is caused by Plum pox virus (PPV, genus Potyvirus), transmitted by at least twenty different aphid species in a non persistent manner. Due to the inefficiency of insecticides and the very rare sources of resistance against the virus available in the host species, prophylactic disease control is based on the removal of the diseased trees in the orchards. Thus, a very good knowledge of the dispersal function of sharka is crucial for building epidemiological models and optimizing the strategies of surveillance and control. We adapted the methodology published by Gibson (1997) based on a Markov chain Monte Carlo (MCMC) algorithm in order to estimate sharka dispersal function from the maps of 157 contiguous peach orchards reporting the exact location and the sanitary status (asymptomatic/symptomatic) of each of the trees during six consecutive years. An estimation method based on the Gibbs sampling algorithm was developed taking into account the specificities of the dataset (more than two dates of observation, annual removal of diseased trees). This estimation algorithm was validated on simulated data and was proved to be more powerful and better adapted to large datasets than the