How to recalculate bone transforms/constraints to originate from armature origin point?
I found a solution for myself. However, it regrettably does use bpy.ops. If anyone has a solution which does not use bpy.ops I'd gladly take a look at that answer.
Here is the blend file for anyone who wants to take a look.