As discussed on Slack with Tim Jenness, we should
- convert SkyCoord objects to ICRs (via .transform_to("icrs"))
- interpret all other coordinate objects as already being in the right frame and just extract their angles directly
It'd be ideal if we could make this work on arrays/vectors as well.
This will be especially valuable for a C++ version of obs_metadata's ObservationInfo.