def hba_dipole_etrs(self, field_name):
"""Return a list of all ETRS dipole coordinates for a given HBA antenna field
Args:
field_name (str): Field name (e.g. 'CS001HBA0')
Returns:
array: array of ETRS coordinates
Example:
>>> import lofarantpos.db
>>> import numpy
>>> if float(".".join(numpy.__version__.split('.')[:2]))>=1.14: numpy.set_printoptions(legacy=True)
>>> db = lofarantpos.db.LofarAntennaDatabase()
>>> db.hba_dipole_etrs("IE613HBA")[:5]
array([[ 3801679.57332..., -528959.80788..., 5076969.80405...],
[ 3801680.56727..., -528959.55814..., 5076969.08837...],
[ 3801681.56122..., -528959.3083985 , 5076968.37269...],
[ 3801682.55516..., -528959.05865..., 5076967.65701...],
[ 3801679.71139..., -528961.02799..., 5076969.57003...]])
"""
return geo.transform(
self.hba_dipole_pqr(field_name),
numpy.zeros(3),
self.pqr_to_etrs[field_name]) + \
self.phase_centres[field_name][numpy.newaxis, :]
def hba_dipole_etrs(self, field_name):
"""Return a list of all ETRS dipole coordinates for a given HBA antenna field
Args:
field_name (str): Field name (e.g. 'CS001HBA0')
Returns:
array: array of ETRS coordinates
Example:
>>> import lofarantpos.db
>>> import numpy
>>> db = lofarantpos.db.LofarAntennaDatabase()
>>> db.hba_dipole_etrs("IE613HBA")[:5]
array([[3801679.57332033, -528959.80788382, 5076969.80405122],
[3801680.56726901, -528959.55814198, 5076969.08837304],
[3801681.56121763, -528959.30839824, 5076968.37269509],
[3801682.55516625, -528959.0586545 , 5076967.65701715],
[3801679.7113895 , -528961.02799576, 5076969.57003303]])
"""
return geo.transform(
self.hba_dipole_pqr(field_name),
numpy.zeros(3),
self.pqr_to_etrs[field_name]) + \
self.phase_centres[field_name][numpy.newaxis, :]
def antenna_pqr(self, field_name):
"""Return a list of all PQR antenna coordinates for a given antenna field
Args:
field_name (str): Field name (e.g. 'CS001HBA0')
Returns:
array: array of PQR coordinates
"""
return geo.transform(self.antenna_etrs(field_name),
self.phase_centres[field_name],
self.pqr_to_etrs[field_name].T)