Sky

sky = Sky(ra_range, dec_range, z_range)

Function	Result
sky.boxsize	Tuple $\in \mathbb(R)^3$; size of bounding box of the volume in the sky
sky.left	Tuple $\in \mathbb(R)^3$; left-bottom corner of the bounding box
sky.right	Tuple $\in \mathbb(R)^3$; upper-right corner of the bounding box



In [6]:

    
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import mockgallib as mock



In [5]:

    
# Sky
ra_range  = [30.1, 38.9]
dec_range = [-6.0, -4.1]
z_range   = [0.38889, 1.21239]
sky = mock.Sky(ra_range, dec_range, z_range)
r_range   = sky.r;

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')

n = 20
ras  = np.linspace(ra_range[0], ra_range[1], n)
decs = np.linspace(dec_range[0], dec_range[1], n)
rs   = np.linspace(r_range[0], r_range[1], 2)


def draw_line(r_list, ra_list, dec_list):
    x3 = [ sky.compute_x(r, ra, dec) for r in r_list for ra in ra_list for dec in dec_list ]

    x = [ x_ for (x_, y_, z_) in x3 ]
    y = [ y_ for (x_, y_, z_) in x3 ]
    z = [ z_ for (x_, y_, z_) in x3 ]
    
    ax.plot(x,y,z, color='red')

draw_line([r_range[0]], ras, [dec_range[0]])
draw_line([r_range[0]], ras, [dec_range[1]])
draw_line([r_range[1]], ras, [dec_range[0]])
draw_line([r_range[1]], ras, [dec_range[1]])

draw_line([r_range[0]], [ra_range[0]], decs)
draw_line([r_range[0]], [ra_range[1]], decs)
draw_line([r_range[1]], [ra_range[0]], decs)
draw_line([r_range[1]], [ra_range[1]], decs)

draw_line(rs, [ra_range[0]], [dec_range[0]])
draw_line(rs, [ra_range[1]], [dec_range[0]])
draw_line(rs, [ra_range[0]], [dec_range[1]])
draw_line(rs, [ra_range[1]], [dec_range[1]])

print('Close the window to continue.')
plt.show();









    



Close the window to continue.



In [1]:

    
%%html
<style>table {float:left}</style>