Added new function to generate matrix of (from, to) counts by station