In [2]:
import numpy as np
import astropy.units as u
from mmtwfs.wfs import cart2pol, pol2cart
In [3]:
# 2019-01-25 data
center = (251, 267)
p10az = (177, 254)
p10el = (239, 347)
sysrot = (234-90) * u.deg
rot = -46.93 * u.deg
az_parity = -1
el_parity = -1
pix_size = 0.135
In [4]:
# 2019-01-26 data
center = (241, 265)
p10az = (185, 318)
p10el = (289, 327)
sysrot = (234-90) * u.deg
rot = 0.0 * u.deg
az_parity = -1
el_parity = -1
pix_size = 0.135
In [6]:
total_rot = sysrot - rot
dx = p10el[0] - center[0]
dy = p10el[1] - center[1]
dr, phi = cart2pol([dx, -dy])
derot_phi = phi + total_rot.to(u.rad).value
az, el = pol2cart([dr, derot_phi])
az *= az_parity * pix_size
el *= el_parity * pix_size
az, el
Out[6]:
In [ ]: