Restructured TS map computation

As close to final as it could get today.
Would it be much work to put this radius computation into a function in `test_statistic.py`?
(doesn't need to be public with a docstring for now).