<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class="">Hi John,</div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">I solved a similar problem using the RA, DEC, and time from the header (and knowing the telescopes location) then just using the RA-LST method.  Here’s some (pseudo)code which calculates the HA for the mid point of the exposure (adjust as needed for start time if you like):</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><div class=""><font face="Monaco" class=""># Location of the observation</font></div><div class=""><font face="Monaco" class="">loc = EarthLocation(lat=cfg_loc['latitude’], lon=cfg_loc['longitude’], height=cfg_loc['elevation’])</font></div><div class=""><span style="font-family: Monaco;" class=""># Read the start time and exposure time of the image from the header</span></div><div class=""><font face="Monaco" class="">starttime = Time(time.strptime(header['DATE-OBS’], </font><span style="font-family: Monaco;" class="">'%Y-%m-%dT%H:%M:%S'), location=loc)</span></div><div class=""><font face="Monaco" class="">exptime = TimeDelta(float(header['EXPTIME']), format='sec')</font></div><div class=""><font face="Monaco" class=""># estimate the exposure mid time</font></div><div class=""><font face="Monaco" class="">midtime = starttime + exptime/2.0</font></div><div class=""><font face="Monaco" class=""># determine the sidereal time at the mid point of the exposure</font></div><div class=""><font face="Monaco" class="">sidereal = midtime.sidereal_time('apparent')</font></div><div class=""><font face="Monaco" class=""># read the position of the field from the header (assumes J2000)</font></div><div class=""><font face="Monaco" class="">header_pointing = SkyCoord(ra=float(header['RA'])*u.degree, </font><span style="font-family: Monaco;" class="">dec=float(header['DEC'])*u.degree)</span></div><div class=""><span style="font-family: Monaco;" class="">## To calculate HA, we need to precess to the current equinox</span></div><div class=""><font face="Monaco" class="">## otherwise the RA - LST method will be off.</font></div><div class=""><font face="Monaco" class="">FK5_Jnow = FK5(equinox=midtime)</font></div><div class=""><font face="Monaco" class="">HA = header_pointing.transform_to(FK5_Jnow).ra.to(u.hourangle) - sidereal</font></div></div><div class=""><br class=""></div><div class="">This was good enough for my purposes, but did seem a bit cumbersome.</div><div class=""><br class=""></div><div class="">-Josh</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><br class=""><div><blockquote type="cite" class=""><div class="">On Nov 4, 2016, at 2:04 PM, John K. Parejko <<a href="mailto:parejkoj@uw.edu" class="">parejkoj@uw.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="">Hello astropythians,<br class=""><br class="">I’m trying to work around a distinct lack of UT1 and LST in some telescope file headers, and was wondering if there’s a straight-forward astropy way to get a source’s Hour Angle from the alt/az, ra/dec, and telescope location. I can build an altaz SkyCoord, but I’m not sure what to do with it from there:<br class=""><br class="">loc = coord.EarthLocation.from_geodetic(long, lat, elevation)<br class="">sky = coord.SkyCoord(az, alt, location=loc, frame='altaz', unit='deg’)<br class=""><br class="">I can’t transform this to ICRS:<br class=""><br class="">ipdb> sky.transform_to(coord.ICRS)<br class="">*** TypeError: unsupported operand type(s) for -: 'NoneType' and ‘float’<br class=""><br class="">and there’s not an “hour_angle” method that I can see in either of those.<br class=""><br class="">Per the internet[1], one should be able to get HA from alt/az and dec directly (modulo refraction corrections):<br class=""><br class="">HA = asin(- sin(az)*cos(alt) / cos(dec))<br class=""><br class="">I can do that, but I was hoping for a nice and tidy astropy method.<br class=""><br class="">Thanks,<br class="">John<br class=""><br class="">1: <a href="http://star-www.st-and.ac.uk/~fv/webnotes/chapter7.htm" class="">http://star-www.st-and.ac.uk/~fv/webnotes/chapter7.htm</a><br class=""><br class="">--<br class="">*************************<br class="">John Parejko<br class=""><a href="mailto:parejkoj@uw.edu" class="">parejkoj@uw.edu</a><br class="">http://staff.washington.edu/parejkoj/<br class="">Department of Physics and Astronomy<br class="">University of Washington<br class="">Seattle, WA<br class="">**************************<br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class=""><br class="">_______________________________________________<br class="">AstroPy mailing list<br class="">AstroPy@scipy.org<br class="">https://mail.scipy.org/mailman/listinfo/astropy<br class=""></div></div></blockquote></div><br class=""></body></html>