when using od_aggregate, st_join sometimes joins 2 polygons to an o or d, due to the fact that aggzones are overlapping, and the o or d resides in both.
the function should throw an error when aggzones are intersecting.
I'll supply a reprex if requested