How to calculate a best-fitting Euler pole from a fault trace?

Howdy,

I’m currently working on a detailed kinematic-palinspastic model of the North American Cordillera. In Alaska and western Canada there are many strike-slip faults (e.g. the Tintina fault system) that transect and offset accreted terranes and parts of the former continental margin.

As can be seen in the above screenshot, the traces of the strike-slip faults tend to follow arc segments of circles with various radii and center coordinates. Therefore, rotating one side of such a fault around the best-fit Euler pole of the fault trace, to bring piercing points into alignment, would closely approximate the pre-offset configuration (aside from some slight kinks associated with local zones of transpression and transtension).

Is there a tool, GPlates or otherwise, that would allow me to calculate a best-fitting Euler pole location and circle radius from the trace of a strike-slip fault (or a segment of such a fault)? Such a tool would save me a lot of time and manual adjustment. And if GPlates doesn’t currently have such a tool, I think it should be added in a future version.