Skip to content
Snippets Groups Projects
Commit 08090e13 authored by Taya Snijder's avatar Taya Snijder
Browse files

update

parent f4bd89b9
No related branches found
No related tags found
1 merge request!194Resolve L2SS-479 "09 dec 2021 convert pointing to delays"
...@@ -7,10 +7,11 @@ def subtract(a, b): ...@@ -7,10 +7,11 @@ def subtract(a, b):
return numpy.array([x - y for x, y in zip(a, b)]) return numpy.array([x - y for x, y in zip(a, b)])
class tempname: class delay_calculator:
def __init__(self, itrf: list([float])): def __init__(self, itrf: list([float])):
""" Return a measure object, configured for the specified location. """ """ Create a measure object, configured for the specified terrestrial location. """
measure = casacore.measures.measures() measure = casacore.measures.measures()
frame_location = measure.position("ITRF", *[f"{x}m" for x in itrf]) frame_location = measure.position("ITRF", *[f"{x}m" for x in itrf])
...@@ -55,29 +56,7 @@ class tempname: ...@@ -55,29 +56,7 @@ class tempname:
reference_dir_vector = self.get_direction_vector(pointing) reference_dir_vector = self.get_direction_vector(pointing)
# # compute the delays for an antennas w.r.t. the reference position # # compute the delays for an antennas w.r.t. the reference position
antenna_relative_itrf = [subtract(pos, reference_itrf) for pos in antenna_itrf] antenna_relative_itrf = [subtract(pos, self.reference_itrf) for pos in antenna_itrf]
delays = [self._get_delay(reference_dir_vector, relative_itrf) for relative_itrf in antenna_relative_itrf] delays = [self._get_delay(reference_dir_vector, relative_itrf) for relative_itrf in antenna_relative_itrf]
return delays return delays
if __name__ == '__main__':
# # create a frame tied to the reference position
reference_itrf = [3826577.066, 461022.948, 5064892.786] # CS002LBA, in ITRF2005 epoch 2012.5
d = tempname(reference_itrf)
# # set the timestamp to solve for
timestamp = datetime.datetime(2021,1,1,0,0,5)
d.set_measure_time(timestamp)
# compute the delays for an antennas w.r.t. the reference position
antenna_itrf = [[3826923.546, 460915.441, 5064643.489]] # CS001LBA, in ITRF2005 epoch 2012.5
# # obtain the direction vector for a specific pointing
direction = "J2000","0deg","0deg"
# calculate the delays based on the set reference position, the set time and now the set direction and antenna positions.
delays = d.convert(direction, antenna_itrf)
# print the delays
import pprint
pprint.pprint(delays)
from delays import *
if __name__ == '__main__':
# # create a frame tied to the reference position
reference_itrf = [3826577.066, 461022.948, 5064892.786] # CS002LBA, in ITRF2005 epoch 2012.5
d = delay_calculator(reference_itrf)
# # set the timestamp to solve for
timestamp = datetime.datetime(2021,1,1,0,0,5)
d.set_measure_time(timestamp)
# compute the delays for an antennas w.r.t. the reference position
antenna_itrf = [[3826923.546, 460915.441, 5064643.489]] # CS001LBA, in ITRF2005 epoch 2012.5
# # obtain the direction vector for a specific pointing
direction = "J2000","0deg","0deg"
# calculate the delays based on the set reference position, the set time and now the set direction and antenna positions.
delays = d.convert(direction, antenna_itrf)
# print the delays
import pprint
pprint.pprint(delays)
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment