the situation is pretty problematic, I hope someday scipy will add a python package for rotating and interpolating quaternions, up to now: