In [1]:

    

%matplotlib inline

import numpy as np
import matplotlib.pyplot as plt
from astropy.time import Time
from astropy.coordinates import SkyCoord, get_body


    

In [2]:

    

mars = get_body('mars', Time('2018-08-04 15:30'))
betAqr = SkyCoord.from_name('bet Aqr')


    

In [3]:

    

distance_deg = mars.separation(betAqr)
distance_deg.deg


    

    
Out[3]:





26.500244810281522

In [4]:

    

distance_pixels = 454
(distance_deg / distance_pixels * (640, 480)).deg


    

    
Out[4]:





array([ 37.3571733 ,  28.01787998])

In [5]:

    

camera_fov_long = 18.5
camera_fov_short = 14
(camera_fov_long * 2, camera_fov_short * 2)


    

    
Out[5]:





(37.0, 28)