ZDES Function
Equivalent to running a P-Median algorithm for the data zones assigned to each of the M regions and then summing the weighted distances (i.e. mean travel cost minimised)
Works well because each change to region boundaries only involves re-computing 2 P-Medians from starting locations that are very close to the optimal