Modified data fetching code in both Python and R.  Not tested properly, and DO need to do this to make sure that results make sense with fictitious data.
Major topics: filtering duplicates in Cad, buffering Camden boundary to avoid 'out of region' errors.
R code now computes STIK function for CAD data but it doesn't seem to make any sense.
New Python code not yet used, but hopefully much more efficient (and correct?).