  Spatial spread of infectious diseases among populations via the mobility of
humans is highly stochastic and heterogeneous. Accurate forecast/mining of the
spread process is often hard to be achieved by using statistical or mechanical
models. Here we propose a new reverse problem, which aims to identify the
stochastically spatial spread process itself from observable information
regarding the arrival history of infectious cases in each subpopulation. We
solved the problem by developing an efficient optimization algorithm based on
dynamical programming, which comprises three procedures: i, anatomizing the
whole spread process among all subpopulations into disjoint componential
patches; ii, inferring the most probable invasion pathways underlying each
patch via maximum likelihood estimation; iii, recovering the whole process by
assembling the invasion pathways in each patch iteratively, without burdens in
parameter calibrations and computer simulations. Based on the entropy theory,
we introduced an identifiability measure to assess the difficulty level that an
invasion pathway can be identified. Results on both artificial and empirical
metapopulation networks show the robust performance in identifying actual
invasion pathways driving pandemic spread.
