def get_circle(ra, dec, rad_deg, n=100):
th = np.linspace(np.radians(rad_deg), np.radians(rad_deg), n)
phi = np.linspace(0, 2 * np.pi, n)
v = Vector3D.createThetaPhi(th, phi)
v.rotatey(np.pi / 2. - np.radians(dec))
v.rotatez(np.radians(ra))
return np.degrees(v.lon()), np.degrees(v.lat())
jint.compute()
sys.exit(0)
#z = np.ones(shape=(360,360))
#hc = Histogram2D(theta_edges,phi_edges)
for i0, th in enumerate(theta):
jtot = integrate(lambda t: jspline(t) * np.sin(t),
np.radians(theta_edges[i0]),
np.radians(theta_edges[i0 + 1]), 100)
# jval = jspline(np.radians(th))*costh_width[i0]
v = Vector3D.createThetaPhi(np.radians(th), np.radians(phi))
v.rotate(yaxis)
lat = np.degrees(v.lat())
lon = np.degrees(v.phi())
src_msk = len(lat) * [True]
if not sources is None:
for k in range(len(v.lat())):
p = Vector3D(v._x[:, k])
sep = np.degrees(p.separation(sources))
imin = np.argmin(sep)
minsep = sep[imin]
def get_circle(ra,dec,rad_deg,n=100):
th = np.linspace(np.radians(rad_deg),
np.radians(rad_deg),n)
phi = np.linspace(0,2*np.pi,n)
v = Vector3D.createThetaPhi(th,phi)
v.rotatey(np.pi/2.-np.radians(dec))
v.rotatez(np.radians(ra))
return np.degrees(v.lon()), np.degrees(v.lat())
jint.compute()
sys.exit(0)
#z = np.ones(shape=(360,360))
#hc = Histogram2D(theta_edges,phi_edges)
for i0, th in enumerate(theta):
jtot = integrate(lambda t: jspline(t)*np.sin(t),
np.radians(theta_edges[i0]),
np.radians(theta_edges[i0+1]),100)
# jval = jspline(np.radians(th))*costh_width[i0]
v = Vector3D.createThetaPhi(np.radians(th),np.radians(phi))
v.rotate(yaxis)
lat = np.degrees(v.lat())
lon = np.degrees(v.phi())
src_msk = len(lat)*[True]
if not sources is None:
for k in range(len(v.lat())):
p = Vector3D(v._x[:,k])
sep = np.degrees(p.separation(sources))
imin = np.argmin(sep)
minsep = sep[imin]
