The equatorialToHorizontal method looks great, but I want the inverse as well - given local horizontal coordinates (e.g. Altitude/Azimuth from an instrumented telescope), calculate the equatorial coordinates (RA/Dec), based on the site's location and time.
(I saw on @eteq mention on another issue possibly redirecting efforts to astropy, see also astropy/astropy-api#6 about adding similar functionality to astropy)