In [45]:

    

import numpy as np
import astropy.units as u


    

In [69]:

    

pixscale = 0.09639 * 0.873


    

In [116]:

    

cen = np.array([425, 435])
elp10 = np.array([503, 347])
azp10 = np.array([511, 519])


    

In [71]:

    

elp10, azp10


    

    
Out[71]:





(array([ 78, -88]), array([86, 84]))

In [72]:

    

d1 = np.sqrt(elp10.dot(elp10)) * pixscale
d1


    

    
Out[72]:





9.8952303695365256

In [73]:

    

d2 = np.sqrt(azp10.dot(azp10)) * pixscale
d2


    

    
Out[73]:





10.11603213855525

In [74]:

    

ave = (d1 + d2)/2
ave


    

    
Out[74]:





10.005631254045888

In [68]:

    

10/ave * 0.883


    

    
Out[68]:





0.87250866820321082

In [42]:

    

pixscale


    

    
Out[42]:





0.08511237

In [47]:

    

a1 = np.arctan2(76, -88) * u.radian


    

In [48]:

    

a1.to(u.deg)


    

    
Out[48]:





$139.18492 \; \mathrm{{}^{\circ}}$

In [50]:

    

a2 = np.arctan2(84, 84) * u.radian


    

In [51]:

    

a2.to(u.deg)


    

    
Out[51]:





$45 \; \mathrm{{}^{\circ}}$

In [61]:

    

from mmtwfs.wfs import pol2cart, cart2pol


    

In [243]:

    

xref = cen[0]
yref = cen[1]
xc = elp10[0]
yc = elp10[1]
dx = -(xc - xref) * pixscale
dy = (yc - yref) * pixscale


    

In [244]:

    

rotation = -225 * u.deg


    

In [245]:

    

total_rotation = rotation.to(u.radian).value


    

In [246]:

    

dr, phi = cart2pol([dx, dy])


    

In [247]:

    

derot_phi = phi - total_rotation


    

In [248]:

    

az, el = pol2cart([dr, derot_phi])


    

In [249]:

    

az, el


    

    
Out[249]:





(-0.59501953763472659, 9.8773243247364775)

In [253]:

    

(747 * u.mm).to(u.imperial.inch)


    

    
Out[253]:





$29.409449 \; \mathrm{inch}$

In [ ]: