def export(start_datetime, end_datetime, num_samples):
dt = (end_datetime - start_datetime) / num_samples
for i in range(num_samples):
t = start_datetime + i * dt
ts = calendar.timegm(t.utctimetuple())
moon.compute(t)
ecliptic_pos = Ecliptic(moon)
positions.append(
dict(timestamp=ts, lon=ecliptic_pos.lon, lat=ecliptic_pos.lat))
def _get_radec_from_ecliptic(self, elong, elat):
# convert via pyephem
try:
ec = Ecliptic(elong, elat)
# check for B name
if "B" in self.name:
epoch = "1950"
else:
epoch = "2000"
eq = Equatorial(ec, epoch=str(epoch))
raj = np.double(eq.ra)
decj = np.double(eq.dec)
except TypeError:
msg = "WARNING: Cannot fine sky location coordinates "
msg += "for PSR {0}. ".format(self.name)
msg += "Setting values to 0.0"
logger.warning(msg)
raj = 0.0
decj = 0.0
return raj, decj
def _get_radec_from_ecliptic(self, elong, elat):
# convert via pyephem
try:
ec = Ecliptic(elong, elat)
# check for B name
if 'B' in self.name:
epoch = '1950'
else:
epoch = '2000'
eq = Equatorial(ec, epoch=str(epoch))
raj = np.double(eq.ra)
decj = np.double(eq.dec)
# TODO: should use logging here
except TypeError:
msg = 'WARNING: Cannot fine sky location coordinates '
msg += 'for PSR {0}. '.format(self.name)
msg += 'Setting values to 0.0'
print(msg)
raj = 0.0
decj = 0.0
return raj, decj